Source author record

Gary J. Ferland

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

31 published item(s)

preprint2024arXiv

AGN STORM 2: V. Anomalous Behavior of the CIV Light Curve in Mrk 817

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response.

preprint2022arXiv

Multiple locations of near-infrared coronal lines in NGC 5548

We present the first intensive study of the variability of the near-infrared coronal lines in an active galactic nucleus (AGN). We use data from a one-year long spectroscopic monitoring campaign with roughly weekly cadence on NGC 5548 to study the variability in both emission line fluxes and profile shapes. We find that in common with many AGN coronal lines, those studied here are both broader than the low-ionisaton forbidden lines and blueshifted relative to them, with a stratification that implies an origin in an outflow interior to the standard narrow line region. We observe for the first time [S VIII] and [Si VI] coronal line profiles that exhibit broad wings in addition to narrow cores, features not seen in either [S IX] or [Si X]. These wings are highly variable, whereas the cores show negligible changes. The differences in both the profile shapes and variability properties of the different line components indicate that there are at least two coronal line regions in AGN. We associate the variable, broad wings with the base of an X-ray heated wind evaporated from the inner edge of the dusty torus. The coronal line cores may be formed at several locations interior to the narrow line region: either along this accelerating, clumpy wind or in the much more compact outflow identified with the obscurer and so emerging on scales similar to the outer accretion disc and broad line region.

preprint2022arXiv

The Paschen Jump as a Diagnostic of the Diffuse Nebular Continuum Emission in Active Galactic Nuclei

Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 Å is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope STIS spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multi-component spectral models over the range $6800-9700$ Å and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from $\sim$10% to 50% at 8000 Å, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe II emission, high-order Balmer line emission, and other components.

preprint2020arXiv

Fe III emission in quasars: evidence for a dense turbulent medium

Recent improvements to atomic energy-level data allow, for the first time, accurate predictions to be made for the Fe III line emission strengths in the spectra of luminous, $L_\text{bol}=10^{46}-10^{48}$ erg/s, Active Galactic Nuclei. The Fe III emitting gas must be primarily photoionized, consistent with observations of line reverberation. We use CLOUDY models exploring a wide range of parameter space, together with 26,500 rest-frame ultraviolet spectra from the Sloan Digital Sky Survey, to constrain the physical conditions of the line emitting gas. The observed Fe III emission is best accounted for by dense ($n_H=10^{14}$ cm$^{-3}$) gas which is microturbulent, leading to smaller line optical depths and fluorescent excitation. Such high density gas appears to be present in the central regions of the majority of luminous quasars. Using our favoured model, we present theoretical predictions for the relative strengths of the Fe III UV34 $λλ$1895,1914,1926 multiplet. This multiplet is blended with the Si III] $λ$1892 and C III] $λ$1909 emission lines and an accurate subtraction of UV34 is essential when using these lines to infer information about the physics of the broad line region in quasars.

preprint2020arXiv

Newly Improved Ionization Corrections for the Neutral Interstellar Medium: Enabling Accurate Abundance Determinations in Star-forming Galaxies throughout the Universe

Studies measuring the chemical abundances of the neutral gas in star-forming galaxies (SFGs) require ionization correction factors (ICFs) to accurately measure their metal contents. In the work presented here we calculate newly improved ICFs for a sample of SFGs. These new corrections include both the contaminating ionized gas along the line of sight (ICF$_{\rm ionized}$) and unaccounted higher ionization stages in the neutral gas (ICF$_{\rm neutral}$). We make use of recently acquired spectroscopic observations taken with the Cosmic Origins Spectrograph (COS) on board Hubble to measure column densities for Fe II and Fe III. Using the Fe III/Fe II ratios as well as other physical properties (i.e. $\log$[L$_{\rm UV}$], $N$(H I), T, and $Z$) we generate ad-hoc photoionization models with CLOUDY to quantify the corrections required for each of the targets. We identify a luminosity threshold of $\log$[L$_{\rm UV}$]$\sim$ 40.75 erg s$^{-1}$ above which the ICF$_{\rm neutral}$ values for nitrogen are relatively higher (ICF$_{\rm neutral}=0.05$-0.7) than those for the rest of the elements (ICF$_{\rm neutral}\sim 0.01$). This behavior indicates that for the high UV luminosity objects, N II is found in non-negligible quantities in the neutral gas, making these ICF$_{\rm neutral}$ corrections critical for determining the true abundances in the interstellar medium. In addition, we calculate ICFs from a uniform grid of models covering a wide range of physical properties typically observed in studies of SFGs and extragalactic H II regions. We provide the community with tabulated ICF values for the neutral gas abundances measured from a variety of environments and applicable to chemical studies of the high redshift universe.

preprint2016arXiv

Interpreting the Ionization Sequence in Star-Forming Galaxy Emission-Line Spectra

High ionization star forming (SF) galaxies are easily identified with strong emission line techniques such as the BPT diagram, and form an obvious ionization sequence on such diagrams. We use a locally optimally emitting cloud model to fit emission line ratios that constrain the excitation mechanism, spectral energy distribution, abundances and physical conditions along the star-formation ionization sequence. Our analysis takes advantage of the identification of a sample of pure star-forming galaxies, to define the ionization sequence, via mean field independent component analysis. Previous work has suggested that the major parameter controlling the ionization level in SF galaxies is the metallicity. Here we show that the observed SF- sequence could alternatively be interpreted primarily as a sequence in the distribution of the ionizing flux incident on gas spread throughout a galaxy. Metallicity variations remain necessary to model the SF-sequence, however, our best models indicate that galaxies with the highest and lowest observed ionization levels (outside the range -0.37 < log [O III]/H\b{eta} < -0.09) require the variation of an additional physical parameter other than metallicity, which we determine to be the distribution of ionizing flux in the galaxy.

preprint2016arXiv

Ultraviolet emission lines of Si II in quasars --- investigating the "Si II disaster"

The observed line intensity ratios of the Si II 1263 and 1307 Å multiplets to that of Si II 1814\,Å in the broad line region of quasars are both an order of magnitude larger than the theoretical values. This was first pointed out by Baldwin et al. (1996), who termed it the "Si II disaster", and it has remained unresolved. We investigate the problem in the light of newly-published atomic data for Si II. Specifically, we perform broad line region calculations using several different atomic datasets within the CLOUDY modeling code under optically thick quasar cloud conditions. In addition, we test for selective pumping by the source photons or intrinsic galactic reddening as possible causes for the discrepancy, and also consider blending with other species. However, we find that none of the options investigated resolves the Si II disaster, with the potential exception of microturbulent velocity broadening and line blending. We find that a larger microturbulent velocity ($\sim 500 \rm \, kms^{-1}$) may solve the Si II disaster through continuum pumping and other effects. The CLOUDY models indicate strong blending of the Si II 1307 Å multiplet with emission lines of O I, although the predicted degree of blending is incompatible with the observed 1263/1307 intensity ratios. Clearly, more work is required on the quasar modelling of not just the Si II lines but also nearby transitions (in particular those of O I) to fully investigate if blending may be responsible for the Si II disaster.

preprint2015arXiv

Atomic data for Zn II - Improving Spectral Diagnostics of Chemical Evolution in High-redshift Galaxies

Damped Lyman-alpha (DLA) and sub-DLA absorbers in quasar spectra provide the most sensitive tools for measuring element abundances of distant galaxies. Estimation of abundances from absorption lines depends sensitively on the accuracy of the atomic data used. We have started a project to produce new atomic spectroscopic parameters for optical/UV spectral lines using state-of-the-art computer codes employing very broad configuration interaction basis. Here we report our results for Zn II, an ion used widely in studies of the interstellar medium (ISM) as well as DLA/sub-DLAs. We report new calculations of many energy levels of Zn II, and the line strengths of the resulting radiative transitions. Our calculations use the configuration interaction approach within a numerical Hartree-Fock framework. We use both non-relativistic and quasi-relativistic one-electron radial orbitals. We have incorporated the results of these atomic calculations into the plasma simulation code Cloudy, and applied them to a lab plasma and examples of a DLA and a sub-DLA. Our values of the Zn II λλ 2026, 2062 oscillator strengths are higher than previous values by 0.10 dex. Cloudy calculations for representative absorbers with the revised Zn atomic data imply ionization corrections lower than calculated before by 0.05 dex. The new results imply Zn metallicities should be lower by 0.1 dex for DLAs and by 0.13-0.15 dex for sub-DLAs than in past studies. Our results can be applied to other studies of Zn II in the Galactic and extragalactic ISM.

preprint2015arXiv

Herschel dust emission as a probe of starless cores mass: MCLD 123.5+24.9 of the Polaris Flare

We present newly processed archival Herschel images of molecular cloud MCLD 123.5+24.9 in the Polaris Flare. This cloud contains five starless cores. Using the spectral synthesis code Cloudy, we explore uncertainties in the derivation of column densities, hence, masses of molecular cores from Herschel data. We first consider several detailed grain models that predict far-IR grain opacities. Opacities predicted by the models differ by more than a factor of two, leading to uncertainties in derived column densities by the same factor. Then we consider uncertainties associated with the modified blackbody fitting process used by observers to estimate column densities. For high column density clouds (N(H) $\gg$ 10$^{22}$ cm$^{-2}$), this fitting technique can underestimate column densities by about a factor of three. Finally, we consider the virial stability of the five starless cores in MCLD 123.5+24.9. All of these cores appear to have strongly sub-virial masses, assuming, as we argue, that $^{13}$CO line data provide reliable estimates of velocity dispersions. Evidently, they are not self-gravitating, so it is no surprise that they are starless.

preprint2015arXiv

Strong variability of the coronal line region in NGC 5548

We present the second extensive study of the coronal line variability in an active galaxy. Our data set for the well-studied Seyfert galaxy NGC 5548 consists of five epochs of quasi-simultaneous optical and near-infrared spectroscopy spanning a period of about five years and three epochs of X-ray spectroscopy overlapping in time with it. Whereas the broad emission lines and hot dust emission varied only moderately, the coronal lines varied strongly. However, the observed high variability is mainly due to a flux decrease. Using the optical [FeVII] and X-ray OVII emission lines we estimate that the coronal line gas has a relatively low density of n~10^3/cm^3 and a relatively high ionisation parameter of log U~1. The resultant distance of the coronal line gas from the ionising source of about eight light years places this region well beyond the hot inner face of the dusty torus. These results imply that the coronal line region is an independent entity. We find again support for the X-ray heated wind scenario of Pier & Voit; the increased ionising radiation that heats the dusty torus also increases the cooling efficiency of the coronal line gas, most likely due to a stronger adiabatic expansion. The much stronger coronal line variability of NGC 5548 relative to that of NGC 4151 can also be explained within this picture. NGC 5548 has much stronger coronal lines relative to the low ionisation lines than NGC 4151 indicating a stronger wind, in which case a stronger adiabatic expansion of the gas and so fading of the line emission is expected.

preprint2015arXiv

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

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

preprint2015arXiv

Variability of the coronal line region in NGC 4151

We present the first extensive study of the coronal line variability in an active galaxy. Our data set for the nearby source NGC 4151 consists of six epochs of quasi-simultaneous optical and near-infrared spectroscopy spanning a period of about eight years and five epochs of X-ray spectroscopy overlapping in time with it. None of the coronal lines showed the variability behaviour observed for the broad emission lines and hot dust emission. In general, the coronal lines varied only weakly, if at all. Using the optical [Fe VII] and X-ray O VII emission lines we estimate that the coronal line gas has a relatively low density of n~10^3 cm^-3 and a relatively high ionisation parameter of log U~1. The resultant distance of the coronal line gas from the ionising source is about two light years, which puts this region well beyond the hot inner face of the obscuring dusty torus. The high ionisation parameter implies that the coronal line region is an independent entity rather than part of a continuous gas distribution connecting the broad and narrow emission line regions. We present tentative evidence for the X-ray heated wind scenario of Pier & Voit. We find that the increased ionising radiation that heats the dusty torus also increases the cooling efficiency of the coronal line gas, most likely due to a stronger adiabatic expansion.

preprint2014arXiv

Line-driven radiative outflows in luminous quasars

An analysis of ~19500 narrow (<200 km/s) CIV 1548.2,1550.8 absorbers in ~34000 Sloan Digital Sky Survey quasar spectra is presented. The statistics of the number of absorbers as a function of outflow-velocity shows that in approximately two-thirds of outflows, with multiple CIV absorbers present, absorbers are line-locked at the 500 km/s velocity separation of the CIV absorber doublet; appearing as 'triplets' in the quasar spectra. Line-locking is an observational signature of radiative line driving in outflowing material, where the successive shielding of 'clouds' of material in the outflow locks the clouds together in outflow velocity. Line-locked absorbers are seen in both broad absorption line quasars (BALs) and non-BAL quasars with comparable frequencies and with velocities out to at least 20000 km/s. There are no detectable differences in the absorber properties and the dust content of single CIV doublets and line-locked CIV doublets. The gas associated with both single and line-locked CIV absorption systems includes material with a wide range of ionization potential (14-138 eV). Both single and line-locked CIV absorber systems show strong systematic trends in their ionization as a function of outflow velocity, with ionization decreasing rapidly with increasing outflow velocity. Initial simulations, employing CLOUDY, demonstrate that a rich spectrum of line-locked signals at various velocities may be expected due to significant opacities from resonance lines of Li-, He- and H-like ions of O, C and N, along with contributions from HeII and HI resonance lines. The simulations confirm that line driving can be the dominant acceleration mechanism for clouds with N(HI) ~ 10^19 cm^-2.

preprint2013arXiv

Atomic data for S II - Toward Better Diagnostics of Chemical Evolution in High-redshift Galaxies

Absorption-line spectroscopy is a powerful tool used to estimate element abundances in the nearby as well as distant universe. The accuracy of the abundances thus derived is, naturally, limited by the accuracy of the atomic data assumed for the spectral lines. We have recently started a project to perform the new extensive atomic data calculations used for optical/UV spectral lines in the plasma modeling code Cloudy using state-of-the-art quantal calculations. Here we demonstrate our approach by focussing on S II, an ion used to estimate metallicities for Milky Way interstellar clouds as well as distant damped Lyman-alpha (DLA) and sub-DLA absorber galaxies detected in the spectra of quasars and gamma-ray bursts (GRBs). We report new extensive calculations of a large number of energy levels of S II, and the line strengths of the resulting radiative transitions. Our calculations are based on the configuration interaction approach within a numerical Hartree-Fock framework, and utilize both non-ralativistic and quasirelativistic one-electron radial orbitals. The results of these new atomic calculations are then incorporated into Cloudy and applied to a lab plasma, and a typical DLA, for illustrative purposes. The new results imply relatively modest changes (~0.04 dex) to the metallicities estimated from S II in past studies. These results will be readily applicable to other studies of S II in the Milky Way and other galaxies.

preprint2013arXiv

Classification and analysis of emission-line galaxies using mean field independent component analysis

We present an analysis of the optical spectra of narrow emission-line galaxies, based on mean field independent component analysis (MFICA). Samples of galaxies were drawn from the Sloan Digital Sky Survey (SDSS) and used to generate compact sets of `continuum' and `emission-line' component spectra. These components can be linearly combined to reconstruct the observed spectra of a wider sample of galaxies. Only 10 components - five continuum and five emission line - are required to produce accurate reconstructions of essentially all narrow emission-line galaxies; the median absolute deviations of the reconstructed emission-line fluxes, given the signal-to-noise ratio (S/N) of the observed spectra, are 1.2-1.8 sigma for the strong lines. After applying the MFICA components to a large sample of SDSS galaxies we identify the regions of parameter space that correspond to pure star formation and pure active galactic nucleus (AGN) emission-line spectra, and produce high S/N reconstructions of these spectra. The physical properties of the pure star formation and pure AGN spectra are investigated by means of a series of photoionization models, exploiting the faint emission lines that can be measured in the reconstructions. We are able to recreate the emission line strengths of the most extreme AGN case by assuming the central engine illuminates a large number of individual clouds with radial distance and density distributions, f(r) ~ r^gamma and g(n) ~ n^beta, respectively. The best fit is obtained with gamma = -0.75 and beta = -1.4. From the reconstructed star formation spectra we are able to estimate the starburst ages. These preliminary investigations serve to demonstrate the success of the MFICA-based technique in identifying distinct emission sources, and its potential as a tool for the detailed analysis of the physical properties of galaxies in large-scale surveys.

preprint2013arXiv

Interpreting the Ionization Sequence in AGN Emission-Line Spectra

We investigate the physical cause of the great range in the ionization level seen in the spectra of narrow lined active galactic nuclei (AGN). Mean field independent component analysis identifies examples of individual SDSS galaxies whose spectra are not dominated by emission due to star formation (SF), which we designate as AGN. We assembled high S/N ratio composite spectra of a sequence of these AGN defined by the ionization level of their narrow-line regions (NLR), extending down to very low-ionization cases. We used a local optimally emitting cloud (LOC) model to fit emission-line ratios in this AGN sequence. These included the weak lines that can be measured only in the co-added spectra, providing consistency checks on strong line diagnostics. After integrating over a wide range of radii and densities our models indicate that the radial extent of the NLR is the major parameter in determining the position of high to moderate ionization AGN along our sequence, providing a physical interpretation for their systematic variation. Higher ionization AGN contain optimally emitting clouds that are more concentrated towards the central continuum source than in lower ionization AGN. Our LOC models indicate that for the objects that lie on our AGN sequence, the ionizing luminosity is anticorrelated with the NLR ionization level, and hence anticorrelated with the radial concentration and physical extent of the NLR. A possible interpretation that deserves further exploration is that the ionization sequence might be an age sequence where low ionization objects are older and have systematically cleared out their central regions by radiation pressure. We consider that our AGN sequence instead represents a mixing curve of SF and AGN spectra, but argue that while many galaxies do have this type of composite spectra, our AGN sequence appears to be a special set of objects with negligible SF excitation.

preprint2012arXiv

Herschel observations of extended atomic gas in the core of the Perseus cluster

We present Herschel observations of the core of the Perseus cluster of galaxies. The brightest cluster galaxy, NGC 1275, is surrounded by a network of filaments previously imaged extensively in Hα and CO. In this work, we report detections of FIR lines with Herschel. All but one of the lines are spatially extended, with the [CII] line emission extending up to 25 kpc from the core. There is spatial and kinematical correlation among [CII], Hα and CO, which gives us confidence to model the different components of the gas with a common heating model. With the help of FIR continuum Herschel measurements, together with a suite of coeval radio, submm and infrared data, we performed a SED fitting of NGC 1275 using a model that contains contributions from dust emission as well as synchrotron AGN emission. The data indicate a low dust emissivity index, beta ~ 1, a total dust mass close to 10^7 solar mass, a cold dust component with temperature 38 \pm 2 K and a warm dust component with temperature of 116 \pm 9 K. The FIR-derived star formation rate (SFR) is 24 \pm 1 solar mass per yr, in close agreement with the FUV-derived SFR. We investigated in detail the source of the Herschel FIR and Hα emissions emerging from a core region 4 kpc in radius. Based on simulations conducted using the radiative transfer code, CLOUDY, a heating model comprising old and young stellar populations is sufficient to explain these observations. We have also detected [CII] in three well-studied regions of the filaments. We find a [OI]/[CII] ratio about 1 dex smaller than predicted by the otherwise functional Ferland (2009) model. The line ratio suggests that the lines are optically thick, as is typical of galactic PDRs, and implies that there is a large reservoir of cold atomic gas. [abridged]

preprint2012arXiv

Intervening BLR Clouds' Effects on Optical/UV Spectrum

Recent x-ray observations of Mrk 766 suggest that broad emission line region clouds cross our line of sight and produce variable x-ray absorption. Here we investigate what optical/ultraviolet spectroscopic features would be produced by such "Intervening BLR Clouds" (IBC) crossing our line of sight to the accretion disk, the source of the optical/UV continuum. Although the emission spectrum produced by intervening clouds is identical to the standard BLR model, they may produce absorption features on the optical or UV continuum. Single clouds will have little effect on the optical/UV spectrum because BLR clouds are likely to be much smaller than the accretion disk. This is unlike the X-ray case, where the radiation source is considerably smaller. However, an ensemble of intervening clouds will produce spectroscopic features in the FUV including a strong depression between the Lyman limit and Lyα. The amount of the depression will indicate the line-of-sight covering factor of clouds, an unknown quantity that is important for the ionization of the intergalactic medium and the energy budget of AGN. Comparison with observations suggests that the SED of Mrk 766 may be affected by intervening BLR clouds and IBC may exist in most of AGNs.

preprint2012arXiv

Two-Component Structure of the Hbeta Broad-Line Region in Quasars. I. Evidence from Spectral Principal Component Analysis

We report on a spectral principal component analysis (SPCA) of a sample of 816 quasars, selected to have small Fe II velocity shifts with spectral coverage in the rest wavelength range 3500--5500 Å. The sample is explicitly designed to mitigate spurious effects on SPCA induced by Fe II velocity shifts. We improve the algorithm of SPCA in the literature and introduce a new quantity, \emph{the fractional-contribution spectrum}, that effectively identifies the emission features encoded in each eigenspectrum. The first eigenspectrum clearly records the power-law continuum and very broad Balmer emission lines. Narrow emission lines dominate the second eigenspectrum. The third eigenspectrum represents the Fe II emission and a component of the Balmer lines with kinematically similar intermediate velocity widths. Correlations between the weights of the eigenspectra and parametric measurements of line strength and continuum slope confirm the above interpretation for the eigenspectra. Monte Carlo simulations demonstrate the validity of our method to recognize cross talk in SPCA and firmly rule out a single-component model for broad Hbeta. We also present the results of SPCA for four other samples that contain quasars in bins of larger Fe II velocity shift; similar eigenspectra are obtained. We propose that the Hbeta-emitting region has two kinematically distinct components: one with very large velocities whose strength correlates with the continuum shape, and another with more modest, intermediate velocities that is closely coupled to the gas that gives rise to Fe II emission.

preprint2011arXiv

Critical Differences and Clues in Eta Car's 2009 Event

We monitored Eta Carinae with HST WFPC2 and Gemini GMOS throughout the 2009 spectroscopic event, which was expected to differ from its predecessor in 2003 (Davidson et al. 2005). Here we report major observed differences between events, and their implications. Some of these results were quite unexpected. (1) The UV brightness minimum was much deeper in 2009. This suggests that physical conditions in the early stages of an event depend on different parameters than the "normal" inter-event wind. Extra mass ejection from the primary star is one possible cause. (2) The expected He II 4687 brightness maximum was followed several weeks later by another. We explain why this fact, and the timing of the 4687 maxima, strongly support a "shock breakup" hypothesis for X-ray and 4687 behavior as proposed 5-10 years ago. (3) We observed a polar view of the star via light reflected by dust in the Homunculus nebula. Surprisingly, at that location the variations of emission-line brightness and Doppler velocities closely resembled a direct view of the star; which should not have been true for any phenomena related to the orbit. This result casts very serious doubt on all the proposed velocity interpretations that depend on the secondary star's orbital motion. (4) Latitude-dependent variations of H I, He I and Fe II features reveal aspects of wind behavior during the event. In addition, we discuss implications of the observations for several crucial unsolved problems.

preprint2011arXiv

Ion-by-ion Cooling efficiencies

We present ion-by-ion cooling efficiencies for low-density gas. We use Cloudy (ver. 08.00) to estimate the cooling efficiencies for each ion of the first 30 elements (H-Zn) individually. We present results for gas temperatures between 1e4 and 1e8K, assuming low densities and optically thin conditions. When nonequilibrium ionization plays a significant role the ionization states deviate from those that obtain in collisional ionization equilibrium (CIE), and the local cooling efficiency at any given temperature depends on specific non-equilibrium ion fractions. The results presented here allow for an efficient estimate of the total cooling efficiency for any ionic composition. We also list the elemental cooling efficiencies assuming CIE conditions. These can be used to construct CIE cooling efficiencies for non-solar abundance ratios, or to estimate the cooling due to elements not explicitly included in any nonequilibrium computation. All the computational results are listed in convenient online tables.

preprint2011arXiv

N II 5668-5712, a New Class of Spectral Features in Eta Carinae

We report on the N II 5668-5712 emission and absorption lines in the spectrum of Eta Carinae. Spectral lines of the stellar wind regions can be classified into four physically distinct categories: 1) low-excitation emission such as H I and Fe II, 2) higher excitation He I features, 3) the N II lines discussed in this paper, and 4) He II emission. These categories have different combinations of radial velocity behavior, excitation processes, and dependences on the secondary star. The N II lines are the only known features that originate in "normal" undisturbed zones of the primary wind but depend primarily on the location of the hot secondary star. N II probably excludes some proposed models, such as those where He I lines originate in the secondary star's wind or in an accretion disk.

preprint2011arXiv

Powering The Intra-cluster Filaments in Cool-Core Clusters of Galaxies

The first radio surveys of the sky discovered that some large clusters of galaxies contained powerful sources of synchrotron emission. Optical images showed that long linear filaments with bizarre emission-line spectra permeated the intra-cluster medium. Recent observations in the infrared and radio show that these filaments have very strong emission lines of molecular hydrogen and carbon monoxide. The mass of molecular material is quite large, the gas is quite warm, and the filaments have not formed stars despite their ~Gyr age. I will discuss the general astrophysical context of large clusters of galaxies and how large masses of molecular gas can be heated to produce what we observe. The unique properties of the filaments are a result of the unique environment. Magnetically confined molecular filaments are surrounded by the hot intra-cluster medium. Thermal particles with keV energies enter atomic and molecular regions and produce a shower of secondary nonthermal electrons. These secondaries collisionally heat, excite, dissociate, and ionize the cool gas. While ionization is dominated by these secondary particles, recombination is controlled by charge exchange, which produces the unusual optical emission line spectrum. I will describe some of the physical processes that are unique to this environment and outline some of the atomic physics issues.

preprint2011arXiv

Understanding strong molecular hydrogen emission in astronomical environments

Here I describe recent studies of objects with molecular hydrogen emission that is strong relative to other spectral lines. Large telescopes and fast spectrometers have made the 2 μm window accessible even for relatively faint objects. I summarize several environments where strong molecular hydrogen 2.121 micron emission is observed. The line is hard to excite due to its large excitation potential, and is most emissive in regions that have temperatures that are nearly high enough to dissociate molecular hydrogen. I outline several case studies. In the Helix planetary nebula strong emission is produced by rapidly flowing molecular gas that is exposed to an intense ionizing radiation field. This advective production of molecular hydrogen is a fundamentally non-equilibrium process. In the filaments surrounding brightest cluster galaxies in cool core clusters ionizing particles penetrate into magnetically confined molecular cores and excite the gas. Finally, I outline ongoing work on the Crab Nebula, where the first complete maps of molecular emission have only recently been completed. Both ionizing particles and high-energy photons may be important. Finally I speculate on the origin of the correlation between H2 / H I intensity ratios and other properties found in Active Galaxies. This is suggestive of a hardening of the radiation field along the Eigenvector 1 sequence. In all of this work I take the approach of understanding H2 emission along with emission from low and moderate ionization species, a necessary step if we are to really understand the context in which molecular hydrogen emission forms.

preprint2010arXiv

A Bare Molecular Cloud at z~0.45

Several neutral species (MgI, SiI, CaI, FeI) have been detected in a weak MgII absorption line system (W_r(2796)~0.15 Angstroms) at z~0.45 along the sightline toward HE0001-2340. These observations require extreme physical conditions, as noted in D'Odorico (2007). We place further constraints on the properties of this system by running a wide grid of photoionization models, determining that the absorbing cloud that produces the neutral absorption is extremely dense (~100-1000/cm^3), cold (<100 K), and has significant molecular content (~72-94%). Structures of this size and temperature have been detected in Milky Way CO surveys, and have been predicted in hydrodynamic simulations of turbulent gas. In order to explain the observed line profiles in all neutral and singly ionized chemical transitions, the lines must suffer from unresolved saturation and/or the absorber must partially cover the broad emission line region of the background quasar. In addition to this highly unusual cloud, three other ordinary weak MgII clouds (within densities of ~0.005/cm^3 and temperatures of ~10000K) lie within 500 km/s along the same sightline. We suggest that the "bare molecular cloud", which appears to reside outside of a galaxy disk, may have had in situ star formation and may evolve into an ordinary weak MgII absorbing cloud.

preprint2010arXiv

Spitzer reveals what's behind Orion's Bar

We present Spitzer Space Telescope observations of 11 regions SE of the Bright Bar in the Orion Nebula, along a radial from the exciting star theta1OriC, extending from 2.6 to 12.1'. Our Cycle 5 programme obtained deep spectra with matching IRS short-high (SH) and long-high (LH) aperture grid patterns. Most previous IR missions observed only the inner few arcmin. Orion is the benchmark for studies of the ISM particularly for elemental abundances. Spitzer observations provide a unique perspective on the Ne and S abundances by virtue of observing the dominant ionization states of Ne (Ne+, Ne++) and S (S++, S3+) in Orion and H II regions in general. The Ne/H abundance ratio is especially well determined, with a value of (1.01+/-0.08)E-4. We obtained corresponding new ground-based spectra at CTIO. These optical data are used to estimate the electron temperature, electron density, optical extinction, and the S+/S++ ratio at each of our Spitzer positions. That permits an adjustment for the total gas-phase S abundance because no S+ line is observed by Spitzer. The gas-phase S/H abundance ratio is (7.68+/-0.30)E-6. The Ne/S abundance ratio may be determined even when the weaker hydrogen line, H(7-6) here, is not measured. The mean value, adjusted for the optical S+/S++ ratio, is Ne/S = 13.0+/-0.6. We derive the electron density versus distance from theta1OriC for [S III] and [S II]. Both distributions are for the most part decreasing with increasing distance. A dramatic find is the presence of high-ionization Ne++ all the way to the outer optical boundary ~12' from theta1OriC. This IR result is robust, whereas the optical evidence from observations of high-ionization species (e.g. O++) at the outer optical boundary suffers uncertainty because of scattering of emission from the much brighter inner Huygens Region.

preprint2010arXiv

Structure and Feedback in 30 Doradus I: Observations

We have completed a a new optical imaging and spectrophotometric survey of a 140 x 80 pc$^2$ region of 30 Doradus centered on R136, covering key optical diagnostic emission lines including \Ha, \Hb, \Hg, [O III] $λλ$4363, 4959, 5007, [N II] $λλ$6548, 6584, [S II] $λλ$6717, 6731 [S III] $λ$6312 and in some locations [S III] $λ$9069. We present maps of fluxes and intensity ratios for these lines, and catalogs of isolated ionizing stars, elephant-trunk pillars, and edge-on ionization fronts. The final science-quality spectroscopic data products are available to the public. Our analysis of the new data finds that, while stellar winds and supernovae undoubtedly produce shocks and are responsible for shaping the nebula, there are no global spectral signatures to indicate that shocks are currently an important source of ionization. We conclude that the considerable region covered by our survey is well described by photoionization from the central cluster where the ionizing continuum is dominated by the most massive O stars. We show that if 30 Dor were viewed at a cosmological distance, its integrated light would be dominated by its extensive regions of lower surface-brightness rather than by the bright, eye-catching arcs.

preprint2009arXiv

High-excitation emission lines near eta Carinae, and its likely companion star

In order to study the distribution of gas and ionizing radiation around eta Car and their implications for its likely companion star, we have examined high-excitation emission lines of [Ne III], [Fe III], etc., in spectra obtained with the HST/STIS instrument during 1998-2004. Our principal results, some of them unexpected, are: (1) The high-excitation fluxes varied systematically and non-trivially throughout eta Car's 5.5-year spectroscopic cycle. Instead of rising to a plateau after the 1998 "event," they changed continuously with a maximum in mid-cycle. (2) At one significant location a brief, strong secondary maximum occurred just before the 2003.5 spectroscopic event. (3) These emission lines are strongly concentrated at the "Weigelt knots" several hundred AU northwest of the star. With less certainty, [Ne III] appears to be somewhat more concentrated than [Fe III]. (4) A faster, blue-shifted component of each feature appears concentrated near the star and elongated perpendicular to the system's bipolar axis. This structure may be related to the equatorial outflow and/or to dense material known to exist along our line of sight to the star. (5) Using the photoionization program Cloudy, we estimated the range of parameters for the hot secondary star that would give satisfactory high-excitation line ratios in the ejecta. Teff = 39000 K and L = 400000 Lsun, for example, would be satisfactory. The allowed region in parameter space is wider (and mostly less luminous) than some previous authors suggested.

preprint2009arXiv

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

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

preprint1999arXiv

Emission-line Helium Abundances in Highly Obscured Nebulae

This paper outlines a way to determine the ICF using only infrared data. We identify four line pairs, [NeIII] 36$\micron$/[NeII] 12.8$\micron$, [NeIII]~15.6$\micron$ /[NeII] 12.8$\micron$, [ArIII] 9$\micron$/[ArII] 6.9$\micron$, and [ArIII] 21$\micron$/[ArII] 6.9$\micron$, that are sensitive to the He ICF. This happens because the ions cover a wide range of ionization, the line pairs are not sensitive to electron temperature, they have similar critical densities, and are formed within the He$^+$/H$^+$ region of the nebula. We compute a very wide range of photoionization models appropriate for galactic HII regions. The models cover a wide range of densities, ionization parameters, stellar temperatures, and use continua from four very different stellar atmospheres. The results show that each line pair has a critical intensity ratio above which the He ICF is always small. Below these values the ICF depends very strongly on details of the models for three of the ratios, and so other information would be needed to determine the helium abundance. The [Ar III] 9$\micron$/[ArII] 6.9$\micron$ ratio can indicate the ICF directly due to the near exact match in the critical densities of the two lines. Finally, continua predicted by the latest generation of stellar atmospheres are sufficiently hard that they routinely produce significantly negative ICFs.