Search

We describe the implementation of multigrid solvers in the Athena++ adaptive mesh refinement (AMR) framework and their application to the solution of the Poisson equation for self-gravity. The new solvers are built on top of the AMR hierarchy and TaskList framework of Athena++ for efficient parallelization. We adopt a conservative formulation for the Laplacian operator that avoids artificial accelerations at level boundaries. Periodic, fixed, and zero-gradient boundary conditions are implemented, as well as open boundary conditions based on a multipole expansion. Hybrid parallelization using both MPI and OpenMP is adopted, and we present results of tests demonstrating the accuracy and scaling of the methods. On a uniform grid we show multigrid significantly outperforms methods based on FFTs, and requires only a small fraction of the compute time required by the (highly optimized) magnetohydrodynamic solver in Athena++. As a demonstration of the capabilities of the methods, we present the results of a test calculation of magnetized protostellar collapse on an adaptive mesh., Comment: 28 pages, 13 figures, submitted to AAS Journals

We present global 3D radiation magnetohydrodynamic simulations of accretion onto a 6.62 solar-mass black hole, with quasi-steady-state accretion rates reaching 0.016–0.9 times the critical accretion rate, which is defined as the accretion rate for powering the Eddington luminosity, assuming a 10% radiative efficiency, in three different runs. The simulations show no sign of thermal instability over hundreds of thermal timescales at 10 r g. The energy dissipation occurs close to the mid-plane in the near-critical runs and near the disk surface in the low–accretion rate run. The total radiative luminosity inside ∼20 r g is about 1%–30% of the Eddington limit, with radiative efficiencies of about 6% and 3%, respectively, in the sub- and near-critical accretion regimes. In both cases, self-consistent turbulence generated by the magnetorotational instability leads to angular momentum transfer, and the disk is supported by magnetic pressure. Outflows from the central low-density funnel, with a terminal velocity of ∼0.1c, are seen only in the near-critical runs. We conclude that these magnetic pressure–dominated disks are thermally stable and thicker than the α disk, and that the effective temperature profiles are much flatter than those in the α disks. The magnetic pressures of these disks are comparable within an order of magnitude to the previous analytical magnetic pressure–dominated disk model.

We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new GPU-accelerated version of the Athena++ AMR code, and span more than 6 orders of magnitude in radius, reaching scales similar to the black hole horizon. The key physical ingredients are radiative cooling and a phenomenological heating model. We find that the accretion flow takes the form of multiphase gas at radii less than about a kpc. The cold gas accretion includes two dynamically distinct stages: the typical disk stage in which the cold gas resides in a rotationally supported disk and relatively rare chaotic stages ($\lesssim 10\%$ of the time) in which the cold gas inflows via chaotic streams. Though cold gas accretion dominates the time-averaged accretion rate at intermediate radii, accretion at the smallest radii is dominated by hot virialized gas at most times. The accretion rate scales with radius as $\dot{M}\propto r^{1/2}$ when hot gas dominates and we obtain $\dot{M}\simeq10^\mathrm{-4}-10^\mathrm{-3}\,M_\odot\,\mathrm{yr^{-1}}$ near the event horizon, similar to what is inferred from EHT observations. The orientation of the cold gas disk can differ significantly on different spatial scales. We propose a subgrid model for accretion in lower-resolution simulations in which the hot gas accretion rate is suppressed relative to the Bondi rate by $\sim (r_\mathrm{g}/r_{\rm Bondi})^{1/2}$. Our results can also provide more realistic initial conditions for simulations of black hole accretion at the event horizon scale., Comment: 28 pages, 25 figures, 2 tables, accepted for publication in ApJ

Magnetic fields on the order of 100 $\mu$G observed in young supernova remnants cannot be amplified by shock compression alone. To investigate the amplification caused by turbulent dynamo, we perform three-dimensional MHD simulations of the interaction between shock wave and inhomogeneous density distribution with a shallow spectrum in the preshock medium. The postshock turbulence is mainly driven by the strongest preshock density contrast and follows the Kolmogorov scaling. The resulting turbulence amplifies the postshock magnetic field. The time evolution of the magnetic fields agrees with the prediction of the nonlinear turbulent dynamo theory in Xu & Lazarian (2016). When the initial weak magnetic field is perpendicular to the shock normal, the maximum amplification of the field's strength achieves a factor of $\approx200$, which is twice larger than that for a parallel shock. We find that the perpendicular shock exhibits a smaller turbulent Alfv\'en Mach number in the vicinity of the shock front than the parallel shock. However, the strongest magnetic field has a low volume filling factor and is limited by the turbulent energy due to the reconnection diffusion taking place in a turbulent and magnetized fluid. The magnetic field strength averaged along the $z$-axis is reduced by a factor $\gtrsim10$. We decompose the turbulent velocity and magnetic field into solenoidal and compressive modes. The solenoidal mode is dominant and evolves to follow the Kolmogorov scaling, even though the preshock density distribution has a shallow spectrum. When the preshock density distribution has a Kolmogorov spectrum, the turbulent velocity's compressive component increases., Comment: 19 pages, 15 figures, accepted for publication in ApJ

We present cosmological hydrodynamic simulations of a quasar-mass halo ($M_{\rm halo} \approx 10^{12.5}\,{\rm M}_{\odot}$ at z=2) that for the first time resolve gas transport down to the inner 0.1 pc surrounding the central massive black hole. We model a multi-phase interstellar medium including stellar feedback by supernovae, stellar winds, and radiation, and a hyper-Lagrangian refinement technique increasing the resolution dynamically approaching the black hole. We do not include black hole feedback. We show that the sub-pc inflow rate (1) can reach ~6 M$_{\odot}$yr$^{-1}$ roughly in steady state during the epoch of peak nuclear gas density (z~2), sufficient to power a luminous quasar, (2) is highly time variable in the pre-quasar phase, spanning 0.001-10 M$_{\odot}$yr$^{-1}$ on Myr timescales, and (3) is limited to short (~2 Myr) active phases (0.01-0.1 M$_{\odot}$yr$^{-1}$) followed by longer periods of inactivity at lower nuclear gas density and late times (z~1), owing to the formation of a hot central cavity. Inflowing gas is primarily cool, rotational support dominates over turbulence and thermal pressure, and star formation can consume as much gas as provided by inflows across 1 pc - 10 kpc. Gravitational torques from multi-scale stellar non-axisymmetries dominate angular momentum transport over gas self-torquing and pressure gradients, with accretion weakly dependent on black hole mass. Sub-pc inflow rates correlate with nuclear (but decouple from global) star formation and can exceed the Eddington rate by x10. The black hole can move ~10 pc from the galaxy center on ~0.1 Myr. Accreting gas forms pc-scale, rotationally supported, obscuring structures often misaligned with the galaxy-scale disk. These simulations open a new avenue to investigate black hole-galaxy co-evolution., Comment: 34 pages, 23 figures, ApJ accepted

We use 2D (axisymmetric) and 3D hydrodynamic simulations to study Bondi-Hoyle-Lyttleton (BHL) accretion with and without transverse upstream gradients. We mainly focus on the regime of high (upstream) Mach number, weak upstream gradients and small accretor size, which is relevant to neutron star (NS) accretion in wind-fed Supergiant X-ray binaries (SgXBs). We present a systematic exploration of the flow in this regime. When there are no upstream gradients, the flow is always stable regardless of accretor size or Mach number. For finite upstream gradients, there are three main types of behavior: stable flow (small upstream gradient), turbulent unstable flow without a disk (intermediate upstream gradient), and turbulent flow with a disk-like structure (relatively large upstream gradient). When the accretion flow is turbulent, the accretion rate decreases non-convergently as the accretor size decreases. The flow is more prone to instability and the disk is less likely to form than previously expected; the parameters of most observed SgXBs place them in the regime of a turbulent, disk-less accretion flow. Among the SgXBs with relatively well-determined parameters, we find OAO 1657-415 to be the only one that is likely to host a persistent disk (or disk-like structure); this finding is consistent with observations., Comment: 25 pages, 28 figures, 2 tables. Accepted by MNRAS

In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^{+}$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We carry out simulations of MHD turbulence in molecular clouds including heating and cooling, and post-process them to study H$_2$ line emission and hot-gas chemistry, particularly the formation of CH$^+$. We explore multiple magnetic field strengths and equations of state. We use a new H$_2$ cooling function for $n_{\rm H} \leq 10^5\,{\rm cm}^{-3}$, $T\leq 5000\,{\rm K}$, and variable H$_2$ fraction. We make two important simplifying assumptions: (i) the ${\rm H}_2/{\rm H}$ fraction is fixed everywhere, and (ii) we exclude from our analysis regions where the ion-neutral drift velocity is calculated to be greater than 5 km/s. Our models produce H$_2$ emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic r.m.s. magnetic field strengths ($\approx 10$ $\mu$G) and velocity dispersions, we reproduce observed CH$^+$ abundances. These findings contrast with those of Valdivia et al. (2017). Comparison of predicted dust polarization with observations by {\it Planck} suggests that the mean field $\gtrsim 5 \mu$G, so that the turbulence is sub-Alfv\'enic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH$^+$ abundance., Comment: 20 pages, 18 figures, accepted for publication to MNRAS

The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both non-relativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available., 50 pages, 41 figures, accepted for publication in American Astronomical Society journals

We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $\sim 10^{-8} M_\odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation., Comment: Accepted by ApJL. 14 pages, 8 figures. Animations available at https://smressle.bitbucket.io/animations.html

We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion onto Sagittarius A* via the magnetized winds of the orbiting Wolf-Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach $\approx$ 300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g., plasma $\beta$ $\sim$ $10^6$), flux freezing/compression in the inflowing gas amplifies the field to $\beta$ $\sim$ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g., density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of i) strong magnetic fields that are amplified by flux freezing/compression, and ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulations is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*., Comment: Accepted by MNRAS; 22 pages; 24 figures, comments welcome

We present results of the first global magnetohydrodynamic (MHD) simulations of accretion disks fed by Roche lobe overflow, including vertical stratification, in order to investigate the roles of spiral shocks, magnetorotational instability (MRI), and the accretion stream on disk structure and evolution. Our models include a simple treatment of gas thermodynamics, with orbital Mach numbers at the inner edge of the disk $M_{\rm in}$ of 5 and 10. We find mass accretion rates to vary considerably on all time scales, with only the Mach 5 model reaching a clear quasi-stationary state. For Mach 10, the model undergoes an outside-in magnetically-driven accretion event occurring on a time scale of $\sim10$ orbital periods of the binary. Both models exhibit spiral shocks inclined with respect to the binary plane, with their position and inclination changing rapidly. However, the time-averaged location of these shocks in the equatorial plane is well-fit by simple linear models. MRI turbulence in the disk generates toroidal magnetic field patterns (butterfly diagrams) that are in some cases irregular, perhaps due to interaction with spiral structure. While many of our results are in good agreement with local studies, we find some features (most notably those related to spiral shocks) can only be captured in global models such as studied here. Thus, while global studies remain computationally expensive -- even as idealized models -- they are essential (along with more sophisticated treatment of radiation transport and disk thermodynamics) for furthering our understanding of accretion in binary systems., Comment: 30 pages, 13 figures; Accepted for publication in ApJ

Highly magnetized accretion disks are present in high-mass X-ray binaries (HMXBs). A potential mechanism to explain the transition between the High/Soft and Low/Hard states observed in HMXBs can be attributed to fast magnetic reconnection induced in the turbulent corona. In this work, we present results of global general relativistic MHD (GRMHD) simulations of accretion disks around black holes that show that fast reconnection events can naturally arise in the coronal region of these systems in presence of turbulence triggered by MHD instabilities, indicating that such events can be a potential mechanism to explain the transient non-thermal emission in HMXBs. To find the zones of fast reconnection, we have employed an algorithm to identify the presence of current sheets in the turbulent regions and computed statistically the magnetic reconnection rates in these locations obtaining average reconnection rates consistent with the predictions of the theory of turbulence-induced fast reconnection.

The largest ever simulation of astrophysical turbulence substantially improves our understanding of how energy injection on large interstellar scales governs how stars form on small scales.

We investigate the properties of the hydrodynamic flow around eccentric protoplanets and compare them with the often assumed case of a circular orbit. To this end, we perform a set of 3D hydrodynamic simulations of protoplanets with small eccentricities ($e\leq 0.1$). We adopt an isothermal equation of state and concentrate resolution on the protoplanet to investigate flows down to the scale of the protoplanet's circumplanetary disk (CPD). We find enhanced prograde rotation exterior to the CPD for low planet masses undergoing subsonic eccentric motion. If the eccentricity is made large enough to develop a bow shock, this trend reverses and rotation becomes increasingly retrograde. The instantaneous eccentric flow field is dramatically altered compared to circular orbits. Whereas the latter exhibit a generic pattern of polar inflow and midplane outflow, the flow geometry depends on orbital phase in the eccentric case. For even the modest eccentricities tested here, the dominant source of inflow can come from the midplane instead of the poles. We find that the amount of inflow and outflow increases for higher $e$ and lower protoplanet masses, thereby recycling more gas through the planet's Bondi radius. These increased fluxes may increase the pebble accretion rate for eccentric planets up to several times that of the circular orbit rate. In response to eccentric motion, the structure and rotation of the planet's bound CPD remains unchanged. Because the CPD regulates the eventual accretion of gas onto the planet, we predict little change to the gas accretion rates between eccentric and circular planets.

Gravitationally focused wind accretion in binary systems consisting of an evolved star with a gaseous envelope and a compact accreting companion is a possible mechanism to explain mass transfer in symbiotic binaries. We study the mass accretion around the secondary caused by the strong wind from the primary late-type component using global three-dimensional hydrodynamic numerical simulations during quiescence and outburst stages. In particular, the dependence of the mass accretion rate on the mass-loss rate, wind parameters and phases of wind outburst development is considered. For a typical wind from an asymptotic giant branch star with a mass-loss rate of 1e-6 Msun/year and wind speeds of 20-50 km/s, the mass transfer through a focused wind results in efficient infall on to the secondary. Accretion rates onto the secondary of 5-20 per cent of the mass-loss from the primary are obtained during quiescence and outburst periods where the wind velocity and mass-loss rates are varied, about 20-50 per cent larger than in the standard Bondi-Hoyle-Lyttleton approximation. This mechanism could be an important method for explaining observed accretion luminosities and periodic modulations in the accretion rates for a broad range of interacting binary systems., Comment: 11 pages, 8 figures, to be published in MNRAS

Recent high resolution observations unveil ring structures in circumstellar disks. The origin of these rings has been widely investigated under various theoretical scenarios. In this work we perform global 3D non-ideal MHD simulations including effects from both Ohmic resistivity and ambipolar diffusion (AD) to model the HL Tau disk. The non-ideal MHD diffusion profiles are calculated based on the global dust evolution calculation including sintering effects. Disk ionization structure changes dramatically across the snow line due to the change of dust size distribution close to snow line of major volatiles. We find that accretion is mainly driven by disk wind. Gaps and rings can be quickly produced from different accretion rates across snow line. Furthermore, ambipolar diffusion (AD) leads to highly preferential accretion at midplane, followed by magnetic reconnection. This results a local zone of decretion that drains of mass in the field reconnection area, which leaves a gap and an adjacent ring just outside it. Overall, under the favorable condition, both snow lines and non-ideal MHD effects can lead to gaseous gaps and rings in protoplanetary disks., 13 pages, 10 figures, published on ApJ

Gaseous disks have been proposed as a mechanism for facilitating mergers of binary black holes. We explore circumbinary disk systems to determine the evolution of the central binary. To do so, we perform 3D, hydrodynamic, locally isothermal simulations of circumbinary disks on a Cartesian grid. We focus on binaries of equal mass ratios on fixed circular orbits. To investigate the orbital evolution of the binary, we examine the various torques exerted on the system. For the case where the disk plane and binary orbital plane are aligned, we find that the total torque is positive so that the semi-major axis of the binary increases. For the misaligned case, we run simulations with the binary orbital plane and disk midplane misaligned by 45\degree and find the same results - the binary grows. The timescale for the circumbinary disk to realign to the plane of the binary is consistent with the global viscous timescale of the disk., Accepted to ApJ. 13 pages, 19 figures

We investigate numerical convergence in simulations of magnetically arrested disks around spinning black holes. Using the general-relativistic magnetohydrodynamics code Athena++, we study the same system at four resolutions (up to an effective 512 by 256 by 512 cells) and with two different spatial reconstruction algorithms. The accretion rate and general large-scale structure of the flow agree across the simulations. This includes the amount of magnetic flux accumulated in the saturated state and the ensuing suppression of the magnetorotational instability from the strong field. The energy of the jet and the efficiency with which spin energy is extracted via the Blandford-Znajek process also show convergence. However the spatial structure of the jet shows variation across the set of grids employed, as do the Lorentz factors. Small-scale features of the turbulence, as measured by correlation lengths, are not fully converged. Despite convergence of a number of aspects of the flow, modeling of synchrotron emission shows that variability is not converged and decreases with increasing resolution even at our highest resolutions., Comment: 19 pages, 18 figures, to be published in ApJ

The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfv\'en waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($\lambda \ll \rm pc$), low cosmic ray number density ($n_{\rm CR}/n_{\rm ISM} \sim 10^{-9}$), and weak CR anisotropy ($\sim v_A/c$) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic-particle-in-cell (MHD-PIC) method, which has unique advantages to alleviate these issues, we conduct one-dimensional simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the $\delta f$ method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum, equally shared between left and right handed Alfv\'en modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires pitch angles of most CRs to efficiently cross $90^\circ$, and can be captured in simulations with relatively high wave amplitude and/or high spatial resolution. We attribute this crossing to non-linear wave-particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate our methodology will open up opportunities for future investigations that incorporate additional physics., Comment: 34 pages, 28 figures, submitted to ApJ

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., Comment: 22 pages, 20 figures, accepted by MNRAS

The dynamical evolution of protoplanetary disks is of key interest for building a comprehensive theory of planet formation and to explain the observational properties of these objects. Using the magnetohydrodynamics code Athena++, with an isothermal shearing box setup, we study the boundary between the active and dead zone, where the accretion rate changes and mass can accumulate. We quantify how the turbulence level is affected by the presence of a non uniform ohmic resistivity in the radial-x direction that leads to a region of inhibited turbulence (or dead zone). Comparing the turbulent activity to that of ideal simulations, the turbulence inhibited area shows density fluctuations and magnetic activity at its boundaries, driven by energy injection from the active (ideal) zone boundaries. We find magnetic dissipation to be significantly stronger in the ideal regions, and the turbulence penetration through the boundary of the dead zone is determined by the value of the resistivity itself, through the ohmic dissipation process, though the thickness of the transition does not play a significant role in changing the dissipation. We investigate the 1D spectra along the shearing direction: magnetic spectra appear flat at large scales both in ideal as well as resistive simulations, though a Kolmogorov scaling over more than one decade persists in the dead zone, suggesting the turbulent cascade is determined by the hydrodynamics of the system: MRI dynamo action is inhibited where sufficiently high resistivity is present., 11 pages 12 figures

The dynamics of solid bodies in protoplanetary disks are subject to the properties of any underlying gas turbulence. Turbulence driven by disk self-gravity shows features distinct from those driven by the magnetorotational instability (MRI). We study the dynamics of solids in gravito-turbulent disks with two-dimensional (in the disk plane), hybrid (particle and gas) simulations. Gravito-turbulent disks can exhibit stronger gravitational stirring than MRI-active disks, resulting in greater radial diffusion and larger eccentricities and relative speeds for large particles (those with dimensionless stopping times $t_{stop} \Omega > 1$, where $\Omega$ is the orbital frequency). The agglomeration of large particles into planetesimals by pairwise collisions is therefore disfavored in gravito-turbulent disks. However, the relative speeds of intermediate-size particles $t_{stop} \Omega \sim 1$ are significantly reduced as such particles are collected by gas drag and gas gravity into coherent filament-like structures with densities high enough to trigger gravitational collapse. First-generation planetesimals may form via gravitational instability of dust in marginally gravitationally unstable gas disks., Comment: MNRAS accepted

Previous studies of the nonlinear regime of the magnetorotational instability in one particular type of shearing box model -- unstratified with no net magnetic flux -- find that without explicit dissipation (viscosity and resistivity) the saturation amplitude decreases with increasing numerical resolution. We show that this result is strongly dependent on the vertical aspect ratio of the computational domain $L_z/L_x$. When $L_z/L_x\lesssim 1$, we recover previous results. However, when the vertical domain is extended $L_z/L_x \gtrsim 2.5$, we find the saturation level of the stress is greatly increased (giving a ratio of stress to pressure $\alpha \gtrsim 0.1$), and moreover the results are independent of numerical resolution. Consistent with previous results, we find that saturation of the MRI in this regime is controlled by a cyclic dynamo which generates patches of strong toroidal field that switches sign on scales of $L_x$ in the vertical direction. We speculate that when $L_z/L_x\lesssim 1$, the dynamo is inhibited by the small size of the vertical domain, leading to the puzzling dependence of saturation amplitude on resolution. We show that previous toy models developed to explain the MRI dynamo are consistent with our results, and that the cyclic pattern of toroidal fields observed in stratified shearing box simulations (leading to the so-called butterfly diagram) may also be related. In tall boxes the saturation amplitude is insensitive to whether or not explicit dissipation is included in the calculations, at least for large magnetic Reynolds and Prandtl number. Finally, we show MRI turbulence in tall domains has a smaller critical $\rm{Pm}_c$, and an extended lifetime compared to $L_z/L_x\lesssim 1$ boxes., Comment: MNRAS in press; paper with higher resolution figures can be found at http://www.astro.princeton.edu/~jmshi/tallbox/ms_tallbox.pdf

The nonlinear evolution of the Kelvin-Helmholtz instability is a popular test for code verification. To date, most Kelvin-Helmholtz problems discussed in the literature are ill-posed: they do not converge to any single solution with increasing resolution. This precludes comparisons among different codes and severely limits the utility of the Kelvin-Helmholtz instability as a test problem. The lack of a reference solution has led various authors to assert the accuracy of their simulations based on ad-hoc proxies, e.g., the existence of small-scale structures. This paper proposes well-posed Kelvin-Helmholtz problems with smooth initial conditions and explicit diffusion. We show that in many cases numerical errors/noise can seed spurious small-scale structure in Kelvin-Helmholtz problems. We demonstrate convergence to a reference solution using both Athena, a Godunov code, and Dedalus, a pseudo-spectral code. Problems with constant initial density throughout the domain are relatively straightforward for both codes. However, problems with an initial density jump (which are the norm in astrophysical systems) exhibit rich behavior and are more computationally challenging. In the latter case, Athena simulations are prone to an instability of the inner rolled-up vortex; this instability is seeded by grid-scale errors introduced by the algorithm, and disappears as resolution increases. Both Athena and Dedalus exhibit late-time chaos. Inviscid simulations are riddled with extremely vigorous secondary instabilities which induce more mixing than simulations with explicit diffusion. Our results highlight the importance of running well-posed test problems with demonstrated convergence to a reference solution. To facilitate future comparisons, we include the resolved, converged solutions to the Kelvin-Helmholtz problems in this paper in machine-readable form., Comment: Reference solution snapshots can be found at https://doi.org/10.5281/zenodo.5803759

Recent observations have revealed that the remnants of stellar-coalescence transients are bipolar. This raises the questions of how these bipolar morphologies arise and what they teach us about the mechanisms of mass ejection during stellar mergers and common-envelope phases. In this paper, we analyze hydrodynamic simulations of the lead-in to binary coalescence, a phase of unstable Roche lobe overflow that takes the binary from the Roche limit separation to the engulfment of the more compact accretor within the envelope of the extended donor. As mass transfer runs away at increasing rates, gas trails away from the binary. Contrary to previous expectations, early mass loss from the system remains bound to the binary and forms a circumbinary torus. Later ejecta, generated as the accretor grazes the surface of the donor, have very different morphologies and are unbound. These two components of mass loss from the binary interact as later, higher-velocity ejecta collide with the circumbinary torus formed by earlier mass loss. Unbound ejecta are redirected toward the poles, and escaping material creates a bipolar outflow. Our findings show that the transition from bound to unbound ejecta from coalescing binaries can explain the bipolar nature of their remnants, with implications for our understanding of the origin of bipolar remnants of stellar-coalescence transients and, perhaps, some preplanetary nebulae., ApJ accepted

We perform relativistic hydrodynamic simulations of internal shocks formed in microquasar jets by continuous variation of the bulk Lorentz factor, in order to investigate the internal shock model. We consider one-, two-, and flicker noise 20-mode variability. We observe emergence of a forward-reverse shock structure for each peak of the Lorentz factor modulation. The high pressure in the shocked layer launches powerful outflows perpendicular to the jet beam into the ambient medium. These outflows dominate the details of the jet's kinetic energy thermalization. They are responsible for mixing between the jet and surrounding medium and generate powerful shocks in the latter. These results do not concur with the popular picture of well-defined internal shells depositing energy as they collide within the confines of the jet, in fact collisions between internal shells themselves are quite rare in our continuous formulation of the problem. For each of our simulations, we calculate the internal energy deposited in the system, the "efficiency" of this deposition (defined as the ratio of internal to total flow energy), and the maximum temperature reached in order to make connections to emission mechanisms. We probe the dependence of these diagnostics on the Lorentz factor variation amplitudes, modulation frequencies, as well as the initial density ratio between the jet and the ambient medium., Comment: 18 pages, 15 figures; accepted for publication in MNRAS

The observed rotation measures (RMs) towards the galactic centre magnetar and towards Sagittarius A* provide a strong constraint on MHD models of the galactic centre accretion flow, probing distances from the black hole separated by many orders of magnitude. We show, using 3D simulations of accretion via magnetized stellar winds of the Wolf-Rayet stars orbiting the black hole, that the large, time-variable RM observed for the pulsar PSR J1745-2900 can be explained by magnetized wind-wind shocks of nearby stars in the clockwise stellar disc. In the same simulation, both the total X-ray luminosity integrated over 2-10$''$, the time variability of the magnetar's dispersion measure, and the RM towards Sagittarius A* are consistent with observations. We argue that (in order for the large RM of the pulsar to not be a priori unlikely) the pulsar should be on an orbit that keeps it near the clockwise disc of stars. We present a 2D RM map of the central 1/2 parsec of the galactic centre that can be used to test our models. Our simulations predict that Sgr A* is typically accreting a significantly ordered magnetic field that ultimately could result in a strongly magnetized flow with flux threading the horizon at $\sim$ 10$\%$ of the magnetically arrested limit., Comment: Accepted into MNRAS Letters. Comments welcome

Luminous red nova transients, presumably from stellar coalescence, exhibit long-term precursor emission over hundreds of binary orbits, leading to impulsive outbursts with durations similar to a single orbital period. In an effort to understand these signatures, we present and analyze a hydrodynamic model of unstable mass transfer from a giant-star donor onto a more compact accretor in a binary system. Our simulation begins with mass transfer at the Roche limit separation and traces a phase of runaway decay leading up to the plunge of the accretor within the envelope of the donor. We characterize the fluxes of mass and angular momentum through the system and show that the orbital evolution can be reconstructed from measurements of these quantities. The morphology of outflow from the binary changes significantly as the binary orbit tightens. At wide separations, a thin stream of relatively high-entropy gas trails from the outer Lagrange points. As the orbit tightens, the orbital motion desynchronizes from the donor's rotation, and low-entropy ejecta trace a broad fan of largely ballistic trajectories. An order-of-magnitude increase in mass ejection rate accompanies the plunge of the accretor with the envelope of the donor. We argue that this transition marks the precursor-to-outburst transition observed in stellar coalescence transients., Comment: Revised following peer-review. ApJ accepted. Animated version of Figure 5 will be available via the Journal's online publication

We present Athena++ grid-based, hydrodynamic simulations of accretion onto Sagittarius A* via the stellar winds of the $\sim 30$ Wolf-Rayet stars within the central parsec of the galactic center. These simulations span $\sim$ 4 orders of magnitude in radius, reaching all the way down to 300 gravitational radii of the black hole, $\sim 32$ times further in than in previous work. We reproduce reasonably well the diffuse thermal X-ray emission observed by Chandra in the central parsec. The resulting accretion flow at small radii is a superposition of two components: 1) a moderately unbound, sub-Keplerian, thick, pressure-supported disc that is at most (but not all) times aligned with the clockwise stellar disc, and 2) a bound, low-angular momentum inflow that proceeds primarily along the southern pole of the disc. We interpret this structure as a natural consequence of a few of the innermost stellar winds dominating accretion, which produces a flow with a broad distribution of angular momentum. Including the star S2 in the simulation has a negligible effect on the flow structure. Extrapolating our results from simulations with different inner radii, we find an accretion rate of $\sim$ a few $\times 10^{-8} M_\odot$/yr at the horizon scale, consistent with constraints based on modeling the observed emission of Sgr A*. The flow structure found here can be used as more realistic initial conditions for horizon scale simulations of Sgr A*., Comment: 20 pages, 19 figures. See http://w.astro.berkeley.edu/~sressler/animations.html for animations. Accepted by MNRAS

We perform global three-dimensional MHD simulations of unstratified accretion disks in cataclysmic variables (CVs). By including mass inflow via an accretion stream, we are able to evolve the disk to a steady state. We investigate the relative importance of spiral shocks and the magnetorotational instability (MRI) in driving angular momentum transport and how each depend on the geometry and strength of seed magnetic field and the Mach number of the disk (where Mach number is ratio of the azimuthal velocity and the sound speed of gas). We use a locally isothermal equation of state and adopt temperature profiles that are consistent with CV disk observations. Our results indicate that the relative importance of spiral shocks and MRI in driving angular momentum transport is controlled by the gas Mach number and the seed magnetic field strength. MRI and spiral shocks provide comparable efficiency of angular momentum transport when the disk Mach number is around 10 and the seed magnetic field has plasma $\beta=400$ (where $\beta$ is ratio of gas pressure and magnetic pressure). The MRI dominates whenever the seed field strength, or the disk Mach number, is increased. Among all of our simulations, the effective viscosity parameter $\alpha_{eff} \sim 0.016-0.1$ after MRI saturates and the disk reaches steady state. Larger values of $\alpha_{eff}$ are favored when the seed magnetic field has vertical components or the flow has stronger magnetization ($1/\beta$). Our models all indicate that the role of MRI in driving angular momentum transport thus mass accretion in CV disks is indispensable, especially in cool disks with weak spiral shocks., Comment: 22 pages, 15 figures

We report new global ideal MHD simulations for thin accretion disks (with thermal scale height H/R=0.1 and 0.05) threaded by net vertical magnetic fields. Our computations span three orders of magnitude in radius, extend all the way to the pole, and are evolved for more than one viscous time over the inner decade in radius. Static mesh refinement is used to properly resolve MRI. We find that:(1) inward accretion occurs mostly in the upper magnetically dominated regions of the disk, similar to the predictions from some previous analytical work and the "coronal accretion" in previous GRMHD simulations. Rapid inflow in the upper layers combined with slow outflow at the midplane creates strong $R\phi$ and $z\phi$ stresses in the mean field; the vertically integrated $\alpha\sim 0.5-1$ when the initial field has $\beta_{0}=10^3$ at the midplane. (2) A quasi-static global field geometry is established in which flux transport by inflows at the surface is balanced by turbulent diffusion. The field is strongly pinched inwards at the surface. A steady-state advection-diffusion model, with turbulent magnetic Prandtl number of order unity, reproduces this geometry well. (3) Weak unsteady disk winds are launched at $z/R\sim1$ with the Alfven radius $R_{A}/R_{0}\sim3$. Although the wind is episodic, the time averaged properties are well described by steady wind theory. Wind is not efficient at transporting angular momentum. Even with $\beta_{0}=10^3$, only 5% of the angular momentum transport is driven by torque from the wind, and the wind mass flux from the inner decade of radius is only $\sim$ 0.4% of the mass accretion rate. With weaker fields or thinner disks, the wind contributes even less. (4) Most of the disk accretion is driven by the $R\phi$ stress from the MRI and global magnetic fields. Our simulations have many applications to astrophysical accretion disk systems., Comment: Two additional figures, accepted by the AAS Journals

The interaction of a satellite with a gaseous disk results in the excitation of spiral density waves which remove angular momentum from the orbit. In addition, if the orbit is not coplanar with the disk, three-dimensional effects will excite bending and eccentricity waves. We perform three-dimensional hydrodynamic simulations to study nonlinear disk-satellite interaction in inviscid protoplanetary disks for a variety of orbital inclinations from $0^\circ$ to $180^\circ$. It is well known that three-dimensional effects are important even for zero inclination. In this work we (1) show that for planets with small inclinations (as in the Solar system), effects such as the total torque and migration rate strongly depend on the inclination and are significantly different (about 2.5 times smaller) from the two-dimensional case, (2) give formulae for the migration rate, inclination damping, and precession rate of planets with different inclination angles in disk with different scale heights, and (3) present the observational signatures of a planet on an inclined orbit with respect to the protoplanetary disk. For misaligned planets we find good agreement with linear theory in the limit of small inclinations, and with dynamical friction estimates for intermediate inclinations. We find that in the latter case, the dynamical friction force is not parallel to the relative planetary velocity. Overall, the derived formulae will be important for studying exoplanets with obliquity., Comment: 12 pages, 9 figures, accepted by MNRAS

The early optical afterglow emission of several gamma-ray bursts (GRBs) shows a high linear polarization degree (PD) of tens of percent, suggesting an ordered magnetic field in the emission region. The light curves are consistent with being of a reverse shock (RS) origin. However, the magnetization parameter, $\sigma$, of the outflow is unknown. If $\sigma$ is too small, an ordered field in the RS may be quickly randomized due to turbulence driven by various perturbations so that the PD may not be as high as observed. Here we use the "Athena++" relativistic MHD code to simulate a relativistic jet with an ordered magnetic field propagating into a clumpy ambient medium, with a focus on how density fluctuations may distort the ordered magnetic field and reduce PD in the RS emission for different $\sigma$ values. For a given density fluctuation, we discover a clear power-law relationship between the relative PD reduction and the $\sigma$ value of the outflow. Such a relation may be applied to estimate $\sigma$ of the GRB outflows using the polarization data of early afterglows., Comment: 6 pages, 4 figures, matching the publication version of ApJL

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., 13 pages, 12 figures, submitted to ApJ

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_\bullet-T_{\rm x}$, followed by $M_\bullet-L_{\rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_\bullet-\sigma_{\rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_\bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions)., Comment: 40 pages, 27 figures, 3 tables; ApJ accepted version (minor revision) - Reviewer: "The authors are to be congratulated on a very interesting paper that deserves to be read by anyone interested in the connection between the SMBH and host galaxy."

We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5\times 10^8M_{\odot}$ black hole with accretion rates varying from $\sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $\sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $\sim 0.1-0.4c$. When the accretion rate is smaller than $\sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $\sim 5\%-7\%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1\%$. We always find the kinetic energy luminosity associated with the outflow is only $\sim 15\%-30\%$ of the radiative luminosity. The mass flux lost in the outflow is $\sim 15\%-50\%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks., 21 pages, 17 figures, submitted to ApJ, movies available here https://goo.gl/5yR3ig

Circumplanetary disks (CPDs) control the growth of planets, supply material for satellites to form, and provide observational signatures of young forming planets. We have carried out two dimensional hydrodynamical simulations with radiative cooling to study CPDs, and suggested a new mechanism to drive the disk accretion. Two spiral shocks are present in CPDs, excited by the central star. We find that spiral shocks can at least contribute to, if not dominate the angular momentum transport and energy dissipation in CPDs. Meanwhile, dissipation and heating by spiral shocks have a positive feedback on shock-driven accretion itself. As the disk is heated up by spiral shocks, the shocks become more open, leading to more efficient angular momentum transport. This shock driven accretion is, on the other hand, unsteady on a timescale of months/years due to production and destruction of vortices in disks. After being averaged over time, a quasi-steady accretion is reached from the planet's Hill radius all the way to the planet surface, and the disk $\alpha$-coefficient characterizing angular momentum transport due to spiral shocks is $\sim$0.001-0.02. The disk surface density ranges from 10 to 1000 g cm$^{-2}$ in our simulations, which is at least 3 orders of magnitude smaller than the "minimum mass sub-nebula" model used to study satellite formation; instead it is more consistent with the "gas-starved" satellite formation model. Finally, we calculate the millimeter flux emitted by CPDs at ALMA and EVLA wavelength bands and predict the flux for several recently discovered CPD candidates, which suggests that ALMA is capable of discovering these accreting CPDs., Comment: 19 pages, 16 figures, and 3 tables, ApJ in press

To understand the conditions under which dense, molecular gas is able to form within a galaxy, we post-process a series of three-dimensional galactic-disk-scale simulations with ray-tracing based radiative transfer and chemical network integration to compute the equilibrium chemical and thermal state of the gas. In performing these simulations we vary a number of parameters, such as the ISRF strength, vertical scale height of stellar sources, cosmic ray flux, to gauge the sensitivity of our results to these variations. Self-shielding permits significant molecular hydrogen (H2) abundances in dense filaments around the disk midplane, accounting for approximately ~10-15% of the total gas mass. Significant CO fractions only form in the densest, n>~10^3 cm^-3, gas where a combination of dust, H2, and self-shielding attenuate the FUV background. We additionally compare these ray-tracing based solutions to photochemistry with complementary models where photo-shielding is accounted for with locally computed prescriptions. With some exceptions, these local models for the radiative shielding length perform reasonably well at reproducing the distribution and amount of molecular gas as compared with a detailed, global ray tracing calculation. Specifically, an approach based on the Jeans Length with a T=40K temperature cap performs the best in regards to a number of different quantitative measures based on the H2 and CO abundances., 21 Pages, 15 figures. Submitted to MNRAS. Comments welcome

We present results from the first global 3D MHD simulations of accretion disks in Cataclysmic Variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared to that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component to the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states and both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. 1) Linear wave dispersion relation fits the pitch angles of spiral density waves very well. 2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. 3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate alpha_{eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. 4) Using the mass accretion rate measured in our simulations to directly measure alpha defined in standard thin-disk theory, we find 0.02 < alpha_{eff} < 0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae (DNe). In this regime the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime, and therefore future models with more realistic thermodynamics and non-ideal MHD are warranted., 23 pages, 15 figures, accepted by ApJ

Protoplanetary disks (PPDs) are widely believed to be turbulent as a result of the magnetorotational instability (MRI). We perform magnetohydrodynamical simulations of PPDs that for the first time, take into account both Ohmic resistivity and ambipolar diffusion in a self-consistent manner. We show that in the inner region of PPDs that corresponds the habitable zone, the MRI is completely suppressed due to the interplay between magnetic field and ambipolar diffusion. The gas in this region is laminar throughout the entire vertical extent of the disk. Instead of MRI-driven accretion, a strong magnetocentrifugal wind is launched that efficiently carries away disk angular momentum. A physical wind geometry requires the presence of a strong current layer that is offset from the disk midplane where horizontal magnetic fields flip. We show that the entire accretion flow proceeds through this strong current layer. The non-turbulent nature of the gas flow strongly favors the habitable zone as the site for planetesimal formation, and has important implications for their subsequent growth into terrestrial planets.

In a hot, dilute, magnetized accretion flow, the electron mean-free path can be much greater than the Larmor radius, thus thermal conduction is anisotropic and along magnetic field lines. In this case, if the temperature decreases outward, the flow may be subject to a buoyancy instability (the magnetothermal instability, or MTI). The MTI amplifies the magnetic field, and aligns field lines with the radial direction. If the accretion flow is differentially rotating, the magnetorotational instability (MRI) may also be present. Using two-dimensional, time-dependent magnetohydrodynamic simulations, we investigate the interaction between these two instabilities. We use global simulations that span over two orders of magnitude in radius, centered on the region around the Bondi radius where the infall time of gas is longer than the growth time of both the MTI and MRI. Significant amplification of the magnetic field is produced by both instabilities, although we find that the MTI primarily amplifies the radial component, and the MRI primarily the toroidal component, of the field, respectively. Most importantly, we find that if the MTI can amplify the magnetic energy by a factor $F_t$, and the MRI by a factor $F_r$, then when the MTI and MRI are both present, the magnetic energy can be amplified by a factor of $F_t \cdot F_r$. We therefore conclude that amplification of the magnetic energy by the MTI and MRI operates independently. We also find that the MTI contributes to the transport of angular momentum, because radial motions induced by the MTI increase the Maxwell (by amplifying the magnetic field) and Reynolds stresses. Finally, we find that thermal conduction decreases the slope of the radial temperature profile. The increased temperature near the Bondi radius decreases the mass accretion rate.

The X-ray signal from hydrogen-rich supernovae (SNe II) in the first tens to hundreds of days after the shock breakout encodes important information about the circumstellar material (CSM) surrounding their progenitors before explosion. In this study, we describe a way to generate the SN II X-ray light curves from hydrodynamical simulations performed with the code Athena++, using the X-ray package XSPEC. In addition, we employ a radiation diffusion hydrodynamic code SNEC for generating the optical light curves in different bands. In this numerical setup, we model the X-ray and optical emission from a set of progenitor models, consisting of either two (red supergiant + low density steady wind), or three (red supergiant + dense CSM + low density steady wind) components. We vary the density in the wind and the slope in the CSM to see how these parameters influence the resulting X-ray and optical light curves. Among our models, we identify one that is able to roughly reproduce both optical and X-ray data of the well observed SN 2013ej. In order to achieve this, the slope of the dense CSM in this model should be steeper than the one of a steady wind ($\rho\propto r^{-2}$), and closer to $\rho\propto r^{-5}$. On the other hand, we show that too steep and extended CSM profiles may produce excessive X-ray emission in the first few tens of days, up to a few orders of magnitude larger than observed. We conclude that ability to reproduce the observed X-ray signal from SNe~II together with their optical light curves is crucial in establishing the validity of different CSM models., Comment: 17 pages, 18 figures

Magnetohydrodynamic instabilities play an important role in accretion disks systems. Besides the well-known effects of the magnetorotational instability (MRI), the Parker-Rayleigh-Taylor instability (PRTI) also arises as an important mechanism to help in the formation of the coronal region around an accretion disk and in the production of magnetic reconnection events similar to those occurring in the solar corona. In this work, we have performed three-dimensional magnetohydrodynamical (3D-MHD) shearing-box numerical simulations of accretion disks with an initial stratified density distribution and a strong azimuthal magnetic field with a ratio between the thermal and magnetic pressures of the order of unity. This study aimed at verifying the role of these instabilities in driving fast magnetic reconnection in turbulent accretion disk/corona systems. As we expected, the simulations showed an initial formation of large-scale magnetic loops due to the PRTI followed by the development of a nearly steady-state turbulence driven by both instabilities. In this turbulent environment, we have employed an algorithm to identify the presence of current sheets produced by the encounter of magnetic flux ropes of opposite polarity in the turbulent regions of both the corona and the disk. We computed the magnetic reconnection rates in these locations obtaining average reconnection velocities in Alfv\'en speed units of the order of $0.13 \pm 0.09$ in the accretion disk and $0.17 \pm 0.10$ in the coronal region (with mean peak values of order of $0.2$), which are consistent with the predictions of the theory of turbulence-induced fast reconnection., Comment: 23 pages, 14 figures. Accepted for publication in ApJ

The nonlinear saturation of the magnetorotational instability (MRI) is best studied through numerical MHD simulations. Recent results of simulations that adopt the local shearing box approximation, and fully global models that follow the entire disk, are described. Outstanding issues remain, such as a first-principles understanding of the dynamo processes that control saturation with no net magnetic flux. Important directions for future work include a better understanding of basic plasma processes, such as reconnection, dissipation, and particle acceleration, in the MHD turbulence driven by the MRI.

We examine the effects of density stratification on magnetohydrodynamic turbulence driven by the magnetorotational instability in local simulations that adopt the shearing box approximation. Our primary result is that, even in the absence of explicit dissipation, the addition of vertical gravity leads to convergence in the turbulent energy densities and stresses as the resolution increases, contrary to results for zero net flux, unstratified boxes. The ratio of total stress to midplane pressure has a mean of ~0.01, although there can be significant fluctuations on long (>~50 orbit) timescales. We find that the time averaged stresses are largely insensitive to both the radial or vertical aspect ratio of our simulation domain. For simulations with explicit dissipation, we find that stratification extends the range of Reynolds and magnetic Prandtl numbers for which turbulence is sustained. Confirming the results of previous studies, we find oscillations in the large scale toroidal field with periods of ~10 orbits and describe the dynamo process that underlies these cycles., 13 pages, 18 figures, submitted to ApJ

We perform two-dimensional and three-dimensional radiation hydrodynamic simulations to study cold clouds accelerated by radiation pressure on dust in the environment of rapidly star-forming galaxies dominated by infrared flux. We utilize the reduced speed of light approximation to solve the frequency-averaged, time-dependent radiative transfer equation. We find that radiation pressure is capable of accelerating the clouds to hundreds of kilometers per second while remaining dense and cold, consistent with observations. We compare these results to simulations where acceleration is provided by entrainment in a hot wind, where the momentum injection of the hot flow is comparable to the momentum in the radiation field. We find that the survival time of the cloud accelerated by the radiation field is significantly longer than that of a cloud entrained in a hot outflow. We show that the dynamics of the irradiated cloud depends on the initial optical depth, temperature of the cloud, and the intensity of the flux. Additionally, gas pressure from the background may limit cloud acceleration if the density ratio between the cloud and background is $\lesssim 10^{2}$. In general, a 10 pc-scale optically thin cloud forms a pancake structure elongated perpendicular to the direction of motion, while optically thick clouds form a filamentary structure elongated parallel to the direction of motion. The details of accelerated cloud morphology and geometry can also be affected by other factors, such as the cloud lengthscale, the reduced speed of light approximation, spatial resolution, initial cloud structure, and the dimensionality of the run, but these have relatively little affect on the cloud velocity or survival time., 22 pages, 21 figures, 2 table, Accepted for publication in ApJ

We present a set of numerical simulations addressing the effects of magnetic field strength and orientation on the flow-driven formation of molecular clouds. Fields perpendicular to the flows sweeping up the cloud can efficiently prevent the formation of massive clouds but permit the build-up of cold, diffuse filaments. Fields aligned with the flows lead to substantial clouds, whose degree of fragmentation and turbulence strongly depends on the background field strength. Adding a random field component leads to a "selection effect" for molecular cloud formation: high column densities are only reached at locations where the field component perpendicular to the flows is vanishing. Searching for signatures of colliding flows should focus on the diffuse, warm gas, since the cold gas phase making up the cloud will have lost the information about the original flow direction because the magnetic fields redistribute the kinetic energy of the inflows., Comment: 11 pages, 9 figures, accepted by ApJ

We describe a numerical algorithm based on Godunov methods for integrating the equations of compressible magnetohydrodynamics (MHD) in multidimensions. It combines a simple, dimensionally-unsplit integration method with the constrained transport (CT) discretization of the induction equation to enforce the divergence-free constraint. We present the results of a series of fully three-dimensional tests which indicate the method is second-order accurate for smooth solutions in all MHD wave families, and captures shocks, contact and rotational discontinuities well. However, it is also more diffusive than other more complex unsplit integrators combined with CT. Thus, the primary advantage of the method is its simplicity. It does not require a characteristic tracing step to construct interface values for the Riemann solver, it is straightforward to extend with additional physics, and it is suitable for use with nested and adaptive meshes. The method is implemented as one of two dimensionally unsplit MHD integrators in the Athena code, which is freely available for download from the web.

We study the effects of anisotropic thermal conduction on magnetized spherical accretion flows using global axisymmetric MHD simulations. In low collisionality plasmas, the Bondi spherical accretion solution is unstable to the magnetothermal instability (MTI). The MTI grows rapidly at large radii where the inflow is subsonic. For a weak initial field, the MTI saturates by creating a primarily radial magnetic field, i.e., by aligning the field lines with the background temperature gradient. The saturation is quasilinear in the sense that the magnetic field is amplified by a factor of $\sim 10-30$ independent of the initial field strength (for weak fields). In the saturated state, the conductive heat flux is much larger than the convective heat flux, and is comparable to the field-free (Spitzer) value (since the field lines are largely radial). The MTI by itself does not appreciably change the accretion rate $\dot M$ relative to the Bondi rate $\dot M_B$. However, the radial field lines created by the MTI are amplified by flux freezing as the plasma flows in to small radii. Oppositely directed field lines are brought together by the converging inflow, leading to significant resistive heating. When the magnetic energy density is comparable to the gravitational potential energy density, the plasma is heated to roughly the virial temperature; the mean inflow is highly subsonic; most of the energy released by accretion is transported to large radii by thermal conduction; and the accretion rate $\dot M \ll \dot M_B$. The predominantly radial magnetic field created by the MTI at large radii in spherical accretion flows may account for the stable Faraday rotation measure towards Sgr A* in the Galactic Center.

We study the probability distribution function (PDF) of the mass density in simulations of supersonic turbulence with properties appropriate for molecular clouds. For this study we use Athena, a new higher-order Godunov code. We find there are surprisingly similar relationships between the mean of the time-averaged PDF and the turbulent Mach number for driven hydrodynamic and strong-field MHD turbulence. There is, however, a large scatter about these relations, indicating a high level of temporal and spatial variability in the PDF. Thus, the PDF of the mass density is unlikely to be a good measure of magnetic field strength. We also find the PDF of decaying MHD turbulence deviates from the mean-Mach relation found in the driven case. This implies that the instantaneous Mach number alone is not enough to determine the statistical properties of turbulence that is out of equilibrium. The scatter about the mean-Mach relation for driven turbulence, along with the large departure of decaying turbulence PDFs from those of driven turbulence, may illuminate one factor contributing to the large observed cloud-to-cloud variation in the star formation rate per solar mass., Comment: 4 pages, 2 figures, published in ApJL; corrected sign in eqn. (4)