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.

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.

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

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{\rho}}/\partial{\theta}|$ 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$.

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(\xi,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(\xi,t)$ starts to decrease with $\xi$ for $\xi \gtrsim 1$ as shown by others, this decrease in our calculations is relatively gradual and could be non-monotonic as $M(\xi,t)$ can increase by a factor of few for $\xi \approx 10-1000$; 2) at these same $\xi$ 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,\xi)$ by over two orders of magnitude for $\xi \gtrsim 100$. The dynamical consequence of result (1) is that line driving can be important for $\xi$ 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 $\xi\gtrsim100$ in tenuous plasma.

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

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.

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.

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.

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

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.

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}\beta$.

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.

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.

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

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

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.