Search NASA⌕ Search

Engineering topics

Waters, Tim

Publications and source records attributed to Waters, Tim.

Photoionization Calculations of the Radiation Force Due To Spectral Lines in AGNs

One of the main mechanisms that could drive mass outflows in active galactic nuclei (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, as well as 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 = 104 per cu.cm) and column density (N ≤ 1017 per sq.cm), 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 ξ gsim 1 as shown by others, this decrease in our calculations is relatively gradual and could be nonmonotonic as M(ξ, t) can increase by a factor of a few for ξ ≈ 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-local thermodynamic equilibrium effects reduce M(t, ξ) by over two orders of magnitude for ξ gsim 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 ξ gsim 100 in tenuous plasma.

Dannen, Randall C.↗

Non-isobaric Thermal Instability

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. 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. Here, 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 that is triggered when the cloud suddenly "lands" on the stable cold branch of the equilibrium curve.

79 ASTRONOMY AND ASTROPHYSICS↗

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

Magnetic, radiation pressure, and thermal driving are the three mechanisms capable of launching accretion disc 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++, here 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 disc, which is less conducive to acceleration. Finally, 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.

79 ASTRONOMY AND ASTROPHYSICS↗

Synthetic Absorption Lines for a Clumpy Medium: A Spectral Signature for Cloud Acceleration in AGN?

There is increasing evidence that the highly ionized multiphase components of AGN disc winds may be due to thermal instability. The ions responsible for forming the observed X-ray absorption lines may only exist in relatively cool 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 warm gas becomes highly mixed but not depleted. Prior to complete disruption, clouds that 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 that 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.

line: formation – galaxies: nuclei – quasars: ↗

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 blueshifted excess on both the echo image and line profile when seen from lower inclinations (i≲45{sup ∘}). This effect may be observable in low ionization lines such as Hβ.

79 ASTRONOMY AND ASTROPHYSICS↗