Showing 1–35 of 35 results for author: Waters, T

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

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

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

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

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

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

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

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

Submitted 24 June, 2024; originally announced June 2024.

A Stellar Dynamical Mass Measurement of the Supermassive Black Hole in NGC 3258

Authors: Thomas K. Waters, Kayhan Gültekin, Karl Gebhardt, Neil Nagar, Vanessa Ávila

Abstract: We present a stellar dynamical mass measurement of the supermassive black hole in the elliptical (E1) galaxy NGC 3258. Our findings are based on Integral Field Unit spectroscopy from the Multi Unit Spectroscopic Explorer observations in Narrow Field Mode with adaptive optics and in the MUSE Wide Field Mode, from which we extract kinematic information by fitting the Ca II and Mg $b$ triplets, respe… ▽ More We present a stellar dynamical mass measurement of the supermassive black hole in the elliptical (E1) galaxy NGC 3258. Our findings are based on Integral Field Unit spectroscopy from the Multi Unit Spectroscopic Explorer observations in Narrow Field Mode with adaptive optics and in the MUSE Wide Field Mode, from which we extract kinematic information by fitting the Ca II and Mg $b$ triplets, respectively. Using axisymmetric, three-integral Schwarzschild orbit-library models, we fit the observed line-of-sight velocity distributions to infer the supermassive black hole mass, the $H$-band mass-to-light ratio, the asymptotic circular velocity, and the dark matter halo scale radius of the galaxy. We report a black hole mass of $(2.2 \pm 0.2)\times10^9 \ \rm M_{\scriptscriptstyle\odot}$ at an assumed distance of $31.9 \ \rm Mpc$. This value is in close agreement with a previous measurement from ALMA CO observations. The consistency between these two measurements provides strong support for both the gas dynamical and stellar dynamical methods. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 18 pages, 12 figures, 1 table

Position dependent radiation fields near accretion disks

Authors: Kara Smith, Daniel Proga, Randall Dannen, Sergei Dyda, Tim Waters

Abstract: In disk wind models for active galactic nuclei (AGN) outflows, high-energy radiation poses a significant problem wherein the gas can become overionized, effectively disabling what is often inferred to be the largest force acting on the gas: the radiation force due to spectral line opacity. Calculations of this radiation force depend on the magnitude of ionizing radiation, which can strongly depend… ▽ More In disk wind models for active galactic nuclei (AGN) outflows, high-energy radiation poses a significant problem wherein the gas can become overionized, effectively disabling what is often inferred to be the largest force acting on the gas: the radiation force due to spectral line opacity. Calculations of this radiation force depend on the magnitude of ionizing radiation, which can strongly depend on the position above a disk where the radiation is anisotropic. As our first step to quantify the position and direction dependence of the radiation field, we assumed free streaming of photons and computed energy distributions of the mean intensity and components of flux as well as energy-integrated quantities such as mean photon energy. We find a significant dependence of radiation field properties on position, but this dependence is not necessarily the same for different field quantities. A key example is that the mean intensity is much softer than the radial flux at many points near the disk. Because the mean intensity largely controls ionization, this softening decreases the severity of the overionization problem. The position dependence of mean intensity implies the position dependence of gas opacity, which we illustrate by computing the radiation force a fluid element feels in an accelerating wind. We find that in a vertical accelerating flow, the force due to radiation is not parallel to the radiation flux. This misalignment is due to the force's geometric weighting by both the velocity field's directionality and the position dependence of the mean intensity. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 14 pages (13 without bibliography), 5 figures, Accepted for publication in ApJ

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

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

Abstract: An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different tim… ▽ More An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response. △ Less

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

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

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

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

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

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

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

Coordinated X-ray and UV absorption within the accretion disk wind of the active galactic nucleus PG 1126-041

Authors: M. Giustini, P. Rodríguez Hidalgo, J. N. Reeves, G. Matzeu, V. Braito, M. Eracleous, G. Chartas, N. Schartel, C. Vignali, P. B. Hall, T. Waters, G. Ponti, D. Proga, M. Dadina, M. Cappi, G. Miniutti, L. de Vries

Abstract: Accretion disk winds launched close to supermassive black holes (SMBHs) are a viable mechanism to provide feedback between the SMBH and the host galaxy. We aim to characterize the X-ray properties of the inner accretion disk wind of the nearby active galactic nucleus (AGN) PG 1126-041, and to study its connection with the ultraviolet (UV)-absorbing wind. We perform spectroscopic analysis of eight… ▽ More Accretion disk winds launched close to supermassive black holes (SMBHs) are a viable mechanism to provide feedback between the SMBH and the host galaxy. We aim to characterize the X-ray properties of the inner accretion disk wind of the nearby active galactic nucleus (AGN) PG 1126-041, and to study its connection with the ultraviolet (UV)-absorbing wind. We perform spectroscopic analysis of eight XMM-Newton observations of PG 1126-041 taken between 2004 and 2015, using both phenomenological models and the most advanced accretion disk wind models available. For half of the dataset, we can compare the X-ray analysis results with the results of quasi-simultaneous, high-resolution spectroscopic UV observations taken with the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope. The X-ray spectra of PG 1126-041 are complex and absorbed by ionized material which is highly variable on multiple time scales, sometimes as short as 11 days. Accretion disk wind models can account for most of the X-ray spectral complexity of PG 1126-041, with the addition of massive clumps, represented by a partially covering absorber. Variations in column density ($N_H \sim 5-20 \times 10^{22}$ cm$^{-2}$) of the partially covering absorber drive the observed X-ray spectral variability of PG 1126-041. The absorption from the X-ray partially covering gas and from the blueshifted C IV troughs appear to vary in a coordinated way. The line of sight toward PG 1126-041 offers a privileged view through a highly dynamic nuclear wind originating on inner accretion disk scales, making the source a very promising candidate for future detailed studies of the physics of accretion disk winds around SMBHs. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: 24 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 679, A73 (2023)

On the transition from efficient to inefficient line-driving in irradiated flows

Authors: Randall Dannen, Daniel Proga, Tim Waters

Abstract: Observations of ionized AGN outflows have provided compelling evidence that the radiation field transfers both momentum and energy to the plasma. At parsec scale distances in AGN, energy transfer can dominate, in which case the only force needed to launch an outflow is that from gas pressure. Much closer to the black hole, gravity dominates thermal energy due to insufficient heating by the radiati… ▽ More Observations of ionized AGN outflows have provided compelling evidence that the radiation field transfers both momentum and energy to the plasma. At parsec scale distances in AGN, energy transfer can dominate, in which case the only force needed to launch an outflow is that from gas pressure. Much closer to the black hole, gravity dominates thermal energy due to insufficient heating by the radiation and the gas is in the so-called `cold' regime. Only magnetic or radiation forces can then lead to outflow, but it is unclear at what temperature and ionization state the radiation force weakens, as these properties depend on the spectral energy distribution (SED). In this work, we survey the parameter space of radiation forces due to spectral lines resulting from blackbody SEDs with varying temperatures in the range $\sim 10^4 - 10^6$~K to identify the radiation temperature at which line-driving begins to lose efficiency. We find that the temperature $\lesssim4\times10^5$~K marks the transition to inefficient line driving. We also self-consistently compute the heating and cooling balance to estimate the gas temperature, so that our parameter survey covers the transition where thermal driving goes from negligible to comparable to line driving. We summarize a large set of hydrodynamical simulations of radial flows to illustrate how the wind properties change during the transition and the dependence of these properties on the assumed SED and governing flow parameters. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 17 pages, 6 figures

The saturation mechanism of thermal instability

Authors: Tim Waters, Daniel Proga

Abstract: The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing recent developments on the linear and nonlinear theory of TI, here we derive several general… ▽ More The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates: whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of the condensation. For isobaric evolution, this 'pressure reversal' occurs nearly simultaneously for every fluid element in the condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave-mediated. We show that both the contraction rebound event and the subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal nonequilibrium. △ Less

Submitted 25 August, 2023; v1 submitted 31 March, 2023; originally announced March 2023.

Comments: Accepted version to appear as a contribution to the research topic "Thermal Imbalance and Multiphase Plasmas Across Scales: From the Solar Corona to the Intracluster Medium" (21 pages, 2 figures)

AGN STORM 2. III. A NICER view of the variable X-ray obscurer in Mrk 817

Authors: Ethan R. Partington, Edward M. Cackett, Erin Kara, Gerard A. Kriss, Aaron J. Barth, Gisella De Rosa, Y. Homayouni, Keith Horne, Hermine Landt, Abderahmen Zoghbi, Rick Edelson, Nahum Arav, Benjamin D. Boizelle, Misty C. Bentz, Michael S. Brotherton, Doyee Byun, Elena Dalla Bonta, Maryam Dehghanian, Pu Du, Carina Fian, Alexei V. Filippenko, Jonathan Gelbord, Michael R. Goad, Diego H. Gonzalez Buitrago, Catherine J. Grier , et al. (22 additional authors not shown)

Abstract: The AGN STORM 2 collaboration targeted the Seyfert 1 galaxy Mrk 817 for a year-long multiwavelength, coordinated reverberation mapping campaign including HST, Swift, XMM-Newton, NICER, and ground-based observatories. Early observations with NICER and XMM revealed an X-ray state ten times fainter than historical observations, consistent with the presence of a new dust-free, ionized obscurer. The fo… ▽ More The AGN STORM 2 collaboration targeted the Seyfert 1 galaxy Mrk 817 for a year-long multiwavelength, coordinated reverberation mapping campaign including HST, Swift, XMM-Newton, NICER, and ground-based observatories. Early observations with NICER and XMM revealed an X-ray state ten times fainter than historical observations, consistent with the presence of a new dust-free, ionized obscurer. The following analysis of NICER spectra attributes variability in the observed X-ray flux to changes in both the column density of the obscurer by at least one order of magnitude ($N_\mathrm{H}$ ranges from $2.85\substack{+0.48\\ -0.33} \times 10^{22}\text{ cm}^{-2}$ to $25.6\substack{+3.0\\ -3.5} \times 10^{22} \text{ cm}^{-2}$) and the intrinsic continuum brightness (the unobscured flux ranges from $10^{-11.8}$ to $10^{-10.5}$ erg s$^{-1}$ cm$^{-2}$ ). While the X-ray flux generally remains in a faint state, there is one large flare during which Mrk 817 returns to its historical mean flux. The obscuring gas is still present at lower column density during the flare but it also becomes highly ionized, increasing its transparency. Correlation between the column density of the X-ray obscurer and the strength of UV broad absorption lines suggests that the X-ray and UV continua are both affected by the same obscuration, consistent with a clumpy disk wind launched from the inner broad line region. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: 19 pages, 9 figures, Accepted to ApJ

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

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

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

Submitted 22 February, 2023; originally announced February 2023.

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

Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching

Authors: Thomas K. Waters, Colton Peterson, Razieh Emami, Xuejian Shen, Lars Hernquist, Randall Smith, Mark Vogelsberger, Charles Alcock, Grant Tremblay, Matthew Liska, John C. Forbes, Jorge Moreno

Abstract: We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner… ▽ More We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback. △ Less

Submitted 5 January, 2024; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 12 figures

Thermal instability in radiation hydrodynamics: instability mechanisms, position-dependent S-curves, and attenuation curves

Authors: Daniel Proga, Tim Waters, Sergei Dyda, Zhaohuan Zhu

Abstract: Local thermal instability can plausibly explain the formation of multiphase gas in many different astrophysical environments, but the theory is only well understood in the optically thin limit of the equations of radiation hydrodynamics (RHD). Here we lay groundwork for transitioning from this limit to a full RHD treatment assuming a gray opacity formalism. We consider a situation where the gas be… ▽ More Local thermal instability can plausibly explain the formation of multiphase gas in many different astrophysical environments, but the theory is only well understood in the optically thin limit of the equations of radiation hydrodynamics (RHD). Here we lay groundwork for transitioning from this limit to a full RHD treatment assuming a gray opacity formalism. We consider a situation where the gas becomes thermally unstable due to the hardening of the radiation field when the main radiative processes are free-free cooling and Compton heating. We identify two ways in which this can happen: (i) when the Compton temperature increases with time, through a rise in either the intensity or energy of a hard X-ray component; and (ii) when attenuation reduces the flux of the thermal component so that the Compton temperature increases with depth through the slab. Both ways likely occur in the broad line region of active galactic nuclei where columns of gas can be ionization bounded. In such instances where attenuation is significant, thermal equilibrium solution curves become position-dependent and it no longer suffices to assess the stability of an irradiated column of gas at all depths using a single equilibrium curve. We demonstrate how to analyze a new equilibrium curve -- the attenuation curve -- for this purpose, and we show that by Field's instability criterion, a negative slope along this curve indicates that constant density slabs are thermally unstable whenever the gas temperature increases with depth. △ Less

Submitted 13 August, 2022; v1 submitted 8 July, 2022; originally announced July 2022.

Comments: 10 pages, 1 figure, version slightly revised and accepted to ApJ Letters

Dynamical thermal instability in highly supersonic outflows

Authors: Tim Waters, Daniel Proga, Randall Dannen, Sergei Dyda

Abstract: Acceleration can change the ionization of X-ray irradiated gas to the point that the gas becomes thermally unstable. Cloud formation, the expected outcome of thermal instability (TI), will be suppressed in a dynamic flow, however, due to the stretching of fluid elements that accompanies acceleration. It is therefore unlikely that cloud formation occurs during the launching phase of a supersonic ou… ▽ More Acceleration can change the ionization of X-ray irradiated gas to the point that the gas becomes thermally unstable. Cloud formation, the expected outcome of thermal instability (TI), will be suppressed in a dynamic flow, however, due to the stretching of fluid elements that accompanies acceleration. It is therefore unlikely that cloud formation occurs during the launching phase of a supersonic outflow. In this paper, we show that the most favorable conditions for dynamical TI in highly supersonic outflows are found at radii beyond the acceleration zone, where the growth rate of entropy modes is set by the linear theory rate for a static plasma. This finding implies that even mildly relativistic outflows can become clumpy, and we explicitly demonstrate this using hydrodynamical simulations of ultrafast outflows. We describe how the continuity and heat equations can be used to appreciate another impediment (beside mode disruption due to the stretching) to making an outflow clumpy: background flow conditions may not allow the plasma to enter a TI zone in the first place. The continuity equation reveals that both impediments are in fact tightly coupled, yet one is easy to overcome. Namely, time variability in the radiation field is found to be a robust means of placing gas in a TI zone. We further show how the ratio of the dynamical and thermal timescales enters linear theory; the heat equation reveals how this ratio depends on the two processes that tend to remove gas from a TI zone -- adiabatic cooling and heat advection. △ Less

Submitted 8 April, 2022; v1 submitted 14 November, 2021; originally announced November 2021.

Comments: Accepted version to appear in ApJ (now 20 pages - Section 2 is somewhat extended/reorganized; Section 4 includes a discussion of the effects of cosmic rays). Simulations viewable at https://trwaters.github.io/dynamicalTI/

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

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

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

Submitted 12 May, 2021; originally announced May 2021.

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

On Synthetic Absorption Line Profiles of Thermally Driven Winds from Active Galactic Nuclei

Authors: Shalini Ganguly, Daniel Proga, Tim Waters, Randall C. Dannen, Sergei Dyda, Margherita Giustini, Timothy Kallman, John Raymond, Jon Miller, Paola Rodriguez Hidalgo

Abstract: The warm absorbers observed in more than half of all nearby active galactic nuclei (AGN) are tracers of ionized outflows located at parsec scale distances from the central engine. If the smallest inferred ionization parameters correspond to plasma at a few $10^4$~K, then the gas undergoes a transition from being bound to unbound provided it is further heated to $\sim 10^6$~K at larger radii. Danne… ▽ More The warm absorbers observed in more than half of all nearby active galactic nuclei (AGN) are tracers of ionized outflows located at parsec scale distances from the central engine. If the smallest inferred ionization parameters correspond to plasma at a few $10^4$~K, then the gas undergoes a transition from being bound to unbound provided it is further heated to $\sim 10^6$~K at larger radii. Dannen et al. recently discovered that under these circumstances, thermally driven wind solutions are unsteady and even show very dense clumps due to thermal instability. To explore the observational consequences of these new wind solutions, we compute line profiles based on the one-dimensional simulations of Dannen et al. We show how the line profiles from even a simple steady state wind solution depend on the ionization energy (IE) of absorbing ions, which is a reflection of the wind ionization stratification. To organize the diversity of the line shapes, we group them into four categories: weak Gaussians, saturated boxy profiles with and without an extended blue wing, and broad weak profiles. The lines with profiles in the last two categories are produced by ions with the highest IE that probe the fastest regions. Their maximum blueshifts agree with the highest flow velocities in thermally unstable models, both steady state and clumpy versions. In contrast, the maximum blueshifts of the most high IE lines in thermally stable models can be less than half of the actual solution velocities. Clumpy solutions can additionally imprint distinguishable absorption troughs at widely separated velocities. △ Less

Submitted 15 May, 2021; v1 submitted 11 March, 2021; originally announced March 2021.

Comments: Accepted for publication in ApJ, 20 pages, 15 figures, 3 tables, videos and FITS files of line profiles available on: http://www.physics.unlv.edu/astro/clumpywindsims-lps.html

Multiphase AGN winds from X-ray irradiated disk atmospheres

Authors: Tim Waters, Daniel Proga, Randall Dannen

Abstract: The mechanism of thermal driving for launching mass outflows is interconnected with classical thermal instability (TI). In a recent paper, we demonstrated that as a result of this interconnectedness, radial wind solutions of X-ray heated flows are prone to becoming clumpy. In this paper, we first show that the Bernoulli function determines whether or not the entropy mode can grow due to TI in dyna… ▽ More The mechanism of thermal driving for launching mass outflows is interconnected with classical thermal instability (TI). In a recent paper, we demonstrated that as a result of this interconnectedness, radial wind solutions of X-ray heated flows are prone to becoming clumpy. In this paper, we first show that the Bernoulli function determines whether or not the entropy mode can grow due to TI in dynamical flows. Based on this finding, we identify a critical `unbound' radius beyond which TI should accompany thermal driving. Our numerical disk wind simulations support this result and reveal that clumpiness is a consequence of buoyancy disrupting the stratified structure of steady state solutions. Namely, instead of a smooth transition layer separating the highly ionized disk wind from the cold phase atmosphere below, hot bubbles formed from TI rise up and fragment the atmosphere. These bubbles first appear within large scale vortices that form below the transition layer, and they result in the episodic production of distinctive cold phase structures referred to as irradiated atmospheric fragments (IAFs). Upon interacting with the wind, IAFs advect outward and develop extended crests. The subsequent disintegration of the IAFs takes place within a turbulent wake that reaches high elevations above the disk. We show that this dynamics has the following observational implications: dips in the absorption measure distribution are no longer expected within TI zones and there can be a less sudden desaturation of X-ray absorption lines such as \OVIII as well as multiple absorption troughs in \FeXXVK. △ Less

Submitted 17 June, 2021; v1 submitted 22 January, 2021; originally announced January 2021.

Comments: Published in ApJ. Simulations viewable at https://trwaters.github.io/multiphaseAGNdiskwinds/ Companion paper arXiv:2103.06497

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

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

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

Submitted 1 October, 2020; originally announced October 2020.

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

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

Clumpy AGN outflows due to thermal instability

Authors: Randall Dannen, Daniel Proga, Tim Waters, Sergei Dyda

Abstract: One of the main mechanisms that could drive mass outflows on parsec scales in AGN is thermal driving. The same X-rays that ionize and heat the plasma are also expected to make it thermally unstable. Indeed, it has been proposed that the observed clumpiness in AGN winds is caused by thermal instability (TI). While many studies employing time-dependent numerical simulations of AGN outflows have incl… ▽ More One of the main mechanisms that could drive mass outflows on parsec scales in AGN is thermal driving. The same X-rays that ionize and heat the plasma are also expected to make it thermally unstable. Indeed, it has been proposed that the observed clumpiness in AGN winds is caused by thermal instability (TI). While many studies employing time-dependent numerical simulations of AGN outflows have included the necessary physics for TI, none have so far managed to produce clumpiness. Here we present the first such clumpy wind simulations in 1-D and 2-D, obtained by simulating parsec scale outflows irradiated by an AGN. By combining an analysis of our extensive parameter survey with physical arguments, we show that the lack of clumps in previous numerical models can be attributed to the following three effects: (i) insufficient radiative heating or other physical processes that prevent the outflowing gas from entering the TI zone; (ii) the stabilizing effect of stretching (due to rapid radial acceleration) in cases where the gas enters the TI zone; and (iii) a flow speed effect: in circumstances where stretching is inefficient, the flow can still be so fast that it passes through the TI zone too quickly for perturbations to grow. Besides these considerations, we also find that a necessary condition to trigger TI in an outflow is for the pressure ionization parameter to decrease along a streamline once gas enters a TI zone. △ Less

Submitted 9 April, 2020; v1 submitted 31 December, 2019; originally announced January 2020.

Comments: 8 pages, 4 figures, 1 table, movies and additional information available on our project webpage: http://www.physics.unlv.edu/astro/clumpywindsims.html

The AGN broad line region as a clumpy turbulent outflow: a physical basis for LOC modeling

Authors: Tim Waters, Hui Li

Abstract: Many studies have considered the roles of clouds, outflows, and turbulence in producing the broad emission and absorption lines in the spectra of active galactic nuclei (AGNs). However, these are often treated as separate or even competing models. Here, we consider the possibility that AGN clouds are condensations formed within the thermally unstable zones of outflows and then compare the typical… ▽ More Many studies have considered the roles of clouds, outflows, and turbulence in producing the broad emission and absorption lines in the spectra of active galactic nuclei (AGNs). However, these are often treated as separate or even competing models. Here, we consider the possibility that AGN clouds are condensations formed within the thermally unstable zones of outflows and then compare the typical sizes of such condensations with the injection scale $k_0^{-1} \sim L_0$ of turbulence, where $L_0$ is assumed to be the scale height of a representative global outflow model. We find that for broad line region (BLR) parameters, clouds are many orders of magnitude smaller than $L_0$ and this has the following implication: BLR cloud dynamics can be modeled using a local approximation through the use of multiphase turbulence simulations of X-ray irradiated plasmas. We present the first such 3D local clumpy turbulent outflow simulations. We show that the condensations share the same type of selection effects characterizing the locally optimally emitting cloud (LOC) scenario, thereby offering a physical interpretation for the LOC model and accounting for its almost uncanny successes. The ubiquitous presence of emission line regions in AGNs can be simply explained as the natural outcome of there being a multiphase interval of $k$-space within the inertial range of a turbulent cascade. △ Less

Submitted 6 December, 2019; originally announced December 2019.

Comments: Submitted to ApJ Letters

Outflows from inflows: the nature of Bondi-like accretion

Authors: Tim Waters, Aycin Aykutalp, Daniel Proga, Jarrett Johnson, Hui Li, Joseph Smidt

Abstract: The classic Bondi solution remains a common starting point both for studying black hole growth across cosmic time in cosmological simulations and for smaller scale simulations of AGN feedback. In nature, however, there will be inhomogenous distributions of rotational velocity and density along the outer radius ($R_o$) marking the sphere of influence of a black hole. While there have been many stud… ▽ More The classic Bondi solution remains a common starting point both for studying black hole growth across cosmic time in cosmological simulations and for smaller scale simulations of AGN feedback. In nature, however, there will be inhomogenous distributions of rotational velocity and density along the outer radius ($R_o$) marking the sphere of influence of a black hole. While there have been many studies of how the Bondi solution changes with a prescribed angular momentum boundary condition, they have all assumed a constant density at $R_o$. In this Letter, we show that a non-uniform density at $R_o$ causes a meridional flow and due to conservation of angular momentum, the Bondi solution qualitatively changes into an inflow-outflow solution. Using physical arguments, we analytically identify the critical logarithmic density gradient $|\partial{\lnρ}/\partialθ|$ above which this change of the solution occurs. For realistic $R_o$, this critical gradient is less than 0.01 and tends to 0 as $R_o \rightarrow \infty$. We show using numerical simulations that, unlike for solutions with an imposed rotational velocity, the accretion rate for solutions under an inhomogenous density boundary condition remains constant at nearly the Bondi rate $\dot{M}_B$, while the outflow rate can greatly exceed $\dot{M}_B$. △ Less

Submitted 21 October, 2019; v1 submitted 2 October, 2019; originally announced October 2019.

Comments: Accepted by MNRAS Letters. Animations viewable at https://trwaters.github.io/OutflowsFromInflows/

Space Telescope and Optical Reverberation Mapping Project. VIII. Time Variability of Emission and Absorption in NGC 5548 Based on Modeling the Ultraviolet Spectrum

Authors: G. A. Kriss, G. De Rosa, J. Ely, B. M. Peterson, J. Kaastra, M. Mehdipour, G. J. Ferland, M. Dehghanian, S. Mathur, R. Edelson, K. T. Korista, N. Arav, A. J. Barth, M. C. Bentz, W. N. Brandt, D. M. Crenshaw, E. Dalla Bontà, K. D. Denney, C. Done, M. Eracleous, M. M. Fausnaugh, E. Gardner, M. R. Goad, C. J. Grier, Keith Horne , et al. (142 additional authors not shown)

Abstract: We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-correcte… ▽ More We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$α$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$α$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state. △ Less

Submitted 12 July, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

Cloud coalescence: a dynamical instability affecting multiphase environments

Authors: Tim Waters, Daniel Proga

Abstract: The mass and size distributions are the key characteristics of any astrophysical objects, including the densest clumps comprising the cold phase of multiphase environments. In our recent papers, we showed how individual clouds of various sizes form and evolve in AGN. In particular, we showed that large clouds undergo damped oscillations as a response to their formation process. Here we followup th… ▽ More The mass and size distributions are the key characteristics of any astrophysical objects, including the densest clumps comprising the cold phase of multiphase environments. In our recent papers, we showed how individual clouds of various sizes form and evolve in AGN. In particular, we showed that large clouds undergo damped oscillations as a response to their formation process. Here we followup this investigation, addressing how different size clouds interact. We find that smaller clouds become trapped in the advective flows generated by larger clouds. The explanation for this behavior leads to a rather remarkable conclusion: even in the absence of gravity, complexes of clouds are dynamically unstable. In an idealized environment (e.g., one free of turbulence and magnetic fields) a perfectly symmetric arrangement of static clouds will remain static, but any small spatial perturbation will lead to all clouds coalescing into a single, large cloud, given enough time. Using numerical simulations, we investigate the main factors that determine the rate of coalescence. Besides the cloud separation distance, we find that the transient response of clouds to a disturbance is the primary factor. Turbulent motions in the flow can easily suppress this tendency for spatially well-separated clouds to coalesce, so it is as yet unclear if this phenomenon can occur in nature. Nevertheless, this work casts strong doubts on a recent hypothesis that large clouds are prone to fragmentation. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: Accepted to ApJ Letters. View simulations at http://www.physics.unlv.edu/astro/WP19bSims.html

Photoionization calculations of the radiation force due to spectral lines in AGNs

Authors: Randall C. Dannen, Daniel Proga, Timothy R. Kallman, Tim Waters

Abstract: One of the main mechanisms that could drive mass outflows in AGNs is radiation pressure due to spectral lines. Although straightforward to understand, the actual magnitude of the radiation force is challenging to compute because the force depends on the physical conditions in the gas, and the strength, spectral energy distribution (SED), and geometry of the radiation field. We present results from… ▽ More One of the main mechanisms that could drive mass outflows in AGNs is radiation pressure due to spectral lines. Although straightforward to understand, the actual magnitude of the radiation force is challenging to compute because the force depends on the physical conditions in the gas, and the strength, spectral energy distribution (SED), and geometry of the radiation field. We present results from our photoionization and radiation transfer calculations of the force multiplier, $M(ξ,t)$, using the same radiation field to compute the gas photoionization and thermal balance. We assume low gas density ($n = 10^4~\rm{cm^{-3}}$) and column density ($N \leq 10^{17}~\rm{ cm^{-2}}$), a Boltzmann distribution for the level populations, and the Sobolev approximation. Here, we describe results for two SEDs corresponding to an unobscured and obscured AGN in NGC 5548. Our main results are the following: 1) although $M(ξ,t)$ starts to decrease with $ξ$ for $ξ\gtrsim 1$ as shown by others, this decrease in our calculations is relatively gradual and could be non-monotonic as $M(ξ,t)$ can increase by a factor of few for $ξ\approx 10-1000$; 2) at these same $ξ$ for which the multiplier is higher than in previous calculations, the gas is thermally unstable by the isobaric criterion; 3) non-LTE effects reduce $M(t,ξ)$ by over two orders of magnitude for $ξ\gtrsim 100$. The dynamical consequence of result (1) is that line driving can be important for $ξ$ as high as $1000$ when the LTE approximation holds, while result (2) provides a natural cloud formation mechanism that may account for the existence of narrow line regions. Result (3) suggests that line driving may not be important for $ξ\gtrsim100$ in tenuous plasma. △ Less

Submitted 19 July, 2019; v1 submitted 4 December, 2018; originally announced December 2018.

Comments: 13 pages, 6 figures, 1 table, project website: http://www.physics.unlv.edu/astro/xstartables.html

Non-isobaric Thermal Instability

Authors: Tim Waters, Daniel Proga

Abstract: Multiphase media have very complex structure and evolution. Accurate numerical simulations are necessary to make advances in our understanding of this rich physics. Because simulations can capture both the linear and nonlinear evolution of perturbations with a relatively wide range of sizes, it is important to thoroughly understand the stability of condensation and acoustic modes between the two e… ▽ More Multiphase media have very complex structure and evolution. Accurate numerical simulations are necessary to make advances in our understanding of this rich physics. Because simulations can capture both the linear and nonlinear evolution of perturbations with a relatively wide range of sizes, it is important to thoroughly understand the stability of condensation and acoustic modes between the two extreme wavelength limits of isobaric and isochoric instability as identified by Field (1965). Partially motivated by a recent suggestion that large non-isobaric clouds can `shatter' into tiny cloudlets, we revisit the linear theory to survey all possible regimes of thermal instability. We uncover seven regimes in total, one of which allows three unstable condensation modes. Using the code Athena++, we determine the numerical requirements to properly evolve small amplitude perturbations of the entropy mode into the nonlinear regime. Our 1D numerical simulations demonstrate that for a typical AGN cooling function, the nonlinear evolution of a single eigenmode in an isobarically unstable plasma involves increasingly larger amplitude oscillations in cloud size, temperature and density as the wavelength increases. Such oscillations are the hallmark behavior of non-isobaric multiphase gas dynamics and may be observable as correlations between changes in brightness and the associated periodic redshifts and blueshifts in systems that can be spatially resolved. Intriguingly, we discuss regimes and derive characteristic cloud sizes for which the saturation process giving rise to these oscillations can be so energetic that the cloud may indeed break apart. However, we dub this process `splattering' instead of `shattering', as it is a different fragmentation mechanism triggered when the cloud suddenly `lands' on the stable cold branch of the equilibrium curve. △ Less

Submitted 8 March, 2019; v1 submitted 26 October, 2018; originally announced October 2018.

Comments: Accepted to ApJ. Simulations can be viewed at http://www.physics.unlv.edu/astro/WP19aSims.html

Magnetothermal disk winds in X-ray binaries: poloidal magnetic fields suppress thermal winds

Authors: Tim Waters, Daniel Proga

Abstract: Magnetic, radiation pressure, and thermal driving are the three mechanisms capable of launching accretion disk winds. In X-ray binaries, radiation pressure is often not significant, as in many systems the luminosity is too low for driving due to continuum transitions yet too high for driving due to line transitions. This leaves thermal and magnetic driving as the contender launching mechanisms in… ▽ More Magnetic, radiation pressure, and thermal driving are the three mechanisms capable of launching accretion disk winds. In X-ray binaries, radiation pressure is often not significant, as in many systems the luminosity is too low for driving due to continuum transitions yet too high for driving due to line transitions. This leaves thermal and magnetic driving as the contender launching mechanisms in these systems. Using Athena++, we perform axisymmetric ideal MHD simulations that include radiative heating and cooling processes appropriate for Compton heated winds to show that the inclusion of magnetic fields into a thermally driven wind has the opposite effect of what one might expect: rather than provide a velocity boost, the thermal wind is suppressed in low plasma beta regions where the field lines are strong enough to reshape the direction of the flow. Our analysis reveals that magneto-centrifugal launching is present but weak, while the reduction in wind velocity is not due to the change in gravitational potential through the magnetically imposed streamline geometry, but rather due to the increased flow tube area just above the surface of the disk, which is less conducive to acceleration. Our results suggest that for magnetothermal wind models to be successful at producing fast dense outflows in low mass X-ray binaries, the winds must be magnetically launched well within the Compton radius. △ Less

Submitted 17 September, 2018; v1 submitted 19 June, 2018; originally announced June 2018.

Comments: 19 pages, updated to match version published in MNRAS

Synthetic absorption lines for a clumpy medium: a spectral signature for cloud acceleration in AGN?

Authors: Tim Waters, Daniel Proga, Randall Dannen, Tim Kallman

Abstract: There is increasing evidence that the highly ionized multiphase components of AGN disk winds may be due to thermal instability. The ions responsible for forming the observed X-ray absorption lines may only exist in relatively cold clumps that can be identified with the so-called 'warm absorbers'. Here we calculate synthetic absorption lines for such warm absorbers from first principles by combinin… ▽ More There is increasing evidence that the highly ionized multiphase components of AGN disk winds may be due to thermal instability. The ions responsible for forming the observed X-ray absorption lines may only exist in relatively cold clumps that can be identified with the so-called 'warm absorbers'. Here we calculate synthetic absorption lines for such warm absorbers from first principles by combining 2D hydrodynamic solutions of a two-phase medium with a dense grid of photoionization models to determine the detailed ionization structure of the gas. Our calculations reveal that cloud disruption, which leads to a highly complicated velocity field (i.e. a clumpy flow), will only mildly affect line shapes and strengths when the cold gas becomes highly mixed but not depleted. Prior to complete disruption, clouds which are optically thin to the driving UV resonance lines will cause absorption at an increasingly blueshifted line of sight velocity as they are accelerated. This behavior will imprint an identifiable signature on the line profile if warm absorbers are enshrouded in an even broader absorption line produced by a high column of intercloud gas. Interestingly, we show that it is possible to develop a spectral diagnostic for cloud acceleration by differencing the absorption components of a doublet line, a result which can be qualitatively understood using a simple partial covering model. Our calculations also permit us to comment on the spectral differences between cloud disruption and ionization changes driven by flux variability. Notably, cloud disruption offers another possibility for explaining absorption line variability. △ Less

Submitted 23 January, 2017; v1 submitted 1 November, 2016; originally announced November 2016.

Comments: 14 pages, 7 figures, to appear in MNRAS

Irradiation of Astrophysical Objects - SED and Flux Effects on Thermally Driven Winds

Authors: Sergei Dyda, Randall Dannen, Tim Waters, Daniel Proga

Abstract: We develop a general method for the self consistent calculation of the hydrodynamics of an astrophysical object irradiated by a radiation field with an arbitrary strength and spectral energy distribution (SED). Using the XSTAR photoionization code, we calculate heating and cooling rates as a function of gas photoionization parameter and temperature for several examples of SEDs: bremsstrahlung, bla… ▽ More We develop a general method for the self consistent calculation of the hydrodynamics of an astrophysical object irradiated by a radiation field with an arbitrary strength and spectral energy distribution (SED). Using the XSTAR photoionization code, we calculate heating and cooling rates as a function of gas photoionization parameter and temperature for several examples of SEDs: bremsstrahlung, blackbody, hard and soft state XRBs, Type 1 and Type 2 AGN. As an application of our method we study the hydrodynamics of 1-dimensional spherical winds heated by a uniform radiation field using the code Athena++. We find that in all cases explored a wind settles into a transonic, steady state. The wind evolves along the radiative heating equilibrium curve until adiabatic cooling effects become important and the flow departs from radiative equilibrium. If the flow is heated very rapidly, for example as in a thermally unstable regime, the corresponding column density of gas is low. Perhaps one of the most intriguing results of our work is the two stage acceleration of the wind that happens when there are two thermally unstable regions and the flux is relatively high. The efficiency with which the radiation field transfers energy to the wind is dependent on the SED of the external source, particularly the relative flux of soft X-rays. These results suggest that detailed photoionization calculations are essential not only to predict spectra but also to properly capture the flow dynamics. △ Less

Submitted 3 February, 2017; v1 submitted 13 October, 2016; originally announced October 2016.

Comments: 14 pages, 8 figures

On the efficient acceleration of clouds in active galactic nuclei

Authors: Tim Waters, Daniel Proga

Abstract: In the broad line region of AGN, acceleration occurs naturally when a cloud condenses out of the hot confining medium due to the increase in line opacity as the cloud cools. However, acceleration by radiation pressure is not very efficient when the flux is time-independent, unless the flow is one-dimensional. Here we explore how acceleration is affected by a time-varying flux, as AGN are known to… ▽ More In the broad line region of AGN, acceleration occurs naturally when a cloud condenses out of the hot confining medium due to the increase in line opacity as the cloud cools. However, acceleration by radiation pressure is not very efficient when the flux is time-independent, unless the flow is one-dimensional. Here we explore how acceleration is affected by a time-varying flux, as AGN are known to be highly variable. If the period of flux oscillations is longer than the thermal timescale, we expect the gas to cool during the low flux state, and therefore line opacity should quickly increase. The cloud will receive a small kick due to the increased radiation force. We perform hydrodynamical simulations using Athena to confirm this effect and quantify its importance. We find that despite the flow becoming turbulent in 2D due to hydrodynamic instabilities, a 20% modulation of the flux leads to a net increase in acceleration --- by more than a factor of 2 --- in both 1D and 2D. We show that this is sufficient to produce the observed line widths, although we only consider optically thin clouds. We discuss the implications of our results for photoionization modeling and reverberation mapping. △ Less

Submitted 6 March, 2016; originally announced March 2016.

Comments: Submitted to MNRAS Letters. Simulations can be viewed online at http://www.physics.unlv.edu/astro/wp16sims.html

Reverberation Mapping of the Broad Line Region: application to a hydrodynamical line-driven disk wind solution

Authors: Tim Waters, Amit Kashi, Daniel Proga, Michael Eracleous, Aaron J. Barth, Jenny Greene

Abstract: The latest analysis efforts in reverberation mapping are beginning to allow reconstruction of echo images (or velocity-delay maps) that encode information about the structure and kinematics of the broad line region (BLR) in active galactic nuclei (AGNs). Such maps can constrain sophisticated physical models for the BLR. The physical picture of the BLR is often theorized to be a photoionized wind l… ▽ More The latest analysis efforts in reverberation mapping are beginning to allow reconstruction of echo images (or velocity-delay maps) that encode information about the structure and kinematics of the broad line region (BLR) in active galactic nuclei (AGNs). Such maps can constrain sophisticated physical models for the BLR. The physical picture of the BLR is often theorized to be a photoionized wind launched from the AGN accretion disk. Previously we showed that the line-driven disk wind solution found in an earlier simulation by Proga and Kallman is virialized over a large distance from the disk. This finding implies that, according to this model, black hole masses can be reliably estimated through reverberation mapping techniques. However, predictions of echo images expected from line-driven disk winds are not available. Here, after presenting the necessary radiative transfer methodology, we carry out the first calculations of such predictions. We find that the echo images are quite similar to other virialized BLR models such as randomly orbiting clouds and thin Keplerian disks. We conduct a parameter survey exploring how echo images, line profiles, and transfer functions depend on both the inclination angle and the line opacity. We find that the line profiles are almost always single peaked, while transfer functions tend to have tails extending to large time delays. The outflow, despite being primarily equatorially directed, causes an appreciable blue-shifted excess on both the echo image and line profile when seen from lower inclinations ($i \lesssim 45^\circ$). This effect may be observable in low ionization lines such as $\rm{H}β$. △ Less

Submitted 23 April, 2016; v1 submitted 20 January, 2016; originally announced January 2016.

Comments: Accepted version (ApJ): added Sec 4.5; improved Sec 2.2

Cloud Formation and Acceleration in a Radiative Environment

Authors: Daniel Proga, Tim Waters

Abstract: In a radiatively heated and cooled medium, the thermal instability is a plausible mechanism for forming clouds, while the radiation force provides a natural acceleration, especially when ions recombine and opacity increases. Here we extend Field's theory to self-consistently account for a radiation force resulting from bound-free and bound-bound transitions in the optically thin limit. We present… ▽ More In a radiatively heated and cooled medium, the thermal instability is a plausible mechanism for forming clouds, while the radiation force provides a natural acceleration, especially when ions recombine and opacity increases. Here we extend Field's theory to self-consistently account for a radiation force resulting from bound-free and bound-bound transitions in the optically thin limit. We present physical arguments for clouds to be significantly accelerated by a radiation force due to lines during a nonlinear phase of the instability. To qualitatively illustrate our main points, we perform both one and two-dimensional (1-D/2-D) hydrodynamical simulations that allow us to study the nonlinear outcome of the evolution of thermally unstable gas subjected to this radiation force. Our 1-D simulations demonstrate that the thermal instability can produce long-lived clouds that reach a thermal equilibrium between radiative processes and thermal conduction, while the radiation force can indeed accelerate the clouds to supersonic velocities. However, our 2-D simulations reveal that a single cloud with a simple morphology cannot be maintained due to destructive processes, triggered by the Rayleigh-Taylor instability and followed by the Kelvin-Helmholtz instability. Nevertheless, the resulting cold gas structures are still significantly accelerated before they are ultimately dispersed. △ Less

Submitted 10 February, 2015; originally announced February 2015.

Comments: 10 pages, 7 figures. Resubmitted to ApJ after addressing referee's comments. Simulations viewable at http://www.physics.unlv.edu/astro/pw15sims.html

Evolution of the L1 halo family in the radial solar sail CRTBP

Authors: Patricia Verrier, Thomas Waters, Jan Sieber

Abstract: We present a detailed investigation of the dramatic changes that occur in the $\mathcal{L}_1$ halo family when radiation pressure is introduced into the Sun-Earth circular restricted three-body problem (CRTBP). This photo-gravitational CRTBP can be used to model the motion of a solar sail orientated perpendicular to the Sun-line. The problem is then parameterized by the sail lightness number, the… ▽ More We present a detailed investigation of the dramatic changes that occur in the $\mathcal{L}_1$ halo family when radiation pressure is introduced into the Sun-Earth circular restricted three-body problem (CRTBP). This photo-gravitational CRTBP can be used to model the motion of a solar sail orientated perpendicular to the Sun-line. The problem is then parameterized by the sail lightness number, the ratio of solar radiation pressure acceleration to solar gravitational acceleration. Using boundary-value problem numerical continuation methods and the AUTO software package (Doedel et al. 1991) the families can be fully mapped out as the parameter $β$ is increased. Interestingly, the emergence of a branch point in the retrograde satellite family around the Earth at $β\approx0.0387$ acts to split the halo family into two new families. As radiation pressure is further increased one of these new families subsequently merges with another non-planar family at $β\approx0.289$, resulting in a third new family. The linear stability of the families changes rapidly at low values of $β$, with several small regions of neutral stability appearing and disappearing. By using existing methods within AUTO to continue branch points and period-doubling bifurcations, and deriving a new boundary-value problem formulation to continue the folds and Krein collisions, we track bifurcations and changes in the linear stability of the families in the parameter $β$ and provide a comprehensive overview of the halo family in the presence of radiation pressure. The results demonstrate that even at small values of $β$ there is significant difference to the classical CRTBP, providing opportunity for novel solar sail trajectories. Further, we also find that the branch points between families in the solar sail CRTBP provide a simple means of generating certain families in the classical case. △ Less

Submitted 28 July, 2014; v1 submitted 15 August, 2013; originally announced August 2013.

Comments: 31 pages, 17 figures, accepted by Celestial Mechanics and Dynamical Astronomy

Parker Winds Revisited: An Extension to Disc Winds

Authors: Timothy R. Waters, Daniel Proga

Abstract: A simple 1D dynamical model of thermally driven disc winds is proposed, based on the results of recent, 2.5D axi-symmetric simulations. Our formulation of the disc wind problem is in the spirit of the original Parker (1958) and Bondi (1952) problems, namely we assume an elementary flow configuration consisting of an outflow following pre-defined trajectories in the presence of a central gravitatin… ▽ More A simple 1D dynamical model of thermally driven disc winds is proposed, based on the results of recent, 2.5D axi-symmetric simulations. Our formulation of the disc wind problem is in the spirit of the original Parker (1958) and Bondi (1952) problems, namely we assume an elementary flow configuration consisting of an outflow following pre-defined trajectories in the presence of a central gravitating point mass. Viscosity and heat conduction are neglected. We consider two different streamline geometries, both comprised of straight lines in the (x,z)-plane: (i) streamlines that converge to a geometric point located at (x,z)=(0,-d) and (ii) streamlines that emerge at a constant inclination angle from the disc midplane (the x-axis, as we consider geometrically thin accretion discs). The former geometry is commonly used in kinematic models to compute synthetic spectra, while the latter, which exhibits self-similarity, is likely unused for this purpose, although it easily can be with existing kinematic models. We make the case that it should be, i.e. that geometry (ii) leads to transonic wind solutions with substantially different properties owing to its lack of streamline divergence. Both geometries can be used to complement recent efforts to estimate photoevaporative mass loss rates from protoplanetary discs. Pertinent to understanding our disc wind results, which are also applicable to X-ray binaries and active galactic nuclei, is a focused discussion on lesser known properties of classic Parker wind solutions. We find that the parameter space corresponding to decelerating Parker wind solutions is made larger due to rotation and leads instead to disc wind solutions that always accelerate after the bulk velocity is slowed to a minimum value. Surprisingly, Keplerian rotation may allow for two different transonic wind solutions for the same physical conditions. △ Less

Submitted 31 July, 2012; originally announced July 2012.

Comments: 30 pages, 13 figures, accepted for publication in MNRAS

Parker Winds Revisited: An Extension to Disk Winds

Authors: Tim Waters, Daniel Proga

Abstract: A simple, one-dimensional dynamical model of thermally driven disk winds, one in the spirit of the original Parker (1958) model, is presented. We consider two different axi-symmetric streamline geometries: geometry (i) is commonly used in kinematic models to compute synthetic spectra, while geometry (ii), which exhibits self-similarity and more closely resembles the geometry found by many numerica… ▽ More A simple, one-dimensional dynamical model of thermally driven disk winds, one in the spirit of the original Parker (1958) model, is presented. We consider two different axi-symmetric streamline geometries: geometry (i) is commonly used in kinematic models to compute synthetic spectra, while geometry (ii), which exhibits self-similarity and more closely resembles the geometry found by many numerical simulations of disk winds, is likely unused for this purpose - although it easily can be with existing kinematic models. We make the case that it should be, i.e. that geometry (ii) leads to transonic wind solutions with substantially different properties. △ Less

Submitted 4 January, 2012; originally announced January 2012.

Comments: 2 pages, 1 figure. Conference proceedings, AGN Winds in Charleston