Search

Cooling flows are common in galaxy clusters which have cool cores. The soft X-ray emission below 1 keV from the flows is mostly absorbed by cold dusty gas within the central cooling sites. Further evidence for this process is presented here through a more detailed analysis of the nearby Centaurus cluster and some additional clusters. Predictions of JWST near and mid-infrared spectra from cooling gas are presented. [NeVI] emission at 7.65 micron should be an important diagnostic of gas cooling between 6 and 1.5 times 10^5 K. The emerging overall picture of hidden cooling flows is explored. The efficiency of AGN feedback in reducing the total cooling rate in cool cores is shown to be above 50 percent for many clusters but is rarely above 90 per cent. The reduction is mostly in outer gas. Cooling dominates in elliptical galaxies and galaxy groups which have mass flow rates below about 15M/yr and in some massive clusters where rates can exceed 1000M/yr., Comment: 17 pages, 49 figures submitted to MNRAS

High resolution X-ray spectra reveal hidden cooling flows depositing cold gas at the centres of massive nearby early-type galaxies with little sign of normal star formation. Optical observations are revealing that a bottom-heavy Initial Mass Function is common within the inner kpc of similar galaxies. We revive the possibility that a low-mass star formation mode is operating due to the high thermal pressure in the cooling flow, thus explaining the accumulation of low-mass stars. We further explore whether such a mode operated in early, high-redshift galaxies and has sporadically continued to the present day. The idea links observed distant galaxies with black holes which are ultramassive for their stellar mass, nearby red nuggets and massive early-type galaxies. Nearby elliptical galaxies may be red but they are not dead., Comment: 4 pages, 2 Figures, Accepted for publication in MNRAS

Near IR spectroscopic reverberation of Active Galactic Nuclei (AGN) potentially allows the infrared (IR) broad line region (BLR) to be reverberated alongside the disc and dust continua, while the spectra can also reveal details of dust astro-chemistry. Here, we describe results of a short pilot study (17 near-IR spectra over a 183 d period) for Mrk 509. The spectra give a luminosity-weighted dust radius of $\langle R_{\mathrm{d,lum}} \rangle = 186 \pm 4$ light-days for blackbody (large grain dust), consistent with previous (photometric) reverberation campaigns, whereas carbon and silicate dust give much larger radii. We develop a method of calibrating spectral data in objects where the narrow lines are extended beyond the slit width. We demonstrate this by showing our resultant photometric band light curves are consistent with previous results, with a hot dust lag at >40 d in the K band, clearly different from the accretion disc response at <20 d in the z band. We place this limit of 40 d by demonstrating clearly that the modest variability that we do detect in the H and K band does not reverberate on time-scales of less than 40 d. We also extract the Pa$\beta$ line light curve, and find a lag which is consistent with the optical BLR H$\beta$ line of $\sim$70-90 d. This is important as direct imaging of the near-IR BLR is now possible in a few objects, so we need to understand its relation to the better studied optical BLR., Comment: 16 pages, 12 figures

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

Recently, we have uncovered Hidden Cooling Flows (HCF) in the X-ray spectra of the central Brightest Galaxies of 11 clusters, 1 group and 2 elliptical galaxies. Here we report such flows in a further 15 objects, consisting of 8 clusters, 3 groups, 3 ellipticals and 1 Red Nugget. The mass cooling rates are about 1 Msun/yr in the ellipticals, 2 to 20 Msun/yr in the groups and 20 to 100 Msun/yr in regular clusters. The Red Nugget, MRK1216, has an HCF of 10 Msun/yr. We review the fate of the cooled gas and investigate how some of it might accrete onto the central black hole. The gas is likely to be very cold and to have fragmented into low mass stars and smaller objects before being swallowed whole, with little luminous output. If such a scenario is correct and operates at a few Msun/yr then such objects may host the fastest growing black holes in the low redshift Universe. We briefly discuss the relevance of HCF to the growth of early galaxies and black holes., Comment: 15 pages, 18 Figures, accepted for publication in MNRAS

We have recently uncovered Hidden Cooling Flows (HCFs) in the XMM RGS spectra of 3 clusters of galaxies, Centaurus, Perseus and A1835. Here we search for them in a wider sample of objects: the X-ray brightest group NGC5044; 4 moderate X-ray luminosity clusters Sersic 159, A262, A2052 and RXJ0821; and 3 high X-ray luminosity clusters RXJ1532, MACS 1931 and the Phoenix cluster. Finally we examine two Virgo elliptical galaxies, M49 and M84. All statistically allow the addition of an HCF. We find a significant detection of an HCF in 6 clusters and 2 elliptical galaxies. The hidden mass cooling rates are 5 to 40 Solar masses per year for the normal clusters, 1000 Solar masses per year or more for the extreme clusters and 1 to 2 Solar masses per year for the elliptical galaxies. We discuss the implications of the results for the composition of the innermost parts of the massive host galaxies and look forward to future observations., Comment: 10 Pages, 13 figures, 3 tables, MNRAS submitted

Examination of emission lines in high-velocity resolution optical spectra of the Orion Nebula confirms that the velocity component on the red wing of the main ionization front emission line is due to backscattering in the Photon Dominated Region. This scattered light component has a weak wavelength dependence that is consistent with either general interstellar medium particles or particles in the foreground of the Orion Nebula Cluster. An anomalous line-broadening component that has been known for 60+ years is characterized in unprecedented detail. Although this extra broadening may be due to turbulence along the line-of-sight of our spectra, we explore the possibility that it is due to \alf\ waves in conditions where the ratio of magnetic and thermal energies are about equal and constant throughout the ionized gas.

Radiation pressure-driven outflows from luminous accreting supermassive black holes are an important part of active galactic nucleus (AGN) feedback. The effective Eddington limit, based on absorption of radiation by dust, not electron scattering, is revealed in the plane of AGN absorption column density $N_{\mathrm{H}}$ as a function of Eddington fraction $\lambda_{\mathrm{Edd}} = L_{\mathrm{bol}}/L_{\mathrm{Edd}}$, where a lack of objects is seen in the region where the effective limit is exceeded. Here, we conduct radiation simulation using the CLOUDY code to deduce the radiative force applied onto dusty gas at the nucleus and compare to the gravitational force to reveal the outflow region and its boundary with long-lived absorption clouds. We also investigate how the outflow condition is affected by various AGN and dust properties and distribution. As expected, the dust abundance has the largest effect on the $N_{\mathrm{H}} - \lambda_{\mathrm{Edd}}$ diagram since the higher the abundance, the more effective the radiative feedback, while the impact of the inner radius of the dusty gas shell, the shell width and the AGN spectral shape are relatively negligible. The presence of other central masses, such as a nuclear star cluster, can also make the feedback less effective. The AGN spectral energy distribution depends on the mass of the black hole and its spin. Though the effects of the AGN SED on the diagram are relatively small, the fraction of ionizing ultraviolet (UV) photons from the blackbody accretion disc is affected more by black hole mass than spin, and can influence the efficiency of radiation pressure., Comment: Accepted for publication in MNRAS; 11 pages, 7 figures

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, Comment: 11 pages. 7 figures, 2 tables MNRAS in press

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., Comment: 40 pages, 6 figures, 18 tables, Accepted for publication in ApJ

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., Comment: 16 pages, 6 figures, accepted to MNRAS

The partition function, $U$, the number of available states in an atom or molecules, is crucial for understanding the physical state of any astrophysical system in thermodynamic equilibrium. There are surprisingly few {\em useful} discussions of the partition function's numerical value. Textbooks often define $U$; some give tables of representative values, while others do a deep dive into the theory of a dense plasma. Most say that it depends on temperature, atomic structure, density, and that it diverges, that is, it goes to infinity, at high temperatures, but few give practical examples. We aim to rectify this. We show that there are two limits, one and two-electron (or closed-shell) systems like H or He, and species with a complicated electronic structure like C, N, O, and Fe. The high-temperature divergence does not occur for one and two-electron systems in practical situations since, at high temperatures, species are collisionally ionized to higher ionization stages and are not abundant. The partition function is then close to the statistical weight of the ground state. There is no such simplification for many-electron species. $U$ is temperature-sensitive across the range of temperatures where an ion is abundant but remains finite at even the highest practical temperatures. The actual value depends on highly uncertain truncation theories in high-density plasmas. We show that there are various theories for continuum lowering but that they are not in good agreement. This remains a long-standing unsolved problem., Comment: Tutorial introduction to the partition function in astrophysics

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., Comment: ApJ in press, comments welcome

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

The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, a.k.a BCGs (brightest cluster galaxies), is key to our understanding of active galactic nucleus feedback, and of the impact of environmental and local effects on star formation. We investigate the contribution of the thermal radiation from the cooling flow surrounding BCGs to the excitation of the filaments. We explore the effects of small levels of extra-heating (turbulence), and of metallicity, on the optical and infrared lines. Using the Cloudy code, we model the photoionization and photodissociation of a slab of gas of optical depth AV{\leq}30mag at constant pressure, in order to calculate self-consistently all of the gas phases, from ionized gas to molecular gas. The ionizing source is the EUV and soft X-ray radiation emitted by the cooling gas. We test these models comparing their predictions to the rich multi-wavelength observations, from optical to submillimeter. These models reproduce most of the multi-wavelength spectra observed in the nebulae surrounding the BCGs, not only the LINER-like optical diagnostics: [O iii]{\lambda} 5007 {\AA}/H\b{eta}, [N ii]{\lambda} 6583 {\AA}/H{\alpha} and ([S ii]{\lambda} 6716 {\AA}+[S ii]{\lambda} 6731 {\AA})/H{\alpha} but also the infrared emission lines from the atomic gas. The modeled ro-vib H2 lines also match observations, which indicates that near and mid-IR H2 lines are mostly excited by collisions between H2 molecules and secondary electrons produced naturally inside the cloud by the interaction between the X-rays and the cold gas in the filament. However, there is still some tension between ionized and molecular line tracers (i.e. CO), which requires to optimize the cloud structure and the density of the molecular zone., Comment: 16 pages, 9 figures, 3 tables. Accepted for publication in A&A

We have extended the membership and determined the 3-D structure of the large (0.19 pc) HH~269 sequence of shocks in the Orion Nebula. All of the components lie along a track that is highly tilted to the plane-of-the-sky and emerge from within the Orion-S embedded molecular cloud. Their source is probably either the highly obscured mm 9 source associated with a high N2H+ density core (more likely) or the more distant star COUP 632 (less likely). The former must be located in the Photon Dominated Region (PDR) underlying the ionized surface of the Orion South Cloud, while the latter would be embedded within the cloud. The flows seem to be episodic, with intervals of 1900 to 2600 years or 700 to 2600 years if COUP 632 is the source., Comment: 7 pages, 3 figures

The cosmic-ray ionization rate ($\zeta$, 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, $\lambda$ Cep, X Per, HD 73882, HD 154368, Cyg OB2 5, Cyg OB2 12. The value of $\zeta$ 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$ ($\zeta$(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 $\zeta$ is nearly 1 dex smaller for sightlines toward Cyg OB2 12. We estimate the average value of $\zeta$(H$_2$)= (95.69 $\pm$ 46.56) $\times$ 10$^{-16}$ s$^{-1}$ for models without PAHs., Comment: Accepted for publication in ApJ, 31 pages, 22 figures. A typo in the title is corrected

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., Comment: 21pages, 7 figures Accepted Astrophysical Journal

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., Comment: 9 pages, 2 figures

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{\alpha}. 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, Comment: accepted for publication in ApJ

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.

Dissociation of molecular hydrogen by secondary electrons produced by cosmic ray or X-ray ionization plays a crucial role in the chemistry of the densest part of molecular clouds. Here we study the effect of the mean kinetic energy of secondary electrons on this process. We compare predictions using a range of secondary electron energies and predictions of the cross-sections with the values in the UMIST database. We find that the predicted column densities change by nearly one dex., Comment: corrected minor typo

X-ray missions with microcalorimeter technology will resolve spectral features with unprecedented detail. In this work, we improve the H-like K$\alpha$ energies for elements between 6 $\leq Z \leq$ 30 for the release version of the spectral simulation code Cloudy to match laboratory energies. We update the ionization potential ($I_{\rm ion}$) for these elements and add a fourth-order polynomial to the level energy difference. This brings the release version of Cloudy into a near-perfect agreement with NIST. The updated energies are $\sim$ 15-4000 times more precise than that of the current release version of Cloudy (C17.02). These new changes will be a part of the next update to the release version, C17.03.

We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$\beta$, C IV, and Ly$\alpha$ broad emission lines. We find an extended disk-like H$\beta$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$\alpha$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$\alpha$ emission arising at smaller radii than the H$\beta$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$\beta$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole., Comment: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

The Soft X-ray Spectrometer (SXS) on board \textit{Hitomi}, with the unprecedented resolving power of R$\sim$1250, allowed the detection of members of the Fe XXV \ka\ complex emission spectra from the center of the Perseus Cluster. In this paper, we introduce a novel method of measuring the column density using the optically thin (Case A) to optically thick (Case B) transition for one- and two-electron systems. We compare the Fe XXV K$\alpha$ line ratios computed with CLOUDY with that from the \textit{Hitomi} observations in the outer region of the Perseus core using collision strengths from different atomic datasets, and obtain good agreement. We also show the effect of turbulence on Fe XXV K$\alpha$ line ratios and interplay between column density and metallicity. Besides, we discuss the atomic number dependence of transition probabilities for allowed and non-allowed transitions, which causes the highly charged He-like systems, such as Fe XXV, to behave fundamentally differently from He I.

We discuss the importance of Fe$^{23+}$ in determining the line intensities of the Fe XXV K$\alpha$ complex in an optically thick cloud, and investigate the prediction of Liedahl (2005) on Resonance Auger Destruction (RAD) with CLOUDY. Although initially motivated by the Perseus cluster, our calculations are extended to the wide range of column densities encountered in astronomy. A Fe XXV line photon can change/lose its identity upon absorption by three-electron iron as a result of "line interlocking". This may lead to the autoionization of the absorbing ion, ultimately destroying the Fe XXV K$\alpha$ photon by RAD. Out of the four members in the Fe XXV K$\alpha$ complex, a significant fraction of the x line photons is absorbed by Fe$^{23+}$ and destroyed, causing the x line intensity to decrease. For example, at a hydrogen column density of 10$^{25}$ cm$^{-2}$, $\sim$ 32\% of x photons are destroyed due to RAD while w is mostly unaffected. The line intensity of y is slightly ($\leq$2\%) reduced. z is not directly affected by RAD, but the contrasting behavior between z and x line intensities points towards the possible conversion of a tiny fraction ($\sim$ 2\%) of x photons into z photons. The change in line intensities due to Electron Scattering Escape (ESE) off fast thermal electrons is also discussed.

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., Comment: Accepted for publication in ApJ

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.

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., Comment: 19 pages, 9 figures, ApJ in press

The gamma-ray burst (GRB) afterglows provide an unique opportunity to study the interstellar medium (ISM) of star-forming galaxies at high-redshift. The GRB-DLAs (damped Lyman-$\alpha$ absorbers) contain a large neutral hydrogen column density, N(H I), and are observed against the GRB afterglow. A large fraction of GRB-DLAs show presence of molecular hydrogen (H$_2$) which is an indicator of star-formation. Hence it is important to study those GRB-DLAs which have H$_2$ lines to decipher and understand their physical conditions. The GRB-DLAs 121024A and 120815A, situated at redshift 2.30 and 2.36 respectively, are two such important H$_2$-bearing GRB-DLAs. Besides H$_2$, these two GRB-DLAs also show many metal lines. In this work we have carried out a detail numerical study on the H$_2$ lines, as well as on those metal lines, in GRB-DLAs 121024A and 120815A self-consistently. We use the spectral synthesis code CLOUDY for this study. This modeling helps us to determine the underlying physical conditions which give rise such lines and hence to understand these two GRB-DLAs in much more detail than any other previous investigation. We find that the hydrogen densities for these two H$_2$-bearing DLAs are $\geq$ 60 cm$^{-3}$. Moreover our study infers that the linear sizes are $\leq 17.7$ pc for these two GRB-DLAs, and the mean gas temperatures averaged over the cloud thickness, are $\leq$ 140 K., Comment: 10 pages, 5 figures, 6 tables. Accepted for publication in MNRAS

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., Comment: 12 pages, 5 figures. This paper was presented as an invited lecture at the 12th Serbian Conference on Spectral Line Shapes in Astrophysics in Vrdnik, Serbia

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.

The 180-day Space Telescope and Optical Reverberation Mapping campaign on NGC 5548 discovered an anomalous period, the broad-line region (BLR) holiday, in which the emission lines decorrelated from the continuum variations. This is important since the correlation between the continuum-flux variations and the emission-line response is the basic assumption for black hole (BH) mass determinations through reverberation mapping. During the BLR holiday, the high-ionization intrinsic absorption lines also decorrelated from the continuum as a result of variable covering factor of the line of sight (LOS) obscurer. The emission lines are not confined to the LOS, so this does not explain the BLR holiday. If the LOS obscurer is a disk wind, its streamlines must extend down to the plane of the disk and the base of the wind would lie between the BH and the BLR, forming an equatorial obscurer. This obscurer can be transparent to ionizing radiation, or can be translucent, blocking only parts of the SED, depending on its density. An emission-line holiday is produced if the wind density increases only slightly above its transparent state. Both obscurers are parts of the same wind, so they can have associated behavior in a way that explains both holidays. A very dense wind would block nearly all ionizing radiation, producing a Seyfert 2 and possibly providing a contributor to the changing-look AGN phenomenon. Disk winds are very common and we propose that the equatorial obscurers are too, but mostly in a transparent state., Comment: Accepted for publication in ApJL

We have recently initiated the first spectroscopic dust reverberation programme on active galactic nuclei (AGN) in the near-infrared. Spectroscopy enables measurement of dust properties, such as flux, temperature and covering factor, with higher precision than photometry. In particular, it enables measurement of both luminosity-based dust radii and dust response times. Here we report results from a one-year campaign on NGC 5548. The hot dust responds to changes in the irradiating flux with a lag time of ~70 light-days, similar to what was previously found in photometric reverberation campaigns. The mean and rms spectra are similar, implying that the same dust component dominates both the emission and the variations. The dust lag time is consistent with the luminosity-based dust radius only if we assume a wavelength-independent dust emissivity-law, i.e. a blackbody, which is appropriate for grains of large sizes (of a few microns). For such grains the dust temperature is ~1450 K. Therefore, silicate grains have most likely evaporated and carbon is the main chemical component. But the hot dust is not close to its sublimation temperature, contrary to popular belief. This is further supported by our observation of temperature variations largely consistent with a heating/cooling process. Therefore, the inner dust-free region is enlarged and the dusty torus rather a "dusty wall", whose inner radius is expected to be luminosity-invariant. The dust-destruction mechanism that enlarges the dust-free region seems to partly affect also the dusty region. We observe a cyclical decrease in dust mass with implied dust reformation times of ~5-6 months., Comment: 22 pages, 13 figures, accepted for publication in MNRAS

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

We have used existing optical emission and absorption lines, [C II] emission lines, and H I absorption lines to create a new model for a Central Column of material near the Trapezium region of the Orion Nebula. This was necessary because recent high spectral resolution spectra of optical emission lines and imaging spectra in the [C II] 158 micron line have shown that there are new velocity systems associated with the foreground Veil and the material lying between Theta 1 Ori C and the Main Ionization Front of the nebula. When a family of models generated with the spectral synthesis code Cloudy were compared with the surface brightness of the emission lines and strengths of the Veil absorption lines seen in the Trapezium stars, distances from Theta 1 Ori C, were derived, with the closest, highest ionization layer being 1.3 pc. The line of sight distance of this layer is comparable with the size of the inner Huygens Region in the plane of the sky. These layers are all blueshifted with respect to the Orion Nebula Cluster of stars, probably because of the pressure of a hot central bubble created by Theta 1 Ori C's stellar wind. We find velocity components that are ascribed to both sides of this bubble. Our analysis shows that the foreground [C II] 158 micron emission is part of a previously identified layer that forms a portion of a recently discovered expanding shell of material covering most of the larger Extended Orion Nebula., Comment: 13 pages, 9 figures, accepted for publication in the Astrophysical Journal

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)

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

We present improved fits to our treatment of suppression of dielectronic recombination at intermediate densities. At low densities, most recombined excited states eventually decay to the ground state, and therefore the total dielectronic recombination rate to all levels is preserved. At intermediate densities, on the other hand, collisions can lead to ionization of higher-lying excited states, thereby suppressing the dielectronic recombination rate. The improved suppression factors presented here, although highly approximate, allow summed recombination rate coefficients to be used to intermediate densities. There have been several technical improvements to our previously presented fits. For H- through B-like ions the activation log densities have been adjusted to better reproduce existing data. For B-, C-, Al-, and Si-like ions secondary autoionization is now included. The treatment of density discontinuity in electron excitations out of ground state H-, He-, and Ne-like ions has been improved. These refined dielectronic recombination suppression factors are used in the most recent version of the plasma simulation code Cloudy. We show how the ionization and emission spectrum change when this physics is included. Although these suppression factors improve the treatment of intermediate densities, they are highly approximate and are not a substitution for a complete collisional-radiative model of the ionization balance.

In this paper we report the experimental implementation of a theoretically-proposed technique for creating a photoionized plasma in the laboratory using X-ray line radiation. Using a Sn laser-plasma to irradiate an Ar gas target, the photoionization parameter, \xi = 4{\pi}F/Ne, reached values of order 50 erg cm/s, where F is the radiation flux in erg/cm2/s. The significance of this is that this technique allows us to mimic effective spectral radiation temperatures in excess of 1 keV. We show that our plasma starts to be collisionally dominated before the peak of the X-ray drive. However, the technique is extendable to higher energy laser systems to create plasmas with parameters relevant to benchmarking codes used to model astrophysical objects., Comment: 10 pages; 8 figures; Physical Review E, in press

The requirement of intermediate line component in the recently observed spectra of several AGNs points to possibility of the existence of a physically separate region between broad line region (BLR) and narrow line region (NLR). In this paper we explore the emission from intermediate line region (ILR) by using the photoionization simulations of the gas clouds distributed radially from the AGN center. The gas clouds span distances typical for BLR, ILR and NLR, and the appearance of dust at the sublimation radius is fully taken into account in our model. Single cloud structure is calculated under the assumption of the constant pressure. We show that the slope of the power law cloud density radial profile does not affect the existence of ILR in major types of AGN. We found that the low ionization iron line, Fe~II, appears to be highly sensitive for the presence of dust and therefore becomes potential tracer of dust content in line emitting regions. We show that the use of disk-like cloud density profile computed at the upper part of the accretion disc atmosphere reproduces the observed properties of the line emissivities. In particular, the distance of H${\beta}$ line inferred from our model agrees with that obtained from the reverberation mapping studies in Sy1 galaxy NGC 5548., Comment: 15 pages, 13 figures

In this paper we explore the intermediate line region (ILR) by using the photoionisation simulations of the gas clouds present at different radial distances from the center, corresponding to the locations from BLR out to NLR in four types of AGNs. We let for the presence of dust whenever conditions allow for dust existence. All spectral shapes are taken from the recent multi-wavelength campaigns. The cloud density decreases with distance as a power law. We found that the slope of the power law density profile does not affect the line emissivity radial profiles of major emission lines: H${\beta}$, He~II, Mg~II, C~III] ~and [O~III]. When the density of the cloud at the sublimation radius is as high as 10$^{11.5}$ cm$^{-3}$, the ILR should clearly be seen in the observations independently of the shape of the illuminating radiation. Moreover, our result is valid for low ionization nuclear emission regions of active galaxies., Comment: 8 pages, 2 figures, Accepted for publication in the Journal Frontiers in Astronomy and Space Sciences

The gas kinetic temperature ($T_K$) of various interstellar environments is often inferred from observations that can deduce level populations of atoms, ions, or molecules using spectral line observations; H I 21 cm is perhaps the most widely used with a long history. Usually the H I 21 cm line is assumed to be in thermal equilibrium and the populations are given by the Boltzmann distribution. A variety of processes, many involving Lyman alpha ($Ly\alpha$), can affect the 21 cm line. Here we show how this is treated in the spectral simulation code Cloudy, and present numerical simulations of environments where this temperature indicator is used, with a detailed treatment of the physical processes that determine level populations within $H^0$. We discuss situations where this temperature indicator traces $T_K$, cases where they fail, as well as the effects of $Ly\alpha$ pumping on the 21 cm spin temperature. We also show that the $Ly\alpha$ excitation temperature rarely traces the gas kinetic temperature., Comment: Accepted for publication in ApJ

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

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

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly., Comment: ApJ in press. Replaced with the accepted version

The common assumption that Theta-1-Ori C is the dominant ionizing source for the Orion Nebula is critically examined. This assumption underlies much of the existing analysis of the nebula. In this paper we establish through comparison of the relative strengths of emission lines with expectations from Cloudy models and through the direction of the bright edges of proplyds that Theta-2-Ori-A, which lies beyond the Bright Bar, also plays an important role. Theta-1-Ori-C does dominate ionization in the inner part of the Orion Nebula, but outside of the Bright Bar as far as the southeast boundary of the Extended Orion Nebula, Theta-2-Ori-A is the dominant source. In addition to identifying the ionizing star in sample regions, we were able to locate those portions of the nebula in 3-D. This analysis illustrates the power of MUSE spectral imaging observations in identifying sources of ionization in extended regions., Comment: 7 pages with 5 figures

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548., Comment: 22 pages, 13 figures, accepted to ApJ

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