Francisco Guzmán

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.

In this paper, we discuss atomic processes modifying the soft X-ray spectra 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 _mod).” If 0 ≤ f _mod ≤ 1, the line is suppressed, which could be the case in both collisionally ionized and photoionized systems. If f _mod ≥ 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 Medium Energy Grating (MEG) spectrum of V1223 Sgr with a combination of 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.

Future microcalorimeter X-ray observations will resolve spectral features in unmatched detail. Understanding of line formation processes in 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 the 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.

Understanding the Fe II emission from active galactic nuclei (AGNs) 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 data sets with radiative and electron collisional rates have become available. We have incorporated these into the spectral synthesis code Cloudy and examined predictions using a new generation of AGN spectral energy distribution (SED), which indicates that the ultraviolet (UV) emission can be quite different depending on the data set utilized. The Smyth et al. data set 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 ≈100 km s⁻¹ reproduces the observed shape and strength of the so-called Fe II “UV bump.” Comparing our predictions to the observed Fe II template, we derive a typical cloud density of 10¹¹ cm⁻³ and photon flux of 10²⁰ cm⁻² s⁻¹, 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)) ∼ 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.

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 that 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, in order 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 He ii or FeKα. We propose that this is the situation in most of the AGNs. Case 2 winds, which 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 a H⁺ ionization front and block much of the Lyman continuum. The results show that the presence of these 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.

The Soft X-ray Spectrometer on board Hitomi, with the unprecedented resolving power of R ∼ 1250, allowed the detection of members of the Fe xxv Kα 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 Hitomi observations in the outer region of the Perseus core using collision strengths from different atomic data sets, and obtain good agreement. We also show the effect of turbulence on Fe xxv Kα line ratios and interplay between column density and metallicity. Additionally, we discuss the atomic number dependence of transition probabilities for allowed and unallowed transitions, which causes highly charged He-like systems, such as Fe xxv, to behave fundamentally differently from He i.

We discuss the importance of Fe²³⁺ in determining the line intensities of the Fe xxv Kα complex in an optically thick cloud, and investigate resonant 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 are absorbed by Fe²³⁺ and destroyed, causing the x line intensity to decrease. For example, at a hydrogen column density of 10²⁵ cm⁻², ∼32% of x photons are destroyed due to RAD while w is mostly unaffected. The line intensity of y is slightly (≤2%) reduced. z is not directly affected by RAD, but the contrasting behavior between z and x line intensities points toward the possible conversion of a tiny fraction (∼2%) of x photons into z photons. The change in line intensities due to electron scattering escape off fast thermal electrons is also discussed.

In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper showed that changes in the LOS obscurer produces the decorrelation between the absorption lines and the continuum. A second study showed that the base of the wind shields the broad emission-line region (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 spectral energy distribution (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 Hubble Space Telescope 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.

The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is one of the most intensive multiwavelength AGN monitoring campaigns 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,” suggesting that unobserved changes to the ionizing continuum were present. To understand this remarkable phenomenon and to obtain an independent assessment of the ionizing continuum variations, 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 this cycle, variations in this obscurer’s line-of-sight covering factor modify the soft X-ray continuum, changing the ionization of helium. Ionizing radiation produced by recombining helium then affects the level of ionization of some ions seen by the Hubble Space Telescope. In particular, high-ionization species are affected by changes in the obscurer covering factor, which does not affect the optical or UV continuum, and thus appear as uncorrelated changes, a “holiday.” It is likely that any other model that selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission-line behavior observed.