Showing 1–50 of 208 results for author: Ferland, G J

Unveiling the Cosmic Chemistry: Revisiting the Mass-Metallicity Relation with JWST/NIRSpec at 4 < z < 10

Authors: Arnab Sarkar, Priyanka Chakraborty, Mark Vogelsberger, Michael McDonald, Paul Torrey, Alex M. Garcia, Gourav Khullar, Gary J. Ferland, William Forman, Scott Wolk, Benjamin Schneider, Mark Bautz, Eric Miller, Catherine Grant, John ZuHone

Abstract: We present star formation rates (SFR), the mass-metallicity relation (MZR), and the SFR-dependent MZR across redshifts 4 to 10 using 81 star-forming galaxies observed by the JWST NIRSpec employing both low-resolution PRISM and medium-resolution gratings, including galaxies from the JADES GOODS-N and GOODS-S fields, the JWST-PRIMAL Legacy Survey, and additional galaxies from the literature in Abell… ▽ More We present star formation rates (SFR), the mass-metallicity relation (MZR), and the SFR-dependent MZR across redshifts 4 to 10 using 81 star-forming galaxies observed by the JWST NIRSpec employing both low-resolution PRISM and medium-resolution gratings, including galaxies from the JADES GOODS-N and GOODS-S fields, the JWST-PRIMAL Legacy Survey, and additional galaxies from the literature in Abell 2744, SMACS-0723, RXJ2129, BDF, COSMOS, and MACS1149 fields. These galaxies span a 3 dex stellar mass range of $10^7 < M_{\ast}/M_{\odot} < 10^{10}$, with an average SFR of $7.2 \pm 1.2 M_{\odot} {\rm yr}^{-1}$ and an average metallicity of $12+{\rm log(O/H)} = 7.91 \pm 0.08$. Our findings align with previous observations up to $z=8$ for the MZR and indicate no deviation from local universe FMR up to this redshift. Beyond $z=8$, we observe a significant deviation $\sim 0.27$ dex) in FMR, consistent with recent JWST findings. We also integrate CEERS (135 galaxies) and JADES (47 galaxies) samples with our data to study metallicity evolution with redshift in a combined sample of 263 galaxies, revealing a decreasing metallicity trend with a slope of $0.067 \pm 0.013$, consistent with IllustrisTNG and EAGLE, but contradicts with FIRE simulations. We introduce an empirical mass-metallicity-redshift (MZ-$z$ relation): $12+{\rm log(O/H)}=6.29 + 0.237 \times{\rm log}(M_{\ast}/M_{\odot}) - 0.06 \times (1+z)$, which accurately reproduces the observed trends in metallicity with both redshift and stellar mass. This trend underscores the ``Grand Challenge'' in understanding the factors driving high-redshift galactic metallicity trends, such as inflow, outflow, and AGN/stellar feedback -- and emphasizes the need for further investigations with larger samples and enhanced simulations. △ Less

Submitted 16 August, 2024; v1 submitted 15 August, 2024; originally announced August 2024.

Comments: 24 Pages, 9 Figures. Submitted to ApJ. Comments are welcome!

Recent Update of Gas-Phase Chemical Reactions and Molecular Lines of TiO in Cloudy

Authors: Gargi Shaw, Gary J. Ferland, Phillip Stancil, Ryan Porter

Abstract: We present our current update on the gas-phase chemical reactions and spectral lines of TiO in the spectral synthesis code Cloudy. For this purpose, we have added 229 Ti-related reactions in the chemical network. In addition, we consider 230 fine-structure energy levels, the corresponding 223 radiative transitions, and 444 collisional transitions with ortho and para H$_2$ and predict 66 TiO lines.… ▽ More We present our current update on the gas-phase chemical reactions and spectral lines of TiO in the spectral synthesis code Cloudy. For this purpose, we have added 229 Ti-related reactions in the chemical network. In addition, we consider 230 fine-structure energy levels, the corresponding 223 radiative transitions, and 444 collisional transitions with ortho and para H$_2$ and predict 66 TiO lines. We perform spectroscopic simulations of TiO emission from the circumstellar region of the oxygen-rich red supergiant VY Canis Majoris to validate our update. Our model reproduces the observed TiO column density. This update is helpful in modeling dust-free astrophysical environments where Ti is in the gas phase and TiO can form. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: 13 pages, 4 figures, accepted for publication in RMxAA

AGN STORM 2: VIII. Investigating the Narrow Absorption Lines in Mrk 817 Using HST-COS Observations

Authors: Maryam Dehghanian, Nahum Arav, Gerard A. Kriss, Missagh Mehdipour, Doyee Byun, Gwen Walker, Mayank Sharma, Aaron J. Barth, Misty C. Bentz, Benjamin D. Boizelle, Michael S. Brotherton, Edward M. Cackett, Elena Dalla Bonta, Gisella De Rosa, Gary J. Ferland, Carina Fian, Alexei V. Filippenko, Jonathan Gelbord, Michael R. Goad, Keith Horne, Yasaman Homayouni, Dragana Ilic, Michael D. Joner, Erin A. Kara, Shai Kaspi , et al. (17 additional authors not shown)

Abstract: We observed the Seyfert 1 galaxy Mrk817 during an intensive multi-wavelength reverberation mapping campaign for 16 months. Here, we examine the behavior of narrow UV absorption lines seen in HST/COS spectra, both during the campaign and in other epochs extending over 14 years. We conclude that while the narrow absorption outflow system (at -3750 km/s with FWHM=177 km/s) responds to the variations… ▽ More We observed the Seyfert 1 galaxy Mrk817 during an intensive multi-wavelength reverberation mapping campaign for 16 months. Here, we examine the behavior of narrow UV absorption lines seen in HST/COS spectra, both during the campaign and in other epochs extending over 14 years. We conclude that while the narrow absorption outflow system (at -3750 km/s with FWHM=177 km/s) responds to the variations of the UV continuum as modified by the X-ray obscurer, its total column density (logNH =19.5 cm-2) did not change across all epochs. The adjusted ionization parameter (scaled with respect to the variations in the Hydrogen ionizing continuum flux) is log UH =-1.0. The outflow is located at a distance smaller than 38 parsecs from the central source, which implies a hydrogen density of nH > 3000 cm-3. The absorption outflow system only covers the continuum emission source and not the broad emission line region, which suggests that its transverse size is small (< 1e16 cm), with potential cloud geometries ranging from spherical to elongated along the line of sight. △ Less

Submitted 8 July, 2024; v1 submitted 4 July, 2024; originally announced July 2024.

Comments: 18 pages, 12 Figures, Accepted for publication in ApJ

AGN STORM 2: IX. Studying the Dynamics of the Ionized Obscurer in Mrk 817 with High-resolution X-ray Spectroscopy

Authors: Fatima Zaidouni, Erin Kara, Peter Kosec, Missagh Mehdipour, Daniele Rogantini, Gerard A. Kriss, Ehud Behar, Jelle Kaastra, Aaron J. Barth, Edward M. Cackett, Gisella De Rosa, Yasaman Homayouni, Keith Horne, Hermine Landt, Nahum Arav, Misty C. Bentz, Michael S. Brotherton, Elena Dalla Bontà, Maryam Dehghanian, Gary J. Ferland, Carina Fian, Jonathan Gelbord, Michael R. Goad, Diego H. González Buitrago, Catherine J. Grier , et al. (23 additional authors not shown)

Abstract: We present the results of the XMM-Newton and NuSTAR observations taken as part of the ongoing, intensive multi-wavelength monitoring program of the Seyfert 1 galaxy Mrk 817 by the AGN Space Telescope and Optical Reverberation Mapping 2 (AGN STORM 2) Project. The campaign revealed an unexpected and transient obscuring outflow, never before seen in this source. Of our four XMM-Newton/NuSTAR epochs,… ▽ More We present the results of the XMM-Newton and NuSTAR observations taken as part of the ongoing, intensive multi-wavelength monitoring program of the Seyfert 1 galaxy Mrk 817 by the AGN Space Telescope and Optical Reverberation Mapping 2 (AGN STORM 2) Project. The campaign revealed an unexpected and transient obscuring outflow, never before seen in this source. Of our four XMM-Newton/NuSTAR epochs, one fortuitously taken during a bright X-ray state has strong narrow absorption lines in the high-resolution grating spectra. From these absorption features, we determine that the obscurer is in fact a multi-phase ionized wind with an outflow velocity of $\sim$5200 km s$^{-1}$, and for the first time find evidence for a lower ionization component with the same velocity observed in absorption features in the contemporaneous HST spectra. This indicates that the UV absorption troughs may be due to dense clumps embedded in diffuse, higher ionization gas responsible for the X-ray absorption lines of the same velocity. We observe variability in the shape of the absorption lines on timescales of hours, placing the variable component at roughly 1000 $R_g$ if attributed to transverse motion along the line of sight. This estimate aligns with independent UV measurements of the distance to the obscurer suggesting an accretion disk wind at the inner broad line region. We estimate that it takes roughly 200 days for the outflow to travel from the disk to our line of sight, consistent with the timescale of the outflow's column density variations throughout the campaign. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Consequences of a low-mass, high-pressure, star formation mode in early galaxies

Authors: A. C. Fabian, J. S. Sanders, G. J. Ferland, B. R. McNamara, C. Pinto, S. A. Walker

Abstract: 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 therma… ▽ More 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. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 4 pages, 2 Figures, Accepted for publication in MNRAS

The first spectroscopic IR reverberation programme on Mrk 509

Authors: J. A. J. Mitchell, M. J. Ward, D. Kynoch, J. V. Hernández Santisteban, K. Horne, J. -U. Pott, J. Esser, P. Mercatoris, C. Packham, G. J. Ferland, A. Lawrence, T. Fischer, A. J. Barth, C. Villforth, H. Winkler

Abstract: 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… ▽ More 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$β$ line light curve, and find a lag which is consistent with the optical BLR H$β$ 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. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 16 pages, 12 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 529, Issue 4, April 2024, Pages 4824 - 4839

The 23.01 release of Cloudy

Authors: Chamani M. Gunasekera, Peter A. M. van Hoof, Marios Chatzikos, Gary J. Ferland

Abstract: We announce the C23.01 update of Cloudy. This corrects a simple coding error, present since $\sim$ 1990, in one routine that required a conversion from the line-center to the mean normalization of the Ly$α$ optical depth. This affects the destruction of H I Ly$α$ by background opacities. Its largest effect is upon the Ly$α$ intensity in high-ionization dusty clouds, where the predicted intensity i… ▽ More We announce the C23.01 update of Cloudy. This corrects a simple coding error, present since $\sim$ 1990, in one routine that required a conversion from the line-center to the mean normalization of the Ly$α$ optical depth. This affects the destruction of H I Ly$α$ by background opacities. Its largest effect is upon the Ly$α$ intensity in high-ionization dusty clouds, where the predicted intensity is now up to three times stronger. Other properties that depend on Ly$α$ destruction, such as grain infrared emission, change in response. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 4 pages, 1 figure

AGN STORM 2. VI. Mapping Temperature Fluctuations in the Accretion Disk of Mrk 817

Authors: Jack M. M. Neustadt, Christopher S. Kochanek, John Montano, Jonathan Gelbord, Aaron J. Barth, Gisella De Rosa, Gerard A. Kriss, Edward M. Cackett, Keith Horne, Erin A. Kara, Hermine Landt, Hagai Netzer, Nahum Arav, Misty C. Bentz, Elena Dalla Bonta, Maryam Dehghanian, Pu Du, Rick Edelson, Gary J. Ferland, Carina Fian, Travis Fischer, Michael R. Goad, Diego H. Gonzalez Buitrago, Varoujan Gorjian, Catherine J. Grier , et al. (27 additional authors not shown)

Abstract: We fit the UV/optical lightcurves of the Seyfert 1 galaxy Mrk 817 to produce maps of the accretion disk temperature fluctuations $δT$ resolved in time and radius. The $δT$ maps are dominated by coherent radial structures that move slowly ($v \ll c$) inwards and outwards, which conflicts with the idea that disk variability is driven only by reverberation. Instead, these slow-moving temperature fluc… ▽ More We fit the UV/optical lightcurves of the Seyfert 1 galaxy Mrk 817 to produce maps of the accretion disk temperature fluctuations $δT$ resolved in time and radius. The $δT$ maps are dominated by coherent radial structures that move slowly ($v \ll c$) inwards and outwards, which conflicts with the idea that disk variability is driven only by reverberation. Instead, these slow-moving temperature fluctuations are likely due to variability intrinsic to the disk. We test how modifying the input lightcurves by smoothing and subtracting them changes the resulting $δT$ maps and find that most of the temperature fluctuations exist over relatively long timescales ($\sim$100s of days). We show how detrending AGN lightcurves can be used to separate the flux variations driven by the slow-moving temperature fluctuations from those driven by reverberation. We also simulate contamination of the continuum emission from the disk by continuum emission from the broad line region (BLR), which is expected to have spectral features localized in wavelength, such as the Balmer break contaminating the $U$ band. We find that a disk with a smooth temperature profile cannot produce a signal localized in wavelength and that any BLR contamination should appear as residuals in our model lightcurves. Given the observed residuals, we estimate that only $\sim$20% of the variable flux in the $U$ and $u$ lightcurves can be due to BLR contamination. Finally, we discus how these maps not only describe the data, but can make predictions about other aspects of AGN variability. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 23 pages, 18 figures, submitting to ApJ, comments welcome

Generation of photoionized plasmas in the laboratory of relevance to accretion-powered x-ray sources using keV line radiation

Authors: D. Riley, R. L. Singh, S White, M. Charlwood, D. Bailie, C. Hyland, T. Audet, G. Sarri, B. Kettle, G. Gribakin, S. J. Rose, E. G. Hill, G. J. Ferland, R. J. R. Williams, F. P. Keenan

Abstract: 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… ▽ More 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 ξ > 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. △ Less

Submitted 19 March, 2024; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: High Energy Density Physics, in press; 21 pages, 10 figures

The 2023 release of Cloudy

Authors: Marios Chatzikos, Stefano Bianchi, Francesco Camilloni, Priyanka Chakraborty, Chamani M. Gunasekera, Francisco Guzmán, Jonathan S. Milby, Arnab Sarkar, Gargi Shaw, Peter A. M. van Hoof, Gary J. Ferland

Abstract: We describe the 2023 release of the spectral synthesis code Cloudy. Since the previous major release, migrations of our online services motivated us to adopt git as our version control system. This change alone led us to adopt an annual release scheme, accompanied by a short release paper, the present being the inaugural. Significant changes to our atomic and molecular data have improved the accur… ▽ More We describe the 2023 release of the spectral synthesis code Cloudy. Since the previous major release, migrations of our online services motivated us to adopt git as our version control system. This change alone led us to adopt an annual release scheme, accompanied by a short release paper, the present being the inaugural. Significant changes to our atomic and molecular data have improved the accuracy of Cloudy predictions: we have upgraded our instance of the Chianti database from version 7 to 10; our H- and He-like collisional rates to improved theoretical values; our molecular data to the most recent LAMDA database, and several chemical reaction rates to their most recent UDfA and KiDA values. Finally, we describe our progress on upgrading Cloudy's capabilities to meet the requirements of the X-ray microcalorimeters aboard the upcoming XRISM and Athena missions, and outline future development that will make Cloudy of use to the X-ray community. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: 18 pages, 10 figures, accepted for publication in RMxAA

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

Authors: Y. Homayouni, Gerard A. Kriss, Gisella De Rosa, Rachel Plesha, Edward M. Cackett, Michael R. Goad, Kirk T. Korista, Keith Horne, Travis Fischer, Tim Waters, Aaron J. Barth, Erin A. Kara, Hermine Landt, Nahum Arav, Benjamin D. Boizelle, Misty C. Bentz, Michael S. Brotherton, Doron Chelouche, Elena Dalla Bonta, Maryam Dehghanian, Pu Du, Gary J. Ferland, Carina Fian, Jonathan Gelbord, Catherine J. Grier , et al. (27 additional authors not shown)

Abstract: 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 tim… ▽ More 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. △ Less

Submitted 5 January, 2024; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: 20 pages, 8 figures, Accepted for publication in ApJ

AGN STORM 2. IV. Swift X-ray and ultraviolet/optical monitoring of Mrk 817

Authors: Edward M. Cackett, Jonathan Gelbord, Aaron J. Barth, Gisella De Rosa, Rick Edelson, Michael R. Goad, Yasaman Homayouni, Keith Horne, Erin A. Kara, Gerard A. Kriss, Kirk T. Korista, Hermine Landt, Rachel Plesha, Nahum Arav, Misty C. Bentz, Benjamin D. Boizelle, Elena Dalla Bonta, Maryam Dehghanian, Fergus Donnan, Pu Du, Gary J. Ferland, Carina Fian, Alexei V. Filippenko, Diego H. Gonzalez Buitrago, Catherine J. Grier , et al. (26 additional authors not shown)

Abstract: The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitori… ▽ More The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitoring shows that Mrk 817 was in a significantly fainter state than in previous observations, with only a brief flare where it reached prior flux levels. The X-ray spectrum is heavily obscured. The UV/optical light curves show significant variability throughout the campaign and are well correlated with one another, but uncorrelated with the X-rays. Combining the Swift UV/optical light curves with Hubble UV continuum light curves, we measure interband continuum lags, $τ(λ)$, that increase with increasing wavelength roughly following $τ(λ) \propto λ^{4/3}$, the dependence expected for a geometrically thin, optically thick, centrally illuminated disk. Modeling of the light curves reveals a period at the beginning of the campaign where the response of the continuum is suppressed compared to later in the light curve - the light curves are not simple shifted and scaled versions of each other. The interval of suppressed response corresponds to a period of high UV line and X-ray absorption, and reduced emission line variability amplitudes. We suggest that this indicates a significant contribution to the continuum from the broad line region gas that sees an absorbed ionizing continuum. △ Less

Submitted 26 September, 2023; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: 20 pages, 13 figures, 3 tables, accepted for publication in ApJ

Hidden Cooling Flows in Clusters of Galaxies III: Accretion onto the Central Black Hole

Authors: A. C. Fabian, J. S. Sanders, G. J. Ferland, B. R. McNamara, C. Pinto, S. A. Walker

Abstract: 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 re… ▽ More 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. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 15 pages, 18 Figures, accepted for publication in MNRAS

AGN STORM 2: II. Ultraviolet Observations of Mrk817 with the Cosmic Origins Spectrograph on the Hubble Space Telescope

Authors: Y. Homayouni, Gisella De Rosa, Rachel Plesha, Gerard A. Kriss, Aaron J. Barth, Edward M. Cackett, Keith Horne, Erin A. Kara, Hermine Landt, Nahum Arav, Benjamin D. Boizelle, Misty C. Bentz, Thomas G. Brink, Michael S. Brotherton, Doron Chelouche, Elena Dalla Bonta, Maryam Dehghanian, Pu Du, Gary J. Ferland, Laura Ferrarese, Carina Fian, Alexei V. Filippenko, Travis Fischer, Ryan J. Foley, Jonathan Gelbord , et al. (40 additional authors not shown)

Abstract: We present reverberation mapping measurements for the prominent ultraviolet broad emission lines of the active galactic nucleus Mrk817 using 165 spectra obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our ultraviolet observations are accompanied by X-ray, optical, and near-infrared observations as part of the AGN Space Telescope and Optical Reverberation Mapping Progra… ▽ More We present reverberation mapping measurements for the prominent ultraviolet broad emission lines of the active galactic nucleus Mrk817 using 165 spectra obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our ultraviolet observations are accompanied by X-ray, optical, and near-infrared observations as part of the AGN Space Telescope and Optical Reverberation Mapping Program 2 (AGN STORM 2). Using the cross-correlation lag analysis method, we find significant correlated variations in the continuum and emission-line light curves. We measure rest-frame delayed responses between the far-ultraviolet continuum at 1180 A and Ly$α$ $\lambda1215$ A ($10.4_{-1.4}^{+1.6}$ days), N V $\lambda1240$ A ($15.5_{-4.8}^{+1.0}$days), SiIV + OIV] $\lambda1397$ A ($8.2_{-1.4}^{+1.4}$ days), CIV $\lambda1549$ A ($11.8_{-2.8}^{+3.0}$ days), and HeII $\lambda1640$ A ($9.0_{-1.9}^{+4.5}$ days) using segments of the emission-line profile that are unaffected by absorption and blending, which results in sampling different velocity ranges for each line. However, we find that the emission-line responses to continuum variations are more complex than a simple smoothed, shifted, and scaled version of the continuum light curve. We also measure velocity-resolved lags for the Ly$α$, and CIV emission lines. The lag profile in the blue wing of Ly$α$ is consistent with virial motion, with longer lags dominating at lower velocities, and shorter lags at higher velocities. The CIV lag profile shows the signature of a thick rotating disk, with the shortest lags in the wings, local peaks at $\pm$ 1500 $\rm km\,s^{-1}$, and a local minimum at line center. The other emission lines are dominated by broad absorption lines and blending with adjacent emission lines. These require detailed models, and will be presented in future work. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: Submitted to ApJ. 25 pages, 8 figures, and 6 tables

Testing AGN outflow and accretion models with CIV and HeII emission line demographics in z=2 quasars

Authors: Matthew J. Temple, James H. Matthews, Paul C. Hewett, Amy L. Rankine, Gordon T. Richards, Manda Banerji, Gary J. Ferland, Christian Knigge, Matthew Stepney

Abstract: Using 190,000 spectra from the seventeenth data release of the Sloan Digital Sky Survey, we investigate the ultraviolet emission line properties in z=2 quasars. Specifically, we quantify how the shape of CIV 1549A and the equivalent width (EW) of HeII 1640A depend on the black hole mass and Eddington ratio inferred from MgII 2800A. Above L/L_Edd>0.2, there is a strong mass dependence in both CIV b… ▽ More Using 190,000 spectra from the seventeenth data release of the Sloan Digital Sky Survey, we investigate the ultraviolet emission line properties in z=2 quasars. Specifically, we quantify how the shape of CIV 1549A and the equivalent width (EW) of HeII 1640A depend on the black hole mass and Eddington ratio inferred from MgII 2800A. Above L/L_Edd>0.2, there is a strong mass dependence in both CIV blueshift and HeII EW. Large CIV blueshifts are observed only in regions with both high mass and high accretion rate. Including X-ray measurements for a subsample of 5,000 objects, we interpret our observations in the context of AGN accretion and outflow mechanisms. The observed trends in HeII and 2 keV strength are broadly consistent with theoretical QSOSED models of AGN spectral energy distributions (SEDs) for low spin black holes, where the ionizing SED depends on the accretion disc temperature and the strength of the soft excess. High spin models are not consistent with observations, suggesting SDSS quasars at z=2 may in general have low spins. We find a dramatic switch in behaviour at L/L_Edd<0.2: the ultraviolet emission properties show much weaker trends, and no longer agree with QSOSED predictions, hinting at changes in the structure of the broad line region. Overall the observed emission line trends are generally consistent with predictions for radiation line driving where quasar outflows are governed by the SED, which itself results from the accretion flow and hence depends on both the SMBH mass and accretion rate. △ Less

Submitted 27 April, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: Resubmitted to MNRAS. Minor edits to text; figures updated; licence updated

Journal ref: MNRAS 523, 646 (2023)

Hidden Cooling Flows in Clusters of Galaxies II: A Wider Sample

Authors: A. C. Fabian, J. S. Sanders, G. J. Ferland, B. R. McNamara, C. Pinto, S. A. Walker

Abstract: 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 tw… ▽ More 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. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: 10 Pages, 13 figures, 3 tables, MNRAS submitted

Backscattering and Line Broadening in Orion

Authors: C. R. O'Dell, G. J. Ferland, J. E. Mendez-Delgado

Abstract: 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… ▽ More 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. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Radiation pressure-driven outflows from dusty AGN

Authors: N. Arakawa, A. C. Fabian, G. J. Ferland, W. Ishibashi

Abstract: 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… ▽ More 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 $λ_{\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}} - λ_{\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. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS; 11 pages, 7 figures

Creating a CLOUDY Compatible Database with CHIANTI version 10 Data

Authors: Chamani M. Gunasekera, Marios Chatzikos, Gary J. Ferland

Abstract: Atomic and molecular data are required to conduct detailed calculations of microphysical processes performed by CLOUDY to predict the spectra of a theoretical model. One of the three databases CLOUDY currently utilizes is CHIANTI version 7.1. CHIANTI version 10.0.1 is available, but its format has changed. CLOUDY is incompatible with the newer version. We have developed a script to convert the ver… ▽ More Atomic and molecular data are required to conduct detailed calculations of microphysical processes performed by CLOUDY to predict the spectra of a theoretical model. One of the three databases CLOUDY currently utilizes is CHIANTI version 7.1. CHIANTI version 10.0.1 is available, but its format has changed. CLOUDY is incompatible with the newer version. We have developed a script to convert the version 10.0.1 database into its version 7.1 format so that CLOUDY does not have to change every time there is a new CHIANTI version with an evolved format. This study outlines the steps taken by the script for this version format change. We have also found a modest number of significant changes to spectral line intensities/luminosities calculated by CLOUDY with the adoption of CHIANTI version 10.0.1. These changes are a result of improvements to collision strength data. △ Less

Submitted 2 December, 2022; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 16 pages, 7 figures

Journal ref: Astronomy 2022, 1, 255-270

Multiple locations of near-infrared coronal lines in NGC 5548

Authors: Daniel Kynoch, Hermine Landt, Maryam Dehghanian, Martin J. Ward, Gary J. Ferland

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 20 pages, 16 figures and 6 tables. Accepted for publication in MNRAS

Hidden Cooling Flows in Clusters of Galaxies

Authors: A. C. Fabian, G. J. Ferland, J. S. Sanders, B. R. McNamara, C. Pinto, S. A. Walker

Abstract: 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 c… ▽ More 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, △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: 11 pages. 7 figures, 2 tables MNRAS in press

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

Authors: Gargi Shaw, G. J. Ferland, M. Chatzikos

Abstract: 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 involvi… ▽ More 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. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 40 pages, 6 figures, 18 tables, Accepted for publication in ApJ

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

Authors: D. Souropanis, A. Chiotellis, P. Boumis, M. Chatzikos, S. Akras, L. Piersanti, A. J. Ruiter, G. J. Ferland

Abstract: 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 th… ▽ More 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. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 16 pages, 6 figures, accepted to MNRAS

A Practical Guide to the Partition Function of Atoms and Ions

Authors: P. Alimohamadi, G. J. Ferland

Abstract: 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 theo… ▽ More 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. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: Tutorial introduction to the partition function in astrophysics

Dense Molecular Clouds in the Crab Supernova Remnant

Authors: Alwyn Wootten, Rory O. Bentley, J. Baldwin, F. Combes, A. C. Fabian, G. J. Ferland, E. Loh, P. Salome, C. N. Shingledecker, A. Castro-Carrizo

Abstract: 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… ▽ More 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. △ Less

Submitted 10 November, 2021; originally announced November 2021.

Comments: ApJ in press, comments welcome

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

Authors: Hengxiao Guo, Aaron J. Barth, Kirk T. Korista, Michael R. Goad, Edward M. Cackett, Misty C. Bentz, William N. Brandt, D. Gonzalez-Buitrago, Gary J. Ferland, Jonathan M. Gelbord, Luis C. Ho, Keith Horne, Michael D. Joner, Gerard A. Kriss, Ian McHardy, Missagh Mehdipour, Daeseong Park, Raymond Remigio, Vivian U, Marianne Vestergaard

Abstract: 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… ▽ More 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. △ Less

Submitted 8 January, 2022; v1 submitted 4 November, 2021; originally announced November 2021.

Comments: 17 pages, 6 figures, and 3 tables; ApJ in press

H, He-like recombination spectra IV; H, He-like recombination spectra IV: clarification and refinement of methodology for $l$-changing collisions

Authors: N. R. Badnell, F. Guzmán, S. Brodie, R. J. R. Williams, P. A. M. van Hoof, M. Chatzikos, G. J. Ferland

Abstract: 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 Gu… ▽ More 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án \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. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: accepted in MNRAS

High-ionization emission line ratios from quasar broad line regions: metallicity or density?

Authors: Matthew J. Temple, Gary J. Ferland, Amy L. Rankine, Marios Chatzikos, Paul C. Hewett

Abstract: 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 f… ▽ More 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. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 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

Journal ref: MNRAS 505, 3247-3259 (2021)

AGN STORM 2: I. First results: A Change in the Weather of Mrk 817

Authors: Erin Kara, Missagh Mehdipour, Gerard A. Kriss, Edward M. Cackett, Nahum Arav, Aaron J. Barth, Doyee Byun, Michael S. Brotherton, Gisella De Rosa, Jonathan Gelbord, Juan V. Hernandez Santisteban, Chen Hu, Jelle Kaastra, Hermine Landt, Yan-Rong Li, Jake A. Miller, John Montano, Ethan Partington, Jesus Aceituno, Jin-Ming Bai, Dongwei Bao, Misty C. Bentz, Thomas G. Brink, Doron Chelouche, Yong-Jie Chen , et al. (47 additional authors not shown)

Abstract: We present the first results from the ongoing, intensive, multi-wavelength monitoring program of the luminous Seyfert 1 galaxy Mrk 817. While this AGN was, in part, selected for its historically unobscured nature, we discovered that the X-ray spectrum is highly absorbed, and there are new blueshifted, broad and narrow UV absorption lines, which suggest that a dust-free, ionized obscurer located at… ▽ More We present the first results from the ongoing, intensive, multi-wavelength monitoring program of the luminous Seyfert 1 galaxy Mrk 817. While this AGN was, in part, selected for its historically unobscured nature, we discovered that the X-ray spectrum is highly absorbed, and there are new blueshifted, broad and narrow UV absorption lines, which suggest that a dust-free, ionized obscurer located at the inner broad line region partially covers the central source. Despite the obscuration, we measure UV and optical continuum reverberation lags consistent with a centrally illuminated Shakura-Sunyaev thin accretion disk, and measure reverberation lags associated with the optical broad line region, as expected. However, in the first 55 days of the campaign, when the obscuration was becoming most extreme, we observe a de-coupling of the UV continuum and the UV broad emission line variability. The correlation recovers in the next 42 days of the campaign, as Mrk 817 enters a less obscured state. The short CIV and Ly alpha lags suggest that the accretion disk extends beyond the UV broad line region. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: 28 pages, 14 figures, submitted to ApJ. Comments welcome

Excitation mechanisms in the intracluster filaments surrounding Brightest Cluster Galaxies

Authors: F. L. Polles, P. Salomé, P. Guillard, B. Godard, G. Pineau des Forêts, V. Olivares, R. S. Beckmann, R. E. A. Canning, F. Combes, Y. Dubois, A. C. Edge, A. C. Fabian, G. J. Ferland, S. L. Hamer, M. D. Lehnert

Abstract: 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… ▽ More 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]λ 5007 Å/H\b{eta}, [N ii]λ 6583 Å/Hα and ([S ii]λ 6716 Å+[S ii]λ 6731 Å)/Hα 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. △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 16 pages, 9 figures, 3 tables. Accepted for publication in A&A

Deciphering the 3-D Orion Nebula-IV: The HH~269 flow emerges from the Orion-S Embedded Molecular Cloud

Authors: C. R. O'Dell, N. P. Abel, G. J. Ferland

Abstract: 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)… ▽ More 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. △ Less

Submitted 30 January, 2021; originally announced February 2021.

Comments: 7 pages, 3 figures

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

Authors: Gargi Shaw, G. J. Ferland

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

Submitted 12 January, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication in ApJ, 31 pages, 22 figures. A typo in the title is corrected

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

Authors: C. R. O'Dell, G. J. Ferland, N. P. Abel

Abstract: 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… ▽ More 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. △ Less

Submitted 7 January, 2021; originally announced January 2021.

Comments: 21pages, 7 figures Accepted Astrophysical Journal

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

Authors: C. R. O'Dell, N. P. Abel, G. J. Ferland

Abstract: 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,… ▽ More 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. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: 9 pages, 2 figures

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

Authors: M. Dehghanian, G. J. Ferland, B. M. Peterson, G. A. Kriss, K. T. Korista, M. R. Goad, M. Chatzikos, M. C. Bentz, F. Guzman, M. Mehdipour, G. De Rosa

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

Submitted 17 November, 2020; originally announced November 2020.

Comments: accepted for publication in ApJ

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

Authors: A. Sarkar, G. J. Ferland, M. Chatzikos, F. Guzmán, P. A. M. van Hoof, R. T. Smyth, C. A. Ramsbottom, F. P. Keenan, C. P. Ballance

Abstract: 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 av… ▽ More 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. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Cosmic ray dissociation of molecular hydrogen and dense cloud chemistry

Authors: Gargi Shaw, G. J. Ferland, S. Ploeckinger

Abstract: 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… ▽ More 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. △ Less

Submitted 30 October, 2020; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: corrected minor typo

Journal ref: 2020RNAAS...4...78S

Cloudy in the microcalorimeter era: improved energies for K$α$ transitions

Authors: P. Chakraborty, G. J. Ferland, S. Bianchi, M. Chatzikos

Abstract: X-ray missions with microcalorimeter technology will resolve spectral features with unprecedented detail. In this work, we improve the H-like K$α$ 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… ▽ More X-ray missions with microcalorimeter technology will resolve spectral features with unprecedented detail. In this work, we improve the H-like K$α$ 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. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Space Telescope and Optical Reverberation Mapping Project. XII. Broad-Line Region Modeling of NGC 5548

Authors: P. R. Williams, A. Pancoast, T. Treu, B. J. Brewer, B. M. Peterson, A. J. Barth, M. A. Malkan, G. De Rosa, Keith Horne, G. A. Kriss, N. Arav, M. C. Bentz, E. M. Cackett, E. Dalla Bontà, M. Dehghanian, C. Done, G. J. Ferland, C. J. Grier, J. Kaastra, E. Kara, C. S. Kochanek, S. Mathur, M. Mehdipour, R. W. Pogge, D. Proga , et al. (133 additional authors not shown)

Abstract: 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$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas traje… ▽ More 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$β$, C IV, and Ly$α$ broad emission lines. We find an extended disk-like H$β$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$α$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$α$ emission arising at smaller radii than the H$β$ 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$β$ 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. △ Less

Submitted 1 October, 2020; originally announced October 2020.

Comments: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

X-ray spectroscopy in the microcalorimeter era II: A new Diagnostic on Column Density from the Case A to B Transition in H- and He-like Iron

Authors: P. Chakraborty, M. Chatzikos, F. Guzmán, Y. Su, G. J. Ferland

Abstract: 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-ele… ▽ More 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$α$ 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$α$ 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. △ Less

Submitted 31 July, 2020; originally announced July 2020.

X-ray spectroscopy in the microcalorimeter era I: Effects of Fe XXIV Resonance Auger Destruction on Fe XXV K$α$ spectra

Authors: P. Chakraborty, G. J. Ferland, M. Chatzikos, F. Guzmán, Y. Su

Abstract: We discuss the importance of Fe$^{23+}$ in determining the line intensities of the Fe XXV K$α$ 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 c… ▽ More We discuss the importance of Fe$^{23+}$ in determining the line intensities of the Fe XXV K$α$ 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$α$ photon by RAD. Out of the four members in the Fe XXV K$α$ 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. △ Less

Submitted 30 July, 2020; originally announced July 2020.

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

Authors: Matthew J. Temple, Gary J. Ferland, Amy L. Rankine, Paul C. Hewett, N. R. Badnell, C. P. Ballance, G. Del Zanna, R. P. Dufresne

Abstract: 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… ▽ More 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. △ Less

Submitted 30 June, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 12 pages, 11 figures, plus appendices

Journal ref: MNRAS 496, 2565 (2020)

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

Authors: M. Dehghanian, G. J. Ferland, G. A. Kriss, B. M. Peterson, K. T. Korista, M. R. Goad, M. Chatzikos, F. Guzman, G. de Rosa, M. Mehdipour, J. Kaastra, S. Mathur, M. Vestergaard, D. Proga, T. Waters, M. C. Bentz, S. Bisogni, W. N. Brandt, E. Dalla Bont`a, M. M. Fausnaugh, J. M. Gelbord, Keith Horne, I. M. McHardy, R. W. Pogge, D. A. Starkey

Abstract: 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… ▽ More 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. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: Accepted for publication in ApJ

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

Authors: C. R. O'Dell, N. P. Abel, G. J. Ferland

Abstract: 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… ▽ More 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. △ Less

Submitted 3 March, 2020; originally announced March 2020.

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

Authors: Keith Horne, G. De Rosa, B. M. Peterson, A. J. Barth, J. Ely, M. M. Fausnaugh, G. A. Kriss, L. Pei, S. M. Adams, M. D. Anderson, P. Arevalo, T G. Beatty, V. N. Bennert, M. C. Bentz, A. Bigley, S. Bisogni, G. A. Borman, T. A. Boroson, M. C. Bottorff, W. N. Brandt, A. A. Breeveld, M. Brotherton, J. E. Brown, J. S. Brown, E. M. Cackett , et al. (133 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 27 November, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: 19 pages, 9 figures, ApJ in press

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

Authors: Gargi Shaw, G. J. Ferland

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

Submitted 1 March, 2020; originally announced March 2020.

Comments: 10 pages, 5 figures, 6 tables. Accepted for publication in MNRAS

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

Authors: Svea Hernandez, Alessandra Aloisi, Bethan L. James, Gary J. Ferland, Andrew J. Fox, Monica Tosi, Jason Tumlinson

Abstract: 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 io… ▽ More 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. △ Less

Submitted 18 February, 2020; originally announced February 2020.

Comments: 20 pages, 8 tables, 8 Figures, Accepted for Publication in ApJ

Current and future development of the photoionization code Cloudy

Authors: P. A. M. van Hoof, G. C. Van de Steene, F. Guzmán, M. Dehghanian, M. Chatzikos, G. J. Ferland

Abstract: 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… ▽ More 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. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: 12 pages, 5 figures. This paper was presented as an invited lecture at the 12th Serbian Conference on Spectral Line Shapes in Astrophysics in Vrdnik, Serbia

Journal ref: Contributions of the Astronomical Observatory Skalnaté Pleso, 2020, vol. 50, pp 32-43

2018 NASA Laboratory Astrophysics Workshop: Scientific Organizing Committee Report

Authors: N. Brickhouse, G. J. Ferland, S. Milam, E. Sciamma-O'Brien, A. Smale, A. Spyrou, P. C. Stancil, L. Storrie-Lombardi, G. M. Wahlgren

Abstract: 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 Worksh… ▽ More 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. △ Less

Submitted 28 January, 2020; originally announced January 2020.

A wind-based unification model for NGC 5548: spectral holidays, non-disk emission, and implications for changing-look quasars

Authors: M. Dehghanian, G. J. Ferland, B. M. Peterson, G. A. Kriss, K. T. Korista, M. Chatzikos, F. Guzman, N. Arav, G. De Rosa, M. R. Goad, M. Mehdipour, P. A. M. Van Hoof

Abstract: 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 thr… ▽ More 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. △ Less

Submitted 22 August, 2019; v1 submitted 20 August, 2019; originally announced August 2019.

Comments: Accepted for publication in ApJL