BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daniel Proga

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-ph.COastro-phastro-ph.SRastro-ph.IMastro-ph.EP

Catalog footprint

What is connected

44works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

44 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the wind-driven relaxation cycle in accretion disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The role of failed accretion disk winds in active galactic nuclei

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Effects of Opacity Temperature Dependence on Radiatively Accelerated Clouds

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.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}β$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coronae and Winds from Irradiated Disks in X-ray binaries

X-ray and UV line emission in X-ray binaries can be accounted for by a hot corona. Such a corona forms through irradiation of the outer disk by radiation produced in the inner accretion flow. The same irradiation can produce a strong outflow from the disk at sufficiently large radii. Outflowing gas has been recently detected in several X-ray binaries via blue-shifted absorption lines. However, the causal connection between winds produced by irradiation and the blue-shifted absorption lines is problematic, particularly in the case of GRO J1655-40. Observations of this source imply wind densities about two orders of magnitude higher than theoretically predicted. This discrepancy does not mean that these `thermal disk-winds' cannot explain blue-shifted absorption in other systems, nor that they are unimportant as a sink of matter. Motivated by the inevitability of thermal disk-winds and wealth of data taken with current observatories such as Chandra, XMM-Newton and Suzaku, as well as the future AstroH mission, we decided to investigate the requirements to produce very dense winds. Using physical arguments, hydrodynamical simulations and absorption line calculations, we found that modification of the heating and cooling rates by a factor of a few results in an increase of the wind density of up to an order of magnitude and the wind velocity by a factor of about two. Therefore, the mass loss rate from the disk can be one, if not even two orders of magnitude higher than the accretion rate onto the central object. Such a high mass loss rate is expected to destabilize the disk and perhaps provides a mechanism for state change.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Parsec-scale accretion and winds irradiated by a quasar

We present numerical simulations of properties of a parsec-scale torus exposed to illumination by the central black hole in an active galaxy (AGN). Our physical model allows to investigate the balance between the formation of winds and accretion simultaneously. Radiation-driven winds are allowed by taking into account radiation pressure due to UV and IR radiation along with X-ray heating and dust sublimation. Accretion is allowed through angular momentum transport and the solution of the equations of radiation hydrodynamics. Our methods adopt flux-limited diffusion radiation-hydrodynamics for the dusty, infrared pressure driven part of the flow, along with X-ray heating and cooling. Angular momentum transport in the accreting part of the flow is modeled using effective viscosity. Our results demonstrate that radiation pressure on dust can play an important role in shaping AGN obscuration. For example, when the luminosity illuminating the torus exceeds $L>0.01\,L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington luminosity, we find no episodes of sustained disk accretion because radiation pressure does not allow a disk to form. Despite the absence of the disk accretion, the flow of gas to smaller radii still proceeds at a rate $10^{-4}-10^{-1}\,M_\odot\,{\rm yr}^{-1}$ through the capturing of the gas from the hot evaporative flow, thus providing a mechanism to deliver gas from a radiation-pressure dominated torus to the inner accretion disk. As $L/L_{\rm edd}$ increases, larger radiation input leads to larger torus aspect ratios and increased obscuration of the central black hole. We also find the important role of the X-ray heated gas in shaping of the obscuring torus.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Line-driven Disk Winds in Active Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer

Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGN) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman (2004). To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line-driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Virialization of Disk Winds: Implications for the Black Hole Mass Estimates in AGN

Estimating the mass of a supermassive black hole (SMBH) in an active galactic nucleus (AGN) usually relies on the assumption that the broad line region (BLR) is virialized. However, this assumption seems invalid in BLR models that consists of an accretion disk and its wind. The disk is likely Keplerian and therefore virialized. However, the wind material must, beyond a certain point, be dominated by an outward force that is stronger than gravity. Here, we analyze hydrodynamic simulations of four different disk winds: an isothermal wind, a thermal wind from an X-ray heated disk, and two line-driven winds, one with and the other without X-ray heating and cooling. For each model, we check whether gravity governs the flow properties, by computing and analyzing the volume-integrated quantities that appear in the virial theorem: internal, kinetic, and gravitational energies, We find that in the first two models, the winds are non-virialized whereas the two line-driven disk winds are virialized up to a relatively large distance. The line-driven winds are virialized because they accelerate slowly so that the rotational velocity is dominant and the wind base is very dense. For the two virialized winds, the so-called projected virial factor scales with inclination angle as $1/ \sin^2{i}$. Finally, we demonstrate that an outflow from a Keplerian disk becomes unvirialized more slowly when it conserves the gas specific angular momentum -- as in the models considered here, than when it conserves the angular velocity -- as in the so-called magneto-centrifugal winds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effects of irradiation on the cloud evolution in active galactic nuclei

We report on the first phase of our study of cloud irradiation. We study irradiation by means of numerical, two-dimensional time-dependent radiation-hydrodynamic simulations of a cloud irradiated by a strong radiation. We adopt a very simple treatment of the opacity, neglect photoionization and gravity, and instead focus on assessing the role of the type and magnitude of the opacity on the cloud evolution. Our main result is that even relatively dense clouds that are radiatively heated (i.e., with significant absorption opacity) do not move as a whole instead they undergo a very rapid and major evolution in its shape, size and physical properties. In particular, the cloud and its remnants become optical thin within less than one sound crossing time and before they can travel over a significant distance (a distance of a few radii of the initial cloud). We also found that a cloud can be accelerated as a whole under quite extreme conditions, e.g., the opacity must be dominated by scattering. However, the acceleration due to the radiation force is relatively small and unless the cloud is optically thin the cloud quickly changes its size and shape. We discuss implications for the modelling and interpetation broad line regions of active galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gas Accretion onto a Supermassive Black Hole: a step to model AGN feedback

We study the gas accretion onto a supermassive black hole (SMBH) using the 3D SPH code GADGET-3 on scales of 0.1-200 pc. First we test our code with spherically symmetric, adiabatic Bondi accretion problem. We find that our simulation can reproduce the expected Bondi accretion flow very well for a limited amount of time until the effect of outer boundary starts to be visible. We also find artificial heating of gas near the inner accretion boundary due to the artificial viscosity of SPH. Second, we implement radiative cooling and heating due to X-rays, and examine the impact of thermal feedback by the central X-ray source. The accretion flow roughly follows the Bondi solution for low central X-ray luminosities, however, the flow starts to exhibit non-spherical fragmentation due to thermal instability for a certain range of central L_X, and a strong overall outflow develops for greater L_X. The cold gas develops filamentary structures that fall into the central SMBH, whereas the hot gas tries to escape through the channels in-between the cold filaments. Such fragmentation of accreting gas can assist in the formation of clouds around AGN, induce star-formation, and contribute to the observed variability of narrow-line regions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multiphase, non-spherical gas accretion onto a black hole

(Abridged) We investigate non-spherical behavior of gas accreting onto a central supermassive black hole performing simulations using the SPH code GADGET-3 including radiative cooling and heating by the central X-ray source. As found in earlier 1D studies, our 3D simulations show that the accretion mode depends on the X-ray luminosity (L_X) for a fixed density at infinity and accretion efficiency. In the low L_X limit, gas accretes in a stable, spherically symmetric fashion. In the high L_X limit, the inner gas is significantly heated up and expands, reducing the central mass inflow rate. The expanding gas can turn into a strong enough outflow capable of expelling most of the gas at larger radii. For some intermediate L_X, the accretion flow becomes unstable developing prominent non-spherical features, the key reason for which is thermal instability (TI) as shown by our analyses. Small perturbations of the initially spherically symmetric accretion flow that is heated by the intermediate L_X quickly grow to form cold and dense clumps surrounded by overheated low density regions. The cold clumps continue their inward motion forming filamentary structures; while the hot infalling gas slows down because of buoyancy and can even start outflowing through the channels in between the filaments. We found that the ratio between the mass inflow rates of the cold and hot gas is a dynamical quantity depending on several factors: time, spatial location, and L_X; and ranges between 0 and 4. We briefly discuss astrophysical implications of such TI-driven fragmentation of accreting gas on the formation of clouds in narrow and broad line regions of AGN, the formation of stars, and the observed variability of the AGN luminiosity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the diversity and complexity of absorption line profiles produced by outflows in Active Galactic Nuclei

Understanding the origin of AGN absorption line profiles and their diversity could help to explain the physical structure of the accretion flow, and also to assess the impact of accretion on the evolution of the AGN host galaxies. Here we present our first attempt to systematically address the issue of the origin of the complexities observed in absorption profiles. Using a simple method, we compute absorption line profiles against a continuum point source for several simulations of accretion disk winds. We investigate the geometrical, ionization, and dynamical effects on the absorption line shapes. We find that significant complexity and diversity of the absorption line profile shapes can be produced by the non-monotonic distribution of the wind velocity, density, and ionization state. Non-monotonic distributions of such quantities are present even in steady-state, smooth disk winds, and naturally lead to the formation of multiple and detached absorption troughs. These results demonstrate that the part of a wind where an absorption line is formed is not representative of the entire wind. Thus, the information contained in the absorption line is incomplete if not even insufficient to well estimate gross properties of the wind such as the total mass and energy fluxes. In addition, the highly dynamical nature of certain portions of disk winds can have important effects on the estimates of the wind properties. For example, the mass outflow rates can be off up to two orders of magnitude with respect to estimates based on a spherically symmetric, homogeneous, constant velocity wind.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On time variability and other complications in studying the UV broad absorption lines of quasars: results from numerical simulations of radiation driven disk winds

We review the main results from axisymmetric, time-dependent hydrodynamical simulations of radiation driven disk winds in AGN. We illustrate the capability of such simulations to provide useful insights into the three domains of observational astronomy: spectroscopy, time-variability, and imaging. Specifically, the synthetic line profiles predicted by the simulations resemble the broad absorption lines observed in quasars. The intrinsically time dependent nature of radiation driven disk winds that have been predicted by the simulations can be supported by a growing number of the observed dramatic variability in the UV absorption lines. And finally, the intensity maps predicted by the simulations give physical and geometrical justification to the phenomenologically deduced fact that a proper interpretation of the observed line absorption requires the wind covering factor to be considered as being partial, inhomogeneous, and velocity dependent.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Parker Winds Revisited: An Extension to Disc Winds

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Regulation of Black Hole Winds and Jets Across the Mass Scale

We present a study of the mechanical power generated by both winds and jets across the black hole mass scale. We begin with the study of ionized X-ray winds and present a uniform analysis using Chandra grating spectra. The high quality grating spectra facilitate the characterization of the outflow velocity, ionization and column density of the absorbing gas. We find that the kinetic power of the winds scales with increasing bolometric luminosity as log(L_wind) \propto (1.58 \pm 0.07) log(L_Bol). This means that SMBH may be more efficient than stellar-mass black holes in launching winds. In addition, the simplicity of the scaling may suggest common driving mechanisms across the mass scale. For comparison, we next examine jet production, estimating jet power based on the energy required to inflate local bubbles. The jet relation is log(L_Jet)\propto (1.18\pm0.24) log(L_Bol). The energetics of the bubble associated with Cygnus X-1 are particularly difficult to determine, and the bubble could be a background SNR. If we exclude Cygnus X-1, then the jets follow a consistent relation to the winds within errors but with a higher normalization, log(L_Jet) \propto (1.34 \pm 0.50) log(L_Bol). The formal consistency in the wind and jet scaling relations suggests that a common launching mechanism may drive both flows; magnetic processes are viable possibilities. We also examine winds with especially high velocities, v > 0.01c. These ultra-fast outflows tend to resemble the jets more than the winds, indicating we may be observing a regime in which winds become jets. This study allows for the total power from black hole accretion, both mechanical and radiative, to be characterized in a simple manner and suggests a possible connection between winds and jets. Finally, we find at low Eddington fractions, the jet power is dominant, and at high Eddington fractions the wind power is dominant.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CIV Emission and the Ultraviolet through X-ray Spectral Energy Distribution of Radio-Quiet Quasars

In the restframe UV, two of the parameters that best characterize the range of emission-line properties in quasar broad emission-line regions are the equivalent width and the blueshift of the CIV line relative to the quasar rest frame. We explore the connection between these emission-line properties and the UV through X-ray spectral energy distribution (SED) for radio-quiet (RQ) quasars. Our sample consists of a heterogeneous compilation of 406 quasars from the Sloan Digital Sky Survey and Palomar-Green survey that have well-measured CIV emission-line and X-ray properties (including 164 objects with measured Gamma). We find that RQ quasars with both strong CIV emission and small CIV blueshifts can be classified as "hard-spectrum" sources that are (relatively) strong in the X-ray as compared to the UV. On the other hand, RQ quasars with both weak CIV emission and large CIV blueshifts are instead "soft-spectrum" sources that are (relatively) weak in the X-ray as compared to the UV. This work helps to further bridge optical/soft X-ray "Eigenvector 1" relationships to the UV and hard X-ray. Based on these findings, we argue that future work should consider systematic errors in bolometric corrections (and thus accretion rates) that are derived from a single mean SED. Detailed analysis of the CIV emission line may allow for SED-dependent corrections to these quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SPH Simulations of Black Hole Accretion: A Step to Model Black Hole Feedback in Galaxies

(Abridged) We test how accurately the smoothed particle hydrodynamics (SPH) numerical technique can follow spherically-symmetric Bondi accretion. Using the 3D SPH code GADGET-3, we perform simulations of gas accretion onto a central supermassive black hole (SMBH) of mass 10^8 M_sun within the radial range of 0.1 - 200 pc. We carry out simulations without and with radiative heating by a central X-ray corona and radiative cooling. For an adiabatic case, the radial profiles of hydrodynamical properties match the Bondi solution, except near the inner and outer radius of the computational domain. We find that adiabatic Bondi accretion can be reproduced for durations of a few dynamical times at the Bondi radius, and for longer times if the outer radius is increased. With radiative heating and cooling included, the spherically accreting gas takes a longer time to reach a steady-state than the adiabatic Bondi accretion runs, and in some cases does not reach a steady-state even within several hundred dynamical times. We find that artificial viscosity in the GADGET code causes excessive heating near the inner radius, making the thermal properties of the gas inconsistent with a physical solution. This overheating occurs typically only in the supersonic part of the flow, so that it does not affect the mass accretion rate. We see that increasing the X-ray luminosity produces a lower central mass inflow rate, implying that feedback due to radiative heating is operational in our simulations. With a sufficiently high X-ray luminosity, the inflowing gas is radiatively heated up, and an outflow develops. We conclude that the SPH simulations can capture the gas dynamics needed to study radiative feedback provided artificial viscosity alters only highly supersonic part of the inflow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Variable X-ray absorption in the mini-BAL QSO PG 1126-041

X-ray studies of active galactic nuclei (AGN) with powerful nuclear winds are important for constraining the physics of the inner accretion/ejection flow around supermassive black holes (SMBHs) and for understanding the impact of such winds on the AGN environment. Our main scientific goal is to constrain the properties of the circum-nuclear matter close to the SMBH in the mini-broad absorption line quasar (mini-BAL QSO) PG 1126-041 using a multi-epoch observational campaign with XMM-Newton. We performed temporally resolved X-ray spectroscopy and simultaneous UV and X-ray photometry on the most complete set of observations and on the deepest X-ray exposure of a mini-BAL QSO ever. We found complex X-ray spectral variability on time scales of both months and hours, which is best reproduced by means of variable massive ionized absorbers along the line of sight. As a consequence, the observed optical-to-X-ray spectral index is found to be variable with time. In the highest signal-to-noise observation we detected highly ionized X-ray absorbing material outflowing much faster (v ~ 16500 km/s) than the UV absorbing one (v ~ 5000 km/s). This highly ionized absorber is found to be variable on very short (a few kiloseconds) time scales. Our findings are qualitatively consistent with line-driven accretion disk winds scenarios. Our observations have opened the time-resolved X-ray spectral analysis field for mini-BAL QSOs. Only with future deep studies will we be able to map the dynamics of the inner flow and understand the physics of AGN winds and their impact on the environment.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Variable X-ray absorption in the mini-broad absorption line quasar PG 1126-041

We present the results of a multi-epoch observational campaign on the mini-broad absorption line quasar (mini-BAL QSO) PG 1126-041 performed with XMM-Newton from 2004 to 2009. Time-resolved X-ray spectroscopy and simultaneous UV and X-ray photometry were performed on the most complete set of observations and on the deepest X-ray exposure of a mini-BAL QSO to date. Complex X-ray spectral variability, found on time scales of both months and hours, is best reproduced by means of variable and massive ionized absorbers along the line of sight. In the highest signal-to-noise observation, highly-ionized X-ray absorbing material outflowing much faster than the UV absorbing one is detected. This highly-ionized absorber is found to be variable on very short time scales of a few hours.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large-Scale Outflows from AGN: a link between central black holes and galaxies

We summarize the results from numerical simulations of mass outflows from AGN. We focus on simulations of outflows driven by radiation from large-scale inflows. We discuss the properties of these outflows in the context of the so-called AGN feedback problem. Our main conclusion is that this type of outflows are efficient in removing matter but inefficient in removing energy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Momentum Driving: which physical processes dominate AGN feedback?

The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $η=dot{M_outflow}/dot{M_accretion}= (2 ε_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $η$, and many prior treatments have essentially assumed that $η=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $ε_w \lesssim 10^{-4}$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Outflows from AGN: their Impact on Spectra and the Environment

We present a brief summary of the main results from our multi-dimensional, time-dependent simulations of gas dynamics in AGN. We focus on two types of outflows powered by radiation emitted from the AGN: disk winds and winds driven from large-scale inflows. We show spectra predicted by the simulations and discuss their relevance to observations of broad- and narrow-line regions of the AGN. We finish with a few remarks on whether these outflows can have a significant impact on their environment and host galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The black holes of radio galaxies during the "Quasar Era": Masses, accretion rates, and evolutionary stage

We present an analysis of the AGN broad-line regions of 6 powerful radio galaxies at z>~2 (HzRGs) with rest-frame optical imaging spectroscopy obtained at the VLT. All galaxies have luminous (L(H-alpha)=few x 10^44 erg s^-1), spatially unresolved H-alpha line emission with FWHM>= 10,000 km s^-1 at the position of the nucleus, suggesting their AGN are powered by supermassive black holes with masses of few x 10^9 M_sun and accretion luminosities of a few percent of the Eddington luminosity. In two galaxies we also detect the BLRs in H-beta, suggesting relatively low extinction of A_V~1 mag, which agrees with constraints from X-ray observations. By relating black hole and bulge mass, we find a possible offset towards higher black-hole masses of at most ~0.6 dex relative to nearby galaxies at a given host mass, although each individual galaxy is within the scatter of the local relationship. If not entirely from systematic effects, this would then suggest that the masses of the host galaxies have increased by at most a factor ~4 since z~2 relative to the black-hole masses, perhaps through accretion of satellite galaxies or because of a time lag between star formation in the host galaxy and AGN fueling. We also compare the radiative and mechanical energy output (from jets) of our targets with predictions of recent models of "synthesis" or "grand unified" AGN feedback, which postulate that AGN with similar radiative and mechanical energy output rates to those found in our HzRGs may be nearing the end of their period of active growth. We discuss evidence that they may reach this stage at the same time as their host galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Magnetized Accretion Flows: Effects of Gas Pressure

We study how axisymmetric magnetohydrodynamical (MHD) accretion flows depend on gamma adiabatic index in the polytropic equation of state. This work is an extension of Moscibrodzka & Proga (2008), where we investigated the gamma dependence of 2-D Bondi-like accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if simulations for various gamma can give us insights into to the problem of various modes of accretion observed in several types of accretion systems such as black hole binaries (BHB), active galactic nuclei (AGN), and gamma-ray bursts (GRBs). We find that for gamma >~ 4/3, the fast rotating flow forms a thick torus that is supported by rotation and gas pressure. As shown before for gamma=5/3, such a torus produces a strong, persistent bipolar outflow that can significantly reduce the polar funnel accretion of a slowly rotating flow. For low gamma, close to 1, the torus is thin and is supported by rotation. The thin torus produces an unsteady outflow which is too weak to propagate throughout the polar funnel inflow. Compared to their HD counterparts, the MHD simulations show that the magnetized torus can produce an outflow and does not exhibit regular oscillations. Generally, our simulations demonstrate how the torus thickness affects the outflow production. They also support the notion that the geometrical thickness of the torus correlates with the power of the torus outflow. Our results, applied to observations, suggest that the torus ability to radiatively cool and become thin can correspond to a suppression of a jet as observed in the BHB during a transition from a hard/low to soft/high spectral state and a transition from a quiescent to hard/low state in AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Feedback Efficiency of Active Galactic Nuclei

(Abridged) We measure and analyze the energy, momentum, and mass feedback efficiencies due to radiation from AGN in relatively large scale outflows. Our measurements are based on the axisymmetric and time-dependent radiation-hydrodynamical simulations recently presented in Kurosawa & Proga. These simulations follow dynamics of gas under the influences of the AGN radiation and the gravity of the central 10^8 Msun black hole on scales from ~0.01 to ~10pc. We compare our Mdot_a-rho_o relation with that predicted by the Bondi accretion model. For high luminosities comparable to the Eddington limit, the power-law fit (Mdot_a ~ rho_o^q) to our models yields q=~0.5 instead of q=1.0 which is predicted by the Bondi model. This difference is caused by the outflows which are important for the overall mass budget at high luminosities. The maximum momentum and mass feedback efficiencies found in our models are ~10^-2 and ~10^-1, respectively. However, the outflows are much less important energetically: their thermal and kinetic powers in units of the radiative luminosity are ~10^-5 and ~10^-4, respectively. The efficiencies do not increase monotonically with the accretion luminosity but rather peak around the Eddington limit beyond which a steady state disk-wind-like solution exists. Our energy feedback efficiencies are significantly lower than 0.05, which is required in some cosmological and galaxy merger simulations. The low feedback efficiencies found here could have significant implications on the mass growth of super massive black holes in the early universe. We stress however that we have not considered the innermost parts of the accretion and outflow where radiation and matter interact most strongly. The feedback from this region could have efficiencies significantly above the low values found here.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the large-scale outflows in active galactic nuclei: consequences of coupling the mass-supply rate and accretion luminosity

(Abridged) We present two-dimensional hydrodynamical simulations of slowly rotating gas that is under the influence of the gravity of a super massive black hole and is irradiated by a thin UV accretion disc and a spherical X-ray corona. We calculate the accretion luminosity of a system based on the accretion-rate which is assumed to be equal to the mass-supply rate at the radius of ~10^-2 pc. For the models with high temperature gas at large radii (~10 pc) and high luminosities, we find a strong correlation between the mass-outflow rate (Mdot_out) and the luminosity (L). The power law index (q) describing the Mdot_out-L relation is q=2.0(+/-0.1), which is very similar to that for radiation-driven stellar and disc wind models. More surprisingly, for high density at large radii, we find steady state solutions with the accretion luminosity exceeding the Eddington limit. The super-Eddington accretion proceeds in the equatorial region and is possible because the radiation flux from the disc is significantly reduced in the equatorial direction due to the geometrical foreshortening effect. In all models, an outflow is driven from an inflow with sub-Keplerian rotation. For high temperature at large radii, the inflow occurs over a wide range of the polar angles, and the outflow occurs in a relatively narrow polar cone. However, for the super-Eddington cases with low temperature at large radii, the inflow persists only very close to the equatorial plane, resembling a thin accretion disc, while the outflow arises in a wide range of radii and polar angles. The geometry of this extreme inflow-outflow solution is very similar to a radiation-driven wind from a luminous Keplerian accretion disc.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Radiation-Driven Outflows in Active Galactic Nuclei

We review the results from multi-dimensional, time-dependent simulations of gas dynamics in AGN. We will focus on two types of outflows powered by radiation emitted from the AGN central engine: (i) outflows driven from the innermost part of an accretion disk and (2) outflows driven from a large-scale inflow that is likely the main supplier of material to the central engine. We discuss the relevance of both types of outflows to the so-called AGN feedback problem. However, the AGN feedback should not be considered separately from the AGN physics. Therefore, we also discuss the issue whether the properties of the same outflows are consistent with the gas properties in broad- and narrow-line regions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Time Evolution of the 3-D Accretion Flows: Effects of the Adiabatic Index and Outer Boundary Condition

We study a slightly rotating accretion flow onto a black hole, using the fully three dimensional (3-D)numerical simulations. We consider hydrodynamics of an inviscid flow, assuming a spherically symmetric density distribution at the outer boundary and a small, latitude-dependent angular momentum. We investigate the role of the adiabatic index and gas temperature, and the flow behaviour due to non-axisymmetric effects. Our 3-D simulations confirm axisymmetric results: the material that has too much angular momentum to be accreted forms a thick torus near the equator and the mass accretion rate is lower than the Bondi rate. In our previous study of the 3-D accretion flows, for gamma=5/3, we found that the inner torus precessed, even for axisymmetric conditions at large radii. The present study shows that the inner torus precesses also for other values of the adiabatic index: gamma=4/3, 1.2 and 1.01. However, the time for the precession to set increases with decreasing gamma. In particular, for gamma=1.01 we find that depending on the outer boundary conditions, the torus may shrink substantially due to the strong inflow of the non-rotating matter and the precession will have insufficient time to develop. On the other hand, if the torus is supplied by the continuous inflow of the rotating material from the outer radii, its inner parts will eventually tilt and precess, as it was for the larger gamma's.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Anisotropic winds from close-in extra-solar planets

We present two-dimensional hydrodynamic models of thermally driven winds from highly irradiated, close-in extra-solar planets. We adopt a very simple treatment of the radiative heating processes at the base of the wind, and instead focus on the differences between the properties of outflows in multidimensions in comparison to spherically symmetric models computed with the same methods. For hot (T > 2 x 10^{4} K) or highly ionized gas, we find strong (supersonic) polar flows are formed above the planet surface which produce weak shocks and outflow on the night-side. In comparison to a spherically symmetric wind with the same parameters, the sonic surface on the day-side is much closer to the planet surface in multidimensions, and the total mass loss rate is reduced by almost a factor of four. We also compute the steady-state structure of interacting planetary and stellar winds. Both winds end in a termination shock, with a parabolic contact discontinuity which is draped over the planet separating the two shocked winds. The planetary wind termination shock and the sonic surface in the wind are well separated, so that the mass loss rate from the planet is essentially unaffected. However, the confinement of the planetary wind to the small volume bounded by the contact discontinuity greatly enhances the column density close to the planet, which might be important for the interpretation of observations of absorption lines formed by gas surrounding transiting planets.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Three-Dimensional Simulations of Dynamics of Accretion Flows Irradiated by a Quasar

We study the axisymmetric and non-axisymmetric, time-dependent hydrodynamics of gas that is under the influence of the gravity of a super massive black hole (SMBH) and the radiation force produced by a radiatively efficient flow accreting onto the SMBH. We have considered two cases: (1) the formation of an outflow from the accretion of the ambient gas without rotation and (2) that with weak rotation. The main goals of this study are: (1) to examine if there is a significant difference between the models with identical initial and boundary conditions but in different dimensionality (2-D and 3-D), and (2) to understand the gas dynamics in AGN. Our 3-D simulations of a non-rotating gas show small yet noticeable non-axisymmetric small-scale features inside the outflow. The outflow as a whole and the inflow do not seem to suffer from any large-scale instability. In the rotating case, the non-axisymmetric features are very prominent, especially in the outflow which consists of many cold dense clouds entrained in a smoother hot flow. The 3-D outflow is non-axisymmetric due to the shear and thermal instabilities. In both 2-D and 3-D simulations, gas rotation increases the outflow thermal energy flux, but reduces the outflow mass and kinetic energy fluxes. Rotation also leads to time variability and fragmentation of the outflow in the radial and latitudinal directions. The collimation of the outflow is reduced in the models with gas rotation. The time variability in the mass and energy fluxes is reduced in the 3-D case because of the outflow fragmentation in the azimuthal direction. The virial mass estimated from the kinematics of the dense cold clouds found in our 3-D simulations of rotating gas underestimates the actual mass used in the simulations by about 40 %. (Abbreviated)

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Dynamics of Accretion Flows Irradiated by a Quasar

We present the results from axisymmetric time-dependent HD calculations of gas flows which are under the influence of gravity of a black hole in quasars. We assume that the flows are non-rotating and exposed to quasar radiation. We take into account X-ray heating and the radiation force due to electron scattering and spectral lines. To compute the radiation field, we consider a standard accretion disk as a source of UV photons and a spherical central object as a source of X-rays. The gas temperature and ionization state in the flow are calculated self-consistently from the photoionization and heating rate of the central object. We find that for a 10e8 MSUN black hole with an accretion luminosity of 0.6 of the Eddington luminosity the flow settles into a steady state and has two components: (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the disk rotational axis. The inflow is a realization of a Bondi-like accretion flow. The second component is an example of a non-radial accretion flow which becomes an outflow once it is pushed close to the rotational axis where thermal expansion and the radiation pressure accelerate it outward. Our main result is that the existence of the above two flow components is robust to the outer boundary conditions and the geometry and spectral energy distribution of the radiation field. However, the flow properties are not robust. In particular, the outflow power and collimation is higher for the radiation dominated by the UV/disk emission than for the radiation dominated by the X-ray/central engine emission. Our most intriguing result is that a very narrow outflow driven by radiation pressure on lines can carry more energy and mass than a broad outflow driven by thermal expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Chemical Abundances in AGN Environment: X-Ray/UV Campaign on the MRK 279 Outflow

We present the first reliable determination of chemical abundances in an AGN outflow. The abundances are extracted from the deep and simultaneous FUSE and HST/STIS observations of Mrk 279. This data set is exceptional for its high signal-to-noise, unblended doublet troughs and little Galactic absorption contamination. These attributes allow us to solve for the velocity-dependent covering fraction, and therefore obtain reliable column densities for many ionic species. For the first time we have enough such column densities to simultaneously determine the ionization equilibrium and abundances in the flow. Our analysis uses the full spectral information embedded in these high-resolution data. Slicing a given trough into many independent outflow elements yields the extra constraints needed for a physically meaningful abundances determination. We find that relative to solar the abundances in the Mrk 279 outflow are (linear scaling): carbon 2.2+/-0.7, nitrogen 3.5+/-1.1 and oxygen 1.6+/-0.8. Our UV-based photoionization and abundances results are in good agreement with the independent analysis of the simultaneous Mrk 279 X-ray spectra. This is the best agreement between the UV and X-ray analyses of the same outflow to date.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Dynamics of radiatively inefficient flows accreting onto radiatively efficient black hole objects

I present results from numerical simulations of gas dynamics outside luminous accretion disks in active galactic nuclei. The gas, gravitationally captured by a super massive black hole, can be driven away by the energy and momentum of the radiation emitted during black hole accretion. Assuming axisymmetry, I study how the mass accretion and outflow rates, and the flow dynamics respond to changes in radiation heating relative to radiation pressure. I find that for a 10^8 MSUN black hole with the accretion luminosity of 0.6 of the Eddington luminosity the flow settles into a steady state and has two components: (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the disk rotational axis. The inflow is a realization of a Bondi-like accretion flow. The second component is an example of a non-radial accretion flow which becomes an outflow once it is pushed close to the rotational axis where thermal expansion and radiation pressure accelerate it outward. The main result of this preliminary work is that although the above two-component solution is robust, its properties are sensitive to the geometry and spectral energy distribution of the radiation field.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The late time evolution of Gamma-Ray Bursts: ending hyperaccretion and producing flares

We consider the properties of a hyperaccretion model for gamma-ray bursts (GRBs) at the late time when the mass supply rate is expected to decrease with time. We point out that the region in the vicinity of the accretor and the accretor itself can play an important role in determining the rate and time behavior of the accretion and ultimately the energy output. Motivated by numerical simulations and theoretical results, we conjecture that the energy release can be repeatedly stopped and then restarted by the magnetic flux accumulated around the accretor. We propose that the episode or episodes when the accretion resumes correspond to X-ray flares discovered recently in a number of GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Radiation driven winds from luminous accretion disks

We study the 2-D, time-dependent hydrodynamics of radiation-driven winds from accretion disks in which the radiation force is mediated by spectral lines. If the dominant contribution to the total radiation field comes from the disk, then we find the outflow is intrinsically unsteady and characterised by large amplitude velocity and density changes. Both infall and outflow can occur in different regions of the wind at the same time. On the other hand, if the total luminosity of the system is dominated by the central star, then the outflow is steady. In either case, we find the 2-D structure of the wind consists of a dense, slow outflow, typically confined to angles within about 45$^o$ of the equatorial plane, that is bounded on the polar side by a high -velocity, lower density stream. The flow geometry is controlled largely by the geometry of the radiation field. Global properties e.g., the total mass loss rate and terminal velocity depend more on the system luminosity and are insensitive to geometry. Matter is fed into the fast wind from within a few stellar radii of the central star. Our solutions agree qualitatively with the kinematics of outflows in CV systems inferred from spectroscopic observations. We predict that low luminosity systems may display unsteady behavior in wind-formed spectral lines. Our study also has application to winds from active galactic nuclei and from high mass YSOs.

0 citations0 reviewsNo code linkedOpen access