sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Ferland, G. J.

Normalized to: Ferland, G.

213 article(s) in total. 636 co-authors, from 1 to 40 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 2126946
Hyper-massive Black Holes have Faint Broad and Narrow Emission Lines
Comments: 20 pages, 26 figures, to be published in MNRAS
Submitted: 2020-07-02
The EUV provides most of the ionization that creates the high equivalent width (EW) broad and narrow emission lines (BELs, NELs) of quasars. Spectra of Hypermassive Schwarzschild black holes (HMBHs, $M_{BH} \geq 10^{10} M_{\odot}$) with $\alpha$-discs, decline rapidly in the EUV suggesting much lower EWs. Model spectra for black holes of mass $10^{6}-10^{12} M_{\odot}$ and accretion rates $0.03 \leq L_{bol}/L_{edd} \leq 1.0$ were input to the CLOUDY photoionization code. BELs become $\sim$100 times weaker in EW from $M_{BH} \sim 10^8 M_{\odot}$ to $M_{BH} \sim 10^{10} M_{\odot}$. The high ionization BELs (O VI 1034 $\overset{\circ}{\mathrm {A}}$, C IV 1549 $\overset{\circ}{\mathrm {A}}$, He II 1640 $\overset{\circ}{\mathrm {A}}$) decline in EW from ($M_{BH} \geq 10^6 M_{\odot}$, reproducing the Baldwin effect, but regain EW for $M_{BH} \geq 10^{10} M_{\odot}$). The low ionization lines (MgII 2798 $\overset{\circ}{\mathrm {A}}$, H$\beta$ 4861 $\overset{\circ}{\mathrm {A}}$ and H$\alpha$ 6563 $\overset{\circ}{\mathrm {A}}$) remain weak. Lines for maximally spinning HMBHs behave similarly. Line ratio diagrams for the BELs show that high OVI/H$\beta$ and low CIV/H$\alpha$ may pick out HMBH, although OVI is often hard to observe. In NEL BPT diagrams HMBHs lie among star-forming regions, except for highly spinning, high accretion rate HMBHs. In summary, the BELs expected from HMBHs would be hard to detect using the current optical facilities. From 100 to $10^{12} M_{\odot}$, the emission lines used to detect AGN only have high EW in the $10^6 - 10^9 M_{\odot}$ window, where most AGN are found. This selection effect may be distorting reported distributions of $M_{BH}$.
[2]  [pdf] - 2124923
Fe III emission in quasars: evidence for a dense turbulent medium
Comments: 12 pages, 11 figures, plus appendices
Submitted: 2020-06-15, last modified: 2020-06-30
Recent improvements to atomic energy-level data allow, for the first time, accurate predictions to be made for the Fe III line emission strengths in the spectra of luminous, $L_\text{bol}=10^{46}-10^{48}$ erg/s, Active Galactic Nuclei. The Fe III emitting gas must be primarily photoionized, consistent with observations of line reverberation. We use CLOUDY models exploring a wide range of parameter space, together with 26,500 rest-frame ultraviolet spectra from the Sloan Digital Sky Survey, to constrain the physical conditions of the line emitting gas. The observed Fe III emission is best accounted for by dense ($n_H=10^{14}$ cm$^{-3}$) gas which is microturbulent, leading to smaller line optical depths and fluorescent excitation. Such high density gas appears to be present in the central regions of the majority of luminous quasars. Using our favoured model, we present theoretical predictions for the relative strengths of the Fe III UV34 $\lambda\lambda$1895,1914,1926 multiplet. This multiplet is blended with the Si III] $\lambda$1892 and C III] $\lambda$1909 emission lines and an accurate subtraction of UV34 is essential when using these lines to infer information about the physics of the broad line region in quasars.
[3]  [pdf] - 2112760
Space Telescope and Optical Reverberation Mapping Project. XI. Disk-wind characteristics and contributions to the very broad emission lines of NGC 5548
Comments: Accepted for publication in ApJ
Submitted: 2020-06-11
In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper showed that changes in the LOS obscurer produce the decorrelation between the absorption lines and the continuum. A second study showed that the base of the wind shields the BLR, leading to the emission-line decorrelation. In that study, we proposed the wind is normally transparent with no effect on the spectrum. Changes in the wind properties alter its shielding and affect the SED striking the BLR, producing the observed decorrelations. In this work, we investigate the impact of a translucent wind on the emission lines. We simulate the obscuration using XMM-Newton, NuSTAR, and HST observations to determine the physical characteristics of the wind. We find that a translucent wind can contribute a part of the He II and Fe K? emission. It has a modest optical depth to electron scattering, which explains the fainter far-side emission in the observed velocity delay maps. The wind produces the very broad base seen in the UV emission lines and may also be present in the Fe K? line. Our results highlight the importance of accounting for the effects of such winds in the analysis of the physics of the central engine.
[4]  [pdf] - 2105774
A molecular absorption line survey toward the AGN of Hydra-A
Comments: Accepted for publication in MNRAS
Submitted: 2020-05-20, last modified: 2020-05-28
We present Atacama Large Millimeter/submillimeter Array observations of the brightest cluster galaxy Hydra-A, a nearby ($z=0.054$) giant elliptical galaxy with powerful and extended radio jets. The observations reveal CO(1-0), CO(2-1), $^{13}$CO(2-1), CN(2-1), SiO(5-4), HCO$^{+}$(1-0), HCO$^{+}$(2-1), HCN(1-0), HCN(2-1), HNC(1-0) and H$_{2}$CO(3-2) absorption lines against the galaxy's bright and compact active galactic nucleus. These absorption features are due to at least 12 individual molecular clouds which lie close to the centre of the galaxy and have velocities of approximately $-50$ to $+10$ km/s relative to its recession velocity, where positive values correspond to inward motion. The absorption profiles are evidence of a clumpy interstellar medium within brightest cluster galaxies composed of clouds with similar column densities, velocity dispersions and excitation temperatures to those found at radii of several kpc in the Milky Way. We also show potential variation in a $\sim 10$ km/s wide section of the absorption profile over a two year timescale, most likely caused by relativistic motions in the hot spots of the continuum source which change the background illumination of the absorbing clouds.
[5]  [pdf] - 2093325
State-of-the-art AGN SEDs for Photoionization Models: BLR Predictions Confront the Observations
Comments: 6 pages, 2 figures. Accepted for publication in MNRAS
Submitted: 2020-04-24
The great power offered by photoionization models of Active Galactic Nuclei (AGN) emission-line regions has long been mitigated by the fact that very little is known about the spectral energy distribution (SED) between the Lyman limit, where intervening absorption becomes a problem, and 0.3 keV, where soft x-ray observations become possible. The emission lines themselves can, to some degree, be used to probe the SED, but only in the broadest terms. This paper employs a new generation of theoretical SEDs which are internally self-consistent, energy-conserving, and tested against observations, to infer properties of the emission-line regions. The SEDs are given as a function of the Eddington ratio, allowing emission-line correlations to be investigated on a fundamental basis. We apply the simplest possible tests, based on the foundations of photoionization theory, to investigate the implications for the geometry of the emission-line region. The SEDs become more far-ultraviolet bright as the Eddington ratio increases, so the equivalent widths of recombination lines should also become larger, an effect which we quantify. The observed lack of correlation between Eddington ratio and equivalent width shows that the cloud covering factor must decrease as Eddington ratio increases. This would be consistent with recent models proposing that the broad-line region is a failed dusty wind off the accretion disc.
[6]  [pdf] - 2058005
Space Telescope and Optical Reverberation Mapping Project. IX. Velocity-Delay Maps for Broad Emission Lines in NGC 5548
Horne, Keith; De Rosa, G.; Peterson, B. M.; Barth, A. J.; Ely, J.; Fausnaugh, M. M.; Kriss, G. A.; Pei, L.; Adams, S. M.; Anderson, M. D.; Arevalo, P.; Beatty, T G.; Bennert, V. N.; Bentz, M. C.; Bigley, A.; Bisogni, S.; Borman, G. A.; Boroson, T. A.; Bottorff, M. C.; Brandt, W. N.; Breeveld, A. A.; Brotherton, M.; Brown, J. E.; Brown, J. S.; Cackett, E. M.; Canalizo, G.; Carini, M. T.; Clubb, K. I.; Comerford, J. M.; Corsini, E. M.; Crenshaw, D. M.; Croft, S.; Croxall, K. V.; Bonta, E. Dalla; Deason, A. J.; Dehghanian, M.; De Lorenzo-Caceres, A.; Denney, K. D.; Dietrich, M.; Done, C.; Edelson, R.; Efimova, N. V.; Eracleous, M.; Evans, P. A.; Ferland, G. J.; Filippenko, A. V.; Flatland, K.; Fox, O. D.; Gardner, E.; Gates, E. L.; Gehrels, N.; Geier, S.; Gelbord, J. M.; Goad, M. R.; Gonzalez, L.; Gorjian, V.; Greene, J. E.; Grier, C. J.; Grupe, D.; Gupta, A.; Hall, P. B.; Henderson, C. B.; Hicks, S.; Holmbeck, E.; Holoien, T. W. -S.; Hutchison, T.; Im, M.; Jensen, J. J.; Johnson, C. A.; Joner, M. D.; Jones, J.; Kaastra, J.; Kaspi, S.; Kelly, P. L.; Kennea, J. A.; Kim, M.; Kim, S.; Kim, S. C.; King, A.; Klimanov, S. A.; Kochanek, C. S.; Korista, K. T.; Krongold, Y.; Lau, M. W.; Lee, J. C.; Leonard, D. C.; Li, Miao; Lira, P.; Lochhaas, C.; Ma, Zhiyuan; MacInnis, F.; Malkan, M. A.; Manne-Nicholas, E. R.; Mathur, S.; Mauerhan, J. C.; McGurk, R.; Hardy, I. M. Mc; Montuori, C.; Morelli, L.; Mosquera, A.; Mudd, D.; Mueller-Sanchez, F.; Nazarov, S. V.; Netzer, H.; Norris, R. P.; Nousek, J. A.; Nguyen, M. L.; Ochner, P.; Okhmat, D. N.; Pancoast, A.; Papadakis, I.; Parks, J. R.; Penny, M. T.; Pizzella, A.; Pogge, R. W.; Poleski, R.; Pott, J. -U.; Proga, D.; Rafter, S. E.; Rix, H. -W.; Runnoe, J.; Saylor, D. A.; Schimoia, J. S.; Schnuelle, K.; Scott, B.; Sergeev, S. G.; Shappee, B. J.; Shivvers, I.; Siegel, M.; Simonian, G. V.; Siviero, A.; Skielboe, A.; Somers, G.; Spencer, M.; Starkey, D.; Stevens, D. J.; Strauss, M. A.; Sung, H. -I.; Tayar, J.; Teems, K. G.; Tejos, N.; Treu, T.; Turner, C. S.; Uttley, P.; Van Saders, J .; Vestergaard, M.; Vican, L.; Villanueva, S.; Villforth, C.; Weiss, Y.; Woo, J. -H.; Yan, H.; Young, S.; Yuk, H.; Zakamska, N. L.; Zheng, W.; Zhu, W.; Zu, Y.
Comments: 19 pages, 9 figures, to be submitted to ApJ
Submitted: 2020-03-03
We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.
[7]  [pdf] - 2059269
Deciphering the 3-D Orion Nebula-I: Expanding Shells in the Huygens Region
Submitted: 2020-03-03
Based on imaging and spectroscopic data, we develop a 3-D model for the Huygens Region of the Orion Nebula. Theta1OriC , the hottest star in the Trapezium, is surrounded by a wind-blown Central Bubble that opens SW into the Extended Orion Nebula. Outside of this feature lies a layer of ionized gas at about 0.4 pc from Theta1OriC. Both of these features are moving rapidly away from Theta1OriC with an expansion age for the Central Bubble of only 15,000 yrs.
[8]  [pdf] - 2065455
Physical conditions in two high-redshift H$_2$-bearing GRB-DLAs, 120815A and 121024A
Comments: 10 pages, 5 figures, 6 tables. Accepted for publication in MNRAS
Submitted: 2020-03-01
The gamma-ray burst (GRB) afterglows provide an unique opportunity to study the interstellar medium (ISM) of star-forming galaxies at high-redshift. The GRB-DLAs (damped Lyman-$\alpha$ absorbers) contain a large neutral hydrogen column density, N(H I), and are observed against the GRB afterglow. A large fraction of GRB-DLAs show presence of molecular hydrogen (H$_2$) which is an indicator of star-formation. Hence it is important to study those GRB-DLAs which have H$_2$ lines to decipher and understand their physical conditions. The GRB-DLAs 121024A and 120815A, situated at redshift 2.30 and 2.36 respectively, are two such important H$_2$-bearing GRB-DLAs. Besides H$_2$, these two GRB-DLAs also show many metal lines. In this work we have carried out a detail numerical study on the H$_2$ lines, as well as on those metal lines, in GRB-DLAs 121024A and 120815A self-consistently. We use the spectral synthesis code CLOUDY for this study. This modeling helps us to determine the underlying physical conditions which give rise such lines and hence to understand these two GRB-DLAs in much more detail than any other previous investigation. We find that the hydrogen densities for these two H$_2$-bearing DLAs are $\geq$ 60 cm$^{-3}$. Moreover our study infers that the linear sizes are $\leq 17.7$ pc for these two GRB-DLAs, and the mean gas temperatures averaged over the cloud thickness, are $\leq$ 140 K.
[9]  [pdf] - 2072698
Newly Improved Ionization Corrections for the Neutral Interstellar Medium: Enabling Accurate Abundance Determinations in Star-forming Galaxies throughout the Universe
Comments: 20 pages, 8 tables, 8 Figures, Accepted for Publication in ApJ
Submitted: 2020-02-18
Studies measuring the chemical abundances of the neutral gas in star-forming galaxies (SFGs) require ionization correction factors (ICFs) to accurately measure their metal contents. In the work presented here we calculate newly improved ICFs for a sample of SFGs. These new corrections include both the contaminating ionized gas along the line of sight (ICF$_{\rm ionized}$) and unaccounted higher ionization stages in the neutral gas (ICF$_{\rm neutral}$). We make use of recently acquired spectroscopic observations taken with the Cosmic Origins Spectrograph (COS) on board Hubble to measure column densities for Fe II and Fe III. Using the Fe III/Fe II ratios as well as other physical properties (i.e. $\log$[L$_{\rm UV}$], $N$(H I), T, and $Z$) we generate ad-hoc photoionization models with CLOUDY to quantify the corrections required for each of the targets. We identify a luminosity threshold of $\log$[L$_{\rm UV}$]$\sim$ 40.75 erg s$^{-1}$ above which the ICF$_{\rm neutral}$ values for nitrogen are relatively higher (ICF$_{\rm neutral}=0.05$-0.7) than those for the rest of the elements (ICF$_{\rm neutral}\sim 0.01$). This behavior indicates that for the high UV luminosity objects, N II is found in non-negligible quantities in the neutral gas, making these ICF$_{\rm neutral}$ corrections critical for determining the true abundances in the interstellar medium. In addition, we calculate ICFs from a uniform grid of models covering a wide range of physical properties typically observed in studies of SFGs and extragalactic H II regions. We provide the community with tabulated ICF values for the neutral gas abundances measured from a variety of environments and applicable to chemical studies of the high redshift universe.
[10]  [pdf] - 2048644
Current and future development of the photoionization code Cloudy
Comments: 12 pages, 5 figures. This paper was presented as an invited lecture at the 12th Serbian Conference on Spectral Line Shapes in Astrophysics in Vrdnik, Serbia
Submitted: 2020-02-13
The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical state of the gas and predicts the spectrum that it emits. Cloudy is continually being developed to improve the treatment of the microphysical processes and the database of fundamental data that it uses. In this paper we will discuss how we are developing the code to improve our high-density predictions by implementing better collisional-radiative models for all ions. We will also briefly discuss the experimental mode in Cloudy to model gas that is not in steady-state equilibrium and present a preliminary model of recombining gas in a planetary nebula that is on the cooling track. We finish with a short discussion of how we are speeding up the code by using parallelization.
[11]  [pdf] - 2040327
2018 NASA Laboratory Astrophysics Workshop: Scientific Organizing Committee Report
Submitted: 2020-01-28
This report provides detailed findings on the critical laboratory astrophysics data needs that are required to maximize the scientific return for NASA's current and near-term planned astrophysics missions. It also provides prioritized rankings on said laboratory astrophysics data, generally by waveband. The Report is based on community input gathered at the 2018 NASA Laboratory Astrophysics Workshop (LAW) from presentations, from discussions during workshop breakout sessions, and from other solicited input deemed appropriate by the Scientific Organizing Committee (SOC) obtained prior to and after the meeting. Hence, the Report is a direct reflection of the spirit and participant make-up of LAW 2018. The Report also outlines specific opportunities and threats facing NASA's Laboratory Astrophysics Program, and articulates concrete actions by which the Agency can capitalize on the opportunities and mitigate the challenges. The Report was prepared by the SOC, with help from some invited speakers, and input and review from community members.
[12]  [pdf] - 1966978
A wind-based unification model for NGC 5548: spectral holidays, non-disk emission, and implications for changing-look quasars
Comments: Accepted for publication in ApJL
Submitted: 2019-08-20, last modified: 2019-08-22
The 180-day Space Telescope and Optical Reverberation Mapping campaign on NGC 5548 discovered an anomalous period, the broad-line region (BLR) holiday, in which the emission lines decorrelated from the continuum variations. This is important since the correlation between the continuum-flux variations and the emission-line response is the basic assumption for black hole (BH) mass determinations through reverberation mapping. During the BLR holiday, the high-ionization intrinsic absorption lines also decorrelated from the continuum as a result of variable covering factor of the line of sight (LOS) obscurer. The emission lines are not confined to the LOS, so this does not explain the BLR holiday. If the LOS obscurer is a disk wind, its streamlines must extend down to the plane of the disk and the base of the wind would lie between the BH and the BLR, forming an equatorial obscurer. This obscurer can be transparent to ionizing radiation, or can be translucent, blocking only parts of the SED, depending on its density. An emission-line holiday is produced if the wind density increases only slightly above its transparent state. Both obscurers are parts of the same wind, so they can have associated behavior in a way that explains both holidays. A very dense wind would block nearly all ionizing radiation, producing a Seyfert 2 and possibly providing a contributor to the changing-look AGN phenomenon. Disk winds are very common and we propose that the equatorial obscurers are too, but mostly in a transparent state.
[13]  [pdf] - 1945881
The first spectroscopic dust reverberation programme on active galactic nuclei: the torus in NGC 5548
Comments: 22 pages, 13 figures, accepted for publication in MNRAS
Submitted: 2019-08-05
We have recently initiated the first spectroscopic dust reverberation programme on active galactic nuclei (AGN) in the near-infrared. Spectroscopy enables measurement of dust properties, such as flux, temperature and covering factor, with higher precision than photometry. In particular, it enables measurement of both luminosity-based dust radii and dust response times. Here we report results from a one-year campaign on NGC 5548. The hot dust responds to changes in the irradiating flux with a lag time of ~70 light-days, similar to what was previously found in photometric reverberation campaigns. The mean and rms spectra are similar, implying that the same dust component dominates both the emission and the variations. The dust lag time is consistent with the luminosity-based dust radius only if we assume a wavelength-independent dust emissivity-law, i.e. a blackbody, which is appropriate for grains of large sizes (of a few microns). For such grains the dust temperature is ~1450 K. Therefore, silicate grains have most likely evaporated and carbon is the main chemical component. But the hot dust is not close to its sublimation temperature, contrary to popular belief. This is further supported by our observation of temperature variations largely consistent with a heating/cooling process. Therefore, the inner dust-free region is enlarged and the dusty torus rather a "dusty wall", whose inner radius is expected to be luminosity-invariant. The dust-destruction mechanism that enlarges the dust-free region seems to partly affect also the dusty region. We observe a cyclical decrease in dust mass with implied dust reformation times of ~5-6 months.
[14]  [pdf] - 1938579
Constraining cold accretion onto supermassive black holes: molecular gas in the cores of eight brightest cluster galaxies revealed by joint CO and CN absorption
Comments: Accepted for publication in MNRAS
Submitted: 2019-07-31
To advance our understanding of the fuelling and feedback processes which power the Universe's most massive black holes, we require a significant increase in our knowledge of the molecular gas which exists in their immediate surroundings. However, the behaviour of this gas is poorly understood due to the difficulties associated with observing it directly. We report on a survey of 18 brightest cluster galaxies lying in cool cores, from which we detect molecular gas in the core regions of eight via carbon monoxide (CO), cyanide (CN) and silicon monoxide (SiO) absorption lines. These absorption lines are produced by cold molecular gas clouds which lie along the line of sight to the bright continuum sources at the galaxy centres. As such, they can be used to determine many properties of the molecular gas which may go on to fuel supermassive black hole accretion and AGN feedback mechanisms. The absorption regions detected have velocities ranging from -45 to 283 km s$^{-1}$ relative to the systemic velocity of the galaxy, and have a bias for motion towards the host supermassive black hole. We find that the CN N = 0 - 1 absorption lines are typically 10 times stronger than those of CO J = 0 - 1. This is due to the higher electric dipole moment of the CN molecule, which enhances its absorption strength. In terms of molecular number density CO remains the more prevalent molecule with a ratio of CO/CN $\sim 10$, similar to that of nearby galaxies. Comparison of CO, CN and HI observations for these systems shows many different combinations of these absorption lines being detected.
[15]  [pdf] - 1953637
Space Telescope and Optical Reverberation Mapping Project. VIII. Time Variability of Emission and Absorption in NGC 5548 Based on Modeling the Ultraviolet Spectrum
Kriss, G. A.; De Rosa, G.; Ely, J.; Peterson, B. M.; Kaastra, J.; Mehdipour, M.; Ferland, G. J.; Dehghanian, M.; Mathur, S.; Edelson, R.; Korista, K. T.; Arav, N.; Barth, A. J.; Bentz, M. C.; Brandt, W. N.; Crenshaw, D. M.; Bontà, E. Dalla; Denney, K. D.; Done, C.; Eracleous, M.; Fausnaugh, M. M.; Gardner, E.; Goad, M. R.; Grier, C. J.; Horne, Keith; Kochanek, C. S.; Mchardy, I. M.; Netzer, H.; Pancoast, A.; Pei, L.; Pogge, R. W.; Proga, D.; Silva, C.; Tejos, N.; Vestergaard, M.; Adams, S. M.; Anderson, M. D.; Arévalo, P.; Beatty, T G.; Behar, E.; Bennert, V. N.; Bianchi, S.; Bigley, A.; Bisogni, S.; Boissay-Malaquin, R.; Borman, G. A.; Bottorff, M. C.; Breeveld, A. A.; Brotherton, M.; Brown, J. E.; Brown, J. S.; Cackett, E. M.; Canalizo, G.; Cappi, M.; Carini, M. T.; Clubb, K. I.; Comerford, J. M.; Coker, C. T.; Corsini, E. M.; Costantini, E.; Croft, S.; Croxall, K. V.; Deason, A. J.; De Lorenzo-Cáceres, A.; De Marco, B.; Dietrich, M.; Di Gesu, L.; Ebrero, J.; Evans, P. A.; Filippenko, A. V.; Flatland, K.; Gates, E. L.; Gehrels, N.; Geier, Z S.; Gelbord, J. M.; Gonzalez, L.; Gorjian, V.; Grupe, D.; Gupta, A.; Hall, P. B.; Henderson, C. B.; Hicks, S.; Holmbeck, E.; Holoien, T. W. -S.; Hutchison, T. A.; Im, M.; Jensen, J. J.; Johnson, C. A.; Joner, M. D.; Kaspi, S.; Kelly, B. C.; Kelly, P. L.; Kennea, J. A.; Kim, M.; Kim, S. C.; Kim, S. Y.; King, A.; Klimanov, S. A.; Krongold, Y.; Lau, M. W.; Lee, J. C.; Leonard, D. C.; Li, Miao; Lira, P.; Lochhaas, C.; Ma, Zhiyuan; Macinnis, F.; Malkan, M. A.; Manne-Nicholas, E. R.; Matt, G.; Mauerhan, J. C.; Mcgurk, R.; Montuori, C.; Morelli, L.; Mosquera, A.; Mudd, D.; Müller-Sánchez, F.; Nazarov, S. V.; Norris, R. P.; Nousek, J. A.; Nguyen, M. L.; Ochner, P.; Okhmat, D. N.; Paltani, S.; Parks, J. R.; Pinto, C.; Pizzella, A.; Poleski, R.; Ponti, G.; Pott, J. -U.; Rafter, S. E.; Rix, H. -W.; Runnoe, J.; Saylor, D. A.; Schimoia, J. S.; Schnülle, K.; Scott, B.; Sergeev, S. G.; Shappee, B. J.; Shivvers, I.; Siegel, M.; Simonian, G. V.; Siviero, A.; Skielboe, A.; Somers, G.; Spencer, M.; Starkey, D.; Stevens, D. J.; Sung, H. -I.; Tayar, J.; Teems, K. G.; Treu, T.; Turner, C. S.; Uttley, P.; Van Saders, J .; Vican, L.; Villforth, C.; Villanueva, S.; Walton, D. J.; Waters, T.; Weiss, Y.; Woo, J. -H.; Yan, H.; Yuk, H.; Zheng, W.; Zhu, W.; Zu, Y.
Comments: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly
Submitted: 2019-07-08, last modified: 2019-07-12
We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.
[16]  [pdf] - 1953577
The structure of the Orion Nebula in the direction of Theta 1 Ori C
Comments: 13 pages, 9 figures, accepted for publication in the Astrophysical Journal
Submitted: 2019-06-18
We have used existing optical emission and absorption lines, [C II] emission lines, and H I absorption lines to create a new model for a Central Column of material near the Trapezium region of the Orion Nebula. This was necessary because recent high spectral resolution spectra of optical emission lines and imaging spectra in the [C II] 158 micron line have shown that there are new velocity systems associated with the foreground Veil and the material lying between Theta 1 Ori C and the Main Ionization Front of the nebula. When a family of models generated with the spectral synthesis code Cloudy were compared with the surface brightness of the emission lines and strengths of the Veil absorption lines seen in the Trapezium stars, distances from Theta 1 Ori C, were derived, with the closest, highest ionization layer being 1.3 pc. The line of sight distance of this layer is comparable with the size of the inner Huygens Region in the plane of the sky. These layers are all blueshifted with respect to the Orion Nebula Cluster of stars, probably because of the pressure of a hot central bubble created by Theta 1 Ori C's stellar wind. We find velocity components that are ascribed to both sides of this bubble. Our analysis shows that the foreground [C II] 158 micron emission is part of a previously identified layer that forms a portion of a recently discovered expanding shell of material covering most of the larger Extended Orion Nebula.
[17]  [pdf] - 1897896
Space Telescope and Optical Reverberation Mapping Project. X. Understanding the Absorption-Line Holiday in NGC 5548
Comments: This is the version accepted by ApJ. Will be published shortly
Submitted: 2018-12-30, last modified: 2019-04-24
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.
[18]  [pdf] - 1853150
Unlocking the Capabilities of Future High-Resolution X-ray Spectroscopy Missions Through Laboratory Astrophysics
Comments: Science white paper submitted to the Astro2020 Decadal Survey
Submitted: 2019-03-19
Thanks to high-resolution and non-dispersive spectrometers onboard future X-ray missions such as XRISM and Athena, we are finally poised to answer important questions about the formation and evolution of galaxies and large-scale structure. However, we currently lack an adequate understanding of many atomic processes behind the spectral features we will soon observe. Large error bars on parameters as critical as transition energies and atomic cross sections can lead to unacceptable uncertainties in the calculations of e.g., elemental abundance, velocity, and temperature. Unless we address these issues, we risk limiting the full scientific potential of these missions. Laboratory astrophysics, which comprises theoretical and experimental studies of the underlying physics behind observable astrophysical processes, is therefore central to the success of these missions.
[19]  [pdf] - 1864050
H-, He-like recombination spectra III: $n$-changing collisions in highly-excited Rydberg states and their impact on the radio, IR and optical recombination lines
Comments: 16 pages (accepted by MNRAS, replaced)
Submitted: 2019-03-13, last modified: 2019-03-18
At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact parameter calculations or semi-empirical fits and the rate coefficients agree within a factor of two. We show in this paper that these uncertainties cause errors lower than 5% in the emission of radio recombination lines (RRL) of most ionized plasmas of typical nebulae. However, in special circumstances where the transitions between Rydberg levels are amplified by maser effects, the errors can increase up to 20%. We present simulations of the optical depth and H$n\alpha$ line emission of Active Galactic Nuclei (AGN) Broad Line Regions (BLRs) and the Orion Nebula Blister to showcase our findings.
[20]  [pdf] - 1850848
The physics and astrophysics of X-ray outflows from Active Galactic Nuclei
Comments: A Science White Paper submitted to the Astro2020 Decadal Survey
Submitted: 2019-03-15
The highly energetic outflows from Active Galactic Nuclei detected in X-rays are one of the most powerful mechanisms by which the central supermassive black hole (SMBH) interacts with the host galaxy. The last two decades of high resolution X-ray spectroscopy with XMM and Chandra have improved our understanding of the nature of these outflowing ionized absorbers and we are now poised to take the next giant leap with higher spectral resolution and higher throughput observatories to understand the physics and impact of these outflows on the host galaxy gas. The future studies on X-ray outflows not only have the potential to unravel some of the currently outstanding puzzles in astronomy, such as the physical basis behind the MBH$-\sigma$ relation, the cooling flow problem in intra-cluster medium (ICM), and the evolution of the quasar luminosity function across cosmic timescales, but also provide rare insights into the dynamics and nature of matter in the immediate vicinity of the SMBH. Higher spectral resolution ($\le 0.5$ eV at $1$ keV) observations will be required to identify individual absorption lines and study the asymmetries and shifts in the line profiles revealing important information about outflow structures and their impact. Higher effective area ($\ge 1000 \rm \,cm^{2}$) will be required to study the outflows in distant quasars, particularly at the quasar peak era (redshift $1\le z\le 3$) when the AGN population was the brightest. Thus, it is imperative that we develop next generation X-ray telescopes with high spectral resolution and high throughput for unveiling the properties and impact of highly energetic X-ray outflows. A simultaneous high resolution UV + X-ray mission will encompass the crucial AGN ionizing continuum, and also characterize the simultaneous detections of UV and X-ray outflows, which map different spatial scales along the line of sight.
[21]  [pdf] - 1983820
Ubiquitous cold and massive filaments in cool core clusters
Comments: 33 Pages, 17 Figures, submitted to A&A
Submitted: 2019-02-25
Multi-phase filamentary structures around Brightest Cluster Galaxies are likely a key step of AGN-feedback. We observed molecular gas in 3 cool cluster cores: Centaurus, Abell S1101, and RXJ1539.5 and gathered ALMA and MUSE data for 12 other clusters. Those observations show clumpy, massive and long, 3--25 kpc, molecular filaments, preferentially located around the radio bubbles inflated by the AGN (Active Galactic Nucleus). Two objects show nuclear molecular disks. The optical nebula is certainly tracing the warm envelopes of cold molecular filaments. Surprisingly, the radial profile of the H$\alpha$/CO flux ratio is roughly constant for most of the objects, suggesting that (i) between 1.2 to 7 times more cold gas could be present and (ii) local processes must be responsible for the excitation. Projected velocities are between 100--400 km s$^{-1}$, with disturbed kinematics and sometimes coherent gradients. This is likely due to the mixing in projection of several thin unresolved filaments. The velocity fields may be stirred by turbulence induced by bubbles, jets or merger-induced sloshing. Velocity and dispersions are low, below the escape velocity. Cold clouds should eventually fall back and fuel the AGN. We compare the filament's radial extent, r$_{fil}$, with the region where the X-ray gas can become thermally unstable. The filaments are always inside the low-entropy and short cooling time region, where t$_{cool}$/t$_{ff}$<20 (9 of 13 sources). The range t$_{cool}$/t$_{ff}$, 8-23 at r$_{fil}$, is likely due to (i) a more complex gravitational potential affecting the free-fall time (e.g., sloshing, mergers); (ii) the presence of inhomogeneities or uplifted gas in the ICM, affecting the cooling time. For some of the sources, r$_{fil}$ lies where the ratio of the cooling time to the eddy-turnover time, t$_{cool}$/t$_{eddy}$, is approximately unity.
[22]  [pdf] - 1834292
Deep and narrow CO absorption revealing molecular clouds in the Hydra-A brightest cluster galaxy
Comments: 11 pages, 8 figures, accepted to MNRAS
Submitted: 2019-02-05
Active galactic nuclei play a crucial role in the accretion and ejection of gas in galaxies. Although their outflows are well studied, finding direct evidence of accretion has proved very difficult and has so far been done for very few sources. A promising way to study the significance of cold accretion is by observing the absorption of an active galactic nucleus's extremely bright radio emission by the cold gas lying along the line-of-sight. As such, we present ALMA CO(1-0) and CO(2-1) observations of the Hydra-A brightest cluster galaxy (z=0.054) which reveal the existence of cold, molecular gas clouds along the line-of-sight to the galaxy's extremely bright and compact mm-continuum source. They have apparent motions relative to the central supermassive black hole of between -43 and -4 km s$^{-1}$ and are most likely moving along stable, low ellipticity orbits. The identified clouds form part of a $\sim$$10^{9}$ $\text{M}_{\odot}$, approximately edge-on disc of cold molecular gas. With peak CO(2-1) optical depths of $\tau$=0.88 $^{+0.06}_{-0.06}$, they include the narrowest and by far the deepest absorption of this type which has been observed to date in a brightest cluster galaxy. By comparing the relative strengths of the lines for the most strongly absorbing region, we are able to estimate a gas temperature of $42^{+25}_{-11}$ K and line-of-sight column densities of $N_{CO}=2^{+3}_{-1}\times 10 ^{17} cm^{-2}$ and $N_{ H_{2} }=7^{+10}_{-4}\times 10 ^{20} cm^{-2}$.
[23]  [pdf] - 1818759
The First Swift Intensive AGN Accretion Disk Reverberation Mapping Survey
Comments: 25 pages, 9 tables, 6 figures, published in ApJ. The paper has been revised to conform to the published version. Please note that the data (Table 2) are available at the Digital Repository at the University of Maryland,
Submitted: 2018-11-19, last modified: 2019-01-20
Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled $\sim$200-350 times in 0.3-10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation analysis of these datasets yields three main results: 1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. 2) The UV/optical interband lags are generally consistent with $ \tau \propto \lambda^{4/3} $ as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. 3) The one exception is the U band lags, which are on average a factor of $\sim$2.2 larger than predicted from the surrounding band data and fits. This excess appears due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by javelin vs. more standard lag measurement techniques, with javelin reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein.
[24]  [pdf] - 1798384
Fe II velocity shifts in optical spectra of type 1 AGN
Comments: 8 pages, 2 figures, 1 table, presented at conference "Revisiting narrow-line Seyfert 1 galaxies and their place in the Universe"; PoS(NLS1-2018)007
Submitted: 2018-12-14
Here we present a critical view of practical problems in the analysis of optical Fe II emission lines of Type 1 AGN spectra. Besides the very complex and unclear physical interpretation of the Fe II template shape, there are other issues that might affect the results of Fe II contribution as well (like the S/N, the AGN continuum component modelling, the complex shapes of other broad and narrow emission lines that are blending in the same part of the spectrum, the galactic host stellar component that could lead to mimicking of the Fe II template shape on some parts of spectrum, etc.). In this paper we concentrate particularly on the claims of possibility that in some objects Fe II could be strongly shifted to the red. We examine the effects in the fitting procedure that could artificially lead to such results.
[25]  [pdf] - 1709479
Production of photoionized plasmas in the laboratory using X-ray line radiation
Comments: 10 pages; 8 figures; Physical Review E, in press
Submitted: 2018-05-15
In this paper we report the experimental implementation of a theoretically-proposed technique for creating a photoionized plasma in the laboratory using X-ray line radiation. Using a Sn laser-plasma to irradiate an Ar gas target, the photoionization parameter, \xi = 4{\pi}F/Ne, reached values of order 50 erg cm/s, where F is the radiation flux in erg/cm2/s. The significance of this is that this technique allows us to mimic effective spectral radiation temperatures in excess of 1 keV. We show that our plasma starts to be collisionally dominated before the peak of the X-ray drive. However, the technique is extendable to higher energy laser systems to create plasmas with parameters relevant to benchmarking codes used to model astrophysical objects.
[26]  [pdf] - 1660848
Intermediate-line Emission in AGNs: The Effect of Prescription of the Gas Density
Comments: 15 pages, 13 figures
Submitted: 2018-02-28, last modified: 2018-04-04
The requirement of intermediate line component in the recently observed spectra of several AGNs points to possibility of the existence of a physically separate region between broad line region (BLR) and narrow line region (NLR). In this paper we explore the emission from intermediate line region (ILR) by using the photoionization simulations of the gas clouds distributed radially from the AGN center. The gas clouds span distances typical for BLR, ILR and NLR, and the appearance of dust at the sublimation radius is fully taken into account in our model. Single cloud structure is calculated under the assumption of the constant pressure. We show that the slope of the power law cloud density radial profile does not affect the existence of ILR in major types of AGN. We found that the low ionization iron line, Fe~II, appears to be highly sensitive for the presence of dust and therefore becomes potential tracer of dust content in line emitting regions. We show that the use of disk-like cloud density profile computed at the upper part of the accretion disc atmosphere reproduces the observed properties of the line emissivities. In particular, the distance of H${\beta}$ line inferred from our model agrees with that obtained from the reverberation mapping studies in Sy1 galaxy NGC 5548.
[27]  [pdf] - 1588686
On the Intermediate Line Region in AGNs
Comments: 8 pages, 2 figures, Accepted for publication in the Journal Frontiers in Astronomy and Space Sciences
Submitted: 2017-09-21
In this paper we explore the intermediate line region (ILR) by using the photoionisation simulations of the gas clouds present at different radial distances from the center, corresponding to the locations from BLR out to NLR in four types of AGNs. We let for the presence of dust whenever conditions allow for dust existence. All spectral shapes are taken from the recent multi-wavelength campaigns. The cloud density decreases with distance as a power law. We found that the slope of the power law density profile does not affect the line emissivity radial profiles of major emission lines: H${\beta}$, He~II, Mg~II, C~III] ~and [O~III]. When the density of the cloud at the sublimation radius is as high as 10$^{11.5}$ cm$^{-3}$, the ILR should clearly be seen in the observations independently of the shape of the illuminating radiation. Moreover, our result is valid for low ionization nuclear emission regions of active galaxies.
[28]  [pdf] - 1587005
Carbon and Nitrogen Abundance Ratio In the Broad Line Region of Tidal Disruption Events
Comments: 10 pages, 5 figures
Submitted: 2017-08-11, last modified: 2017-08-16
The rest-frame UV spectra of three recent tidal disruption events (TDEs), ASASSN-14li, PTF15af and iPTF16fnl display strong nitrogen emission lines but weak or undetectable carbon lines. In these three objects, the upper limits of the C III] $\lambda 1908$/N III] $\lambda 1750$ ratio are about two orders of magnitude lower than those of quasars, suggesting a high abundance ratio of [N/C]. With detailed photoionization simulations, we demonstrate that $C^{2+}$ and $N^{2+}$ are formed in the same zone, so the C III]/N III] ratio depends only moderately on the physical conditions in the gas and weakly on the shape of the ionizing continuum. There are smaller than $0.5$ dex variations in the line ratio over wide ranges of gas densities and ionization parameters at a given metallicity. This allows a robust estimate of the relative abundance ratio nitrogen to carbon. We derive a relative abundance ratio of [N/C]$>1.5$ for ASASSN-14li, and even higher for PTF15af and iPTF16fnl. This suggests that the broad line region in those TDE sources is made of nitrogen-enhanced core material that falls back at later times. Based on stellar evolution models, the lower limit of the disrupted star should be larger than 0.6M$_{\odot}$. The chemical abundance of the line emitting gas gives a convincing evidence that the flares origin from stellar tidal disruptions. The coincidence of the weakness of the X-ray emission with the strong broad absorption lines in PTF15af, iPTF16fnl and the strong X-ray emission without such lines in ASASSN-li14 are in analogy to quasars with and without broad absorption lines.
[29]  [pdf] - 1582748
Space Telescope and Optical Reverberation Mapping Project. VII. Understanding the UV anomaly in NGC 5548 with X-Ray Spectroscopy
Mathur, S.; Gupta, A.; Page, K.; Pogge, R. W.; Krongold, Y.; Goad, M. R.; Adams, S. M.; Anderson, M. D.; Arevalo, P.; Barth, A. J.; Bazhaw, C.; Beatty, T. G.; Bentz, M. C.; Bigley, A.; Bisogni, S.; Borman, G. A.; Boroson, T. A.; Bottorff, M. C.; Brandt, W. N.; Breeveld, A. A.; Brown, J. E.; Brown, J. S.; Cackett, E. M.; Canalizo, G.; Carini, M. T.; Clubb, K. I.; Comerford, J. M.; Coker, C. T.; Corsini, E. M.; Crenshaw, D. M.; Croft, S.; Croxall, K. V.; Bonta, E. Dalla; Deason, A. J.; Denney, K. D.; De Lorenzo-Caceres, A.; De Rosa, G.; Dietrich, M.; Edelson, R.; Ely, J.; Eracleous, M.; Evans, P. A.; Fausnaugh, M. M.; Ferland, G. J.; Filippenko, A. V.; Flatland, K.; Fox, O. D.; Gates, E. L.; Gehrels, N.; Geier, S.; Gelbord, J. M.; Gorjian, V.; Greene, J. E.; Grier, C. J.; Grupe, D.; Hall, P. B.; Henderson, C. B.; Hicks, S.; Holmbeck, E.; Holoien, T. W. -S.; Horenstein, D.; Horne, Keith; Hutchison, T.; Im, M.; Jensen, J. J.; Johnson, C. A.; Joner, M. D.; Jones, J.; Kaastra, J.; Kaspi, S.; Kelly, B. C.; Kelly, P. L.; Kennea, J. A.; Kim, M.; Kim, S.; Kim, S. C.; King, A.; Klimanov, S. A.; Kochanek, C. S.; Korista, K. T.; Kriss, G. A.; Lau, M. W.; Lee, J. C.; Leonard, D. C.; Li, M.; Lira, P.; Ma, Z.; MacInnis, F.; Manne-Nicholas, E. R.; Malkan, M. A.; Mauerhan, J. C.; McGurk, R.; McHardy, I. M.; Montouri, C.; Morelli, L.; Mosquera, A.; Mudd, D.; Muller-Sanchez, F.; Musso, R.; Nazarov, S. V.; Netzer, H.; Nguyen, M. L.; Norris, R. P.; Nousek, J. A.; Ochner, P.; Okhmat, D. N.; Ou-Yang, B.; Pancoast, A.; Papadakis, I.; Parks, J. R.; Pei, L.; Peterson, B. M.; Pizzella, A.; Poleski, R.; Pott, J. -U.; Rafter, S. E.; Rix, H. -W.; Runnoe, J.; Saylor, D. A.; Schimoia, J. S.; Schnülle, K.; Sergeev, S. G.; Shappee, B. J.; Shivvers, I.; Siegel, M.; Simonian, G. V.; Siviero, A.; Skielboe, A.; Somers, G.; Spencer, M.; Starkey, D.; Stevens, D. J.; Sung, H. -I.; Tayar, J.; Tejos, N.; Turner, C. S.; Uttley, P.; Van Saders, J.; Vestergaard, M.; Vican, L.; VillanuevaJr., S.; Villforth, C.; Weiss, Y.; Woo, J. -H.; Yan, H.; Young, S.; Yuk, H.; Zheng, W.; Zhu, W.; Zu, Y.
Comments: ApJ in press. Replaced with the accepted version
Submitted: 2017-04-20, last modified: 2017-08-01
During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.
[30]  [pdf] - 1592077
A Magnetically-Supported Photodissociation Region in M17
Comments: Published as 2007, ApJ,658,1119
Submitted: 2006-11-27, last modified: 2017-08-01
The southwestern (SW) part of the Galactic H II region M17 contains an obscured ionization front that is most easily seen at infrared and radio wavelengths. It is nearly edge-on, thus offering an excellent opportunity to study the way in which the gas changes from fully ionized to molecular as radiation from the ionizing stars penetrates into the gas. M17 is also one of the very few H II regions for which the magnetic field strength can be measured in the photodissociation region ( PDR) that forms the interface between the ionized and molecular gas. Here we model an observed line of sight through the gas cloud, including the H+, H0 (PDR), and molecular layers, in a fully self-consistent single calculation. An interesting aspect of the M17 SW bar is that the PDR is very extended. We show that the strong magnetic field that is observed to be present inevitably leads to a very deep PDR, because the structure of the neutral and molecular gas is dominated by magnetic pressure, rather than by gas pressure, as previously had been supposed.We also show that a wide variety of observed facts can be explained if a hydrostatic geometry prevails, in which the gas pressure from an inner X-ray hot bubble and the outward momentum of the stellar radiation field compress the gas and its associated magnetic field in the PDR, as has already been shown to occur in the Orion Nebula. The magnetic field compression may also amplify the local cosmic-ray density. The pressure in the observed magnetic field balances the outward forces, suggesting that the observed geometry is a natural consequence of the formation of a star cluster within a molecular cloud.
[31]  [pdf] - 1584071
The 2017 Release of Cloudy
Comments: Revista Mexicana accepted, to appear in Vol. 53, N{\deg}2, October 2017
Submitted: 2017-05-30, last modified: 2017-06-28
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.
[32]  [pdf] - 1581198
Atomic Data Revisions for Transitions Relevant to Observations of Interstellar, Circumgalactic, and Intergalactic Matter
Comments: Final version published in the Astrophysical Journal Supplement Series. Includes 1 machine-readable table
Submitted: 2017-01-17, last modified: 2017-06-16
Measurements of element abundances in galaxies from astrophysical spectroscopy depend sensitively on the atomic data used. With the goal of making the latest atomic data accessible to the community, we present a compilation of selected atomic data for resonant absorption lines at wavelengths longward of 911.753 {\AA} (the \ion{H}{1} Lyman limit), for key heavy elements (heavier than atomic number 5) of astrophysical interest. In particular, we focus on the transitions of those ions that have been observed in the Milky Way interstellar medium (ISM), the circumgalactic medium (CGM) of the Milky Way and/or other galaxies, and the intergalactic medium (IGM). We provide wavelengths, oscillator strengths, associated accuracy grades, and references to the oscillator strength determinations. We also attempt to compare and assess the recent oscillator strength determinations. For about 22\% of the lines that have updated oscillator strength values, the differences between the former values and the updated ones are $\gtrsim$~0.1 dex. Our compilation will be a useful resource for absorption line studies of the ISM, as well as studies of the CGM and IGM traced by sight lines to quasars and gamma-ray bursts. Studies (including those enabled by future generations of extremely large telescopes) of absorption by galaxies against the light of background galaxies will also benefit from our compilation.
[33]  [pdf] - 1584262
The validity of 21 cm spin temperature as a kinetic temperature indicator in atomic and molecular gas
Comments: Accepted for publication in ApJ
Submitted: 2017-06-05
The gas kinetic temperature ($T_K$) of various interstellar environments is often inferred from observations that can deduce level populations of atoms, ions, or molecules using spectral line observations; H I 21 cm is perhaps the most widely used with a long history. Usually the H I 21 cm line is assumed to be in thermal equilibrium and the populations are given by the Boltzmann distribution. A variety of processes, many involving Lyman alpha ($Ly\alpha$), can affect the 21 cm line. Here we show how this is treated in the spectral simulation code Cloudy, and present numerical simulations of environments where this temperature indicator is used, with a detailed treatment of the physical processes that determine level populations within $H^0$. We discuss situations where this temperature indicator traces $T_K$, cases where they fail, as well as the effects of $Ly\alpha$ pumping on the 21 cm spin temperature. We also show that the $Ly\alpha$ excitation temperature rarely traces the gas kinetic temperature.
[34]  [pdf] - 1582973
Emission line ratios of Fe III as astrophysical plasma diagnostics
Comments: Accepted for publication in the ApJ
Submitted: 2017-04-29
Recent state-of-the-art calculations of A-values and electron impact excitation rates for Fe III are used in conjunction with the Cloudy modeling code to derive emission line intensity ratios for optical transitions among the fine-structure levels of the 3d$^6$ configuration. A comparison of these with high resolution, high signal-to-noise spectra of gaseous nebulae reveals that previous discrepancies found between theory and observation are not fully resolved by the latest atomic data. Blending is ruled out as a likely cause of the discrepancies, because temperature- and density-independent ratios (arising from lines with common upper levels) match well with those predicted by theory. For a typical nebular plasma with electron temperature $T_{\rm e} = 9000$ K and electron density $\rm N_{e}=10^4 \, cm^{-3}$, cascading of electrons from the levels $\rm ^3G_5$, $\rm ^3G_4$ and $\rm ^3G_3$ plays an important role in determining the populations of lower levels, such as $\rm ^3F_4$, which provide the density diagnostic emission lines of Fe III, such as $\rm ^5D_4$ - $\rm ^3F_4$ at 4658 \AA. Hence further work on the A-values for these transitions is recommended, ideally including measurements if possible. However, some Fe III ratios do provide reliable $N_{\rm e}$-diagnostics, such as 4986/4658. The Fe III cooling function calculated with Cloudy using the most recent atomic data is found to be significantly greater at $T_e$ $\simeq$ 30000 K than predicted with the existing Cloudy model. This is due to the presence of additional emission lines with the new data, particularly in the 1000--4000 \AA\ wavelength region.
[35]  [pdf] - 1576461
Thermodynamically-consistent semi-classical $\ell$-changing rates
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
Submitted: 2017-04-27
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 $\Delta\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.
[36]  [pdf] - 1574689
Swift monitoring of NGC 4151: Evidence for a Second X-ray/UV Reprocessing
Comments: 14 pages, 6 figures, 7 tables. To appear in ApJ April 2017 issue. This version incorporates minor revisions to conform to published version and corrects the links to the author names above
Submitted: 2017-03-20, last modified: 2017-04-04
Swift monitoring of NGC 4151 with ~6 hr sampling over a total of 69 days in early 2016 is used to construct light curves covering five bands in the X-rays (0.3-50 keV) and six in the ultraviolet (UV)/optical (1900-5500 A). The three hardest X-ray bands (>2.5 keV) are all strongly correlated with no measurable interband lag while the two softer bands show lower variability and weaker correlations. The UV/optical bands are significantly correlated with the X-rays, lagging ~3-4 days behind the hard X-rays. The variability within the UV/optical bands is also strongly correlated, with the UV appearing to lead the optical by ~0.5-1 day. This combination of >~3 day lags between the X-rays and UV and <~1 day lags within the UV/optical appears to rule out the "lamp-post" reprocessing model in which a hot, X-ray emitting corona directly illuminates the accretion disk, which then reprocesses the energy in the UV/optical. Instead, these results appear consistent with the Gardner & Done picture in which two separate reprocessings occur: first, emission from the corona illuminates an extreme-UV-emitting toroidal component that shields the disk from the corona; this then heats the extreme-UV component which illuminates the disk and drives its variability.
[37]  [pdf] - 1549550
Which Stars are Ionizing the Orion Nebula ?
Comments: 7 pages with 5 figures
Submitted: 2017-03-17
The common assumption that Theta-1-Ori C is the dominant ionizing source for the Orion Nebula is critically examined. This assumption underlies much of the existing analysis of the nebula. In this paper we establish through comparison of the relative strengths of emission lines with expectations from Cloudy models and through the direction of the bright edges of proplyds that Theta-2-Ori-A, which lies beyond the Bright Bar, also plays an important role. Theta-1-Ori-C does dominate ionization in the inner part of the Orion Nebula, but outside of the Bright Bar as far as the southeast boundary of the Extended Orion Nebula, Theta-2-Ori-A is the dominant source. In addition to identifying the ionizing star in sample regions, we were able to locate those portions of the nebula in 3-D. This analysis illustrates the power of MUSE spectral imaging observations in identifying sources of ionization in extended regions.
[38]  [pdf] - 1550468
Space Telescope and Optical Reverberation Mapping Project. V. Optical Spectroscopic Campaign and Emission-Line Analysis for NGC 5548
Pei, L.; Fausnaugh, M. M.; Barth, A. J.; Peterson, B. M.; Bentz, M. C.; De Rosa, G.; Denney, K. D.; Goad, M. R.; Kochanek, C. S.; Korista, K. T.; Kriss, G. A.; Pogge, R. W.; Bennert, V. N.; Brotherton, M.; Clubb, K. I.; Bontà, E. Dalla; Filippenko, A. V.; Greene, J. E.; Grier, C. J.; Vestergaard, M.; Zheng, W.; Adams, Scott M.; Beatty, Thomas G.; Bigley, A.; Brown, Jacob E.; Brown, Jonathan S.; Canalizo, G.; Comerford, J. M.; Coker, Carl T.; Corsini, E. M.; Croft, S.; Croxall, K. V.; Deason, A. J.; Eracleous, Michael; Fox, O. D.; Gates, E. L.; Henderson, C. B.; Holmbeck, E.; Holoien, T. W. -S.; Jensen, J. J.; Johnson, C. A.; Kelly, P. L.; Kim, S.; King, A.; Lau, M. W.; Li, Miao; Lochhaas, Cassandra; Ma, Zhiyuan; Manne-Nicholas, E. R.; Mauerhan, J. C.; Malkan, M. A.; McGurk, R.; Morelli, L.; Mosquera, Ana; Mudd, Dale; Sanchez, F. Muller; Nguyen, M. L.; Ochner, P.; Ou-Yang, B.; Pancoast, A.; Penny, Matthew T.; Pizzella, A.; Poleski, Radosław; Runnoe, Jessie; Scott, B.; Schimoia, Jaderson S.; Shappee, B. J.; Shivvers, I.; Simonian, Gregory V.; Siviero, A.; Somers, Garrett; Stevens, Daniel J.; Strauss, M. A.; Tayar, Jamie; Tejos, N.; Treu, T.; Van Saders, J.; Vican, L.; Villanueva, S.; Yuk, H.; Zakamska, N. L.; Zhu, W.; Anderson, M. D.; Arévalo, P.; Bazhaw, C.; Bisogni, S.; Borman, G. A.; Bottorff, M. C.; Brandt, W. N.; Breeveld, A. A.; Cackett, E. M.; Carini, M. T.; Crenshaw, D. M.; De Lorenzo-Cáceres, A.; Dietrich, M.; Edelson, R.; Efimova, N. V.; Ely, J.; Evans, P. A.; Ferland, G. J.; Flatland, K.; Gehrels, N.; Geier, S.; Gelbord, J. M.; Grupe, D.; Gupta, A.; Hall, P. B.; Hicks, S.; Horenstein, D.; Horne, Keith; Hutchison, T.; Im, M.; Joner, M. D.; Jones, J.; Kaastra, J.; Kaspi, S.; Kelly, B. C.; Kennea, J. A.; Kim, M.; Kim, S. C.; Klimanov, S. A.; Lee, J. C.; Leonard, D. C.; Lira, P.; MacInnis, F.; Mathur, S.; McHardy, I. M.; Montouri, C.; Musso, R.; Nazarov, S. V.; Netzer, H.; Norris, R. P.; Nousek, J. A.; Okhmat, D. N.; Papadakis, I.; Parks, J. R.; Pott, J. -U.; Rafter, S. E.; Rix, H. -W.; Saylor, D. A.; Schnülle, K.; Sergeev, S. G.; Siegel, M.; Skielboe, A.; Spencer, M.; Starkey, D.; Sung, H. -I.; Teems, K. G.; Turner, C. S.; Uttley, P.; Villforth, C.; Weiss, Y.; Woo, J. -H.; Yan, H.; Young, S.; Zu, Y.
Comments: 22 pages, 13 figures, accepted to ApJ
Submitted: 2017-02-03
We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.
[39]  [pdf] - 1548696
Testing atomic collision theory with the two-photon continuum of astrophysical nebulae
Comments: 7 pages, 6 fig, 2 tables, accepted in MNRAS
Submitted: 2017-01-26
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.
[40]  [pdf] - 1500509
Structure and Physical Conditions in the Huygens Region of the Orion Nebula
Submitted: 2016-10-20
HST images, MUSE maps of emission-lines, and an atlas of high velocity resolution emission-line spectra have been used to establish for the firrst time correlations of the electron temperature, electron density, radial velocity, turbulence, and orientation within the main ionization front of the nebula. From the study of the combined properties of multiple features, it is established that variations in the radial velocity are primarily caused by the photo-evaporating ionization front being viewed at different angles. There is a progressive increase of the electron temperature and density with decreasing distance from the dominant ionizing star Theta1 Ori C. The product of these characteristics (NexTe) is the most relevant parameter in modeling a blister-type nebula like the Huygens Region, where this quantity should vary with the surface brightness in Halpha. Several lines of evidence indicate that small-scale structure and turbulence exists down to the level of our resolution of a few arcseconds. Although photo-evaporative ow must contribute at some level to the well-known non-thermal broadening of the emission lines, comparison of quantitative predictions with the observed optical line widths indicate that it is not the major additive broad- ening component. Derivation of Te values for H+ from radio+optical and optical-only ionized hydro- gen emission showed that this temperature is close to that derived from [Nii] and that the transition from the well-known at extinction curve that applies in the Huygens Region to a more normal steep extinction curve occurs immediately outside of the Bright Bar feature of the nebula.
[41]  [pdf] - 1483511
H, He-like recombination spectra II: $l$-changing collisions for He Rydberg states
Comments: 9 pages, $ figures, 2 tables, accepted for publication in MNRAS
Submitted: 2016-08-31
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.
[42]  [pdf] - 1472985
Evidence for Fluorescent Fe II Emission from Extended Low Ionization Outflows in Obscured Quasars
Comments: 24 pages, 6 figures
Submitted: 2016-08-31
Recent studies have shown that outflows in at least some broad absorption line (BAL) quasars are extended well beyond the putative dusty torus. Such outflows should be detectable in obscured quasars. We present four WISE selected infrared red quasars with very strong and peculiar ultraviolet Fe ii emission lines: strong UV Fe II UV arising from transitions to ground/low excitation levels, and very weak Fe II at wavelengths longer than 2800 {\AA}. The spectra of these quasars display strong resonant emission lines, such as C IV, Al III and Mg II but sometimes, a lack of non-resonant lines such as C III], S III and He II. We interpret the Fe II lines as resonantly scattered light from the extended outflows that are viewed nearly edge-on, so that the accretion disk and broad line region are obscured by the dusty torus, while the extended outflows are not. We show that dust free gas exposed to strong radiation longward of 912 {\AA} produces Fe II emission very similar to that observed. The gas is too cool to collisionally excite Fe II lines, accounting for the lack of optical emission. The spectral energy distribution from the UV to the mid-infrared can be modeled as emission from a clumpy dusty torus, with UV emission being reflected/scattered light either by the dusty torus or the outflow. Within this scenario, we estimate a minimum covering factor of the outflows from a few to 20% for the Fe II scattering region, suggesting that Fe II BAL quasars are at a special stage of quasar evolution.
[43]  [pdf] - 1515729
The intermediate line region in active galactic nuclei
Comments: 9 pages, 6 figures, accepted for publication in ApJ
Submitted: 2016-06-01, last modified: 2016-08-16
We show that the recently observed suppression of the gap between the broad line region (BLR) and the narrow line region (NLR) in some AGN can be fully explained by an increase of the gas density in the emitting region. Our model predicts the formation of the intermediate line region (ILR) that is observed in some Seyfert galaxies by the detection of emission lines with intermediate velocity full width half maximum (FWHM) $\sim$ 700 - 1200 km s$^{-1}$. These lines are believed to be originating from an ILR located somewhere between the BLR and NLR. As it was previously proved, the apparent gap is assumed to be caused by the presence of dust beyond the sublimation radius. Our computations with the use of {\sc cloudy} photoionization code, show that the differences in the shape of spectral energy distribution (SED) from the central region of AGN, do not diminish the apparent gap in the line emission in those objects. A strong discontinuity in the line emission vs radius exists for all lines at the dust sublimation radius. However, increasing the gas density to $\sim$ 10$^{11.5}$ cm$^{-3}$ at the sublimation radius provides the continuous line emission vs radius and fully explains the recently observed lack of apparent gap in some AGN. We show that such a high density is consistent with the density of upper layers of an accretion disk atmosphere. Therefore, the upper layers of the disk atmosphere can give rise to the formation of observed emission line clouds.
[44]  [pdf] - 1422401
HST imaging of the dusty filaments and nucleus swirl in NGC4696 at the centre of the Centaurus Cluster
Comments: 7 pages, 10 figures, MNRAS in press
Submitted: 2016-06-08
Narrow-band HST imaging has resolved the detailed internal structure of the 10 kpc diameter H alpha+[NII] emission line nebulosity in NGC4696, the central galaxy in the nearby Centaurus cluster, showing that the dusty, molecular, filaments have a width of about 60pc. Optical morphology and velocity measurements indicate that the filaments are dragged out by the bubbling action of the radio source as part of the AGN feedback cycle. Using the drag force we find that the magnetic field in the filaments is in approximate pressure equipartition with the hot gas. The filamentary nature of the cold gas continues inward, swirling around and within the Bondi accretion radius of the central black hole, revealing the magnetic nature of the gas flows in massive elliptical galaxies. HST imaging resolves the magnetic, dusty, molecular filaments at the centre of the Centaurus cluster to a swirl around and within the Bondi radius.
[45]  [pdf] - 1405609
The abundance discrepancy factor and t^2 in nebulae: are non-thermal electrons the culprits?
Comments: to appear in Revista Mexicana de Astronomia y Astrofisica
Submitted: 2016-05-11
Photoionization produces supra-thermal electrons, electrons with much more energy than is found in a thermalized gas at electron temperatures characteristic of nebulae. The presence of these high energy electrons may solve the long-standing t^2/ADF puzzle, the observations that abundances obtained from recombination and collisionally excited lines do not agree, and that different temperature indicators give different results, if they survive long enough to affect diagnostic emission lines. The presence of these non-Maxwellian distribution electrons is usually designated by the term kappa. Here we use well-established methods to show that the distance over which heating rates change are much longer than the distance supra thermal electrons can travel, and that the timescale to thermalize these electrons are much shorter than the heating or cooling timescales. These estimates establish that supra thermal electrons will have disappeared into the Maxwellian velocity distribution long before they affect the collisionally excited forbidden and recombination lines that are used for deriving abundances relative to hydrogen. The electron velocity distribution in nebulae should be closely thermal.
[46]  [pdf] - 1433034
Ultraviolet emission lines of Si II in quasars --- investigating the "Si II disaster"
Comments: Accepted for publication in ApJ
Submitted: 2016-04-25
The observed line intensity ratios of the Si II 1263 and 1307 \AA\ multiplets to that of Si II 1814\,\AA\ in the broad line region of quasars are both an order of magnitude larger than the theoretical values. This was first pointed out by Baldwin et al. (1996), who termed it the "Si II disaster", and it has remained unresolved. We investigate the problem in the light of newly-published atomic data for Si II. Specifically, we perform broad line region calculations using several different atomic datasets within the CLOUDY modeling code under optically thick quasar cloud conditions. In addition, we test for selective pumping by the source photons or intrinsic galactic reddening as possible causes for the discrepancy, and also consider blending with other species. However, we find that none of the options investigated resolves the Si II disaster, with the potential exception of microturbulent velocity broadening and line blending. We find that a larger microturbulent velocity ($\sim 500 \rm \, kms^{-1}$) may solve the Si II disaster through continuum pumping and other effects. The CLOUDY models indicate strong blending of the Si II 1307 \AA\ multiplet with emission lines of O I, although the predicted degree of blending is incompatible with the observed 1263/1307 intensity ratios. Clearly, more work is required on the quasar modelling of not just the Si II lines but also nearby transitions (in particular those of O I) to fully investigate if blending may be responsible for the Si II disaster.
[47]  [pdf] - 1397025
H, He-like recombination spectra I: $l$-changing collisions for hydrogen
Comments: 6 pages, 3 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2016-04-13
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.
[48]  [pdf] - 1382240
The Lexington Benchmarks for Numerical Simulations of Nebulae
Comments: The Analysis of Emission Lines, Proc. STScI Symposium May 16--18, 1994, Eds. R.E. Williams & M. Livio, Cambridge University Press, p. 83
Submitted: 2016-03-29
We present the results of a meeting on numerical simulations of ionized nebulae held at the University of Kentucky in conjunction with the celebration of the 70th birthdays of Profs. Donald Osterbrock and Michael Seaton.
[49]  [pdf] - 1422270
Space Telescope and Optical Reverberation Mapping Project. IV. Anomalous behavior of the broad ultraviolet emission lines in NGC 5548
Comments: 12 pages, 2 figures, accepted for publication in ApJ, Friday 25th March 2016. A movie of the anomalous emission-line behavior can be found in the ancilliary documents
Submitted: 2016-03-29
During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation.
[50]  [pdf] - 1378803
Interpreting the Ionization Sequence in Star-Forming Galaxy Emission-Line Spectra
Comments: 41 pages, 17 figures, 9 tables, accepted to MNRAS
Submitted: 2016-01-13
High ionization star forming (SF) galaxies are easily identified with strong emission line techniques such as the BPT diagram, and form an obvious ionization sequence on such diagrams. We use a locally optimally emitting cloud model to fit emission line ratios that constrain the excitation mechanism, spectral energy distribution, abundances and physical conditions along the star-formation ionization sequence. Our analysis takes advantage of the identification of a sample of pure star-forming galaxies, to define the ionization sequence, via mean field independent component analysis. Previous work has suggested that the major parameter controlling the ionization level in SF galaxies is the metallicity. Here we show that the observed SF- sequence could alternatively be interpreted primarily as a sequence in the distribution of the ionizing flux incident on gas spread throughout a galaxy. Metallicity variations remain necessary to model the SF-sequence, however, our best models indicate that galaxies with the highest and lowest observed ionization levels (outside the range -0.37 < log [O III]/H\b{eta} < -0.09) require the variation of an additional physical parameter other than metallicity, which we determine to be the distribution of ionizing flux in the galaxy.
[51]  [pdf] - 1334062
The influence of soft spectral components on the structure and stability of warm absorbers in AGN
Comments: 16 pages, 10 figures, Published in MNRAS. Typos corrected, references elaborated and journal reference added. The definitive version is available at
Submitted: 2012-01-25, last modified: 2015-12-29
The radiation from the central regions of active galactic nuclei, including that from the accretion disk surrounding the black hole, is likely to peak in the extreme ultraviolet $\sim 13 -100$ eV. However, due to Galactic absorption, we are limited to constrain the physical properties, i.e. the black hole mass and the accretion rate, from what observations we have below $\sim 10$ eV or above $\sim 100$ eV. In this paper we predict the thermal and ionization states of warm absorbers as a function of the shape of the unobservable continuum. In particular we model an accretion disk at $kT_{in} \sim 10$ eV and a {\it soft excess} at $kT_{se} \sim 150$ eV. The warm absorber, which is the highly ionized gas along the line of sight to the continuum, shows signatures in the $\sim 0.3 - 2$ keV energy range consisting of numerous absorption lines and edges of various ions, some of the prominent ones being H- and He-like oxygen, neon, magnesium and silicon. We find that the properties of the warm absorber are significantly influenced by the changes in the temperature of the accretion disk, as well as by the strength of the {\it soft excess}, as they affect the optical depth particularly for iron and oxygen. These trends may help develop a method of characterising the shape of the unobservable continuum and the occurrence of warm absorbers.
[52]  [pdf] - 1323563
Ultraviolet emission lines of Si II in cool star and solar spectra
Comments: Accepted for publication in MNRAS
Submitted: 2015-10-29
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.
[53]  [pdf] - 1323518
Evidence for Photoionization Driven Broad Absorption Line Variability
Comments: 41 pages, 15 figures, accepted to ApJ
Submitted: 2015-10-16
We present a qualitative analysis of the variability of quasar broad absorption lines using the large multi-epoch spectroscopic dataset of the Sloan Digital Sky Survey Data Release 10. We confirm that variations of absorption lines are highly coordinated among different components of the same ion or the same absorption component of different ions for C IV, Si IV and N V. Furthermore, we show that the equivalent widths of the lines decrease or increase statistically when the continuum brightens or dims. This is further supported by the synchronized variations of emission and absorption line equivalent width, when the well established intrinsic Baldwin effect for emission lines is taken into account. We find that the emergence of an absorption component is usually accompanying with dimming of the continuum while the disappearance of an absorption line component with brightening of the continuum. This suggests that the emergence or disappearance of a C IV absorption component is only the extreme case, when the ionic column density is very sensitive to continuum variations or the continuum variability amplitude is larger. These results support the idea that absorption line variability is driven mainly by changes in the gas ionization in response to continuum variations, that the line-absorbing gas is highly ionized, and in some extreme cases, too highly ionized to be detected in UV absorption lines. Due to uncertainties in the spectroscopic flux calibration, we cannot quantify the fraction of quasars with asynchronized continuum and absorption line variations.
[54]  [pdf] - 1300329
Strong variability of the coronal line region in NGC 5548
Comments: 10 pages, 6 figures; accepted for publication in MNRAS. arXiv admin note: substantial text overlap with arXiv:1501.02928
Submitted: 2015-09-19
We present the second extensive study of the coronal line variability in an active galaxy. Our data set for the well-studied Seyfert galaxy NGC 5548 consists of five epochs of quasi-simultaneous optical and near-infrared spectroscopy spanning a period of about five years and three epochs of X-ray spectroscopy overlapping in time with it. Whereas the broad emission lines and hot dust emission varied only moderately, the coronal lines varied strongly. However, the observed high variability is mainly due to a flux decrease. Using the optical [FeVII] and X-ray OVII emission lines we estimate that the coronal line gas has a relatively low density of n~10^3/cm^3 and a relatively high ionisation parameter of log U~1. The resultant distance of the coronal line gas from the ionising source of about eight light years places this region well beyond the hot inner face of the dusty torus. These results imply that the coronal line region is an independent entity. We find again support for the X-ray heated wind scenario of Pier & Voit; the increased ionising radiation that heats the dusty torus also increases the cooling efficiency of the coronal line gas, most likely due to a stronger adiabatic expansion. The much stronger coronal line variability of NGC 5548 relative to that of NGC 4151 can also be explained within this picture. NGC 5548 has much stronger coronal lines relative to the low ionisation lines than NGC 4151 indicating a stronger wind, in which case a stronger adiabatic expansion of the gas and so fading of the line emission is expected.
[55]  [pdf] - 1281002
The Nature and Frequency of Outflows from Stars in the Central Orion Nebula Cluster
Comments: 152 pages, 46 figures, 7 tables. Accepted by AJ
Submitted: 2015-08-13
Recent Hubble Space Telescope images have allowed the determination with unprecedented accuracy of motions and changes of shocks within the inner Orion Nebula. These originate from collimated outflows from very young stars, some within the ionized portion of the nebula and others within the host molecular cloud. We have doubled the number of Herbig-Haro objects known within the inner Orion Nebula. We find that the best-known Herbig-Haro shocks originate from a relatively few stars, with the optically visible X-ray source COUP 666 driving many of them. While some isolated shocks are driven by single collimated outflows, many groups of shocks are the result of a single stellar source having jets oriented in multiple directions at similar times. This explains the feature that shocks aligned in opposite directions in the plane of the sky are usually blue shifted because the redshifted outflows pass into the optically thick Photon Dominated Region behind the nebula. There are two regions from which optical outflows originate for which there are no candidate sources in the SIMBAD data base.
[56]  [pdf] - 1259914
Herschel dust emission as a probe of starless cores mass: MCLD 123.5+24.9 of the Polaris Flare
Comments: ApJ, Accepted. Minor typographical errors corrected and figures 6 & 7 updated in v3
Submitted: 2015-04-23, last modified: 2015-06-22
We present newly processed archival Herschel images of molecular cloud MCLD 123.5+24.9 in the Polaris Flare. This cloud contains five starless cores. Using the spectral synthesis code Cloudy, we explore uncertainties in the derivation of column densities, hence, masses of molecular cores from Herschel data. We first consider several detailed grain models that predict far-IR grain opacities. Opacities predicted by the models differ by more than a factor of two, leading to uncertainties in derived column densities by the same factor. Then we consider uncertainties associated with the modified blackbody fitting process used by observers to estimate column densities. For high column density clouds (N(H) $\gg$ 10$^{22}$ cm$^{-2}$), this fitting technique can underestimate column densities by about a factor of three. Finally, we consider the virial stability of the five starless cores in MCLD 123.5+24.9. All of these cores appear to have strongly sub-virial masses, assuming, as we argue, that $^{13}$CO line data provide reliable estimates of velocity dispersions. Evidently, they are not self-gravitating, so it is no surprise that they are starless.
[57]  [pdf] - 1245847
Stout: Cloudy's Atomic and Molecular Database
Comments: ApJ in press
Submitted: 2015-06-04
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.
[58]  [pdf] - 1045085
Atomic Data and Spectral Models for FeII
Comments: 36 pages, 12 figures, 9 tables, accepted for publication in the ApJ
Submitted: 2015-05-20
We present extensive calculations of radiative transition rates and electron impact collision strengths for Fe II. The data sets involve 52 levels from the $3d\,^7$, $3d\,^64s$, and $3d\,^54s^2$ configurations. Computations of $A$-values are carried out with a combination of state-of-the-art multiconfiguration approaches, namely the relativistic Hartree--Fock, Thomas--Fermi--Dirac potential, and Dirac--Fock methods; while the $R$-matrix plus intermediate coupling frame transformation, Breit--Pauli $R$-matrix and Dirac $R$-matrix packages are used to obtain collision strengths. We examine the advantages and shortcomings of each of these methods, and estimate rate uncertainties from the resulting data dispersion. We proceed to construct excitation balance spectral models, and compare the predictions from each data set with observed spectra from various astronomical objects. We are thus able to establish benchmarks in the spectral modeling of [Fe II] emission in the IR and optical regions as well as in the UV Fe II absorption spectra. Finally, we provide diagnostic line ratios and line emissivities for emission spectroscopy as well as column densities for absorption spectroscopy. All atomic data and models are available online and through the AtomPy atomic data curation environment.
[59]  [pdf] - 1043201
Atomic data for Zn II - Improving Spectral Diagnostics of Chemical Evolution in High-redshift Galaxies
Comments: accepted The Astrophysical Journal
Submitted: 2015-04-07
Damped Lyman-alpha (DLA) and sub-DLA absorbers in quasar spectra provide the most sensitive tools for measuring element abundances of distant galaxies. Estimation of abundances from absorption lines depends sensitively on the accuracy of the atomic data used. We have started a project to produce new atomic spectroscopic parameters for optical/UV spectral lines using state-of-the-art computer codes employing very broad configuration interaction basis. Here we report our results for Zn II, an ion used widely in studies of the interstellar medium (ISM) as well as DLA/sub-DLAs. We report new calculations of many energy levels of Zn II, and the line strengths of the resulting radiative transitions. Our calculations use the configuration interaction approach within a numerical Hartree-Fock framework. We use both non-relativistic and quasi-relativistic one-electron radial orbitals. We have incorporated the results of these atomic calculations into the plasma simulation code Cloudy, and applied them to a lab plasma and examples of a DLA and a sub-DLA. Our values of the Zn II {\lambda}{\lambda} 2026, 2062 oscillator strengths are higher than previous values by 0.10 dex. Cloudy calculations for representative absorbers with the revised Zn atomic data imply ionization corrections lower than calculated before by 0.05 dex. The new results imply Zn metallicities should be lower by 0.1 dex for DLAs and by 0.13-0.15 dex for sub-DLAs than in past studies. Our results can be applied to other studies of Zn II in the Galactic and extragalactic ISM.
[60]  [pdf] - 1043170
Bootstrapping dielectronic recombination from second-row elements and the Orion Nebula
Comments: Ap.J. At Press
Submitted: 2015-03-19
Dielectronic recombination (DR) is the dominant recombination process for most heavy elements in photoionized clouds. Accurate DR rates for a species can be predicted when the positions of autoionizing states are known. Unfortunately such data are not available for most third and higher-row elements. This introduces an uncertainty that is especially acute for photoionized clouds, where the low temperatures mean that DR occurs energetically through very low-lying autoionizing states. This paper discusses S$^{2+} \rightarrow$ S$^+$ DR, the process that is largely responsible for establishing the [S~III]/[S~II] ratio in nebulae. We derive an empirical rate coefficient using a novel method for second-row ions, which do have accurate data. Photoionization models are used to reproduce the [O~III] / [O~II] / [O~I] / [Ne~III] intensity ratios in central regions of the Orion Nebula. O and Ne have accurate atomic data and can be used to derive an empirical S$^{2+} \rightarrow$ S$^+$ DR rate coefficient at $\sim 10^{4}$~K. We present new calculations of the DR rate coefficient for S$^{2+} \rightarrow$ S$^+$ and quantify how uncertainties in the autoionizing level positions affect it. The empirical and theoretical results are combined and we derive a simple fit to the resulting rate coefficient at all temperatures for incorporation into spectral synthesis codes. This method can be used to derive empirical DR rates for other ions, provided that good observations of several stages of ionization of O and Ne are available.
[61]  [pdf] - 972606
Accurate determination of the free-free Gaunt factor. II - relativistic Gaunt factors
Comments: 8 pages, 8 figures, 2 tables
Submitted: 2015-02-24
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.
[62]  [pdf] - 1223992
Variability of the coronal line region in NGC 4151
Comments: 12 pages, 9 figures; accepted by MNRAS
Submitted: 2015-01-13
We present the first extensive study of the coronal line variability in an active galaxy. Our data set for the nearby source NGC 4151 consists of six epochs of quasi-simultaneous optical and near-infrared spectroscopy spanning a period of about eight years and five epochs of X-ray spectroscopy overlapping in time with it. None of the coronal lines showed the variability behaviour observed for the broad emission lines and hot dust emission. In general, the coronal lines varied only weakly, if at all. Using the optical [Fe VII] and X-ray O VII emission lines we estimate that the coronal line gas has a relatively low density of n~10^3 cm^-3 and a relatively high ionisation parameter of log U~1. The resultant distance of the coronal line gas from the ionising source is about two light years, which puts this region well beyond the hot inner face of the obscuring dusty torus. The high ionisation parameter implies that the coronal line region is an independent entity rather than part of a continuous gas distribution connecting the broad and narrow emission line regions. We present tentative evidence for the X-ray heated wind scenario of Pier & Voit. We find that the increased ionising radiation that heats the dusty torus also increases the cooling efficiency of the coronal line gas, most likely due to a stronger adiabatic expansion.
[63]  [pdf] - 1319517
Collisional excitation of [C II], [O I] and CO in Massive Galaxies
Comments: 17 pages, submitted to MNRAS
Submitted: 2015-01-06
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 \mu m , [O I]63 \mu 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.
[64]  [pdf] - 903870
Implications of Coronal Line Emission in NGC 4696
Comments: 11 pages, 8 figures, 6 tables; accepted for publication in MNRAS
Submitted: 2014-10-16
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] {\lambda}5303 / [Fe X] {\lambda}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] {\lambda}6375, but not [Fe XIV] {\lambda}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.
[65]  [pdf] - 1216496
Line-driven radiative outflows in luminous quasars
Comments: 19 pages, 14 figures, 5 tables, accepted to MNRAS (updated to include correct Fig. 6 and 8)
Submitted: 2014-08-25, last modified: 2014-08-28
An analysis of ~19500 narrow (<200 km/s) CIV 1548.2,1550.8 absorbers in ~34000 Sloan Digital Sky Survey quasar spectra is presented. The statistics of the number of absorbers as a function of outflow-velocity shows that in approximately two-thirds of outflows, with multiple CIV absorbers present, absorbers are line-locked at the 500 km/s velocity separation of the CIV absorber doublet; appearing as 'triplets' in the quasar spectra. Line-locking is an observational signature of radiative line driving in outflowing material, where the successive shielding of 'clouds' of material in the outflow locks the clouds together in outflow velocity. Line-locked absorbers are seen in both broad absorption line quasars (BALs) and non-BAL quasars with comparable frequencies and with velocities out to at least 20000 km/s. There are no detectable differences in the absorber properties and the dust content of single CIV doublets and line-locked CIV doublets. The gas associated with both single and line-locked CIV absorption systems includes material with a wide range of ionization potential (14-138 eV). Both single and line-locked CIV absorber systems show strong systematic trends in their ionization as a function of outflow velocity, with ionization decreasing rapidly with increasing outflow velocity. Initial simulations, employing CLOUDY, demonstrate that a rich spectrum of line-locked signals at various velocities may be expected due to significant opacities from resonance lines of Li-, He- and H-like ions of O, C and N, along with contributions from HeII and HI resonance lines. The simulations confirm that line driving can be the dominant acceleration mechanism for clouds with N(HI) ~ 10^19 cm^-2.
[66]  [pdf] - 863133
Accurate determination of the free-free Gaunt factor; I - non-relativistic Gaunt factors
Comments: 10 pages, 5 tables, 3 figures. Fixed typo in Eq. 14
Submitted: 2014-07-18, last modified: 2014-08-17
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.
[67]  [pdf] - 831976
On the Observability of Optically Thin Coronal Hyperfine Structure Lines
Comments: 48 pages; 13 figures; 9 tables; additional images available at the online journal
Submitted: 2014-06-03
We present Cloudy calculations for the intensity of coronal hyperfine lines in various environments. We model indirect collisional and radiative transitions, and quantify the collisionally-excited line emissivity in the density-temperature phase-space. As an observational aid, we also express the emissivity in units of the continuum in the 0.4--0.7 keV band. For most hyperfine lines, knowledge of the X-ray surface brightness and the plasma temperature is sufficient for rough estimates. We find that the radiation fields of both Perseus A and Virgo A can enhance the populations of highly ionized species within 1 kpc. They can also enhance line emissivity within the cluster core. This could have implications for the interpretation of spectra around bright AGN. We find the intensity of the $^{57}$Fe XXIV {\lambda}3.068 mm to be about two orders of magnitude fainter than previously thought, at about 20 {\mu}K. Comparably bright lines may be found in the infrared. Finally, we find the intensity of hyperfine lines in the Extended Orion Nebula to be low, due to the shallow sightline. Observations of coronal hyperfine lines will likely be feasible with the next generation of radio and sub-mm telescopes.
[68]  [pdf] - 821717
Radiative Cooling II: Effects of Density and Metallicity
Comments: Published by MNRAS, online tables can be found at
Submitted: 2014-03-12, last modified: 2014-05-08
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.
[69]  [pdf] - 1180630
Atomic data for S II - Toward Better Diagnostics of Chemical Evolution in High-redshift Galaxies
Comments: Accepted for publication in The Astrophysical Journal; 34 pages, 10 figures
Submitted: 2013-11-11
Absorption-line spectroscopy is a powerful tool used to estimate element abundances in the nearby as well as distant universe. The accuracy of the abundances thus derived is, naturally, limited by the accuracy of the atomic data assumed for the spectral lines. We have recently started a project to perform the new extensive atomic data calculations used for optical/UV spectral lines in the plasma modeling code Cloudy using state-of-the-art quantal calculations. Here we demonstrate our approach by focussing on S II, an ion used to estimate metallicities for Milky Way interstellar clouds as well as distant damped Lyman-alpha (DLA) and sub-DLA absorber galaxies detected in the spectra of quasars and gamma-ray bursts (GRBs). We report new extensive calculations of a large number of energy levels of S II, and the line strengths of the resulting radiative transitions. Our calculations are based on the configuration interaction approach within a numerical Hartree-Fock framework, and utilize both non-ralativistic and quasirelativistic one-electron radial orbitals. The results of these new atomic calculations are then incorporated into Cloudy and applied to a lab plasma, and a typical DLA, for illustrative purposes. The new results imply relatively modest changes (~0.04 dex) to the metallicities estimated from S II in past studies. These results will be readily applicable to other studies of S II in the Milky Way and other galaxies.
[70]  [pdf] - 753970
Effects of External Radiation Fields on Line Emission - Application to Star-forming Regions
Comments: 12 pages, 7 figures, accepted for publication to ApJ
Submitted: 2013-10-29
A variety of astronomical environments contain clouds irradiated by a combination of isotropic and beamed radiation fields. For example, molecular clouds may be irradiated by the isotropic cosmic microwave background (CMB), as well as by a nearby active galactic nucleus (AGN). These radiation fields excite atoms and molecules and produce emission in different ways. We revisit the escape probability theorem and derive a novel expression that accounts for the presence of external radiation fields. We show that when the field is isotropic the escape probability is reduced relative to that in the absence of external radiation. This is in agreement with previous results obtained under ad hoc assumptions or with the two-level system, but can be applied to complex many-level models of atoms or molecules. This treatment is in the development version of the spectral synthesis code Cloudy. We examine the spectrum of a Spitzer cloud embedded in the local interstellar radiation field, and show that about 60 percent of its emission lines are sensitive to background subtraction. We argue that this geometric approach could provide an additional tool toward understanding the complex radiation fields of starburst galaxies.
[71]  [pdf] - 1180192
Interpreting the Ionization Sequence in AGN Emission-Line Spectra
Comments: 57 pages; 18 figures, accepted by MNRAS
Submitted: 2013-10-23
We investigate the physical cause of the great range in the ionization level seen in the spectra of narrow lined active galactic nuclei (AGN). Mean field independent component analysis identifies examples of individual SDSS galaxies whose spectra are not dominated by emission due to star formation (SF), which we designate as AGN. We assembled high S/N ratio composite spectra of a sequence of these AGN defined by the ionization level of their narrow-line regions (NLR), extending down to very low-ionization cases. We used a local optimally emitting cloud (LOC) model to fit emission-line ratios in this AGN sequence. These included the weak lines that can be measured only in the co-added spectra, providing consistency checks on strong line diagnostics. After integrating over a wide range of radii and densities our models indicate that the radial extent of the NLR is the major parameter in determining the position of high to moderate ionization AGN along our sequence, providing a physical interpretation for their systematic variation. Higher ionization AGN contain optimally emitting clouds that are more concentrated towards the central continuum source than in lower ionization AGN. Our LOC models indicate that for the objects that lie on our AGN sequence, the ionizing luminosity is anticorrelated with the NLR ionization level, and hence anticorrelated with the radial concentration and physical extent of the NLR. A possible interpretation that deserves further exploration is that the ionization sequence might be an age sequence where low ionization objects are older and have systematically cleared out their central regions by radiation pressure. We consider that our AGN sequence instead represents a mixing curve of SF and AGN spectra, but argue that while many galaxies do have this type of composite spectra, our AGN sequence appears to be a special set of objects with negligible SF excitation.
[72]  [pdf] - 1180003
Cold gas dynamics in Hydra-A: evidence for a rotating disk
Comments: Accepted for publication in MNRAS on October 9 2013
Submitted: 2013-10-16
We present multi-frequency observations of the radio galaxy Hydra-A (3C218) located in the core of a massive, X-ray luminous galaxy cluster. IFU spectroscopy is used to trace the kinematics of the ionised and warm molecular hydrogen which are consistent with a ~ 5 kpc rotating disc. Broad, double-peaked lines of CO(2-1), [CII]157 $\mu$m and [OI]63 $\mu$m are detected. We estimate the mass of the cold gas within the disc to be M$_{gas}$ = 2.3 $\pm$ 0.3 x 10$^9$ M$_{\odot}$. These observations demonstrate that the complex line profiles found in the cold atomic and molecular gas are related to the rotating disc or ring of gas. Finally, an HST image of the galaxy shows that this gas disc contains a substantial mass of dust. The large gas mass, SFR and kinematics are consistent with the levels of gas cooling from the ICM. We conclude that the cold gas originates from the continual quiescent accumulation of cooled ICM gas. The rotation is in a plane perpendicular to the projected orientation of the radio jets and ICM cavities hinting at a possible connection between the kpc-scale cooling gas and the accretion of material onto the black hole. We discuss the implications of these observations for models of cold accretion, AGN feedback and cooling flows.
[73]  [pdf] - 1172680
Long Term Spectral Evolution of Tidal Disruption Candidates Selected by Strong Coronal Lines
Comments: 30 Pages, 10 Figures, 2 Tables, Accepted for Publication in ApJ
Submitted: 2013-07-11
We present results of follow-up optical spectroscopic observations of seven rare, extreme coronal line emitting galaxies reported by Wang et al. (2012) with Multi-Mirror Telescope (MMT). Large variations in coronal lines are found in four objects, making them strong candidates of tidal disruption events (TDE). For the four TDE candidates, all the coronal lines with ionization status higher than [Fe VII] disappear within 5-9 years. The [Fe VII] faded by a factor of about five in one object (J0952+2143) within 4 years, whereas emerged in other two without them previously. A strong increment in the [O III] flux is observed, shifting the line ratios towards the loci of active galactic nucleus on the BPT diagrams. Surprisingly, we detect a non-canonical [O III]5007/[O III]4959 2 in two objects, indicating a large column density of O$^{2+}$ and thus probably optical thick gas. This also requires a very large ionization parameter and relatively soft ionizing spectral energy distribution (e.g. blackbody with $T < 5\times 10^4$ K). Our observations can be explained as echoing of a strong ultraviolet to soft X-ray flare caused by tidal disruption events, on molecular clouds in the inner parsecs of the galactic nuclei. Re-analyzing the SDSS spectra reveals double-peaked or strongly blue-shouldered broad lines in three of the objects, which disappeared in the MMT spectra in two objects, and faded by a factor of ten in 8 years in the remaining object with a decrease in both the line width and centroid offset. We interpret these broad lines as arising from decelerating biconical outflows. Our results demonstrate that the signatures of echoing can persist for as long as ten years, and can be used to probe the gas environment in the quiescent galactic nuclei.
[74]  [pdf] - 1172469
Detecting the rapidly expanding outer shell of the Crab Nebula: where to look
Comments: Accepted by ApJ
Submitted: 2013-07-03
We present a range of steady-state photoionization simulations, corresponding to different assumed shell geometries and compositions, of the unseen postulated rapidly expanding outer shell to the Crab Nebula. The properties of the shell are constrained by the mass that must lie within it, and by limits to the intensities of hydrogen recombination lines. In all cases the photoionization models predict very strong emission from high ionization lines that will not be emitted by the Crab's filaments, alleviating problems with detecting these lines in the presence of light scattered from brighter parts of the Crab. The NIR [Ne VI] $\lambda$7.652 $\mu$m line is a particularly good case; it should be dramatically brighter than the optical lines commonly used in searches. The C IV $\lambda1549\AA$ doublet is predicted to be the strongest absorption line from the shell, which is in agreement with HST observations. We show that the cooling timescale for the outer shell is much longer than the age of the Crab, due to the low density. This means that the temperature of the shell will actually "remember" its initial conditions. However, the recombination time is much shorter than the age of the Crab, so the predicted level of ionization should approximate the real ionization. In any case, it is clear that IR observations present the best opportunity to detect the outer shell and so guide future models that will constrain early events in the original explosion.
[75]  [pdf] - 1166087
Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--III:Tangential Motions using AstroDrizzle Images
Submitted: 2013-04-18
We have been able to compare with astrometric precision AstroDrizzle processed images of NGC 6720 (the Ring Nebula) made using two cameras on the Hubble Space Telescope. The time difference of the observations was 12.925 yrs. This large time-base allowed determination of tangential velocities of features within this classic planetary nebula. Individual features were measured in [N II] images as were the dark knots seen in silhouette against background nebular [O III] emission. An image magnification and matching technique was also used to test the accuracy of the usual assumption of homologous expansion. We found that homologous expansion does apply, but the rate of expansion is greater along the major axis of the nebula, which is intrinsically larger than the minor axis. We find that the dark knots expand more slowly that the nebular gas, that the distance to the nebula is 720 pc +/-30%, and the dynamic age of the Ring Nebula is about 4000 yrs. The dynamic age is in agreement with the position of the central star on theoretical curves for stars collapsing from the peak of the Asymptotic Giant Branch to being white dwarfs.
[76]  [pdf] - 1165443
Corrigendum: Improved He~I Emissivities in the Case B Approximation
Comments: Accepted to MNRAS Letters. 2 pages, 1 figure, 1 table, 1 supplemental table
Submitted: 2013-03-20
Errors were discovered in our recent calculations of HeI emissivities (Porter et al. 2012; arXiv:1206.4115). This corrigendum discusses the errors and reports corrected emissivities.
[77]  [pdf] - 1165141
Suppression of Dielectronic Recombination Due to Finite Density Effects
Comments: main article (27 pages, 7 figures); supporting information (37 pages, 46 figures)
Submitted: 2013-03-10
We have developed a general model for determining density-dependent effective dielectronic recombination (DR) rate coefficients in order to explore finite-density effects on the ionization balance of plasmas. Our model consists of multiplying by a suppression factor those highly-accurate total zero-density DR rate coefficients which have been produced from state-of-the-art theoretical calculations and which have been benchmarked by experiment. The suppression factor is based-upon earlier detailed collision-radiative calculations which were made for a wide range of ions at various densities and temperatures, but used a simplified treatment of DR. A general suppression formula is then developed as a function of isoelectronic sequence, charge, density, and temperature. These density-dependent effective DR rate coefficients are then used in the plasma simulation code Cloudy to compute ionization balance curves for both collisionally ionized and photoionized plasmas at very low (ne = 1 cm^-3) and finite (ne=10^10 cm^-3) densities. We find that the denser case is significantly more ionized due to suppression of DR, warranting further studies of density effects on DR by detailed collisional-radiative calculations which utilize state-of-the-art partial DR rate coefficients. This is expected to impact the predictions of the ionization balance in denser cosmic gases such as those found in nova and supernova shells, accretion disks, and the broad emission line regions in active galactic nuclei.
[78]  [pdf] - 1508116
Expanded Iron UTA spectra -- probing the thermal stability limits in AGN clouds
Comments: ApJ in press
Submitted: 2013-03-01
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.
[79]  [pdf] - 628502
The 2013 Release of Cloudy
Comments: Major review article to be published in RevMexAA. 27 pages, 18 figures
Submitted: 2013-02-18
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.
[80]  [pdf] - 1508114
Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--II:Physical Conditions
Comments: Accepted for publication in the Astronomical Journal
Submitted: 2013-01-28
We have performed a detailed analysis of the electron temperature and density in the the Ring Nebula using the calibrated HST WFC3 images described in the preceding paper. The electron temperature (Te) determined from [N II] and [O III] rises slightly and monotonically towards the central star. The observed equivalent width (EW) in the central region indicates that Te rises as high as 13000 K. In contrast, the low EW's in the outer regions are largely due to scattered diffuse Galactic radiation by dust. The images allowed determination of unprecedented small scale variations in Te. These variations indicate that the mean square area temperature fluctuations are significantly higher than expected from simple photoionization. The power producing these fluctuations occurs at scales of less than 3.5E15 cm. This scale length provides a strong restriction on the mechanism causing the large t^2 values observed.
[81]  [pdf] - 1159325
Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--I: Structure and Evolution
Comments: Accepted for publication by the Astronomical Journal
Submitted: 2013-01-28
We have performed a detailed analysis of the Ring Nebula (NGC 6720) using HST WFC3 images and derived a new 3-D model. Existing high spectral resolution spectra played an important supplementary role in our modeling. It is shown that the Main Ring of the nebula is an ionization-bounded irregular non-symmetric disk with a central cavity and perpendicular extended lobes pointed almost towards the observer. The faint outer halos are determined to be fossil radiation, i.e. radiation from gas ionized in an earlier stage of the nebula when it was not ionization bounded. The narrow-band WFC3 filters that isolate some of the emission-lines are affected by broadening on their short wavelength side and all the filters were calibrated using ground-based spectra. The filter calibration results are presented in an appendix.
[82]  [pdf] - 1159271
Classification and analysis of emission-line galaxies using mean field independent component analysis
Comments: MNRAS accepted. 29 pages, 24 figures, 3 tables
Submitted: 2013-01-24
We present an analysis of the optical spectra of narrow emission-line galaxies, based on mean field independent component analysis (MFICA). Samples of galaxies were drawn from the Sloan Digital Sky Survey (SDSS) and used to generate compact sets of `continuum' and `emission-line' component spectra. These components can be linearly combined to reconstruct the observed spectra of a wider sample of galaxies. Only 10 components - five continuum and five emission line - are required to produce accurate reconstructions of essentially all narrow emission-line galaxies; the median absolute deviations of the reconstructed emission-line fluxes, given the signal-to-noise ratio (S/N) of the observed spectra, are 1.2-1.8 sigma for the strong lines. After applying the MFICA components to a large sample of SDSS galaxies we identify the regions of parameter space that correspond to pure star formation and pure active galactic nucleus (AGN) emission-line spectra, and produce high S/N reconstructions of these spectra. The physical properties of the pure star formation and pure AGN spectra are investigated by means of a series of photoionization models, exploiting the faint emission lines that can be measured in the reconstructions. We are able to recreate the emission line strengths of the most extreme AGN case by assuming the central engine illuminates a large number of individual clouds with radial distance and density distributions, f(r) ~ r^gamma and g(n) ~ n^beta, respectively. The best fit is obtained with gamma = -0.75 and beta = -1.4. From the reconstructed star formation spectra we are able to estimate the starburst ages. These preliminary investigations serve to demonstrate the success of the MFICA-based technique in identifying distinct emission sources, and its potential as a tool for the detailed analysis of the physical properties of galaxies in large-scale surveys.
[83]  [pdf] - 1158674
The Nature of the H2-Emitting Gas in the Crab Nebula
Comments: 51 pages, 15 figures, accepted for publication in MNRAS, revised Figure 12 results unchanged
Submitted: 2012-12-21, last modified: 2013-01-01
Understanding how molecules and dust might have formed within a rapidly expanding young supernova remnant is important because of the obvious application to vigorous supernova activity at very high redshift. In previous papers, we found that the H2 emission is often quite strong, correlates with optical low-ionization emission lines, and has a surprisingly high excitation temperature. Here we study Knot 51, a representative, bright example, for which we have available long slit optical and NIR spectra covering emission lines from ionized, neutral, and molecular gas, as well as HST visible and SOAR Telescope NIR narrow-band images. We present a series of CLOUDY simulations to probe the excitation mechanisms, formation processes and dust content in environments that can produce the observed H2 emission. We do not try for an exact match between model and observations given Knot 51's ambiguous geometry. Rather, we aim to explain how the bright H2 emission lines can be formed from within the volume of Knot 51 that also produces the observed optical emission from ionized and neutral gas. Our models that are powered only by the Crab's synchrotron radiation are ruled out because they cannot reproduce the strong, thermal H2 emission. The simulations that come closest to fitting the observations have the core of Knot 51 almost entirely atomic with the H2 emission coming from just a trace molecular component, and in which there is extra heating. In this unusual environment, H2 forms primarily by associative detachment rather than grain catalysis. In this picture, the 55 H2-emitting cores that we have previously catalogued in the Crab have a total mass of about 0.1 M_sun, which is about 5% of the total mass of the system of filaments. We also explore the effect of varying the dust abundance. We discuss possible future observations that could further elucidate the nature of these H2 knots.
[84]  [pdf] - 1158313
Radiative cooling in collisionally and photo ionized plasmas
Submitted: 2012-12-05
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.
[85]  [pdf] - 1157772
A transition mass for black holes to show broad emission lines
Comments: 7 pages, 5 figures. Submitted to MNRAS
Submitted: 2012-11-14
Although the super-massive (AGN) and stellar mass (XRBs) black holes have many properties in common, the broad emission lines (BELs) are exclusively signatures of the AGN. Based on the detection of these lines from SDSS data bases, there seems to be no AGN with mass M_BH \lesssim 10^5 M_sun. In this paper we investigate if such low mass black holes are really non-existent or they are undetected because the BELs in them are not produced efficiently. Using the ionizing spectral energy distribution for a wide range of black hole mass, 10 - 10^9 M_sun, spanning XRBs to AGN, we calculate the equivalent widths (EWs) of ultraviolet and optical lines Ly\alpha 1216 \AA, H\beta 4861 \AA, CIV 1549 \AA and MgII 2798 \AA. The LOC (locally optimally emitting cloud) model has been used to describe the broad emission line region (BELR) for the calculations. We find that the hardening of the SED shape with decreasing mass do not decrease the BEL EWs. However, finite size of the BELR, as measured by the line widths, which is controlled by the mass of the black hole, regulates the production of these emission lines. There seems to be a peak in the EWs of the emission lines for typical AGN black holes of ~ 10^8 M_sun, below which the lines become intrinsically fainter with a sharp fall-off below ~ 10^6 M_sun. This may be the cause of the absence of low mass AGN in SDSS.
[86]  [pdf] - 1152149
Two-Component Structure of the Hbeta Broad-Line Region in Quasars. I. Evidence from Spectral Principal Component Analysis
Comments: 22 pages, 17 figures, accepted for publication in The Astrophysical Journal
Submitted: 2012-10-15
We report on a spectral principal component analysis (SPCA) of a sample of 816 quasars, selected to have small Fe II velocity shifts with spectral coverage in the rest wavelength range 3500--5500 \AA. The sample is explicitly designed to mitigate spurious effects on SPCA induced by Fe II velocity shifts. We improve the algorithm of SPCA in the literature and introduce a new quantity, \emph{the fractional-contribution spectrum}, that effectively identifies the emission features encoded in each eigenspectrum. The first eigenspectrum clearly records the power-law continuum and very broad Balmer emission lines. Narrow emission lines dominate the second eigenspectrum. The third eigenspectrum represents the Fe II emission and a component of the Balmer lines with kinematically similar intermediate velocity widths. Correlations between the weights of the eigenspectra and parametric measurements of line strength and continuum slope confirm the above interpretation for the eigenspectra. Monte Carlo simulations demonstrate the validity of our method to recognize cross talk in SPCA and firmly rule out a single-component model for broad Hbeta. We also present the results of SPCA for four other samples that contain quasars in bins of larger Fe II velocity shift; similar eigenspectra are obtained. We propose that the Hbeta-emitting region has two kinematically distinct components: one with very large velocities whose strength correlates with the continuum shape, and another with more modest, intermediate velocities that is closely coupled to the gas that gives rise to Fe II emission.
[87]  [pdf] - 1151032
Pumping up the [N I] nebular lines
Comments: ApJ in press. 8 pages of main paper plus 11 pages of appendices, with 13 figures and 12 tables
Submitted: 2012-08-31
The optical [N I] doublet near 5200 {\AA} 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{\beta} 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.
[88]  [pdf] - 1150594
Herschel observations of extended atomic gas in the core of the Perseus cluster
Comments: 21 pages, 16 figures and 10 tables. Accepted for publication in MNRAS
Submitted: 2012-08-08
We present Herschel observations of the core of the Perseus cluster of galaxies. The brightest cluster galaxy, NGC 1275, is surrounded by a network of filaments previously imaged extensively in H{\alpha} and CO. In this work, we report detections of FIR lines with Herschel. All but one of the lines are spatially extended, with the [CII] line emission extending up to 25 kpc from the core. There is spatial and kinematical correlation among [CII], H{\alpha} and CO, which gives us confidence to model the different components of the gas with a common heating model. With the help of FIR continuum Herschel measurements, together with a suite of coeval radio, submm and infrared data, we performed a SED fitting of NGC 1275 using a model that contains contributions from dust emission as well as synchrotron AGN emission. The data indicate a low dust emissivity index, beta ~ 1, a total dust mass close to 10^7 solar mass, a cold dust component with temperature 38 \pm 2 K and a warm dust component with temperature of 116 \pm 9 K. The FIR-derived star formation rate (SFR) is 24 \pm 1 solar mass per yr, in close agreement with the FUV-derived SFR. We investigated in detail the source of the Herschel FIR and H{\alpha} emissions emerging from a core region 4 kpc in radius. Based on simulations conducted using the radiative transfer code, CLOUDY, a heating model comprising old and young stellar populations is sufficient to explain these observations. We have also detected [CII] in three well-studied regions of the filaments. We find a [OI]/[CII] ratio about 1 dex smaller than predicted by the otherwise functional Ferland (2009) model. The line ratio suggests that the lines are optically thick, as is typical of galactic PDRs, and implies that there is a large reservoir of cold atomic gas. [abridged]
[89]  [pdf] - 1124436
Hydrogen Two-Photon Continuum Emission from the Horseshoe Filament in NGC 1275
Comments: 9 pages, 8 figures. Accepted for publication in MNRAS
Submitted: 2012-06-27
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.
[90]  [pdf] - 1124219
Improved He I Emissivities in the Case B Approximation
Comments: Accepted to MNRAS Letters; 4 pages, 4 figures, 2 tables, 1 supplemental file
Submitted: 2012-06-18
We update our prior work on the case B collisional-recombination spectrum of He I to incorporate \textit{ab initio} photoionisation cross-sections. This large set of accurate, self-consistent cross-sections represents a significant improvement in He I emissivity calculations because it largely obviates the piecemeal nature that has marked all modern works. A second, more recent set of \textit{ab initio} cross-sections is also available, but we show that those are less consistent with bound-bound transition probabilities than our adopted set. We compare our new effective recombination coefficients with our prior work and our new emissivities with those by other researchers, and we conclude with brief remarks on the effects of the present work on the He I error budget. Our calculations cover temperatures $5000 \le T_e \le 25000$ K and densities $10^1 \le n_e \le 10^{14}$ cm$^{-3}$. Full results are available online.
[91]  [pdf] - 545510
Intervening BLR Clouds' Effects on Optical/UV Spectrum
Comments: 22 pages, 11 figures
Submitted: 2012-05-25
Recent x-ray observations of Mrk 766 suggest that broad emission line region clouds cross our line of sight and produce variable x-ray absorption. Here we investigate what optical/ultraviolet spectroscopic features would be produced by such "Intervening BLR Clouds" (IBC) crossing our line of sight to the accretion disk, the source of the optical/UV continuum. Although the emission spectrum produced by intervening clouds is identical to the standard BLR model, they may produce absorption features on the optical or UV continuum. Single clouds will have little effect on the optical/UV spectrum because BLR clouds are likely to be much smaller than the accretion disk. This is unlike the X-ray case, where the radiation source is considerably smaller. However, an ensemble of intervening clouds will produce spectroscopic features in the FUV including a strong depression between the Lyman limit and Ly{\alpha}. The amount of the depression will indicate the line-of-sight covering factor of clouds, an unknown quantity that is important for the ionization of the intergalactic medium and the energy budget of AGN. Comparison with observations suggests that the SED of Mrk 766 may be affected by intervening BLR clouds and IBC may exist in most of AGNs.
[92]  [pdf] - 574508
Multiphase Signatures of AGN Feedback in Abell 2597
Comments: 16 pages, 11 figures; Accepted for publication in MNRAS
Submitted: 2012-05-10
We present new Chandra X-ray observations of the brightest cluster galaxy (BCG) in the cool core cluster Abell 2597. The data reveal an extensive kpc-scale X-ray cavity network as well as a 15 kpc filament of soft-excess gas exhibiting strong spatial correlation with archival VLA radio data. In addition to several possible scenarios, multiwavelength evidence may suggest that the filament is associated with multiphase (10^3 - 10^7 K) gas that has been entrained and dredged-up by the propagating radio source. Stemming from a full spectral analysis, we also present profiles and 2D spectral maps of modeled X-ray temperature, entropy, pressure, and metal abundance. The maps reveal an arc of hot gas which in projection borders the inner edge of a large X-ray cavity. Although limited by strong caveats, we suggest that the hot arc may be (a) due to a compressed rim of cold gas pushed outward by the radio bubble or (b) morphologically and energetically consistent with cavity-driven active galactic nucleus (AGN) heating models invoked to quench cooling flows, in which the enthalpy of a buoyant X-ray cavity is locally thermalized as ambient gas rushes to refill its wake. If confirmed, this would be the first observational evidence for this model.
[93]  [pdf] - 574507
Residual Cooling and Persistent Star Formation amid AGN Feedback in Abell 2597
Comments: 20 pages, 7 figures; Accepted for publication in MNRAS
Submitted: 2012-05-10
New Chandra X-ray and Herschel FIR observations enable a multiwavelength study of active galactic nucleus (AGN) heating and intracluster medium (ICM) cooling in the brightest cluster galaxy of Abell 2597. The new Chandra observations reveal the central < 30 kiloparsec X-ray cavity network to be more extensive than previously thought, and associated with enough enthalpy to theoretically inhibit the inferred classical cooling flow. Nevertheless, we present new evidence, consistent with previous results, that a moderately strong residual cooling flow is persisting at 4%-8% of the classically predicted rates in a spatially structured manner amid the feedback-driven excavation of the X-ray cavity network. New Herschel observations are used to estimate warm and cold dust masses, a lower-limit gas-to-dust ratio, and a star formation rate consistent with previous measurements. The cooling time profile of the ambient X-ray atmosphere is used to map the locations of the observational star formation entropy threshold as well as the theoretical thermal instability threshold. Both lie just outside the < 30 kpc central region permeated by X-ray cavities, and star formation as well as ionized and molecular gas lie interior to both. The young stars are distributed in an elongated region that is aligned with the radio lobes, and their estimated ages are both younger and older than the X-ray cavity network, suggesting both jet-triggered as well as persistent star formation over the current AGN feedback episode. Bright X-ray knots that are coincident with extended Ly-alpha and FUV continuum filaments motivate a discussion of structured cooling from the ambient hot atmosphere along a projected axis that is perpendicular to X-ray cavity and radio axis. We conclude that the cooling ICM is the dominant contributor of the cold gas reservoir fueling star formation and AGN activity in the Abell 2597 BCG.
[94]  [pdf] - 1116459
The Curious Conundrum Regarding Sulfur Abundances In Planetary Nebulae
Comments: 40 pages, 14 figures, 3 tables. Accepted for publication in the Astrophysical Journal
Submitted: 2012-02-07
Sulfur abundances derived from optical emission line measurements and ionization correction factors in planetary nebulae are systematically lower than expected for the objects' metallicities. We have carefully considered a large range of explanations for this "sulfur anomaly", including: (1) correlations between the size of the sulfur deficit and numerous nebular and central star properties; (2) ionization correction factors which under-correct for unobserved ions; (3) effects of dielectronic recombination on the sulfur ionization balance; (4) sequestering of S into dust and/or molecules; and (5) excessive destruction of S or production of O by AGB stars. It appears that all but the second scenario can be ruled out. However, we find evidence that the sulfur deficit is generally reduced but not eliminated when S^+3 abundances determined directly from IR measurements are used in place of the customary sulfur ionization correction factor. We tentatively conclude that the sulfur anomaly is caused by the inability of commonly used ICFs to properly correct for populations of ionization stages higher than S^+2.
[95]  [pdf] - 1116294
Star formation in self-gravitating disks in active galactic nuclei. II. Episodic formation of broad line regions
Comments: 27 pages in emulateapj style, 10 figures, Accepted by ApJ 2012
Submitted: 2012-01-31
(abridged) We study the consequence of star formation (SF) in an self-gravity dominated accretion disk in quasars. The warm skins of the SF disk are governed by the radiation from the inner part of the accretion disk to form Compton atmosphere (CAS). The CAS are undergoing four phases to form broad line regions. Phase I is the duration of pure accumulation supplied by the SF disk. During phase II clouds begin to form due to line cooling and sink to the SF disk. Phase III is a period of preventing clouds from sinking to the SF disk through dynamic interaction between clouds and the CAS. Finally, phase IV is an inevitable collapse of the entire CAS through line cooling. This CAS evolution drives the episodic appearance of BLRs. Geometry and dynamics of BLRs can be self-consistently derived from the thermal instability of the CAS during phases II and III by linear analysis. The metallicity gradient of SF disk gives rise to different properties of clouds from outer to inner part of BLRs. We find that clouds have column density N_H < 10^22cm^{-2} in the metal-rich regions whereas they have N_H > 10^22 cm^{-2} in the metal-poor regions. The metal-rich clouds compose the high ionization line (HIL) regions whereas the metal-poor clouds are in low ionization line (LIL) regions. Metal-rich clouds in HIL regions will be blown away by radiation pressure, forming the observed outflows. The LIL regions are episodic due to the mass cycle of clouds with the CAS in response to continuous injection by the SF disk, giving rise to different types of AGNs. Based on SDSS quasar spectra, we identify a spectral sequence in light of emission line equivalent width from Phase I to IV. A key phase in the episodic appearance of the BLRs is bright type II AGNs with no or only weak BLRs. We discuss observational implications and tests of the theoretical predictions of this model.
[96]  [pdf] - 486726
The Impact of Recent Advances in Laboratory Astrophysics on our Understanding of the Cosmos
Comments: 61 pages. 11 figures; to appear in Reports on Progress in Physics
Submitted: 2011-12-09
An emerging theme in modern astrophysics is the connection between astronomical observations and the underlying physical phenomena that drive our cosmos. Both the mechanisms responsible for the observed astrophysical phenomena and the tools used to probe such phenomena - the radiation and particle spectra we observe - have their roots in atomic, molecular, condensed matter, plasma, nuclear and particle physics. Chemistry is implicitly included in both molecular and condensed matter physics. This connection is the theme of the present report, which provides a broad, though non-exhaustive, overview of progress in our understanding of the cosmos resulting from recent theoretical and experimental advances in what is commonly called laboratory astrophysics. This work, carried out by a diverse community of laboratory astrophysicists, is increasingly important as astrophysics transitions into an era of precise measurement and high fidelity modeling.
[97]  [pdf] - 1092216
H2 Temperatures in the Crab Nebula
Comments: 16 pages, 5 figures
Submitted: 2011-12-07
We used K-band spectra to measure the H2 excitation temperatures in six molecular knots associated with the filaments in the Crab Nebula. The temperatures are quite high - in the range T ~ 2000-3000K, just below the H2 dissociation temperature. This is the temperature range over which the H2 1-0 S(1) line at 2.121\mum has its maximum emissivity per unit mass, so there may be many additional H2 cores with lower temperatures that are too faint to detect. We also measured the electron density in adjacent ionized gas, which on the assumption of gas pressure balance indicates densities in the molecular region n_mol ~ 20,000 H baryons cm-3, although this really is just a lower limit since the H2 gas may be confined by other means. The excited region may be just a thin skin on a much more extensive blob of molecular gas that does not have the correct temperature and density to be as easily detectable. At the opposite extreme, the observed knots could consist of a fine mist of molecular gas in which we are detecting essentially all of the H2. Future CO observations could distinguish between these two cases. The Crab filaments serve as the nearby laboratories for understanding the very much larger filamentary structures that have formed in the intracluster medium of cool-core galaxy clusters.
[98]  [pdf] - 1092012
Ion-by-ion Cooling efficiencies
Comments: Submitted to ApJS. Electronic data available at
Submitted: 2011-11-29
We present ion-by-ion cooling efficiencies for low-density gas. We use Cloudy (ver. 08.00) to estimate the cooling efficiencies for each ion of the first 30 elements (H-Zn) individually. We present results for gas temperatures between 1e4 and 1e8K, assuming low densities and optically thin conditions. When nonequilibrium ionization plays a significant role the ionization states deviate from those that obtain in collisional ionization equilibrium (CIE), and the local cooling efficiency at any given temperature depends on specific non-equilibrium ion fractions. The results presented here allow for an efficient estimate of the total cooling efficiency for any ionic composition. We also list the elemental cooling efficiencies assuming CIE conditions. These can be used to construct CIE cooling efficiencies for non-solar abundance ratios, or to estimate the cooling due to elements not explicitly included in any nonequilibrium computation. All the computational results are listed in convenient online tables.
[99]  [pdf] - 423198
How to Make an Atomic Blog in Your Own Kitchen. Summary of the Workshop: Uncertainties in Atomic Data and How They Propagate in Chemical Abundances
Comments: Summary contribution; 5 pages, no figures. The complete proceedings will be published on
Submitted: 2011-10-09
This workshop brought together scientists (including atomic physicists, theoretical astrophysicists and astronomers) concerned with the completeness and accuracy of atomic data for astrophysical applications. The topics covered in the workshop included the evaluation of uncertainties in atomic data, the propagation of such uncertainties in chemical abundances, and the feedback between observations and calculations. On a different level, we also discussed communication issues such as how to ensure that atomic data are correctly understood and used, and which forum is the best one for a fluid interaction between all communities involved in the production and use of atomic data. This paper reports on the discussions held during the workshop and introduces AstroAtom, a blog created as a platform for timely and open discussions on the needs and concerns over atomic data, and their effects on astronomical research. The complete proceedings will be published on
[100]  [pdf] - 1084115
Powering The Intra-cluster Filaments in Cool-Core Clusters of Galaxies
Comments: Atomic processes in plasmas - Proceedings of the 17th International Conference on Atomic processes in plasmas (2011) Edited by: KM Aggarwal and SFC Shearer (AIP)
Submitted: 2011-09-15
The first radio surveys of the sky discovered that some large clusters of galaxies contained powerful sources of synchrotron emission. Optical images showed that long linear filaments with bizarre emission-line spectra permeated the intra-cluster medium. Recent observations in the infrared and radio show that these filaments have very strong emission lines of molecular hydrogen and carbon monoxide. The mass of molecular material is quite large, the gas is quite warm, and the filaments have not formed stars despite their ~Gyr age. I will discuss the general astrophysical context of large clusters of galaxies and how large masses of molecular gas can be heated to produce what we observe. The unique properties of the filaments are a result of the unique environment. Magnetically confined molecular filaments are surrounded by the hot intra-cluster medium. Thermal particles with keV energies enter atomic and molecular regions and produce a shower of secondary nonthermal electrons. These secondaries collisionally heat, excite, dissociate, and ionize the cool gas. While ionization is dominated by these secondary particles, recombination is controlled by charge exchange, which produces the unusual optical emission line spectrum. I will describe some of the physical processes that are unique to this environment and outline some of the atomic physics issues.
[101]  [pdf] - 1084001
The curious conundrum regarding sulfur and oxygen abundances in planetary nebulae
Comments: 2 pages, 1 figure, proceedings of the IAU Symposium No. 283, Planetary Nebulae: An Eye to the Future
Submitted: 2011-09-12
We carefully consider numerous explanations for the sulfur abundance anomaly in planetary nebulae. No one rationale appears to be satisfactory, and we suggest that the ultimate explanation is likely to be a heretofore unidentified feature of the nebular gas which significantly impacts the sulfur ionization correction factor.
[102]  [pdf] - 1083385
Herschel observations of the Centaurus cluster - the dynamics of cold gas in a cool core
Comments: 17 pages, 9 figures, 5 tables, accepted for publication in MNRAS
Submitted: 2011-08-13
Brightest cluster galaxies (BCGs) in the cores of galaxy clusters have distinctly different properties from other low redshift massive ellipticals. The majority of the BCGs in cool-core clusters show signs of active star formation. We present observations of NGC 4696, the BCG of the Centaurus galaxy cluster, at far-infrared (FIR) wavelengths with the Herschel space telescope. Using the PACS spectrometer, we detect the two strongest coolants of the interstellar medium, CII at 157.74 micron and OI at 63.18 micron, and in addition NII at 121.90 micron. The CII emission is extended over a region of 7 kpc with a similar spatial morphology and kinematics to the optical H-alpha emission. This has the profound implication that the optical hydrogen recombination line, H-alpha, the optical forbidden lines, NII 6583 Angstrom, the soft X-ray filaments and the far-infrared CII line all have the same energy source. We also detect dust emission using the PACS and SPIRE photometers at all six wavebands. We perform a detailed spectral energy distribution fitting using a two-component modified black-body function and find a cold 19 K dust component with mass 1.6x10^6 solar mass and a warm 46 K dust component with mass 4.0x10^3 solar mass. The total FIR luminosity between 8 micron and 1000 micron is 7.5x10^8 solar luminosity, which using Kennicutt relation yields a low star formation rate of 0.13 solar mass per yr. This value is consistent with values derived from other tracers, such as ultraviolet emission. Combining the spectroscopic and photometric results together with optical H-alpha, we model emitting clouds consisting of photodissociation regions (PDRs) adjacent to ionized regions. We show that in addition to old and young stellar populations, there is another source of energy, such as cosmic rays, shocks or reconnection diffusion, required to excite the H-alpha and CII filaments.
[103]  [pdf] - 389465
Understanding strong molecular hydrogen emission in astronomical environments
Comments: Narrow-Line Seyfert 1 Galaxies and Their Place in the Universe (NLS1 Milan, Italy)
Submitted: 2011-07-22
Here I describe recent studies of objects with molecular hydrogen emission that is strong relative to other spectral lines. Large telescopes and fast spectrometers have made the 2 {\mu}m window accessible even for relatively faint objects. I summarize several environments where strong molecular hydrogen 2.121 micron emission is observed. The line is hard to excite due to its large excitation potential, and is most emissive in regions that have temperatures that are nearly high enough to dissociate molecular hydrogen. I outline several case studies. In the Helix planetary nebula strong emission is produced by rapidly flowing molecular gas that is exposed to an intense ionizing radiation field. This advective production of molecular hydrogen is a fundamentally non-equilibrium process. In the filaments surrounding brightest cluster galaxies in cool core clusters ionizing particles penetrate into magnetically confined molecular cores and excite the gas. Finally, I outline ongoing work on the Crab Nebula, where the first complete maps of molecular emission have only recently been completed. Both ionizing particles and high-energy photons may be important. Finally I speculate on the origin of the correlation between H2 / H I intensity ratios and other properties found in Active Galaxies. This is suggestive of a hardening of the radiation field along the Eigenvector 1 sequence. In all of this work I take the approach of understanding H2 emission along with emission from low and moderate ionization species, a necessary step if we are to really understand the context in which molecular hydrogen emission forms.
[104]  [pdf] - 1078102
A deep spectroscopic study of the filamentary nebulosity in NGC4696, the brightest cluster galaxy in the Centaurus cluster
Comments: 21 pages, 23 figures, accepted for publication in MNRAS
Submitted: 2011-07-20
We present results of deep integral field spectroscopy observations using high resolution optical (4150-7200 A) VIMOS VLT spectra, of NGC 4696, the dominant galaxy in the Centaurus cluster (Abell 3526). After the Virgo cluster, this is the second nearest (z=0.0104) example of a cool core cluster. NGC 4696 is surrounded by a vast, luminous H alpha emission line nebula (L = 2.2 \times 10^40 ergs per second). We explore the origin and excitation of the emission-line filaments and find their origin consistent with being drawn out, under rising radio bubbles, into the intracluster medium as in other similar systems. Contrary to previous observations we do not observe evidence for shock excitation of the outer filaments. Our optical spectra are consistent with the recent particle heating excitation mechanism of Ferland et al.
[105]  [pdf] - 1078063
Star formation in self-gravitating disks in active galactic nuclei. I. Metallicity gradients in broad line regions
Comments: 11 pages, 5 Figures, ApJ, Vol. 737, in press
Submitted: 2011-07-18
It has been suggested that the high metallicity generally observed in active galactic nuclei (AGNs) and quasars originates from ongoing star formation in the self-gravitating part of accretion disks around the supermassive black holes. We designate this region as the star forming (SF) disk, in which metals are produced from supernova explosions (SNexp) while at the same time inflows are driven by SNexp-excited turbulent viscosity to accrete onto the SMBHs. In this paper, an equation of metallicity governed by SNexp and radial advection is established to describe the metal distribution and evolution in the SF disk. We find that the metal abundance is enriched at different rates at different positions in the disk, and that a metallicity gradient is set up that evolves for steady-state AGNs. Metallicity as an integrated physical parameter can be used as a probe of the SF disk age during one episode of SMBH activity. In the SF disk, evaporation of molecular clouds heated by SNexp blast waves unavoidably forms hot gas. This heating is eventually balanced by the cooling of the hot gas, but we show that the hot gas will escape from the SF disk before being cooled, and diffuse into the BLRs forming with a typical rate of $\sim 1\sunmyr$. The diffusion of hot gas from a SF disk depends on ongoing star formation, leading to the metallicity gradients in BLR observed in AGNs. We discuss this and other observable consequences of this scenario.
[106]  [pdf] - 1051548
Structure and Feedback in 30 Doradus II. Structure and Chemical Abundances
Comments: 54 pages, 21 figures, Accepted for publication in the Astrophysical Journal
Submitted: 2011-01-19, last modified: 2011-07-06
We use our new optical-imaging and spectrophotometric survey of key diagnostic emission lines in 30 Doradus, together with CLOUDY photoionization models, to study the physical conditions and ionization mechanisms along over 4000 individual lines of sight at points spread across the face of the extended nebula, out to a projected radius 75 pc from R136 at the center of the ionizing cluster NGC 2070. We focus on the physical conditions, geometry and importance of radiation pressure on a point-by-point basis, with the aim of setting observational constraints on important feedback processes. We find that the dynamics and large scale structure of 30 Dor are set by a confined system of X-ray bubbles in rough pressure equilibrium with each other and with the confining molecular gas. Although the warm (10,000K) gas is photoionized by the massive young stars in NGC 2070, the radiation pressure does not currently play a major role in shaping the overall structure. The completeness of our survey also allows us to create a composite spectrum of 30 Doradus, simulating the observable spectrum of a spatially-unresolved, distant giant extragalactic H II region. We find that the highly simplified models used in the "strong line" abundance technique do in fact reproduce our observed lines strengths and deduced chemical abundances, in spite of the more than one order of magnitude range in the ionization parameter and density of the actual gas in 30 Dor.
[107]  [pdf] - 1077636
Critical Differences and Clues in Eta Car's 2009 Event
Comments: 45 pages, 9 figures, submitted to ApJ
Submitted: 2011-06-29
We monitored Eta Carinae with HST WFPC2 and Gemini GMOS throughout the 2009 spectroscopic event, which was expected to differ from its predecessor in 2003 (Davidson et al. 2005). Here we report major observed differences between events, and their implications. Some of these results were quite unexpected. (1) The UV brightness minimum was much deeper in 2009. This suggests that physical conditions in the early stages of an event depend on different parameters than the "normal" inter-event wind. Extra mass ejection from the primary star is one possible cause. (2) The expected He II 4687 brightness maximum was followed several weeks later by another. We explain why this fact, and the timing of the 4687 maxima, strongly support a "shock breakup" hypothesis for X-ray and 4687 behavior as proposed 5-10 years ago. (3) We observed a polar view of the star via light reflected by dust in the Homunculus nebula. Surprisingly, at that location the variations of emission-line brightness and Doppler velocities closely resembled a direct view of the star; which should not have been true for any phenomena related to the orbit. This result casts very serious doubt on all the proposed velocity interpretations that depend on the secondary star's orbital motion. (4) Latitude-dependent variations of H I, He I and Fe II features reveal aspects of wind behavior during the event. In addition, we discuss implications of the observations for several crucial unsolved problems.
[108]  [pdf] - 1053275
N II 5668-5712, a New Class of Spectral Features in Eta Carinae
Comments: 4 figures, 1 table
Submitted: 2011-04-10, last modified: 2011-05-20
We report on the N II 5668-5712 emission and absorption lines in the spectrum of Eta Carinae. Spectral lines of the stellar wind regions can be classified into four physically distinct categories: 1) low-excitation emission such as H I and Fe II, 2) higher excitation He I features, 3) the N II lines discussed in this paper, and 4) He II emission. These categories have different combinations of radial velocity behavior, excitation processes, and dependences on the secondary star. The N II lines are the only known features that originate in "normal" undisturbed zones of the primary wind but depend primarily on the location of the hot secondary star. N II probably excludes some proposed models, such as those where He I lines originate in the secondary star's wind or in an accretion disk.
[109]  [pdf] - 1076496
The energy source of the filaments around the giant galaxy NGC1275
Comments: 7 pages, 4 figures, MNRAS in press
Submitted: 2011-05-09
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.
[110]  [pdf] - 1053064
A Survey of Molecular Hydrogen in the Crab Nebula
Comments: Accepted for publication in the ApJS
Submitted: 2011-03-30
We have carried out a near-infrared, narrow-band imaging survey of the Crab Nebula, in the H2 2.12 micron and Br-gamma 2.17 micron lines, using the Spartan Infrared camera on the SOAR Telescope. Over a 2.8' x 5.1' area that encompasses about 2/3 of the full visible extent of the Crab, we detect 55 knots that emit strongly in the H2 line. We catalog the observed properties of these knots. We show that they are in or next to the filaments that are seen in optical-passband emission lines. Comparison to HST [S II] and [O III] images shows that the H2 knots are strongly associated with compact regions of low-ionization gas. We also find evidence of many additional, fainter H2 features, both discrete knots and long streamers following gas that emits strongly in [S II]. A pixel-by-pixel analysis shows that about 6 percent of the Crab's projected surface area has significant H2 emission that correlates with [S II] emission. We measured radial velocities of the [S II] lambda6716 emission lines from 47 of the cataloged knots and find that most are on the far (receding) side of the nebula. We also detect Br-gamma emission. It is right at the limit of our survey, and our Br-gamma filter cuts off part of the expected velocity range. But clearly the Br-gamma emission has a quite different morphology than the H2 knots, following the long linear filaments that are seen in H-alpha and in [O III] optical emission lines.
[111]  [pdf] - 1052726
Physical Conditions in Barnard's Loop, Components of the Orion-Eridanus Bubble, and Implications for the WIM Component of the ISM
Submitted: 2011-03-14
We have supplemented existing spectra of Barnard's Loop with high accuracy spectrophotometry of one new position. Cloudy photoionization models were calculated for a variety of ionization parameters and stellar temperatures and compared with the observations. After testing the procedure with recent observations of M43, we establish that Barnard's Loop is photoionized by four candidate ionizing stars, but agreement between the models and observations is only possible if Barnard's Loop is enhanced in heavy elements by about a factor of 1.4. Barnard's Loop is very similar in properties to the brightest components of the Orion-Eridanus Bubble and the Warm Ionized Medium (WIM). We are able to establish models that bound the range populated in low-ionization color-color diagrams (I([SII])/I(H{\alpha}) versus I([NII])/I(H{\alpha})) using only a limited range of ionization parameters and stellar temperatures. Previously established variations in the relative abundance of heavy elements render uncertain the most common method of determining electron temperatures for components of the Orion-Eridanus Bubble and the WIM based on only the I([NII])/I(H{\alpha}) ratio, although we confirm that the lowest surface brightness components of the WIM are on average of higher electron temperature. The electron temperatures for a few high surface brightness WIM components determined by direct methods are comparable to those of classical bright H II regions. In contrast, the low surface brightness HII regions studied by the Wisconsin H{\alpha} Mapper are of lower temperatures than the classical bright HII regions.
[112]  [pdf] - 1040831
Structure and Feedback in 30 Doradus I: Observations
Comments: 42 pages, 16 figures, Accepted for publication in ApJS
Submitted: 2010-09-24
We have completed a a new optical imaging and spectrophotometric survey of a 140 x 80 pc$^2$ region of 30 Doradus centered on R136, covering key optical diagnostic emission lines including \Ha, \Hb, \Hg, [O III] $\lambda\lambda$4363, 4959, 5007, [N II] $\lambda\lambda$6548, 6584, [S II] $\lambda\lambda$6717, 6731 [S III] $\lambda $6312 and in some locations [S III] $\lambda$9069. We present maps of fluxes and intensity ratios for these lines, and catalogs of isolated ionizing stars, elephant-trunk pillars, and edge-on ionization fronts. The final science-quality spectroscopic data products are available to the public. Our analysis of the new data finds that, while stellar winds and supernovae undoubtedly produce shocks and are responsible for shaping the nebula, there are no global spectral signatures to indicate that shocks are currently an important source of ionization. We conclude that the considerable region covered by our survey is well described by photoionization from the central cluster where the ionizing continuum is dominated by the most massive O stars. We show that if 30 Dor were viewed at a cosmological distance, its integrated light would be dominated by its extensive regions of lower surface-brightness rather than by the bright, eye-catching arcs.
[113]  [pdf] - 1034942
Detection of optical coronal emission from 10^6 K gas in the core of the Centaurus cluster
Comments: 11 pages, 13 figures, 2 tables, accepted for publication in MNRAS
Submitted: 2010-09-14
We report a detection (3.5x10^37 \pm 5.6x10^36 ergps) of the optical coronal emission line [Fe X]6374 and upper limits of four other coronal lines using high resolution VIMOS spectra centred on NGC 4696, the brightest cluster galaxy in the Centaurus cluster. Emission from these lines is indicative of gas at temperatures between 1 and 5 million K so traces the interstellar gas in NGC 4696. The rate of cooling derived from the upper limits is consistent with the cooling rate from X-ray observations (~10 solar masses per year) however we detect twice the luminosity expected for [Fe X]6374 emission, at 1 million K, our lowest temperature probe. We suggest this emission is due to the gas being heated rather than cooling out of the intracluster medium. We detect no coronal lines from [Ca XV], which are expected from the 5 million K gas seen near the centre in X-rays with Chandra. Calcium is however likely to be depleted from the gas phase onto dust grains in the central regions of NGC 4696.
[114]  [pdf] - 1034307
Spitzer reveals what's behind Orion's Bar
Comments: 60 pages, 16 figures, 10 tables. MNRAS accepted
Submitted: 2010-08-16
We present Spitzer Space Telescope observations of 11 regions SE of the Bright Bar in the Orion Nebula, along a radial from the exciting star theta1OriC, extending from 2.6 to 12.1'. Our Cycle 5 programme obtained deep spectra with matching IRS short-high (SH) and long-high (LH) aperture grid patterns. Most previous IR missions observed only the inner few arcmin. Orion is the benchmark for studies of the ISM particularly for elemental abundances. Spitzer observations provide a unique perspective on the Ne and S abundances by virtue of observing the dominant ionization states of Ne (Ne+, Ne++) and S (S++, S3+) in Orion and H II regions in general. The Ne/H abundance ratio is especially well determined, with a value of (1.01+/-0.08)E-4. We obtained corresponding new ground-based spectra at CTIO. These optical data are used to estimate the electron temperature, electron density, optical extinction, and the S+/S++ ratio at each of our Spitzer positions. That permits an adjustment for the total gas-phase S abundance because no S+ line is observed by Spitzer. The gas-phase S/H abundance ratio is (7.68+/-0.30)E-6. The Ne/S abundance ratio may be determined even when the weaker hydrogen line, H(7-6) here, is not measured. The mean value, adjusted for the optical S+/S++ ratio, is Ne/S = 13.0+/-0.6. We derive the electron density versus distance from theta1OriC for [S III] and [S II]. Both distributions are for the most part decreasing with increasing distance. A dramatic find is the presence of high-ionization Ne++ all the way to the outer optical boundary ~12' from theta1OriC. This IR result is robust, whereas the optical evidence from observations of high-ionization species (e.g. O++) at the outer optical boundary suffers uncertainty because of scattering of emission from the much brighter inner Huygens Region.
[115]  [pdf] - 1026425
A Sea Change in Eta Carinae
Comments: 11 pages, 3 figures, accepted for publication in ApJL
Submitted: 2010-04-20, last modified: 2010-06-01
Major stellar-wind emission features in the spectrum of Eta Car have recently decreased by factors of order 2 relative to the continuum. This is unprecedented in the modern observational record. The simplest, but unproven, explanation is a rapid decrease in the wind density.
[116]  [pdf] - 1026675
A Bright Molecular Core in a Crab Nebula Filament
Comments: Accepted for publication in ApJ Letters
Submitted: 2010-05-05, last modified: 2010-05-24
In a sub-arcsec near-infrared survey of the Crab Nebula using the new Spartan Infrared Camera, we have found several knots with high surface brightness in the H_2 2.12 micron line and a very large H_2 2.12 micron to Br-gamma ratio. The brightest of these knots has an intensity ratio I(H_2 2.12 micron)/I(Br-gamma) = 18+/-9, which we show sets a lower limit on the ratio of masses in the molecular and recombination (i.e. ionized) zones M_mol / M_rec >/- 0.9, and a total molecular mass within this single knot M_mol >/- 5E-5 M_sun. We argue that the knot discussed here probably is able to emit so strongly in the 2.12 micron line because its physical conditions are better tuned for such emission than is the case in other filaments. It is unclear whether this knot has an unusually large M_mol / M_rec ratio, or if many other Crab filaments also have similar amounts of molecular gas which is not emitting because the physical conditions are not so well tuned.
[117]  [pdf] - 1026739
Herschel images of NGC 6720: H2 formation on dust grains
Comments: Herschel first results letter; 7 pages, 6 figures
Submitted: 2010-05-10
Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots.
[118]  [pdf] - 1026720
Herschel photometry of brightest cluster galaxies in cooling flow clusters
Comments: Accepted for A&A Herschel Special Issue, 7 pages, 3 figures
Submitted: 2010-05-07, last modified: 2010-05-10
The dust destruction timescales in the cores of clusters of galaxies are relatively short given their high central gas densities. However, substantial mid-infrared and sub-mm emission has been detected in many brightest cluster galaxies. In this letter we present Herschel PACS and SPIRE photometry of the brightest cluster galaxy in three strong cooling flow clusters, A1068, A2597 and Zw3146. This photometry indicates that a substantial mass of cold dust is present (>3 x 10^7 Mo) at temperatures significantly lower (20-28K) than previously thought based on limited MIR and/or sub-mm results. The mass and temperature of the dust appear to match those of the cold gas traced by CO with a gas-to-dust ratio of 80-120.
[119]  [pdf] - 1026719
Herschel observations of FIR emission lines in brightest cluster galaxies
Comments: Accepted for the A&A Herschel Special Issue, 5 pages, 2 figures
Submitted: 2010-05-07, last modified: 2010-05-10
The question of how much gas cools in the cores of clusters of galaxies has been the focus of many, multiwavelength studies in the past 30 years. In this letter we present the first detections of the strongest atomic cooling lines, [C II], [O I] and [N I] in two strong cooling flow clusters, A1068 and A2597, using Herschel PACS. These spectra indicate that the substantial mass of cold molecular gas (>10^9 Mo) known to be present in these systems is being irradiated by intense UV radiation, most probably from young stars. The line widths of these FIR lines indicate that they share dynamics similar but not identical to other ionised and molecular gas traced by optical, near-infrared and CO lines. The relative brightness of the FIR lines compared to CO and FIR luminosity is consistent with other star-forming galaxies indicating that the properties of the molecular gas clouds in cluster cores and the stars they form are not unusual. These results provide additional evidence for a reservoir of cold gas that is fed by the cooling of gas in the cores of the most compact clusters and provide important diagnostics of the temperature and density of the dense clouds this gas resides in.
[120]  [pdf] - 1026261
A Bare Molecular Cloud at z~0.45
Comments: 15 pages, 4 figures, 4 tables, ApJ accepted
Submitted: 2010-04-12
Several neutral species (MgI, SiI, CaI, FeI) have been detected in a weak MgII absorption line system (W_r(2796)~0.15 Angstroms) at z~0.45 along the sightline toward HE0001-2340. These observations require extreme physical conditions, as noted in D'Odorico (2007). We place further constraints on the properties of this system by running a wide grid of photoionization models, determining that the absorbing cloud that produces the neutral absorption is extremely dense (~100-1000/cm^3), cold (<100 K), and has significant molecular content (~72-94%). Structures of this size and temperature have been detected in Milky Way CO surveys, and have been predicted in hydrodynamic simulations of turbulent gas. In order to explain the observed line profiles in all neutral and singly ionized chemical transitions, the lines must suffer from unresolved saturation and/or the absorber must partially cover the broad emission line region of the background quasar. In addition to this highly unusual cloud, three other ordinary weak MgII clouds (within densities of ~0.005/cm^3 and temperatures of ~10000K) lie within 500 km/s along the same sightline. We suggest that the "bare molecular cloud", which appears to reside outside of a galaxy disk, may have had in situ star formation and may evolve into an ordinary weak MgII absorbing cloud.
[121]  [pdf] - 902648
High-excitation emission lines near eta Carinae, and its likely companion star
Submitted: 2009-12-05
In order to study the distribution of gas and ionizing radiation around eta Car and their implications for its likely companion star, we have examined high-excitation emission lines of [Ne III], [Fe III], etc., in spectra obtained with the HST/STIS instrument during 1998-2004. Our principal results, some of them unexpected, are: (1) The high-excitation fluxes varied systematically and non-trivially throughout eta Car's 5.5-year spectroscopic cycle. Instead of rising to a plateau after the 1998 "event," they changed continuously with a maximum in mid-cycle. (2) At one significant location a brief, strong secondary maximum occurred just before the 2003.5 spectroscopic event. (3) These emission lines are strongly concentrated at the "Weigelt knots" several hundred AU northwest of the star. With less certainty, [Ne III] appears to be somewhat more concentrated than [Fe III]. (4) A faster, blue-shifted component of each feature appears concentrated near the star and elongated perpendicular to the system's bipolar axis. This structure may be related to the equatorial outflow and/or to dense material known to exist along our line of sight to the star. (5) Using the photoionization program Cloudy, we estimated the range of parameters for the hot secondary star that would give satisfactory high-excitation line ratios in the ejecta. Teff = 39000 K and L = 400000 Lsun, for example, would be satisfactory. The allowed region in parameter space is wider (and mostly less luminous) than some previous authors suggested.
[122]  [pdf] - 18559
Properties of warm absorbers in active galaxies: a systematic stability curve analysis
Comments: The paper contains 10 figures and 2 tables. Published in MNRAS. Typos corrected, references elaborated and journal reference added. The definitive version is available at
Submitted: 2008-11-14, last modified: 2009-11-25
Signatures of warm absorbers are seen in soft X-ray spectra of about half of all Seyfert1 galaxies observed and in some quasars and blazars. We use the thermal equilibrium curve to study the influence of the shape of the ionizing continuum, density and the chemical composition of the absorbing gas on the existence and nature of the warm absorbers. We describe circumstances in which a stable warm absorber can exist as a multiphase medium or one with continuous variation in pressure. In particular we find the following results: i) the warm absorber exists only if the spectral index of the X-ray power-law ionizing continuum $\alpha > 0.2$ and has a multiphase nature if $\alpha \sim 0.8$, which interestingly is the spectral index for most of the observed Seyfert 1 galaxies; ii) thermal and ionization states of highly dense warm absorbers are sensitive to their density if the ionizing continuum is sufficiently soft, i.e. dominated by the ultraviolet iii) absorbing gas with super-Solar metallicity is more likely to have a multiphase nature; iv) the nature of the warm absorber is significantly influenced by the absence of iron and associated elements which are produced in the later stages of star formation history in supernovae of type Ia.
[123]  [pdf] - 902327
Implications of infalling Fe II - emitting clouds in active galactic nuclei: anisotropic properties
Comments: 5 figures, Accepted for publication in ApJ (Letters)
Submitted: 2009-11-06, last modified: 2009-11-19
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.
[124]  [pdf] - 1932174
Laboratory Studies for Planetary Sciences. A Planetary Decadal Survey White Paper Prepared by the American Astronomical Society (AAS) Working Group on Laboratory Astrophysics (WGLA)
Submitted: 2009-10-02
The WGLA of the AAS ( promotes collaboration and exchange of knowledge between astronomy and planetary sciences and the laboratory sciences (physics, chemistry, and biology). Laboratory data needs of ongoing and next generation planetary science missions are carefully evaluated and recommended in this white paper submitted by the WGLA to Planetary Decadal Survey.
[125]  [pdf] - 1003075
Electron impact excitation of OII fine-structure levels
Comments: 2009MNRAS.397..903K
Submitted: 2009-07-23
Effective collision strengths for forbidden transitions among the 5 energetically lowest finestructure levels of O II are calculated in the Breit-Pauli approximation using the R-matrix method. Results are presented for the electron temperature range 100 to 100 000 K. The accuracy of the calculations is evaluated via the use of different types of radial orbital sets and a different configuration expansion basis for the target wavefunctions. A detailed assessment of previous available data is given, and erroneous results are highlighted. Our results reconfirm the validity of the original Seaton and Osterbrock scaling for the optical O II ratio, a matter of some recent controversy. Finally we present plasma diagnostic diagrams using the best collision strengths and transition probabilities.
[126]  [pdf] - 25199
Rotationally Warm Molecular Hydrogen in the Orion Bar
Comments: Accepted for publication in ApJ (34 pages, including 16 figures)
Submitted: 2009-06-12
The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients.
[127]  [pdf] - 1932072
Roles and Needs of Laboratory Astrophysics in NASA's Space and Earth Science Mission
Comments: White Paper submitted by the AAS Working Group on Laboratory Astrophysics (WGLA) to the National Academy of Science, Space Sciences Board, Ad Hoc Committee studying "The Role and Scope of Mission-Enabling Activities in NASA's Space and Earth Science Missions"
Submitted: 2009-03-26
Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The mission enabling impact of laboratory astrophysics ranges from the scientific conception stage for airborne and space-based observatories, all the way through to the scientific return of these missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA. These efforts are necessary for the success of astronomical research being funded by NASA. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics (WGLA), have prepared this White Paper on the laboratory astrophysics infrastructure needed to maximize the scientific return from NASA's space and Earth sciences program.
[128]  [pdf] - 554081
Laboratory Astrophysics and the State of Astronomy and Astrophysics
Comments: Position paper submitted by the AAS Working Group on Laboratory Astrophysics (WGLA) to the State of the Profession (Facilities, Funding and Programs Study Group) of the Astronomy and Astrophysics Decadal Survey (Astro2010)
Submitted: 2009-03-13
Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The impact of laboratory astrophysics ranges from the scientific conception stage for ground-based, airborne, and space-based observatories, all the way through to the scientific return of these projects and missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA, NSF, and DOE. These efforts are necessary for the success of astronomical research being funded by the agencies. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics, have prepared this State of the Profession Position Paper on the laboratory astrophysics infrastructure needed to ensure the advancement of astronomy and astrophysics in the next decade.
[129]  [pdf] - 21919
Magnetic fields, stellar feedback, and the geometry of H II Regions
Comments: 10 pages, to appear in the proceedings of IAU Symposium 259: "Cosmic Magnetic Fields:from Planets, to Stars and Galaxies," edited by K.G. Strassmeier, A.G. Kosovichev, J. Beckmann
Submitted: 2009-03-01
Magnetic pressure has long been known to dominate over gas pressure in atomic and molecular regions of the interstellar medium. Here I review several recent observational studies of the relationships between the H^+, H^0 and H_2 regions in M42 (the Orion complex) and M17. A simple picture results. When stars form they push back surrounding material, mainly through the outward momentum of starlight acting on grains, and field lines are dragged with the gas due to flux freezing. The magnetic field is compressed and the magnetic pressure increases until it is able to resist further expansion and the system comes into approximate magnetostatic equilibrium. Magnetic field lines can be preferentially aligned perpendicular to the long axis of quiescent cloud before stars form. After star formation and pushback occurs ionized gas will be constrained to flow along field lines and escape from the system along directions perpendicular to the long axis. The magnetic field may play other roles in the physics of the H II region and associated PDR. Cosmic rays may be enhanced along with the field and provide additional heating of atomic and molecular material. Wave motions may be associated with the field and contribute a component of turbulence to observed line profiles.
[130]  [pdf] - 21886
New Discoveries in Planetary Systems and Star Formation through Advances in Laboratory Astrophysics
Comments: White paper submitted by the AAS Working Group on Laboratory Astrophysics (WGLA) to the PSF SFP of the Astronomy and Astrophysics Decadal Survey (Astro2010)
Submitted: 2009-02-27
As the panel on Planetary Systems and Star Formation (PSF) is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of solar system bodies (other than the Sun) and extrasolar planets, debris disks, exobiology, the formation of individual stars, protostellar and protoplanetary disks, molecular clouds and the cold ISM, dust, and astrochemistry. Central to the progress in these areas are the corresponding advances in laboratory astro- physics which are required for fully realizing the PSF scientific opportunities in the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics and chemistry which produce the observed spectra and describe the astrophysical processes. We discuss four areas of laboratory astrophysics relevant to the PSF panel: atomic, molecular, solid matter, and plasma physics. Section 2 describes some of the new opportunities and compelling themes which will be enabled by advances in laboratory astrophysics. Section 3 provides the scientific context for these opportunities. Section 4 discusses some experimental and theoretical advances in laboratory astrophysics required to realize the PSF scientific opportunities of the next decade. As requested in the Call for White Papers, we present in Section 5 four central questions and one area with unusual discovery potential. We give a short postlude in Section 6.
[131]  [pdf] - 554074
New Discoveries in Cosmology and Fundamental Physics through Advances in Laboratory Astrophysics
Comments: White paper submitted by the AAS Working Group on Laboratory Astrophysics (WGLA) to the CFP SFP of the Astronomy and Astrophysics Decadal Survey (Astro2010)
Submitted: 2009-02-26
As the Cosmology and Fundamental Physics (CFP) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the early universe, the microwave background, the reionization and galaxy formation up to virialization of protogalaxies, large scale structure, the intergalactic medium, the determination of cosmological parameters, dark matter, dark energy, tests of gravity, astronomically determined physical constants, and high energy physics using astronomical messengers. Central to the progress in these areas are the corresponding advances in laboratory astrophysics which are required for fully realizing the CFP scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics which produce the observed astrophysical processes. The 5 areas of laboratory astrophysics which we have identified as relevant to the CFP panel are atomic, molecular, plasma, nuclear, and particle physics. Here, Section 2 describes some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the CFP scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.
[132]  [pdf] - 554076
New Discoveries in Stars and Stellar Evolution through Advances in Laboratory Astrophysics
Comments: White paper submitted by the AAS Working Group on Laboratory Astrophysics (WGLA) to the SSE SFP of the Astronomy and Astrophysics Decadal Survey (Astro2010)
Submitted: 2009-02-26
As the Stars and Stellar Evolution (SSE) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the Sun as a star, stellar astrophysics, the structure and evolution of single and multiple stars, compact objects, SNe, gamma-ray bursts, solar neutrinos, and extreme physics on stellar scales. Central to the progress in these areas are the corresponding advances in laboratory astrophysics, required to fully realize the SSE scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produces the observed astrophysical processes. The 6 areas of laboratory astrophysics, which we have identified as relevant to the CFP panel, are atomic, molecular, solid matter, plasma, nuclear physics, and particle physics. In this white paper, we describe in Section 2 the scientific context and some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we discuss some of the experimental and theoretical advances in laboratory astrophysics required to realize the SSE scientific opportunities of the next decade. As requested in the Call for White Papers, Section 4 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 5.
[133]  [pdf] - 21859
New Discoveries in the Galactic Neighborhood through Advances in Laboratory Astrophysics
Comments: White paper submitted by the AAS Working Group on Laboratory Astrophysics (WGLA) to the GAN SFP of the Astronomy and Astrophysics Decadal Survey (Astro2010)
Submitted: 2009-02-26
As the Galactic Neighborhood (GAN) panel is fully aware, the next decade will see major advances in our understanding of this area of research. To quote from their charge, these advances will occur in studies of the galactic neighborhood, including the structure and properties of the Milky Way and nearby galaxies, and their stellar populations and evolution, as well as interstellar media and star clusters. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GAN scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics and chemistry that produces the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the GAN panel are atomic, molecular, solid matter, plasma, and nuclear physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GAN scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.
[134]  [pdf] - 554075
New Discoveries in Galaxies across Cosmic Time through Advances in Laboratory Astrophysics
Comments: White paper submitted by the AAS Working Group on Laboratory Astrophysics (WGLA) to the GCT SFP of the Astronomy and Astrophysics Decadal Survey (Astro2010)
Submitted: 2009-02-26
As the Galaxies across Cosmic Time (GCT) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the formation, evolution, and global properties of galaxies and galaxy clusters, as well as active galactic nuclei and QSOs, mergers, star formation rate, gas accretion, and supermassive black holes. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GCT scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produce the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the CFP panel are atomic, molecular, solid matter, plasma, nuclear, and particle physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GCT scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.
[135]  [pdf] - 18705
The location and kinematics of the coronal-line emitting regions in AGN
Comments: 6 Pages, 5 Figures, Accepted for publication in MNRAS Letters
Submitted: 2008-11-19
We use the photoionisation code Cloudy to determine both the location and the kinematics of the optical forbidden, high ionisation line (hereafter, FHIL) emitting gas in the narrow line Seyfert 1 galaxy Ark 564. The results of our models are compared with the observed properties of these emission lines to produce a physical model that is used to explain both the kinematics and the source of this gas. The main features of this model are that the FHIL emitting gas is launched from the putative dusty torus and is quickly accelerated to its terminal velocity of a few hundred km/s. Iron-carrying grains are destroyed during this initial acceleration. This velocity is maintained by a balance between radiative forces and gravity in this super-Eddington source. Eventually the outflow is slowed at large radii by the gravitational forces of and interactions with the host galaxy. In this model, FHIL emission traces the transition between the AGN and bulge zones of influence.
[136]  [pdf] - 315293
Orion's Bar: Physical Conditions across the Definitive H+ / H0 / H2 Interface
Comments: 32 pages, including 10 figures. ApJ (in press). Revised to correct 2 spelling errors
Submitted: 2008-11-07, last modified: 2008-11-14
Previous work has shown the Orion Bar to be an interface between ionized and molecular gas, viewed roughly edge on, which is excited by the light from the Trapezium cluster. Much of the emission from any star-forming region will originate from such interfaces, so the Bar serves as a foundation test of any emission model. Here we combine X-ray, optical, IR and radio data sets to derive emission spectra along the transition from H+ to H0 to H2 regions. We then reproduce the spectra of these layers with a simulation that simultaneously accounts for the detailed microphysics of the gas, the grains, and molecules, especially H2 and CO. The magnetic field, observed to be the dominant pressure in another region of the Orion Nebula, is treated as a free parameter, along with the density of cosmic rays. Our model successfully accounts for the optical, IR and radio observations across the Bar by including a significant magnetic pressure and also heating by an excess density of cosmic rays, which we suggest is due to cosmic rays being trapped in the compressed magnetic field. In the Orion Bar, as we had previously found in M17, momentum carried by radiation and winds from the newly formed stars pushes back and compresses the surrounding gas. There is a rough balance between outward momentum in starlight and the total pressure in atomic and molecular gas surrounding the H+ region. If the gas starts out with a weak magnetic field, the starlight from a newly formed cluster will push back the gas and compress the gas, magnetic field, and cosmic rays until magnetic pressure becomes an important factor.
[137]  [pdf] - 6799
Dielectronic recombination and stability of warm gas in AGN
Comments: 5 pages, 3 figures, Accepted for publication in MNRAS Letters. The definitive version is available at
Submitted: 2007-11-07, last modified: 2008-11-13
Understanding the thermal equilibrium (stability) curve may offer insights into the nature of the warm absorbers often found in active galactic nuclei. Its shape is determined by factors like the spectrum of the ionizing continuum and the chemical composition of the gas. We find that the stability curves obtained under the same set of the above mentioned physical factors, but using recently derived dielectronic recombination rates, give significantly different results, especially in the regions corresponding to warm absorbers, leading to different physical predictions. Using the current rates we find a larger probability of having thermally stable warm absorber at $10^5 \kel$ than previous predictions and also a greater possibility for its multiphase nature. the results obtained with the current dielectronic recombination rate coefficients are more reliable because the warm absorber models along the stability curve have computed coefficient values, whereas previous calculations relied on guessed averages for the same due to lack of available data.
[138]  [pdf] - 18486
Magnetic fields and the location of the PDR
Submitted: 2008-11-12
I review recent studies of the emission-line regions in Orion and M17. Both have similar geometries, a bubble of hot shocked gas surrounding the central star cluster, with H^+, H^0, and H_2 regions, often referred to as H II regions, PDRs, and molecular clouds, forming successive shells on the surface of a molecular cloud. The magnetic fields in the H^0 regions have been measured with 21 cm Zeeman polarization and are found to be 1 -- 2 dex stronger than the field in the diffuse ISM. The regions appear to be in rough hydrostatic equilibrium. The H^+ region is pushed away from the star cluster by starlight radiation pressure. Since most starlight is in ionizing radiation, most of its outward push will act on the H^+ region and then on to the H^0 region. The magnetic pressure in the H^0 region balances the momentum in starlight and together they set the location of the H^0 region. The picture is that, when the star cluster formed, it created a bubble of ionized gas which expanded and compressing surrounding H^0 and H_2 regions. The magnetic field was amplified until its pressure was able to support the momentum in starlight. This offers a great simplification in understanding the underlying physics that establishes parameters for PDR models.
[139]  [pdf] - 18317
Uncertainties in Theoretical HeI Emissivities: HII Regions, Primordial Abundance, and Cosmological Recombination
Comments: 5 pages, 3 figures, accepted to MNRAS Letters
Submitted: 2008-11-07
A number of recent works in astronomy and cosmology have relied upon theoretical He I emissivities, but we know of no effort to quantify the uncertainties in the atomic data. We analyze and assign uncertainties to all relevant atomic data, perform Monte Carlo analyses, and report standard deviations in the line emissivities. We consider two sets of errors, which we call "optimistic" and "pessimistic." We also consider three different conditions, corresponding to prototypical Galactic and extragalactic H II regions and the epoch of cosmological recombination. In the extragalactic H II case, the errors we obtain are comparable to or larger than the errors in some recent $Y_p$ calculations, including those derived from CMB observations. We demonstrate a systematic effect on primordial abundance calculations; this effect cannot be reduced by observing a large number of objects. In the cosmological recombination case, the errors are comparable to many of the effects considered in recent calculations.
[140]  [pdf] - 315261
Fluorescent excitation of Balmer lines in gaseous nebulae: case D
Comments: 49 pages, 14 figures. AAS Latex. Accepted for publication in the Astrophysical Journal
Submitted: 2008-10-30
(abridged) Non-ionizing stellar continua are a source of photons for continuum pumping in the hydrogen Lyman transitions. In the environments where these transitions are optically thick, deexcitation occurs through higher series lines, so that the flux in these lines has a fluorescent contribution in addition to recombination; in particular, Balmer emissivities are systematically enhanced above case B. The effectiveness of such mechanism in HII regions and the adequacy of photoionization models as a tool to study it are the two main focuses of this work. We find that photoionization models of H II regions illuminated by low-resolution population synthesis models significantly overpredict the fluorescent contribution to the Balmer lines. Conversely, photoionization models in which the non-ionizing part of the continuum is omitted or is not transferred underpredict the fluorescent contribution to the Balmer lines, producing a bias of similar amplitude in the opposite direction. In this paper, we carry out realistic estimations of the fluorescent Balmer intensity and discuss the variations to be expected as the simulated observational setup and the stellar population's parameters are varied. In all the cases explored, we find that fluorescent excitation provides a significant contribution. We also show that differential fluorescent enhancement may produce line-of-sight differences in the Balmer decrement, mimicking interstellar extinction. Fluorescent excitation emerges from our study as a small but important mechanism for the enhancement of Balmer lines, which should be taken into account in the abundance analysis of photoionized regions, particularly in the case of high-precision applications such as the determination of primordial helium.
[141]  [pdf] - 17992
Collisional heating as the origin of filament emission in galaxy clusters
Comments: 28 pages, 21 figures, accepted to MNRAS
Submitted: 2008-10-29
[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.
[142]  [pdf] - 17785
The three dimensional dynamic structure of the inner Orion Nebula
Comments: AJ in press, 18 pages, 7 figures (2 in full color)
Submitted: 2008-10-24
The three dimensional structure of the brightest part of the Orion Nebula is assessed in the light of published and new data. We find that the widely accepted model of a concave blister of ionized material needs to be altered in the southwest direction from the Trapezium, where we find that the Orion-S feature is a separate cloud of very optically thick molecules within the body of ionized gas, which is probably the location of the multiple embedded sources that produce the outflows that define the Orion-S star formation region. Evidence for this cloud comes from the presence of H2CO lines in absorption in the radio continuum and discrepancies in the extinction derived from radio-optical and optical only emission. We present an equilibrium Cloudy model of the Orion-S cloud, which successfully reproduces many observed properties of this feature. We also report the discovery of an open-sided shell of [O III] surrounding the Trapezium stars, revealed through emission line ratio images and the onset of radiation shadows beyond some proplyds. We show that the observed properties of the shell are consistent with it being a stationary structure, produced by shock interactions between the ambient nebular gas and the high-velocity wind from theta^1 Ori C. We examine the implications of the recently published evidence for a large blueshifted velocity of theta^1 Ori C with respect to the Orion Molecular Cloud, which could mean that this star has only recently begun to photoionize the Orion Nebula. We show that current observations of the Nebula do not rule out such a possibility, so long as the ionization front has propagated into a pre-existing low-density region. In addition, a young age for the Nebula would help explain the presence of nearby proplyds with a short mass-loss timescale to photoablation.
[143]  [pdf] - 15492
Sensitivity of PDR Calculations to Microphysical Details
Comments: 36 pages, 17 figures, accepted for publication in ApJ
Submitted: 2008-08-19
Our understanding of physical processes in Photodissociation regions or Photon Dominated Regions (PDRs) largely depends on the ability of spectral synthesis codes to reproduce the observed infrared emission-line spectrum. In this paper, we explore the sensitivity of a single PDR model to microphysical details. Our calculations use the Cloudy spectral synthesis code, recently modified to include a wealth of PDR physical processes. We show how the chemical/thermal structure of a PDR, along with the calculated spectrum, changes when the treatment of physical processes such as grain physics and atomic/molecular rates are varied. We find a significant variation in the intensities of PDR emission lines, depending on different treatments of the grain physics. We also show how different combinations of the cosmic-ray ionization rate, inclusion of grain-atom/ion charge transfer, and the grain size distribution can lead to very similar results for the chemical structure. Additionally, our results show the utility of Cloudy for the spectral modeling of molecular environments.
[144]  [pdf] - 12526
Rotational quenching rate coefficients for H_2 in collisions with H_2 from 2 to 10,000 K
Submitted: 2008-05-12
Rate coefficients for rotational transitions in H_2 induced by H_2 impact are presented. Extensive quantum mechanical coupled-channel calculations based on a recently published (H_2)_2 potential energy surface were performed. The potential energy surface used here is presumed to be more reliable than surfaces used in previous work. Rotational transition cross sections with initial levels J <= 8 were computed for collision energies ranging between 0.0001 and 2.5 eV, and the corresponding rate coefficients were calculated for the temperature range 2 < T <10,000 K. In general, agreement with earlier calculations, which were limited to 100-6000 K, is good though discrepancies are found at the lowest and highest temperatures. Low-density-limit cooling functions due to para- and ortho-H_2 collisions are obtained from the collisional rate coefficients. Implications of the new results for non-thermal H_2 rotational distributions in molecular regions are also investigated.
[145]  [pdf] - 11954
Conference Summary: The Cosmic Agitator - Magnetic Fields in the Galaxy
Comments: 9 pages, 2008 Mar 26-29 Conference Summary
Submitted: 2008-04-21
We present a summary of the conference "The Cosmic Agitator: Magnetic Fields in the Galaxy" held in Lexington KY in 2008 Mar 26-29. The presentation draws primarily from material in the slides prepared for the Conference Summary by one of us (Carl Heiles). Interested readers may navigate to the conference web site given in the paper to view the posted presentations in detail.
[146]  [pdf] - 10198
The origin of molecular hydrogen emission in cooling-flow filaments
Comments: 5 pages, 4 figures, accepted to MNRAS Letters
Submitted: 2008-02-18
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.
[147]  [pdf] - 7446
Merged ionization/dissociation fronts in planetary nebulae
Comments: 4 pages, accepted by ApJL, scheduled December 20 issue
Submitted: 2007-11-27
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.
[148]  [pdf] - 6656
On the enhanced cosmic-ray ionization rate in the diffuse cloud towards Zeta Persei
Comments: 25 pages, 4 figures, accepted for publication in ApJ
Submitted: 2007-11-02
The spatial distribution of the cosmic-ray flux is important in understanding the Interstellar Medium (ISM) of the Galaxy. This distribution can be analyzed by studying different molecular species along different sight lines whose abundances are sensitive to the cosmic-ray ionization rate. Recently several groups have reported an enhanced cosmic-ray ionization rate in diffuse clouds compared to the standard value, zeta(average)=2.5e-17 s^-1, measured toward dense molecular clouds. In an earlier work we reported an enhancement in the cosmic ray rate of 20 towards HD185418. McCall et al. have reported enhancements of 48 towards zeta Persei based on the observed abundance of H$_{3}^+ while Le Petit et al. found a cosmic ray enhancement of ~10 to be consistent with their models for this same sight line. Here we revisit zeta Persei and perform a detailed calculation using a self-consistent treatment of the hydrogen chemistry, grain physics, energy and ionization balance, and excitation physics. We show that the value of zeta deduced from the H$_{3}^+$ column density in the diffuse region of the sightline depends strongly on the properties of the grains because they remove free electrons and change the hydrogen chemistry. The observations are largely consistent with a cosmic ray enhancement of 40, with several diagnostics indicating higher values. This underscores the importance of a full treatment of grain physics in studies of interstellar chemistry.
[149]  [pdf] - 89485
Discovery of atomic and molecular mid infra-red emission lines in off-nuclear regions of NGC 1275 and NGC 4696 with the Spitzer Space Telescope
Comments: 15 pages. 13 figures. 12 tables. Accepted for publication in MNRAS. Moderate revision following referee's comments
Submitted: 2007-02-16, last modified: 2007-09-12
We present Spitzer high-resolution spectra of off-nuclear regions in the central cluster galaxies NGC 1275 and NGC 4696 in the Perseus and Centaurus clusters, respectively. Both objects are surrounded by extensive optical emission-line filamentary nebulae, bright outer parts of which are the targets of our observations. The 10-37um spectra show strong pure rotational lines from molecular hydrogen revealing a molecular component to the filaments which has an excitation temperature of ~300-400K. The flux in the 0-0S(1) molecular hydrogen line correlates well with the strength of the optical lines, having about 3 per cent of the Halpha+[NII] emission. The 11.3um PAH feature is seen in some spectra. Emission is also seen from both low and high ionization fine structure lines. Molecular hydrogen cooler than ~400K dominates the mass of the outer filaments; the nebulae are predominantly molecular.
[150]  [pdf] - 1944041
Revisiting He-like X-ray Emission Line Plasma Diagnostics
Comments: 28 pages, 7 figures, accepted to ApJ
Submitted: 2007-04-19
A complete model of helium-like line and continuum emission has been incorporated into the plasma simulation code Cloudy. All elements between He and Zn are treated, any number of levels can be considered, and radiative and collisional processes are included. This includes photoionization from all levels, line transfer including continuum pumping and destruction by background opacities, scattering, and collisional processes. The model is calculated self-consistently along with the ionization and thermal structure of the surrounding nebula. The result is a complete line and continuum spectrum of the plasma. Here we focus on the ions of the He I sequence and reconsider the standard helium-like X-ray diagnostics. We first consider semi-analytical predictions and compare these with previous work in the low-density, optically-thin limit. We then perform numerical calculations of helium-like X-ray emission (such as is observed in some regions of Seyferts) and predict line ratios as a function of ionizing flux, hydrogen density, and column density. In particular, we demonstrate that, in photoionized plasmas, the $R$-ratio, a density indicator in a collisional plasma, depends on the ionization fraction and is strongly affected by optical depth for large column densities. We also introduce the notion that the $R$-ratio is a measure of the incident continuum at UV wavelengths. The $G$-ratio, which is temperature-sensitive in a collisional plasma, is also discussed, and shown to be strongly affected by continuum pumping and optical depth as well. These distinguish a photoionized plasma from the more commonly studied collisional case.
[151]  [pdf] - 316832
Optically Thick [O I] and [C II] Emission Towards NGC 6334 A
Comments: 27 pages, 8 figures, accepted for publication in the Astrophysical Journal
Submitted: 2007-02-09
This work focuses on [O I] and [C II] emission towards NGC 6334 A, an embedded H+ region/PDR only observable at infrared or longer wavelengths. A geometry where nearly all the emission escapes out the side of the cloud facing the stars, such as Orion, is not applicable to this region. Instead, we find the geometry to be one where the H+ region and associated PDR is embedded in the molecular cloud. Constant density PDR calculations are presented which predict line intensities as a function of AV (or N(H)), hydrogen density (nH), and incident UV radiation field (G0). We find that a single component model with AV ~650 mag, nH = 5x10^5 cm-3, and G0 = 7x10^4 reproduces the observed [O I] and [C II] intensities, and that the low [O I] 63 to 146 micron ratio is due to line optical depth effects in the [O I] lines, produced by a large column density of atomic/molecular gas. We find that the effects of a density-law would increase our derived AV, while the effects of an asymmetric geometry would decrease AV, with the two effects largely canceling. We conclude that optically selected H+ regions adjacent to PDRs, such as Orion, likely have a different viewing angle or geometry than similar regions detected through IR observations. Overall, the theoretical calculations presented in this work have utility for any PDR embedded in a molecular cloud.
[152]  [pdf] - 338144
A PDR-Code Comparison Study
Comments: 28 pages, 16 figures, see all benchmark data under:
Submitted: 2007-02-08
We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.
[153]  [pdf] - 88765
Determination of the Physical Conditions of the Knots in the Helix Nebula from Optical and Infrared Observations
Comments: Astronomical Journal, in press. Some figures are shown at reduced resolution. A full resolution version is available at
Submitted: 2007-01-22
[Abridged] We use new HST and archived images to clarify the nature of the knots in the Helix Nebula. We employ published far infrared spectrophotometry and existing 2.12 micron images to establish that the population distribution of the lowest ro-vibrational states of H2 is close to the distribution of a gas in LTE at 988 +- 119 K. We derive a total flux from the nebula in H2 lines and compare this with the power available from the central star for producing this radiation. We establish that neither soft X-rays nor FUV radiation has enough energy to power the H2 radiation, only the stellar EUV radiation shortward of 912 Angstrom does. Advection of material from the cold regions of the knots produces an extensive zone where both atomic and molecular hydrogen are found, allowing the H2 to directly be heated by Lyman continuum radiation, thus providing a mechanism that can explain the excitation temperature and surface brightness of the cusps and tails. New images of the knot 378-801 reveal that the 2.12 micron cusp and tail lie immediately inside the ionized atomic gas zone. This firmly establishes that the "tail" structure is an ionization bounded radiation shadow behind the optically thick core of the knot. A unique new image in the HeII 4686 Angstrom line fails to show any emission from knots that might have been found in the He++ core of the nebula. We also re-examined high signal-to-noise ratio ground-based telescope images of this same inner region and found no evidence of structures that could be related to knots.
[154]  [pdf] - 88632
Quasar Metal Abundances and Host Galaxy Evolution
Comments: 10 pages, 5 figures in 6 files, invited review for the conference "The Central Engines of AGN," Xi'an, China, October 2006. ASP Conference Series, in press
Submitted: 2007-01-17
Quasars signal a unique phase of galaxy evolution -- when massive spheroids are rapidly being assembled, forming stars and growing their central super-massive black holes. Measurements of the metal abundances around quasars provide unique information about these complex evolutionary processes. Here we provide a brief review of the current status and implications of quasar abundance research.
[155]  [pdf] - 86965
He~I Emission in the Orion Nebula and Implications for Primordial Helium Abundance
Comments: accepted for publication in ApJ
Submitted: 2006-11-17
We apply a recently developed theoretical model of helium emission to observations of both the Orion Nebula and a sample of extragalactic H II regions. In the Orion analysis, we eliminate some weak and blended lines and compare theory and observation for our reduced line list. With our best theoretical model we find an average difference between theoretical and observed intensities $<I_{predicted}/I_{observed}-1> = 6.5%$. We argue that both the red and blue ends of the spectrum may have been inadequately corrected for reddening. For the 22 highest quality lines, with $3499 {\AA}\le\lambda\le6678 {\AA}$, our best model predicts observations to an average of 3.8%. We also perform an analysis of the reported observational errors and conclude they have been underestimated. In the extragalactic analysis, we demonstrate the likelihood of a large systematic error in the reported data and discuss possible causes. This systematic error is at least as large as the errors associated with nearly all attempts to calculate the primordial helium abundance from such observations. Our Orion analysis suggests that the problem does not lie in the theoretical models. We demonstrate a correlation between equivalent width and apparent helium abundance of lines from extragalactic sources that is most likely due to underlying stellar absorption. Finally, we present fits to collisionless case-B He I emissivities as well as the relative contributions due to collisional excitations out of the metastable $2s {}^{3}S$ term.
[156]  [pdf] - 85954
The Structure of the Homunculus. II. Modeling the physical conditions in Eta Car's molecular shell
Comments: 18 pages, 5 figures. Accepted by ApJ
Submitted: 2006-10-17
We present models that reproduce the observed double-shell structure of the Homunculus Nebula around eta Carinae, including the stratification of infrared H2 and [FeII] emission seen in data obtained with the Phoenix spectrograph on Gemini South, as well as the corresponding stratified grain temperature seen in thermal-infrared data. Tuning the model to match the observed shell thickness allows us to determine the threshold density which permits survival of H2. An average hydrogen density of n_H=(0.5-1)x10^7 cm-3 in the outer zone is required to allow H2 to exist at all latitudes in the nebula, and for Fe+ to recombine. This gives independent confirmation of the very large mass of the Homunculus, indicating a total of roughly 15--35 Msun (although we note reasons why the lower end of this range is favored). At the interface between the atomic and molecular zones, we predict a sharp drop in the dust temperature, in agreement with the bimodal dust color temperatures observed in the two zones. In the outer molecular shell, the dust temperature drops to nearly the blackbody temperature, and becomes independent of grain size because of self-shielding at shorter UV wavelengths and increased heating at longer wavelengths. This relaxes constraints on large grain sizes suggested by near-blackbody color temperatures. Finally, from the strength of infrared [FeII] emission in the inner shell we find that the gas-phase Fe abundance is roughly solar. This is astonishing in such a dusty object, where one normally expects gaseous iron to be depleted by two orders of magnitude.
[157]  [pdf] - 119174
State-to-state rotational transitions in H$_2$+H$_2$ collisions at low temperatures
Comments: Accepted by J. Chem. Phys
Submitted: 2006-07-18
We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the super-thermal region. The two most recently published potential energy surfaces for the H$_2$-H$_2$ complex, the so-called DJ (Diep and Johnson, 2000) and BMKP (Boothroyd et al., 2002) surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H$_2$+H$_2$ collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H$_2$-H$_2$ interaction. We found significant differences in rotational excitation/de-excitation cross sections computed on the two surfaces in collisions between two para-H$_2$ molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment (Mat\'e et al., 2005) for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.
[158]  [pdf] - 82472
A Cloudy/Xspec Interface
Comments: 13 pages, 1 table, 4 figures, accepted for publication in PASP
Submitted: 2006-06-01
We discuss new functionality of the spectral simulation code CLOUDY which allows the user to calculate grids with one or more initial parameters varied and formats the predicted spectra in the standard FITS format. These files can then be imported into the x-ray spectral analysis software XSPEC and used as theoretical models for observations. We present and verify a test case. Finally, we consider a few observations and discuss our results.
[159]  [pdf] - 81279
Determining the H+ Region/PDR Equation of State in Star Forming Regions
Comments: 21 pages, 7 figures, accepted for publication in ApJ
Submitted: 2006-04-10
The emission line regions of starburst galaxies and active nuclei reveal a wealth of spectroscopic information. A unified picture of the relationship between ionized, atomic, and molecular gas makes it possible to better understand these observations. We performed a series of calculations designed to determine the equation of state, the relationship between density, temperature, and pressure, through emission-line diagnostic ratios that form in the H+ region and PDR. We consider a wide range of physical conditions in the H+ region. We connect the H+ region to the PDR by considering two constant pressure cases, one with no magnetic field and one where the magnetic field overwhelms the thermal pressure. We show that diagnostic ratios can yield the equation of state for single H+ regions adjacent to single PDRs, with the results being more ambiguous when considering observations of entire galaxies. As a test, we apply our calculations to the Orion H+/PDR region behind the Trapezium. We find the ratio of thermal to magnetic pressure in the PDR to be \~1.2. If magnetic and turbulent energy are in equipartition, our results mean the magnetic field is not the cause of the unexplained broadening in M 42, but may significantly affect line broadening in the PDR. Since Orion is often used to understand physical processes in extragalactic environments, our calculations suggest magnetic pressure should be considered in modeling such regions.
[160]  [pdf] - 78746
Physical Conditoins in Orion's Veil II: A Multi-Component Study of the Line of Sight Toward the Trapezium
Comments: 43 pages, 15 figures, to be published in ApJ
Submitted: 2005-12-21
Orion's Veil is an absorbing screen that lies along the line of sight to the Orion H II region. It consists of two or more layers of gas that must lie within a few parsecs of the Trapezium cluster. Our previous work considered the Veil as a whole and found that the magnetic field dominates the energetics of the gas in at least one component. Here we use high-resolution STIS UV spectra that resolve the two velocity components in absorption and determine the conditions in each. We derive a volume hydrogen density, 21 cm spin temperature, turbulent velocity, and kinetic temperature, for each. We combine these estimates with magnetic field measurements to find that magnetic energy significantly dominates turbulent and thermal energies in one component, while the other component is close to equipartition between turbulent and magnetic energies. We observe molecular hydrogen absorption for highly excited v, J levels that are photoexcited by the stellar continuum, and detect blueshifted S III and P III. These ions must arise from ionized gas between the mostly neutral portions of the Veil and the Trapezium and shields the Veil from ionizing radiation. We find that this layer of ionized gas is also responsible for He I absorption in the Veil, which resolves a 40-year-old debate on the origin of He I absorption towards the Trapezium. Finally, we determine that the ionized and mostly atomic layers of the Veil will collide in less than 85,000 years.
[161]  [pdf] - 77167
Physical conditions in the ISM towards HD185418
Comments: 31 pages, accepted for publication in ApJ
Submitted: 2005-10-23
We have developed a complete model of the hydrogen molecule as part of the spectral simulation code Cloudy. Our goal is to apply this to spectra of high-redshift star-forming regions where H2 absorption is seen, but where few other details are known, to understand its implication for star formation. The microphysics of H2 is intricate, and it is important to validate these numerical simulations in better-understood environments. This paper studies a well-defined line-of-sight through the Galactic interstellar medium (ISM) as a test of the microphysics and methods we use. We present a self-consistent calculation of the observed absorption-line spectrum to derive the physical conditions in the ISM towards HD185418, a line-of-sight with many observables. We deduce density, temperature, local radiation field, cosmic ray ionization rate, chemical composition and compare these conclusions with conditions deduced from analytical calculations. We find a higher density, similar abundances, and require a cosmic ray flux enhanced over the Galactic background value, consistent with enhancements predicted by MHD simulations.
[162]  [pdf] - 73956
Molecular hydrogen in the diffuse interstellar medium at high redshift
Comments: 21 pages, 12 figures, submitted to MNRAS
Submitted: 2005-06-23
Physical conditions within DLAs can reveal the star formation history, determine the chemical composition of the associated ISM, and hence document the first steps in the formation of present day galaxies. Here we present calculations that self-consistently determine the gas ionization, level populations (atomic fine-structure levels and rotational levels of H_2), grain physics, and chemistry. We show that for a low-density gas (n< 0.1 cm^-3) the meta-galactic UV background due to quasars is sufficient to maintain H_2 column densities below the detection limit (i.e N(H_2)< 10^14 cm^-2) irrespective of the metallicity and dust content in the gas. Such a gas will have a 21 cm spin temperature in excess of 7000 K and very low C I and C II* column densities for H I column densities typically observed in DLAs. We show that the observed properties of the ~15% per cent of the DLAs that do show detectable H_2 absorption cannot be reproduced with only the quasar dominated meta-galactic UV radiation field. Gas with higher densities (n>10 cm^-3), a moderate radiation field (flux density one to ten times that of the background radiation of the Galactic ISM), the observed range of metallicity and dust-to-gas ratio reproduce all the observed properties of the DLAs that show H_2 absorption lines. This favors the presence of ongoing star formation in DLAs with H_2. The absence of detectable H_2 and C I absorption in a large fraction of DLAs can be explained if they originate either in a low-density gas or in a high-density gas with a large ambient radiation field. The absence of 21 cm absorption and C II* absorption will be consistent with the first possibility. The presence of 21 cm absorption and strong C II* without H_2 and C I absorption will suggest the second alternative. (Abridged)
[163]  [pdf] - 73955
The VLT-UVES survey for molecular hydrogen in high-redshift damped Lyman-alpha systems: Physical conditions in the neutral gas
Comments: 21 pages, 16 figures, accepted for publication in MNRAS
Submitted: 2005-06-23
We study the physical conditions in DLAs, using a sample of 33 systems toward 26 QSOs acquired for a recently completed survey of H_2 by Ledoux et al. (2003). We use the column densities of H_2 in different rotational levels, together with those of C I, C I*, C I**, C II* and singly ionized atomic species to discuss the kinetic temperature, the density of hydrogen and the electronic density in the gas together with the ambient UV radiation field. Detailed comparisons are made between the observed properties in DLAs, the interstellar medium (ISM) of the Galaxy, LMC and SMC. The mean kinetic temperature of the gas corresponding to DLA subcomponents in which H_2 absorption line is detected, derived from the ortho-to-para ratio (153+/-78 K), is higher than that measured in the ISM (77+/-17 K) and the Magellanic clouds (82+/-21 K). Typical pressure in these components (corresponding to T = 100-300 K and n = 10-200 cm^-3), measured using C I fine-structure excitation, are higher than what is measured along ISM sightlines. This is consistent with the corresponding higher values for N(H_2,J=2)/N(H_2,J=0) seen in DLAs. From the column densities of the high-J rotational levels, we derive that the typical radiation field in the H_2 bearing components is of the order of or slightly higher than the mean UV field in the Galactic ISM. Determination of electron density in the gas with H_2 and C I show the ionization rate is similar to that of a cold neutral medium (CNM) in a moderate radiation field. This, together with the fact that we see H_2 in 13-20% of the DLAs, can be used to conclude that DLAs at z>1.9 could contribute as much as 50% star formation rate density seen in Lyman break galaxies (LBGs). (Abridged)
[164]  [pdf] - 73914
The H II Region/PDR Connection: Self-Consistent Calculations of Physical Conditions in Star-Forming Regions
Comments: 60 pages, 35 figures, to be published in the Astrophysical Journal. Version with full resolution is available at
Submitted: 2005-06-21
We have performed a series of calculations designed to reproduce infrared diagnostics used to determine physical conditions in star forming regions. We self-consistently calculate the thermal and chemical structure of an H II region and photodissociation region (PDR) that are in pressure equilibrium. This differs from previous work, which used separate calculations for each gas phase. Our calculations span a wide range of stellar temperatures, gas densities, and ionization parameters. We describe improvements made to the spectral synthesis code Cloudy that made these calculations possible. These include the addition of a molecular network with ~1000 reactions involving 68 molecular species and improved treatment of the grain physics. Data from the Spitzer First Look Survey, along with other archives, are used to derive important physical characteristics of the H II region and PDR. These include stellar temperatures, electron densities, ionization parameters, UV radiation flux, and PDR density. Finally, we calculate the contribution of the H II region to PDR emission line diagnostics, which allows for a more accurate determination of physical conditions in the PDR.
[165]  [pdf] - 72307
A Multi-Instrument Study of the Helix Nebula Knots with the Hubble Space Telescope
Comments: 16 pages, accepted for publication in The Astronomical Journal
Submitted: 2005-04-08
We have conducted a combined observational and theoretical investigation of the ubiquitous knots in the Helix Nebula (NGC 7293). We have constructed a combined hydrodynamic+radiation model for the ionized portion of these knots and have accurately calculated a static model for their molecular regions. Imaging observations in optical emission lines were made with the Hubble Space Telescope's STIS spectrograph, operating in a "slitless" mode, complemented by WFPC2 images in several of the same lines. The NICMOS camera was used to image the knots in molecular hydrogen. These observations, when combined with other studies of molecular hydrogen and CO provide a complete characterization of the knots. They possess dense molecular cores of densities about 1,000,000 per cubic centimeter surrounded on the central star side by a zone of hot molecular hydrogen. The temperature of the molecular hydrogen emitting layer defies explanation either through detailed calculations for radiative equilibrium or for simplistic calculations for shock excitation. Further away from the core is the ionized zone, whose peculiar distribution of emission lines is explained by the expansion effects of material flowing through this region. The shadowed region behind the core is the source of most of the CO emission from the knot and is of the low temperature expected for a radiatively heated molecular region.
[166]  [pdf] - 72180
J-resolved He I emission predictions in the low-density limit
Comments: 31 pages, 2 figures, accepted by ApJ, tentatively scheduled for 20 July 2005, v628 1
Submitted: 2005-04-04
Determinations of the primordial helium abundance are used in precision cosmological tests. These require highly accurate He I recombination rate coefficients. Here we reconsider the formation of He I recombination lines in the low-density limit. This is the simplest case and it forms the basis for the more complex situation where collisions are important. The formation of a recombination line is a two-step process, beginning with the capture of a continuum electron into a bound state and followed by radiative cascade to ground. The rate coefficient for capture from the continuum is obtained from photoionization cross sections and detailed balancing, while radiative transition probabilities determine the cascades. We have made every effort to use today's best atomic data. Radiative decay rates are from Drake's variational calculations, which include QED, fine structure, and singlet-triplet mixing. Certain high-$L$ fine-structure levels do not have a singlet-triplet distinction and the singlets and triplets are free to mix in dipole-allowed radiative decays. We use quantum defect or hydrogenic approximations to include levels higher than those treated in the variational calculations. [Abstract clipped]
[167]  [pdf] - 72015
Physical conditions in the Homunculus
Comments: 9 pages, 3 figs, ASP Conf 332
Submitted: 2005-03-29
Conditions within the Homunculus nebula around Eta Car are determined by many of the same physical processes that occur in molecular clouds in the interstellar medium. But there is one major exception -- we know when the ejection occurred and something about its composition and initial state. The gas was warm, ionized, and dust-free when it was located within the star's atmosphere and it is currently cold, molecular, and dusty. It undertook this transformation in a bit over 150 years. It offers a laboratory for the study of physical processes in a well-constrained environment. We derive a photoionization model of the Homunculus nebula that reproduces many of its observed properties. We conclude by outlining how observations of the Homunculus could address basic problems in the physics of the interstellar medium.
[168]  [pdf] - 71015
Theoretical He I Emissivities in the Case B Approximation
Comments: Accepted to ApJL
Submitted: 2005-02-10
We calculate the He I case B recombination cascade spectrum using improved radiative and collisional data. We present new emissivities over a range of electron temperatures and densities. The differences between our results and the current standard are large enough to have a significant effect not only on the interpretation of observed spectra of a wide variety of objects but also on determinations of the primordial helium abundance.
[169]  [pdf] - 70597
Molecular Hydrogen in Star-forming regions: implementation of its micro-physics in Cloudy
Comments: Accepted for publication in ApJ (58 pages). Version with full resolution is available at
Submitted: 2005-01-22
Much of the baryonic matter in the Universe is in the form of H2 which includes most of the gas in galactic and extragalactic interstellar clouds. Molecular hydrogen plays a significant role in establishing the thermal balance in many astrophysical environments and can be important as a spectral diagnostic of the gas. Modeling and interpretation of observations of such environments requires a quantitatively complete and accurate treatment of H2. Using this micro-physical model of H2, illustrative calculations of prototypical astrophysical environments are presented. This work forms the foundation for future investigations of these and other environments where H2 is an important constituent.
[170]  [pdf] - 70146
Self-Consistent Dynamic Models of Steady Ionization Fronts I. Weak-D and Weak-R Fronts
Comments: ApJ, main journal, in press 1 March 2005, v621 (22 pages, 3 appendices, 21 figures) PDF version with high-res figures available from
Submitted: 2005-01-04
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.
[171]  [pdf] - 69605
Close-coupling calculations of low-energy inelastic and elastic processes in $^4$He collisions with H$_2$: A comparative study of two potential energy surfaces
Comments: 26 pages, 11 figures
Submitted: 2004-12-09
The two most recently published potential energy surfaces (PESs) for the HeH$_2$ complex, the so-called MR (Muchnick and Russek) and BMP (Boothroyd, Martin, and Peterson) surfaces, are quantitatively evaluated and compared through the investigation of atom-diatom collision processes. The BMP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the PES compared to that of the MR surface. We found significant differences in inelastic rovibrational cross sections computed on the two surfaces for processes dominated by large changes in target rotational angular momentum. In particular, the H$_2$($\nu=1, j=0$) total quenching cross section computed using the BMP potential was found to be a factor of $~$1000 larger than that obtained with the MR surface.
[172]  [pdf] - 66246
The Origin of Fe II Emission in AGN
Comments: 23 pages. Accepted by ApJ
Submitted: 2004-07-19
We used a very large set of models of broad emission line (BEL) clouds in AGN to investigate the formation of the observed Fe II emission lines. We show that photoionized BEL clouds cannot produce both the observed shape and observed equivalent width of the 2200-2800A Fe II UV bump unless there is considerable velocity structure corresponding to a microturbulent velocity parameter v_turb > 100 km/s for the LOC models used here. This could be either microturbulence in gas that is confined by some phenomenon such as MHD waves, or a velocity shear such as in the various models of winds flowing off the surfaces of accretion disks. The alternative way that we can find to simultaneously match both the observed shape and equivalent width of the Fe II UV bump is for the Fe II emission to be the result of collisional excitation in a warm, dense gas. Such gas would emit very few lines other than Fe II. However, since the collisionally excited gas would constitute yet another component in an already complicated picture of the BELR, we prefer the model involving turbulence. In either model, the strength of Fe II emission relative to the emission lines of other ions such as Mg II depends as much on other parameters (either v_turb or the surface area of the collisionally excited gas) as it does on the iron abundance. Therefore, the measurement of the iron abundance from the FeII emission in quasars becomes a more difficult problem.
[173]  [pdf] - 64307
Physical Conditions in Orion's Veil
Comments: 45 pages, 20 figures, accepted for publication in ApJ
Submitted: 2004-04-20
Orion's veil consists of several layers of largely neutral gas lying between us and the main ionizing stars of the Orion nebula. It is visible in 21cm H I absorption and in optical and UV absorption lines of H I and other species. Toward the Trapezium, the veil has two remarkable properties, high magnetic field (~100 microGauss) and a surprising lack of molecular hydrogen given its total hydrogen column density. Here we compute photoionization models of the veil to establish its gas density and its distance from the Trapezium. We use a greatly improved model of the hydrogen molecule that determines level populations in ~1e5 rotational/vibrational levels and provides improved estimates of molecular hydrogen destruction via the Lyman-Werner bands. Our best fit photoionization models place the veil 1-3 pc in front of the star at a density of 1e3-1e4 cubic centimeters. Magnetic energy dominates the energy of non-thermal motions in at least one of the 21cm H I velocity components. Therefore, the veil is the first interstellar environment where magnetic dominance appears to exist. We find that the low ratio of molecular to atomic hydrogen (< 1e-4) is a consequence of high UV flux incident upon the veil due to its proximity to the Trapezium stars and the absence of small grains in the region.
[174]  [pdf] - 63753
The Broad Emission Line Regions of Quasars - current status and future prospects
Comments: Invited review at conference "AGN Physics with the Sloan Digital Sky Survey"
Submitted: 2004-03-23
The spectrum emitted by the Broad Emission Line Regions of Active Galactic Nuclei can measure the luminosity of the central object and the chemical composition of the interstellar medium in a very young galaxy, and so constrain the expansion of the universe and the first stages of the evolution of massive galaxies. Here I review some recent developments in our understanding of the nature of the clouds and the interpretation of their spectra.
[175]  [pdf] - 62895
Grain size distributions and photo-electric heating in ionized media
Comments: 12 pages, 3 figures, 8 tables, minor typographical corrections wrt v1
Submitted: 2004-02-16, last modified: 2004-03-23
Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. A detailed code is necessary to model grains in a photo-ionized medium. In this paper we will use the spectral synthesis code Cloudy for this purpose. The grain model of this code has recently undergone a comprehensive upgrade. Part of this is the newly developed hybrid grain charge model, which will be described in detail. This model allows discrete charge states of very small grains to be modelled accurately while simultaneously avoiding the overhead of fully resolving the charge distribution of large grains, thus making the model both accurate and computationally efficient. A comprehensive comparison with the fully resolved charge state models of Weingartner & Draine (2001) shows that the agreement is very satisfactory for realistic size distributions. The effect of the grain size distribution on the line emission from photo-ionized regions is studied by taking standard models for an H II region and a planetary nebula and adding a dust component to the models with varying grain size distributions. A comparison of the models shows that varying the size distribution has a dramatic effect on the emitted spectrum, and affects the ionization balance as well. These results clearly demonstrate that the grain size distribution is an important parameter in photo-ionization models.
[176]  [pdf] - 61915
$^{12}$C/$^{13}$C ratio in planetary nebulae from the IUE archives
Comments: ApJ accepted: 19 pages, 3 Figures, 2 Tables
Submitted: 2004-01-06
We investigated the abundance ratio of $^{12}$C/$^{13}$C in planetary nebulae by examining emission lines arising from \ion{C}{3} $2s2p ^3P_{2,1,0} \to 2s^2 ^1S_0$. Spectra were retrieved from the International Ultraviolet Explorer archives, and multiple spectra of the same object were coadded to achieve improved signal-to-noise. The $^{13}$C hyperfine structure line at 1909.6 \AA was detected in NGC 2440. The $^{12}$C/$^{13}$C ratio was found to be $\sim4.4\pm$1.2. In all other objects, we provide an upper limit for the flux of the 1910 \AA line. For 23 of these sources, a lower limit for the $^{12}$C/$^{13}$C ratio was established. The impact on our current understanding of stellar evolution is discussed. The resulting high signal-to-noise \ion{C}{3} spectrum helps constrain the atomic physics of the line formation process. Some objects have the measured 1907/1909 flux ratio outside the low-electron density theoretical limit for $^{12}$C. A mixture of $^{13}$C with $^{12}$C helps to close the gap somewhat. Nevertheless, some observed 1907/1909 flux ratios still appear too high to conform to the presently predicted limits. It is shown that this limit, as well as the 1910/1909 flux ratio, are predominantly influenced by using the standard partitioning among the collision strengths for the multiplet $^1S_0$--$^3P_J$ according to the statistical weights. A detailed calculation for the fine structure collision strengths between these individual levels would be valuable.
[177]  [pdf] - 61148
The Effects of Low-Temperature Dielectronic Recombination on the Relative Populations of the Fe M-Shell States
Comments: 10 pages, 5 figures (3 color), 2 tables. Accepted for publication in the Astrophysical Journal
Submitted: 2003-11-25
We examine the effects of low-temperature, or $\Delta$$n = 0$, dielectronic recombination (DR) on the ionization balance of the Fe M-Shell (Fe IX through Fe XVI). Since $\Delta$$n = 0$ rates are not available for these ions, we have derived estimates based on the existing rates for the first four ionization states of the CNO sequence and newly calculated rates for L-shell ions of 3rd row elements and Fe. For a range of ionization parameter and column density applicable to the intrinsic absorbers detected in {\it ASCA}, {\it Chandra}, and {\it XMM-Newton} observations of Seyfert galaxies, we generated two grids of photoionization models, with and without DR. The results show that the ionization parameter at which the population of an Fe M-shell ion peakscan increase in some cases by factor $>$ 2 when these rates are included. More importantly, there are dramatic changes in the range in ionization parameter over which individual M-shell ions contain significant fractions of the total Fe (e.g. $>$ 10%) in the plasma. These results may explain the mismatch between the range of Fe ionization states detected in the X-ray spectra of Seyferts, identified by the energies of the M-Shell Unresolved Transition Array, and those predicted by photoionization models of the X-ray absorbers that reproduce lines of second and third row elements. The results suggest that care should be taken in using 3rd and 4th row ions to constrain the physical conditions in photoionized X-ray plasmas until accurate DR rates are available. This underscores the importance of atomic physics in interpreting astronomical spectroscopy.
[178]  [pdf] - 61100
Implications of IR continua for x-ray emission/reflection in AGN
Comments: invited review at X-ray spectroscopy of AGN with Chandra and XMM-Newton, Garching, Dec 3-6, 2001, eds Th. Boller, S. Komossa, S. Kahn, H. Kunieda, MPE Report 279, p. 103
Submitted: 2003-11-22
Observations of infrared emission from AGN show that grains exist over a broad range of distances from the central object, extending to the point where they are destroyed by sublimation. These ~1000 K grains produce much of the 1 micron continuum. In this region closest to the central object there must be a gaseous component associated with the hot grains. This paper employs a state of the art grain model and shows that the gas must be very hot, with temperatures in the neighborhood of 10^6 K. The dusty component has a covering factor of roughly 50% and so this region also reprocesses much of the total x-ray emission. Our explicit models of the IR through x-ray spectral energy distributions allow the x-ray component to be predicted from IR observations. We are creating a grid of such predictions and will make them available as an XSPEC add-in, allowing this spectral component to be included in quantitative modeling of AGN spectra.
[179]  [pdf] - 59879
Grain Size Distributions and Photo-Electric Heating in Ionized Regions
Comments: 6 pages, 6 figures; to appear in "Asymmetrical Planetary Nebulae III" editors M. Meixner, J. Kastner, N. Soker, B. Balick (ASP Conf. Series)
Submitted: 2003-10-08
In this paper we present results obtained with the new grain code in Cloudy which underline the strong effect of photo-electric heating by grains in photo-ionized regions. We study the effect that the distribution of grain sizes has on the magnitude of the effect, and show that this effect is nothing short of dramatic. This makes the grain size distribution an important parameter in modeling of photo-ionized regions such as H II regions and planetary nebulae.
[180]  [pdf] - 58163
The spectrum of a single photoionized cloud
Comments: review from Quasars & Cosmology, ASP 162, Ferland & Baldwin, editors
Submitted: 2003-07-25
The emission-line spectrum of a quasar is most likely emitted by an ensemble of photoionized clouds moving with a variety of velocities, with a range of densities and distances from the central object. The state of the art in this field is to first compute the emission from a single cloud, and then prescribe a mix of clouds to reproduce line intensities, profiles, and reverberation lags. Here I review the parameters that de-fine emission from a single cloud and the outstanding problems in qua-sar emission-line analysis.
[181]  [pdf] - 53634
Temperature Variations from HST Spectroscopy of the Orion Nebula
Comments: MNRAS accepted: 30 pages, 3 Figures, 2 Tables
Submitted: 2002-12-10
We present HST/STIS long-slit spectroscopy of NGC 1976. Our goal is to measure the intrinsic line ratio [O III] 4364/5008 and thereby evaluate the electron temperature (T_e) and the fractional mean-square T_e variation (t_A^2) across the nebula. We also measure the intrinsic line ratio [N II] 5756/6585 in order to estimate T_e and t_A^2 in the N^+ region. The interpretation of the [N II] data is not as clear cut as the [O III] data because of a higher sensitivity to knowledge of the electron density as well as a possible contribution to the [N II] 5756 emission by recombination (and cascading). We present results from binning the data along the various slits into tiles that are 0.5" square (matching the slit width). The average [O III] temperature for our four HST/STIS slits varies from 7678 K to 8358 K; t_A^2 varies from 0.00682 to at most 0.0176. For our preferred solution, the average [N II] temperature for each of the four slits varies from 9133 K to 10232 K; t_A^2 varies from 0.00584 to 0.0175. The measurements of T_e reported here are an average along each line of sight. Therefore, despite finding remarkably low t_A^2, we cannot rule out significantly larger temperature fluctuations along the line of sight. The result that the average [N II] T_e exceeds the average [O III] T_e confirms what has been previously found for Orion and what is expected on theoretical grounds. Observations of the proplyd P159-350 indicate: large local extinction associated; ionization stratification consistent with external ionization by theta^1 Ori C; and indirectly, evidence of high electron density.
[182]  [pdf] - 53198
Observational consequences of fine structure line optical depths on infrared spectral diagnostics
Comments: 16 pages, 4 figures, to be published in the February 2003 issue of the PASP
Submitted: 2002-11-20
It has long been known that infrared fine structure lines of abundant ions, like the [O III] 88 micron line, can become optically thick in H II regions under certain high luminosity conditions. This could mitigate their potential as diagnostic tools, especially if the source is too dusty for optical spectroscopy to otherwise determine the system's parameters. We examined a series of photoionization calculations which were designed to push the nebulae into the limit where many IR lines should be quite optically thick. We find that radiative transfer effects do not significantly change the observed emission line spectrum. This is due to a combination of grain absorption of the hydrogen ionizing continuum and the fact that the correction for stimulated emission in these lines is large. Given these results, and the likelihood that real objects have non-thermal line broadening, it seems unlikely that line optical depth presents a problem in using these lines as diagnostics of the physical conditions or chemical composition.
[183]  [pdf] - 52199
Reliability in the Face of Complexity; The Challenge of High-End Scientific Computing
Comments: 8 pages, PDF only
Submitted: 2002-10-07
This conference centers on photoionized plasmas, and the tools necessary to understand them. One of the major goals of meetings held between developers of plasma codes is to identify sources of differences between various codes, and resolve them. Very often these differences are caused by bugs. This contribution outlines some practical aspects of developing bug-free codes, and is based on my experience creating the plasma simulation code Cloudy. Computers present nothing less than a new way to understand the universe. We can do any numerical experiment we can imagine. But machines are growing faster far more quickly than people are getting smarter. Harnessing the power of future machines in a reliable way will present an ever greater intellectual challenge. More powerful machines allow a more complete and complex simulation of physical processes, and so offers a far greater opportunity for errors. Industrial standards of quality control, and the associated code development techniques, can play a major role in developing reliable simulations.
[184]  [pdf] - 52191
Chemical Abundances in Broad Emission Line Regions: The "Nitrogen-Loud" QSO 0353-383
Comments: 21 pages, including 3 tables and 7 figures. Accepted for publication in ApJ
Submitted: 2002-10-07
The intensity of the strong N V 1240 line relative to C IV 1549 or to He II 1640 has been proposed as an indicator of the metallicity of QSO broad emission line regions, allowing abundance measurements in a large number of QSOs out to the highest redshifts. Previously, it had been shown that the (normally) much weaker lines N III] 1750 and N IV] 1486 could be used in the same way. The redshift 1.96 QSO 0353-383 has long been known to have N III] and N IV] lines that are far stronger relative to Ly-alpha or C IV than in any other QSO. Because in this particular case these intercombination lines can be easily measured, this unusual object provides an ideal opportunity for testing whether the N V line is a valid abundance indicator. Using new observations of Q0353-383 made both with HST in the ultraviolet and from the ground in the visible passband, we find that intensity ratios involving the strengths of N V, N IV] and N III] relative to lines of He, C and O all indicate that nitrogen is overabundant relative to oxygen in Q0353-383 by a factor of ~15 compared to solar ratios. This agreement among the diagnostics supports the use of these lines for measuring BLR chemical abundances. If nitrogen behaves like a secondary element, such that N/O is proportional to O/H, then the extreme nitrogen enhancement in Q0353-383 implies a metallicity of ~15 times the solar value. Even if Q0353-383 represents an extreme outlier in the N/O proportional to O/H relation, the overall metallicity should still be at least five times solar. Unusually high metallicities in Q0353-383 might imply that we caught this object just as the gas-phase metallicity in the central part of its host galaxy has peaked, at a time when the interstellar gas supply is nearly exhausted and hence the fuel source for the central QSO is ready to shut off.
[185]  [pdf] - 51741
The Mass of Quasar Broad Emission Line Regions
Comments: 12 pages, including 1 table and 2 figures. Accepted for publication in ApJ
Submitted: 2002-09-17
We show that the mass of ionized gas in the Broad Line Regions (BLRs) of luminous QSOs is at least several hundred Msun, and probably of order 10^3-10^4 M_sun. BLR mass estimates in several existing textbooks suggest lower values, but pertain to much less luminous Seyfert galaxies or include only a small fraction of the ionized/emitting volume of the BLR. The previous estimates also fail to include the large amounts of BLR gas that emit at low efficiency (in a given line), but that must be present based on reverberation and other studies. Very highly ionized gas, as well as partially ionized and neutral gas lying behind the ionization zones, are likely additional sources of mass within the BLR. The high masses found here imply that the chemical enrichment of the BLR cannot be controlled by mass ejection from one or a few stars. A significant stellar population in the host galaxies must be contributing. Simple scaling arguments based on normal galactic chemical enrichment and solar or higher BLR metallicities show that the minimum mass of the enriching stellar population is of order 10 times the BLR mass, or greater than 10^4-10^5 M_sun. More realistic models of the chemical and dynamical evolution in galactic nuclei suggest that much larger, bulge-size stellar populations are involved.
[186]  [pdf] - 48073
The physical conditions within dense cold clouds in cooling flows II
Comments: 9 pages, 7 figures. Accepted for publication in MNRAS. Minor updates to refer to Puy et al 1999, and to tidy the references
Submitted: 2002-03-05, last modified: 2002-05-13
This is a progress report on our numerical simulations of conditions in the cold cores of cooling flow condensations. The physical conditions in any non-equilibrium plasma are the result of a host of microphysical processes, many involving reactions that are research areas in themselves. We review the dominant physical processes in our previously published simulations, to clarify those issues that have caused confusion in the literature. We show that conditions in the core of an X-ray illuminated cloud are very different from those found in molecular clouds, largely because carbon remains substantially atomic and provides powerful cooling through its far infrared lines. We show how the results of the Opacity Project have had a major impact on our predictions, largely because photoionization cross sections of atoms and first ions are now calculated to be far larger than previously estimated. Finally we show that the predicted conditions are strongly affected by such complexities as microturbulence or the presence of small amounts of dust. Large masses of cold dense gas, in addition to the warmer molecular gas detected recently, could be present in cooling flows.
[187]  [pdf] - 48821
Temperature Variations from HST Imagery and Spectroscopy of NGC 7009
Comments: MNRAS accepted: 26 pages, 4 Figures, 1 Table
Submitted: 2002-04-17
We present new HST/WFPC2 imagery and STIS long-slit spectroscopy of the planetary nebula NGC 7009. The primary goal was to obtain high spatial resolution of the intrinsic line ratio [O III] 4364/5008 and thereby evaluate the electron temperature (T_e) and the fractional mean-square T_e variation (t_A^2) across the nebula. The WFPC2 T_e map is rather uniform; almost all values are between 9000 - 11,000 K, with the higher T_e's closely coinciding with the inner He^++-zone. The results indicate very small values - <~ 0.01 - for t_A^2 throughout. Our STIS data allow an even more direct determination of T_e and t_A^2, albeit for a much smaller area than with WFPC2. We present results from binning the data along the slit into tiles that are 0.5'' square (matching the slit width). The average [O III] temperature using 45 tiles (excluding the central star and STIS fiducial bars) is 10,139 K; t_A^2 is 0.0035. The measurements of T_e reported here are an average along each line of sight. Therefore, despite finding remarkably low t_A^2, we cannot completely rule out temperature fluctuations along the line of sight as the cause of the large abundance discrepancy between heavy element abundances inferred from collisionally excited emission lines compared to those derived from recombination lines.
[188]  [pdf] - 44505
Metallicities and Abundance Ratios from Quasar Broad Emission Lines
Comments: 15 pages, 5 figures in 6 files, in press with ApJ. This version fixes some minor glitches and includes previously omitted references
Submitted: 2001-09-02, last modified: 2001-11-20
The broad emission lines (BELs) of quasars and active galactic nuclei (AGNs) are important diagnostics of the relative abundances and overall metallicity in the gas. Here we present new theoretical predictions for several UV BELs. We focus specifically on the relative nitrogen abundance as a metallicity indicator, based on the expected secondary enrichment of nitrogen at metallicities Z > 0.2 Z_o. Among the lines we consider, N III] 1750/O III] 1664, N V 1240/(C IV 1549 + O VI 1034) and N V/He II 1640 are the most robust diagnostics. We argue, in particular, that the average N V BEL is not dominated by scattered Ly-alpha photons from a broad absorption line wind. We then compare our calculated line ratios with observations from the literature. The results support earlier claims that the gas-phase metallicities near quasars are typically near or several times above the solar value. We conclude that quasar activity is preceded by, or coeval with, an episode of rapid and extensive star formation in the surrounding galactic (or proto-galactic) nuclei. Chemical evolution models of these environments suggest that, to reach Z > Z_o in well-mixed interstellar gas, the star formation must have begun > 10^8 yr before the observed quasar activity.
[189]  [pdf] - 43554
Grains in Photo-Ionized Environments
Comments: To appear in the proceedings of the Lexington Plasma Meeting: 'Spectroscopic Challenges of Photoionized Plasmas', eds. Ferland G.J., Savin D.W., ASP Conference Series, (16 pages). Abstract and minor corrections added in revision
Submitted: 2001-07-10, last modified: 2001-07-12
Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. This includes the diffuse ISM, as well as H II regions, planetary nebulae, and photo-dissociation regions. A detailed code is necessary to model grains in a photo-ionized medium since the interactions of grains with their environment include a host of microphysical processes, and their importance can only be judged by performing a complete simulation. In this paper we will use the spectral synthesis code Cloudy for this purpose. A comprehensive upgrade of the grain model has been recently incorporated in Cloudy, and certain aspects of this upgrade will be discussed. Special emphasis will be on the new grain charge model. We will consider in detail the physics of grains in both ionized and neutral environments, and will present a calculation of photo-electric heating and collisional cooling rates for a range of physical conditions and grain materials and for a range of grain sizes (including a realistic size distribution). We conclude with a brief discussion of the problems currently hampering progress in this field. The new grain model will be used to model the silicate emission in the Ney-Allen nebula, and will help us better understand the nature of the grains in that part of the Orion complex.
[190]  [pdf] - 34390
The Primordial Helium Abundance: Towards Understanding and Removing the Cosmic Scatter in the dY/dZ Relation
Comments: 10 pages, 5 figures, accepted by the ApJ
Submitted: 2000-02-02
We present results from photoionization models of low-metallicity HII regions. These nebulae form the basis for measuring the primordial helium abundance. Our models show that the helium ionization correction factor (ICF) can be non-negligible for nebulae excited by stars with effective temperatures larger than 40,000 K. Furthermore, we find that when the effective temperature rises to above 45,000 K, the ICF can be significantly negative. This result is independent of the choice of stellar atmosphere. However, if an HII region has an [O III] 5007/[O I] 6300 ratio greater than 300, then our models show that, regardless of its metallicity, it will have a negligibly small ICF. A similar, but metallicity dependent, result was found using the [O III] 5007/H$\beta$ ratio. These two results can be used as selection criteria to remove nebulae with potentially non-negligible ICFs. Using our metallicity independent criterion on the data of Izotov & Thuan (1998) results in a 20% reduction of the rms scatter about the best fit $Y-Z$ line. A fit to the selected data results in a slight increase of the value of the primordial helium abundance.
[191]  [pdf] - 34329
Magnetic Confinement, MHD Waves, and Smooth Line Profiles in AGN
Comments: to appear in MNRAS
Submitted: 2000-01-31
In this paper, we show that if the broad line region clouds are in approximate energy equipartition between the magnetic field and gravity, as hypothesized by Rees, there will be a significant effect on the shape and smoothness of broad emission line profiles in active galactic nuclei. Line widths of contributing clouds or flow elements are much wider than their thermal widths, due to the presence of non-dissipative MHD waves, and their collective contribution produce emission line profiles broader and smoother than would be expected if a magnetic field were not present. As an illustration, a simple model of isotropically emitting clouds, normally distributed in velocity, is used to show that smoothness can be achieved for less than 80,000 clouds and may even be as low as a few hundred. We conclude that magnetic confinement has far reaching consequences for observing and modeling active galactic nuclei.
[192]  [pdf] - 33993
Current Development of the Photo-Ionization Code Cloudy
Comments: 4 pages, 1 figure. To appear in ``Cosmic Evolution and Galaxy Formation: Structure, Interactions, and Feedback'', ASP Conference Series, J. Franco, E. Terlevich, O. Lopez-Cruz and I. Aretxaga, eds. (extended version)
Submitted: 2000-01-11
This paper describes the most important differences between Cloudy 90 and Cloudy 94, the most recent upgrade of the code which has now been posted on the WWW. It also describes the upgrade of the grain model in Cloudy which is currently underway.
[193]  [pdf] - 109916
The collision strength of the [Ne V] infrared fine-structure lines
Comments: 4 pages, 1 figure, accepted for publication in A&A Letters
Submitted: 1999-12-13
The calculation of accurate collision strengths for atomic transitions has been a long standing problem in quantitative spectroscopy. Most modern calculations are based on the R-matrix method and problems pertaining to the use of this method have led to a discussion of the accuracy of these results. More in particular, based on an analysis of the spectra of NGC 3918 and NGC 6302, Clegg et al. (1987) and Oliva et al. (1996) have questioned R-matrix calculations for the infrared [Ne V] fine-structure transitions. Using improved flux measurements for the [Ne V] lines, we show that the conclusion that these collision strengths would be too high, is not correct. The discrepancies found by Clegg et al. (1987) can be explained by the inaccuracy of the [Ne V] 342.6 nm flux they adopted. The discrepancies found by Oliva et al. (1996) can be explained by the inaccuracy of the LRS flux for the [Ne V] 14.32 micron line. Based on the data presented in this paper there is no reason to assume that there are any problems with the R-matrix calculations for Ne^4+ of Lennon & Burke (1994). We show that the data are accurate at the 30 % level or better. This confirms the validity of the close coupling method.
[194]  [pdf] - 1943634
Properties of Dust Grains in Planetary Nebulae -- I. The Ionized Region of NGC 6445
Comments: 36 pages, 8 tables, 7 figures, accepted for publication in ApJ
Submitted: 1999-10-21
In this paper we study new infrared spectra of the evolved planetary nebula NGC 6445 obtained with ISO. These data show that the thermal emission from the grains is very cool and has a low flux compared to H beta. A model of the ionized region is constructed, using the photo-ionization code CLOUDY 90.05. Based on this model, we show from depletions in the gas phase elements that little grain destruction can have occurred in the ionized region of NGC 6445. We also argue that dust-gas separation in the nebula is not plausible. The most likely conclusion is that grains are residing inside the ionized region of NGC 6445 and that the low temperature and flux of the grain emission are caused by the low luminosity of the central star and the low optical depth of the grains. This implies that the bulk of the silicon-bearing grains in this nebula were able to survive exposure to hard UV photons for at least several thousands of years, contradicting previously published results. A comparison between optical and infrared diagnostic line ratios gives a marginal indication for the presence of a t^2-effect in the nebula. However, the evidence is not convincing and other explanations for the discrepancy are also plausible. The off-source spectrum taken with ISO-LWS clearly shows the presence of a warm cirrus component with a temperature of 24 K as well as a very cold component with a temperature of 7 K. Since our observation encompasses only a small region of the sky, it is not clear how extended the 7 K component is and whether it contributed significantly to the FIRAS spectrum taken by COBE. Because our line of sight is in the galactic plane, the very cold component could be a starless core.
[195]  [pdf] - 1469801
Emission-line Helium Abundances in Highly Obscured Nebulae
Comments: Accepted by PASP. Scheduled for the October 1999 issue. 11 pages, 5 figures
Submitted: 1999-07-27
This paper outlines a way to determine the ICF using only infrared data. We identify four line pairs, [NeIII] 36$\micron$/[NeII] 12.8$\micron$, [NeIII]~15.6$\micron$ /[NeII] 12.8$\micron$, [ArIII] 9$\micron$/[ArII] 6.9$\micron$, and [ArIII] 21$\micron$/[ArII] 6.9$\micron$, that are sensitive to the He ICF. This happens because the ions cover a wide range of ionization, the line pairs are not sensitive to electron temperature, they have similar critical densities, and are formed within the He$^+$/H$^+$ region of the nebula. We compute a very wide range of photoionization models appropriate for galactic HII regions. The models cover a wide range of densities, ionization parameters, stellar temperatures, and use continua from four very different stellar atmospheres. The results show that each line pair has a critical intensity ratio above which the He ICF is always small. Below these values the ICF depends very strongly on details of the models for three of the ratios, and so other information would be needed to determine the helium abundance. The [Ar III] 9$\micron$/[ArII] 6.9$\micron$ ratio can indicate the ICF directly due to the near exact match in the critical densities of the two lines. Finally, continua predicted by the latest generation of stellar atmospheres are sufficiently hard that they routinely produce significantly negative ICFs.
[196]  [pdf] - 106417
The PG X-ray QSO sample: Links between the UV-X-ray Continuum and Emission Lines
Comments: Contribution to ``Quasars and Cosmology'', A.S.P. Conference Series (1999), eds. G. J. Ferland and J. A. Baldwin (San Francisco: ASP); 12 pages, 4 figures
Submitted: 1999-05-07
The UV to soft X-rays of luminous AGNs dominate their bolometric luminosity, driven by an accretion-powered dynamo at the center. These photons ionize the surrounding gas, thereby providing clues to fueling and exhaust. Two sets of important relationships - neither of them understood - link the continuum and gas properties. (i) Boroson & Green's `eigenvector 1' relationships: Steeper soft X-ray spectra are clearly related to narrower Hbeta emission and stronger optical Fe II emission from the BLR, and weaker [O III] 5007 from the NLR. We show that these relationships extend to UV spectra: narrower C III] 1909, stronger low ionization lines, larger Si III] 1892/C III] 1909 (a density indicator), weaker C IV 1549 but stronger higher-ionization N V 1240. We speculate that high accretion rates are linked to high columns of dense (1e10 - 1e11 cm-3), nitrogen-enhanced, low-ionization gas from nuclear starbursts. Linewidth, inverse Fe II-[O III] and inverse Fe II-C IV relationships hint at the geometrical arrangement of this gas. (ii) The Baldwin effect (inverse equivalent width - luminosity relationships): Our correlation analyses suggest that these are independent of the above eigenvector 1 relationships. The eigenvector 1 relationships can therefore be used in future work, to reduce scatter in the Baldwin relationships, perhaps fulfilling the dream of using the Baldwin effect for cosmological studies.
[197]  [pdf] - 106120
Elemental Abundances in QSOs: Star Formation and Galactic Nuclear Evolution at High Redshifts
Comments: To appear in the 1999 Annual Review of Astronomy and Astrophysics, 47 pages, 13 figures (3-6 and 12 not embedded)
Submitted: 1999-04-16
Quasar (or QSO) elemental abundances provide unique probes of high-redshift star formation and galaxy evolution. There is growing evidence from both the emission and intrinsic absorption lines that QSO environments have roughly solar or higher metallicities out to redshifts >4. The range is not well known, but solar to a few times solar appears to be typical. There is also evidence for higher metallicities in more luminous objects, and for generally enhanced N/C and Fe/alpha abundances compared to solar ratios. These results identify QSOs with vigorous, high-redshift star formation -- consistent with the early evolution of massive galactic nuclei or dense proto-galactic clumps. However, the QSOs offer new constraints. For example, 1) most of the enrichment and star formation must occur before the QSOs "turn on" or become observable, on time scales of ~<1 Gyr at least at the highest redshifts. 2) The tentative result for enhanced Fe/alpha suggests that the first local star formation began at least ~1 Gyr prior to the QSO epoch. 3) The star formation must ultimately be extensive in order to reach high metallicities, i.e. a substantial fraction of the local gas must be converted into stars and stellar remnants. The exact fraction depends on the shape of the initial mass function (IMF). 4) The highest derived metallicities require IMFs that are weighted slightly more toward massive stars than the in solar neighborhood. 5) High metallicities also require deep gravitational potentials. By analogy with the well-known mass--metallicity relation among low-redshift galaxies, metal-rich QSOs should reside in galaxies (or proto-galaxies) that are minimally as massive (or as tightly bound) as our own Milky Way.
[198]  [pdf] - 102451
Physical Conditions in the Orion HII Region
Comments: Invited Review for 1997 Ringberg Orion Complex Meeting
Submitted: 1998-08-11
The Orion Nebula is the defining galactic HII Region, and has been well studied over a broad range of wavelengths. The gas-phase metallicity is 2/3 of solar, close to the star cluster. Grains survive throughout. The Trapezium cluster's radiation field will drive the dusty gas away, pushing it back against the photodissociation region (PDR). The nebula is in a state of energy equipartition between energy input from starlight and energy stored in gas pressure, kinetic outflow, and the magnetic field. There is much that is not clear, however. How can gas flow away from the PDR into the face of a radiation field that produces such effective radiative acceleration? Grains are photoionized by the stellar radiation field, and this charge establishes a Coulomb drag that couples the grains to the gas. The grains drift relative to the gas, with speeds of order a fraction of a km s^-1. The grains can drift an appreciable fraction of the thickness of the ionized larger in a timescale of roughly 10^5 yr. The age of the geometry is unknown since the motion of t1 C Ori is unknown. The likely equilibrium is one where grains drift away from the illuminated face of the HII region and accumulate near the PDR. The ionized gas then expands into the direction of the central star cluster, looses its emission measure, and eventually becomes invisible. The loss of photoelectric grain heating may cause a dramatic change in the temperature, and provide a natural explanation for the temperature fluctuations phenomenon. The role of the magnetic field in controlling properties of the HII region is unclear, although a component of turbulence equal to the Alfven speed has been seen. Much work remains to be done.
[199]  [pdf] - 102403
The PG X-Ray QSO Sample: Links among X-ray, UV & Optical Spectra
Comments: 7 pages with 2 figures. Needs paspconf.sty. To appear in Proc. of the Conference `Structure and Kinematics of Quasar Broad Line Regions' held 1998, March 23-26 in Lincoln, Nebraska (eds. C.M. Gaskell et al.)
Submitted: 1998-08-06
A unique, essentially complete sample of 22 QSOs, with high quality soft X-ray spectra from ROSAT, as well as HST and optical spectrophotometry from below Ly-alpha to above H-alpha, is being used to investigate the relationships among the ionizing continuum and the optical and UV continuum, emission and absorption lines. Here we present a first analysis showing that optical `Eigenvector 1' linking steeper soft X-ray spectra with increasing optical Fe II strength, decreasing [O III] 5007 emission, and narrower BLR H-beta emission, extends to the UV emission lines, and is manifested by weaker C IV 1549 emission, stronger Si III] 1892/C III] 1909 ratio, and narrower C III] 1909 emission. Steeper soft X-ray spectra have been linked to higher L/L_Edd ratios, thus apparently linking BLR densities, high and low ionization gas, and kinematics, to the accretion process.
[200]  [pdf] - 1235229
Luminosity indicators in dusty photoionized environments
Comments: 12 pages 3 figures. Accepted ASP Sept.98
Submitted: 1998-06-09
The luminosity of the central source in ionizing radiation is an essential parameter in a photoionized environment, and one of the most fundamental physical quantities one can measure. We outline a method of determining luminosity for any emission-line region using only infrared data. In dusty environments, grains compete with hydrogen in absorbing continuum radiation. Grains produce infrared emission, and hydrogen produces recombination lines. We have computed a very large variety of photoionization models, using ranges of abundances, grain mixtures, ionizing continua, densities, and ionization parameters. The conditions were appropriate for such diverse objects as H II regions, planetary nebulae, starburst galaxies, and the narrow and broad line regions of active nuclei. The ratio of the total thermal grain emission relative to H$\beta$ (IR/H$\beta$) is the primary indicator of whether the cloud behaves as a classical Str\"{o}mgren sphere (a hydrogen-bounded nebula) or whether grains absorb most of the incident continuum (a dust-bounded nebula). We find two global limits: when $IR/H\beta<100$ infrared recombination lines determine the source luminosity in ionizing photons; when $IR/H\beta\gg100$ the grains act as a bolometer to measure the luminosity.
[201]  [pdf] - 101559
Quasars as Cosmological Probes: The Ionizing Continuum, Gas Metallicity and the EW-L Relation
Comments: 11 pages, 3 figures, AASTeX aas2pp4.sty, accepted for publication in ApJ
Submitted: 1998-05-27
Using a realistic model for line emission from the broad emission line regions of quasars, we are able to reproduce the previously observed correlations of emission-line ratios with the shape of the spectral energy distribution (SED). In agreement with previous studies, we find that the primary driving force behind the Baldwin Effect (EW ~ L^beta, beta < 0) is a global change in the SED with quasar luminosity, in that more luminous quasars must have characteristically softer ionizing continua. This is completely consistent with observations that show correlations between L_uv, alpha_ox, alpha_uvx, line ratios and EWs. However, to explain the complete lack of a correlation in the EW(NV)--L_uv diagram we propose that the more luminous quasars have characteristically larger gas metallicities (Z). As a secondary element, nitrogen's rapidly increasing abundance with increasing Z compensates for the losses in EW(NV) emitted by gas illuminated by softer continua in higher luminosity quasars. A characteristic relationship between Z and L has an impact on the EW--L_uv relations for other lines as well. For a fixed SED, an increasing gas metallicity reduces the EW of the stronger metal lines (the gas cools) and that of Ly_alpha and especially HeII (because of the increasing metal opacity), while the weaker lines (e.g., CIII] 1909) generally respond positively. The interplay between the effects of a changing SED and Z with L results in the observed luminosity dependent spectral variations. All of the resulting dependences on L_uv are within the range of the observed slopes.
[202]  [pdf] - 1469571
The UV-Optical Albedo of Broad Emission Line Clouds
Comments: 15 pages, LATeX (AASTeX aaspp4.sty), 5 postscript figures, accepted for publication in The Astrophysical Journal
Submitted: 1997-10-29
We explore the effective UV-optical albedos of a variety of types of broad emission line clouds, as well as their possible effects on the observed spectra of AGN. An important albedo source in moderately ionized ionization-bounded clouds is that due to neutral hydrogen: Rayleigh scattering of continuum photons off the extreme damping wings of Lya. The photons resulting from this scattering mechanism may contribute significantly to the Lya emission line, especially in the very broad wings. In addition, line photons emitted near 1200 Angstroms (e.g., N V 1240) that stream toward the neutral portion of the cloud may be reflected off this Rayleigh scattering mirror, so that they preferentially escape from the illuminated face. Inclusion of this effect can alter predicted emission line strengths and profiles. In more highly-ionized ionization-bounded clouds, Thompson scattering dominates the UV-optical albedo, but this albedo is lessened by the hydrogen gas opacity. These clouds are most reflective on the long wavelength side of the hydrogen recombination edges; this feature may then alter the shapes of the spectral regions near the recombination edges. We illustrate the effects of gas density and line broadening on the effective albedo. We also discuss the reflection effects of the accretion disk and the ``dusty torus.'' The accretion disk is an effective reflector of UV-optical photons, whether by electron or Rayleigh scattering, and it is possible that we observe a significant fraction of this light from the AGN in reflection. This effect can alter the emission line profiles and even destroy emission at the Lyman jump emitted by broad line clouds. Finally, we discuss the possibility that continuum reflection from broad line clouds is at least in part responsible for the polarized broad absorption line troughs.
[203]  [pdf] - 97344
Physical Conditions of the Coronal Line Region in Seyfert Galaxies
Comments: 18 pages plus 6 postscript figures, to appear in the June ApJS
Submitted: 1997-05-12
The launch of the Infrared Space Observatory and new atomic data have opened a window to the study of high ionization gas in active galactic nuclei (AGN). We present the results of a large number of photoionization simulations of the ``coronal line'' region in AGN, employing new atomic data from the Opacity and Iron Projects. Our grid of line emission spans 8 orders of magnitude in gas density and 14 orders of magnitude in ionizing flux in an effort to identify the optimal conditions in which these lines form. We show that coronal lines form at distances from just outside the broad line region to ~400L_{43.5}^{1/2} pc, in gas with ionization parameter -2.0 < log U(H) < 0.75, corresponding to gas densities of 10$^2$ to 10$^{8.5}$ cm$^{-3}$, with electron temperatures ~12,000K -- 150,000K. A large range of distances from the central source implies significant line width variation among the coronal lines. We identify several line ratios which could be used to measure relative abundances, and we use these to show that the coronal line gas is likely to be dust free.
[204]  [pdf] - 97345
Locally Optimally-emitting Clouds and the Narrow Emission Lines in Seyfert Galaxies
Comments: 29 pages plus 16 figures, accepted for publication in ApJ
Submitted: 1997-05-12
The narrow emission line spectra of active galactic nuclei are not accurately described by simple photoionization models of single clouds. Recent Hubble Space Telescope images of Seyfert 2 galaxies show that these objects are rich with ionization cones, knots, filaments, and strands of ionized gas. Here we extend to the narrow line region the ``locally optimally emitting cloud'' (LOC) model, in which the observed spectra are predominantly determined by powerful selection effects. We present a large grid of photoionization models covering a wide range of physical conditions and show the optimal conditions for producing many of the strongest emission lines. We show that the integrated narrow line spectrum can be predicted by an integration of an ensemble of clouds, and we present these results in the form of diagnostic line ratio diagrams making comparisons with observations. We also predict key diagnostic line ratios as a function of distance from the ionizing source, and compare these to observations. The predicted radial dependence of the [O III]/[O II] ratio may be matched to the observed one in NGC4151, if the narrow line clouds see a more intense continuum than we see. The LOC scenario when coupled with a simple Keplerian gravitational velocity field will quite naturally predict the observed line width versus critical density relationship. The influence of dust within the ionized portion of the clouds is discussed and we show that the more neutral gas is likely to be dusty, although a high ionization dust-free region is most likely present too. This argues for a variety of NLR cloud origins.
[205]  [pdf] - 97224
Do the Broad Emission Line Clouds See the Same Continuum that We See?
Comments: 11 pages AAS-LATeX aaspp4.sty format, including 1 figure; accepted for publication in ApJ
Submitted: 1997-04-27
Recent observations of quasars, Mrk 335 and the HST quasar composite spectrum, have indicated that many of them have remarkably soft ionizing continua (fnu ~ nu^-2.0, 13.6 eV -- 100 eV). We point out that the number of E > 54.4 eV photons is insufficient to create the observed strengths of the He II emission lines. While the numbers of photons which energize C IV 1549 and O VI 1034 are sufficient, even the most efficiently emitting clouds for these two lines must each cover at least 20% -- 40% of the source. If the typical quasar ionizing continuum is indeed this soft, then we must conclude that the broad emission line clouds must see a very different (harder) continuum than we see. The other viable possibility is that the UV -- EUV SED is double-peaked, with the second peaking near 54 eV, its Wien tail the observed soft X-ray excess.
[206]  [pdf] - 96866
Charge Transfer Between Neutral Atoms and Highly Ionized Species: Implications for ISO Observations
Comments: 6 pages, latex (aaspp4.sty), 2 figures. To appear in ApJ Letters
Submitted: 1997-03-17
We estimate rate coefficients for charge transfer between neutral hydrogen and helium and moderate to highly ionized heavy elements. Although charge transfer does not have much influence on hot collisionally ionized plasmas, its effects on photoionized plasmas can be profound. We present several photoionization models which illustrate the significant effect of charge transfer on the far infrared lines detected by ISO.
[207]  [pdf] - 96014
An Atlas of Computed Equivalent Widths of Quasar Broad Emission Lines
Comments: 16 pages, LaTeX (AASTeX aaspp4.sty); to appear in the 1997 ApJS: full contents of the 9 photoionization grids presented in this paper may be found at
Submitted: 1996-11-26
We present graphically the results of several thousand photoionization calculations of broad emission line clouds in quasars, spanning seven orders of magnitude in hydrogen ionizing flux and particle density. The equivalent widths of 42 quasar emission lines are presented as contours in the particle density - ionizing flux plane for a typical incident continuum shape, solar chemical abundances, and cloud column density of $N(H) = 10^{23} cm^{-2}$. Results are similarly given for a small subset of emission lines for two other column densities ($10^{22} cm^{-2}$ and $10^{24} cm^{-2}$), five other incident continuum shapes, and a gas metallicity of 5 \Zsun. These graphs should prove useful in the analysis of quasar emission line data and in the detailed modeling of quasar broad emission line regions. The digital results of these emission line grids and many more are available over the Internet.
[208]  [pdf] - 95664
Accurate Hydrogen Spectral Simulations with a Compact Model Atom
Comments: 18 pages, prepared in MS-Word, Postscript only, 12 Figures, also available at, accepted for publication in the Astrophysical Journal
Submitted: 1996-10-18
Many large scale numerical simulations of astrophysical plasmas must also reproduce the hydrogen ionization and the resulting emission spectrum, in some cases quite accurately. We describe a compact model hydrogen atom that can be readily incorporated into such simulations. It reproduces the recombination efficiency and line spectrum predicted by much larger calculations for a very broad range of densities and temperatures. Uncertainties in hydrogen collision data are the largest source of differences between our compact atom and predictions of more extensive calculations, and underscore the need for accurate atomic data.
[209]  [pdf] - 94987
Physical Conditions in Low Ionization Regions of the Orion Nebula
Comments: 12 pages, latex (aaspp4.sty), 1 figure. To appear in ApJ Letters
Submitted: 1996-07-08
We reexamine the spectroscopic underpinnings of recent suggestions that [O I] and [Fe II] lines from the Orion H II region are produced in gas where the iron-carrying grains have been destroyed and the electron density is surprisingly high. Our new observations show that previous detections of [O I] 5577 were dominated by telluric emission. Our limits are consistent with a moderate density (10^4 cm^{-3}) photoionized gas. We show that a previously proposed model of the Orion H II region reproduces the observed [O I] and [Fe II] spectrum. These lines are fully consistent with formation in a dusty region of moderate density.
[210]  [pdf] - 110234
Atomic Data for Permitted Resonance Lines of Atoms and Ions from H to Si, and S, Ar, Ca and Fe
Comments: 6 pages of text, latex, 1 figure, 4 tables; tables in ASCII format available at or at Accepted by Atomic Data Nucl. Data Tables
Submitted: 1996-04-10
We list vacuum wavelengths, energy levels, statistical weights, transition probabilities and oscillator strengths for permitted resonance spectral lines of all ions of 18 astrophysically important elements (H through Si, S, Ar, Ca, Fe). Using a compilation of experimental energy levels, we derived accurate wavelengths for 5599 lines of 1828 ground-term multiplets which have gf-values calculated in the Opacity Project. We recalculated the Opacity Project multiplet gf-values to oscillator strengths and transition probabilities of individual lines. For completeness, we added 372 resonance lines of NeI, ArI, FeI and FeII ions which are not covered by the Opacity Project. Intercombination and forbidden lines are not included in the present compilation.
[211]  [pdf] - 93864
Atomic Data for Astrophysics. II. New Analytic Fits for Photoionization Cross Sections of Atoms and Ions
Comments: 24 pages including Postscript figures and tables, uses aaspp4.sty, accepted for publication in Astrophysical Journal. Misprint in Eq.(1) is corrected
Submitted: 1996-01-03, last modified: 1996-01-16
We present a complete set of analytic fits to the non-relativistic photoionization cross sections for the ground states of atoms and ions of elements from H through Si, and S, Ar, Ca, and Fe. Near the ionization thresholds, the fits are based on the Opacity Project theoretical cross sections interpolated and smoothed over resonances. At higher energies, the fits reproduce calculated Hartree-Dirac-Slater photoionization cross sections.
[212]  [pdf] - 354164
Locally Optimally Emitting Clouds and the Origin of Quasar Emission Lines
Comments: 9 pages including 1 ps figure. LaTeX format using aaspp4.sty
Submitted: 1995-10-16
The similarity of quasar line spectra has been taken as an indication that the emission line clouds have preferred parameters, suggesting that the environment is subject to a fine tuning process. We show here that the observed spectrum is a natural consequence of powerful selection effects. We computed a large grid of photoionization models covering the widest possible range of cloud gas density and distance from the central continuum source. For each line only a narrow range of density and distance from the continuum source results in maximum reprocessing efficiency, corresponding to ``locally optimally-emitting clouds'' (LOC). These parameters depend on the ionization and excitation potentials of the line, and its thermalization density. The mean QSO line spectrum can be reproduced by simply adding together the full family of clouds, with an appropriate covering fraction distribution. The observed quasar spectrum is a natural consequence of the ability of various clouds to reprocess the underlying continuum, and can arise in a chaotic environment with no preferred pressure, gas density, or ionization parameter.
[213]  [pdf] - 93297
Atomic data for astrophysics. I. Radiative recombination rates for H-like, He-like, Li-like and Na-like ions over a broad range of temperature
Comments: 15 pages, latex using aaspp.sty, 2 tables and 8 figures added as uufiles output. Accepted for publication in Astrophysical Journal Suppl
Submitted: 1995-09-15
We present new calculations and analytic fits to the rates of radiative recombination towards H-like, He-like, Li-like and Na-like ions of all elements from H through Zn (Z=30). The fits are valid over a wide range of temperature, from 3 K to 10**9 K.