Source author record

Tim Waters

Tim Waters appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.HE astro-ph.SR

Catalog footprint

What is connected

9works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2024arXiv

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

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.

preprint2022arXiv

Dynamical thermal instability in highly supersonic outflows

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.

preprint2022arXiv

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

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.

preprint2020arXiv

Clumpy AGN outflows due to thermal instability

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.

preprint2019arXiv

Outflows from inflows: the nature of Bondi-like accretion

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$.

preprint2016arXiv

On the efficient acceleration of clouds in active galactic nuclei

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

preprint2016arXiv

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

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

preprint2015arXiv

Cloud Formation and Acceleration in a Radiative Environment

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

preprint2012arXiv

Parker Winds Revisited: An Extension to Disk Winds

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