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

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

Physics of 1 keV line in X-ray binaries

Authors: Priyanka Chakraborty, Gary Ferland, Andrew Fabian, Arnab Sarkar, Renee Ludlam, Stefano Bianchi, Hayden Hall, Peter Kosec

Abstract: X-ray binaries (XRBs) often exhibit spectral residuals in the 0.5 to 2 keV range, known as the "1 keV residual/1 keV feature", with variable centroid and intensity across different systems. Yet a comprehensive scientific explanation of the variability of the 1 keV feature has remained largely elusive. In this paper, we explain for the first time the origin and variability of the 1 keV feature in X… ▽ More X-ray binaries (XRBs) often exhibit spectral residuals in the 0.5 to 2 keV range, known as the "1 keV residual/1 keV feature", with variable centroid and intensity across different systems. Yet a comprehensive scientific explanation of the variability of the 1 keV feature has remained largely elusive. In this paper, we explain for the first time the origin and variability of the 1 keV feature in XRBs using the spectral synthesis code \textsc{Cloudy}. We constructed line blends for the emission and absorption lines and study the variability of these blends with ionization parameters, temperature, and column density. We conducted a sample study involving five XRBs including two ultraluminous X-ray sources (ULXs): NGC 247 ULX-1, NGC 1313 X-1, a binary X-ray pulsar: Hercules X-1, and two typical low-mass X-ray binaries (LMXBs): Cygnus X-2, and Serpens X-1, providing a comprehensive explanation of the 1 keV feature observed across these targets. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 22 pages, 12 figures, Submitted to MNRAS

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.

JADES -- The small blue bump in GN-z11: insights into the nuclear region of a galaxy at z=10.6

Authors: Xihan Ji, Roberto Maiolino, Gary Ferland, Francesco D'Eugenio, Rachana Bhatawdekar, Stéphane Charlot, Jacopo Chevallard, Mirko Curti, Emma Curtis-Lake, Kevin Hainline, Zhiyuan Ji, Brant Robertson, Bruno Rodríguez Del Pino, Jan Scholtz, Sandro Tacchella, Christina C. Williams, Joris Witstok

Abstract: We report the detection of continuum excess in the rest-frame UV between 3000 Å and 3550 Å in the JWST/NIRSpec spectrum of GN-z11, a galaxy hosting an active galactic nucleus (AGN) at z = 10.603. The shape of the continuum excess resembles a Balmer continuum but has a break around 3546 Å in the rest frame, which is 100 Å bluewards to the Balmer limit at 3646 Å. A Balmer continuum model alone canno… ▽ More We report the detection of continuum excess in the rest-frame UV between 3000 Å and 3550 Å in the JWST/NIRSpec spectrum of GN-z11, a galaxy hosting an active galactic nucleus (AGN) at z = 10.603. The shape of the continuum excess resembles a Balmer continuum but has a break around 3546 Å in the rest frame, which is 100 Å bluewards to the Balmer limit at 3646 Å. A Balmer continuum model alone cannot fit the spectrum, implying a different origin for the continuum excess. The absence of the Balmer jump indicates an electron temperature of $\sim 3\times 10^4$ K, which is significantly higher than the temperature of $T_{e}({\rm O^{2+}}) \approx 1.3\times 10^{4}$ K inferred from [OIII]$λ4363$. The temperature difference must result from mixing of different ionized regions: the Balmer emission mainly arises from dense and hot clouds in the Broad Line Region, close to the accreting black hole, whereas the forbidden lines originate from less dense and colder gas in the host galaxy (although these ionized regions are kinematically similar in GN-z11 due to its small BH mass). We propose a potential explanation for the observed continuum excess to come from a complex of FeII emission, which shows a characteristic jump bluewards to the Balmer limit as previously seen in the spectra of many lower-redshift quasars. Through comparisons with Cloudy models, we show an Fe abundance or an overall metallicity above $\sim 1/3$ solar is likely needed. Besides the FeII emission, part of the small blue bump might also be associated with an OIII Bowen fluorescent line, a line often enhanced in dense AGN-ionized gas. Finally, the spectrum provides further evidence against Wolf-Rayet or massive stars dominating the nebular emission in GN-z11. △ Less

Submitted 20 May, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: 22 pages (including appendix), 18 figures, submitted to MNRAS

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

Two distinct molecular cloud populations detected in massive galaxies

Authors: Tom Rose, B. R. McNamara, F. Combes, A. C. Edge, M. McDonald, Ewan O'Sullivan, H. Russell, A. C. Fabian, G. Ferland, P. Salome, G. Tremblay

Abstract: We present new ALMA observations of CO, CN, CS, HCN and HCO$^{+}$ absorption seen against the bright and compact radio continuum sources of eight massive galaxies. Combined with archival observations, they reveal two distinct populations of molecular clouds, which we identify by combining CO emission and absorption profiles to unambiguously reveal each cloud's direction of motion and likely locati… ▽ More We present new ALMA observations of CO, CN, CS, HCN and HCO$^{+}$ absorption seen against the bright and compact radio continuum sources of eight massive galaxies. Combined with archival observations, they reveal two distinct populations of molecular clouds, which we identify by combining CO emission and absorption profiles to unambiguously reveal each cloud's direction of motion and likely location. In galaxy disks, we see clouds with low velocity dispersions, low line of sight velocities and a lack of any systemic inflow or outflow. In galactic cores, we find high velocity dispersion clouds inflowing at up to 550 km/s. This provides observational evidence in favour of cold accretion onto galactic centres, which likely contributes to the fuelling of active galactic nuclei. We also see a wide range in the CO(2-1)/CO(1-0) ratios of the absorption lines. This is likely the combined effect of hierarchical substructure within the molecular clouds and continuum sources which vary in size with frequency. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: Submitted to MNRAS

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

Revisiting the gas-phase chemical rate coefficients at high temperatures in CLOUDY

Authors: Gargi Shaw, Gary Ferland, M. Chatzikos

Abstract: A two-body gas-phase reaction rate coefficient can be given by the usual Arrhenius-type formula which depends on temperature. The UMIST Database for Astrochemistry is a widely used database for reaction rate coefficients. They provide fittings for coefficients valid over a particular range of temperatures. The permissible upper-temperature limits vary over a wide range: from 100 K to 41000K. A wid… ▽ More A two-body gas-phase reaction rate coefficient can be given by the usual Arrhenius-type formula which depends on temperature. The UMIST Database for Astrochemistry is a widely used database for reaction rate coefficients. They provide fittings for coefficients valid over a particular range of temperatures. The permissible upper-temperature limits vary over a wide range: from 100 K to 41000K. A wide range of temperatures occurs in nature; thus, it requires evaluating the rate coefficients at temperatures outside the range of validity. As a result, a simple extrapolation of the rate coefficients can lead to unphysically large values at high temperatures. These result in unrealistic predictions. Here we present a solution to prevent the gas-phase reaction coefficients from diverging at a very high temperature. We implement this into the spectral synthesis code CLOUDY which operates over a wide range of temperatures from CMB to 10$^{10}$ K subject to different astrophysical environments. △ Less

Submitted 20 July, 2023; originally announced August 2023.

Comments: Submitted, comments welcome. Three pages, one figure

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

Recent updates of gas-phase chemical reactions and molecular lines of SiS in CLOUDY

Authors: Gargi Shaw, Gary Ferland, M. Chatzikos

Abstract: Here we present our current update of CLOUDY on gas-phase chemical reactions for the formation and destruction of the SiS molecule, its energy levels, and collisional rate coefficients with H$_2$, H, and He over a wide range of temperatures. As a result, henceforth the spectral synthesis code CLOUDY predicts SiS line intensities and column densities for various astrophysical environments. Here we present our current update of CLOUDY on gas-phase chemical reactions for the formation and destruction of the SiS molecule, its energy levels, and collisional rate coefficients with H$_2$, H, and He over a wide range of temperatures. As a result, henceforth the spectral synthesis code CLOUDY predicts SiS line intensities and column densities for various astrophysical environments. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: Submitted, comments welcome. Three pages, one figure

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.

Does absorption against AGN reveal supermassive black hole accretion?

Authors: Tom Rose, B. R. McNamara, F. Combes, A. C. Edge, A. C. Fabian, M. Gaspari, H. Russell, P. Salomé, G. Tremblay, G. Ferland

Abstract: Galaxies often contain large reservoirs of molecular gas which shape their evolution. This can be through cooling of the gas -- which leads to star formation, or accretion onto the central supermassive black hole -- which fuels AGN activity and produces powerful feedback. Molecular gas has been detected in early-type galaxies on scales of just a few tens to hundreds of solar masses by searching fo… ▽ More Galaxies often contain large reservoirs of molecular gas which shape their evolution. This can be through cooling of the gas -- which leads to star formation, or accretion onto the central supermassive black hole -- which fuels AGN activity and produces powerful feedback. Molecular gas has been detected in early-type galaxies on scales of just a few tens to hundreds of solar masses by searching for absorption against their compact radio cores. Using this technique, ALMA has found absorption in several brightest cluster galaxies, some of which show molecular gas moving towards their galaxy's core at hundreds of km/s. In this paper we constrain the location of this absorbing gas by comparing each galaxy's molecular emission and absorption. In four galaxies, the absorption properties are consistent with chance alignments between the continuum and a fraction of the molecular clouds visible in emission. In four others, the properties of the absorption are inconsistent with this scenario. In these systems the absorption is likely produced by a separate population of molecular clouds in close proximity to the galaxy core and with high inward velocities and velocity dispersions. We thus deduce the existence of two types of absorber, caused by chance alignments between the radio core and: (i) a fraction of the molecular clouds visible in emission, and (ii) molecular clouds close to the AGN, in the process of accretion. We also present the first ALMA observations of molecular emission in S555, Abell 2390, RXC J1350.3+0940 and RXC J1603.6+1553 -- with the latter three having molecular masses of $>10^{10}$M$_{\odot}$. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: Submitted to MNRAS

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

X-ray spectroscopy in the microcalorimeter era 4: Optical depth effects on the soft X-rays studied with CLOUDY

Authors: Priyanka Chakraborty, Gary Ferland, Marios Chatzikos, Andrew Fabian, Stefano Bianchi, Francisco Guzmán, Yuanyuan Su

Abstract: In this paper, we discuss atomic processes modifying the soft X-ray spectra in from optical depth effects like photoelectric absorption and electron scattering suppressing the soft X-ray lines. We also show the enhancement in soft X-ray line intensities in a photoionized environment via continuum pumping. We quantify the suppression/enhancement by introducing a "line modification factor (… ▽ More In this paper, we discuss atomic processes modifying the soft X-ray spectra in from optical depth effects like photoelectric absorption and electron scattering suppressing the soft X-ray lines. We also show the enhancement in soft X-ray line intensities in a photoionized environment via continuum pumping. We quantify the suppression/enhancement by introducing a "line modification factor ($f_{\rm mod}$)." If 0 $\leq$ $f_{\rm mod}$ $\leq$ 1, the line is suppressed, which could be the case in both collisionally-ionized and photoionized systems. If $f_{\rm mod}$ $\geq$ 1, the line is enhanced, which occurs in photoionized systems. Hybrid astrophysical sources are also very common, where the environment is partly photoionized and partly collisionally-ionized. Such a system is V1223 Sgr, an intermediate polar binary. We show the application of our theory by fitting the first-order Chandra MEG spectrum of V1223 Sgr with a combination of \textsc{Cloudy}-simulated additive cooling-flow and photoionized models. In particular, we account for the excess flux for O~VII, O~VIII, Ne~IX, Ne~X, and Mg~XI lines in the spectrum found in a recent study, which could not be explained with an absorbed cooling-flow model. △ Less

Submitted 5 July, 2022; originally announced July 2022.

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

Self-Consistent Grain Depletions and Abundances II: Effects on strong-line diagnostics of extragalactic H II regions

Authors: Chamani M. Gunasekera, Xihan Ji, Marios Chatzikos, Renbin Yan, Gary Ferland

Abstract: The depletion of elements onto dust grains is characterized using a generalized depletion strength $F_*$ for any sightline, and trend-line parameters $A_X, B_X$ and $z_X$. The parameters $A_X, B_X$ and $z_X$ define the relative depletion pattern, for which values are published in previous works. The present study uses these parameters to calculate post-depleted gas-phase abundances of 15 different… ▽ More The depletion of elements onto dust grains is characterized using a generalized depletion strength $F_*$ for any sightline, and trend-line parameters $A_X, B_X$ and $z_X$. The parameters $A_X, B_X$ and $z_X$ define the relative depletion pattern, for which values are published in previous works. The present study uses these parameters to calculate post-depleted gas-phase abundances of 15 different elements while varying $F_*$ from 0 to 1. An analysis of emergent strong spectral line intensities, obtained by inputting the calculated abundances into a cloudy model, shows that the depletion strength has a non-trivial effect on predicted emission lines and the thermal balance of the ionized cloud. The amount by which elements deplete also affects the coolant abundances in the gas. Furthermore, it was found that each of the parameters - metallicity, ionization parameter U and depletion strength $F_*$ have degenerate effects on the emission-line strengths, and thermal balance of the interstellar medium (ISM). Finally, comparing our results to a sample of H II regions using data obtained from the Mapping Nearby Galaxies at Apache Point Observatory survey (MaNGA) revealed that the best-fit $F_*$ was approximately 0.5. However, this best-fit value does not work well for all metallicities. Removing the sulfur depletion and changing the nitrogen abundance pattern can improve the fit. As a result, extra observational evidence is required to verify the choices of parameters and better constrain the typical depletion strength in galaxies. △ Less

Submitted 27 March, 2023; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 11 pages, 5 figures

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

Gas condensation in Brightest Group Galaxies unveiled with MUSE

Authors: V. Olivares, P. Salome, S. L. Hamer, F. Combes, M. Gaspari, K. Kolokythas, E. O'Sullivan, R. S. Beckmann, A. Babul, F. L. Polles, M. Lehnert, S. I. Loubser, M. Donahue, M. -L. Gendron-Marsolais, P. Lagos, G. Pineau des Forets, B. Godard, T. Rose, G. Tremblay, G. Ferland, P. Guillard

Abstract: The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extendi… ▽ More The origin of the cold gas in central galaxies in groups is still a matter of debate. We present Multi-Unit Spectroscopic Explorer (MUSE) observations of 18 optically selected local Brightest Group Galaxies (BGGs) to study the kinematics and distribution of the optical emission-line gas. MUSE observations reveal a distribution of gas morphologies including ten complex networks of filaments extending up to 10 kpc to two compact (<3 kpc) and five extended (>5 kpc) disk-dominated structures. Some rotating disks show rings and elongated structures arising from the central disk. The kinematics of the stellar component is mainly rotation-dominated, which is very different from the disturbed kinematics and distribution found in the filamentary sources. The ionized gas is kinematically decoupled from the stellar component for most systems, suggesting an external origin for the gas. We also find that the Halpha luminosity correlates with the cold molecular mass. By exploring the thermodynamical properties of the hot atmospheres, we find that the filamentary sources and compact disks are found in systems with small central entropy values and tcool/teddy ratios. This suggests that, like for Brightest Cluster Galaxies in cool core clusters, the ionized gas are likely formed from hot halo gas condensations, consistently with the Chaotic Cold Accretion simulations (as shown via the C-ratio, Tat, and k-plot). We note that gaseous rotating disks are more frequent than in BCGs. An explanation for the origin of the gas in those objects is a contribution to gas fueling by mergers or group satellites, as qualitatively hinted by some sources of the present sample. Nonetheless, we discuss the possibility that some extended disks could also be a transition stage in an evolutionary sequence including filaments, extended disks and compact disks, as described by hot gas condensation models of cooling flows. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 38 pages, and 16 figures. Resubmitted to A&A after referee report

Journal ref: A&A 666, A94 (2022)

Self-consistent grain depletions and abundances I: The Orion Nebula as a test case

Authors: Chamani Gunasekera, Xihan Ji, Marios Chatzikos, Renbin Yan, Gary Ferland

Abstract: Atomic species in the interstellar medium (ISM) transition out of their gas phase mainly by depletion onto dust. In this study, we examine if there is any change to the spectral line ratio predictions from a photoionization model of the Orion H II region when the degree of dust depletions is altered according to the most recently published model. We use equations and parameters published by previo… ▽ More Atomic species in the interstellar medium (ISM) transition out of their gas phase mainly by depletion onto dust. In this study, we examine if there is any change to the spectral line ratio predictions from a photoionization model of the Orion H II region when the degree of dust depletions is altered according to the most recently published model. We use equations and parameters published by previous works, in order to streamline the calculation of depleted abundances within CLOUDY. Our aim is for CLOUDY users to be able to vary the level of depletion using a single parameter in the input file. This makes it possible to explore predictions for a large range of depletions more efficiently. Finally, we discuss the results obtained for a model of the Orion Nebula when the degree of depletions are manipulated in this way. We found that the intensity of line ratios are significantly affected by depletions onto dust grains. Further, we found that adjusting dust abundances along with depletion affects the structure and the overall temperature of the H$^+$ layer across the H II region. △ Less

Submitted 30 March, 2022; v1 submitted 8 January, 2022; originally announced January 2022.

Comments: 8 pages, 5 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 512, Issue 2, May 2022, Pages 2310-2317

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

Evidence of a tidal disruption event in GSN 069 from the abnormal carbon and nitrogen abundance ratio

Authors: Zhenfeng Sheng, Tinggui Wang, Gary Ferland, Xinwen Shu, Chenwei Yang, Ning Jiang, Yang Chen

Abstract: GSN 069 is an ultra-soft X-ray active galactic nucleus that previously exhibited a huge X-ray outburst and a subsequent long-term decay. It has recently presented X-ray quasi-periodic eruptions (QPEs). We report the detection of strong nitrogen lines but weak or undetectable carbon lines in its far ultraviolet spectrum. With a detailed photoionization model, we use the C IV/N IV] ratio and other r… ▽ More GSN 069 is an ultra-soft X-ray active galactic nucleus that previously exhibited a huge X-ray outburst and a subsequent long-term decay. It has recently presented X-ray quasi-periodic eruptions (QPEs). We report the detection of strong nitrogen lines but weak or undetectable carbon lines in its far ultraviolet spectrum. With a detailed photoionization model, we use the C IV/N IV] ratio and other ratios between nitrogen lines to constrain the [C/N] abundance of GSN 069 to be from -3.33 to -1.91. We argue that a partially disrupted red giant star can naturally explain the abnormal C/N abundance in the UV spectrum, while the surviving core orbiting the black hole might produce the QPEs. △ Less

Submitted 12 October, 2021; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 6 pages,5 figures, published in ApJL

Journal ref: ApJL 920 L25 (2021)

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

Cloudy in the microcalorimeter era: improved energies for Si and S K$α$ fluorescence lines

Authors: Francesco Camilloni, Stefano Bianchi, Roberta Amato, Gary Ferland, Victoria Grinberg

Abstract: The upcoming X-ray missions based on the microcalorimeter technology require exquisite precision in spectral simulation codes in order to match the unprecedented spectral resolution. In this work, we improve the fluorescence K$α$ energies for Si II-XI and S II-XIII in the code Cloudy. In particular, we provide here a patch to update the Cloudy fluorescence energy table, originally based on Kaastra… ▽ More The upcoming X-ray missions based on the microcalorimeter technology require exquisite precision in spectral simulation codes in order to match the unprecedented spectral resolution. In this work, we improve the fluorescence K$α$ energies for Si II-XI and S II-XIII in the code Cloudy. In particular, we provide here a patch to update the Cloudy fluorescence energy table, originally based on Kaastra & Mewe (1993), with the laboratory energies measured by Hell et al. (2016). The new Cloudy simulations were used to model the Chandra/HETG spectra of the High Mass X-ray Binary Vela X-1 previously presented in Amato et al. (2021), showing a remarkable agreement and a dramatic improvement with respect to the current release version of Cloudy (C17.02). △ Less

Submitted 5 July, 2021; originally announced July 2021.

Comments: Published on RNAAS

Journal ref: Res. Notes AAS 5 149 (2021)

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

X-ray spectroscopy in the microcalorimeter era III: line formation under Case A, Case B, Case C, and Case D in H- and He-like iron for a photoionized cloud

Authors: Priyanka Chakraborty, Gary Ferland, Marios Chatzikos, Francisco Guzmán, Yuanyuan Su

Abstract: Future microcalorimeter X-ray observations will resolve spectral features in unmatched detail. Understanding the line formation processes in the X-rays deserves much attention. The purpose of this paper is to discuss such processes in the presence of a photoionizing source. Line formation processes in one and two-electron species are broadly categorized into four cases. Case A occurs when the Lyma… ▽ More Future microcalorimeter X-ray observations will resolve spectral features in unmatched detail. Understanding the line formation processes in the X-rays deserves much attention. The purpose of this paper is to discuss such processes in the presence of a photoionizing source. Line formation processes in one and two-electron species are broadly categorized into four cases. Case A occurs when the Lyman line optical depths are very small and photoexcitation does not occur. Line photons escape the cloud without any scattering. Case B occurs when the Lyman-line optical depths are large enough for photons to undergo multiple scatterings. Case C occurs when a broadband continuum source strikes an optically thin cloud. The Lyman lines are enhanced by induced radiative excitation of the atoms/ions by continuum photons, also known as continuum pumping. A fourth less-studied scenario, where the Case B spectrum is enhanced by continuum pumping, is called Case D. Here, we establish the mathematical foundation of Cases A, B, C, and D in an irradiated cloud with Cloudy. We also show the total X-ray emission spectrum for all four cases within the energy range 0.1 - 10 keV at the resolving power of XRISM around 6 keV. Additionally, we show that a combined effect of electron scattering and partial blockage of continuum pumping reduces the resonance line intensities. Such reduction increases with column density and can serve as an important tool to measure the column density/optical depth of the cloud. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: accepted for publication in ApJ, comments are welcome

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