Showing 1–43 of 43 results for author: Williams, R J R

Generation of photoionized plasmas in the laboratory of relevance to accretion-powered x-ray sources using keV line radiation

Authors: D. Riley, R. L. Singh, S White, M. Charlwood, D. Bailie, C. Hyland, T. Audet, G. Sarri, B. Kettle, G. Gribakin, S. J. Rose, E. G. Hill, G. J. Ferland, R. J. R. Williams, F. P. Keenan

Abstract: We describe laboratory experiments to generate X-ray photoionized plasmas of relevance to accretion-powered X-ray sources such as neutron star binaries and quasars, with significant improvements over previous work. A key quantity is referenced, namely the photoionization parameter. This is normally meaningful in an astrophysical steady-state context, but is also commonly used in the literature as… ▽ More We describe laboratory experiments to generate X-ray photoionized plasmas of relevance to accretion-powered X-ray sources such as neutron star binaries and quasars, with significant improvements over previous work. A key quantity is referenced, namely the photoionization parameter. This is normally meaningful in an astrophysical steady-state context, but is also commonly used in the literature as a figure of merit for laboratory experiments that are, of necessity, time-dependent. We demonstrate emission-weighted values of ξ > 50 ergcm/s using laser-plasma X-ray sources, with higher results at the centre of the plasma which are in the regime of interest for several astrophysical scenarios. Comparisons of laboratory experiments with astrophysical codes are always limited, principally by the many orders of magnitude differences in time and spatial scales, but also other plasma parameters. However useful checks on performance can often be made for a limited range of parameters. For example, we show that our use of a keV line source, rather than the quasi-blackbody radiation fields normally employed in such experiments, has allowed the generation of the ratio of inner-shell to outer-shell photoionization expected from a blackbody source with ~keV spectral temperature. We compare calculations from our in-house plasma modelling code with those from Cloudy and find moderately good agreement for the time evolution of both electron temperature and average ionisation. However, a comparison of code predictions for a K-beta argon X-ray spectrum with experimental data reveals that our Cloudy simulation overestimates the intensities of more highly ionised argon species. This is not totally surprising as the Cloudy model was generated for a single set of plasma conditions, while the experimental data are spatially integrated. △ Less

Submitted 19 March, 2024; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: High Energy Density Physics, in press; 21 pages, 10 figures

H, He-like recombination spectra IV; H, He-like recombination spectra IV: clarification and refinement of methodology for $l$-changing collisions

Authors: N. R. Badnell, F. Guzmán, S. Brodie, R. J. R. Williams, P. A. M. van Hoof, M. Chatzikos, G. J. Ferland

Abstract: Precise spectral diagnostic modelling of H~{\sc i} and He~{\sc ii} recombination spectra can constrain theoretical models which describe many astrophysical environments. Simple analytic expressions are of interest for collisional $l$-changing rate coefficients that are used by large-scale population modelling codes. We review, clarify and improve-upon the modified Pengelly \& Seaton formulae of Gu… ▽ More Precise spectral diagnostic modelling of H~{\sc i} and He~{\sc ii} recombination spectra can constrain theoretical models which describe many astrophysical environments. Simple analytic expressions are of interest for collisional $l$-changing rate coefficients that are used by large-scale population modelling codes. We review, clarify and improve-upon the modified Pengelly \& Seaton formulae of Guzmán \etal We show that the recent poor results for it shown by Vrinceanu \etal are due to their misinterpretation of its usage. We also detail efficient numerical algorithms which should enable the full quantum mechanical expression for such rate coefficients to be used much more routinely by modelling codes. We illustrate with some collisional-radiative population modelling for hydrogen. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: accepted in MNRAS

Space Telescope and Optical Reverberation Mapping Project. X. Understanding the Absorption-Line Holiday in NGC 5548

Authors: M. Dehghanian, G. J. Ferland, G. A. Kriss, B. M. Peterson, S. Mathur, M. Mehdipour, F. Guzman, M. Chatzikos, P. A. M. Van Hoof, R. J. R. Williams, N. Arav, A. J. Barth, M. C. Bentz, S. Bisogni, W. N. Brandt, D. M. Crenshaw, E. Dalla Bonta, G. De Rosa, M. M. Fausnaugh, J. M. Gelbord, M. R. Goad, A. Gupta, Keith Horne, J. Kaastra, C. Knigge , et al. (5 additional authors not shown)

Abstract: The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central bla… ▽ More The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central black hole mass in AGNs. The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is the most intensive multi-wavelength AGN monitoring campaign ever. For most of the campaign, the emission-line variations followed changes in the continuum with a time lag, as expected. However, the lines varied independently of the observed UV-optical continuum during a 60 -- 70 day "holiday." To understand this remarkable phenomenon, we study the intrinsic absorption lines present in NGC 5548. We identify a novel cycle that reproduces the absorption line variability and thus identify the physics that allows the holiday to occur. In our model, variations in this obscurer's line-of-sight covering factor modify the soft X-ray continuum. This leads to changes in the ionization of helium gas in the broad-line region. Ionizing radiation produced by recombining helium then affects the ionization of other species as observed during the AGN STORM holiday. It is likely that any other model which selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission line behavior observed. △ Less

Submitted 24 April, 2019; v1 submitted 30 December, 2018; originally announced December 2018.

Comments: This is the version accepted by ApJ. Will be published shortly

The classical D-type expansion of spherical H II regions

Authors: Robin J. R. Williams, Thomas G. Bisbas, Thomas J. Haworth, Jonathan Mackey

Abstract: Recent numerical and analytic work has highlighted some shortcomings in our understanding of the dynamics of H II region expansion, especially at late times, when the H II region approaches pressure equilibrium with the ambient medium. Here we reconsider the idealized case of a constant radiation source in a uniform and spherically symmetric ambient medium, with an isothermal equation of state. A… ▽ More Recent numerical and analytic work has highlighted some shortcomings in our understanding of the dynamics of H II region expansion, especially at late times, when the H II region approaches pressure equilibrium with the ambient medium. Here we reconsider the idealized case of a constant radiation source in a uniform and spherically symmetric ambient medium, with an isothermal equation of state. A thick-shell solution is developed which captures the stalling of the ionization front and the decay of the leading shock to a weak compression wave as it escapes to large radii. An acoustic approximation is introduced to capture the late-time damped oscillations of the H II region about the stagnation radius. Putting these together, a matched asymptotic equation is derived for the radius of the ionization front which accounts for both the inertia of the expanding shell and the finite temperature of the ambient medium. The solution to this equation is shown to agree very well with the numerical solution at all times, and is superior to all previously published solutions. The matched asymptotic solution can also accurately model the variation of H II region radius for a time-varying radiation source. △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: 8 pages, 5 figures. Authors' original version, revised version with minor changes accepted for publication in MNRAS

The 2017 Release of Cloudy

Authors: G. J. Ferland, M. Chatzikos, F. Guzmán, M. L. Lykins, P. A. M. van Hoof, R. J. R. Williams, N. P. Abel, N. R. Badnell, F. P. Keenan, R. L. Porter, P. C. Stancil

Abstract: We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully… ▽ More We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds. △ Less

Submitted 5 October, 2023; v1 submitted 30 May, 2017; originally announced May 2017.

Comments: replace Fig 6 with the published version

Journal ref: 2017RMxAA..53..385F

Thermodynamically-consistent semi-classical $\ell$-changing rates

Authors: R. J. R. Williams, F. Guzmán, N. R. Badnell, P. A. M. van Hoof, M. Chatzikos, G. J. Ferland

Abstract: We compare the results of the semi-classical (SC) and quantum-mechanical (QM) formalisms for angular-momentum changing transitions in Rydberg atom collisions given by Vrinceanu & Flannery, J. Phys. B 34, L1 (2001), and Vrinceanu, Onofrio & Sadeghpour, ApJ 747, 56 (2012), with those of the SC formalism using a modified Monte Carlo realization. We find that this revised SC formalism agrees well with… ▽ More We compare the results of the semi-classical (SC) and quantum-mechanical (QM) formalisms for angular-momentum changing transitions in Rydberg atom collisions given by Vrinceanu & Flannery, J. Phys. B 34, L1 (2001), and Vrinceanu, Onofrio & Sadeghpour, ApJ 747, 56 (2012), with those of the SC formalism using a modified Monte Carlo realization. We find that this revised SC formalism agrees well with the QM results. This provides further evidence that the rates derived from the QM treatment are appropriate to be used when modelling recombination through Rydberg cascades, an important process in understanding the state of material in the early universe. The rates for $Δ\ell=\pm1$ derived from the QM formalism diverge when integrated to sufficiently large impact parameter, $b$. Further to the empirical limits to the $b$ integration suggested by Pengelly & Seaton, MNRAS 127, 165 (1964), we suggest that the fundamental issue causing this divergence in the theory is that it does not fully cater for the finite time taken for such distant collisions to complete. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 12 pages, 4 figures. To be published in J Phys B. Author original version: journal accepted version, with minor changes, will be posted after embargo period

Journal ref: J. Phys. B: At. Mol. Opt. Phys. (2017) 50, 115201

Testing atomic collision theory with the two-photon continuum of astrophysical nebulae

Authors: F. Guzmán, N. R. Badnell, M. Chatzikos, P. A. M. van Hoof, R. J. R. Williams, G. J. Ferland

Abstract: Accurate rates for energy-degenerate l-changing collisions are needed to determine cosmological abundances and recombination. There are now several competing theories for the treatment of this process, and it is not possible to test these experimentally. We show that the H I two-photon continuum produced by astrophysical nebulae is strongly affected by l-changing collisions. We perform an analysis… ▽ More Accurate rates for energy-degenerate l-changing collisions are needed to determine cosmological abundances and recombination. There are now several competing theories for the treatment of this process, and it is not possible to test these experimentally. We show that the H I two-photon continuum produced by astrophysical nebulae is strongly affected by l-changing collisions. We perform an analysis of the different underlying atomic processes and simulate the recombination and two-photon spectrum of a nebula containing H and He. We provide an extended set of effective recombination coefficients and updated l-changing 2s-2p transition rates using several competing theories. In principle, accurate astronomical observations could determine which theory is correct. △ Less

Submitted 26 January, 2017; originally announced January 2017.

Comments: 7 pages, 6 fig, 2 tables, accepted in MNRAS

Journal ref: Mon Not R Astron Soc (2017) 467 (4): 3944-3950

H, He-like recombination spectra II: $l$-changing collisions for He Rydberg states

Authors: F. Guzmán, N. R. Badnell, R. J. R. Williams, P. A. M. van Hoof, M. Chatzikos, G. J. Ferland

Abstract: Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which foc… ▽ More Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $\sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy. △ Less

Submitted 31 August, 2016; originally announced August 2016.

Comments: 9 pages, $ figures, 2 tables, accepted for publication in MNRAS

H, He-like recombination spectra I: $l$-changing collisions for hydrogen

Authors: F. Guzmán, N. R. Badnell, R. J. R. Williams, P. A. M. van Hoof, M. Chatzikos, G. J. Ferland

Abstract: Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe di… ▽ More Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe dipole $l$-changing collisions ($nl\rightarrow nl'=l\pm 1$) between energy-degenerate states within an $n$-shell. The work of Pengelly & Seaton (1964) has, for half-a-century, been considered the definitive study which "solved" the problem. Recent work by Vrinceanu et al.(2012) recommended the use of rate coefficients from a semi-classical approximation which are nearly an order of magnitude smaller than those of Pengelly & Seaton (1964), with the result that significantly higher densities are needed for the $nl$ populations to come into local thermodynamic equilibrium. Here, we compare predicted H~I emissivities from the two works and find widespread differences, of up to $\approx 10$%. This far exceeds the 1% precision required to obtain the primordial He/H abundance ratio from observations so as to constrain Big Bang cosmologies. We recommend using the rate coefficients of Pengelly & Seaton (1964) for $l$-changing collisions, to describe the H recombination spectrum, based-on their quantum mechanical representation of the long-range dipole interaction. △ Less

Submitted 13 April, 2016; originally announced April 2016.

Comments: 6 pages, 3 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Ultraviolet emission lines of Si II in cool star and solar spectra

Authors: Sibasish Laha, Francis P. Keenan, Gary J. ferland, Catherine A. Ramsbottom, Kanti M. Aggarwal, Thomas R. Ayres, Marios Chatzikos, Peter A. M. van Hoof, Robin J. R. Williams

Abstract: Recent atomic physics calculations for Si II are employed within the Cloudy modelling code to analyse Hubble Space Telescope (HST) STIS ultraviolet spectra of three cool stars, Beta-Geminorum, Alpha-Centauri A and B, as well as previously published HST/GHRS observations of Alpha-Tau, plus solar quiet Sun data from the High Resolution Telescope and Spectrograph. Discrepancies found previously betwe… ▽ More Recent atomic physics calculations for Si II are employed within the Cloudy modelling code to analyse Hubble Space Telescope (HST) STIS ultraviolet spectra of three cool stars, Beta-Geminorum, Alpha-Centauri A and B, as well as previously published HST/GHRS observations of Alpha-Tau, plus solar quiet Sun data from the High Resolution Telescope and Spectrograph. Discrepancies found previously between theory and observation for line intensity ratios involving the 3s$^{2}$3p $^{2}$P$_{J}$--3s3p$^{2}$ $^{4}$P$_{J^{\prime}}$ intercombination multiplet of Si II at 2335 Angs are significantly reduced, as are those for ratios containing the 3s$^{2}$3p $^{2}$P$_{J}$--3s3p$^{2}$ $^{2}$D$_{J^{\prime}}$ transitions at 1816 Angs. This is primarily due to the effect of the new Si II transition probabilities. However, these atomic data are not only very different from previous calculations, but also show large disagreements with measurements, specifically those of Calamai et. al. (1993) for the intercombination lines. New measurements of transition probabilities for Si II are hence urgently required to confirm (or otherwise) the accuracy of the recently calculated values. If the new calculations are confirmed, then a long-standing discrepancy between theory and observation will have finally been resolved. However, if the older measurements are found to be correct, then the agreement between theory and observation is simply a coincidence and the existing discrepancies remain. △ Less

Submitted 29 October, 2015; originally announced October 2015.

Comments: Accepted for publication in MNRAS

StarBench: The D-type expansion of an HII region

Authors: T. G. Bisbas, T. J. Haworth, R. J. R. Williams, J. Mackey, P. Tremblin, A. C. Raga, S. J. Arthur, C. Baczynski, J. E. Dale, T. Frostholm, S. Geen, T. Haugboelle, D. Hubber, I. T. Iliev, R. Kuiper, J. Rosdahl, D. Sullivan, S. Walch, R. Wuensch

Abstract: StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase o… ▽ More StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: 20 pages, 12 Figures, 4 Tables. Accepted for publication in MNRAS. Comments are welcome. Participation in future StarBench tests is also welcome

Stout: Cloudy's Atomic and Molecular Database

Authors: M. L. Lykins, G. J. Ferland, R. Kisielius, M. Chatzikos, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams, F. P. Keenan, P. C. Stancil

Abstract: We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures 2.8 K < T < 1e10 K and densities spanning the low to high-density limits. The radiation field between photon energies $10^{-8}$ Ry and 100 MeV is considered… ▽ More We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures 2.8 K < T < 1e10 K and densities spanning the low to high-density limits. The radiation field between photon energies $10^{-8}$ Ry and 100 MeV is considered, along with all atoms and ions of the lightest 30 elements, and ~100 molecules. For ease of maintenance, the data are stored in a format as close as possible to the original data sources. Few data sources include the full range of data we need. We describe how we fill in the gaps in the data or extrapolate rates beyond their tabulated range. We tabulate data sources both for the atomic spectroscopic parameters and for collision data for the next release of Cloudy. This is not intended as a review of the current status of atomic data, but rather a description of the features of the database which we will build upon. △ Less

Submitted 4 June, 2015; originally announced June 2015.

Comments: ApJ in press

Accurate determination of the free-free Gaunt factor. II - relativistic Gaunt factors

Authors: P. A. M. van Hoof, G. J. Ferland, R. J. R. Williams, K. Volk, M. Chatzikos, M. Lykins, R. L. Porter

Abstract: When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter spac… ▽ More When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter space than has ever been produced before. We present tables of the thermally averaged Gaunt factor covering the range log10(gamma^2) = -6 to 10 and log10(u) = -16 to 13 for all atomic numbers Z = 1 through 36. The data were calculated using the relativistic Bethe-Heitler-Elwert (BHE) approximation and were subsequently merged with accurate non-relativistic results in those parts of the parameter space where the BHE approximation is not valid. These data will be incorporated in the next major release of the spectral synthesis code Cloudy. We also produced tables of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to 10 for all values of Z between 1 and 36. All the data presented in this paper are available online. △ Less

Submitted 24 February, 2015; originally announced February 2015.

Comments: 8 pages, 8 figures, 2 tables

Journal ref: MNRAS, 449, 2112 (2015)

Collisional excitation of [C II], [O I] and CO in Massive Galaxies

Authors: R. E. A. Canning, G. J. Ferland, A. C. Fabian, R. M. Johnstone, P. A. M. van Hoof, R. L. Porter, N. Werner, R. J. R. Williams

Abstract: Many massive galaxies at the centres of relaxed galaxy clusters and groups have vast reservoirs of cool (~10,000 K) and cold (~100 K) gas. In many low redshift brightest group and cluster galaxies this gas is lifted into the hot ISM in filamentary structures, which are long lived and are typically not forming stars. Two important questions are how far do these reservoirs cool and if cold gas is ab… ▽ More Many massive galaxies at the centres of relaxed galaxy clusters and groups have vast reservoirs of cool (~10,000 K) and cold (~100 K) gas. In many low redshift brightest group and cluster galaxies this gas is lifted into the hot ISM in filamentary structures, which are long lived and are typically not forming stars. Two important questions are how far do these reservoirs cool and if cold gas is abundant what is the cause of the low star formation efficiency? Heating and excitation of the filaments from collisions and mixing of hot particles in the surrounding X-ray gas describes well the optical and near infra-red line ratios observed in the filaments. In this paper we examine the theoretical properties of dense, cold clouds emitting in the far infra-red and submillimeter through the bright lines of [C II]157 μm , [O I]63 μm and CO, exposed to these energetic ionising particles. While some emission lines may be optically thick we find this is not sufficient to model the emission line ratios. Models where the filaments are supported by thermal pressure support alone also cannot account for the cold gas line ratios but a very modest additional pressure support, either from turbulence or magnetic fields can fit the observed [O I]/[C II] line ratios by decreasing the density of the gas. This may also help stabilise the filaments against collapse leading to the low rates of star formation. Finally we make predictions for the line ratios expected from cold gas under these conditions and present diagnostic diagrams for comparison with further observations. We provide our code as an Appendix. △ Less

Submitted 6 January, 2015; originally announced January 2015.

Comments: 17 pages, submitted to MNRAS

Implications of Coronal Line Emission in NGC 4696

Authors: M. Chatzikos, R. J. R. Williams, G. J. Ferland, R. E. A. Canning, A. C. Fabian, J. S. Sanders, P. A. M. van Hoof, R. M. Johnstone, M. Lykins, R. L. Porter

Abstract: We announce a new facility in the spectral code CLOUDY that enables tracking the evolution of a cooling parcel of gas with time. For gas cooling from temperatures relevant to galaxy clusters, earlier calculations estimated the [Fe XIV] λ5303 / [Fe X] λ6375 luminosity ratio, a critical diagnostic of a cooling plasma, to slightly less than unity. By contrast, our calculations predict a ratio ~3. We… ▽ More We announce a new facility in the spectral code CLOUDY that enables tracking the evolution of a cooling parcel of gas with time. For gas cooling from temperatures relevant to galaxy clusters, earlier calculations estimated the [Fe XIV] λ5303 / [Fe X] λ6375 luminosity ratio, a critical diagnostic of a cooling plasma, to slightly less than unity. By contrast, our calculations predict a ratio ~3. We revisit recent optical coronal line observations along the X-ray cool arc around NGC 4696 by Canning et al. (2011), which detected [Fe X] λ6375, but not [Fe XIV] λ5303. We show that these observations are not consistent with predictions of cooling flow models. Differential extinction could in principle account for the observations, but it requires extinction levels (A_V > 3.625) incompatible with previous observations. The non-detection of [Fe XIV] implies a temperature ceiling of 2.1 million K. Assuming cylindrical geometry and transonic turbulent pressure support, we estimate the gas mass at ~1 million solar masses. The coronal gas is cooling isochorically. We propose that the coronal gas has not condensed out of the intracluster medium, but instead is the conductive or mixing interface between the X-ray plume and the optical filaments. We present a number of emission lines that may be pursued to test this hypothesis and constrain the amount of intermediate temperature gas in the system. △ Less

Submitted 16 October, 2014; originally announced October 2014.

Comments: 11 pages, 8 figures, 6 tables; accepted for publication in MNRAS

Accurate determination of the free-free Gaunt factor; I - non-relativistic Gaunt factors

Authors: P. A. M. van Hoof, R. J. R. Williams, K. Volk, M. Chatzikos, G. J. Ferland, M. Lykins, R. L. Porter, Y. Wang

Abstract: Modern spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt fa… ▽ More Modern spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt factors over a much wider range of parameters than has ever been produced before. We first produced a table of non-averaged Gaunt factors, covering the parameter space log10(epsilon_i) = -20 to +10 and log10(w) = -30 to +25. We then continued to produce a table of thermally averaged Gaunt factors covering the parameter space log10(gamma^2) = -6 to +10 and log10(u) = -16 to +13. Finally we produced a table of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to +10. All the data presented in this paper are available online. △ Less

Submitted 17 August, 2014; v1 submitted 18 July, 2014; originally announced July 2014.

Comments: 10 pages, 5 tables, 3 figures. Fixed typo in Eq. 14

Journal ref: MNRAS, 444, 420 (2014)

Radiative Cooling II: Effects of Density and Metallicity

Authors: Ye Wang, G. J. Ferland, M. L. Lykins, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams

Abstract: This work follows Lykins et al. discussion of classic plasma cooling function at low density and solar metallicity. Here we focus on how the cooling function changes over a wide range of density (n_H<10^12 cm^(-3)) and metallicity (Z<30Z _sun ). We find that high densities enhance the ionization of elements such as hydrogen and helium until they reach local thermodynamic equilibrium. By charge tra… ▽ More This work follows Lykins et al. discussion of classic plasma cooling function at low density and solar metallicity. Here we focus on how the cooling function changes over a wide range of density (n_H<10^12 cm^(-3)) and metallicity (Z<30Z _sun ). We find that high densities enhance the ionization of elements such as hydrogen and helium until they reach local thermodynamic equilibrium. By charge transfer, the metallicity changes the ionization of hydrogen when it is partially ionized. We describe the total cooling function as a sum of four parts: those due to H&He, the heavy elements, electron-electron bremsstrahlung and grains. For the first 3 parts, we provide a low-density limit cooling function, a density dependence function, and a metallicity dependence function. These functions are given with numerical tables and analytical fit functions. For grain cooling, we only discuss in ISM case. We then obtain a total cooling function that depends on density, metallicity and temperature. As expected, collisional de-excitation suppresses the heavy elements cooling. Finally, we provide a function giving the electron fraction, which can be used to convert the cooling function into a cooling rate. △ Less

Submitted 8 May, 2014; v1 submitted 12 March, 2014; originally announced March 2014.

Comments: Published by MNRAS, online tables can be found at http://mnras.oxfordjournals.org/content/suppl/2014/04/28/stu514.DC1

Journal ref: 2014, MNRAS 440, 3100

Expanded Iron UTA spectra -- probing the thermal stability limits in AGN clouds

Authors: G. J. Ferland, R. Kisielius, F. P. Keenan, P. A. M. van Hoof, V. Jonauskas, M. L. Lykins, R. L. Porter, R. J. R. Williams

Abstract: The Fe unresolved transition array (UTAs) produce prominent features in the 15-17?A wavelength range in the spectra of Active Galactic Nuclei (AGN). Here we present new calculations of the energies and oscillator strengths of inner- shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation… ▽ More The Fe unresolved transition array (UTAs) produce prominent features in the 15-17?A wavelength range in the spectra of Active Galactic Nuclei (AGN). Here we present new calculations of the energies and oscillator strengths of inner- shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in Active Galactic Nuclei. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration. △ Less

Submitted 1 March, 2013; originally announced March 2013.

Comments: ApJ in press

The 2013 Release of Cloudy

Authors: G. J. Ferland, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams, N. P. Abel, M. L. Lykins, Gargi Shaw, W. J. Henney, P. C. Stancil

Abstract: This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime.… ▽ More This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered. △ Less

Submitted 18 February, 2013; originally announced February 2013.

Comments: Major review article to be published in RevMexAA. 27 pages, 18 figures

Radiative cooling in collisionally and photo ionized plasmas

Authors: M. L. Lykins, G. J. Ferland, Ryan L. Porter, Peter A. M. van Hoof, R. J. R. Williams, Orly Gnat

Abstract: We discuss recent improvements in the calculation of the radiative cooling in both collisionally and photo ionized plasmas. We are extending the spectral simulation code Cloudy so that as much as possible of the underlying atomic data is taken from external databases, some created by others, some developed by the Cloudy team. This paper focuses on recent changes in the treatment of many stages of… ▽ More We discuss recent improvements in the calculation of the radiative cooling in both collisionally and photo ionized plasmas. We are extending the spectral simulation code Cloudy so that as much as possible of the underlying atomic data is taken from external databases, some created by others, some developed by the Cloudy team. This paper focuses on recent changes in the treatment of many stages of ionization of iron, and discusses its extensions to other elements. The H-like and He-like ions are treated in the iso-electronic approach described previously. Fe II is a special case treated with a large model atom. Here we focus on Fe III through Fe XXIV, ions which are important contributors to the radiative cooling of hot, 1e5 to 1e7 K, plasmas and for X-ray spectroscopy. We use the Chianti atomic database to greatly expand the number of transitions in the cooling function. Chianti only includes lines that have atomic data computed by sophisticated methods. This limits the line list to lower excitation, longer wavelength, transitions. We had previously included lines from the Opacity Project database, which tends to include higher energy, shorter wavelength, transitions. These were combined with various forms of the g-bar approximation, a highly approximate method of estimating collision rates. For several iron ions the two databases are almost entirely complementary. We adopt a hybrid approach in which we use Chianti where possible, supplemented by lines from the Opacity Project for shorter wavelength transitions. The total cooling including the lightest thirty elements differs significantly from some previous calculations. △ Less

Submitted 5 December, 2012; originally announced December 2012.

Pumping up the [N I] nebular lines

Authors: G. J. Ferland, W. J. Henney, C. R. O'Dell, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams

Abstract: The optical [N I] doublet near 5200 Å is anomalously strong in a variety of emission-line objects. We compute a detailed photoionization model and use it to show that pumping by far-ultraviolet (FUV) stellar radiation previously posited as a general explanation applies to the Orion Nebula (M42) and its companion M43; but, it is unlikely to explain planetary nebulae and supernova remnants. Our mode… ▽ More The optical [N I] doublet near 5200 Å is anomalously strong in a variety of emission-line objects. We compute a detailed photoionization model and use it to show that pumping by far-ultraviolet (FUV) stellar radiation previously posited as a general explanation applies to the Orion Nebula (M42) and its companion M43; but, it is unlikely to explain planetary nebulae and supernova remnants. Our models establish that the observed nearly constant equivalent width of [N I] with respect to the dust-scattered stellar continuum depends primarily on three factors: the FUV to visual-band flux ratio of the stellar population; the optical properties of the dust; and the line broadening where the pumping occurs. In contrast, the intensity ratio [N I]/Hβ depends primarily on the FUV to extreme-ultraviolet ratio, which varies strongly with the spectral type of the exciting star. This is consistent with the observed difference of a factor of five between M42 and M43, which are excited by an O7 and B0.5 star respectively. We derive a non-thermal broadening of order 5 km/s for the [N I] pumping zone and show that the broadening mechanism must be different from the large-scale turbulent motions that have been suggested to explain the line-widths in this H II region. A mechanism is required that operates at scales of a few astronomical units, which may be driven by thermal instabilities of neutral gas in the range 1000 to 3000 K. In an appendix, we describe how collisional and radiative processes are treated in the detailed model N I atom now included in the Cloudy plasma code. △ Less

Submitted 31 August, 2012; originally announced September 2012.

Comments: ApJ in press. 8 pages of main paper plus 11 pages of appendices, with 13 figures and 12 tables

Hydrogen Two-Photon Continuum Emission from the Horseshoe Filament in NGC 1275

Authors: R. M. Johnstone, R. E. A. Canning, A. C. Fabian, G. J. Ferland, M. Lykins, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams

Abstract: Far ultraviolet emission has been detected from a knot of Halpha emission in the Horseshoe filament, far out in the NGC 1275 nebula. The flux detected relative to the brightness of the Halpha line in the same spatial region is very close to that expected from Hydrogen two-photon continuum emission in the particle heating model of Ferland et al. (2009) if reddening internal to the filaments is take… ▽ More Far ultraviolet emission has been detected from a knot of Halpha emission in the Horseshoe filament, far out in the NGC 1275 nebula. The flux detected relative to the brightness of the Halpha line in the same spatial region is very close to that expected from Hydrogen two-photon continuum emission in the particle heating model of Ferland et al. (2009) if reddening internal to the filaments is taken into account. We find no need to invoke other sources of far ultraviolet emission such as hot stars or emission lines from CIV in intermediate temperature gas to explain these data. △ Less

Submitted 27 June, 2012; originally announced June 2012.

Comments: 9 pages, 8 figures. Accepted for publication in MNRAS

The energy source of the filaments around the giant galaxy NGC1275

Authors: A. C. Fabian, J. S. Sanders, R. J. R. Williams, A. Lazarian, G. J. Ferland, R. M. Johnstone

Abstract: The brightest galaxy in the nearby Perseus cluster, NGC1275, is surrounded by a network of filaments. These were first observed through their Halpha emission but are now known to have a large molecular component with a total mass approaching 10^11Msun of gas. The filaments are embedded in hot intracluster gas and stretch over 80 kpc. They have an unusual low excitation spectrum which is well model… ▽ More The brightest galaxy in the nearby Perseus cluster, NGC1275, is surrounded by a network of filaments. These were first observed through their Halpha emission but are now known to have a large molecular component with a total mass approaching 10^11Msun of gas. The filaments are embedded in hot intracluster gas and stretch over 80 kpc. They have an unusual low excitation spectrum which is well modelled by collisional heating and ionization by secondary electrons. Here we note that the surface radiative flux from the outer filaments is close to the energy flux impacting on them from particles in the hot gas. We propose that the secondary electrons within the cold filaments, which excite the observed submillimetre through UV emission, are due to the hot surrounding gas efficiently penetrating the cold gas through reconnection diffusion. Some of the soft X-ray emission seen from the filaments is then due to charge exchange, although this is insufficient to account for all the observed X-ray flux. The filaments are complex with multiphase gas. Interpenetration of hot and cold gas leads to the filaments growing in mass, at a rate of up to 100Msunpyr. The lack of soft X-ray cooling emission in cool core clusters is then due to the non-radiative cooling of hot gas on mixing with cold gas around and within the central galaxy. △ Less

Submitted 9 May, 2011; originally announced May 2011.

Comments: 7 pages, 4 figures, MNRAS in press

Implications of infalling Fe II - emitting clouds in active galactic nuclei: anisotropic properties

Authors: Gary J. Ferland, Chen Hu, Jian-Min Wang, Jack A. Baldwin, Ryan L. Porter, Peter A. M. van Hoof, R. J. R. Williams

Abstract: We investigate consequences of the discovery that Fe II emission in quasars, one of the spectroscopic signatures of "Eigenvector 1", may originate in infalling clouds. Eigenvector 1 correlates with the Eddington ratio L/L_Edd so that Fe II/Hbeta increases as L/L_Edd increases. We show that the "force multiplier", the ratio of gas opacity to electron scattering opacity, is ~ 10^3 - 10^4 in Fe II-… ▽ More We investigate consequences of the discovery that Fe II emission in quasars, one of the spectroscopic signatures of "Eigenvector 1", may originate in infalling clouds. Eigenvector 1 correlates with the Eddington ratio L/L_Edd so that Fe II/Hbeta increases as L/L_Edd increases. We show that the "force multiplier", the ratio of gas opacity to electron scattering opacity, is ~ 10^3 - 10^4 in Fe II-emitting gas. Such gas would be accelerated away from the central object if the radiation force is able to act on the entire cloud. As had previously been deduced, infall requires that the clouds have large column densities so that a substantial amount of shielded gas is present. The critical column density required for infall to occur depends on L/L_Edd, establishing a link between Eigenvector 1 and the Fe II/Hbeta ratio. We see predominantly the shielded face of the infalling clouds rather than the symmetric distribution of emitters that has been assumed. The Fe II spectrum emitted by the shielded face is in good agreement with observations thus solving several long-standing mysteries in quasar emission lines. △ Less

Submitted 19 November, 2009; v1 submitted 6 November, 2009; originally announced November 2009.

Comments: 5 figures, Accepted for publication in ApJ (Letters)

Journal ref: Astrophys.J.707:L82-L86,2009

Diffuse continuum transfer in H II regions

Authors: R. J. R. Williams, W. J. Henney

Abstract: We compare the accuracy of various methods for determining the transfer of the diffuse Lyman continuum in HII regions, by comparing them with a high-resolution discrete-ordinate integration. We use these results to suggest how, in multidimensional dynamical simulations, the diffuse field may be treated with acceptable accuracy without requiring detailed transport solutions. The angular distribut… ▽ More We compare the accuracy of various methods for determining the transfer of the diffuse Lyman continuum in HII regions, by comparing them with a high-resolution discrete-ordinate integration. We use these results to suggest how, in multidimensional dynamical simulations, the diffuse field may be treated with acceptable accuracy without requiring detailed transport solutions. The angular distribution of the diffuse field derived from the numerical integration provides insight into the likely effects of the diffuse field for various material distributions. △ Less

Submitted 5 October, 2009; originally announced October 2009.

Comments: 12 pages, 7 figures, to be published in MNRAS

Laboratory Experiments, Numerical Simulations, and Astronomical Observations of Deflected Supersonic Jets: Application to HH 110

Authors: P. Hartigan, J. M. Foster, B. H. Wilde, R. F. Coker, P. A. Rosen, J. F. Hansen, B. E. Blue, R. J. R. Williams, R. Carver, A. Frank

Abstract: Collimated supersonic flows in laboratory experiments behave in a similar manner to astrophysical jets provided that radiation, viscosity, and thermal conductivity are unimportant in the laboratory jets, and that the experimental and astrophysical jets share similar dimensionless parameters such as the Mach number and the ratio of the density between the jet and the ambient medium. Laboratory je… ▽ More Collimated supersonic flows in laboratory experiments behave in a similar manner to astrophysical jets provided that radiation, viscosity, and thermal conductivity are unimportant in the laboratory jets, and that the experimental and astrophysical jets share similar dimensionless parameters such as the Mach number and the ratio of the density between the jet and the ambient medium. Laboratory jets can be studied for a variety of initial conditions, arbitrary viewing angles, and different times, attributes especially helpful for interpreting astronomical images where the viewing angle and initial conditions are fixed and the time domain is limited. Experiments are also a powerful way to test numerical fluid codes in a parameter range where the codes must perform well. In this paper we combine images from a series of laboratory experiments of deflected supersonic jets with numerical simulations and new spectral observations of an astrophysical example, the young stellar jet HH 110. The experiments provide key insights into how deflected jets evolve in 3-D, particularly within working surfaces where multiple subsonic shells and filaments form, and along the interface where shocked jet material penetrates into and destroys the obstacle along its path. The experiments also underscore the importance of the viewing angle in determining what an observer will see. The simulations match the experiments so well that we can use the simulated velocity maps to compare the dynamics in the experiment with those implied by the astronomical spectra. The experiments support a model where the observed shock structures in HH 110 form as a result of a pulsed driving source rather than from weak shocks that may arise in the supersonic shear layer between the Mach disk and bow shock of the jet's working surface. △ Less

Submitted 1 October, 2009; originally announced October 2009.

Comments: Full resolution figures available at http://sparky.rice.edu/~hartigan/pub.html To appear in ApJ

Journal ref: Astrophys.J.705:1073-1094,2009

Collisional heating as the origin of filament emission in galaxy clusters

Authors: G. J. Ferland, A. C. Fabian, N. A. Hatch, R. M. Johnstone, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams

Abstract: [Abridged] Photoionization, whether by starlight or other sources, has difficulty in accounting for the observed spectra of the optical filaments that often surround central galaxies in large clusters. Our first paper examined whether heating by energetic particles or dissipative MHD wave can account for the observations. Here we include atomic and low-ionization regions. The model of the hydrog… ▽ More [Abridged] Photoionization, whether by starlight or other sources, has difficulty in accounting for the observed spectra of the optical filaments that often surround central galaxies in large clusters. Our first paper examined whether heating by energetic particles or dissipative MHD wave can account for the observations. Here we include atomic and low-ionization regions. The model of the hydrogen atom, along with all elements of the H-like iso-electronic sequence, is now fully nl-resolved. We show how the predicted HI spectrum differs from the pure recombination case. The second update is to the rates for H^0 - H2 inelastic collisions. We now use the values computed by Wrathmall et al. The rates are often much larger. We calculate the chemistry, ionization, temperature, gas pressure, and emission-line spectrum for a wide range of gas densities and collisional heating rates. We assume that the filaments are magnetically confined and free to move along field lines so that the gas pressure is equal to that of the surrounding hot gas. A mix of clouds, some being dense and cold and others hot and tenuous, can exist. The observed spectrum will be the integrated emission from clouds with different densities and temperatures but the same pressure. We assume that the gas filling factor is given by a power law in density. The power-law index is set by matching the observed intensities of IR H2 lines relative to optical HI lines. We conclude that the filaments are heated by ionizing particles, either conducted in from surrounding regions or produced in situ by processes related to MHD waves. △ Less

Submitted 29 October, 2008; originally announced October 2008.

Comments: 28 pages, 21 figures, accepted to MNRAS

The origin of molecular hydrogen emission in cooling-flow filaments

Authors: G. J. Ferland, A. C. Fabian, N. A. Hatch, R. M. Johnstone, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams

Abstract: The optical filaments found in many cooling flows in galaxy clusters consist of low density ($\sim 10^3 \pcc$) cool ($\sim 10^3$ K) gas surrounded by significant amounts of cosmic-ray and magnetic-field energy. Their spectra show anomalously strong low-ionization and molecular emission lines when compared with galactic molecular clouds exposed to ionizing radiation such as the Orion complex. Pre… ▽ More The optical filaments found in many cooling flows in galaxy clusters consist of low density ($\sim 10^3 \pcc$) cool ($\sim 10^3$ K) gas surrounded by significant amounts of cosmic-ray and magnetic-field energy. Their spectra show anomalously strong low-ionization and molecular emission lines when compared with galactic molecular clouds exposed to ionizing radiation such as the Orion complex. Previous studies have shown that the spectra cannot be produced by O-star photoionization. Here we calculate the physical conditions in dusty gas that is well shielded from external sources of ionizing photons and is energized either by cosmic rays or dissipative MHD waves. Strong molecular hydrogen lines, with relative intensities similar to those observed, are produced. Selection effects introduced by the microphysics produce a correlation between the \htwo line upper level energy and the population temperature. These selection effects allow a purely collisional gas to produce \htwo emission that masquerades as starlight-pumped \htwo but with intensities that are far stronger. This physics may find application to any environment where a broad range of gas densities or heating rates occur. △ Less

Submitted 18 February, 2008; originally announced February 2008.

Comments: 5 pages, 4 figures, accepted to MNRAS Letters

Merged ionization/dissociation fronts in planetary nebulae

Authors: William J. Henney, R. J. R. Williams, Gary J. Ferland, Gargi Shaw, C. R. O'Dell

Abstract: The hydrogen ionization and dissociation front around an ultraviolet radiation source should merge when the ratio of ionizing photon flux to gas density is sufficiently low and the spectrum is sufficiently hard. This regime is particularly relevant to the molecular knots that are commonly found in evolved planetary nebulae, such as the Helix Nebula, where traditional models of photodissociation… ▽ More The hydrogen ionization and dissociation front around an ultraviolet radiation source should merge when the ratio of ionizing photon flux to gas density is sufficiently low and the spectrum is sufficiently hard. This regime is particularly relevant to the molecular knots that are commonly found in evolved planetary nebulae, such as the Helix Nebula, where traditional models of photodissociation regions have proved unable to explain the high observed luminosity in H_2 lines. In this paper we present results for the structure and steady-state dynamics of such advection-dominated merged fronts, calculated using the Cloudy plasma/molecular physics code. We find that the principal destruction processes for H_2 are photoionization by extreme ultraviolet radiation and charge exchange reactions with protons, both of which form H_2^+, which rapidly combines with free electrons to undergo dissociative recombination. Advection moves the dissociation front to lower column densities than in the static case, which vastly increases the heating in the partially molecular gas due to photoionization of He^0, H_2, and H^0. This causes a significant fraction of the incident bolometric flux to be re-radiated as thermally excited infrared H_2 lines, with the lower excitation pure rotational lines arising in 1000 K gas and higher excitation H_2 lines arising in 2000 K gas, as is required to explain the H_2 spectrum of the Helix cometary knots. △ Less

Submitted 27 November, 2007; v1 submitted 27 November, 2007; originally announced November 2007.

Comments: 4 pages, accepted by ApJL, scheduled December 20 issue

Self-Consistent Dynamic Models of Steady Ionization Fronts I. Weak-D and Weak-R Fronts

Authors: William J. Henney, S. Jane Arthur, Robin J. R. Williams, Gary J. Ferland

Abstract: We present a method for including steady-state gas flows in the plasma physics code Cloudy, which was previously restricted to modeling static configurations. The numerical algorithms are described in detail, together with an example application to plane-parallel ionization-bounded HII regions. As well as providing the foundation for future applications to more complex flows, we find the followi… ▽ More We present a method for including steady-state gas flows in the plasma physics code Cloudy, which was previously restricted to modeling static configurations. The numerical algorithms are described in detail, together with an example application to plane-parallel ionization-bounded HII regions. As well as providing the foundation for future applications to more complex flows, we find the following specific results regarding the effect of advection upon ionization fronts in HII regions: 1. Significant direct effects of advection on the global emission properties occur only when the ionization parameter is lower than is typical for HII regions. 2. The overheating of partially ionized gas in the front is not large, even for supersonic (R-type) fronts. 3. The most significant morphological signature of advective fronts is an electron density spike that occurs at the ionization front. Observational evidence for such a spike is found in images of the Orion bar. 4. Plane-parallel, weak-D fronts are found to show at best a shallow correlation between mean velocity and ionization potential for optical emission lines even when the flow velocity closely approaches the ionized sound speed. △ Less

Submitted 4 January, 2005; originally announced January 2005.

Comments: ApJ, main journal, in press 1 March 2005, v621 (22 pages, 3 appendices, 21 figures) PDF version with high-res figures available from http://www.astrosmo.unam.mx/~w.henney/astro-ph/0501034/

Journal ref: Astrophys.J. 621 (2005) 328-347

Non-spherical evolution of the line-driven wind instability

Authors: E. L. Gomez, R. J. R. Williams

Abstract: In this paper, we study the structure and stability of line driven winds using numerical hydrodynamic simulations. We calculate the radiation force from an explicit non-local solution of the radiation transfer equation, rather than a Sobolev approximation, without restricting the flow to one-dimensional symmetry. We find that the solutions which result have complex and highly variable structures… ▽ More In this paper, we study the structure and stability of line driven winds using numerical hydrodynamic simulations. We calculate the radiation force from an explicit non-local solution of the radiation transfer equation, rather than a Sobolev approximation, without restricting the flow to one-dimensional symmetry. We find that the solutions which result have complex and highly variable structures, including dense condensations which we compare to observed variable absorption features in the spectra of early type stars. △ Less

Submitted 30 May, 2003; originally announced May 2003.

Comments: Accepted by MNRAS, 17 pages, 15 postscript figures, 4 large figures appear as jpeg files

Journal ref: Mon.Not.Roy.Astron.Soc. 344 (2003) 725

Instabilities in two-fluid magnetized media with inter-component drift

Authors: P. V. Tytarenko, R. J. R. Williams, S. A. E. G. Falle

Abstract: We analyse the stability of a magnetized medium consisting of a neutral fluid and a fluid of charged particles, coupled to each other through a drag force and exposed to differential body forces (for example, as the result of radiation forces on one phase). We consider a uniform equilibrium and simple model input physics, but do not arbitrarily restrict the relative orientations of the magnetic… ▽ More We analyse the stability of a magnetized medium consisting of a neutral fluid and a fluid of charged particles, coupled to each other through a drag force and exposed to differential body forces (for example, as the result of radiation forces on one phase). We consider a uniform equilibrium and simple model input physics, but do not arbitrarily restrict the relative orientations of the magnetic field, slip velocity and wave vector of the disturbance. We find several instabilities and classify these in terms of wave resonances. We briefly apply our results to the structure of SiO maser regions appearing in the winds from late-type stars. △ Less

Submitted 31 July, 2002; originally announced July 2002.

Comments: 20 pages, 5 figures, MNRAS in press

Journal ref: Mon.Not.Roy.Astron.Soc. 337 (2002) 117

Resolved shocks in clumpy media

Authors: R. J. R. Williams, J. E. Dyson

Abstract: We study the structure of shocks in clumpy media, using a multifluid formalism. As expected, shocks broaden as they weaken: for sufficiently weak shocks, no viscous subshock appears in the structure. This has significant implications for the survival of dense clouds in regions overrun by shocks in a wide range of astrophysical circumstances, from planetary nebulae to the nuclei of starburst gala… ▽ More We study the structure of shocks in clumpy media, using a multifluid formalism. As expected, shocks broaden as they weaken: for sufficiently weak shocks, no viscous subshock appears in the structure. This has significant implications for the survival of dense clouds in regions overrun by shocks in a wide range of astrophysical circumstances, from planetary nebulae to the nuclei of starburst galaxies. △ Less

Submitted 20 November, 2001; originally announced November 2001.

Comments: 9 pages, 4 figures, Latex, MN macros, in press with MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 333 (2002) 1

Hydrodynamics of photoionized columns in the Eagle Nebula, M 16

Authors: R. J. R. Williams, D. Ward-Thompson, A. P. Whitworth

Abstract: We present hydrodynamical simulations of the formation, structure and evolution of photoionized columns, with parameters based on those observed in the Eagle Nebula. On the basis of these simulations we argue that there is no unequivocal evidence that the dense neutral clumps at heads of the columns were cores in the pre-existing molecular cloud. In our simulations, a variety of initial conditio… ▽ More We present hydrodynamical simulations of the formation, structure and evolution of photoionized columns, with parameters based on those observed in the Eagle Nebula. On the basis of these simulations we argue that there is no unequivocal evidence that the dense neutral clumps at heads of the columns were cores in the pre-existing molecular cloud. In our simulations, a variety of initial conditions leads to the formation and maintenance of near-equilibrium columns. Therefore, it is likely that narrow columns will often occur in regions with large-scale inhomogeneities, but that observations of such columns can tell us little about the processes by which they formed. The manner in which the columns in our simulations develop suggests that their evolution may result in extended sequences of radiation-induced star formation. △ Less

Submitted 16 July, 2001; originally announced July 2001.

Comments: 12 pages, 9 figures, Latex, MN macros, in press with MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 327 (2001) 788

Magnetic ionization fronts II: Jump conditions for oblique magnetization

Authors: R. J. R. Williams, J. E. Dyson, T. W. Hartquist

Abstract: We present the jump conditions for ionization fronts with oblique magnetic fields. The standard nomenclature of R- and D-type fronts can still be applied, but in the case of oblique magnetization there are fronts of each type about each of the fast- and slow-mode speeds. As an ionization front slows, it will drive first a fast- and then a slow-mode shock into the surrounding medium. Even for rat… ▽ More We present the jump conditions for ionization fronts with oblique magnetic fields. The standard nomenclature of R- and D-type fronts can still be applied, but in the case of oblique magnetization there are fronts of each type about each of the fast- and slow-mode speeds. As an ionization front slows, it will drive first a fast- and then a slow-mode shock into the surrounding medium. Even for rather weak upstream magnetic fields, the effect of magnetization on ionization front evolution can be important. [Includes numerical MHD models and an application to observations of S106.] △ Less

Submitted 7 January, 2000; originally announced January 2000.

Comments: 9 pages, 10 figures, Latex, to be published in MNRAS

Symbiotic starburst-black hole AGN -- I. Isothermal hydrodynamics of the mass-loaded ISM

Authors: R. J. R. Williams, A. C. Baker, J. J. Perry

Abstract: Compelling evidence associates the nuclei of active galaxies and massive starbursts. The symbiosis between a compact nuclear starburst stellar cluster and a massive black hole can self-consistently explain the properties of active nuclei. The young stellar cluster has a profound effect on the most important observable properties of active galaxies through its gravity, and by mass injection throu… ▽ More Compelling evidence associates the nuclei of active galaxies and massive starbursts. The symbiosis between a compact nuclear starburst stellar cluster and a massive black hole can self-consistently explain the properties of active nuclei. The young stellar cluster has a profound effect on the most important observable properties of active galaxies through its gravity, and by mass injection through stellar winds, supernovae and stellar collisions. Mass injection generates a nuclear ISM which flows under gravitational and radiative forces until it leaves the nucleus or is accreted onto the black hole or accretion disc. The radiative force exerted by the black hole--accretion disc radiation field is not spherically symmetric. This results in complex flows in which regions of inflow can coexist with high Mach number outflowing winds and hydrodynamic jets. We present two-dimensional hydrodynamic models of such nISM flows, which are highly complex and time variable. Shocked shells, jets and explosive bubbles are produced, with bipolar winds driving out from the nucleus. Our results graphically illustrate why broad emission line studies have consistently failed to identify any simple, global flow geometry. The real structure of the flows is _inevitably_ yet more complex. △ Less

Submitted 14 July, 1999; originally announced July 1999.

Comments: 51 pages, 85 postscript figures, Latex, using MNRAS macros, to be published in MNRAS. Postscript will full resolution pictures and mpeg simulations available via http://ast.leeds.ac.uk/~rjrw/agn.html

Jet-cloud interations and the brightening of the narrow line region in Seyfert galaxies

Authors: W. Steffen, J. L. Gomez, A. C. Raga, R. J. R. Williams

Abstract: We study the kinematical and brightness evolution of emission line clouds in the narrow line region (NLR) of Seyfert galaxies during the passage of a jet. We derive a critical density above which a cloud remains radiative after compression by the jet cocoon. The critical density depends mainly on the cocoon pressure. Super-critical clouds increase in emission line brightness, while sub-critical… ▽ More We study the kinematical and brightness evolution of emission line clouds in the narrow line region (NLR) of Seyfert galaxies during the passage of a jet. We derive a critical density above which a cloud remains radiative after compression by the jet cocoon. The critical density depends mainly on the cocoon pressure. Super-critical clouds increase in emission line brightness, while sub-critical clouds generally are highly overheated reducing their luminosity below that of the inter-cloud medium. Due to the pressure stratification in the bow-shock of the jet, a cylindrical structure of nested shells develops around the jet. The most compact and brightest compressed clouds surround the cloud-free channel of the radio jet. To support our analytical model we present a numerical simulation of a supersonic jet propagating into a clumpy NLR. The position-velocity diagram of the simulated H_alpha emission shows total line widths of the order of 500 km/s with large-scale variations in the radial velocities of the clouds due to the stratified pressure in the bow-shock region of the jet. Most of the luminosity is concentrated in a few dense clouds surrounding the jet. These morphological and kinematic signatures are all found in the well observed NLR of NGC1068 and other Seyfert galaxies. △ Less

Submitted 16 October, 1997; originally announced October 1997.

Comments: 11 pages, 3 figures, accepted for publication in The Astrophysical Journal Letters

Jet cocoons and the formation of Narrow Line Clouds in Seyfert galaxies

Authors: W. Steffen, J. L. Gomez, R. J. R. Williams, A. C. Raga, A. Pedlar

Abstract: We present non-adiabatic hydrodynamic simulations of a supersonic light jet propagating into a fully ionized medium of uniform density on a scale representative of the narrow line region (NLR) in Seyfert galaxies with associated radio jets. In this regime the cooling distance of the swept up gas in the bowshock of the jet is of the same order as the transverse extent of the jet bowshock, as oppo… ▽ More We present non-adiabatic hydrodynamic simulations of a supersonic light jet propagating into a fully ionized medium of uniform density on a scale representative of the narrow line region (NLR) in Seyfert galaxies with associated radio jets. In this regime the cooling distance of the swept up gas in the bowshock of the jet is of the same order as the transverse extent of the jet bowshock, as opposed to the more extreme regimes found for more powerful adiabatic large scale jets or the slow galactic jets which have been simulated previously. We calculate the emissivity for the H_alpha line and radio synchrotron emission. We find that the structure of the line emitting cold envelope of the jet cocoon is strongly dependent on the non-stationary dynamics of the jet head as it propagates through the ambient medium. We observe the formation of cloud-like high density regions which we associate with NLR clouds and filaments. We find that some of these clouds might be partially neutral and represent sites of jet induced star formation. The calculated H_alpha flux and the spectral line width are consistent with NLR observations. The simulation of the radio-optical emission with radiative cooling confirms the basic result of the geometric bowshock model developed by Taylor et al. (1989) that the start of noticeable optical line emission can be significantly offset from the hotspot of the radio emission. However, the time-dependent nature of the jet dynamics implies significant differences from their geometric bowshock model. △ Less

Submitted 5 December, 1996; originally announced December 1996.

Comments: 11 pages, 5 figures, accepted for publication in MNRAS

Highly supersonic motions within the outer features of the Eta Carinae nebulosity

Authors: J. Meaburn, P. Boumis, J. R. Walsh, W. Steffen, A. J. Holloway, R. J. R Williams, M. Bryce

Abstract: Spatially resolved H$α$ and [N{\sc ii}] line profiles have been obtained over striking features in the outer regions of the Eta Carinae nebula. The highly irregular outer `shell' of low ionization, [N{\sc ii}] bright, knots has been shown to exhibit radial velocities between -1200 and +300 km/s with respect to the systemic radial velocity, over its perimeter. Furthermore, several knots have been… ▽ More Spatially resolved H$α$ and [N{\sc ii}] line profiles have been obtained over striking features in the outer regions of the Eta Carinae nebula. The highly irregular outer `shell' of low ionization, [N{\sc ii}] bright, knots has been shown to exhibit radial velocities between -1200 and +300 km/s with respect to the systemic radial velocity, over its perimeter. Furthermore, several knots have been found which appear to emit only in the H$α$ line with radial velocities up to -1450 km/s. However, the most intriguing features are a narrow `spike' that projects through this outer shell and a faint `arc' of emission that extends well beyond it. The `spike', which exhibits a change of radial velocity along its length, could be a narrow collimated jet with an outflow velocity of > 1000 km/s. In one interpretation the `arc' is modelled by a conical outflow and mechanisms for generating [N{\sc ii}] emission from highly supersonic gas are also discussed. △ Less

Submitted 14 June, 1996; originally announced June 1996.

Comments: MNRAS style file, 6 pages. 7 eps figures are not included but are available individually at http://axp2.ast.man.ac.uk:8000/Preprints.html

Clumpy ultracompact HII regions -- III. Cometary morphologies around stationary stars

Authors: R. J. R. Williams, J. E. Dyson, M. P. Redman

Abstract: Cometary ultracompact HII regions have been modelled as the interaction of the hypersonic wind from a moving star with the molecular cloud which surrounds the star. We here show that a similar morphology can ensue even if the star is stationary with respect to the cloud material. We assume that the HII region is within a stellar wind bubble which is strongly mass loaded: the cometary shape resul… ▽ More Cometary ultracompact HII regions have been modelled as the interaction of the hypersonic wind from a moving star with the molecular cloud which surrounds the star. We here show that a similar morphology can ensue even if the star is stationary with respect to the cloud material. We assume that the HII region is within a stellar wind bubble which is strongly mass loaded: the cometary shape results from a gradient in the distribution of mass loading sources. This model circumvents problems associated with the necessarily high spatial velocities of stars in the moving star models. △ Less

Submitted 1 March, 1996; originally announced March 1996.

Comments: 6 pages with MN LaTeX macros, includes 5 postscript figures. Also available at http://axp2.ast.man.ac.uk:8000/Preprints.html

Clumpy ultracompact HII regions -- II. Cores, spheres and shells from subsonic flows

Authors: M. P. Redman, R. J. R. Williams, J. E. Dyson

Abstract: We have modelled ultracompact HII regions (UCHIIR) in terms of steady subsonic ionized flows in a clumpy medium. Mass loss from neutral clumps allows the regions to be long-lived. We examine the form of global flows for different dependences of the volume mass injection rate, $\dot{q}$, on radius and Mach number, and describe the solutions in detail. We find that three observed UCHIIR morphologi… ▽ More We have modelled ultracompact HII regions (UCHIIR) in terms of steady subsonic ionized flows in a clumpy medium. Mass loss from neutral clumps allows the regions to be long-lived. We examine the form of global flows for different dependences of the volume mass injection rate, $\dot{q}$, on radius and Mach number, and describe the solutions in detail. We find that three observed UCHIIR morphologies are reproduced with these models. Mach number independent flows that include a radial variation can give centre- brightened core--halo morphologies. Mach number dependent flows reproduce naturally the uniform UCHIIR morphology. In a hybrid model, including subsonic and supersonic flows, we allow a supersonic wind to shock in the ionized region. The ionized subsonic gas has a high density and so dominates the emission. The shell produced has a velocity structure very different from that of fully supersonic models. Several morphologies of spherical UCHIIR can be understood in terms of these various models; however, kinematic data are crucial as a discriminant between them. △ Less

Submitted 1 March, 1996; originally announced March 1996.

Comments: 6 pages with MN LaTeX macros, includes 12 postscript figures. Also available at http://axp2.ast.man.ac.uk:8000/Preprints.html

Breaking the sound barrier in recombination fronts

Authors: R. J. R. Williams, J. E. Dyson

Abstract: We exploit a generic instability in the integration of steady, sonic, near-isothermal flows to find the complete transition diagram for recombination fronts (for a model system of equations). The instability requires the integration of the flow equations for speeds between the isothermal and adiabatic sound speeds to be performed with particular care. As a result of this, the previous work of Ne… ▽ More We exploit a generic instability in the integration of steady, sonic, near-isothermal flows to find the complete transition diagram for recombination fronts (for a model system of equations). The instability requires the integration of the flow equations for speeds between the isothermal and adiabatic sound speeds to be performed with particular care. As a result of this, the previous work of Newman & Axford on the structure of recombination fronts neglected an important class of solution, that of transonic fronts; our method is readily extensible to a more complete treatment of the ionization structure. Future papers will apply these results in models of the structure of ultracompact HII regions. △ Less

Submitted 18 December, 1995; originally announced December 1995.

Comments: Accepted for publication in MNRAS. LaTeX, 7 pages, 6 Postscript figures. Also available at http://axp2.ast.man.ac.uk:8000/Preprints.html

Clumpy Ultracompact HII Regions I: Fully Supersonic Wind-blown Models

Authors: J. E. Dyson, R. J. R. Williams, M. P. Redman

Abstract: We propose that a significant fraction of the ultracompact HII regions found in massive star-forming clouds are the result of the interaction of the wind and ionizing radiation from a young massive star with the clumpy molecular cloud gas in its neighbourhood. Distributed mass loading in the flow allows the compact nebulae to be long-lived. In this paper, we discuss a particularly simple case, i… ▽ More We propose that a significant fraction of the ultracompact HII regions found in massive star-forming clouds are the result of the interaction of the wind and ionizing radiation from a young massive star with the clumpy molecular cloud gas in its neighbourhood. Distributed mass loading in the flow allows the compact nebulae to be long-lived. In this paper, we discuss a particularly simple case, in which the flow in the HII region is everywhere supersonic. The line profiles predicted for this model are highly characteristic, for the case of uniform mass loading. We discuss briefly other observational diagnostics of these models. △ Less

Submitted 18 July, 1995; originally announced July 1995.

Comments: To appear in Monthly Notices of the Royal Astronomical Society. 5 pages LaTeX (uses mn.sty and epsf.sty macros) + 4 PS figures. Also available via http://axp2.ast.man.ac.uk:8000/Preprints.html