No Arabic abstract
We explore the effect of momentum-driven winds representing radiation pressure driven outflows from accretion onto supermassive black holes in a set of numerical hydrodynamical simulations. We explore two matched sets of cosmological zoom-in runs of 24 halos with masses ~$10^{12.0}-10^{13.4}$ M_sun run with two different feedback models. Our `NoAGN model includes stellar feedback via UV heating, stellar winds and supernovae, photoelectric heating and cosmic X-ray background heating from a meta-galactic background. Our fiducial `MrAGN model is identical except that it also includes a model for black hole seeding and accretion, as well as heating and momentum injection associated with the radiation from black hole accretion. Our MrAGN model launches galactic outflows which result in both `ejective feedback - the outflows themselves which drive gas out of galaxies - and `preventative feedback, which suppresses the inflow of new and recycling gas. As much as 80 % of outflowing galactic gas can be expelled, and accretion can be suppressed by as much as a factor of 30 in the MrAGN runs when compared with the NoAGN runs. The histories of NoAGN galaxies are recycling-dominated, with ~70% of material that leaves the galaxy eventually returning, and the majority of outflowing gas re-accretes on 1 Gyr timescales without AGN feedback. Outflowing gas in the MrAGN runs has higher characteristic velocity (500 - 1,000 km/s versus 100-300 km/s for outflowing NoAGN gas) and travels as far as a few Mpcs. Only ~10% of ejected material is re-accreted in the MrAGN galaxies.
We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $rsim R_{rm B}$, where the density follows $rho propto (1+R_{rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $rho propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($dot{M}propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $dot{M}_{rm B}$ The net accretion rate depends on the viscous strength, following $dot{M}/dot{M}_{rm B}propto (alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $dot{M}_{rm B}/dot{M}_{rm Edd}< 10^{-3}$ having minimal radiation cooling, where $dot{M}_{rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $dot{M}/dot{M}_{rm B} sim 10^{-3}-10^{-2}$. This rate is as low as $dot{M}/dot{M}_{rm Edd} sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.
We want to test if self-similar magneto-hydrodynamic (MHD) accretion-ejection models can explain the observational results for accretion disk winds in BHBs. In our models, the density at the base of the outflow, from the accretion disk, is not a free parameter, but is determined by solving the full set of dynamical MHD equations without neglecting any physical term. Different MHD solutions were generated for different values of (a) the disk aspect ratio ($varepsilon$) and (b) the ejection efficiency ($p$). We generated two kinds of MHD solutions depending on the absence (cold solution) or presence (warm solution) of heating at the disk surface. The cold MHD solutions are found to be inadequate to account for winds due to their low ejection efficiency. The warm solutions can have sufficiently high values of $p (gtrsim 0.1)$ which is required to explain the observed physical quantities in the wind. The heating (required at the disk surface for the warm solutions) could be due to the illumination which would be more efficient in the Soft state. We found that in the Hard state a range of ionisation parameter is thermodynamically unstable, which makes it impossible to have any wind at all, in the Hard state. Our results would suggest that a thermo-magnetic process is required to explain winds in BHBs.
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.
Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a subgrid model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($lsim 10^{11.5} msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations.
We explore the poloidal structure of two-dimensional (2D) MHD winds in relation to their potential association with the X-ray warm absorbers (WAs) and the highly-ionized ultra-fast outflows (UFOs) in AGN, in a single unifying approach. We present the density $n(r,theta)$, ionization parameter $xi(r,theta)$, and velocity structure $v(r,theta)$ of such ionized winds for typical values of their fluid-to-magnetic flux ratio, $F$, and specific angular momentum, $H$, for which wind solutions become super-Alfvenic. We explore the geometrical shape of winds for different values of these parameters and delineate the values that produce the widest and narrowest opening angles of these winds, quantities necessary in the determination of the statistics of AGN obscuration. We find that winds with smaller $H$ show a poloidal geometry of narrower opening angles with their Alfven surface at lower inclination angles and therefore they produce the highest line of sight (LoS) velocities for observers at higher latitudes with the respect to the disk plane. We further note a physical and spatial correlation between the X-ray WAs and UFOs that form along the same LoS to the observer but at different radii, $r$, and distinct values of $n$, $xi$ and $v$ consistent with the latest spectroscopic data of radio-quiet Seyfert galaxies. We also show that, at least in the case of 3C 111, the winds pressure is sufficient to contain the relativistic plasma responsible for its radio emission. Stratified MHD disk-winds could therefore serve as a unique means to understand and unify the diverse AGN outflows.