Showing 1–41 of 41 results for author: Chatzikos, M

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

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

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

HST viewing of spectacular star-forming trails behind ESO 137-001

Authors: William Waldron, Ming Sun, Rongxin Luo, Sunil Laudari, Marios Chatzikos, Suresh Sivanandam, Jeffrey D. P. Kenney, Pavel Jachym, G. Mark Voit, Megan Donahue, Matteo Fossati

Abstract: We present the results from the HST WFC3 and ACS data on an archetypal galaxy undergoing ram pressure stripping (RPS), ESO 137-001, in the nearby cluster Abell 3627. ESO 137-001 is known to host a prominent stripped tail detected in many bands from X-rays, Halpha to CO. The HST data reveal significant features indicative of RPS such as asymmetric dust distribution and surface brightness as well as… ▽ More We present the results from the HST WFC3 and ACS data on an archetypal galaxy undergoing ram pressure stripping (RPS), ESO 137-001, in the nearby cluster Abell 3627. ESO 137-001 is known to host a prominent stripped tail detected in many bands from X-rays, Halpha to CO. The HST data reveal significant features indicative of RPS such as asymmetric dust distribution and surface brightness as well as many blue young star complexes in the tail. We study the correlation between the blue young star complexes from HST, HII regions from Halpha (MUSE) and dense molecular clouds from CO (ALMA). The correlation between the HST blue star clusters and the HII regions is very good, while their correlation with the dense CO clumps are typically not good, presumably due in part to evolutionary effects. In comparison to the Starburst99+Cloudy model, many blue regions are found to be young (< 10 Myr) and the total star formation (SF) rate in the tail is 0.3 - 0.6 M_Sun/yr for sources measured with ages less than 100 Myr, about 40% of the SF rate in the galaxy. We trace SF over at least 100 Myr and give a full picture of the recent SF history in the tail. We also demonstrate the importance of including nebular emissions and a nebular to stellar extinction correction factor when comparing the model to the broadband data. Our work on ESO 137-001 demonstrates the importance of HST data for constraining the SF history in stripped tails. △ Less

Submitted 9 May, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 22 pages, 20 figures, 4 tables, published in MNRAS 06/2023

Journal ref: MNRAS 522 (2023) 173-194

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

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

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

ESO 137-002: a large spiral undergoing edge-on ram-pressure stripping with little star formation in the tail

Authors: Sunil Laudari, Pavel Jáchym, Ming Sun, Will Waldron, Marios Chatzikos, Jeffrey Kenney, Rongxin Luo, Paul Nulsen, Craig Sarazin, Françoise Combes, Tim Edge, G. Mark Voit, Megan Donahue, Luca Cortese

Abstract: Ram pressure stripping (RPS) is an important mechanism for galaxy evolution. In this work, we present results from HST and APEX observations of one RPS galaxy, ESO 137-002 in the closest rich cluster Abell 3627. The galaxy is known to host prominent X-ray and H$α$ tails. The HST data reveal significant features indicative of RPS in the galaxy, including asymmetric distribution of dust in the galax… ▽ More Ram pressure stripping (RPS) is an important mechanism for galaxy evolution. In this work, we present results from HST and APEX observations of one RPS galaxy, ESO 137-002 in the closest rich cluster Abell 3627. The galaxy is known to host prominent X-ray and H$α$ tails. The HST data reveal significant features indicative of RPS in the galaxy, including asymmetric distribution of dust in the galaxy, dust filaments and dust clouds in ablation generally aligned with the direction of ram pressure, and young star clusters immediately upstream of the residual dust clouds that suggest star formation (SF) triggered by RPS. The distribution of the molecular gas is asymmetric in the galaxy, with no CO upstream and abundant CO downstream and in the inner tail region. A total amount of $\sim 5.5 \times 10^{9}$ M$_\odot$ of molecular gas is detected in the galaxy and its tail. On the other hand, we do not detect any active SF in the X-ray and H$α$ tails of ESO 137-002 with the HST data and place a limit on the SF efficiency in the tail. Hence, if selected by SF behind the galaxy in the optical or UV (e.g., surveys like GASP or using the Galex data), ESO 137-002 will not be considered a ``jellyfish'' galaxy. Thus, galaxies like ESO 137-002 are important for our comprehensive understanding of RPS galaxies and the evolution of the stripped material. ESO 137-002 also presents a great example of an edge-on galaxy experiencing a nearly edge-on RPS wind. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: 19 pages, 16 figures, 5 tables, MNRAS, re-submitted

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

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

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

Submitted 20 August, 2021; originally announced August 2021.

Comments: accepted in MNRAS

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

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

Abstract: The flux ratios of high-ionization lines are commonly assumed to indicate the metallicity of the broad emission line region in luminous quasars. When accounting for the variation in their kinematic profiles, we show that the NV/CIV, (SiIV+OIV])/CIV and NV/Lya line ratios do not vary as a function of the quasar continuum luminosity, black hole mass, or accretion rate. Using photoionization models f… ▽ More The flux ratios of high-ionization lines are commonly assumed to indicate the metallicity of the broad emission line region in luminous quasars. When accounting for the variation in their kinematic profiles, we show that the NV/CIV, (SiIV+OIV])/CIV and NV/Lya line ratios do not vary as a function of the quasar continuum luminosity, black hole mass, or accretion rate. Using photoionization models from CLOUDY , we further show that the observed changes in these line ratios can be explained by emission from gas with solar abundances, if the physical conditions of the emitting gas are allowed to vary over a broad range of densities and ionizing fluxes. The diversity of broad line emission in quasar spectra can be explained by a model with emission from two kinematically distinct regions, where the line ratios suggest that these regions have either very different metallicity or density. Both simplicity and current galaxy evolution models suggest that near-solar abundances, with parts of the spectrum forming in high-density clouds, are more likely. Within this paradigm, objects with stronger outflow signatures show stronger emission from gas which is denser and located closer to the ionizing source, at radii consistent with simulations of line-driven disc-winds. Studies using broad-line ratios to infer chemical enrichment histories should consider changes in density and ionizing flux before estimating metallicities. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 12 pages, 7 figures, accepted for publication in MNRAS. The interactive version of Figs. 3 and 4 will be made available as supplementary online material with the journal, and will be available at https://people.ast.cam.ac.uk/~alr53/CIV_yvar/all_plots.html until August 2021

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

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

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.

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

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

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

Submitted 16 October, 2020; originally announced October 2020.

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

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

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

Submitted 31 July, 2020; originally announced July 2020.

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

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

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

Submitted 30 July, 2020; originally announced July 2020.

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

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

Abstract: In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper… ▽ More In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper showed that changes in the LOS obscurer produce the decorrelation between the absorption lines and the continuum. A second study showed that the base of the wind shields the BLR, leading to the emission-line decorrelation. In that study, we proposed the wind is normally transparent with no effect on the spectrum. Changes in the wind properties alter its shielding and affect the SED striking the BLR, producing the observed decorrelations. In this work, we investigate the impact of a translucent wind on the emission lines. We simulate the obscuration using XMM-Newton, NuSTAR, and HST observations to determine the physical characteristics of the wind. We find that a translucent wind can contribute a part of the He II and Fe K? emission. It has a modest optical depth to electron scattering, which explains the fainter far-side emission in the observed velocity delay maps. The wind produces the very broad base seen in the UV emission lines and may also be present in the Fe K? line. Our results highlight the importance of accounting for the effects of such winds in the analysis of the physics of the central engine. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: Accepted for publication in ApJ

Current and future development of the photoionization code Cloudy

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

Abstract: The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical… ▽ More The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical state of the gas and predicts the spectrum that it emits. Cloudy is continually being developed to improve the treatment of the microphysical processes and the database of fundamental data that it uses. In this paper we will discuss how we are developing the code to improve our high-density predictions by implementing better collisional-radiative models for all ions. We will also briefly discuss the experimental mode in Cloudy to model gas that is not in steady-state equilibrium and present a preliminary model of recombining gas in a planetary nebula that is on the cooling track. We finish with a short discussion of how we are speeding up the code by using parallelization. △ Less

Submitted 13 February, 2020; originally announced February 2020.

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

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

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

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

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

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

Comments: Accepted for publication in ApJL

H-, He-like recombination spectra III: $n$-changing collisions in highly-excited Rydberg states and their impact on the radio, IR and optical recombination lines

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

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

Submitted 18 March, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: 16 pages (accepted by MNRAS, replaced)

Space Telescope and Optical Reverberation Mapping Project. X. Understanding the Absorption-Line Holiday in NGC 5548

Authors: M. Dehghanian, G. J. Ferland, G. A. Kriss, B. M. Peterson, S. Mathur, M. Mehdipour, F. Guzman, M. Chatzikos, P. A. M. Van Hoof, R. J. R. Williams, N. Arav, A. J. Barth, M. C. Bentz, S. Bisogni, W. N. Brandt, D. M. Crenshaw, E. Dalla Bonta, G. De Rosa, M. M. Fausnaugh, J. M. Gelbord, M. R. Goad, A. Gupta, Keith Horne, J. Kaastra, C. Knigge , et al. (5 additional authors not shown)

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

Submitted 24 April, 2019; v1 submitted 30 December, 2018; originally announced December 2018.

Comments: This is the version accepted by ApJ. Will be published shortly

Challenges and Techniques for Simulating Line Emission

Authors: Karen P. Olsen, Andrea Pallottini, Aida Wofford, Marios Chatzikos, Mitchell Revalski, Francisco Guzmán, Gergö Popping, Enrique Vázquez-Semadeni, Georgios E. Magdis, Mark L. A. Richardson, Michaela Hirschmann, William J. Gray

Abstract: Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodol… ▽ More Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale. △ Less

Submitted 24 August, 2018; originally announced August 2018.

Comments: 29 pages, 10 figures, 2 tables. Report of the findings at the Walking the Line workshop held in March 2018. Submitted to the MDPI journal Galaxies

The [CII] 158 micron line emission in high-redshift galaxies

Authors: G. Lagache, M. Cousin, M. Chatzikos

Abstract: The [CII] fine structure transition at 158 microns is the dominant cooling line of cool interstellar gas, and is the brightest of emission lines from star forming galaxies from FIR through meter wavelengths. With the advent of ALMA and NOEMA, capable of detecting [CII]-line emission in high-redshift galaxies, there has been a growing interest in using the [CII] line as a probe of the physical cond… ▽ More The [CII] fine structure transition at 158 microns is the dominant cooling line of cool interstellar gas, and is the brightest of emission lines from star forming galaxies from FIR through meter wavelengths. With the advent of ALMA and NOEMA, capable of detecting [CII]-line emission in high-redshift galaxies, there has been a growing interest in using the [CII] line as a probe of the physical conditions of the gas in galaxies, and as a SFR indicator at z>4. In this paper, we use a semi-analytical model of galaxy evolution (G.A.S.) combined with the code CLOUDY to predict the [CII] luminosity of a large number of galaxies at 4< z<8. At such high redshift, the CMB represents a strong background and we discuss its effects on the luminosity of the [CII] line. We study the LCII-SFR and LCII-Zg relations and show that they do not strongly evolve with redshift from z=4 and to z=8. Galaxies with higher [CII] luminosities tend to have higher metallicities and higher star formation rates but the correlations are very broad, with a scatter of about 0.5 dex for LCII-SFR. Our model reproduces the LCII-SFR relations observed in high-redshift star-forming galaxies, with [CII] luminosities lower than expected from local LCII-SFR relations. Accordingly, the local observed LCII-SFR relation does not apply at high-z. Our model naturally produces the [CII] deficit, which appears to be strongly correlated with the intensity of the radiation field in our simulated galaxies. We then predict the [CII] luminosity function, and show that it has a power law form in the range of LCII probed by the model with a slope alpha=1. The slope is not evolving from z=4 to z=8 but the number density of [CII]-emitters decreases by a factor of 20x. We discuss our predictions in the context of current observational estimates on both the differential and cumulative luminosity functions. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: In press, will appear in Astronomy & Astrophysics

Journal ref: A&A 609, A130 (2018)

The 2017 Release of Cloudy

Authors: G. J. Ferland, M. Chatzikos, F. Guzmán, M. L. Lykins, P. A. M. van Hoof, R. J. R. Williams, N. P. Abel, N. R. Badnell, F. P. Keenan, R. L. Porter, P. C. Stancil

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

Submitted 5 October, 2023; v1 submitted 30 May, 2017; originally announced May 2017.

Comments: replace Fig 6 with the published version

Journal ref: 2017RMxAA..53..385F

Thermodynamically-consistent semi-classical $\ell$-changing rates

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

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

Submitted 27 April, 2017; originally announced April 2017.

Comments: 12 pages, 4 figures. To be published in J Phys B. Author original version: journal accepted version, with minor changes, will be posted after embargo period

Journal ref: J. Phys. B: At. Mol. Opt. Phys. (2017) 50, 115201

Testing atomic collision theory with the two-photon continuum of astrophysical nebulae

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

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

Submitted 26 January, 2017; originally announced January 2017.

Comments: 7 pages, 6 fig, 2 tables, accepted in MNRAS

Journal ref: Mon Not R Astron Soc (2017) 467 (4): 3944-3950

H, He-like recombination spectra II: $l$-changing collisions for He Rydberg states

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

Abstract: Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which foc… ▽ More Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $\sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy. △ Less

Submitted 31 August, 2016; originally announced August 2016.

Comments: 9 pages, $ figures, 2 tables, accepted for publication in MNRAS

H, He-like recombination spectra I: $l$-changing collisions for hydrogen

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

Abstract: Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe di… ▽ More Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe dipole $l$-changing collisions ($nl\rightarrow nl'=l\pm 1$) between energy-degenerate states within an $n$-shell. The work of Pengelly & Seaton (1964) has, for half-a-century, been considered the definitive study which "solved" the problem. Recent work by Vrinceanu et al.(2012) recommended the use of rate coefficients from a semi-classical approximation which are nearly an order of magnitude smaller than those of Pengelly & Seaton (1964), with the result that significantly higher densities are needed for the $nl$ populations to come into local thermodynamic equilibrium. Here, we compare predicted H~I emissivities from the two works and find widespread differences, of up to $\approx 10$%. This far exceeds the 1% precision required to obtain the primordial He/H abundance ratio from observations so as to constrain Big Bang cosmologies. We recommend using the rate coefficients of Pengelly & Seaton (1964) for $l$-changing collisions, to describe the H recombination spectrum, based-on their quantum mechanical representation of the long-range dipole interaction. △ Less

Submitted 13 April, 2016; originally announced April 2016.

Comments: 6 pages, 3 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

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

Authors: Sibasish Laha, Francis P. Keenan, Gary J. ferland, Catherine A. Ramsbottom, Kanti M. Aggarwal, Thomas R. Ayres, Marios Chatzikos, Peter A. M. van Hoof, Robin J. R. Williams

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

Submitted 29 October, 2015; originally announced October 2015.

Comments: Accepted for publication in MNRAS

Deep Chandra, HST-COS, and Megacam Observations of the Phoenix Cluster: Extreme Star Formation and AGN Feedback on Hundred Kiloparsec Scales

Authors: M. McDonald, B. R. McNamara, R. J. van Weeren, D. E. Applegate, M. Bayliss, M. W. Bautz, B. A. Benson, J. E. Carlstrom, L. E. Bleem, M. Chatzikos, A. C. Edge, A. C. Fabian, G. P. Garmire, J. Hlavacek-Larrondo, C. Jones-Forman, A. B. Mantz, E. D. Miller, B. Stalder, S. Veilleux, J. A. Zuhone

Abstract: We present new ultraviolet, optical, and X-ray data on the Phoenix galaxy cluster (SPT-CLJ2344-4243). Deep optical imaging reveals previously-undetected filaments of star formation, extending to radii of ~50-100 kpc in multiple directions. Combined UV-optical spectroscopy of the central galaxy reveals a massive (2x10^9 Msun)), young (~4.5 Myr) population of stars, consistent with a time-averaged s… ▽ More We present new ultraviolet, optical, and X-ray data on the Phoenix galaxy cluster (SPT-CLJ2344-4243). Deep optical imaging reveals previously-undetected filaments of star formation, extending to radii of ~50-100 kpc in multiple directions. Combined UV-optical spectroscopy of the central galaxy reveals a massive (2x10^9 Msun)), young (~4.5 Myr) population of stars, consistent with a time-averaged star formation rate of 610 +/- 50 Msun/yr. We report a strong detection of OVI(1032,1038) which appears to originate primarily in shock-heated gas, but may contain a substantial contribution (>1000 Msun/yr) from the cooling intracluster medium. We confirm the presence of deep X-ray cavities in the inner ~10 kpc, which are amongst the most extreme examples of radio-mode feedback detected to date, implying jet powers of 2-7 x10^45 erg/s. We provide evidence that the AGN inflating these cavities may have only recently transitioned from "quasar-mode" to "radio-mode", and may currently be insufficient to completely offset cooling. A model-subtracted residual X-ray image reveals evidence for prior episodes of strong radio-mode feedback at radii of ~100 kpc, with extended "ghost" cavities indicating a prior epoch of feedback roughly 100 Myr ago. This residual image also exhibits significant asymmetry in the inner ~200 kpc (0.15R500), reminiscent of infalling cool clouds, either due to minor mergers or fragmentation of the cooling ICM. Taken together, these data reveal a rapidly evolving cool core which is rich with structure (both spatially and in temperature), is subject to a variety of highly energetic processes, and yet is cooling rapidly and forming stars along thin, narrow filaments. △ Less

Submitted 24 August, 2015; originally announced August 2015.

Comments: 20 pages, 12 figures. Accepted for publication in ApJ

Stout: Cloudy's Atomic and Molecular Database

Authors: M. L. Lykins, G. J. Ferland, R. Kisielius, M. Chatzikos, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams, F. P. Keenan, P. C. Stancil

Abstract: We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures 2.8 K < T < 1e10 K and densities spanning the low to high-density limits. The radiation field between photon energies $10^{-8}$ Ry and 100 MeV is considered… ▽ More We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures 2.8 K < T < 1e10 K and densities spanning the low to high-density limits. The radiation field between photon energies $10^{-8}$ Ry and 100 MeV is considered, along with all atoms and ions of the lightest 30 elements, and ~100 molecules. For ease of maintenance, the data are stored in a format as close as possible to the original data sources. Few data sources include the full range of data we need. We describe how we fill in the gaps in the data or extrapolate rates beyond their tabulated range. We tabulate data sources both for the atomic spectroscopic parameters and for collision data for the next release of Cloudy. This is not intended as a review of the current status of atomic data, but rather a description of the features of the database which we will build upon. △ Less

Submitted 4 June, 2015; originally announced June 2015.

Comments: ApJ in press

Accurate determination of the free-free Gaunt factor. II - relativistic Gaunt factors

Authors: P. A. M. van Hoof, G. J. Ferland, R. J. R. Williams, K. Volk, M. Chatzikos, M. Lykins, R. L. Porter

Abstract: When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter spac… ▽ More When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter space than has ever been produced before. We present tables of the thermally averaged Gaunt factor covering the range log10(gamma^2) = -6 to 10 and log10(u) = -16 to 13 for all atomic numbers Z = 1 through 36. The data were calculated using the relativistic Bethe-Heitler-Elwert (BHE) approximation and were subsequently merged with accurate non-relativistic results in those parts of the parameter space where the BHE approximation is not valid. These data will be incorporated in the next major release of the spectral synthesis code Cloudy. We also produced tables of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to 10 for all values of Z between 1 and 36. All the data presented in this paper are available online. △ Less

Submitted 24 February, 2015; originally announced February 2015.

Comments: 8 pages, 8 figures, 2 tables

Journal ref: MNRAS, 449, 2112 (2015)

Implications of Coronal Line Emission in NGC 4696

Authors: M. Chatzikos, R. J. R. Williams, G. J. Ferland, R. E. A. Canning, A. C. Fabian, J. S. Sanders, P. A. M. van Hoof, R. M. Johnstone, M. Lykins, R. L. Porter

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

Submitted 16 October, 2014; originally announced October 2014.

Comments: 11 pages, 8 figures, 6 tables; accepted for publication in MNRAS

Accurate determination of the free-free Gaunt factor; I - non-relativistic Gaunt factors

Authors: P. A. M. van Hoof, R. J. R. Williams, K. Volk, M. Chatzikos, G. J. Ferland, M. Lykins, R. L. Porter, Y. Wang

Abstract: Modern spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt fa… ▽ More Modern spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt factors over a much wider range of parameters than has ever been produced before. We first produced a table of non-averaged Gaunt factors, covering the parameter space log10(epsilon_i) = -20 to +10 and log10(w) = -30 to +25. We then continued to produce a table of thermally averaged Gaunt factors covering the parameter space log10(gamma^2) = -6 to +10 and log10(u) = -16 to +13. Finally we produced a table of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to +10. All the data presented in this paper are available online. △ Less

Submitted 17 August, 2014; v1 submitted 18 July, 2014; originally announced July 2014.

Comments: 10 pages, 5 tables, 3 figures. Fixed typo in Eq. 14

Journal ref: MNRAS, 444, 420 (2014)

On the Observability of Optically Thin Coronal Hyperfine Structure Lines

Authors: Marios Chatzikos, Gary Ferland, Robin Williams, Andy Fabian

Abstract: We present Cloudy calculations for the intensity of coronal hyperfine lines in various environments. We model indirect collisional and radiative transitions, and quantify the collisionally-excited line emissivity in the density-temperature phase-space. As an observational aid, we also express the emissivity in units of the continuum in the 0.4--0.7 keV band. For most hyperfine lines, knowledge of… ▽ More We present Cloudy calculations for the intensity of coronal hyperfine lines in various environments. We model indirect collisional and radiative transitions, and quantify the collisionally-excited line emissivity in the density-temperature phase-space. As an observational aid, we also express the emissivity in units of the continuum in the 0.4--0.7 keV band. For most hyperfine lines, knowledge of the X-ray surface brightness and the plasma temperature is sufficient for rough estimates. We find that the radiation fields of both Perseus A and Virgo A can enhance the populations of highly ionized species within 1 kpc. They can also enhance line emissivity within the cluster core. This could have implications for the interpretation of spectra around bright AGN. We find the intensity of the $^{57}$Fe XXIV λ3.068 mm to be about two orders of magnitude fainter than previously thought, at about 20 μK. Comparably bright lines may be found in the infrared. Finally, we find the intensity of hyperfine lines in the Extended Orion Nebula to be low, due to the shallow sightline. Observations of coronal hyperfine lines will likely be feasible with the next generation of radio and sub-mm telescopes. △ Less

Submitted 3 June, 2014; originally announced June 2014.

Comments: 48 pages; 13 figures; 9 tables; additional images available at the online journal

Journal ref: 2014 ApJ, 787, 96

Effects of External Radiation Fields on Line Emission - Application to Star-forming Regions

Authors: Marios Chatzikos, Gary Ferland, Robin Williams, Peter van Hoof, Ryan Porter

Abstract: A variety of astronomical environments contain clouds irradiated by a combination of isotropic and beamed radiation fields. For example, molecular clouds may be irradiated by the isotropic cosmic microwave background (CMB), as well as by a nearby active galactic nucleus (AGN). These radiation fields excite atoms and molecules and produce emission in different ways. We revisit the escape probabilit… ▽ More A variety of astronomical environments contain clouds irradiated by a combination of isotropic and beamed radiation fields. For example, molecular clouds may be irradiated by the isotropic cosmic microwave background (CMB), as well as by a nearby active galactic nucleus (AGN). These radiation fields excite atoms and molecules and produce emission in different ways. We revisit the escape probability theorem and derive a novel expression that accounts for the presence of external radiation fields. We show that when the field is isotropic the escape probability is reduced relative to that in the absence of external radiation. This is in agreement with previous results obtained under ad hoc assumptions or with the two-level system, but can be applied to complex many-level models of atoms or molecules. This treatment is in the development version of the spectral synthesis code Cloudy. We examine the spectrum of a Spitzer cloud embedded in the local interstellar radiation field, and show that about 60 percent of its emission lines are sensitive to background subtraction. We argue that this geometric approach could provide an additional tool toward understanding the complex radiation fields of starburst galaxies. △ Less

Submitted 29 October, 2013; originally announced October 2013.

Comments: 12 pages, 7 figures, accepted for publication to ApJ

Journal ref: 2013, ApJ, 779, 122

Chandra Observation of Abell 2065: An Unequal Mass Merger?

Authors: Marios Chatzikos, Craig L. Sarazin, Joshua C. Kempner

Abstract: We present an analysis of a 41 ks Chandra observation of the merging cluster Abell 2065 with the ACIS-I detector. Previous observations with ROSAT and ASCA provided evidence for an ongoing merger, but also suggested that there were two surviving cooling cores, which were associated with the two cD galaxies in the center of the cluster. The Chandra observation reveals only one X-ray surface brigh… ▽ More We present an analysis of a 41 ks Chandra observation of the merging cluster Abell 2065 with the ACIS-I detector. Previous observations with ROSAT and ASCA provided evidence for an ongoing merger, but also suggested that there were two surviving cooling cores, which were associated with the two cD galaxies in the center of the cluster. The Chandra observation reveals only one X-ray surface brightness peak, which is associated with the more luminous, southern cD galaxy. The gas related with that peak is cool and displaced slightly from the position of the cD. The data suggest that this cool material has formed a cold front. On the other hand, in the higher spatial resolution Chandra image, the second feature to the north is not associated with the northern cD; rather, it appears to be a trail of gas behind the main cD. We argue that only one of the two cooling cores has survived the merger, although it is possible that the northern cD may not have possessed a cool core prior to the merger. We use the cool core survival to constrain the kinematics of the merger and we find an upper limit of ~< 1900 km/s for the merger relative velocity. A surface brightness discontinuity is found at ~140 kpc from the southern cD; the Mach number for this feature is ${\cal M} = 1.66^{+0.24}_{-0.32}$, although its nature (shock or cold front) is not clear from the data. We argue that Abell 2065 is an example of an unequal mass merger. The more massive southern cluster has driven a shock into the ICM of the infalling northern cluster, which has disrupted the cool core of the latter, if one existed originally. We estimate that core crossing occurred a few hundred Myr ago, probably for the first time. △ Less

Submitted 8 February, 2006; originally announced February 2006.

Comments: 15 pages, 10 figures, ApJ in press

Journal ref: Astrophys.J. 643 (2006) 751-763