Showing 1–50 of 150 results for author: Proga, D

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.

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

Time-Dependent AGN Disc Winds I -- X-ray Irradiation

Authors: Sergei Dyda, Shane W. Davis, Daniel Proga

Abstract: We study AGN line driven disc winds using time-dependent radiation hydrodynamics. The key criterion for determining wind launching is the coupling strength of the UV radiation field via the spectral lines of the gas. The strength of these lines in turn relies crucially on the gas ionization state, determined by the local X-ray intensity. We consider a suite of models where the central ionizing rad… ▽ More We study AGN line driven disc winds using time-dependent radiation hydrodynamics. The key criterion for determining wind launching is the coupling strength of the UV radiation field via the spectral lines of the gas. The strength of these lines in turn relies crucially on the gas ionization state, determined by the local X-ray intensity. We consider a suite of models where the central ionizing radiation is affected by scattering, absorption and re-emission by the intervening gas. In a pure attenuation model, the disc launches an episodic wind, as previous studies have shown. Including scattering or re-emission tends to weaken the wind, lowering the mass flux and outflow velocity and if sufficiently dominant, suppressing the outflow entirely. However, the exponential nature of radiative attenuation means only a modest, factor of a few, increase in the absorption cross section can overcome the wind suppression due to scattering and re-emission. We find mass outflow rates of $\sim 20\%$ or more of the assumed inflow rate through the disk, indicating that radiation driven winds may significantly alter the structure of the accretion flow. The winds also supply a large, time-varying column of material above the nominal constant disk scale height, which will determine the geometry of reprocessed emission from the central source. Our results suggest the need for accurate photoionization modeling, radiation transport as well as accretion disc physics, to study their effects on the AGN disc winds △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 13 pages, 4 figures

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

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

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

Submitted 2 October, 2023; originally announced October 2023.

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

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)

The Black Hole Candidate Swift J1728.9$-$3613 and the Supernova Remnant G351.9$-$0.9

Authors: Mayura Balakrishnan, Paul A. Draghis, Jon M. Miller, Joe Bright, Robert Fender, Mason Ng, Edward Cackett, Andrew Fabian, Kip Kuntz, James C. A. Miller-Jones, Daniel Proga, Paul S. Ray, John Raymond, Mark Reynolds, Abderahmen Zoghbi

Abstract: A number of neutron stars have been observed within the remnants of the core-collapse supernova explosions that created them. In contrast, black holes are not yet clearly associated with supernova remnants. Indeed, some observations suggest that black holes are ``born in the dark'', i.e. without a supernova explosion. Herein, we present a multi-wavelength analysis of the X-ray transient Swift J172… ▽ More A number of neutron stars have been observed within the remnants of the core-collapse supernova explosions that created them. In contrast, black holes are not yet clearly associated with supernova remnants. Indeed, some observations suggest that black holes are ``born in the dark'', i.e. without a supernova explosion. Herein, we present a multi-wavelength analysis of the X-ray transient Swift J1728.9$-$3613, based on observations made with Chandra, ESO-VISTA, MeerKAT, NICER, NuSTAR, Swift, and XMM-Newton. Three independent diagnostics indicate that the system likely harbors a black hole primary. Infrared imaging signals a massive companion star that is broadly consistent with an A or B spectral type. Most importantly, the X-ray binary lies within the central region of the catalogued supernova remnant G351.9$-$0.9. Our deep MeerKAT image at 1.28~GHz signals that the remnant is in the Sedov phase; this fact and the non-detection of the soft X-ray emission expected from such a remnant argue that it lies at a distance that could coincide with the black hole. Utilizing a formal measurement of the distance to Swift J1728.9$-$3613 ($d = 8.4\pm 0.8$ kpc), a lower limit on the distance to G351.9$-$0.9 ($d \geq 7.5$ kpc), and the number and distribution of black holes and supernova remnants within the Milky Way, extensive simulations suggest that the probability of a chance superposition is $<1.7\%$ ($99.7\%$ credible interval). The discovery of a black hole within a supernova remnant would support numerical simulations that produce black holes and remnants, and thus provide clear observational evidence of distinct black hole formation channels. We discuss the robustness of our analysis and some challenges to this interpretation. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 25 pages, 10 figures, accepted for publication to ApJ

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

Supermassive Black Hole Winds in X-rays -- SUBWAYS. I. Ultra-fast outflows in QSOs beyond the local Universe

Authors: G. A. Matzeu, M. Brusa, G. Lanzuisi, M. Dadina, S. Bianchi, G. Kriss, M. Mehdipour, E. Nardini, G. Chartas, R. Middei, E. Piconcelli, V. Gianolli, A. Comastri, A. L. Longinotti, Y. Krongold, F. Ricci, P. O. Petrucci, F. Tombesi, A. Luminari, L. Zappacosta, G. Miniutti, M. Gaspari, E. Behar, M. Bischetti, S. Mathur , et al. (26 additional authors not shown)

Abstract: We present a new X-ray spectroscopic study of $22$ luminous ($2\times10^{45}\lesssim L_{\rm bol}\rm /erg\,s^{-1} \lesssim 2\times10^{46}$) active galactic nuclei (AGNs) at intermediate-redshift ($0.1 \lesssim z \lesssim 0.4$), as part of the SUpermassive Black hole Winds in the x-rAYS (SUBWAYS) sample, mostly composed of quasars (QSOs) and type\,1 AGN. Here, 17 targets were observed with \textit{X… ▽ More We present a new X-ray spectroscopic study of $22$ luminous ($2\times10^{45}\lesssim L_{\rm bol}\rm /erg\,s^{-1} \lesssim 2\times10^{46}$) active galactic nuclei (AGNs) at intermediate-redshift ($0.1 \lesssim z \lesssim 0.4$), as part of the SUpermassive Black hole Winds in the x-rAYS (SUBWAYS) sample, mostly composed of quasars (QSOs) and type\,1 AGN. Here, 17 targets were observed with \textit{XMM-Newton} between 2019--2020 and the remaining 5 are from previous observations. The aim of this large campaign ($1.45\,\rm Ms$ duration) is to characterise the various manifestations of winds in the X-rays driven from supermassive black holes in AGN. In this paper we focus on the search and characterization of ultra-fast outflows (UFOs), which are typically detected through blueshifted absorption troughs in the Fe\,K band ($E>7\,\rm keV$). By following Monte Carlo procedures, we confirm the detection of absorption lines corresponding to highly ionised iron (e.g., Fe\,\textsc{xxv}\,H$α$, Fe\,\textsc{xxvi}\,Ly$α$) in 7/22 sources at the $\gtrsim95\%$ confidence level (for each individual line). The global combined probability of such absorption features in the sample is $>99.9\%$. The SUBWAYS campaign extends at higher luminosity and redshifts than previous local studies on Seyferts, obtained using \xmm and \suzaku observations. We find a UFO detection fraction of $\sim30\%$ on the total sample that is in agreement with the previous findings. This work independently provides further support for the existence of highly-ionised matter propagating at mildly relativistic speed ($\gtrsim0.1c$) in a considerable fraction of AGN over a broad range of luminosities, which is expected to play a key role in the self-regulated AGN feeding-feedback cycle, as also supported by hydrodynamical multiphase simulations. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 30 pages, 10 figures, accepted for publication in A&A

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/

A Spectroscopic Angle on Central Engine Size Scales in Accreting Neutron Stars

Authors: Nicolas Trueba, J. M. Miller, A. C. Fabian, J. Kaastra, T. Kallman, A. Lohfink, R. M. Ludlam, D. Proga, J. Raymond, C. Reynolds, M. Reynolds, A. Zoghbi

Abstract: Analyses of absorption from disk winds and atmospheres in accreting compact objects typically treat the central emitting regions in these systems as point sources relative to the absorber. This assumption breaks down if the absorbing gas is located within $few \times 1000\cdot GM/{c}^{2}$, in which case a small component of the absorber's Keplerian motion contributes to the velocity-width of absor… ▽ More Analyses of absorption from disk winds and atmospheres in accreting compact objects typically treat the central emitting regions in these systems as point sources relative to the absorber. This assumption breaks down if the absorbing gas is located within $few \times 1000\cdot GM/{c}^{2}$, in which case a small component of the absorber's Keplerian motion contributes to the velocity-width of absorption lines. Here, we demonstrate how this velocity-broadening effect can be used to constrain the sizes of central engines in accreting compact objects via a simple geometric relationship, and develop a method for modeling this effect. We apply this method on the Chandra/HETG spectra of three ultra-compact and short period neutron star X-ray binaries in which evidence of gravitationally redshifted absorption, owing to an inner-disk atmosphere, has recently been reported. The significance of the redshift is above $5σ$ for XTE J1710$-$281 (this work) and 4U 1916$-$053, and is inconsistent with various estimates of the relative radial velocity of each binary. For our most sensitive spectrum (XTE J1710$-$281), we obtain a 1$σ$ upper bound of 310 $\text{km}$ $\text{s}^{-1}$ on the magnitude of this geometric effect and a central engine of size ${R}_{CE} < 60 ~ GM/{c}^{2}$ (or, $< 90 ~ GM/{c}^{2}$ at the $3σ$ level). These initial constraints compare favorably to those obtained via microlensing in quasars and approach the sensitivity of constraints via relativistic reflection in neutron stars. This sensitivity will increase with further exposures, as well as the launch of future microcalorimeter and grating missions. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: Accepted to ApJ, 32 Pages (3 Appendices), 17 Figures

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

Ionized outflows from active galactic nuclei as the essential elements of feedback

Authors: Sibasish Laha, Christopher S. Reynolds, James Reeves, Gerard Kriss, Matteo Guainazzi, Randall Smith, Sylvain Veilleux, Daniel Proga

Abstract: Outflows from active galactic nuclei (AGN) are one of the fundamental mechanisms by which the central supermassive black hole interacts with its host galaxy. Detected in $\ge 50\%$ of nearby AGN, these outflows have been found to carry kinetic energy that is a significant fraction of AGN power, and thereby give negative feedback to their host galaxies. To understand the physical processes that reg… ▽ More Outflows from active galactic nuclei (AGN) are one of the fundamental mechanisms by which the central supermassive black hole interacts with its host galaxy. Detected in $\ge 50\%$ of nearby AGN, these outflows have been found to carry kinetic energy that is a significant fraction of AGN power, and thereby give negative feedback to their host galaxies. To understand the physical processes that regulate them, it is important to have a robust estimate of their physical and dynamical parameters. In this review we summarize our current understanding on the physics of the ionized outflows detected in absorption in the UV and X-ray wavelength bands. We discuss the most relevant observations and our current knowledge and uncertainties in the measurements of the outflow parameters. We also discuss their origin and acceleration mechanisms. The commissioning and concept studies of large telescope missions with high resolution spectrographs in UV/optical and X-rays along with rapid advancements in simulations offer great promise for discoveries in this field over the next decade. △ Less

Submitted 12 December, 2020; originally announced December 2020.

Comments: An invited Review Article published in Nature Astronomy: https://www.nature.com/articles/s41550-020-01255-2

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

A Redshifted Inner Disk Atmosphere and Transient Absorbers in the Ultra-Compact Neutron Star X-ray Binary 4U 1916-053

Authors: Nicolas Trueba, J. M. Miller, A. C. Fabian, J. Kaastra, T. Kallman, A. Lohfink, D. Proga, J. Raymond, C. Reynolds, M. Reynolds, A. Zoghbi

Abstract: The very small accretion disks in ultra-compact X-ray binaries (UCXBs) are special laboratories in which to study disk accretion and outflows. We report on three sets of new (250 ks total) and archival (50 ks) Chandra/HETG observations of the "dipping" neutron-star X-ray binary 4U 1916$-$053, which has an orbital period of $P\simeq 50$~minutes. We find that the bulk of the absorption in all three… ▽ More The very small accretion disks in ultra-compact X-ray binaries (UCXBs) are special laboratories in which to study disk accretion and outflows. We report on three sets of new (250 ks total) and archival (50 ks) Chandra/HETG observations of the "dipping" neutron-star X-ray binary 4U 1916$-$053, which has an orbital period of $P\simeq 50$~minutes. We find that the bulk of the absorption in all three spectra originates in a disk atmosphere that is redshifted by $v\simeq 220-290$ $\text{km}$ $\text{s}^{-1}$, corresponding to the gravitational redshift at radius of $R \sim 1200$ $GM/{c}^{2}$. This shift is present in the strongest, most highly ionized lines (Si XIV and Fe XXVI), with a significance of 5$σ$. Absorption lines observed during dipping events (typically associated with the outermost disk) instead display no velocity shifts and serve as a local standard of rest, suggesting that the redshift is intrinsic to an inner disk atmosphere and not due to radial motion in the galaxy or a kick. In two spectra, there is also evidence of a more strongly redshifted component that would correspond to a disk atmosphere at $R \sim 70$ $GM/{c}^{2}$; this component is significant at the 3$σ$ level. Finally, in one spectrum, we find evidence of disk wind with a blue shift of $v = {-1700}^{+1700}_{-1200}$ $\text{km}$ $\text{s}^{-1}$. If real, this wind would require magnetic driving. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 15 Pages, 6 Figures, Accepted to ApJ Letters

An Obscured, Seyfert-2-like State of the Stellar-mass Black Hole GRS 1915+105 Caused by Failed Disk Winds

Authors: J. M. Miller, A. Zoghbi, J. Raymond, M. Balakrishnan, L. Brenneman, E. Cackett, P. Draghis, A. C. Fabian, E. Gallo, J. Kaastra, T. Kallman, E. Kammoun, S. E. Motta, D. Proga, M. T. Reynolds, N. Trueba

Abstract: We report on Chandra gratings spectra of the stellar-mass black hole GRS 1915+105 obtained during a novel, highly obscured state. As the source entered this state, a dense, massive accretion disk wind was detected through strong absorption lines. Photionization modeling indicates that it must originate close to the central engine, orders of magnitude from the outer accretion disk. Strong, nearly s… ▽ More We report on Chandra gratings spectra of the stellar-mass black hole GRS 1915+105 obtained during a novel, highly obscured state. As the source entered this state, a dense, massive accretion disk wind was detected through strong absorption lines. Photionization modeling indicates that it must originate close to the central engine, orders of magnitude from the outer accretion disk. Strong, nearly sinusoidal flux variability in this phase yielded a key insight: the wind is blue-shifted when its column density is relatively low, but red-shifted as it approaches the Compton-thick threshold. At no point does the wind appear to achieve the local escape velocity; in this sense, it is a "failed wind." Later observations suggest that the disk ultimately fails to keep even the central engine clear of gas, leading to heavily obscured and Compton-thick states characterized by very strong Fe K emission lines. Indeed, these later spectra are successfully described using models developed for obscured AGN. We discuss our results in terms the remarkable similarity of GRS 1915+105 deep in its "obscured state" to Seyfert-2 and Compton-thick AGN, and explore how our understanding of accretion and obscuration in massive black holes is impacted by our observations. △ Less

Submitted 14 July, 2020; originally announced July 2020.

Comments: Submitted to ApJ

A hard look at local, optically-selected, obscured Seyfert galaxies

Authors: E. S. Kammoun, J. M. Miller, M. Koss, K. Oh, A. Zoghbi, R. F. Mushotzky, D. Barret, E. Behar, W. N. Brandt, L. W. Brenneman, J. S. Kaastra, A. M. Lohfink, D. Proga, D. Stern

Abstract: We study the X-ray spectra of a sample of 19 obscured, optically-selected Seyfert galaxies (Sy 1.8, 1.9 and 2) in the local universe ($d \leq 175$~Mpc), drawn from the CfA Seyfert sample. Our analysis is driven by the high sensitivity of NuSTAR in the hard X-rays, coupled with soft X-ray spectra using XMM-Newton, Chandra, Suzaku, and Swift/XRT. We also analyze the optical spectra of these sources… ▽ More We study the X-ray spectra of a sample of 19 obscured, optically-selected Seyfert galaxies (Sy 1.8, 1.9 and 2) in the local universe ($d \leq 175$~Mpc), drawn from the CfA Seyfert sample. Our analysis is driven by the high sensitivity of NuSTAR in the hard X-rays, coupled with soft X-ray spectra using XMM-Newton, Chandra, Suzaku, and Swift/XRT. We also analyze the optical spectra of these sources in order to obtain accurate mass estimates and Eddington fractions. We employ four different models to analyze the X-ray spectra of these sources, which all result in consistent results. We find that 79-90 % of the sources are heavily obscured with line-of-sight column density $N_{\rm H} > 10^{23}~\rm cm^{-2}$. We also find a Compton-thick ($N_{\rm H} > 10^{24}~\rm cm^{-2}$) fraction of $37-53$ %. These results are consistent with previous estimates based on multi-wavelength analyses. We find that the fraction of reprocessed to intrinsic emission is positively correlated with $N_{\rm H}$ and negatively correlated with the intrinsic, unabsorbed, X-ray luminosity (in agreement with the Iwasawa-Taniguchi effect). Our results support the hypothesis that radiation pressure regulates the distribution of the circumnuclear material. △ Less

Submitted 24 August, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

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

Swift Spectroscopy of the Accretion Disk Wind in the Black Hole GRO J1655-40

Authors: M. Balakrishnan, J. M. Miller, N. Trueba, M. Reynolds, J. Raymond, D. Proga, A. C. Fabian, T. Kallman, J. Kaastra

Abstract: Chandra obtained two High Energy Transmission Grating (HETG) spectra of the stellar-mass black hole GRO J1655-40 during its 2005 outburst, revealing a rich and complex disk wind. Soon after its launch, the Neil Gehrels Swift Observatory began monitoring the same outburst. Some X-ray Telescope (XRT) observations were obtained in a mode that makes it impossible to remove strong Mn calibration lines,… ▽ More Chandra obtained two High Energy Transmission Grating (HETG) spectra of the stellar-mass black hole GRO J1655-40 during its 2005 outburst, revealing a rich and complex disk wind. Soon after its launch, the Neil Gehrels Swift Observatory began monitoring the same outburst. Some X-ray Telescope (XRT) observations were obtained in a mode that makes it impossible to remove strong Mn calibration lines, so the Fe K-alpha line region in the spectra was previously neglected. However, these lines enable a precise calibration of the energy scale, facilitating studies of the absorption-dominated disk wind and its velocity shifts. Here, we present fits to 15 Swift/XRT spectra, revealing variability and evolution in the outflow. The data strongly point to a magnetically driven disk wind: both the higher velocity (e.g., v ~ 10^4 km/s) and lower velocity (e.g., v ~ 10^3 km/s) wind components are typically much faster than is possible for thermally driven outflows (v < 200 km/s), and photoionization modeling yields absorption radii that are two orders of magnitude below the Compton radius that defines the typical inner extent of thermal winds. Moreover, correlations between key wind parameters yield an average absorption measure distribution (AMD) that is consistent with magnetohydrodynamic wind models. We discuss our results in terms of recent observational and theoretical studies of black hole accretion disks and outflows, and future prospects. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: Accepted for publication in ApJ

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

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

Abstract: We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10… ▽ More We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon. △ Less

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

Comments: 19 pages, 9 figures, ApJ in press

The role of failed accretion disk winds in active galactic nuclei

Authors: Margherita Giustini, Daniel Proga

Abstract: Both observational and theoretical evidence point at outflows originating from accretion disks as fundamental ingredients of active galactic nuclei (AGN). These outflows can have more than one component, for example an unbound supersonic wind and a failed wind (FW). The latter is a prediction of the simulations of radiation-driven disk outflows which show that the former is accompanied by an inner… ▽ More Both observational and theoretical evidence point at outflows originating from accretion disks as fundamental ingredients of active galactic nuclei (AGN). These outflows can have more than one component, for example an unbound supersonic wind and a failed wind (FW). The latter is a prediction of the simulations of radiation-driven disk outflows which show that the former is accompanied by an inner failed component, where the flow struggles to escape from the strong gravitational pull of the supermassive black hole. This FW component could provide a physical framework to interpret various phenomenological components of AGN. Here we briefly discuss a few of them: the broad line region, the X-ray obscurer, and the X-ray corona. △ Less

Submitted 19 February, 2020; v1 submitted 18 February, 2020; originally announced February 2020.

Comments: 5 pages no figures. Proceedings of the IAU Symposium No. 356, "Nuclear Activity in Galaxies Across Cosmic Time" (Addis Ababa, 7-11 Oct 2019)

Journal ref: Proc. IAU 15 (2019) 82-86

On the wind-driven relaxation cycle in accretion disks

Authors: Shalini Ganguly, Daniel Proga

Abstract: A disk wind can cause perturbations that propagate throughout the disk via diffusive processes. On reaching the inner disk, these perturbations can change the disk luminosity, which in turn, can change the wind mass loss rate, $\dot{M}_w$. It has been argued that this so-called "wind driven relaxation cycle" might explain the observed variability in some disk accreting objects. Here, we study the… ▽ More A disk wind can cause perturbations that propagate throughout the disk via diffusive processes. On reaching the inner disk, these perturbations can change the disk luminosity, which in turn, can change the wind mass loss rate, $\dot{M}_w$. It has been argued that this so-called "wind driven relaxation cycle" might explain the observed variability in some disk accreting objects. Here, we study the response of the innermost mass accretion rate $\dot{M}_a$ to the loss of matter at different rates and radii. We allow the wind launching radius, $R_L$, to scale with $\dot{M}_a$. We computed a grid of time-dependent models for various $\dot{M}_w$-$\dot{M}_a$ and $R_{L}$-$\dot{M}_a$ dependencies. We find that the disk behaviour significantly differs for the 'variable $R_L$' case compared to the 'fixed $R_L$' case. In particular, much stronger winds are required to destabilize the disk in the former than the latter case. However, the $\dot{M}_a$ amplitude does not grow significantly even for unstable cases because the oscillations saturate at a low level either due to disk depletion or due to the wind being launched at very small radii, or both. This result implies that disk winds are unlikely to be responsible for state transitions as those require large changes in the inner disk. Despite modest changes at the inner disk regions, the disk surface density at large radii can vary with a large amplitude, i.e., from 0 to a few factors of the steady state value. This dramatic variation of the outer disk could have observable consequences. △ Less

Submitted 11 January, 2020; originally announced January 2020.

Comments: 10 pages, 7 figures, 2 tables; 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

Effects of Opacity Temperature Dependence on Radiatively Accelerated Clouds

Authors: Sergei Dyda, Daniel Proga, Christopher S Reynolds

Abstract: We study how different opacity-temperature scalings affect the dynamical evolution of irradiated gas clouds using time-dependent, radiation-hydrodynamics (rad-HD) simulations. When clouds are optically thick, the bright side heats up and expands, accelerating the cloud via the rocket effect. Clouds that become more optically thick as they heat accelerate $\sim 35\%$ faster than clouds that become… ▽ More We study how different opacity-temperature scalings affect the dynamical evolution of irradiated gas clouds using time-dependent, radiation-hydrodynamics (rad-HD) simulations. When clouds are optically thick, the bright side heats up and expands, accelerating the cloud via the rocket effect. Clouds that become more optically thick as they heat accelerate $\sim 35\%$ faster than clouds that become optically thin. An enhancement of $\sim 85\%$ in the acceleration can be achieved by having a broken powerlaw opacity profile, which allows the evaporating gas driving the cloud to become optically thin and not attenuate the driving radiation flux. We find that up to $\sim 2\%$ of incident radiation is re-emitted by accelerating clouds, which we estimate as the contribution of a single accelerating cloud to an emission or absorption line. Re-emission is suppressed by "bumps" in the opacity-temperature relation since these decrease the opacity of the hot, evaporating gas, primarily responsible for the re-radiation. If clouds are optically thin, they heat nearly uniformly, expand and form shocks. This triggers the Richtmyer-Meshkov instability, leading to cloud disruption and dissipation on thermal time-scales. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 11 pages, 5 figures

A Comprehensive Chandra Study of the Disk Wind in the Black Hole Candidate 4U 1630-472

Authors: Nicolas Trueba, J. M. Miller, J. Kaastra, A. Zoghbi, A. C. Fabian T. Kallman, D. Proga, J. Raymond

Abstract: The mechanisms that drive disk winds are a window into the physical processes that underlie the disk. Stellar-mass black holes are an ideal setting in which to explore these mechanisms, in part because their outbursts span a broad range in mass accretion rate. We performed a spectral analysis of the disk wind found in six Chandra/HETG observations of the black hole candidate 4U~1630$-$472, coverin… ▽ More The mechanisms that drive disk winds are a window into the physical processes that underlie the disk. Stellar-mass black holes are an ideal setting in which to explore these mechanisms, in part because their outbursts span a broad range in mass accretion rate. We performed a spectral analysis of the disk wind found in six Chandra/HETG observations of the black hole candidate 4U~1630$-$472, covering a range of luminosities over two distinct spectral states. We modeled both wind absorption and extended wind re-emission components using PION, a self-consistent photoionized absorption model. In all but one case, two photoionization zones were required in order to obtain acceptable fits. Two independent constraints on launching radii, obtained via the ionization parameter formalism and the dynamical broadening of the re-emission, helped characterize the geometry of the wind. The innermost wind components ($r \simeq {10}^{2-3}$ $GM/{c}^{2}$) tend towards small volume filling factors, high ionization, densities up to $n \simeq {10}^{15-16} {\text{cm}}^{-3}$, and outflow velocities of $\sim 0.003c$. These small launching radii and large densities require magnetic driving, as they are inconsistent with numerical and analytical treatments of thermally driven winds. Outer wind components ($r \simeq {10}^{5}$ $GM/{c}^{2}$) are significantly less ionized and have filling factors near unity. Their larger launching radii, lower densities ($n \simeq {10}^{12} {\text{cm}}^{-3}$), and outflow velocities ($\sim 0.0007c$) are nominally consistent with thermally driven winds. The overall wind structure suggests that these components may also be part of a broader MHD outflow and perhaps better described as magneto-thermal hybrid winds. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 27 pages, 13 figures

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

A hard look at NGC 5347: revealing a nearby Compton-thick AGN

Authors: E. S. Kammoun, J. M. Miller, A. Zoghbi, K. Oh, M. Koss, R. F. Mushotzky, L. W. Brenneman, W. N. Brandt, D. Proga, A. M. Lohfink, J. S. Kaastra, D. Barret, E. Behar, D. Stern

Abstract: Current measurements show that the observed fraction of Compton-thick (CT) AGN is smaller than the expected values needed to explain the cosmic X-ray background. Prior fits to the X-ray spectrum of the nearby Seyfert-2 galaxy NGC 5347 ($z=0.00792,\, D =35.5 \rm ~Mpc $) have alternately suggested a CT and Compton-thin source. Combining archival data from $Suzaku$, $Chandra$, and - most importantly… ▽ More Current measurements show that the observed fraction of Compton-thick (CT) AGN is smaller than the expected values needed to explain the cosmic X-ray background. Prior fits to the X-ray spectrum of the nearby Seyfert-2 galaxy NGC 5347 ($z=0.00792,\, D =35.5 \rm ~Mpc $) have alternately suggested a CT and Compton-thin source. Combining archival data from $Suzaku$, $Chandra$, and - most importantly - new data from $NuSTAR$, and using three distinct families of models, we show that NGC 5347 is an obscured CTAGN ($N_{\rm H} > 2.23\times 10^{24}~\rm cm^{-2}$). Its 2-30~keV spectrum is dominated by reprocessed emission from distant material, characterized by a strong Fe K$α$ line and a Compton hump. We found a large equivalent width of the Fe K$α$ line ($\rm EW = 2.3 \pm 0.3$ keV) and a high intrinsic-to-observed flux ratio ($\sim 100$). All of these observations are typical for bona fide CTAGN. We estimate a bolometric luminosity of $L_{\rm bol} \simeq 0.014 \pm 0.005~L_{\rm Edd.}$. The $Chandra$ image of NGC 5347 reveals the presence of extended emission dominating the soft X-ray spectrum ($E < 2\,\rm keV$), which coincides with the [O III] emission detected in the $Hubble ~Space~ Telescope$ images. Comparison to other CTAGN suggests that NGC 5347 is broadly consistent with the average properties of this source class. We simulated $XRISM$ and $Athena$/X-IFU spectra of the source, showing the potential of these future missions in identifying CTAGN in the soft X-rays. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: Accepted for publication in ApJ

A global view of the inner accretion and ejection flow around super massive black holes: radiation driven accretion disk winds in a physical context

Authors: Margherita Giustini, Daniel Proga

Abstract: Understanding the physics and geometry of accretion and ejection around super massive black holes (SMBHs) is important to understand the evolution of active galactic nuclei (AGN) and therefore of the large scale structures of the Universe. We aim at providing a simple, coherent, and global view of the sub-parsec accretion and ejection flow in AGN with varying Eddington ratio, $\dot{m}$, and black… ▽ More Understanding the physics and geometry of accretion and ejection around super massive black holes (SMBHs) is important to understand the evolution of active galactic nuclei (AGN) and therefore of the large scale structures of the Universe. We aim at providing a simple, coherent, and global view of the sub-parsec accretion and ejection flow in AGN with varying Eddington ratio, $\dot{m}$, and black hole mass, $M_{BH}$. We made use of theoretical insights, results of numerical simulations, as well as UV and X-ray observations to review the inner regions of AGN by including different accretion and ejection modes, with special emphasis on the role of radiation in driving powerful accretion disk winds from the inner regions around the central SMBH. We propose five $\dot{m}$ regimes where the physics of the inner accretion and ejection flow around SMBHs is expected to change, and that correspond observationally to quiescent and inactive galaxies; low luminosity AGN (LLAGN); Seyferts and mini-broad absorption line quasars (mini-BAL QSOs); narrow line Seyfert 1 galaxies (NLS1s) and broad absorption line quasars (BAL QSOs); and super-Eddington sources. We include in this scenario radiation-driven disk winds, which are strong in the high $\dot{m}$, large $M_{BH}$ regime, and possibly present but likely weak in the moderate $\dot{m}$, small $M_{BH}$ regime. A great diversity of the accretion/ejection flows in AGN can be explained to a good degree by varying just two fundamental properties: the Eddington ratio $\dot{m}$ and the black hole mass $M_{BH}$, and by the inclusion of accretion disk winds that can naturally be launched by the radiation emitted from luminous accretion disks. △ Less

Submitted 7 November, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: Matches the published version after language editing; some references added

Journal ref: A&A 630, A94 (2019)

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

Astro2020 Science White Paper: Do Supermassive Black Hole Winds Impact Galaxy Evolution?

Authors: F. Tombesi, M. Cappi, F. Carrera, G. Chartas, K. Fukumura, M. Guainazzi, D. Kazanas, G. Kriss, D. Proga, T. J. Turner, Y. Ueda, S. Veilleux, M. Brusa, M. Gaspari

Abstract: Powerful winds driven by supermassive black holes (SMBHs) are likely the main mechanism through which SMBHs regulate their own growth and influence the host galaxy evolution. However, their origin and their capability to impact the large-scale environment are still highly debated. Fundamental results will come from high-energy and spatial resolution X-ray observatories. Powerful winds driven by supermassive black holes (SMBHs) are likely the main mechanism through which SMBHs regulate their own growth and influence the host galaxy evolution. However, their origin and their capability to impact the large-scale environment are still highly debated. Fundamental results will come from high-energy and spatial resolution X-ray observatories. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Science White Paper Submitted for the Astro2020 Decadal Survey on Astronomy and Astrophysics

The physics and astrophysics of X-ray outflows from Active Galactic Nuclei

Authors: Sibasish Laha, Randall Smith, Panayiotis Tzanavaris, Tim Kallman, Sylvain Veilleux, Francesco Tombesi, Gerard Kriss, Matteo Guainazzi, Massimo Gaspari, Jelle Kaastra, Alex Markowitz, Mike Crenshaw, Ehud Behar, Keigo Fukumura, Anna Lia Longinotti, Agata Rozanska, Jacobo Ebrero, Gary Ferland, Claudio Ricci, Chris Done, Daniel Proga, Mitchell Revalski, Andrey Vayner

Abstract: The highly energetic outflows from Active Galactic Nuclei detected in X-rays are one of the most powerful mechanisms by which the central supermassive black hole (SMBH) interacts with the host galaxy. The last two decades of high resolution X-ray spectroscopy with XMM and Chandra have improved our understanding of the nature of these outflowing ionized absorbers and we are now poised to take the n… ▽ More The highly energetic outflows from Active Galactic Nuclei detected in X-rays are one of the most powerful mechanisms by which the central supermassive black hole (SMBH) interacts with the host galaxy. The last two decades of high resolution X-ray spectroscopy with XMM and Chandra have improved our understanding of the nature of these outflowing ionized absorbers and we are now poised to take the next giant leap with higher spectral resolution and higher throughput observatories to understand the physics and impact of these outflows on the host galaxy gas. The future studies on X-ray outflows not only have the potential to unravel some of the currently outstanding puzzles in astronomy, such as the physical basis behind the MBH$-σ$ relation, the cooling flow problem in intra-cluster medium (ICM), and the evolution of the quasar luminosity function across cosmic timescales, but also provide rare insights into the dynamics and nature of matter in the immediate vicinity of the SMBH. Higher spectral resolution ($\le 0.5$ eV at $1$ keV) observations will be required to identify individual absorption lines and study the asymmetries and shifts in the line profiles revealing important information about outflow structures and their impact. Higher effective area ($\ge 1000 \rm \,cm^{2}$) will be required to study the outflows in distant quasars, particularly at the quasar peak era (redshift $1\le z\le 3$) when the AGN population was the brightest. Thus, it is imperative that we develop next generation X-ray telescopes with high spectral resolution and high throughput for unveiling the properties and impact of highly energetic X-ray outflows. A simultaneous high resolution UV + X-ray mission will encompass the crucial AGN ionizing continuum, and also characterize the simultaneous detections of UV and X-ray outflows, which map different spatial scales along the line of sight. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: A Science White Paper submitted to the Astro2020 Decadal Survey

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

Effects of Radiation Field Geometry on Line Driven Disc Winds

Authors: Sergei Dyda, Daniel Proga

Abstract: We study line driven winds for models with different radial intensity profiles: standard Shakura-Sunyaev radiating thin discs, uniform intensity discs and truncated discs where driving radiation is cutoff at some radius. We find that global outflow properties depend primarily on the total system luminosity but truncated discs can launch outflows with $\sim 2$ times higher mass flux and… ▽ More We study line driven winds for models with different radial intensity profiles: standard Shakura-Sunyaev radiating thin discs, uniform intensity discs and truncated discs where driving radiation is cutoff at some radius. We find that global outflow properties depend primarily on the total system luminosity but truncated discs can launch outflows with $\sim 2$ times higher mass flux and $\sim 50\%$ faster outflow velocity than non-truncated discs with the same total radiation flux. Streamlines interior to the truncation radius are largely unaffected and carry the same momentum flux as non-truncated models whereas those far outside the truncation radius effectively carry no outflow because the local radiation force is too weak to lift matter vertically away from the disc. Near the truncation radius the flow becomes more radial, due to the loss of pressure/radiation support from gas/radiation at larger radii. These models suggest that line driven outflows are sensitive to the geometry of the radiation field driving them, motivating the need for self-consistent disc/wind models. △ Less

Submitted 7 August, 2018; originally announced August 2018.

Comments: 10 pages, 8 figures

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

Time Dependent Radiation Driven Winds

Authors: Sergei Dyda, Daniel Proga

Abstract: We study temporal variability of radiation driven winds using one dimensional, time dependent simulations and an extension of the classic theory of line driven winds developed by Castor Abbott and Klein. We drive the wind with a sinusoidally varying radiation field and find that after a relaxation time, determined by the propagation time for waves to move out of the acceleration zone of the wind,… ▽ More We study temporal variability of radiation driven winds using one dimensional, time dependent simulations and an extension of the classic theory of line driven winds developed by Castor Abbott and Klein. We drive the wind with a sinusoidally varying radiation field and find that after a relaxation time, determined by the propagation time for waves to move out of the acceleration zone of the wind, the solution settles into a periodic state. Winds driven at frequencies much higher than the dynamical frequency behave like stationary winds with time averaged radiation flux whereas winds driven at much lower frequencies oscillate between the high and low flux stationary states. Most interestingly, we find a resonance frequency near the dynamical frequency which results in velocity being enhanced or suppressed by a factor comparable to the amplitude of the flux variation. Whether the velocity is enhanced or suppressed depends on the relative phase between the radiation and the dynamical variables. These results suggest that a time-varying radiation source can induce density and velocity perturbations in the acceleration zones of line driven winds. △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 7 pages, 3 figures

Non-Axisymmetric Line Driven Disc Winds II - Full Velocity Gradient

Authors: Sergei Dyda, Daniel Proga

Abstract: We study non-axisymetric features of 3D line driven winds in the Sobolev approximation, where the optical depth is calculated using the full velocity gradient. We find that non-axisymmetric density features, so called clumps, form primarily at the base of the wind on super-Sobolev length scales. The density of clumps differs by a factor of $\sim 3$ from the azimuthal average, the magnitude of thei… ▽ More We study non-axisymetric features of 3D line driven winds in the Sobolev approximation, where the optical depth is calculated using the full velocity gradient. We find that non-axisymmetric density features, so called clumps, form primarily at the base of the wind on super-Sobolev length scales. The density of clumps differs by a factor of $\sim 3$ from the azimuthal average, the magnitude of their velocity dispersion is comparable to the flow velocity and they produce $\sim 20\%$ variations in the column density. Clumps may be observable because differences in density produce enhancements in emission and absorption profiles or through their velocity dispersion which enhances line broadening. △ Less

Submitted 10 February, 2018; originally announced February 2018.

Comments: 12 pages, 9 figures, 2 tables

Non-Axisymmetric Line Driven Disc Winds I - Disc Perturbations

Authors: Sergei Dyda, Daniel Proga

Abstract: We study mass outflows driven from accretion discs by radiation pressure due to spectral lines. To investigate non-axisymmetric effects, we use the Athena++ code and develop a new module to account for radiation pressure driving. In 2D, our new simulations are consistent with previous 2D axisymmetric solutions by Proga et al. who used the Zeus 2D code. Specifically, we find that the disc winds are… ▽ More We study mass outflows driven from accretion discs by radiation pressure due to spectral lines. To investigate non-axisymmetric effects, we use the Athena++ code and develop a new module to account for radiation pressure driving. In 2D, our new simulations are consistent with previous 2D axisymmetric solutions by Proga et al. who used the Zeus 2D code. Specifically, we find that the disc winds are time dependent, characterized by a dense stream confined to $\sim 45^{\circ}$ relative to the disc midplane and bounded on the polar side by a less dense, fast stream. Introducing a vertical, $φ$-dependent, subsonic velocity perturbation in the disc midplane does not change the overall character of the solution but global outflow properties such as the mass, momentum and kinetic energy fluxes are altered by up to 100%. Non-axisymmetric density structures develop and persist mainly at the base of the wind. They are relatively small, and their densities can be a few times higher that the azimuthal average. The structure of the non-axisymmetric and axisymmetric solutions differ also in other ways. Perhaps most importantly from the observational point of view are the differences in the so called clumping factors, that serve as a proxy for emissivity due to two body processes. In particular, the spatially averaged clumping factor over the entire fast stream, while it is of a comparable value in both solutions, it varies about 10 times faster in the non-axisymmetric case. △ Less

Submitted 10 February, 2018; v1 submitted 21 October, 2017; originally announced October 2017.

Comments: 11 pages, 7 figures

Thermal Disk Winds in X-ray Binaries: Realistic Heating and Cooling Rates Give Rise to Slow, but Massive Outflows

Authors: Nick Higginbottom, Daniel Proga, Christian Knigge, Knox S. Long

Abstract: A number of X-ray binaries exhibit clear evidence for the presence of disk winds in the high/soft state. A promising driving mechanism for these outflows is mass loss driven by the thermal expansion of X-ray heated material in the outer disk atmosphere. Higginbottom \& Proga recently demonstrated that the properties of thermally-driven winds depend critically on the shape of the thermal equilibriu… ▽ More A number of X-ray binaries exhibit clear evidence for the presence of disk winds in the high/soft state. A promising driving mechanism for these outflows is mass loss driven by the thermal expansion of X-ray heated material in the outer disk atmosphere. Higginbottom \& Proga recently demonstrated that the properties of thermally-driven winds depend critically on the shape of the thermal equilibrium curve, since this determines the thermal stability of the irradiated material. For a given spectral energy distribution, the thermal equilibrium curve depends on exact balance between the various heating and cooling mechanisms at work. Most previous work on thermally-driven disk winds relied an analytical approximation to these rates. Here, we use the photoionization code \textsc{cloudy} to generate realistic heating and cooling rates which we then use in a 2.5D hydrodynamic model computed in ZEUS to simulate thermal winds in a typical black-hole X-ray binary. We find that these heating and cooling rates produce a significantly more complex thermal equilibrium curve, with dramatically different stability properties. The resulting flow, calculated in the optically thin limit, is qualitatively different from flows calculated using approximate analytical rates. Specifically, our thermal disk wind is much denser and slower, with a mass-loss rate that is a factor of two higher and characteristic velocities that are a factor of three lower. The low velocity of the flow -- $v_{max} \simeq 200$~km~s$^{-1}$ -- may be difficult to reconcile with observations. However, the high mass-loss rate -- 15$\times$ the accretion rate -- is promising, since it has the potential to destabilize the disk. Thermally-driven disk winds may therefore provide a mechanism for state changes. △ Less

Submitted 28 December, 2016; originally announced December 2016.

Comments: 10 Page, 8 figures, Accepted for publication in ApJ

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

Disk-Wind Connection During the Heartbeats of GRS 1915+105

Authors: Abderahmen Zoghbi, J. M. Miller, A. L. King, M. C. Miller, D. Proga, T. Kallman, A. C. Fabian, F. A. Harrison, J. Kaastra, J. Raymond, C. S. Reynolds, S. E. Boggs, F. E. Christensen, W. Craig, C. J. Hailey, D. Stern, W. W. Zhang

Abstract: Disk and wind signatures are seen in the soft state of Galactic black holes, while the jet is seen in the hard state. Here we study the disk-wind connection in the $ρ$ class of variability in GRS 1915+105 using a joint NuSTAR-Chandra observation. The source shows 50 sec limit cycle oscillations. By including new information provided by the reflection spectrum, and using phase-resolved spectroscopy… ▽ More Disk and wind signatures are seen in the soft state of Galactic black holes, while the jet is seen in the hard state. Here we study the disk-wind connection in the $ρ$ class of variability in GRS 1915+105 using a joint NuSTAR-Chandra observation. The source shows 50 sec limit cycle oscillations. By including new information provided by the reflection spectrum, and using phase-resolved spectroscopy, we find that the change in the inner disk inferred from the blackbody emission is not matched by reflection measurements. The latter is almost constant, independent of the continuum model. The two radii are comparable only if the disk temperature color correction factor changes, an effect that could be due to the changing opacity of the disk caused by changes in metal abundances. The disk inclination is similar to that inferred from the jet axis, and oscillates by ~10 deg. The simultaneous Chandra data show the presence of two wind components with velocities between 500-5000 km/s, and possibly two more with velocities reaching 20,000 km/s (~0.06 c). The column densities are ~5e22 cm$^{-2}$. An upper limit to the wind response time of 2 sec is measured, implying a launch radius of <6e10 cm. The changes in wind velocity and absorbed flux require the geometry of the wind to change during the oscillations, constraining the wind to be launched from a distance of 290 - 1300 rg from the black hole. Both datasets support fundamental model predictions in which a bulge originates in the inner disk and moves outward as the instability progresses. △ Less

Submitted 18 October, 2016; originally announced October 2016.

Comments: Accepted for publication in ApJ