BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Yan-Fei Jiang

Yan-Fei Jiang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-ph.SRastro-ph.IMastro-ph.COphysics.comp-phastro-ph.EPphysics.plasm-ph

Catalog footprint

What is connected

31works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

Cosmic-Ray Transport in Varying Galactic Environments

We study the propagation of mildly-relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta+21 to post-process three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfvén speed, and diffusion relative to the Alfvén waves, with the diffusion coefficient set by the balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the post-processed simulations to make predictions regarding potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extra-planar region (|z|>0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CR Driven Multi-phase Gas Formed via Thermal Instability

Cosmic rays (CRs) are an important energy source in the circum-galactic medium (CGM) that impact the multi-phase gas structure and dynamics. We perform two-dimensional CR-magnetohydrodynamic simulations to investigate the role of CRs in accelerating multi-phase gas formed via thermal instability. We compare outflows driven by CRs to those driven by a hot wind with equivalent momentum. We find that CRs driven outflow produces lower density contrast between cold and hot gas due to non-thermal pressure support, and yields a more filamentary cloud morphology. While entrainment in a hot wind can lead to cold gas increasing due to efficient cooling, CRs tend to suppress cold gas growth. The mechanism of this suppression depends on magnetic field strength, with CRs either reducing cooling or shredding the clouds by differential acceleration. Despite the suppression of cold gas growth, CRs are able to launch the cold clouds to observed velocities without rapid destruction. The dynamical interaction between CRs ad multi-phase gas is also sensitive to the magnetic field strength. In relatively strong fields, the CRs are more important for direct momentum input to cold gas. In relatively weak fields, the CRs impact gas primarily by heating, which modifies gas pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multi-group Radiation Magneto-hydrodynamics based on Discrete Ordinates including Compton Scattering

We present a formulation and numerical algorithm to extend the scheme for grey radiation magneto-hydrodynamics (MHD) developed by Jiang (2021) to include the frequency dependence via the multi-group approach. The entire frequency space can be divided into arbitrary number of groups in the lab frame, and we follow the time dependent evolution of frequency integrated specific intensities along discrete rays inside each group. Spatial transport of photons is done in the lab frame while all the coupling terms are solved in the fluid rest frame. Lorentz transformation is used to connect different frames. Radiation transport equation is solved fully implicitly in time while the MHD equations are evolved explicitly so that time step is not limited by the speed of light. A finite volume approach is used for transport in both spatial and frequency spaces to conserve radiation energy density and momentum. The algorithm includes photon absorption, electron scattering as well as Compton scattering, which is calculated by solving the Kompaneets equation. The algorithm is accurate for a wide range of optical depth conditions and can handle both radiation pressure and gas pressure dominated flows. It works for both Cartesian and curvilinear coordinate systems with adaptive mesh refinement. We provide a variety of test problems including radiating sphere, shadow test, absorption of a moving gas, Bondi type flows as well as a collection of test problems for thermal and bulk Compton scattering. We also discuss examples where frequency dependence can make a big difference compared with the grey approach.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical Simulations of Convective 3-Dimensional Red Supergiant Envelopes

We explore the three-dimensional properties of convective, luminous ($L\approx10^{4.5}-10^{5}L_\odot$), Hydrogen-rich envelopes of Red Supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry using $\texttt{Athena++}$. These computations comprise $\approx30\%$ of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer $\approx 3M_\odot$ of the simulated $M\approx15M_\odot$ stars. This work reveals a radius, $R_\mathrm{corr}$, around which the nature of the convection changes. For $r>R_\mathrm{corr}$, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well-studied in less luminous stars, but in RSGs, the near- (or above-) Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outwards rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This "halo" of material will impact predictions for both shock breakout and early lightcurves of Type II-P supernovae. Inside of $R_\mathrm{corr}$, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length-theory (MLT). Radiation pressure provides $\approx1/3$ of the support against gravity in this region. Our comparisons to MLT suggest a mixing length of $α=3-4$, consistent with the sizes of convective plumes seen in the simulations. The temporal variability of these 3D models is mostly on the timescale of the convective plume lifetimes ($\approx300$ days), with amplitudes consistent with those observed photometrically.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Radiative Relativistic Magnetohydrodynamic Simulations of Neutron Star Column Accretion in Cartesian Geometry

High luminosity accretion onto a strongly magnetized neutron star results in a radiation pressure dominated, magnetically confined accretion column. We investigate the dynamics of these columns using two-dimensional radiative relativistic magnetohydrodynamic simulations, restricting consideration to modest accretion rates where the height of the column is low enough that Cartesian geometry can be employed. The column structure is dynamically maintained through high-frequency oscillations of the accretion shock at $\simeq 10-25$~kHz. These oscillations arise because it is necessary to redistribute the power released at the accretion shock through bulk vertical motions, both to balance the cooling and to provide vertical pressure support against gravity. Sideways cooling always dominates the loss of internal energy. In addition to the vertical oscillations, photon bubbles form in our simulations and add additional spatial complexity to the column structure. They are not themselves responsible for the oscillations, and they do not appear to affect the oscillation period. However, they enhance the vertical transport of radiation and increase the oscillation amplitude in luminosity. The time-averaged column structure in our simulations resembles the trends in standard 1D stationary models, the main difference being that the time-averaged height of the shock front is lower because of the higher cooling efficiency of the 2D column shape.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shock Breakout in 3-Dimensional Red Supergiant Envelopes

Using Athena++, we perform 3D Radiation-Hydrodynamic calculations of the radiative breakout of the shock wave in the outer envelope of a red supergiant (RSG) which has suffered core collapse and will become a Type IIP supernova. The intrinsically 3D structure of the fully convective RSG envelope yields key differences in the brightness and duration of the shock breakout (SBO) from that predicted in a 1D stellar model. First, the lower-density `halo' of material outside of the traditional photosphere in 3D models leads to a shock breakout at lower densities than 1D models. This would prolong the duration of the shock breakout flash at any given location on the surface to $\approx$1-2 hours. However, we find that the even larger impact is the intrinsically 3D effect associated with large-scale fluctuations in density that cause the shock to break out at different radii at different times. This substantially prolongs the SBO duration to $\approx$3-6 hours and implies a diversity of radiative temperatures, as different patches across the stellar surface are at different stages of their radiative breakout and cooling at any given time. These predicted durations are in better agreement with existing observations of SBO. The longer durations lower the predicted luminosities by a factor of 3-10 ($L_\mathrm{bol}\sim10^{44}\mathrm{erg\ s^{-1}}$), and we derive the new scalings of brightness and duration with explosion energies and stellar properties. These intrinsically 3D properties eliminate the possibility of using observed rise times to measure the stellar radius via light-travel time effects.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stochastic Low Frequency Variability in 3-Dimensional Radiation Hydrodynamical Models of Massive Star Envelopes

Increasing main sequence stellar luminosity with stellar mass leads to the eventual dominance of radiation pressure in stellar envelope hydrostatic balance. As the luminosity approaches the Eddington limit, additional instabilities (beyond conventional convection) can occur. These instabilities readily manifest in the outer envelopes of OB stars, where the opacity increase associated with iron yields density and gas pressure inversions in 1D models. Additionally, recent photometric surveys (e.g. TESS) have detected excess broadband low frequency variability in power spectra of OB star lightcurves, called stochastic low frequency variability (SLFV). This motivates our novel 3D Athena++ radiation hydrodynamical (RHD) simulations of two 35$\,$M$_\odot$ star envelopes (the outer $\approx$15$\%$ of the stellar radial extent), one on the zero-age main sequence and the other in the middle of the main sequence. Both models exhibit turbulent motion far above and below the conventional iron opacity peak convection zone (FeCZ), obliterating any ``quiet" part of the near-surface region and leading to velocities at the photosphere of 10-100$\,$km$\,$s$^{-1}$, directly agreeing with spectroscopic data. Surface turbulence also produces SLFV in model lightcurves with amplitudes and power-law slopes that are strikingly similar to those of observed stars. The characteristic frequencies associated with SLFV in our models are comparable to the thermal time in the FeCZ ($\approx$3-7$\,$days$^{-1}$). These simulations, which have no free parameters, are directly validated by observations and, though more models are needed, we remain optimistic that 3D RHD models of main sequence O star envelopes exhibit SLFV originating from the FeCZ.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Cosmic Ray Staircase: the Outcome of the Cosmic Ray Acoustic Instability

Recently, cosmic rays (CRs) have emerged as a leading candidate for driving galactic winds. Small-scale processes can dramatically affect global wind properties. We run two-moment simulations of CR streaming to study how sound waves are driven unstable by phase-shifted CR forces and CR heating. We verify linear theory growth rates. As the sound waves grow non-linear, they steepen into a quasi-periodic series of propagating shocks; the density jumps at shocks create CR bottlenecks. The depth of a propagating bottleneck depends on both the density jump and its velocity; ΔP_c is smaller for rapidly moving bottlenecks. A series of bottlenecks creates a CR staircase structure, which can be understood from a convex hull construction. The system reaches a steady state between growth of new perturbations, and stair mergers. CRs are decoupled at plateaus, but exert intense forces and heating at stair jumps. The absence of CR heating at plateaus leads to cooling, strong gas pressure gradients and further shocks. If bottlenecks are stationary, they can drastically modify global flows; if their propagation times are comparable to dynamical times, their effects on global momentum and energy transfer are modest. The CR acoustic instability is likely relevant in thermal interfaces between cold and hot gas, as well as galactic winds. Similar to increased opacity in radiative flows, the build-up of CR pressure due to bottlenecks can significantly increase mass outflow rates, by up to an order of magnitude. It seeds unusual forms of thermal instability, and the shocks could have distinct observational signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fluid Simulations of Cosmic Ray Modified Shocks

We consider cosmic ray (CR) modified shocks with both streaming and diffusion in the two-fluid description. Previously, numerical codes were unable to incorporate streaming in this demanding regime, and have never been compared against analytic solutions. First, we find a new analytic solution highly discrepant in acceleration efficiency from the standard solution. It arises from bi-directional streaming of CRs away from the subshock, similar to a Zeldovich spike in radiative shocks. Since fewer CRs diffuse back upstream, this results in a much lower acceleration efficiency, typically $\sim 10\%$ as opposed to $\sim 50\%$ found in previous analytic work. At Mach number $\gtrsim 10$, the new solution bifurcates into 3 branches, with efficient, intermediate and inefficient CR acceleration. Our two-moment code (Jiang & Oh 2018) accurately recovers these solutions across the entire parameter space probed, with no ad hoc closure relations. For generic initial conditions, the inefficient branch is the most robust and preferred solution. The intermediate branch is unstable, while the efficient branch appears only when the inefficient branch is not allowed (for CR dominated or high plasma $β$ shocks). CR modified shocks have very long equilibration times ($\sim 1000$ diffusion time) required to develop the precursor, which must be resolved by $\gtrsim 10$ cells for convergence. Non-equilibrium effects, poor resolution and obliquity of the magnetic field all reduce CR acceleration efficiency. Shocks in galaxy scale simulations will generally contribute little to CR acceleration without a subgrid prescription.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori's Accretion Disk and Observational Signatures of Magnetic Fields

FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori's accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori's inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z$\sim$R, similar to the "surface accretion" mechanism in previous locally-isothermal MHD simulations. (2) A moderate disk wind is launched in the vertical field simulations with a terminal speed of $\sim$300-500 km/s and a mass loss rate of 1-10\% the disk accretion rate, which is consistent with observations. Disk wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disk photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disk wind, the disk photosphere temperature is lower than that predicted by the thin disk theory, and the previously inferred disk accretion rate may be lower than the real accretion rate by a factor of $\sim$2-3.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Opacity Driven Convection and Variability in Accretion Disks around Supermassive Black Holes

We study the structure of accretion disks around supermassive black holes in the radial range $30\sim 100$ gravitational radii, using a three dimensional radiation magneto-hydrodynamic simulation. For typical conditions in this region of Active Galactic Nuclei (AGN), the Rosseland mean opacity is expected to be larger than the electron scattering value. We show that the iron opacity bump causes the disk to be convective unstable. Turbulence generated by convection puffs up the disk due to additional turbulent pressure support and enhances the local angular momentum transport. This also results in strong fluctuations in surface density and heating of the disk. The opacity drops with increasing temperature and convection is suppressed. The disk cools down and the whole process repeats again. This causes strong oscillations of the disk scale height and luminosity variations by more than a factor of $\approx 3-6$ over a few years' timescale. Since the iron opacity bump will move to different locations of the disk for black holes with different masses and accretion rates, we suggest that this is a physical mechanism that can explain the variability of AGN with a wide range of amplitudes over a time scale of years to decades.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time Dependent Radiation Hydrodynamics on a Moving Mesh

We describe the structure and implementation of a radiation hydrodynamic solver for MANGA, the moving-mesh hydrodynamics module of the large-scale parallel code, Charm N-body GrAvity solver (ChaNGa). We solve the equations of time dependent radiative transfer using a reduced speed of light approximation following the algorithm of Jiang et al (2014). By writing the radiative transfer equations as a generalized conservation equation, we solve the transport part of these equations on an unstructured Voronoi mesh. We then solve the source part of the radiative transfer equations following Jiang et al (2014) using an implicit solver, and couple this to the hydrodynamic equations. The use of an implicit solver ensure reliable convergence and preserves the conservation properties of these equations even in situations where the source terms are stiff due to the small coupling timescales between radiation and matter. We present the results of a limited number of test cases (energy conservation, momentum conservation, dynamic diffusion, linear waves, crossing beams, and multiple shadows) to show convergence with analytic results and numerical stability. We also show that it produces qualitatively the correct results in the presence of multiple sources in the optically thin case.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Global Radiation Magneto-hydrodynamic Simulations of Sub-Eddington Accretion Disks around Supermassive Black Holes

We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5\times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner \approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Iron Opacity Bump Changes the Stability and Structure of Accretion Disks in Active Galactic Nuclei

Accretion disks around supermassive black holes have regions where the Rosseland mean opacity can be much larger than the electron scattering opacity primarily due to the large number of bound-bound transitions in iron. We study the effects of this iron opacity "bump" on the thermal stability and vertical structure of radiation pressure dominated accretion disks, utilizing three dimensional radiation magneto-hydrodynamic simulations in the local shearing box approximation. The simulations self-consistently calculate the heating due to MHD turbulence caused by magneto-rotational instability and radiative cooling by using the radiative transfer module based on a variable Eddington tensor in Athena. For a $5\times 10^8$ solar mass black hole with $\sim 3\%$ of the Eddington luminosity, a model including the iron opacity bump maintains its structure for more than $10$ thermal times without showing significant signs of thermal runaway. In contrast, if only electron scattering and free-free opacity are included as in the standard thin disk model, the disk collapses on the thermal time scale. The difference is caused by a combination of 1) an anti-correlation between the total optical depth and the temperature, and 2) enhanced vertical advective energy transport. These results suggest that the iron opacity bump may have a strong impact on the stability and structure of AGN accretion disks, and may contribute to a dependence of AGN properties on metallicity. Since this opacity is relevant primarily in UV emitting regions of the flow, it may help to explain discrepancies between observation and theory that are unique to AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Prompt Radiation and Mass Outflows from the Stream-Stream Collisions of Tidal Disruption Events

Stream-stream collisions play an important role for the circularization of highly eccentric streams resulting from tidal disruption events (TDEs). We perform three dimensional radiation hydrodynamic simulations to show that stream collisions can contribute significant optical and ultraviolet light to the flares produced by TDEs, and can sometimes explain the majority of the observed emission. Our simulations focus on the region near the radiation pressure dominated shock produced by a collision and track how the kinetic energy of the stream is dissipated by the associated shock. When the mass flow rate of the stream $\dot{M}$ is a significant fraction of the Eddington accretion rate, $\gtrsim2\%$ of the initial kinetic energy is converted to radiation directly as a result of the collision. In this regime, the collision redistributes the specific kinetic energy into the downstream gas and more than $16\%$ of the mass can become unbound. The fraction of unbound gas decreases rapidly as $\dot{M}$ drops significantly below the Eddington limit, with no unbound gas being produced when $\dot{M}$ drops to $1\%$ of Eddington; we find however that the radiative efficiency increases slightly to $\lesssim 8\%$ in these low $\dot{M}$ cases. The effective radiation temperature and size of the photosphere is determined by the stream velocity and $\dot{M}$, which we find to be a few times $10^4$~K and $10^{14}$~cm in our calculations, comparable to the inferred values of some TDE candidates. The photosphere size is directly proportional to $\dot{M}$, which can explain the rapidly changing photosphere sizes seen in TDE candidates such as PS1-10jh.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Towards Self-Consistent Modelling of the Sgr A* Accretion Flow: Linking Theory and Observation

The interplay between supermassive black holes (SMBHs) and their environments is believed to command an essential role in galaxy evolution. The majority of these SMBHs are in the radiative inefficient accretion phase where this interplay remains elusive, but suggestively important, due to few observational constraints. To remedy this, we directly fit 2-D hydrodynamic simulations to Chandra observations of Sgr A* with Markov Chain Monte Carlo sampling, self-consistently modelling the 2-D inflow-outflow solution for the first time. We find the temperature and density at flow onset are consistent with the origin of the gas in the stellar winds of massive stars in the vicinity of Sgr A*. We place the first observational constraints on the angular momentum of the gas and estimate the centrifugal radius, r$_c$ $\approx$ 0.056 r$_b$ $\approx8\times10^{-3}$ pc, where r$_b$ is the Bondi radius. Less than 1\% of the inflowing gas accretes onto the SMBH, the remainder being ejected in a polar outflow. For the first time, we decouple the quiescent point-like emission from the spatially extended flow. We find this point-like emission, accounting for $\sim4$\% of the quiescent flux, is spectrally too steep to be explained by unresolved flares, nor Bremsstrahlung, but is likely a combination of a relatively steep synchrotron powerlaw and the high energy tail of inverse-Compton emission. With this self-consistent model of the accretion flow structure, we make predictions for the flow dynamics and discuss how future X-ray spectroscopic observations can further our understanding of the Sgr A* accretion flow.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Local Radiation Hydrodynamic Simulations of Massive Star Envelopes at the Iron Opacity Peak

We perform three dimensional radiation hydrodynamic simulations of the structure and dynamics of radiation dominated envelopes of massive stars at the location of the iron opacity peak. One dimensional hydrostatic calculations predict an unstable density inversion at this location, whereas our simulations reveal a complex interplay of convective and radiative transport whose behavior depends on the ratio of the photon diffusion time to the dynamical time. The latter is set by the ratio of the optical depth per pressure scale height, $τ_0$, to $τ_c=c/c_g$, where $c_g \approx$ 50 km/s is the isothermal sound speed in the gas alone. When $τ_0 \gg τ_c$, convection reduces the radiation acceleration and removes the density inversion. The turbulent energy transport in the simulations agrees with mixing length theory and provides its first numerical calibration in the radiation dominated regime. When $τ_0 \ll τ_c$, convection becomes inefficient and the turbulent energy transport is negligible. The turbulent velocities exceed $c_g$, driving shocks and large density fluctuations that allow photons to preferentially diffuse out through low-density regions. However, the effective radiation acceleration is still larger than the gravitational acceleration so that the time average density profile contains a modest density inversion. In addition, the simulated envelope undergoes large-scale oscillations with periods of a few hours. The turbulent velocity field may affect the broadening of spectral lines and therefore stellar rotation measurements in massive stars, while the time variable outer atmosphere could lead to variations in their mass loss and stellar radius.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Global Three Dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks

We study super-Eddington accretion flows onto black holes using a global three dimensional radiation magneto-hydrodynamical simulation. We solve the time dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ~220L_edd/c^2 and forms a radiation driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ~20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ~10L_edd. This yields a radiative efficiency ~4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed published results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation

(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the Athena code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers. The method is ideally suited for problems in which characteristic velocities are non-relativistic, but still within a few percent or more of the speed of light. The method is an intermediate step towards algorithms for fully relativistic flows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radiation Magneto-hydrodynamic Simulations of the Formation of Hot Accretion Disk Coronae

A new mechanism to form a magnetic pressure supported, high temperature corona above the photosphere of an accretion disk is explored using three dimensional radiation magneto-hydrodynamic (MHD) simulations. The thermal properties of the disk are calculated self-consistently by balancing radiative cooling through the surfaces of the disk with heating due to dissipation of turbulence driven by magneto-rotational instability (MRI). As has been noted in previous work, we find the dissipation rate per unit mass increases dramatically with height above the mid-plane, in stark contrast to the alpha-disk model which assumes this quantity is a constant. Thus, we find that in simulations with a low surface density (and therefore a shallow photosphere), the fraction of energy dissipated above the photosphere is significant (about 3.4% in our lowest surface density model), and this fraction increases as surface density decreases. When a significant fraction of the accretion energy is dissipated in the optically thin photosphere, the gas temperature increases substantially and a high temperature, magnetic pressure supported corona is formed. The volume-averaged temperature in the disk corona is more than 10 times larger than at the disk mid-plane. Moreover, gas temperature in the corona is strongly anti-correlated with gas density, which implies the corona formed by MRI turbulence is patchy. This mechanism to form an accretion disk corona may help explain the observed relation between the spectral index and luminosity from AGNs, and the soft X-ray excess from some AGNs. It may also be relevant to spectral state changes in X-ray binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Thermal Stability of Radiation Dominated Accretion Disks

We study the long-term thermal stability of radiation dominated disks in which the vertical structure is determined self-consistently by the balance of heating due to dissipation of MHD turbulence driven by the magneto-rotational instability (MRI), and cooling due to radiation emitted at the photosphere. The calculations adopt the local shearing box approximation, and utilize the recently developed radiation transfer module in the Athena MHD code based on a variable Eddington tensor rather than an assumed local closure. After saturation of the MRI, in many cases the disk maintains a steady vertical structure for many thermal times. However, in every case in which the box size in the horizontal directions is at least one pressure scale height, fluctuations associated with MRI turbulence and dynamo action in the disk eventually trigger a thermal runaway which causes the disk to either expand or contract until the calculation must be terminated. During runaway, the dependence of the heating and cooling rates on total pressure satisfy the simplest criterion for classical thermal instability. We identify several physical reasons why the thermal runaway observed in our simulations differ from the standard alpha disk model, for example the advection of radiation contributes a non-negligible fraction to the vertical energy flux at the largest radiation pressure, most of the dissipation does not happen in the disk mid-plane, and the change of dissipation scale height with mid-plane pressure is slower than the change of density scale height. We discuss how and why our results differ from those published previously. Such thermal runaway behavior might have important implications for interpreting temporal variability in observed systems, but fully global simulations are required to study the saturated state before detailed predictions can be made.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Saturation of the MRI in Strongly Radiation Dominated Accretion Disks

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches 1% of the rest mass energy density of the gas), we find Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of 8 compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find radiative shear viscosity is negligible compared to the Maxwell and Reynolds stress in the flow. This may have important implications for the structure of radiation dominated accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor

We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusion-like approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiation-modified linear acoustic and magnetosonic waves, the structure of radiation modified shocks, and two-dimensional tests of photon bubble instability and the ablation of dense clouds by an intense radiation field. These tests cover a very wide range of regimes, including both optically thick and thin flows, and ratios of the radiation to gas pressure of at least 10^{-4} to 10^{4}. Across most of the parameter space, we find the method is accurate. However, the tests also reveal there are regimes where the method needs improvement, for example when both the radiation pressure and absorption opacity are very large. We suggest modifications to the algorithm that will improve accuracy in this case. We discuss the advantages of this method over those based on flux-limited diffusion. In particular, we find the method is not only substantially more accurate, but often no more expensive than the diffusion approximation for our intended applications.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Way to Conserve Total Energy for Eulerian Hydrodynamic Simulations with Self-Gravity

We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the the self-gravitational force is given, so the change in total energy can be written in conservative form. Numerical experiments with the code Athena show that the total energy is indeed conserved with our new algorithm and the new algorithm is second order accurate. We have performed a set of tests that show the numerical errors in the traditional, non-conservative algorithm can affect the dynamics of the system. The new algorithm only requires one extra solution of the Poisson equation, as compared to the traditional algorithm which includes self-gravity as a source term. If the Poisson solver takes a negligible fraction of the total simulation time, such as when FFTs are used, the new algorithm is almost as efficient as the original method. This new algorithm is useful in Eulerian hydrodynamic simulations with self-gravity, especially when results are sensitive to small energy errors, as for radiation pressure dominated flow.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Radiation Transfer Solver for Athena using Short Characteristics

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-LTE effects. The module is based on well-known and well-tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator-splitting and we describe this approach in detail. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself, and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate, and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad-hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g. LTE grey) problems can be comparable to or less than a single timestep of Athena's MHD integrators, and only few times more expensive than that for more general problems. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-linear Evolution of Rayleigh-Taylor Instability in a Radiation Supported Atmosphere

The non-linear regime of Rayleigh-Taylor instability (RTI) in a radiation supported atmosphere, consisting of two uniform fluids with different densities, is studied numerically. We perform simulations using our recently developed numerical algorithm for multi-dimensional radiation hydrodynamics based on a variable Eddington tensor as implemented in Athena, focusing on the regime where scattering opacity greatly exceeds absorption opacity. We find that the radiation field can reduce the growth and mixing rate of RTI, but this reduction is only significant when radiation pressure significantly exceeds gas pressure. Small scale structures are also suppressed in this case. In the non-linear regime, dense fingers sink faster than rarefied bubbles can rise, leading to asymmetric structures about the interface. By comparing the calculations that use a variable Eddington tensor (VET) versus the Eddington approximation, we demonstrate that anisotropy in the radiation field can affect the non-linear development of RTI significantly. We also examine the disruption of a shell of cold gas being accelerated by strong radiation pressure, motivated by models of radiation driven outflows in ultraluminous infrared galaxies. We find that when the growth rate of RTI is smaller than acceleration time scale, the amount of gas that would be pushed away by the radiation field is reduced due to RTI.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Black Hole Mass and Bulge Luminosity for Low-mass Black Holes

We study the scaling between bulge magnitude and central black hole (BH) mass in galaxies with virial BH masses < 10^6 solar mass. Based on careful image decomposition of a snapshot Hubble Space Telescope I-band survey, we found that these BHs are found predominantly in galaxies with pseudobulges. Here we show that the \mbulge\ relation for the pseudobulges at low mass is significantly different from classical bulges with BH masses >10^7 solar mass. Specfically, bulges span a much wider range of bulge luminosity, and on average the luminosity is larger, at fixed black hole mass. The trend holds both for the active galaxies from Bentz et al. and the inactive sample of Gultekin et al. and cannot be explained by differences in stellar populations, as it persists when we use dynamical bulge masses. Put another way, the ratio between bulge and BH mass is much larger than $\sim 1000$ for our sample. This is consistent with recent suggestions that black hole mass does not scale with the pseudobulge luminosity. The low-mass scaling relations appear to flatten, consistent with predictions from Volonteri & Natarajan for massive seed BHs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Host Galaxies of Low-mass Black Holes

Using HST observations of 147 host galaxies of low-mass black holes (BHs), we systematically study the structures and scaling relations of these active galaxies. Our sample is selected to have central BHs with virial masses ~10^5-10^6 solar mass. The host galaxies have total I-band magnitudes of -23.2<M_I<-18.8 mag and bulge magnitudes of -22.9<M_I<-16.1 mag. Detailed bulge-disk-bar decompositions with GALFIT show that 93% of the galaxies have extended disks, 39% have bars and 5% have no bulges at all at the limits of our observations. Based on the Sersic index and bulge-to-total ratio, we conclude that the majority of the galaxies with disks are likely to contain pseudobulges and very few of these low-mass BHs live in classical bulges. The fundamental plane of our sample is offset from classical bulges and ellipticals in a way that is consistent with the scaling relations of pseudobulges. The sample has smaller velocity dispersion at fixed luminosity in the Faber-Jackson plane, compared with classical bulges and elliptical galaxies. The galaxies without disks are structurally more similar to spheroidals than to classical bulges according to their positions in the fundamental plane, especially the Faber-Jackson projection. Overall, we suggest that BHs with mass < 10^6 solar mass live in galaxies that have evolved secularly over the majority of their history. A classical bulge is not a prerequisite to host a black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Evolution of Wide Binary Stars

We study the orbital evolution of wide binary stars in the solar neighborhood due to gravitational perturbations from passing stars. We include the effects of the Galactic tidal field and continue to follow the stars after they become unbound. For a wide variety of initial semi-major axes and formation times, we find that the number density (stars per unit logarithmic interval in projected separation) exhibits a minimum at a few times the Jacobi radius r_J, which equals 1.7 pc for a binary of solar-mass stars. The density peak interior to this minimum arises from the primordial distribution of bound binaries, and the exterior density, which peaks at \sim 100--300 pc separation, arises from formerly bound binaries that are slowly drifting apart. The exterior peak gives rise to a significant long-range correlation in the positions and velocities of disk stars that should be detectable in large astrometric surveys such as GAIA that can measure accurate three-dimensional distances and velocities.

0 citations0 reviewsNo code linkedOpen access