BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Yan-Fei Jiang

Yan-Fei Jiang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRastro-ph.GAastro-ph.IMastro-ph.COphysics.comp-phastro-ph.EPphysics.plasm-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 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
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
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
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

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