Search

The flux ratios of high-ionization lines are commonly assumed to indicate the metallicity of the broad emission line region in luminous quasars. When accounting for the variation in their kinematic profiles, we show that the NV/CIV, (SiIV+OIV])/CIV and NV/Lya line ratios do not vary as a function of the quasar continuum luminosity, black hole mass, or accretion rate. Using photoionization models from CLOUDY , we further show that the observed changes in these line ratios can be explained by emission from gas with solar abundances, if the physical conditions of the emitting gas are allowed to vary over a broad range of densities and ionizing fluxes. The diversity of broad line emission in quasar spectra can be explained by a model with emission from two kinematically distinct regions, where the line ratios suggest that these regions have either very different metallicity or density. Both simplicity and current galaxy evolution models suggest that near-solar abundances, with parts of the spectrum forming in high-density clouds, are more likely. Within this paradigm, objects with stronger outflow signatures show stronger emission from gas which is denser and located closer to the ionizing source, at radii consistent with simulations of line-driven disc-winds. Studies using broad-line ratios to infer chemical enrichment histories should consider changes in density and ionizing flux before estimating metallicities., 12 pages, 7 figures, accepted for publication in MNRAS. The interactive version of Figs. 3 and 4 will be made available as supplementary online material with the journal, and will be available at https://people.ast.cam.ac.uk/~alr53/CIV_yvar/all_plots.html until August 2021

We present results of time-series analysis of the first year of the Fairall 9 intensive disc-reverberation campaign. We used Swift and the Las Cumbres Observatory global telescope network to continuously monitor Fairall 9 from X-rays to near-infrared at a daily to subdaily cadence. The cross-correlation function between bands provides evidence for a lag spectrum consistent with the τ ∝ λ4/3 scaling expected for an optically thick, geometrically thin blackbody accretion disc. Decomposing the flux into constant and variable components, the variable component’s spectral energy distribution is slightly steeper than the standard accretion disc prediction. We find evidence at the Balmer edge in both the lag and flux spectra for an additional bound-free continuum contribution that may arise from reprocessing in the broad-line region. The inferred driving light curve suggests two distinct components, a rapidly variable (100 d) component with an opposite lag to the reverberation signal.

At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact parameter calculations or semi-empirical fits and the rate coefficients agree within a factor of two. We show in this paper that these uncertainties cause errors lower than 5% in the emission of radio recombination lines (RRL) of most ionized plasmas of typical nebulae. However, in special circumstances where the transitions between Rydberg levels are amplified by maser effects, the errors can increase up to 20%. We present simulations of the optical depth and H$n\alpha$ line emission of Active Galactic Nuclei (AGN) Broad Line Regions (BLRs) and the Orion Nebula Blister to showcase our findings., Comment: 16 pages (accepted by MNRAS, replaced)

Absorption and emission lines of the iron-peak species Fe II are prominent in the infrared, optical and ultraviolet spectra of a myriad of astrophysical sources, requiring extensive and highly reliable sets of atomic structure and collisional data for an accurate quantitative analysis. However, comparisons among existing calculations reveal large discrepancies in the effective collision strengths, often up to factors of three, highlighting the need for further steps towards new converged calculations. Here we report a new 20 configuration, 6069 level atomic structure model, calculated using the multi-configurational Dirac-Fock method. Collision strengths and effective collision strengths are presented, for a wide range of temperatures of astrophysical relevance, from substantial 262 level and 716 level Dirac R-matrix calculations, plus a 716 level Breit-Pauli R-matrix calculation. Convergence of the scattering calculations is discussed, and results are critically compared with existing data in the literature, providing us with error estimates for our data. As a consequence, we assign an uncertainty of ±15% to relevant forbidden and allowed transitions encompassed within a 50 level subset of the 716 level Dirac R-matrix dataset. To illustrate the implications of our new datasets for the analysis of astronomical observations of Fe II, they are incorporated into the Cloudy modelling code, sample Fe II spectra are generated and compared.

Based on prior precision photometry and cluster age analysis, the bright star GSC 00154−01819 is a possible young pre-main sequence member of the Rosette cluster, NGC 2244. As part of a comprehensive study of the large-scale structure of the Rosette and its excitation by the cluster stars, we noted this star as a potential backlight for a probe of the interstellar medium and extinction along the sight line towards a distinctive nebular feature projected on to the cluster centre. New high-resolution spectra of the star were taken with the University College London Echelle Spectrograph of the AAT. They reveal that rather than being a reddened spectral type B or A star within the Mon OB2 association, it is a nearby, largely unreddened, solar twin of spectral type G2V less than 180 Myr old. It is about 219 pc from the Sun with a barycentric radial velocity of +14.35 ± 1.99 km s−1. The spectrum of the Rosette behind it and along this line of sight shows a barycentric radial velocity of +26.0 ± 2.4 km s−1 in H α, and a full width at half-maximum velocity dispersion of 61.94 ± 1.38 km s−1.

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., 9 pages, $ figures, 2 tables, accepted for publication in MNRAS

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., 41 pages, 17 figures, 9 tables, accepted to MNRAS

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., 10 pages, 6 figures; accepted for publication in MNRAS. arXiv admin note: substantial text overlap with arXiv:1501.02928

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

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., 10 pages, 5 tables, 3 figures. Fixed typo in Eq. 14

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, Published by MNRAS, online tables can be found at http://mnras.oxfordjournals.org/content/suppl/2014/04/28/stu514.DC1

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., Comment: MNRAS accepted. 29 pages, 24 figures, 3 tables

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

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