Source author record

G. J. Ferland

G. J. Ferland appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph.HE astro-ph.IM physics.atom-ph astro-ph physics.space-ph physics.plasm-ph physics.acc-ph astro-ph.EP hep-ex hep-ph nucl-ex nucl-th physics.chem-ph

Catalog footprint

What is connected

65works

16topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A recent update of gas-phase chemical reactions and molecular lines in CLOUDY: its effects on millimeter and sub-millimeter molecular line predictions

Here we present our current updates of the gas-phase chemical reaction rates and molecular lines in the spectral synthesis code CLOUDY, and its implications in spectroscopic modelling of various astrophysical environments. We include energy levels, radiative and collisional rates for HF, CF$^+$, HC$_3$N, ArH$^+$, HCl, HCN, CN, CH, and CH$_2$. Simultaneously, we expand our molecular network involving these molecules. For this purpose, we have added 561 new reactions and have updated the existing 165 molecular reaction rates involving these molecules. As a result, CLOUDY now predicts all the lines arising from these nine molecules. In addition, we also update H$_2$--H$_2$ collisional data up to rotational levels $J$=31 for $v$=0. We demonstrate spectroscopic simulations of these molecules for a few astrophysical environments. Our existing model for globules in the Crab nebula successfully predicts the observed column density of ArH$^+$. Our model predicts a detectable amount of HeH$^+$, OH$^+$, and CH$^+$ for the Crab nebula. We also model the ISM towards HD185418, W31C, NGC 253, and our predictions match with most of the observed column densities within the observed error bars. Very often molecular lines trace various physical conditions. Hence, this update will be very supportive for spectroscopic modelling of various astrophysical environments, particularly involving sub-millimeter and mid-infrared observations using ALMA and JWST, respectively.

preprint2022arXiv

Hidden Cooling Flows in Clusters of Galaxies

The radiative cooling time of the hot gas at the centres of cool cores in clusters of galaxies drops down to 10 million years and below. The observed mass cooling rate of such gas is very low, suggesting that AGN feedback is very tightly balanced or that the soft X-ray emission from cooling is somehow hidden from view. We use an intrinsic absorption model in which the cooling and coolest gas are closely interleaved to search for hidden cooling flows in the Centaurus, Perseus and A1835 clusters of galaxies. We find hidden mass cooling rates of between 10 to 500 Msunpyr as the cluster mass increases, with the absorbed emission emerging in the Far Infrared band. Good agreement is found between the hidden cooling rate and observed FIR luminosity in the Centaurus Cluster. The limits on the other two clusters allow for considerable hidden cooling. The implied total mass of cooled gas is much larger than the observed molecular masses. We discuss its fate including possible further cooling and collapse into undetected very cold clouds, low mass stars and substellar objects,

preprint2022arXiv

Linking the properties of accreting white dwarfs with the ionization state of their ambient medium

Steadily accreting white dwarfs (WDs) are efficient sources of ionization and thus, are able to create extended ionized nebulae in their vicinity. These nebulae represent ideal tools for the detection of accreting WDs, given that in most cases the source itself is faint. In this work, we combine radiation transfer simulations with known H and He accreting WD models, providing for the first time the ionization state and the emission line spectra of the formed nebulae as a function of the WD mass, the accretion rate and the chemical composition of the accreted material. We find that the nebular optical line fluxes and radial extent vary strongly with the WD's accretion properties, peaking in systems with WD masses of 0.8 - 1.2 $\rm~M_{\odot}$. Projecting our results on the 'BPT' diagnostic diagrams, we show that accreting WDs nebulae possess characteristics distinct from those of H II-like regions, while they share similar line ratios with the galactic low-ionization emission-line regions. Finally, we compare our results to the relevant constraints imposed by the lack of ionized nebulae in the vicinity of supersoft X-ray sources (SSSs) and Type Ia supernova remnants - sources which are related to steadily accreting WDs. The large discrepancies uncovered by our comparison rule out any steadily accreting WD as a potential progenitor of the studied remnants and additionally require the ambient medium around the SSSs to be less dense than 0.2 $\rm~cm^{-3}$. We discuss possible alternatives that could bridge the incompatibility between the theoretical expectations and the relevant observations.

preprint2021arXiv

Deciphering the 3-D Orion Nebula-III: Structure on the NE boundary of the Orion-S Embedded Molecular Cloud

We have extended the work of Papers I and II of this series to determine at higher spatial resolution the properties of the embedded Orion-S Molecular Cloud that lies within the ionized cavity of the Orion Nebula and of the thin ionized layer that lies between the Cloud and the observer. This was done using existing and new [NII](658.3 nm) and [OIII] (500.7 nm) spectra that map the central region of the Orion Nebula (the Huygens Region). However, it remains unclear how the surface brightness of the ionized layer on the Orion-S Molecular Cloud and that of a foreground Nearer Ionized Layer are linked, as the observations show they must be. It is shown that the Cloud modifies the outer parts of the Huygens Region in the direction of the extended hot X-ray gas.

preprint2021arXiv

Dense Molecular Clouds in the Crab Supernova Remnant

Molecular emission was imaged with ALMA from numerous components near and within bright H2-emitting knots and absorbing dust globules in the Crab Nebula. These observations provide a critical test of how energetic photons and particles produced in a young supernova remnant interact with gas, cleanly differentiating between competing models. The four fields targeted show contrasting properties but within them, seventeen distinct molecular clouds are identified with CO emission; a few also show emission from HCO+, SiO and/or SO. These observations are compared with Cloudy models of these knots. It has been suggested that the Crab filaments present an exotic environment in which H2 emission comes from a mostly-neutral zone probably heated by cosmic rays produced in the supernova surrounding a cool core of molecular gas. Our model is consistent with the observed CO J=3-2 line strength. These molecular line emitting knots in the Crab present a novel phase of the ISM representative of many important astrophysical environments.

preprint2021arXiv

Role of Polycyclic Aromatic Hydrocarbons on the Cosmic-Ray ionization rate in the Galaxy

The cosmic-ray ionization rate ($ζ$, s$^{-1}$) plays an important role in the interstellar medium. It controls ion-molecular chemistry and provides a source of heating. Here we perform a grid of calculations using the spectral synthesis code CLOUDY along nine sightlines towards, HD 169454, HD 110432, HD 204827, $λ$ Cep, X Per, HD 73882, HD 154368, Cyg OB2 5, Cyg OB2 12. The value of $ζ$ is determined by matching the observed column densities of H$_3^+$ and H$_2$. The presence of polycyclic aromatic hydrocarbons (PAHs) affects the free electron density, which changes the H$_3^+$ density and the derived ionization rate. PAHs are ubiquitous in the Galaxy, but there are also regions where PAHs do not exist. Hence, we consider clouds with a range of PAH abundances and show their effects on the H$_3^+$ abundance. We predict an average cosmic-ray ionization rate for H$_2$ ($ζ$(H$_2$))= (7.88 $\pm$ 2.89) $\times$ 10$^{-16}$ s$^{-1}$ for models with average Galactic PAHs abundances, (PAH/H =10$^{-6.52}$), except Cyg OB2 5 and Cyg OB2 12. The value of $ζ$ is nearly 1 dex smaller for sightlines toward Cyg OB2 12. We estimate the average value of $ζ$(H$_2$)= (95.69 $\pm$ 46.56) $\times$ 10$^{-16}$ s$^{-1}$ for models without PAHs.

preprint2020arXiv

2018 NASA Laboratory Astrophysics Workshop: Scientific Organizing Committee Report

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.

preprint2020arXiv

Current and future development of the photoionization code Cloudy

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.

preprint2020arXiv

Deciphering the 3-D Orion Nebula-I: Expanding Shells in the Huygens Region

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.

preprint2020arXiv

Deciphering the 3-D Orion Nebula-II: A low-ionization region of multiple velocity components southwest of Theta1OriC confounds interpretation of low velocity resolution studies of temperature, density, and abundance

We establish that there are two velocity systems along lines-of-sight that contribute to the emission-line spectrum of the the brightest parts of the Orion Nebula. These overlie the Orion-S embedded molecular cloud southwest of the dominant ionizing star (Theta1OriC). Examination of 10x10" samples of high spectral resolution emission-line spectra of this region reveals it to be of low ionization, with velocities and ionization different from the central part of the Nebula. These properties jeopardize earlier determinations of abundance and physical conditions since they indicate that this region is much more complex than has been assumed in analyzing earlier spectroscopic studies and argue for use of very high spectral resolution or known simple regions in future studies.

preprint2020arXiv

Improved Fe II emission line models for AGN using new atomic datasets

Understanding the Fe II emission from Active Galactic Nuclei (AGN) has been a grand challenge for many decades. The rewards from understanding the AGN spectra would be immense, involving both quasar classification schemes such as "Eigenvector 1" and tracing the chemical evolution of the cosmos. Recently, three large Fe II atomic datasets with radiative and electron collisional rates have become available. We have incorporated these into the spectral synthesis code Cloudy and examine predictions using a new generation of AGN Spectral Energy Distribution (SED), which indicates that the UV emission can be quite different depending on the dataset utilized. The Smyth et al dataset better reproduces the observed Fe II template of the I ZW 1 Seyfert galaxy in the UV and optical regions, and we adopt these data. We consider both thermal and microturbulent clouds and show that a microturbulence of $\approx$ 100 km/s reproduces the observed shape and strength of the so-called Fe II "UV bump". Comparing our predictions with the observed Fe II template, we derive a typical cloud density of $10^{11}$ cm$^{-3}$ and photon flux of $10^{20}$ cm$^{-2}$ s$^{-1}$, and show that these largely reproduce the observed Fe II emission in the UV and optical. We calculate the $I$(Fe II)/$I$(Mg II) emission-line intensity ratio using our best-fitting model and obtain log($I$(Fe II)/$I$(Mg II)) $\sim$ 0.7, suggesting many AGNs have a roughly solar Fe/Mg abundance ratio. Finally, we vary the Eddington ratio and SED shape as a step in understanding the Eigenvector 1 correlation.

preprint2020arXiv

Physical conditions in two high-redshift H$_2$-bearing GRB-DLAs, 120815A and 121024A

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-$α$ 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.

preprint2020arXiv

Space Telescope and Optical Reverberation Mapping Project. IX. Velocity-Delay Maps for Broad Emission Lines in NGC 5548

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.

preprint2020arXiv

Space Telescope and Optical Reverberation Mapping Project. XI. Disk-wind characteristics and contributions to the very broad emission lines of NGC 5548

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.

preprint2020arXiv

Space Telescope and Optical Reverberation Mapping Project. XIII. An Atlas of UV and X-ray Spectroscopic Signatures of theDisk Wind in NGC 5548

The unusual behavior of the spectral lines of NGC5548 during the STORM campaign demonstrated a missing piece in the structure of AGNs. For a two-month period in the middle of the campaign, the spectral lines showed a deficit in flux and a reduced response to the variations of the UV continuum. This was the first time that this behavior was unequivocally observed in an AGN. Our previous papers explained this as being due to a variable disk-wind which acts as a shield and alters the SED. Here we use Cloudy to create an atlas of photoionization models for a variety of disk-winds to study their effects on the SED. We show that the winds have three different cases: Case 1 winds are transparent, fully ionized and have minimal effects on the intrinsic SED, although they can produce some line emission, especially HeIIor FeKα. We propose that this is the situation in most of the AGNs. Case 2 winds have a He++-He+ionization-front, block part of the XUV continuum but transmit much of the Lyman continuum. They lead to the observed abnormal behavior. Case 3 winds have H+ionization-front and block much of the Lyman continuum. The results show that the presence of the winds has important effects on the spectral lines of AGNs. They will thus have an effect on the measurements of the black hole mass and the geometry of the AGN. This atlas of spectral simulations can serve as a guide to future reverberation campaigns

preprint2016arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

HST imaging of the dusty filaments and nucleus swirl in NGC4696 at the centre of the Centaurus Cluster

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.

preprint2016arXiv

Space Telescope and Optical Reverberation Mapping Project. IV. Anomalous behavior of the broad ultraviolet emission lines in NGC 5548

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.

preprint2016arXiv

Structure and Physical Conditions in the Huygens Region of the Orion Nebula

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.

preprint2016arXiv

The abundance discrepancy factor and t^2 in nebulae: are non-thermal electrons the culprits?

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.

preprint2016arXiv

The intermediate line region in active galactic nuclei

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Bootstrapping dielectronic recombination from second-row elements and the Orion Nebula

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.

preprint2015arXiv

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

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

preprint2015arXiv

Stout: Cloudy's Atomic and Molecular Database

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.

preprint2015arXiv

The Nature and Frequency of Outflows from Stars in the Central Orion Nebula Cluster

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.

preprint2014arXiv

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

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.

preprint2014arXiv

Implications of Coronal Line Emission in NGC 4696

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

preprint2014arXiv

Radiative Cooling II: Effects of Density and Metallicity

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.

preprint2013arXiv

Cold gas dynamics in Hydra-A: evidence for a rotating disk

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 $μ$m and [OI]63 $μ$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.

preprint2013arXiv

Corrigendum: Improved He~I Emissivities in the Case B Approximation

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.

preprint2013arXiv

Detecting the rapidly expanding outer shell of the Crab Nebula: where to look

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] $λ$7.652 $μ$m line is a particularly good case; it should be dramatically brighter than the optical lines commonly used in searches. The C IV $\lambda1549Å$ 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.

preprint2013arXiv

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

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.

preprint2013arXiv

Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--I: Structure and Evolution

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.

preprint2013arXiv

Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--II:Physical Conditions

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.

preprint2013arXiv

Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--III:Tangential Motions using AstroDrizzle Images

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.

preprint2013arXiv

Suppression of Dielectronic Recombination Due to Finite Density Effects

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.

preprint2013arXiv

The 2013 Release of Cloudy

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.

preprint2013arXiv

The Nature of the H2-Emitting Gas in the Crab Nebula

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Improved He I Emissivities in the Case B Approximation

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.

preprint2012arXiv

Multiphase Signatures of AGN Feedback in Abell 2597

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.

preprint2012arXiv

Pumping up the [N I] nebular lines

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

preprint2012arXiv

Residual Cooling and Persistent Star Formation amid AGN Feedback in Abell 2597

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.

preprint2012arXiv

Star formation in self-gravitating disks in active galactic nuclei. II. Episodic formation of broad line regions

(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.

preprint2012arXiv

The Curious Conundrum Regarding Sulfur Abundances In Planetary Nebulae

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.

preprint2011arXiv

A deep spectroscopic study of the filamentary nebulosity in NGC4696, the brightest cluster galaxy in the Centaurus cluster

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.

preprint2011arXiv

A Survey of Molecular Hydrogen in the Crab Nebula

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.

preprint2011arXiv

H2 Temperatures in the Crab Nebula

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.

preprint2011arXiv

Physical Conditions in Barnard's Loop, Components of the Orion-Eridanus Bubble, and Implications for the WIM Component of the ISM

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α) versus I([NII])/I(Hα)) 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α) 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α Mapper are of lower temperatures than the classical bright HII regions.

preprint2011arXiv

Star formation in self-gravitating disks in active galactic nuclei. I. Metallicity gradients in broad line regions

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.

preprint2011arXiv

The curious conundrum regarding sulfur and oxygen abundances in planetary nebulae

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.

preprint2011arXiv

The energy source of the filaments around the giant galaxy NGC1275

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.

preprint2011arXiv

The Impact of Recent Advances in Laboratory Astrophysics on our Understanding of the Cosmos

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.

preprint2010arXiv

A Bright Molecular Core in a Crab Nebula Filament

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.

preprint2010arXiv

Detection of optical coronal emission from 10^6 K gas in the core of the Centaurus cluster

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.

preprint2010arXiv

Herschel images of NGC 6720: H2 formation on dust grains

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.

preprint2010arXiv

Herschel observations of FIR emission lines in brightest cluster galaxies

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.

preprint2010arXiv

Herschel photometry of brightest cluster galaxies in cooling flow clusters

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.

preprint2009arXiv

Electron impact excitation of OII fine-structure levels

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.

preprint2008arXiv

Fluorescent excitation of Balmer lines in gaseous nebulae: case D

(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.

preprint2008arXiv

Orion's Bar: Physical Conditions across the Definitive H+ / H0 / H2 Interface

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.

preprint2007arXiv

A PDR-Code Comparison Study

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.

preprint2007arXiv

Optically Thick [O I] and [C II] Emission Towards NGC 6334 A

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.

G. J. Ferland

What is connected

Connect this record

See the researcher in context

Building this map preview

65 published item(s)

A recent update of gas-phase chemical reactions and molecular lines in CLOUDY: its effects on millimeter and sub-millimeter molecular line predictions

Hidden Cooling Flows in Clusters of Galaxies

Linking the properties of accreting white dwarfs with the ionization state of their ambient medium

Deciphering the 3-D Orion Nebula-III: Structure on the NE boundary of the Orion-S Embedded Molecular Cloud

Dense Molecular Clouds in the Crab Supernova Remnant

Role of Polycyclic Aromatic Hydrocarbons on the Cosmic-Ray ionization rate in the Galaxy

2018 NASA Laboratory Astrophysics Workshop: Scientific Organizing Committee Report

Current and future development of the photoionization code Cloudy

Deciphering the 3-D Orion Nebula-I: Expanding Shells in the Huygens Region

Deciphering the 3-D Orion Nebula-II: A low-ionization region of multiple velocity components southwest of Theta1OriC confounds interpretation of low velocity resolution studies of temperature, density, and abundance

Improved Fe II emission line models for AGN using new atomic datasets

Physical conditions in two high-redshift H$_2$-bearing GRB-DLAs, 120815A and 121024A

Space Telescope and Optical Reverberation Mapping Project. IX. Velocity-Delay Maps for Broad Emission Lines in NGC 5548

Space Telescope and Optical Reverberation Mapping Project. XI. Disk-wind characteristics and contributions to the very broad emission lines of NGC 5548

Space Telescope and Optical Reverberation Mapping Project. XIII. An Atlas of UV and X-ray Spectroscopic Signatures of theDisk Wind in NGC 5548

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

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

HST imaging of the dusty filaments and nucleus swirl in NGC4696 at the centre of the Centaurus Cluster

Space Telescope and Optical Reverberation Mapping Project. IV. Anomalous behavior of the broad ultraviolet emission lines in NGC 5548

Structure and Physical Conditions in the Huygens Region of the Orion Nebula

The abundance discrepancy factor and t^2 in nebulae: are non-thermal electrons the culprits?

The intermediate line region in active galactic nuclei

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

Bootstrapping dielectronic recombination from second-row elements and the Orion Nebula

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

Stout: Cloudy's Atomic and Molecular Database

The Nature and Frequency of Outflows from Stars in the Central Orion Nebula Cluster

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

Implications of Coronal Line Emission in NGC 4696

Radiative Cooling II: Effects of Density and Metallicity

Cold gas dynamics in Hydra-A: evidence for a rotating disk

Corrigendum: Improved He~I Emissivities in the Case B Approximation

Detecting the rapidly expanding outer shell of the Crab Nebula: where to look

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

Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--I: Structure and Evolution

Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--II:Physical Conditions

Studies of NGC 6720 with Calibrated HST WFC3 Emission-Line Filter Images--III:Tangential Motions using AstroDrizzle Images

Suppression of Dielectronic Recombination Due to Finite Density Effects

The 2013 Release of Cloudy

The Nature of the H2-Emitting Gas in the Crab Nebula

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

Improved He I Emissivities in the Case B Approximation

Multiphase Signatures of AGN Feedback in Abell 2597

Pumping up the [N I] nebular lines

Residual Cooling and Persistent Star Formation amid AGN Feedback in Abell 2597

Star formation in self-gravitating disks in active galactic nuclei. II. Episodic formation of broad line regions

The Curious Conundrum Regarding Sulfur Abundances In Planetary Nebulae

A deep spectroscopic study of the filamentary nebulosity in NGC4696, the brightest cluster galaxy in the Centaurus cluster

A Survey of Molecular Hydrogen in the Crab Nebula

H2 Temperatures in the Crab Nebula

Physical Conditions in Barnard's Loop, Components of the Orion-Eridanus Bubble, and Implications for the WIM Component of the ISM

Star formation in self-gravitating disks in active galactic nuclei. I. Metallicity gradients in broad line regions

The curious conundrum regarding sulfur and oxygen abundances in planetary nebulae

The energy source of the filaments around the giant galaxy NGC1275

The Impact of Recent Advances in Laboratory Astrophysics on our Understanding of the Cosmos

A Bright Molecular Core in a Crab Nebula Filament

Detection of optical coronal emission from 10^6 K gas in the core of the Centaurus cluster

Herschel images of NGC 6720: H2 formation on dust grains

Herschel observations of FIR emission lines in brightest cluster galaxies

Herschel photometry of brightest cluster galaxies in cooling flow clusters

Electron impact excitation of OII fine-structure levels

Fluorescent excitation of Balmer lines in gaseous nebulae: case D

Orion's Bar: Physical Conditions across the Definitive H+ / H0 / H2 Interface

A PDR-Code Comparison Study

Optically Thick [O I] and [C II] Emission Towards NGC 6334 A