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

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.

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

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

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.

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.

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.

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

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

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.

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

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.

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

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

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

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

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

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

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.

NGC 2617 is a nearby ($z\sim 0.01$) active galaxy that recently switched from being a Seyfert 1.8 to be a Seyfert 1.0. At the same time, it underwent a strong increase of X-ray flux by one order of magnitude with respect to archival measurements. We characterise the X-ray spectral and timing properties of NGC 2617 with the aim of studying the physics of a changing-look active galactic nucleus (AGN). We performed a comprehensive timing and spectral analysis of two XMM-Newton pointed observations spaced by one month, complemented by archival quasi-simultaneous INTEGRAL observations. We found that, to the first order, NGC 2617 looks like a type 1 AGN in the X-ray band and, with the addition of a modest reflection component, its continuum can be modelled well either with a power law plus a phenomenological blackbody, a partially covered power law, or a double Comptonisation model. Independent of the continuum adopted, in all three cases a column density of a few $10^{23}$ cm$^{-2}$ of neutral gas covering 20-40\% of the continuum source is required by the data. Most interestingly, absorption structures due to highly ionised iron have been detected in both observations with a redshift of about $0.1c$ with respect to the systemic redshift of the host galaxy. The redshifted absorber can be ascribed to a failed wind/aborted jets component, to gravitational redshift effects, and/or to matter directly falling towards the central supermassive black hole. In either case, we are probing the innermost accretion flow around the central supermassive black hole of NGC 2617 and might be even watching matter in a direct inflow towards the black hole itself.

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.

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

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

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.

We report on a 120 ks Chandra/HETG spectrum of the black hole GRS 1915+105. The observation was made during an extended and bright soft state in June, 2015. An extremely rich disk wind absorption spectrum is detected, similar to that observed at lower sensitivity in 2007. The very high resolution of the third-order spectrum reveals four components to the disk wind in the Fe K band alone; the fastest has a blue-shift of v = 0.03c. Broadened re-emission from the wind is also detected in the first-order spectrum, giving rise to clear accretion disk P Cygni profiles. Dynamical modeling of the re-emission spectrum gives wind launching radii of r ~ 10^(2-4) GM/c^2. Wind density values of n ~ 10^(13-16) cm^-3 are then required by the ionization parameter formalism. The small launching radii, high density values, and inferred high mass outflow rates signal a role for magnetic driving. With simple, reasonable assumptions, the wind properties constrain the magnitude of the emergent magnetic field to B ~ 10^(3-4) Gauss if the wind is driven via magnetohydrodynamic (MHD) pressure from within the disk, and B ~ 10^(4-5) Gauss if the wind is driven by magnetocentrifugal acceleration. The MHD estimates are below upper limits predicted by the canonical alpha-disk model (Shakura & Sunyaev 1973). We discuss these results in terms of fundamental disk physics and black hole accretion modes.

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.

We present an analysis of ionized X-ray disk winds observed in the Fe K band of four stellar-mass black holes observed with Chandra, including 4U 1630-47, GRO J1655-40, H 1743-322, and GRS 1915+105. High-resolution photoionization grids were generated in order to model the data. Third-order gratings spectra were used to resolve complex absorption profiles into atomic effects and multiple velocity components. The Fe XXV line is found to be shaped by contributions from the intercombination line (in absorption), and the Fe XXVI line is detected as a spin-orbit doublet. The data require 2-3 absorption zones, depending on the source. The fastest components have velocities approaching or exceeding 0.01c, increasing mass outflow rates and wind kinetic power by orders of magnitude over prior single-zone models. The first-order spectra require re-emission from the wind, broadened by a degree that is loosely consistent with Keplerian orbital velocities at the photoionization radius. This suggests that disk winds are rotating with the orbital velocity of the underlying disk, and provides a new means of estimating launching radii -- crucial to understanding wind driving mechanisms. Some aspects of the wind velocities and radii correspond well to the broad-line region (BLR) in active galactic nuclei, suggesting a physical connection. We discuss these results in terms of prevalent models for disk wind production and disk accretion itself, and implications for massive black holes in active galactic nuclei.

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.

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.

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.

We report on a Chandra/HETG X-ray spectrum of the black hole candidate MAXI J1305-704. A rich absorption complex is detected in the Fe L band, including density-sensitive lines from Fe XX, XXI, and XXII. Spectral analysis over three bands with photoionization models generally requires a gas density of n > 1 E+17 cm^-3. Assuming a luminosity of L = 1 E+37 erg/s, fits to the 10-14 A band constrain the absorbing gas to lie within r = 3.9(7) E+3 km from the central engine, or about r = 520 +/- 90 (M/5 Msun) r_g, where r_g = GM/c^2. At this distance from the compact object, gas in Keplerian orbits should have a gravitational red-shift of z = v/c ~ 3 +/- 1 E-3 (M/5 Msun), and any tenuous inflowing gas should have a free-fall velocity of v/c ~ 6 +/- 1 E-2 (M/5 Msun)^1/2. The best-fit single-zone photoionization models measure a red-shift of v/c = 2.6-3.2 E-3. Models with two zones provide significantly improved fits; the additional zone is measured to have a red-shift of v/c =4.6-4.9 E-2 (models including two zones suggest slightly different radii and may point to lower densities). Thus, the shifts are broadly consistent with the photoionization radius. The results may be explained in terms of a "failed wind" like those predicted in some numerical simulations. We discuss our results in the context of accretion flows across the mass scale, and the potential role of failed winds in black hole state transitions.

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.

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.

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.

Understanding the astrophysics of feedback in active galactic nuclei (AGN) is key to understanding the growth and co-evolution of supermassive black holes and galaxies. AGN-driven winds/outflows are potentially the most effective way of transporting energy and momentum from the nuclear scales to the host galaxy, quenching star formation by sweeping away the gas reservoir. Key questions in this field are: 1) how do accretion disks around black holes launch winds/outflows, and how much energy do these carry? 2) How are the energy and metals accelerated in winds/outflows transferred and deposited into the circumgalactic medium? X-ray observations are a unique way to address these questions because they probe the phase of the outflows which carries most of the kinetic energy. We show how a high throughput, high spectral resolution instrument like the X-ray Integral Field Unit (X-IFU) on Athena+ will allow us to address these questions by determining the physical parameters (ionization state, density, temperature, abundances, velocities, geometry, etc.) of the outflows on a dynamical time-scale, in a broad sample of nearby bright AGN. The X-IFU will also allow direct spectral imaging of the impact of these winds on the host galaxy for local AGN, forming a template for understanding AGN at higher redshifts where wind shocks cannot be resolved.

(Abridged) We study stability of gas accretion in Active Galactic Nuclei. Our grid based simulations cover a radial range from 0.1 to 200 pc. Here, as in previous studies by our group, we include gas radiative cooling as well as heating by a sub-Eddington X-ray source near the central supermassive black hole of 10^8 M_{\odot}. Our theoretical estimates and simulations show that for the X-ray luminosity L_X \sim 0.008 L_{Edd}, the gas is thermally and convectivelly unstable within the computational domain. In the simulations, we observe that very tiny fluctuations in an initially smooth, spherically symmetric, accretion flow, grow first linearly and then non-linearly. Consequently, an initially one-phase flow relatively quickly transitions into a two-phase/cold-hot accretion flow. For L_X = 0.015 L_{Edd} or higher, the cold clouds continue to accrete but in some regions of the hot phase, the gas starts to move outward. For L_X < 0.015 L_{Edd}, the cold phase contribution to the total mass accretion rate only moderately dominates over the hot phase contribution. This result might have some consequences for cosmological simulations of the so-called AGN feedback problem. Our simulations confirm the previous results of Barai et al. (2012) who used smoothed particle hydrodynamic simulations to tackle the same problem. However here, because we use a grid based code to solve equations in 1-D and 2-D, we are able to follow the gas dynamics at much higher spacial resolution and for longer time in comparison to the 3-D SPH simulations. One of new features revealed by our simulations is that the cold condensations in the accretion flow initially form long filaments, but at the later times, those filaments may break into smaller clouds advected outwards within the hot outflow. These simulations may serve as an attractive model for the so-called Narrow Line Region in AGN.

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.

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.

Ultraviolet and X-ray observations show evidence of outflowing gas around many active galactic nuclei. Some of these outflows may be driven off gas infalling towards the central black hole. We perform radiative transfer calculations to compute the gas ionization state and X-ray spectra for two- and three-dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By comparison with observations, our results can be used to test the theoretical models and guide future numerical simulations. We predict both absorption and emission features, most of which are formed in a polar funnel of outflowing gas. This outflow causes strong absorption for observer orientation angles of < 35 degrees. Particularly in 3D, the strength of this absorption varies significantly for different lines-of-sight owing to the fragmentary structure of the gas flow. Although infalling material occupies a large fraction of the simulation volume, we do not find that it imprints strong absorption features since the ionization state is very high. Thus, an absence of observed inflow absorption features does not exclude the models. The main spectroscopic consequence of the infalling gas is a scattered continuum component that partially re-fills the absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the outflow is predicted to give rise to several emission features for all observer orientations. For the hydrodynamical simulations considered we find both ionization states and column densities for the outflowing gas that are too high to be quantitatively consistent with well-observed X-ray absorption systems. Nevertheless, our results are qualitatively encouraging and further exploration of the model parameter space is warranted. (Abridged.)

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.

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

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.

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.

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.

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.

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.

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.

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

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.

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 $\eta=dot{M_outflow}/dot{M_accretion}= (2 \epsilon_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 $\eta$, and many prior treatments have essentially assumed that $\eta=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 $\epsilon_w \lesssim 10^{-4}$.

We perform multi-dimensional radiative transfer simulations to compute spectra for a hydrodynamical simulation of a line-driven accretion disk wind from an active galactic nucleus. The synthetic spectra confirm expectations from parameterized models that a disk wind can imprint a wide variety of spectroscopic signatures including narrow absorption lines, broad emission lines and a Compton hump. The formation of these features is complex with contributions originating from many of the different structures present in the hydrodynamical simulation. In particular, spectral features are shaped both by gas in a successfully launched outflow and in complex flows where material is lifted out of the disk plane but ultimately falls back. We also confirm that the strong Fe Kalpha line can develop a weak, red-skewed line wing as a result of Compton scattering in the outflow. In addition, we demonstrate that X-ray radiation scattered and reprocessed in the flow has a pivotal part in both the spectrum formation and determining the ionization conditions in the wind. We find that scattered radiation is rather effective in ionizing gas which is shielded from direct irradiation from the central source. This effect likely makes the successful launching of a massive disk wind somewhat more challenging and should be considered in future wind simulations.

We present the results of hydrodynamical simulations of the disk photosphere irradiated by strong X-rays produced in the inner most part of the disk. As expected, the irradiation heats the photosphere and drives a thermal wind. To apply our results to the well-studied X-ray transient source GRO J1655-40, we adopted the observed mass of its black hole, and the observed properties of its X-ray radiation. To compare the results with the observations, we also computed transmitted X-ray spectra based on the wind solution. Our main finding is: the density of the fast moving part of the wind is more than one order of magnitude lower than that inferred from the observations. Consequently, the model fails to predict spectra with line absorption as strong and as blueshifted as those observed. However, despite the thermal wind being weak and Compton thin, the ratio between the mass-loss rate and the mass accretion rate is about seven. This high ratio is insensitive to the accretion luminosity, in the limit of lower luminosities. Most of the mass is lost from the disk between 0.07 and 0.2 of the Compton radius. We discovered that beyond this range the wind solution is self-similar. In particular, soon after it leaves the disk, the wind flows at a constant angle with respect to the disk. Overall, the thermal winds generated in our comprehensive simulations do not match the wind spectra observed in GRO J1655-40. This supports the conclusion of Miller et al. and Kallman et al. that the wind in GRO J1655-40, and possibly other X-ray transients, may be driven by magnetic processes. This in turn implies that the disk wind carries even more material than our simulations predict and as such has a very significant impact on the accretion disk structure and dynamics.

We find, from high-resolution hydro simulations, that winds from AGN effectively heat the inner parts (~100 pc) of elliptical galaxies, reducing infall to the central SMBH; and radiative (photoionization and X-ray) heating reduces cooling flows at the kpc scale. Including both types of feedback with (peak) efficiencies of 3 10^{-4} < epsilon_mech < 10^{-3} and of epsilon_rad ~10^{-1.3} respectively, produces systems having duty-cycles, central SMBH masses, X-ray luminosities, optical light profiles, and E+A spectra in accord with the broad suite of modern observations of massive elliptical systems. Our main conclusion is that mechanical feedback (including all three of energy, momentum and mass) is necessary but the efficiency, based on several independent arguments must be a factor of 10 lower than is commonly assumed. Bursts are frequent at z>1 and decline in frequency towards the present epoch as energy and metal rich gas are expelled from the galaxies into the surrounding medium. For a representative galaxy of final stellar mass ~3 10^{11} Msun, roughly 3 10^{10} Msun of recycled gas has been added to the ISM since z~2 and, of that, roughly 63% has been expelled from the galaxy, 19% has been converted into new metal rich stars in the central few hundred parsecs, and 2% has been added to the central SMBH, with the remaining 16% in the form hot X-ray emitting ISM. The bursts occupy a total time of ~170 Myr, which is roughly 1.4% of the available time. Of this time, the central SMBH would be seen as an UV or optical source for ~45% and ~71$% of the time, respectively. Restricting to the last 8.5 Gyr, the burst occupy ~44 Myr, corresponding to a fiducial duty-cycle of ~5 10^{-3}.

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.

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.

We present the first results from a long (496 ks) Chandra High Energy Transmission Grating observation of the intermediate polar EX Hydrae. In addition to the narrow emission lines from the cooling post-shock gas, for the first time we have detected a broad component in some of the X-ray emission lines, namely O VIII 18.97, Mg XII 8.42, Si XIV 6.18, and Fe XVII 16.78. The broad and narrow components have widths of ~ 1600 km s^-1 and ~ 150 km s^-1, respectively. We propose a scenario where the broad component is formed in the pre-shock accretion flow, photoionized by radiation from the post-shock flow. Because the photoionized region has to be close to the radiation source in order to produce strong photoionized emission lines from ions like O VIII, Fe XVII, Mg XII, and Si XIV, our photoionization model constrains the height of the standing shock above the white dwarf surface. Thus, the X-ray spectrum from EX Hya manifests features of both magnetic and non-magnetic cataclysmic variables.

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.

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

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.

(Abridged) The two high-redshift gamma-ray bursts, GRB 080913 at z=6.7 and GRB 090423 at z=8.3, recently detected by Swift appear as intrinsically short, hard GRBs. They could have been recognized by BATSE as short/hard GRBs should they have occurred at z <= 1. We perform a more thorough investigation on two physically distinct types (Type I/II) of cosmological GRBs and their observational characteristics. We reiterate the definitions of Type I/II GRBs and review the observational criteria and their physical motivations. Contrary to the traditional approach of assigning the physical category based on the gamma-ray properties (duration, hardness, and spectral lag), we take an alternative approach to define the Type I and Type II Gold Samples using several criteria that are more directly related to the GRB progenitors, and study the properties of the two Gold Samples and compare them with the traditional long/soft and short/hard samples. We find that the Type II Gold Sample reasonably tracks the long/soft population, although it includes several intrinsically short (shorter than 1s in the rest frame) GRBs. The Type I Gold Sample only has 5 GRBs, 4 of which are not strictly short but have extended emission. Other short/hard GRBs detected in the Swift era represent the BATSE short/hard sample well, but it is unclear whether all of them belong to Type I. We suggest that some (probably even most) high-luminosity short/hard GRBs instead belong to Type II. We suggest that GRB 080913 and GRB 090423 are more likely Type II events. We re-emphasize the importance of invoking multiple observational criteria, and cautiously propose an operational procedure to infer the physical origin of a given GRB with available multiple observational criteria, with various caveats laid out.

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.

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

The importance of the radiative feedback from SMBHs at the centers of elliptical galaxies is not in doubt, given the well established relations among electromagnetic output, black hole mass and galaxy optical luminosity. In addition, feedback due to mechanical and thermal deposition of energy from jets and winds emitted by the accretion disk around the central SMBH is also expected to occur. In this paper we improve and extend the accretion and feedback physics explored in our previous papers to include also a physically motivated mechanical feedback. We study the evolution of an isolated elliptical galaxy with the aid of a high-resolution 1-D hydrodynamical code, where the cooling and heating functions include photoionization and Compton effects, and restricting to models which include only radiative or only mechanical feedback. We confirm that for Eddington ratios above 0.01 both the accretion and radiative output are forced by feedback effects to be in burst mode, so that strong intermittencies are expected at early times, while at low redshift the explored models are characterized by smooth, very sub-Eddington mass accretion rates punctuated by rare outbursts. However, the explored models always fail some observational tests. If we assume the high mechanical efficiency of 10^{-2.3}, we find that most of the gas is ejected from the galaxy, the resulting X-ray luminosity is far less than is typically observed and little SMBH growth occurs. But models with low enough mechanical efficiency to accomodate satisfactory SMBH growth tend to allow too strong cooling flows and leave galaxies at z=0 with E+A spectra more frequently than is observed. We conclude that both types of feedback are required. Models with combined feedback are explored in a forthcoming paper [abridged]

An extraordinary recent development in astrophysics was the discovery of the fossil relationship between central black hole mass and the stellar mass of galactic bulges. The physical process underpinning this relationship has become known as feedback. The Chandra X-ray Observatory was instrumental in realizing the physical basis for feedback, by demonstrating a tight coupling between the energy released by supermassive black holes and the gaseous structures surrounding them. This white paper discusses how a great leap forward in X-ray collecting area and spectral resolution will allow a qualitatively new way of studying how feedback from black holes influenced the growth of structure.

The large scale structure of the present Universe is determined by the growth of dark matter density fluctuations and by the dynamical action of dark energy and dark matter. While much progress has been made in recent years in constraining the cosmological parameters, and in reconstructing the evolution in the large--scale structure of the dark matter distribution, we still lack an understanding of the evolution of the baryonic component of the Universe. Located at nodes of the cosmic web, clusters of galaxies are the largest collapsed structures in the Universe with total masses up to 10$^{15}$ M$_{\sun}$. Over 80% of their mass resides in the form of dark matter. The remaining mass is composed of baryons, most of which (about 85%) is a diffuse, hot plasma that radiates primarily in X-rays. X-ray observations of the evolving cluster population provide a unique opportunity to address such open and fundamental questions as: How do hot diffuse baryons dynamically evolve in dark matter potentials? How and when was the excess energy which we observe in the intergalactic medium generated? What is the cosmic history of heavy-element production and circulation? Our current knowledge comes primarily from detailed studies of clusters in the relatively nearby Universe (z$<$0.5). Major advances will come from high throughput, high spectral and spatial resolution X-ray observations that measure the thermodynamical properties and metal content of the first low mass clusters emerging at z $\sim$ 2 and directly trace their evolution into today's massive clusters.

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)

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.

abridged: We use XSCORT, together with the hydrodynamic accretion disc wind simulation from Proga & Kallman (2004), to calculate the impact that the accretion disk wind has on the X-ray spectrum from a 1E8 solar mass black hole Active Galactic Nuclei (AGN) accreting at 0.5 L/L_Edd. The properties of the resulting spectra depend on viewing angle and clearly reflect the distinct regions apparent in the original hydrodynamic simulation. Very equatorial lines-of-sight (l.o.s) are dominated by Compton scattering and nearly-neutral absorption. Polar l.o.s result in largely featureless spectra. Finally, l.o.s that intersect the transition region between these extremes have a wide range of absorption features imprinted on the spectrum. Both polar and transition region l.o.s produce spectra that show highly-ionized, blue-shifted, Fe absorption features that are qualitatively similar to features observed in the X-ray spectra of a growing number of AGN. The spectra presented here clearly demonstrate that current simulations of line driven AGN accretion disk winds cannot reproduce the smooth soft X-ray excess. Furthermore, they predict that high accretion rate (L/L_Edd) AGN are likely to be strongly affected by obscuration, in sharp contrast to the clean picture that is generally assumed, based on the observed relation between the opening angle of the molecular torus and AGN luminosity.

In the frame of the collapsar model for long gamma ray bursts (GRBs), we investigate the formation of a torus around a spinning BH and we check what rotational properties a progenitor star must have in order to sustain torus accretion over relatively long activity periods. We also study the time evolution of the BH spin parameter. We take into account the coupling between BH mass, its spin parameter and the critical specific angular momentum of accreting gas, needed for the torus to form. The large BH spin reduces the critical angular momentum which in turn can increase the GRB duration with respect to the Schwarzschild BH case. We quantify this effect and estimate the GRB durations in three cases: when a hyper accreting torus operates or a BH spins very fast or both. We show under what conditions a given progenitor star produces a burst that can last as short as several seconds and as long as several hundred of seconds. Our models indicate that it is possible for a single collapse to produce three kinds of jets: (1) a very short, lasting between a fraction of a second and a few seconds, 'precursor' jet, powered only by a hyper accreting torus before the BH spins up, (2) an 'early' jet, lasting several tens of seconds and powered by both hyper accretion and BH rotation, and (3) a 'late' jet, powered only by the spinning BH.

We have calculated a series of models of outflows from the obscuring torus in active galactic nuclei (AGN). Our modeling assumes that the inner face of a rotationally supported torus is illuminated and heated by the intense X-rays from the inner accretion disk and black hole. As a result of such heating a strong biconical outflow is observed in our simulations. We calculate 3-dimensional hydrodynamical models, assuming axial symmetry, and including the effects of X-ray heating, ionization, and radiation pressure. We discuss the behavior of a large family of these models, their velocity fields, mass fluxes and temperature, as functions of the torus properties and X-ray flux. Synthetic warm absorber spectra are calculated, assuming pure absorption, for sample models at various inclination angles and observing times. We show that these models have mass fluxes and flow speeds which are comparable to those which have been inferred from observations of Seyfert 1 warm absorbers, and that they can produce rich absorption line spectra.

We report on the fourth phase of our study of slightly rotating accretion flows onto black holes. The main new element of this study is that we used fully three dimensional (3-D) numerical simulations. We consider hydrodynamics of inviscid accretion flows. We assume a spherically symmetric density distribution at the outer boundary, but brake the flow symmetry by introducing a small, latitude-dependent angular momentum. We also consider cases where angular momentum at large radii is latitude- and azimuth-dependent. For the latitude-dependent angular momentum, 3-D simulations confirm axisymmetric results: the material that has too much angular momentum to be accreted forms a thick torus near the equator. Consequently, accretion proceeds only through the polar funnel, and the mass accretion rate through the funnel is constrained by the size and shape of the torus, not by the outer conditions. In 3-D simulations, we found that the torus precesses, even for axisymmetric conditions at large radii. For the latitude and azimuth-dependent angular momentum, the non-rotating gas near the equator can also significantly affect the evolution of the rotating gas. In particular, it may prevent the formation of a proper torus (i.e. its closing, in the azimuthal direction). In such models, the mass accretion rate is only slightly less than the corresponding Bondi rate.

We report on time-dependent axisymmetric simulations of an X-ray excited flow from a parsec-scale, rotating, cold torus around an active galactic nucleus. Our simulations account for radiative heating and cooling and radiation pressure force. The simulations follow the development of a broad bi-conical outflow induced mainly by X-ray heating. We compute synthetic spectra predicted by our simulations. The wind characteristics and the spectra support the hypothesis that a rotationally supported torus can serve as the source of a wind which is responsible for the warm absorber gas observed in the X-ray spectra of many Seyfert galaxies.

We report on next phase of our study of rotating accretion flows onto black holes. We consider hydrodynamical (HD) accretion flows with a spherically symmetric density distribution at the outer boundary but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum. We study accretion flows by means of numerical two-dimensional, axisymmetric, HD simulations for variety of the adiabatic index, $\gamma$ and the gas temperature at infinity, $c_\infty$. Our work is an extension of work done by Proga & Begelman who consider models for only $\gamma=5/3$. Our main result is that the flow properties such as the topology of the sonic surface and time behavior strongly depend on $\gamma$ but little on $c_\infty$. In particular, for $1 < \gamma < 5/3$, the mass accretion rate shows large amplitude, slow time-variability which is a result of mixing between slow and fast rotating gas. This temporal behavior differs significantly from that in models with $\gamma\simless 5/3$ where the accretion rate is relatively constant and from that in models with $\gamma\simgreat 1$ where the accretion exhibits small amplitude quasi-periodic oscillations. The key parameter responsible for the differences is the sound speed of the accretion flow which in turn determines whether the flow is dominated by gas pressure, radiation pressure or rotation. Despite these differences the time-averaged mass accretion rate in units of the corresponding Bondi rate is a weak function of $\gamma$ and $c_\infty$.

Sgr A* is a source of strongly variable emission in several energy bands. It is generally agreed that this emission comes from the material surrounding the black hole which is either falling in or flowing out. The activity must be driven by accretion but the character of accretion flow in this object is an open question. We suggest that the inflow is dominated by the relatively low angular momentum material originating in one of the nearby group of stars. Such material flows in directly towards the black hole up to the distance of order of ten Schwarzschild radii or less, where it hits the angular momentum barrier which leads naturally to a flow variability. We study both the analytical and the numerical solutions for the flow dynamics, and we analyze the radiation spectra in both cases using the Monte Carlo code to simulate the synchrotron, bremsstrahlung and the Compton scattering. Our model roughly reproduces the broad band spectrum of Sgr A* and its variability if we allow for a small fraction of energy to be converted to non-thermal population of electrons. It is also consistent (for a range of viewing angles) with the strong constraints on the amount of circumnuclear material imposed by the measurements of the Faraday rotation.

We study the axisymmetric, time-dependent hydrodynamics of rotating flows that are under the influence of supermassive black hole gravity and radiation from an accretion disk surrounding the black hole. This work is an extension of the earlier work presented by Proga, where nonrotating flows were studied. Here, we consider effects of rotation, a position-dependent radiation temperature, density at large radii, and uniform X-ray background radiation. As in the non-rotating case, the rotating flow settles into a configuration with two components (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the pole. However, with rotation the flow does not always reach a steady state. In addition, rotation reduces the outflow collimation and the outward flux of mass and kinetic energy. Moreover rotation increases the outward flux of the thermal energy and can lead to fragmentation and time-variability of the outflow. We also show that a position-dependent radiation temperature can significantly change the flow solution. In particular, the inflow in the equatorial region can be replaced by a thermally driven outflow. Generally, as it have been discussed and shown in the past, we find that self-consistently determined preheating/cooling from the quasar radiation can significantly reduce the rate at which the central BH is fed with matter. However, our results emphasize also a little appreciated feature. Namely, quasar radiation drives a non-spherical, multi-temperature and very dynamic flow. These effects become dominant for luminosities in excess of 0.01 of the Eddington luminosity.

The collapsar model is the most promising scenario to explain the huge release of energy associated with long duration gamma-ray-bursts (GRBs). Within this scenario GRBs are believed to be powered by accretion through a rotationally support torus or by fast rotation of a compact object. In both cases then, rotation of the progenitor star is one of the key properties because it must be high enough for the torus to form, the compact object to rotate very fast, or both. Here, we check what rotational properties a progenitor star must have in order to sustain torus accretion over relatively long activity periods as observed in most GRBs. We show that simple, often cited, estimates of the total mass available for torus formation and consequently the duration of a GRB are only upper limits. We revise these estimates by taking into account the long term effect that as the compact object accretes the minimum specific angular momentum needed for torus formation increases. This in turn leads to a smaller fraction of the stellar envelope that can form a torus. We demostrate that this effect can lead to a significant, an order of magnidute, reduction of the total energy and overall duration of a GRB event. This of course can be mitigated by assuming that the progenitor star rotates faster then we assumed. However, our assumed rotation is already high compared to observational and theoretical constraints. We also discuss implications of our result.

We present three-dimensional (3-D) hydrodynamical simulations of gas flows in the vicinity of an active galactic nucleus (AGN) powered by a precessing accretion disk. We consider the effects of the radiation force from such a disk on its environment on a relatively large scale (up to ~10 pc. We implicitly include the precessing disk by forcing the disk radiation field to precess around a symmetry axis with a given period ($P$) and a tilt angle ($\Theta$). We study time evolution of the flows irradiated by the disk, and investigate basic dependencies of the flow morphology, mass flux, angular momentum on different combinations of $\Theta$ and $P$. We find the gas flow settles into a configuration with two components, (1) an equatorial inflow and (2) a bipolar inflow/outflow with the outflow leaving the system along the poles (the directions of disk normals). However, the flow does not always reach a steady state. We find that the maximum outflow velocity and the kinetic outflow power at the outer boundary can be reduced significantly with increasing $\Theta$. We also find that of the mass inflow rate across the inner boundary does not change significantly with increasing $\Theta$. (Abbreviated)

Numerical simulations of MHD accretion flows in the vicinity of a supermasssive black hole provide important insights to the problem of why and how systems -- such as the Galactic Center -- are underluminous and variable. To access applicability of such flows to real objects, we examine the dynamical MHD studies with computations of the time dependent radiation spectra predicted by the simulations. We apply Monte Carlo methods to calculate spectra predicted by the time-dependent model of an axisymmetic MHD flow accreting onto a black hole presented by Proga and Begelman. Our calculations show that variability in an accretion flow is not always reflected in the corresponding spectra, at least not in all wavelengths. We find no one-to-one correspondence between the accretion state and the predicted spectrum

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.

I present results from magnetohydrodynamic (MHD) simulations of a gaseous envelope collapsing onto a black hole. These results support the notion that the collapsar model is one of most promising scenarios to explain the huge release of energy in a matter of seconds associated with Gamma Ray Bursts (GRB). Additionally, the MHD simulations show that at late times, when the mass supply rate is expected to decrease, the region in the vicinity of the black hole can play an important role in determining the rate of accretion, its time behaviour, and ultimately the energy output. In particular, the magnetic flux accumulated around the black hole can repeatedly stop and then restart the energy release. As proposed by Proga and Zhang, the episode or episodes of reoccurring of accretion processes can correspond to X-ray flares discovered recently in a number of GRBs.

We present results from multi-epoch spectral analysis of XMM-Newton and Chandra observations of the mini broad absorption line (BAL) quasar PG 1115+080. This is one of the few X-ray detected mini-BAL quasars to date that is bright enough in the X-ray band, mostly due to large gravitational-lensing magnifications, to allow in-depth spectral analysis. The present XMM-Newton observations of PG 1115+080 have provided the highest signal-to-noise X-ray spectra of a mini-BAL quasar obtained to date. By modeling the spectra of PG 1115+080 we have obtained constraints on the column density and ionization state of its outflowing absorbing gas. A comparison between these constraints over several epochs indicates significant variability in the properties of the outflowing absorbers in PG 1115+080. The depths of the high-energy broad absorption features in PG 1115+080 show a significant decrease between the first two observation epochs separated by a rest-frame timescale of ~ 1 year. This variability supports the intrinsic nature of these absorbers. Assuming the interpretation that the high-energy absorption features arise from highly ionized Fe XXV we constrain the fraction of the total bolometric energy released by quasars PG 1115+080 and APM 08279+5225 into the IGM in the form of kinetic energy to be epsilon_k = 0.64(-0.40,+0.52) (68% confidence), and epsilon_k =0.09(-0.05,+0.07), respectively. According to recent theoretical studies this range of efficiencies is large enough to influence significantly the formation of the host galaxy and to regulate the growth of the central black hole.

I present a brief review of theory of winds in active galactic nuclei (AGN). Magnetic, radiation, and thermal driving likely operate in AGN. In many cases, it is difficult to distinguish, both from observational and theoretical point of view, which of these wind driving mechanisms dominates in producing winds. Therefore, I focus on specific theoretical predictions which could help to improve our understanding of the physics of AGN winds.

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.

We present low-resolution ultraviolet spectra of 14 low redshift (z<0.8) quasars observed with HST/STIS as part of a Snap project to understand the relationship between quasar outflows and luminosity. By design, all observations cover the CIV emission line. Nine of the quasars are from the Hamburg-ESO catalog, three are from the Palomar-Green catalog, and one is from the Parkes catalog. The sample contains a few interesting quasars including two broad absorption line (BAL) quasars (HE0143-3535, HE0436-2614), one quasar with a mini-BAL (HE1105-0746), and one quasar with associated narrow absorption (HE0409-5004). These BAL quasars are among the brightest known (though not the most luminous) since they lie at z<0.8. We compare the properties of these BAL quasars to the z<0.5 Palomar-Green and z>1.4 Large Bright Quasar samples. By design, our objects sample luminosities in between these two surveys, and our four absorbed objects are consistent with the v ~ L^0.62 relation derived by Laor & Brandt (2002). Another quasar, HE0441-2826, contains extremely weak emission lines and our spectrum is consistent with a simple power-law continuum. The quasar is radio-loud, but has a steep spectral index and a lobe-dominated morphology, which argues against it being a blazar. The unusual spectrum of this quasar resembles the spectra of the quasars PG1407+265, SDSSJ1136+0242, and PKS1004+13 for which several possible explanations have been entertained.

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.

We present a Chandra-LETGS observation of the Seyfert 1 galaxy Mrk 279. This observation was carried out simultaneously with HST-STIS and FUSE, in the context of a multiwavelength study of this source. The Chandra pointings were spread over ten days for a total exposure time of ~360 ks. The spectrum of Mrk279 shows evidence of broad emission features, especially at the wavelength of the OVII triplet. We quantitatively explore the possibility that this emission is produced in the broad line region (BLR). We modeled the broad UV emission lines seen in the FUSE and HST-STIS spectra following the ``locally optimally emitting cloud" approach. We find that the X-ray lines luminosity derived from the best fit BLR model can match the X-ray features, suggesting that the gas producing the UV lines is sufficient to account also for the X-ray emission. The spectrum is absorbed by ionized gas whose total column density is ~5x10^{20} cm^{-2}. The absorption spectrum can be modeled by two distinct gas components (log xi ~ 0.47 and 2.49, respectively) both showing a significant outflow velocity. However, the data allow also the presence of intermediate ionization components. The distribution of the column densities of such extra components as a function of the ionization parameter is not consistent with a continuous, power law-like, absorber, suggesting a complex structure for the gas outflow for Mrk 279 (abridged).

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.

We present our recent results from numerical simulations of a magnetized flow in the vicinity of a black hole in the context of the collapsar model for GRBs. The simulations show that after an initial transient, the flow settles into a complex convolution of several distinct, time-dependent flow components including an accretion torus, its corona and outflow, an inflow and an outflow in the polar funnel. We focus on studying the nature and connection between these components, in particular between the inflows and related outflows. We find that rotational and MHD effects launch, accelerate, and sustain the outflows. We also find that an outflow can be formed even when the collapsing envelope has initially a very weak magnetic field and a very small angular momentum. Our main conclusion is that even for a relatively weak initial magnetic field and a slow rotation, a gravitational collapse of a stellar envelope can lead to formation of a very strong and very fast jet.

Soft X-ray flares were detected to follow the short-duration gamma-ray burst GRB 050724. The temporal properties of the flares suggest that they are likely due to the late time activity of the central engine. We argue that if short GRBs are generated through compact star mergers, as is supported by the recent observations, the jet powering the late X-ray flares must be launched via magnetic processes rather than via neutrino-antineutrino annihilations. As a result, the X-ray flares following short GRBs are expected to be linearly polarized. The argument may also apply to the X-ray flares following long GRBs. Future observations with the upcoming X-ray polarimeters will test this prediction.

Within the standard accretion disk theory for active galactic nuclei (AGN), the observed X-rays are often modeled by Compton up-scattering of ultraviolet (UV) disk photons inside a hot disk corona. Here, we point out that for many AGN, radiation pressure due to the very same UV disk photons can drive a flow from the disk into the corona and couple the processes producing X-rays and UV photons. This coupling could lead to quenching of the disk corona because the regions above the UV disk will be too dense, too opaque, and consequently too cold. We discuss various consequences of this new type of the X-ray/UV coupling on the dynamical and radiative properties of AGN.

We present results from axisymmetric, time-dependent hydrodynamical (HD) and magnetohydrodynamical (MHD) simulations of a gaseous envelope collapsing onto a black hole (BH). We consider gas with so small angular momentum that after an initial transient, the flow in the HD case, accretes directly onto a BH without forming a rotationally support torus. However, in the MHD case even with a very weak initial magnetic field, the flow settles into a configuration with four components: (i) an equatorial inflow, (ii) a bipolar outflow, (iii) polar funnel outflow, and (iv) polar funnel inflow. We focus our analysis on the second flow component of the MHD flow which represents a simple yet robust example of a well-organized inflow/outflow solution to the problem of MHD jet formation. The jet is heavy, highly magnetized, and driven by magnetic and centrifugal forces. A significant fraction of the total energy in the jet is carried out by a large scale magnetic field. We review previous simulations, where specific angular momentum was higher than that assumed here, and conclude that our bipolar outflow develops for a wide range of the properties of the flow near the equator and near the poles. Future work on such a simple inflow/outflow solution will help to pinpoint the key elements of real jets/outflows as well as help to interpret much more complex simulations aimed at studying jet formation and collapse of magnetized envelopes.

We present an analysis of the intrinsic UV absorption in the Seyfert 1 galaxy Mrk 279 based on simultaneous long observations with the Hubble Space Telescope (41 ks) and the Far Ultraviolet Spectroscopic Explorer (91 ks). To extract the line-of-sight covering factors and ionic column densities, we separately fit two groups of absorption lines: the Lyman series and the CNO lithium-like doublets. For the CNO doublets we assume that all three ions share the same covering factors. The fitting method applied here overcomes some limitations of the traditional method using individual doublet pairs; it allows for the treatment of more complex, physically realistic scenarios for the absorption-emission geometry and eliminates systematic errors that we show are introduced by spectral noise. We derive velocity-dependent solutions based on two models of geometrical covering -- a single covering factor for all background emission sources, and separate covering factors for the continuum and emission lines. Although both models give good statistical fits to the observed absorption, we favor the model with two covering factors because: (a) the best-fit covering factors for both emission sources are similar for the independent Lyman series and CNO doublet fits; (b) the fits are consistent with full coverage of the continuum source and partial coverage of the emission lines by the absorbers, as expected from the relative sizes of the nuclear emission components; and (c) it provides a natural explanation for variability in the Ly$\alpha$ absorption detected in an earlier epoch. We also explore physical and geometrical constraints on the outflow from these results.

We review the main results from line-driven (LD) models of winds from accretion disks in cataclysmic variables (CVs). We consider LD disk wind models in the hydrodynamic (HD) and magnetohydrodynamic (MHD limits. We discuss the basic physical conditions needed for a disk wind to exist and the conditions for the wind to be steady or unsteady. We also discuss how the line-driven wind structures revealed in numerical simulations relate to observations. In particular, we present synthetic line profiles predicted by the LD wind models and compare them with observations. Our main conclusion is that, despite some problems, line-driving alone is the most plausible mechanism for driving the CV winds. This conclusion is related to two facts: 1) LD wind models are likely the best studied wind models, analytically and numerically, and 2) it is most likely that the predictive power of the LD wind models is much higher than of any other wind model so far. Preliminary results from LD-MHD wind models confirm that magnetic driving is likely an important element of the wind dynamics. However, magnetic driving does not seem to be necessary to produce a CV wind. The most important issues which need to be addressed by future dynamical models, regardless of driving mechanism, are the effects of the position-dependent photoionization and the dynamical effects in three dimensions.

We present results from our magnetohydrodynamical simulations of accretion flows onto black holes. Our main focus is the interplay between inflows and related outflows. We consider applications of such flows to the Galatic center and low luminosity active galactic nuclei.

We explore consequences of a radiation driven disk wind model for mass outflows from active galactic nuclei (AGN). We performed axisymmetric time-dependent hydrodynamic calculations using the same computational technique as Proga, Stone and Kallman (2000). We test the robustness of radiation launching and acceleration of the wind for relatively unfavorable conditions. In particular, we take into account the central engine radiation as a source of ionizing photons but neglect its contribution to the radiation force. Additionally, we account for the attenuation of the X-ray radiation by computing the X-ray optical depth in the radial direction assuming that only electron scattering contributes to the opacity. Our new simulations confirm the main result from our previous work: the disk atmosphere can 'shield' itself from external X-rays so that the local disk radiation can launch gas off the disk photosphere. We also find that the local disk force suffices to accelerate the disk wind to high velocities in the radial direction. This is true provided the wind does not change significantly the geometry of the disk radiation by continuum scattering and absorption processes; we discuss plausibility of this requirement. Synthetic profiles of a typical resonance ultraviolet line predicted by our models are consistent with observations of broad absorption line (BAL) QSOs.

We investigate the applicability of inhomogeneous absorber models in the formation of AGN outflow absorption-troughs. The models we explore are limited to monotonic gradients of absorbing column densities in front of a finite emission source. Our main finding is that simple power-law and gaussian distributions are hard pressed to fit the Mrk 279 high-quality UV outflow data. An acceptable fit for the O VI troughs can only be obtained by assuming unrealistic optical depth values (upward of 100). The strongest constraints arise from the attempt to fit the Lyman series troughs. In this case it is evident that even allowing for complete freedom of both the power-law exponent and the optical depth as a function of velocity cannot yield an acceptable fit. In contrast, partial covering models do yield good fits for the Lyman series troughs. We conclude that monotonic inhomogeneous absorber models that do not include a sharp edge in the optical depth distribution across the source are not an adequate physical model to explain the trough formation mechanism for the outflow observed in Mrk 279.

One of the main problems in modeling the ionised outflows in Active Galactic Nuclei is the unknown distance of the outflowing wind to the central source. Only if the density is known this distance can be determined through the ionisation parameter. Here we study density diagnostics based upon O V transitions. O V is known to have metastable levels that are density dependent. We study the population of those levels under photoionisation equilibrium conditions and determine for which parameter range they can have a significant population. We find that resonance line trapping plays an important role in reducing the critical densities above which the metastable population becomes important. We investigate the K-shell absorption lines from these metastable levels. Provided that there is a sufficient population of the metastable levels, the corresponding K-shell absorption lines are detectable and are well separated from the main absorption line originating from the ground state. We then present the Chandra LETGS spectrum of the Seyfert 1 galaxy Mrk 279 that may show for the first time the presence of these metastable level absorption lines. A firm identification is not yet possible due to both uncertainties in the observed wavelength of the strongest line as well as uncertainties in the predicted wavelength. If the line is indeed due to absorption from O V, then we deduce a distance to the central source of one light week to a few light months, depending upon the importance of additional heating processes.

We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.

We summarize recent developments in modeling of mass outflows from accretion disks in AGN. We illustrate how the hydrodynamical and magnetohydrodynamical wind structures revealed in numerical simulations relate to observations. Important results from the numerical simulations include: no shielding from X-rays is needed to radiatively launch a wind from the inner disk; shielding is needed to accelerate the wind to high velocities. Disk winds launched magnetically seem to be too dense to be accelerated by line driving. It appears that we have an all or nothing situation for radiation-driving: radiation driving has to be responsible for both launching and accelerating the wind or a wind launched by some other mechanism must have a very finely tuned density in order to be accelerated by radiation driving to velocities as high as measured in BAL QSOs.

We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. We begin the simulations after the 1.7 MSUN iron core of a 25 MSUN presupernova star has collapsed and study the ensuing accretion of the 7 MSUN helium envelope onto the central black hole formed by the collapsed iron core. We consider a spherically symmetric progenitor model, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum and a weak radial magnetic field. Our MHD simulations include a realistic equation of state, neutrino cooling, photodisintegration of helium, and resistive heating. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.

We propose a gamma-ray burst scenario involving relativistic jets dominated by Poynting flux with alternating toroidal magnetic fields. Such a structure may arise naturally if the jet is formed and powered by the accretion flow in the core of a collapsar. We conjecture that the polarity of the toroidal magnetic field changes randomly due to hydromagnetic turbulence driven by the magneto-rotational instability (MRI), with the typical reversal time determined by the time scale for amplifying magnetic fields up to a dynamically important level. Poynting flux-dominated jets with reversing B-fields provide a natural and efficient way to dissipate energy via the reconnection process. Gamma-rays are produced at the spatially separated reconnection sites. In this scenario, the emergent synchrotron radiation can be highly polarized and can form both smooth and spiky light curves. We note the possibility that cold and dense filaments can form in the reconnection zones as the result of thermal instability. One could then explain the production of very hard X-ray spectra, as due to bound-free absorption of the synchrotron radiation.

We report on resonance line profiles predicted by radiation driven disk wind models which extend radially one order of magnitude farther out than our previous models. Our main result is that the inclusion of a disk wind at larger radii changes qualitatively and quantitatively the line profiles predicted by the models. Our new models predict line absorption that is significantly stronger than those predicted by old models. Some of the previous line profiles exhibit a doubled-humped structure near the line center which is now replaced by a more plausible single, blueshifted minimum. We emphasize that the improvements in the shape as well as the strength of the absorption were achieved without changing the gross properties of the wind. In particular, our new models do not predict a higher mass-loss rate than the previous models. The main changes in the line profiles are due to the fact that the ratio between the rotational velocity and poloidal velocity of the wind decreases downstream. The new line profiles reproduce well the line absorption of the nova-like variable, IX~Vel, and promise to reproduce observations of other cataclysmic variables. This success of the radiation driven disk wind model provides an important link between outflows in OB stars and outflows in active galactic nuclei.

We report on the second phase of our study of slightly rotating accretion flows onto black holes. We consider magnetohydrodynamical (MHD) accretion flows with a spherically symmetric density distribution at the outer boundary, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum and a weak radial magnetic field. We study accretion flows by means of numerical 2D, axisymmetric, MHD simulations with and without resistive heating. Our main result is that the properties of the accretion flow depend mostly on an equatorial accretion torus which is made of the material that has too much angular momentum to be accreted directly. The torus accretes, however, because of the transport of angular momentum due to the magnetorotational instability (MRI). Initially, accretion is dominated by the polar funnel, as in the hydrodynamic inviscid case, where material has zero or very low angular momentum. At the later phase of the evolution, the torus thickens towards the poles and develops a corona or an outflow or both. Consequently, the mass accretion through the funnel is stopped. The accretion of rotating gas through the torus is significantly reduced compared to the accretion of non-rotating gas (i.e., the Bondi rate). It is also much smaller than the accretion rate in the inviscid, weakly rotating case.Our results do not change if we switch on or off resistive heating. Overall our simulations are very similar to those presented by Stone, Pringle, Hawley and Balbus despite different initial and outer boundary conditions. Thus, we confirm that MRI is very robust and controls the nature of radiatively inefficient accretion flows.

We study the two-dimensional, time-dependent MHD of radiation-driven winds from luminous accretion disks initially threaded by a purely axial magnetic field. The radiation force is mediated primarily by spectral lines. We use ideal MHD to compute the evolution of Keplerian disks, varying the magnetic field strengths and the luminosity of the disk, the central accreting object or both. We find that the magnetic fields very quickly start deviating from purely axial due to the magnetorotational instability. This leads to fast growth of the toroidal magnetic field. As a result the toroidal field dominates over the poloidal field above the disk and the gradient of the former drives a slow and dense disk outflow, which conserves specific angular momentum. Depending on the strength of the magnetic field relative to the system luminosity the disk wind can be radiation- or MHD driven. The pure radiation-driven wind consists of a dense, slow outflow that is bounded on the polar side by a high-velocity stream. The mass-loss rate is mostly due to the fast stream. As the magnetic field strength increases first the slow part of the flow is affected, namely it becomes denser and slightly faster and begins to dominate the mass-loss rate. In very strong magnetic field or pure MHD cases, the wind consists of only a dense, slow outflow without the presence of the distinctive fast stream so typical to pure radiation-driven winds. Our simulations indicate that winds launched by the magnetic fields are likely to remain dominated by the fields downstream because of their relatively high densities.

We report on the first phase of our study of slightly rotating accretion flows onto black holes. We consider inviscid accretion flows with a spherically symmetric density distribution at the outer boundary, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum. We study accretion flows by means of numerical 2D, axisymmetric, hydrodynamical simulations. Our main result is that the properties of the accretion flow do not depend as much on the outer boundary conditions (i.e., the amount as well as distribution of the angular momentum) as on the geometry of the non-accreting matter. The material that has too much angular momentum to be accreted forms a thick torus near the equator. Consequently, the geometry of the polar region, where material is accreted (the funnel), and the mass accretion rate through it are constrained by the size and shape of the torus. Our results show one way in which the mass accretion rate of slightly rotating gas can be significantly reduced compared to the accretion of non-rotating gas (i.e., the Bondi rate), and set the stage for calculations that will take into account the transport of angular momentum and energy.

We present synthetic line profiles as predicted by the models of 2-D line- driven disk winds due to Proga, Stone & Drew. We compare the model line profiles with HST observations of the cataclysmic variable IX Vel. The model wind consists of a slow outflow that is bounded on the polar side by a fast stream. We find that these two components of the wind produce distinct spectral features. The fast stream produces profiles which show features consistent with observations. These include the appearance of the P-Cygni shape for a range of inclinations, the location of the maximum depth of the absorption component at velocities less than the terminal velocity, and the transition from absorption to emission with increasing inclination. However the model profiles have too little absorption or emission equivalent width. This quantitative difference between our models and observations is not a surprise because the line-driven wind models predict a mass loss rate that is lower than the rate required by the observations. We note that the model profiles exhibit a double-humped structure near the line center which is not echoed in observations. We identify this structure with a non-negligible redshifted absorption which is formed in the slow component of the wind where the rotational velocity dominates over expansion velocity. We conclude that the next generation of disk wind models, developed for application to CVs, needs to yield stronger wind driving out to larger disk radii than do the present models.

To confront the predictions of the most recent line-driven disk wind models with observational evidence, we have obtained HST STIS (1180--1700A) echelle spectra of the nova-like variables IX Vel and V3885 Sgr at three epochs. The targets were observed in timetag mode for ~2000sec on each occasion, allowing us to study the spectral time evolution on timescales down to ~10sec. The mean UV spectra are characterised by the wind signature of broad blueshifted absorption in Lya, NV1240, SiIV1398, CIV1549 and HeII1640. There is very little redshifted emission other than in CIV. Narrow blueshifted absorption dips, superposed on the broad absorption at around -900km/s, accompany periods of well-developed wind activity. The continuum level and mean line profiles vary markedly from observation to observation -- with the wind signatures almost disappearing in one epoch of observation of IX Vel. The strong positive correlation between UV brightness and wind activity predicted by line-driven disk wind models is disobeyed by both binaries. The wind signatures in IX Vel's UV spectrum are revealed to be remarkably steady on timescales ranging from ~10 to ~1000sec. More variability is seen in V3885 Sgr, the binary with the lower opacity outflow. But there is only one epoch in which the line profile changes significantly in ~100sec or less. Narrow absorption dips, when present, show only smooth, small changes in velocity.

Far ultraviolet (FUV) observations are excellent probes of the inner accretion disk, disk outflows, and the mass-accumulating white dwarf in cataclysmic variables. Here we study the contrasting behavior of two canonical dwarf novae in outburst by presenting FUSE FUV (904 -- 1187 A) spectroscopy of U Gem and SS Cyg. We observed each system four times during a single outburst. The outburst peak and early decline spectra of SS Cyg are well fit by models of a steady-state accretion disk and a biconical wind. A broad, blueshifted OVI wind-formed absorption line is the only strong spectral feature. In late outburst decline, OVI and CIII lines are seen as broad emission features and the continuum has flattened. In U Gem, the continua of the optical outburst plateau spectra are plausibly fit by accretion disk model spectra. The spectra also show numerous narrow, low-velocity absorption lines that do not originate in the inner accretion disk. We discuss the line spectra in the context of partial absorption of the FUV continuum by low-velocity, vertically-extended material located at large disk radii. The late outburst decline spectrum of U Gem is dominated by the white dwarf. WD model fits confirm the sub-solar C and super-solar N abundances found in earlier studies.

I review the main results from recent 2-D, time-dependent hydrodynamic models of radiation-driven winds from accretion disks in AGN. I also discuss the physical conditions needed for a disk wind to be shielded from the strong X-rays and to be accelerated to hypersonic velocities. I conclude with a few remarks on winds in hot stars, low mass young stellar objects, cataclysmic variables, low mass X-ray binaries, and galactic black holes and future work.

X-ray heating of the photosphere of an accretion disk is a possible mechanism to produce strong, broad UV emission lines in low mass X-ray binaries (LMXBs). However, detailed photoionization calculations show that this mechanism fails to produce sufficient emission measure. We present the results of hydrodynamical calculations of the disk photosphere irradiated by strong X-rays. We attempt to determine whether LMXBs can harbor significant UV-driven disk winds despite the effects of X-ray ionization. Such winds would be a likely candidate for the site of emission of UV lines and may better explain the observations than the X-ray heated disk photosphere. We find that the local disk radiation cannot launch a wind from the disk because of strong ionizing radiation from the central object. Unphysically high X-ray opacities would be required to shield the UV emitting disk and allow the line force to drive a disk wind. However the same X-ray radiation that inhibits line driving heats the disk and can produce a hot bipolar wind or corona above the disk. To assess the impact of X-ray heating upon driving of a disk wind by the line force in any system with an accretion disk we derive analytic formulae. In particular, we compare results of line-driven disk wind models for accretion disks in LMXBs and active galactic nuclei. The latter show spectral features associated with a strong and fast disk wind. The key parameter determining the role of the line force is not merely the presence of the luminous UV zone in the disk and the presence of the X-rays, but also the distance of this UV zone from the center.

We have obtained FUV (904 -- 1187 \AA) spectra of U Gem in outburst with FUSE. Three of the observations were acquired during the plateau phase of the outburst, while the fourth was obtained during late outburst decline. The plateau spectra have continuum shapes and fluxes that are approximated by steady-state accretion disk model spectra with $\dot{m} \simeq$ 7e-09 M_sun/yr. The spectra also show numerous absorption lines of H I, He II, and two- to five-times ionized transitions of C, N, O, P, S and Si. There are no emission features in the spectra, with the possible exception of a weak feature on the red wing of the O VI doublet. The absorption lines are narrow (FWHM ~ 500 km/s), too narrow to arise from the disk photosphere, and at low velocities (<= 700 km/s). The S VI and O VI doublets are optically thick. The absorption lines in the plateau spectra show orbital variability: in spectra obtained at orbital phases 0.53 -- 0.79, low-ionization absorption lines appear and the central depths of the pre-existing lines increase. The increase in line absorption occurs at the same orbital phases as previously observed EUV and X-ray light curve dips. If the absorbing material is in (near-) Keplerian rotation around the disk, it must be located at large disk radii. The final observation occurred when U Gem was about two magnitudes from optical quiescence. The spectra are dominated by emission from an $\simeq$43,000 K, metal-enriched white dwarf (WD). The inferred radius of the WD is 4.95e+08 cm, close to that observed in quiescence. Allowing for a hot heated region on the surface of the WD improves the fit to the spectrum at short (< 960 \AA) wavelengths.

We analyze optical and ultraviolet observations of the symbiotic binary AG Pegasi acquired during 1992-97. The bolometric luminosity of the hot component declined by a factor of 2-3 from 1980-1985 to 1997. Since 1992, the effective temperature of the hot component may have declined by 10%-20%, but this decline is comparable to the measurement errors. Optical observations of H-beta and He I emission show a clear illumination effect, where high energy photons from the hot component ionize the outer atmosphere of the red giant. Simple illumination models generally account for the magnitude of the optical and ultraviolet emission line fluxes. High ionization emission lines - [Ne V], [Mg V], and [Fe VII] - suggest mechanical heating in the outer portions of the photoionized red giant wind. This emission probably originates in a low density region $\sim$ 30-300 AU from the central binary.

We present the results of time-dependent hydrodynamic calculations of line- driven winds from accretion disks in AGN. To calculate the radiation force, we take into account radiation from the disk and the central engine. The gas temperature and ionization state in the wind are calculated self-consistently. We find that a disk accreting onto a 10^8 MSUN black hole at the rate of 1.8 MSUN/YR can launch a wind at about 10^16 cm from the central engine. The X-rays from the central object are attenuated by the disk atmosphere so they cannot prevent the local disk radiation from pushing matter away from the disk. For a reasonable X-ray opacity, the disk wind can be accelerated by the central UV radiation to velocities of up to 15000 km/s at a distance of about 10^17 cm from the central engine. The covering factor of the disk wind is $\sim 0.2$. The wind is unsteady and consists of an opaque, slow vertical flow near the disk that is bounded on the polar side by a high-velocity stream. A typical column density through the fast stream is a few 10^23 cm^-2. The fast stream contributes nearly 100% to the total wind mass loss rate of 0.5 MSUN/YR.

We study the 2-D, time-dependent hydrodynamics of radiation-driven winds from luminous accretion disks threaded by a strong, large-scale, ordered magnetic field. The radiation force is due to spectral lines and is calculated using a generalized multidimensional formulation of the Sobolev approximation. The effects of the magnetic field are approximated by adding a force that emulates a magnetocentrifugal force. Our approach allows us to calculate disk winds when the magnetic field controls the flow geometry, forces the flow to corotate with the disk, or both. In particular, we calculate models where the lines of the poloidal component of the field are straight and inclined to the disk at a fixed angle. Our numerical calculations show that flows which corotate with the disk have a larger mass loss rate than their counterparts which conserve specific angular momentum. The difference in the mass loss rate between these two types of winds can be several orders of magnitude for low disk luminosities but vanishes for high disk luminosities. Winds which corotate with the disk have much higher velocities than angular momentum conserving winds. Fixing the wind geometry stabilizes winds which are unsteady when the geometry is derived self-consistently. The inclination angle between the poloidal velocity and the normal to the disk midplane is important. Non-zero inclination angles allow the magnetocentrifugal force to increase the mass loss rate for low luminosities, and increase the wind velocity for all luminosities. Our calculations also show that the radiation force can launch winds from magnetized disks. The line force can be essential in producing MHD winds from disks where the thermal energy is too low to launch winds or where the field lines make an angle of < 30^o with respect to the normal to the disk.

We describe an efficient method of calculating the radiation pressure due to spectral lines, including all the terms in the velocity gradient tensor. We apply this method to calculate the two-dimensional, time-dependent structure of winds from luminous disks. Qualitative features of our new models are very similar to those we calculated including only the dominant terms in the tensor (Proga, Stone & Drew 1998, hereafter PSD). In particular, we find that models which displayed unsteady behaviour in PSD are also unsteady with the new method, and gross properties of the winds, such as mass-loss rate and characteristic velocity are not changed by the more accurate approach. The largest change caused by the new method is in the disk-wind opening angle: winds driven only by the disk radiation are more polar with the new method while winds driven by the disk and central object radiation are typically more equatorial. In the closing discussion, we provide further insight into the way the geometry of the radiation field and consequent flow determines the time properties of the flow.

We compare models of line-driven winds from accretion discs and single spherical stars. We look at the problem of scaling mass-loss rates and velocities of stellar and disc winds with model parameters. We find that stellar and disc winds driven by radiation, within the CAK framework, are very similar as far as mass-loss rates and velocities are concerned. Thus we can use analytic results for stellar winds to rescale, in a first order approximation, numerical results for disc winds. We also show how the CAK stellar solutions change when we take into account effects of very low luminosities or line-driving force.

New high resolution spectroscopic and medium resolution spectropolarimetric data, complemented with optical broad and narrow band imaging, of the B[e] star HD 87643 are presented. The spectrum of HD 87643 exhibits the hybrid characteristics well known to be representative of the group of B[e] stars; a fast wind with an expansion velocity in excess of 1000 km/s is measured in the hydrogen and helium lines, while a slower component is traced by lower excitation lines and forbidden lines. Clues to the geometry of the rapidly expanding circumstellar shell are provided by the startling polarization changes across Halpha. Comparison with published schematic calculations indicates that the polarizing material is located in a slowly rotating, expanding disk structure. A hydrodynamical model is then presented whose results are consistent with the original two-wind concept for B[e] stars and exhibits kinematic properties that may well explain the observed spectral features in HD 87643. The model calculations use as input a B star undergoing mass loss, surrounded by an optically thick disk. The resulting configuration consists of a fast polar wind from the star and a slowly expanding disk wind. The model also predicts that the stellar wind at intermediate latitudes is slower and denser than in the polar region.

A radiation-driven disk wind model is proposed that offers great promise of explaining the extreme mass loss signatures of massive young stellar objects (the BN-type objects and more luminous Herbig Be stars). It is argued that the dense low-velocity winds associated with young late-O/early-B stars would be the consequence of continuing optically-thick accretion onto them. The launch of outflow from a Keplerian disk allows wind speeds of about 200 km/s that are substantially less than the escape speed from the stellar surface. The star itself is not required to be a rapid rotator. Disk irradiation is taken into account in the hydrodynamical calculation presented, and identified as an important issue both observationally and from the dynamical point of view.

We describe a non-LTE photoionization code to calculate the wind structure and emergent spectrum of a red giant wind illuminated by the hot component of a symbiotic binary system. We consider spherically symmetric winds with several different velocity and temperature laws and derive predicted line fluxes as a function of the red giant mass loss rate, \mdot. Our models generally match observations of the symbiotic stars EG And and AG Peg for \mdot about 10^{-8} \msunyr to 10^{-7} \msunyr. The optically thick cross- section of the red giant wind as viewed from the hot component is a crucial parameter in these models. Winds with cross-sections of 2--3 red giant radii reproduce the observed fluxes, because the wind density is then high, about 10^9 cm^{-3}. Our models favor winds with acceleration regions that either lie far from the red giant photosphere or extend for 2--3 red giant radii.

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.