Search

Hot accretion flows contain collisionless plasmas that are believed to be capable of accelerating particles to very high energies, as a result of turbulence generated by the magnetorotational instability (MRI). We conduct unstratified shearing-box simulations of the MRI turbulence in ideal magnetohydrodynamic (MHD), and inject energetic (relativistic) test particles in simulation snapshots to conduct a detailed investigation on particle diffusion and stochastic acceleration. We consider different amount of net vertical magnetic flux to achieve different disk magnetizations levels at saturated states, with sufficiently high resolution to resolve the gyro-radii ($R_g$) of most particles. Particles with large $R_g$ ($\gtrsim0.03$ disk scale height $H$) show spatial diffusion coefficients of $\sim30$ and $\sim5$ times Bohm values in the azimuthal and poloidal directions, respectively. We further measure particle momentum diffusion coefficient $D(p)$ by applying the Fokker-Planck equation to particle momentum evolution. For these particles, contribution from turbulent fluctuations scales as $D(p)\propto p$, and shear acceleration takes over when $R_g\gtrsim0.1H$, characterized by $D(p)\propto p^3$. For particles with smaller $R_g$ ($\lesssim0.03H$), their spatial diffusion coefficients roughly scale as $\sim p^{-1}$, and show evidence of $D(p)\propto p^2$ scaling in momentum diffusion but with large uncertainties. We find that multiple effects contribute to stochastic acceleration/deceleration, and the process is also likely affected by intermittency in the MRI turbulence. We also discuss the potential of accelerating PeV cosmic-rays in hot accretion flows around supermassive black holes., Comment: 22 pages, 18 figures, accepted for publication in MNRAS

The vertical shear instability (VSI) is a robust phenomenon in irradiated protoplanetary disks (PPDs). While there is extensive literature on the VSI in the hydrodynamic limit, PPDs are expected to be magnetized and their extremely low ionization fractions imply that non-ideal magneto-hydrodynamic (MHD) effects should be properly considered. To this end, we present linear analyses of the VSI in magnetized disks with Ohmic resistivity. We primarily consider toroidal magnetic fields, which are likely to dominate the field geometry in PPDs. We perform vertically global and radially local analyses to capture characteristic VSI modes with extended vertical structures. To focus on the effect of magnetism, we use a locally isothermal equation of state. We find that magnetism provides a stabilizing effect to dampen the VSI, with surface modes, rather than body modes, being the first to vanish with increasing magnetization. Subdued VSI modes can be revived by Ohmic resistivity, where sufficient magnetic diffusion overcome magnetic stabilization, and hydrodynamic results are recovered. We also briefly consider poloidal fields to account for the magnetorotational instability (MRI), which may develop towards surface layers in the outer parts of PPDs. The MRI grows efficiently at small radial wavenumbers, in contrast to the VSI. When resistivity is considered, we find the VSI dominates over the MRI for Ohmic Els\"{a}sser numbers $\lesssim 0.09$ at plasma beta parameter $\beta_Z \sim 10^4$., Comment: 17 pages, 11 figures, accepted for publication in MNRAS

Convection has been discussed in the field of accretion discs for several decades, both as a means of angular momentum transport and also because of its role in controlling discs' vertical structure via heat transport. If the gas is sufficiently ionized and threaded by a weak magnetic field, convection might interact in non-trivial ways with the magnetorotational instability (MRI). Recently, vertically stratified local simulations of the MRI have reported considerable variation in the angular momentum transport, as measured by the stress to thermal pressure ratio $\alpha$, when convection is thought to be present. Although MRI turbulence can act as a heat source for convection, it is not clear how the instabilities will interact dynamically. Here we aim to investigate the interplay between the two instabilities in controlled numerical experiments, and thus isolate the generic features of their interaction. We perform vertically stratified, 3D MHD shearing box simulations with a perfect gas equation of state with the conservative, finite-volume code PLUTO. We find two characteristic outcomes of the interaction between the two instabilities: straight MRI and MRI/convective cycles, with the latter exhibiting alternating phases of convection-dominated flow (during which the turbulent transport is weak) and MRI-dominated flow. During the latter phase we find that $\alpha$ is enhanced by nearly an order of magnitude, reaching peak values of $\sim 0.08$. In addition, we find that convection in the non-linear phase takes the form of large-scale and oscillatory convective cells. Convection can also help the MRI persist to lower Rm than it would otherwise do. Finally we discuss how our results help interpret simulations of Dwarf Novae., Comment: Accepted for publication in MNRAS (23 pages, 14 figures, 1 table)

We investigate whether the regular Galilean satellites could have formed in the dead zone of a circumplanetary disc. A dead zone is a region of weak turbulence in which the magnetorotational instability (MRI) is suppressed, potentially an ideal environment for satellite formation. With the grid-based hydrodynamic code, FARGO3D, we examine the evolution of a circumplanetary disc model with a dead zone. Material accumulates in the dead zone of the disc leading to a higher total mass and but a similar temperature profile compared to a fully turbulent disc model. The tidal torque increases the rate of mass transport through the dead zone leading to a steady state disc with a dead zone that does not undergo accretion outbursts. We explore a range of disc, dead zone and mass inflow parameters and find that the maximum mass of the disc is around 0.001 MJ . Since the total solid mass of such a disc is much lower, we find that there is not sufficient material in the disc for in situ formation of the Galilean satellites and that external supplement is required., Comment: 10 pages, 6 figures

We explore the pulsationally driven orbital mass ejection mechanism for Be star disc formation using isothermal, 3D magnetohydrodynamic (MHD) and hydrodynamic simulations. Non-radial pulsations are added to a star rotating at 95\% of critical as an inner boundary condition that feeds gas into the domain. In MHD, the initial magnetic field within the star is weak. The hydrodynamics simulation has limited angular momentum transport, resulting in repeating cycles of mass accumulation into a rotationally-supported disc at small radii followed by fall-back onto the star. The MHD simulation, conversely, has efficient (Maxwell $\alpha_{\rm M}$ $\sim$ 0.04) angular momentum transport provided by both of turbulent and coherent magnetic fields; a slowly decreting midplane driven by the magnetorotational instability and a supersonic wind on the surface of the disc driven by global magnetic torques. The angle and time-averaged properties near the midplane agree reasonably well with a 1D viscous decretion disc model with a modified $\tilde\alpha=0.5$, in which the gas transitions from a subsonic thin disc to a supersonic spherical wind at the critical point. 1D models, however, cannot capture the multi-phase decretion/angular structure seen in our simulations. Our results demonstrate that, at least under certain conditions, non-radial pulsations on the surface of a rapidly rotating, weakly magnetized star can drive a Keplerian disc with the basic properties of the viscous decretion disc paradigm, albeit coupled to a laminar wind away from the midplane. Future modeling of Be star discs should consider the possible existence of such a surface wind., Comment: Accepted by MNRAS. 16 pages, 10 figures

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two and three dimensional numerical modelling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical ($z$) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal ($rz$) plane, driving magnetic reconnection when the Lundquist number $S > 3 \times 10^2$, making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at $S = 5 \times 10^3$, self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When non-axisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets further unstable, and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation., 12 pages, 6 figures, Accepted for publication to ApJ Letters

In the present study, we addressed the question of how the growth rate of the magnetorotational instability is modified when the radial component of the stellar dipole magnetic field is taken into account in addition to the vertical component. Considering a fiducial radius in the disc where diamagnetic currents are pronounced, we carried out a linear stability analysis to obtain the growth rates of the magnetorotational instability for various parameters such as the ratio of the radial-to-vertical component and the gradient of the magnetic field, the Alfvenic Mach number and the diamagnetization parameter. Our results show that the interaction between the diamagnetic current and the radial component of the magnetic field increases the growth rate of the magnetorotational instability and generates a force perpendicular to the disc plane which may induce a torque. It is also shown that considering the radial component of the magnetic field and taking into account a radial gradient in the vertical component of the magnetic field causes an increase in the magnitudes of the growth rates of both the axisymmetric ($m=0$) and the non-axisymmetric ($m=1$) modes., Comment: 9 pages, 5 figures, published in MNRAS, 2020, vol. 491, pages 5481-5488

We investigate acceleration and propagation processes of high-energy particles inside hot accretion flows. The magnetorotational instability (MRI) creates turbulence inside accretion flows, which triggers magnetic reconnection and may produce non-thermal particles. They can be further accelerated stochastically by the turbulence. To probe the properties of such relativistic particles, we perform magnetohydrodynamic simulations to obtain the turbulent fields generated by the MRI, and calculate orbits of the high-energy particles using snapshot data of the MRI turbulence. We find that the particle acceleration is described by a diffusion phenomenon in energy space with a diffusion coefficient of the hard-sphere type: $D_\epsilon\propto \epsilon^2$, where $\epsilon$ is the particle energy. Eddies in the largest scale of the turbulence play a dominant role in the acceleration process. On the other hand, the stochastic behaviour in configuration space is not usual diffusion but superdiffusion: the radial displacement increases with time faster than that in the normal diffusion. Also, the magnetic field configuration in the hot accretion flow creates outward bulk motion of high-energy particles. This bulk motion is more effective than the diffusive motion for higher energy particles. Our results imply that typical active galactic nuclei that host hot accretion flows can accelerate CRs up to $\epsilon\sim 0.1-10$ PeV., Comment: 17 pages, 16 figures, Figure 4 is corrected, conclusions are unchanged, published in MNRAS

We explore in situ formation and subsequent evolution of close-in super-Earths and mini-Neptunes. We adopt a steady-state inner protoplanetary gas disc structure that arises from viscous accretion due to the magneto-rotational instability (MRI). We consider the evolution of dust in the inner disc, including growth, radial drift and fragmentation, and find that dust particles that radially drift into the inner disc fragment severely due to the MRI-induced turbulence. This result has two consequences: (1) radial drift of grains within the inner disc is quenched, leading to an enhancement of dust in the inner regions which scales as dust-to-gas-mass-flux-ratio at ~1 AU; (2) however, despite this enhancement, planetesimal formation is impeded by the small grain size. Nevertheless, assuming that planetary cores are present in the inner disc, we then investigate the accretion of atmospheres onto cores and their subsequent photoevaporation. We then compare our results to the observed exoplanet mass-radius relationship. We find that: (1) the low gas surface densities and high temperatures in the inner disc reduce gas accretion onto cores compared to the minimum mass solar nebula, preventing the cores from growing into hot Jupiters, in agreement with the data; (2) however, our predicted envelope masses are still typically larger than observed ones. Finally, we sketch a qualitative picture of how grains may grow and planetesimals form in the inner disc if grain effects on the ionization levels and the MRI and the back-reaction of the dust on the gas (both neglected in our calculations) are accounted for., Comment: Accepted for publication in MNRAS

Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling of such sources by means of general relativistic magnetohydrodynamical (GRMHD) simulations in arbitrary stationary spacetimes. New additions were required to guarantee an accurate evolution of the magnetic field when small and large scales are captured simultaneously. We discuss the adaptive mesh refinement (AMR) techniques employed in BHAC, essential to keep several problems computationally tractable, as well as staggered-mesh-based constrained transport (CT) algorithms to preserve the divergence-free constraint of the magnetic field, including a general class of prolongation operators for face-allocated variables compatible with them. Through several standard tests, we show that the choice of divergence-control method can produce qualitative differences in simulations of scientifically relevant accretion problems. We demonstrate the ability of AMR to reduce the computational costs of accretion simulations while sufficiently resolving turbulence from the magnetorotational instability. In particular, we describe a simulation of an accreting Kerr black hole in Cartesian coordinates using AMR to follow the propagation of a relativistic jet while self-consistently including the jet engine, a problem set up-for which the new AMR implementation is particularly advantageous. The CT methods and AMR strategies discussed here are being employed in the simulations performed with BHAC used in the generation of theoretical models for the Event Horizon Telescope Collaboration., Comment: 21 pages, 16 figures, accepted for publication in Astronomy & Astrophysics. Link to www.bhac.science will soon be active

Many high-energy astrophysical systems contain magnetized collisionless plasmas with relativistic particles, in which turbulence can be driven by an arbitrary mixture of solenoidal and compressive motions. For example, turbulence in hot accretion flows may be driven solenoidally by the magnetorotational instability or compressively by spiral shock waves. It is important to understand the role of the driving mechanism on kinetic turbulence and the associated particle energization. In this work, we compare particle-in-cell simulations of solenoidally driven turbulence with similar simulations of compressively driven turbulence. We focus on plasma that has an initial beta of unity, relativistically hot electrons, and varying ion temperature. Apart from strong large-scale density fluctuations in the compressive case, the turbulence statistics are similar for both drives, and the bulk plasma is described reasonably well by an isothermal equation of state. We find that nonthermal particle acceleration is more efficient when turbulence is driven compressively. In the case of relativistically hot ions, both driving mechanisms ultimately lead to similar power-law particle energy distributions, but over a different duration. In the case of non-relativistic ions, there is significant nonthermal particle acceleration only for compressive driving. Additionally, we find that the electron-to-ion heating ratio is less than unity for both drives, but takes a smaller value for compressive driving. We demonstrate that this additional ion energization is associated with the collisionless damping of large-scale compressive modes via perpendicular electric fields., 29 pages, 28 figures, accepted for publication in ApJ

Aims. Accretion and luminosity bursts triggered by three distinct mechanisms: the magnetorotational instability in the inner disk regions, clump infall in gravitationally fragmented disks and close encounters with an intruder star, were studied to determine the disk kinematic characteristics that can help to distinguish between these burst mechanisms. Methods. Numerical hydrodynamics simulations in the thin-disk limit were employed to model the bursts in disk environments that are expected for each burst mechanism. Results. We found that the circumstellar disks featuring accretion bursts can bear kinematic features that are distinct for different burst mechanisms, which can be useful when identifying the burst origin. The disks in the stellar encounter and clump-infall models are characterized by tens of percent deviations from the Keplerian rotation, whie the disks in the MRI models are characterized only a few percent deviation, which is mostly caused by the gravitational instability that fuels the MRI bursts. Velocity channel maps also show distinct kinks and wiggles, which are caused by gas disk flows that are peculiar to each considered burst mechanism. The deviations of velocity channels in the burst-hosting disks from a symmetric pattern typical of Keplerian disks are strongest for the clump-infall and collision models, and carry individual features that may be useful for the identification of the corresponding burst mechanism. The considered burst mechanisms produce a variety of light curves with the burst amplitudes varying in the \Delta m=2.5-3.7 limits, except for the clump-infall model where \Delta m can reach 5.4, although the derived numbers may be affected by a small sample and boundary conditions. Conclusions. Burst triggering mechanisms are associated with distinct kinematic features in the burst-hosting disks that may be used for their identification. Abridged., Comment: Accepted for publication by Astronomy & Astrophysics

DW Cnc is an intermediate polar which has previously been observed in both high and low states. Observations of the high state of DW Cnc have previously revealed a spin period at ~ 38.6 min, however observations from the 2018/19 low state showed no evidence of the spin period. We present results from our analysis of 12 s cadence photometric data collected by NGTS of DW Cnc during the high state which began in 2019. Following the previously reported suppression of the spin period signal we identify the return of this signal during the high state, consistent with previous observations of it. We identify this as the restarting of accretion during the high state. We further identified three short outbursts lasting ~ 1 d in DW Cnc with a mean recurrence time of ~ 60 d and an amplitude of ~ 1 mag. These are the first outbursts identified in DW Cnc since 2008. Due to the short nature of these events we identify them not as a result of accretion instabilities but instead either from instabilities originating from the interaction of the magnetorotational instability in the accretion disc and the magnetic field generated by the white dwarf or the result of magnetic gating., Comment: Accepted to MNRAS; 8 pages, 4 figues

We investigate the possibility of the growth of magnetorotational instability (MRI) in disks around Class 0 protostars. We construct a disk model and calculate the chemical reactions of neutral and charged atoms, molecules and dust grains to derive the abundance of each species and the ionization degree of the disk. Then, we estimate the diffusion coefficients of non-ideal magnetohydrodynamics effects such as ohmic dissipation, ambipolar diffusion and the Hall effect. Finally, we evaluate the linear growth rate of MRI in each area of the disk. We investigate the effect of changes in the strength and direction of the magnetic field in our disk model and we adopt four different dust models to investigate the effect of dust size distribution on the diffusion coefficients. Our results indicate that an MRI active region possibly exists with a weak magnetic field in a region far from the protostar where the Hall effect plays a role in the growth of MRI. On the other hand, in all models the disk is stable against MRI in the region within $, Accepted for publication in MNRAS. 43 pages, 12 figures

In the early stages of a protoplanetary disk, when its mass is a significant fraction of its star's, turbulence generated by gravitational instability (GI) should feature significantly in the disk's evolution. At the same time, the disk may be sufficiently ionised for magnetic fields to play some role in the dynamics. Though usually neglected, the impact of magnetism on the GI may be critical, with consequences for several processes: the efficiency of accretion, spiral structure formation, fragmentation, and the dynamics of solids. In this paper, we report on global three-dimensional magnetohydrodynamical simulations of a self-gravitating protoplanetary disk using the meshless finite mass (MFM) Lagrangian technique. We confirm that GI spiral waves trigger a dynamo that amplifies an initial magnetic field to nearly thermal amplitudes (plasma beta < 10), an order of magnitude greater than that generated by the magneto-rotational instability alone. We also determine the dynamo's nonlinear back reaction on the gravitoturbulent flow: the saturated state is substantially hotter, with an associated larger Toomre parameter and weaker, more 'flocculent' spirals. But perhaps of greater import is the dynamo's boosting of accretion via a significant Maxwell stress; mass accretion is enhanced by factors of several relative to either pure GI or pure MRI. Our simulations use ideal MHD, an admittedly poor approximation in protoplanetary disks, and thus future studies should explore the full gamut of non-ideal MHD. In preparation for that, we exhibit a small number of Ohmic runs that reveal that the dynamo, if anything, is stronger in a non-ideal environment. This work confirms that magnetic fields are a potentially critical ingredient in gravitoturbulent young disks, possibly controlling their evolution, especially via their enhancement of (potentially episodic) accretion., Happy Chinese new year

Magnetars are highly magnetized neutron stars whose magnetic dipole ranges from $10^{14}$ to $10^{15}$ G. The MRI is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars and form magnetars. This scenario is supported by many local studies showing that magnetic fields could be amplified by the MRI on small scales. However, the efficiency of the MRI at generating a dipole field is still unknown. To answer this question, we study the MRI dynamo in an idealized global model of a fast rotating protoneutron star with differential rotation. We perform 3D incompressible MHD simulations in spherical geometry with explicit diffusivities where the differential rotation is forced at the outer boundary. We vary the initial magnetic field and investigated different magnetic boundary conditions. These simulations were compared to local shearing box simulations. We obtain a self-sustained turbulent MRI-driven dynamo, whose saturated state is independent of the initial magnetic field. The MRI generates a strong turbulent magnetic field of $B \geq 2\times 10^{15}$ G and a non-dominant magnetic dipole, which represents systematically about $5\%$ of the averaged magnetic field strength. Interestingly, this dipole is tilted towards the equatorial plane. We find that local shearing box models can reproduce fairly well several characteristics of global MRI turbulence such as the kinetic and magnetic spectra. The turbulence is nonetheless more vigorous in the local models than in the global ones. Overall, our results support the ability of the MRI to form magnetar-like large-scale magnetic fields. They furthermore predict the presence of a stronger small-scale magnetic field. The resulting magnetic field could be important to power outstanding stellar explosions, such as superluminous supernovae and GRBs., Accepted for publication in A&A, 16 pages, 21 figures. Minor changes after referee report. Addition of a field line figure

High-frequency quasi-periodic oscillations (HFQPOs) observed in the emission of black-hole X-ray binary systems promise insight into strongly curved spacetime. `Discoseismic' oscillations with frequencies set by the intrinsic properties of the central black hole, in particular `trapped inertial waves' (r-modes), offer an attractive explanation for HFQPOs. To produce an observable signature, however, such oscillations must be excited to sufficiently large amplitudes. Turbulence driven by the magnetorotational instability (MRI) fails to provide the necessary amplification, but r-modes may still be excited via interaction with accretion disc warps or eccentricities. We present 3D global hydrodynamic simulations of relativistic accretion discs, which demonstrate for the first time the excitation of trapped inertial waves by an imposed eccentricity in the flow. While the r-modes' saturated state depends on the vertical boundary conditions used in our unstratified, cylindrical framework, their excitation is unambiguous in all runs with eccentricity >0.005 near the ISCO. These simulations provide a proof of concept, demonstrating the robustness of trapped inertial wave excitation in a non-magnetized context. In a companion paper, we explore the competition between this excitation, and damping by magnetohydrodynamic turbulence., 16 pages, 15 figures, accepted in MNRAS. Updated with grammatical and typographical corrections

Aims and Methods. Accretion bursts triggered by the magnetorotational instability (MRI) in the innermost disk regions were studied for protoplanetary gas-dust disks formed from prestellar cores of various mass $M_{\rm core}$ and mass-to-magnetic flux ratio $\lambda$. Numerical magnetohydrodynamics simulations in the thin-disk limit were employed to study the long-term ($\sim 1.0$~Myr) evolution of protoplanetary disks with an adaptive turbulent $\alpha$-parameter, which depends explicitly on the strength of the magnetic field and ionization fraction in the disk. The numerical models also feature the co-evolution of gas and dust, including the back-reaction of dust on gas and dust growth. Results. Dead zone with a low ionization fraction $x, Comment: Accepted for publication in Astronomy & Astrophysics

We investigate sustenance and dependence on magnetic Prandtl number (${\rm Pm}$) for magnetorotational instability (MRI)-driven turbulence in astrophysical Keplerian disks with zero net magnetic flux using standard shearing box simulations. We focus on the turbulence dynamics in Fourier space, capturing specific/noncanonical anisotropy of nonlinear processes due to disk flow shear. This is a new type of nonlinear redistribution of modes over wavevector orientations in Fourier space -- the nonlinear transverse cascade -- which is generic to shear flows and fundamentally different from usual direct/inverse cascade. The zero flux MRI has no exponentially growing modes, so its growth is transient, or nonmodal. Turbulence self-sustenance is governed by constructive cooperation of the transient growth of MRI and the nonlinear transverse cascade. This cooperation takes place at small wavenumbers (on the flow size scales) referred to as the vital area in Fourier space. The direct cascade transfers mode energy from the vital area to larger wavenumbers. At large ${\rm Pm}$, the transverse cascade prevails over the direct one, keeping most of modes' energy contained in small wavenumbers. With decreasing ${\rm Pm}$, however, the action of the transverse cascade weakens and can no longer oppose the action of direct cascade which more efficiently transfers energy to higher wavenumbers, leading to increased resistive dissipation. This undermines the sustenance scheme, resulting in the turbulence decay. Thus, the decay of zero net flux MRI-turbulence with decreasing ${\rm Pm}$ is attributed to topological rearrangement of the nonlinear processes when the direct cascade begins to prevail over the transverse cascade., 26 pages, 16 figures, accepted for publication in the Astrophysical Journal

Previous axisymmetric investigations in two dimensions (2D) have shown that rings and gaps develop naturally in non-ideal magnetohydrodynamic (MHD) disk-wind systems, especially in the presence of ambipolar diffusion (AD). Here we extend these 2D simulations to three dimensions (3D) and find that rings and gaps are still formed in the presence of a moderately strong ambipolar diffusion. The rings and gaps form from the same basic mechanism that was identified in the 2D simulations, namely, the redistribution of the poloidal magnetic flux relative to the disk material as a result of the reconnection of a sharply pinched poloidal magnetic field lines. Thus, the less dense gaps are more strongly magnetized with a large poloidal magnetic field compared to the less magnetized (dense) rings. The rings and gaps start out rather smoothly in 3D simulations that have axisymmetric initial conditions. Non-axisymmetric variations arise spontaneously at later times, but they do not grow to such an extent as to disrupt the rings and gaps. These disk substructures persist to the end of the simulations, lasting up to 3000 orbital periods at the inner edge of the simulated disk. The longevity of the perturbed yet still coherent rings make them attractive sites for trapping large grains that would otherwise be lost to rapid radial migration due to gas drag. As the ambipolar diffusivity decreases, both the disk and the wind become increasingly turbulent, driven by the magnetorotational instability, with tightly-wound spiral arms becoming more prominent in the disk., Comment: 19 pages, 13 figures, accepted by MNRAS

We show that on-off intermittency in solar and stellar cycles is a result of amplitude-phase synchronization in multiscale interactions in solar/stellar dynamos or magnetorotational instability which leads to the formation of kinematic and magnetic coherent structures, and the novel techniques of Lagrangian coherent structures can detect transport barriers and vortices such as magnetic flux tubes/ropes in solar and stellar turbulence with high accuracy.

Geometrically thick tori with constant specific angular momentum have been widely used in the last decades to construct numerical models of accretion flows onto black holes. Such discs are prone to a global non-axisymmetric hydrodynamic instability, known as Papaloizou-Pringle instability (PPI), which can redistribute angular momentum and also lead to an emission of gravitational waves. It is, however, not clear yet how the development of the PPI is affected by the presence of a magnetic field and by the concurrent development of the magnetorotational instability (MRI). We present a numerical analysis using three-dimensional GRMHD simulations of the interplay between the PPI and the MRI considering, for the first time, an analytical magnetized equilibrium solution as initial condition. In the purely hydrodynamic case, the PPI selects as expected the large-scale $m=1$ azimuthal mode as the fastest growing and non-linearly dominant mode. However, when the torus is threaded by a weak toroidal magnetic field, the development of the MRI leads to the suppression of large-scale modes and redistributes power across smaller scales. If the system starts with a significantly excited $m=1$ mode, the PPI can be dominant in a transient phase, before being ultimately quenched by the MRI. Such dynamics may well be important in compact star mergers and tidal disruption events., Comment: 15 pages, 13 figures, submitted to MNRAS

Context. Disks around pre-main-sequence stars evolve over time by turbulent viscous spreading. The main contender to explain the strength of the turbulence is the Magneto-Rotational-Instability (MRI) model, whose efficiency depends on the disk ionization fraction. Aims. We aim at computing self-consistently the chemistry including PAH charge chemistry, the grain charging and an estimate of an effective value of the turbulence alpha parameter in order to find observational signatures of disk turbulence. Methods. We introduced PAH and grain charging physics and their interplay with other gas-phase reactions in the physico-chemical code ProDiMo. Non-ideal magnetohydrodynamics effects such as Ohmic and ambipolar diffusion are parametrized to derive an effective value for the turbulent parameter alpha_eff . We explored the effects of turbulence heating and line broadening on CO isotopologue sub-millimeter lines. Results. The spatial distribution of alpha_eff depends on various unconstrained disk parameters such as the magnetic parameter beta_mag or the cosmic ray density distribution inside the protoplanetary disks. The inner disk midplane shows the presence of the so-called dead-zone where the turbulence is quasi-inexistent. The disk is heated mostly by thermal accommodation on dust grains in the dead-zone, by viscous heating outside the dead-zone up to a few hundred astronomical units, and by chemical heating in the outer disk. The CO rotational lines probe the warm molecular disk layers where the turbulence is at its maximum. However, the effect of turbulence on the CO line profiles is minimal and difficult to distinguish from the thermal broadening. Conclusions. Viscous heating of the gas in the disk midplane outside the dead-zone is efficient. The determination of alpha from CO rotational line observations alone is challenging., accepted to A&A. 29 pages

The Event Horizon Telescope Collaboration, Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope, and high- precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD, and KORAL. Agreement among the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. © 2019. The American Astronomical Society. All rights reserved., R.N. thanks the National Science Foundation (NSF; grants OISE-1743747, AST-1816420) and acknowledges computational support from the NSF via XSEDE resources (grant TG-AST080026N). L.D.Z. acknowledges support from the PRIN-MIUR project Multi-scale Simulations of High-Energy Astrophysical Plasmas (Prot. 2015L5EE2Y) and from the INFN-TEONGRAV initiative. C.J.W. made use of thComision Nacional de Investigacion Cientifica y Tecnologica (CONICYT, ChileComision Nacional de Investigacion Cientifica y Tecnologica (CONICYT, Chilee Extreme Science and Engineering Discovery Environment (XSEDE) Comet at the San Diego Supercomputer Center through allocation AST170012. The H-AMR high-resolution simulation was made possible by NSF PRAC award Nos. 1615281 and OAC-1811605 at the Blue Waters sustained-petascale computing project and supported in part under grant No. NSF PHY-1125915 (PI A. Tchekhovskoy). K.C. and S.M. are supported by the Netherlands Organization for Scientific Research (NWO) VICI grant (No. 639.043.513); M.L. is supported by the NWO Spinoza Prize (PI M.B.M. van der Klis). The HARM-Noble simulations were made possible by NSF PRAC award No. NSF OAC-1515969, OAC-1811228 at the Blue Waters sustained-petascale computing project, and supported in part under grant No. NSF PHY-1125915. The BHAC CKS-GRMHD simulations were performed on the Dutch National Supercomputing cluster Cartesius and are funded by the NWO computing grant 16431. S.C.N. was supported by an appointment to the NASA Postdoctoral Program at the Goddard Space Flight Center administered by USRA through a contract with NASA. Y.M., H.O., O.P., and L.R. acknowledge support from the ERC synergy grant >BlackHoleCam: Imaging the Event Horizon of Black Holes> (grant No. 610058). M.B. acknowledges support from the European Research Council (grant No. 715368-MagBURST) and from the Gauss Centre for Supercomputing e.V. (www.Gauss-centre.eu) for funding this project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (www.lrz.de).P.C.F.was supported by NSF grant AST-1616185 and used resources from the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by NSF grant No. ACI-1548562. Work by P.A. was performed in part under the auspices of the US Department of Energy by Lawrence Livermore National Laboratory under contract DE-AC52-07NA27344. The authors of the present paper further thank the following organizations and programs: the Academy of Finland (projects 274477, 284495, 312496); the Advanced European Network of E-infrastructures for Astronomy with the SKA (AENEAS) project, supported by the European Commission Framework Programme Horizon 2020 Research and Innovation action under grant agreement 731016; the Alexander von Humboldt Stiftung; the Black Hole Initiative at Harvard University, through a grant (60477) from the John Templeton Foundation; the China Scholarship Council; Comision Nacional de Investigacion Cientifica y Tecnologica (CONICYT, Chile, via PIA ACT172033, Fondecyt 1171506, BASAL AFB-170002, ALMA-conicyt 31140007); Consejo Nacional de Ciencia y Tecnologia (CONACYT, Mexico, projects 104497, 275201, 279006, 281692); the Delaney Family via the Delaney Family John A. Wheeler Chair at Perimeter Institute; Direccion General de Asuntos del Personal Academico-Universidad Nacional Autonoma de Mexico (DGAPA-UNAM, project IN112417); the European Research Council Synergy Grant >BlackHoleCam: Imaging the Event Horizon of Black Holes> (grant 610058); the Generalitat Valenciana postdoctoral grant APOSTD/2018/177; the Gordon and Betty Moore Foundation (grants GBMF-3561, GBMF-5278); the Istituto Nazionale di Fisica Nucleare (INFN) sezione di Napoli, iniziative specifiche TEONGRAV; the International Max Planck Research School for Astronomy and Astrophysics at the Universities of Bonn and Cologne; the Jansky Fellowship program of the National Radio Astronomy Observatory (NRAO); the Japanese Government (Monbukagakusho: MEXT) Scholarship; the Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for JSPS Research Fellowship (JP17J08829); the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (CAS, grants QYZDJ-SSW-SLH057, QYZDJ-SSW-SYS008); the Leverhulme Trust Early Career Research Fellowship; the Max-Planck-Gesellschaft (MPG); the Max Planck Partner Group of the MPG and the CAS; the MEXT/JSPS KAKENHI (grants 18KK0090, JP18K13594, JP18K03656, JP18H03721, 18K03709, 18H01245, 25120007); the MIT International Science and Technology Initiatives (MISTI) Funds; the Ministry of Science and Technology (MOST) of Taiwan (105-2112-M-001-025-MY3, 106-2112-M-001-011, 106-2119-M-001-027, 107-2119-M-001-017, 107-2119-M-001-020, and 107-2119-M-110-005); the National Aeronautics and Space Administration (NASA, Fermi Guest Investigator grant 80NSSC17K0649); the National Institute of Natural Sciences (NINS) of Japan; the National Key Research and Development Program of China (grant 2016YFA0400704, 2016YFA0400702); the National Science Foundation (NSF, grants AST-0096454, AST-0352953, AST-0521233, AST-0705062, AST-0905844, AST-0922984, AST-1126433, AST-1140030, DGE-1144085, AST-1207704, AST-1207730, AST-1207752, MRI-1228509, OPP-1248097, AST-1310896, AST-1312651, AST-1337663, AST-1440254, AST-1555365, AST-1715061, AST-1615796, AST-1716327, OISE-1743747, AST-1816420); the Natural Science Foundation of China (grants 11573051, 11633006, 11650110427, 10625314, 11721303, 11725312); the Natural Sciences and Engineering Research Council of Canada (NSERC, including a Discovery Grant and the NSERC Alexander Graham Bell Canada Graduate Scholarships Doctoral Program); the National Youth Thousand Talents Program of China; the National Research Foundation of Korea (the Global PhD Fellowship Grant: grants NRF-2015H1A2A1033752, 2015-R1D1A1A01056807, the Korea Research Fellowship Program: NRF-2015H1D3A1066561); the Netherlands Organization for Scientific Research (NWO) VICI award (grant 639.043. 513) and Spinoza Prize SPI 78-409; the New Scientific Frontiers with Precision Radio Interferometry Fellowship awarded by the South African Radio Astronomy Observatory (SARAO), which is a facility of the National Research Foundation (NRF), an agency of the Department of Science and Technology (DST) of South Africa; the Onsala Space Observatory (OSO) national infrastructure, for the provisioning of its facilities/observational support (OSO receives funding through the Swedish Research Council under grant 2017-00648); the Perimeter Institute for Theoretical Physics (research at Perimeter Institute is supported by the Government of Canada through the Department of Innovation, Science and Economic Development, and by the Province of Ontario through the Ministry of Research, Innovation and Science); the Russian Science Foundation (grant 17-12-01029); the Spanish Ministerio de Economia y Competitividad (grants AYA2015-63939-C2-1-P, AYA2016-80889-P); the State Agency for Research of the Spanish MCIU through the >Center of Excellence Severo Ochoa> award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709); the Toray Science Foundation; the US Department of Energy (USDOE) through the Los Alamos National Laboratory (operated by Triad National Security, LLC, for the National Nuclear Security Administration of the USDOE (Contract 89233218CNA000001)); the Italian Ministero dell'Istruzione Universita e Ricerca through the grant Progetti Premiali 2012-iALMA (CUP C52I13000140001); the European Union's Horizon 2020 research and innovation programme under grant agreement No. 730562 RadioNet; ALMA North America Development Fund; the Academia Sinica; and Chandra TM6-17006X. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), supported by NSF grant ACI-1548562, and CyVerse, supported by NSF grants DBI-0735191, DBI-1265383, and DBI-1743442. XSEDE Stampede2 resource at TACC was allocated through TG-AST170024 and TG-AST080026N. XSEDE JetStream resource at PTI and TACC was allocated through AST170028. The simulations were performed in part on the SuperMUC cluster at the LRZ in Garching, on the LOEWE cluster in CSC in Frankfurt, and on the HazelHen cluster at the HLRS in Stuttgart. This research was enabled in part by support provided by Compute Ontario (http://computeontario.ca), Calcul Quebec (http://www.calculquebec.ca), and Compute Canada (http://www.computecanada.ca).

The evolution of a large-scale poloidal magnetic field in an accretion disc is an important problem because it determines the launching of winds and the feasibility of the magnetorotational instability to generate turbulence or channel flows. Recent studies, both semi-analytical calculations and numerical simulations, have highlighted the crucial role non-ideal MHD effects (Ohmic resistivity, Hall drift and ambipolar diffusion), relevant in the protoplanetary disc context, might play in magnetic flux evolution in the disc. We investigate the flux transport in discs through the use of two one-dimensional semi-analytic models in the vertical direction, exploring regimes where different physical source terms and effects dominate. The governing equations for both models are derived by performing an asymptotic expansion in the limit of a thin disc, with the different regimes isolated through setting the relative order of the leading terms between variables. Flux transport rates and vertical structure profiles are calculated for a range of diffusivities and disc magnetisations. We found that Ohmic and ambipolar diffusivities drive radially outward flux transport with an outwardly inclined field. A wind outflow drives inward flux transport, which is significantly enhanced in the presence of Hall drift in the positive polarity case, $\eta_H (\boldsymbol{B}_z \cdot \boldsymbol{\Omega}) > 0$, an effect which has only been briefly noted before. Coupled only with outward inclination, the Hall effect reduces the flux transport given by a background Ohmic and/or ambipolar diffusivity, but drives no flux transport when it is the only non-ideal effect present., Comment: In press. Accepted by MNRAS

We develop a framework for magnetohydrodynamical (MHD) simulations in a local cylindrical shearing box by extending the formulation of the Cartesian shearing box. We construct shearing-periodic conditions at the radial boundaries of a simulation box from the conservation relations of the basic MHD equations, taking into account the explicit radial dependence of physical quantities. We demonstrate quasi-steady mass accretion, which cannot be handled by the standard Cartesian shearing box model, with an ideal MHD simulation in a vertically unstratified cylindrical shearing box up to 200 rotations. In this demonstrative run we set up (i) net vertical magnetic flux, (ii) a locally isothermal equation of state, and (iii) a sub-Keplerian equilibrium rotation, whereas the sound velocity and the initial Alfven velocity have the same radial dependence as that of the Keplerian velocity. Inward mass accretion is induced to balance with the outward angular momentum flux of the MHD turbulence triggered by the magnetorotational instability in a self-consistent manner. We discuss detailed physical properties of the saturated magnetic field, in comparison to the results of a Cartesian shearing box simulation., 27 pages, 10 figures, accepted by PASJ, movie is available at http://ea.c.u-tokyo.ac.jp/astro/Members/stakeru/research/cylshbx/

The gas temperature in protoplanetary disks (PPDs) is determined by a combination of irradiation heating and accretion heating, with the latter conventionally attributed to turbulent dissipation. However, recent studies have suggested that the inner disk (a few AU) is largely laminar, with accretion primarily driven by magnetized disk winds, as a result of nonideal magnetohydrodynamic (MHD) effects from weakly ionized gas, suggesting an alternative heating mechanism by Joule dissipation. We perform local stratified MHD simulations including all three nonideal MHD effects (ohmic, Hall, and ambipolar diffusion) and investigate the role of Joule heating and the resulting disk vertical temperature profiles. We find that in the inner disk, as ohmic and ambipolar diffusion strongly suppress electrical current around the midplane, Joule heating primarily occurs at several scale heights above the midplane, making the midplane temperature much lower than that with the conventional viscous heating model. Including the Hall effect, Joule heating is enhanced/reduced when the magnetic fields threading the disks are aligned/anti-aligned with the disk rotation, but it is overall ineffective. Our results further suggest that the midplane temperature in the inner PPDs is almost entirely determined by irradiation heating, unless viscous heating can trigger thermal ionization in the disk innermost region to self-sustain magnetorotational instability turbulence., Comment: 23 pages, 11 figures; matches the published version

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

We analyze the sustaining mechanism of nonlinear perturbations/turbulence in spectrally stable smooth shear flows. The essence of the sustenance is a subtle interplay of linear transient growth of Fourier harmonics and nonlinear processes. In spectrally stable shear flows, the transient growth of perturbations is strongly anisotropic in spectral (k-)space. This, in turn, leads to anisotropy of nonlinear processes ink-space and, as a result, the main (new) nonlinear process appears to be not a direct/inverse, but rather a transverse/angular redistribution of harmonics in Fourier space referred to as the nonlinear transverse cascade. It is demonstrated that nonlinear state is sustained owing to the interplay of the linear nonmodal growth and the transverse cascade. The possibility of such course of events has been described ink-space byG. Chagelishvili, J.-P. Zahn, A. Tevzadze and J. Lominadze, A&A, 402, 401 (2003)that reliably exemplifies the well-known bypass scenario of subcritical turbulence in spectrally stable shear flows. We present selected results of the simulations performed in different (HD and MHD; 2D and 3D; plane and Keplerian) shear flows to demonstrate the transverse cascade in action.

Lunar rocks are severely depleted in moderately volatile elements such as Rb, K, and Zn relative to Earth. Identifying the cause of this depletion is important for understanding how the Earth-Moon system evolved in the aftermath of the Moon-forming giant impact. We measured the Rb isotopic compositions of lunar and terrestrial rocks to understand why moderately volatile elements are depleted in the Moon. Combining our new measurements with previous data reveals that the Moon has an 87Rb/85Rb ratio higher than Earth by +0.16 permil. This isotopic composition is consistent with evaporation of Rb into a vapor medium that was ~99 percent saturated. Evaporation under this saturation can also explain the previously documented isotopic fractionations of K, Ga, Cu and Zn of lunar rocks relative to Earth. We show that a possible setting for achieving the same saturation upon evaporation of elements with such diverse volatilities is through viscous drainage of a partially vaporized protolunar disk onto Earth. In the framework of an alpha-disk model, the alpha-viscosity needed to explain the ~99% saturation calculated here is 0.001 to 0.01, which is consistent with a vapor disk where viscosity is controlled by magnetorotational instability., Comment: 28 pages, 8 figures, 2 tables

Magnetorotational instability (MRI) in complex quantum plasma, which comprises electron, ion, and dust grains, is investigated. Electrons are considered to be degenerate while ions and dust grains are taken as nondegenerate. The general dispersion relation for MRI with the effect of spin magnetization is derived using local approximations. The instability criteria for spin quantum plasma are derived using low-frequency approximations. Due to the complex nature of the dispersion relation, MHD assumptions are used to define the instability criteria. The dust mass plays a significant role in both the dispersion and growth rate of the magnetorotational mode. The complex mass of dust modifies the instability criteria. Spin magnetization force arises due to electron Pauli spin matrices and has an important role in low-temperature and highly dense plasmas that are found mainly in white dwarfs, neutron stars, and in core-collapse supernovae. It is found that spin magnetization influences the criteria for instability and hence the properties of MHD waves in the degenerate plasma.

We solve the two-dimensional hydrodynamic equations of hot accretion flow in the presence of the thermal conduction. The flow is assumed to be in steady-state and axisymmetric, and self-similar approximation is adopted in the radial direction. In this hydrodynamic study, we consider the viscous stress tensor to mimic the effects of the magnetorotational instability for driving angular momentum. We impose the physical boundary conditions at both the rotation axis and the equatorial plane and obtain the solutions in the full $ r-\theta $ space. We have found that thermal conduction is indispensable term for investigating the inflow-wind structure of the hot accretion flows with very low mass accretion rates. One of the most interesting results here is that the disc is convectively stable in hot accretion mode and in the presence of the thermal conduction. Furthermore, the properties of wind and also its driving mechanisms are studied. Our analytical results are consistent with previous numerical simulations of hot accretion flow., Comment: 20 pages, 10 figures, accepted for publication in ApJ

We investigate evolution of an accretion disc in binary black hole (BBH) systems and possible electromagnetic counterparts of the gravitational waves from mergers of BBHs. Perna et al. (2016) proposed a novel evolutionary scenario of an accretion disc in BBHs in which a disc eventually becomes "dead", i.e., the magnetorotational instability (MRI) becomes inactive. In their scenario, the dead disc survives until {\it a few seconds before} the merger event. We improve the dead disc model and propose another scenario, taking account of effects of the tidal torque from the companion and the critical ionization degree for MRI activation more carefully. We find that the mass of the dead disc is much lower than that in the Perna's scenario. When the binary separation sufficiently becomes small, the mass inflow induced by the tidal torque reactivates MRI, restarting mass accretion onto the black hole. We also find that this disc "revival" happens {\it more than thousands of years before} the merger. The mass accretion induced by the tidal torque increases as the separation decreases, and a relativistic jet could be launched before the merger. The emissions from these jets are too faint compared to GRBs, but detectable if the merger events happen within $\lesssim 10$ Mpc or if the masses of the black holes are as massive as $\sim 10^5 M_{\odot}$., Comment: 8 pages, 3 figures, accepted for publication in MNRAS

We study the dynamo generation (exponential growth) of large scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal ($x$-$y$) averaging, we demonstrate the presence of large scale fields when either horizontal or vertical ($y$-$z$) averaging is employed. By computing planar averaged fields and power spectra, we find large scale dynamo action in the early MRI growth phase---a previously unidentified feature. Fast growing horizontal low modes and fiducial vertical modes over a narrow range of wave numbers amplify these planar averaged fields in the MRI growth phase, before turbulence sets in. The large scale field growth requires linear fluctuations but not nonlinear turbulence (as defined by mode-mode coupling) and grows as a direct global mode of the MRI. Only by vertical averaging, can it be shown that the growth of horizontal low wavenumber MRI modes directly feed-back to the initial vertical field providing a clue as to why the large scale vertical field sustains against turbulent diffusion in the saturation regime. We compute the terms in the planar averaged mean field equations to identify the individual contributions to large scale field growth for both vertical and horizontal averaging. The large scale fields obtained from such vertical averaging are found to compare well with global cylindrical simulations and quasilinear analytical analysis from a previous study by Ebrahimi \& Blackman. We discuss the potential implications of these new results for understanding large scale MRI dynamo saturation and turbulence., 12 pages, 9 figues, submitted to MNRAS, comments are welcome

The circumnuclear disc (CND) orbiting the Galaxy's central black hole is a reservoir of material that can ultimately provide energy through accretion, or form stars in the presence of the black hole, as evidenced by the stellar cluster that is presently located at the CND's centre. In this paper, we report the results of a computational study of the dynamics of the CND. The results lead us to question two paradigms that are prevalent in previous research on the Galactic Centre. The first is that the disc's inner cavity is maintained by the interaction of the central stellar cluster's strong winds with the disc's inner rim, and second, that the presence of unstable clumps in the disc implies that the CND is a transient feature. Our simulations show that, in the absence of a magnetic field, the interaction of the wind with the inner disc rim actually leads to a filling of the inner cavity within a few orbital time-scales, contrary to previous expectations. However, including the effects of magnetic fields stabilizes the inner disc rim against rapid inward migration. Furthermore, this interaction causes instabilities that continuously create clumps that are individually unstable against tidal shearing. Thus the occurrence of such unstable clumps does not necessarily mean that the disc is itself a transient phenomenon. The next steps in this investigation are to explore the effect of the magnetorotational instability on the disc evolution and to test whether the results presented here persist for longer time-scales than those considered here., 18 pages, 11 figures, accepted by MNRAS

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

Aims. The response of a protoplanetary disk to luminosity bursts of various durations is studied with the purpose to determine the effect of the bursts on the strength and sustainability of gravitational instability in the disk. A special emphasis is paid to the spatial distribution of gas and grown dust (from 1 mm to a few centimetres) during and after the burst. Methods. Numerical hydrodynamics simulations were employed to study the dynamics of gas and dust in the thin-disk limit. Dust-to-gas friction, including back reaction and dust growth, were also considered. Bursts of various durations (from 100 yr to 500 yr) were initiated in accordance with a thermally ignited magnetorotational instability. Luminosity curves for constant- and declining-magnitude bursts were adopted to represent two typical limiting cases for FU Orionis-type eruptions. Results. The short-term effect of the burst is to reduce the strength of gravitational instability by heating and expanding the disk. The longest bursts with durations comparable to the revolution period of the spiral can completely dissolve the original two-armed spiral pattern in the gas disk by the end of the burst, while the shortest bursts only weaken the spiral pattern. The reaction of grown dust to the burst is somewhat different. The spiral-like initial distribution with deep cavities in the inter-armed regions transforms into a ring-like distribution with deep gaps. This transformation is mostly expressed for the longest-duration bursts. The long-term effect of the burst depends on the initial disk conditions at the onset of the burst. In some cases, vigorous disk fragmentation sets in several thousands of years after the burst, which was absent in the model without the burst. Several clumps with masses in the giant-planet mass range form in the outer disk regions. After the disk fragmentation phase, the spatial distribution of grown dust is characterized by multiple sharp rings located from tens to hundreds of astronomical units. The arrangement and sharpness of the rings depends on the strength of dust turbulent diffusion. The wide-orbit rings are likely formed as the result of dust-rich clump dispersal in the preceding gravitationally unstable phase. Conclusions. Luminosity bursts similar in magnitude to FU Orionis-type eruptions can have a profound effect on the dynamics of gas and dust in protoplanetary disks if the burst duration is comparable to, or longer than, the dynamical timescales. In this context, the spatial morphology of the gas-dust disk of V883 Ori, a FU Orionis-like object that is thought to be in the outburst phase for more than a century with an unknown onset date, may be used as test case for the burst models considered in this work. The potential relation of the obtained ring structures to a variety of gaps and rings observed in T Tauri disks remains to be established.

Outflows driven by large-scale magnetic fields likely play an important role in the evolution and dispersal of protoplanetary disks, and in setting the conditions for planet formation. We extend our 2-D axisymmetric non-ideal MHD model of these outflows by incorporating radiative transfer and simplified thermochemistry, with the twin aims of exploring how heating influences wind launching, and illustrating how such models can be tested through observations of diagnostic spectral lines. Our model disks launch magnetocentrifugal outflows primarily through magnetic tension forces, so the mass-loss rate increases only moderately when thermochemical effects are switched on. For typical field strengths, thermochemical and irradiation heating are more important than magnetic dissipation. We furthermore find that the entrained vertical magnetic flux diffuses out of the disk on secular timescales as a result of non-ideal MHD. Through post-processing line radiative transfer, we demonstrate that spectral line intensities and moment-1 maps of atomic oxygen, the HCN molecule, and other species show potentially observable differences between a model with a magnetically driven outflow and one with a weaker, photoevaporative outflow. In particular, the line shapes and velocity asymmetries in the moment-1 maps could enable the identification of outflows emanating from the disk surface., Comment: 35 pages, 20 figures, accepted for publication in ApJ

Turbulence in the protoplanetary disks induces collisions between dust grains, and thus facilitates grain growth. We investigate the two fundamental assumptions of the turbulence in obtaining grain collisional velocities -- the kinetic energy spectrum and the turbulence autocorrelation time -- in the context of the turbulence generated by the magneto-rotational instability (MRI). We carry out numerical simulations of the MRI as well as driven turbulence, for a range of physical and numerical parameters. We find that the convergence of the turbulence $\alpha$-parameter does not necessarily imply the convergence of the energy spectrum. The MRI turbulence is largely solenoidal, for which we observe a persistent kinetic energy spectrum of $k^{-4/3}$. The same is obtained for solenoidal driven turbulence with and without magnetic field, over more than 1 dex near the dissipation scale. This power-law slope appears to be converged in terms of numerical resolution, and to be due to the bottleneck effect. The kinetic energy in the MRI turbulence peaks at the fastest growing mode of the MRI. In contrast, the magnetic energy peaks at the dissipation scale. The magnetic energy spectrum in the MRI turbulence does not show a clear power-law range, and is almost constant over approximately 1 dex near the dissipation scale. The turbulence autocorrelation time is nearly constant at large scales, limited by the shearing timescale, and shows a power-law drop close to $k^{-1}$ at small scales, with a slope steeper than that of the eddy crossing time. The deviation from the standard picture of the Kolmogorov turbulence with the injection scale at the disk scale height can potentially have a significant impact on the grain collisional velocities., Comment: Accepted by ApJ

The prevalence and consequences of convection perpendicular to the plane of accretion discs have been discussed for several decades. Recent simulations combining convection and the magnetorotational instability have given fresh impetus to the debate, as the interplay of the two processes can enhance angular momentum transport, at least in the optically thick outburst stage of dwarf novae. In this paper we seek to isolate and understand the most generic features of disc convection, and so undertake its study in purely hydrodynamical models. First, we investigate the linear phase of the instability, obtaining estimates of the growth rates both semi-analytically, using one-dimensional spectral computations, as well as analytically, using WKBJ methods. Next we perform three-dimensional, vertically stratified, shearing box simulations with the conservative, finite-volume code PLUTO, both with and without explicit diffusion coefficients. We find that hydrodynamic convection can, in general, drive outward angular momentum transport, a result that we confirm with ATHENA, an alternative finite-volume code. Moreover, we establish that the sign of the angular momentum flux is sensitive to the diffusivity of the numerical scheme. Finally, in sustained convection, whereby the system is continuously forced to an unstable state, we observe the formation of various coherent structures, including large- scale and oscillatory convective cells, zonal flows, and small vortices., Accepted for publication in MNRAS (20 pages, 16 figures, 4 tables)

Recent numerical simulations of magnetized accretion discs show that the radial-azimuthal component of the stress tensor due to the magnetorotational instability (MRI) is well represented by a power-law function of the gas pressure rather than a linear relation which has been used in most of the accretion disc studies. The exponent of this power-law function which depends on the net flux of the imposed magnetic field is reported in the range between zero and unity. However, the physical consequences of this power-law stress-pressure relation within the framework of the standard disc model have not been explored so far. In this study, the structure of an accretion disc with a power-law stress-pressure relation is studied using analytical solutions in the steady-state and time-dependent cases. The derived solutions are applicable to different accreting systems, and as an illustrative example, we explore structure of protoplanetary discs using these solutions. We show that the slopes of the radial surface density and temperature distributions become steeper with decreasing the stress exponent. However, if the disc opacity is dominated by icy grains and value of the stress exponent is less than about $0.5$, the surface density and temperature profiles become so steep that make them unreliable. We also obtain analytical solutions for the protoplanetary discs which are irradiated by the host star. Using these solutions, we find that the effect of the irradiation becomes more significant with decreasing the stress exponent., Accepted for publication in MNRAS

MHD turbulence plays a crucial role in the dust dynamics of protoplanetary discs. It affects planet formation, vertical settling and is one possible origin of the large scale axisymmetric structures, such as rings, recently imaged by ALMA and SPHERE. Among the variety of MHD processes, the magnetorotational instability (MRI) has raised particular interest since it provides a source of turbulence and potentially organizes the flow into radial structures. However, the weak ionization of discs prevents the MRI from being excited beyond 1 AU. The strong sedimentation of millimetre dust measured in T-Tauri discs is also in contradiction with predictions based on ideal MRI turbulence. In this paper, we study the effects of non-ideal MHD and winds on the dynamics and sedimentation of dust grains. We consider a weakly ionized plasma subject to ambipolar diffusion characterizing the disc outer regions (>1 AU). For that, we perform numerical MHD simulations in the stratified shearing box, using the PLUTO code. Our simulations show that the mm-cm dust is contained vertically in a very thin layer, with typical heightscale ~0.4 AU at 30 AU, compatible with recent ALMA observations. Horizontally, the grains are trapped within pressure maxima induced by ambipolar diffusion, leading to the formation of dust rings. For micrometer grains, dust and gas scaleheights are similar. In this regime, the settling cannot be explained by a simple 1D diffusion theory but results from a large scale 2D circulation induced by both MHD winds and zonal flows. Overall, our results show that non-ideal MHD effects and their related winds play a major role in shaping the radial and vertical distribution of dust in protoplanetary discs. Leading to substantial accretion efficiency, non-ideal effects also a promising avenue to reconcile the low turbulent activity measured in discs with their relatively high accretion rates., 22 pages, 21 figures, accepted in A&A

We studied dynamical balances in magnetorotational instability (MRI) turbulence with a net vertical field in the shearing box model of disks. Analyzing the turbulence dynamics in Fourier (${\bf k}$-)space, we identified three types of active modes that define turbulence characteristics. These modes have lengths similar to the box size, i.e., lie in the small wavenumber region in Fourier space labeled the vital area and are: (i) the channel mode - uniform in the disk plane with the smallest vertical wavenumber,(ii) the zonal flow mode - azimuthally and vertically uniform with the smallest radial wavenumber and (iii) the rest modes. The rest modes comprise those harmonics in the vital area whose energies reach more than $50 \%$ of the maximum spectral energy. The rest modes individually are not so significant compared to the channel and zonal flow modes, however, the combined action of their multitude is dominant over these two modes. These three mode types are governed by interplay of the linear and nonlinear processes, leading to their interdependent dynamics. The linear processes consist in disk flow nonmodality-modified classical MRI with a net vertical field. The main nonlinear process is transfer of modes over wavevector angles in Fourier space - the transverse cascade. The channel mode exhibits episodic bursts supplied by linear MRI growth, while the nonlinear processes mostly oppose this, draining the channel energy and redistributing it to the rest modes. As for the zonal flow, it does not have a linear source and is fed by nonlinear interactions of the rest modes., 28 pages, 16 figures, published in ApJ

The streaming instability and Kelvin--Helmholtz instability are considered the two major sources causing clumping of dust particles and turbulence in the dust layer of a protoplanetary disk as long as we consider the dead zone where the magneto-rotational instability does not grow. Extensive numerical simulations have been carried out in order to elucidate the condition for the development of particle clumping caused by the streaming instability. In this paper, a set of two parameters suitable for classifying the numerical results is proposed. One is the Stokes number that has been employed in previous works and the other is the dust particle column density that is nondimensionalized using the gas density in the midplane, Keplerian angular velocity, and difference between the Keplerian and gaseous orbital velocities. The magnitude of dust clumping is a measure of the behavior of the dust layer. Using three-dimensional numerical simulations of dust particles and gas based on Athena code v. 4.2, it is confirmed that the magnitude of dust clumping for two disk models are similar if the corresponding sets of values of the two parameters are identical to each other, even if the values of the metallicity (i.e., the ratio of the columns density of the dust particles to that of the gas) are different.

The large population of Earth to super-Earth sized planets found very close to their host stars has motivated consideration of $in$ $situ$ formation models. In particular, Inside-Out Planet Formation is a scenario in which planets coalesce sequentially in the disk, at the local gas pressure maximum near the inner boundary of the dead zone. The pressure maximum arises from a decline in viscosity, going from the active innermost disk (where thermal ionization of alkalis yields high viscosities via the magneto-rotational instability (MRI)) to the adjacent dead zone (where the MRI is quenched). Previous studies of the pressure maximum, based on $\alpha$-disk models, have assumed ad hoc values for the viscosity parameter $\alpha$ in the active zone, ignoring the detailed physics of the MRI. Here we explicitly couple the MRI criteria to the $\alpha$-disk equations, to find steady-state (constant accretion rate) solutions for the disk structure. We consider the effects of both Ohmic and ambipolar resistivities, and find solutions for a range of disk accretion rates ($\dot{M}$ = $10^{-10}$ - $10^{-8}$ ${\rm M}_{\odot}$/yr), stellar masses ($M_{\ast}$ = 0.1 - 1 ${\rm M}_{\odot}$), and fiducial values of the $non$-MRI $\alpha$-viscosity in the dead zone ($\alpha_{\rm {DZ}} = 10^{-5}$ - $10^{-3}$). We find that: (1) A midplane pressure maximum forms radially $outside$ the inner boundary of the dead zone; (2) Hall resistivity dominates near the midplane in the inner disk, which may explain why close-in planets do $not$ form in $\sim$50% of systems; (3) X-ray ionization can be competitive with thermal ionization in the inner disk, because of the low surface density there in steady-state; and (4) our inner disk solutions are viscously unstable to surface density perturbations., Comment: 34 pages, 28 figures, 3 appendices. Accepted by the Astrophysical Journal

No planets exist inside the orbit of Mercury and the terrestrial planets of the solar system exhibit a localized configuration. According to thermal structure calculation of protoplanetary disks, a silicate condensation line (~ 1300 K) is located around 0.1 au from the Sun except for the early phase of disk evolution, and planetesimals could have formed inside the orbit of Mercury. A recent study of disk evolution that includes magnetically driven disk winds showed that the gas disk obtains a positive surface density slope inside ~ 1 au from the central star. In a region with positive midplane pressure gradient, planetesimals undergo outward radial drift. We investigate the radial drift of planetesimals and type I migration of planetary embryos in a disk that viscously evolves with magnetically driven disk winds. We show a case in which no planets remain in the close-in region. Radial drifts of planetesimals are simulated using a recent disk evolution model that includes effects of disk winds. The late stage of planet formation is also examined by performing N-body simulations of planetary embryos. We demonstrate that in the middle stage of disk evolution, planetesimals can undergo convergent radial drift in a magnetorotational instability (MRI)-inactive disk, in which the pressure maximum is created, and accumulate in a narrow ring-like region with an inner edge at ~ 0.7 au from the Sun. We also show that planetary embryos that may grow from the narrow planetesimal ring do not exhibit significant type I migration in the late stage of disk evolution. The origin of the localized configuration of the terrestrial planets of the solar system, in particular the deficit of close-in planets, can be explained by the convergent radial drift of planetesimals in disks with a positive pressure gradient in the close-in region., 5 pages, 4 figures, accepted for publication in A&A Letters

In this paper, regarding the recent detection of significant azimuthal magnetic field in some accretion disks such as protostellar (Donati et al. in Nature 438:466, 2005), the multi-fluid model has been employed to analysis the stability of Keplerian rotational viscous dusty plasma system in a current-free helical magnetic field structure. Using the fluid-Maxwell equations, the general dispersion relation of the excited modes in the system has been obtained by applying the local approximation method in the linear perturbation theory. The typical numerical analysis of the obtained dispersion relation in the high-frequency regime shows that the presence of azimuthal magnetic field component in Keplerian flow has a considerable role in the stability conditions of the system. It also shows that the magnetic field helicity has a stabilization role against the magnetorotational instability (MRI) in the system due to contraction of the unstable wavelength region and decreasing the maximum growth rate of the instability. In this sense, the stabilization role of the viscosity term is more considerable for HMRI (instability in the presence of azimuthal magnetic field component) than the corresponding MRI (instability in the absence of azimuthal magnetic field component). Moreover, considering the discovered azimuthal magnetic field in these systems, the MRI can be arisen in the over-all range of dust grains construction values in contract with traditional MRI. This investigation can greatly contribute to better understanding the physics of some astrophysical phenomena, such as the main source of turbulence and angular momentum transport in protostellar and the other sufficiently ionized astrophysical disks, where the azimuthal magnetic field component in these systems can play a significant role.

We present an approach to turbulence closure based on mixing length theory with three-dimensional fluctuations against a two-dimensional background. This model is intended to be rapidly computable for implementation in stellar evolution software and to capture a wide range of relevant phenomena with just a single free parameter, namely the mixing length. We incorporate magnetic, rotational, baroclinic and buoyancy effects exactly within the formalism of linear growth theories with nonlinear decay. We treat differential rotation effects perturbatively in the corotating frame using a novel controlled approximation which matches the time evolution of the reference frame to arbitrary order. We then implement this model in an efficient open source code and discuss the resulting turbulent stresses and transport coefficients. We demonstrate that this model exhibits convective, baroclinic and shear instabilities as well as the magnetorotational instability (MRI). It also exhibits non-linear saturation behaviour, and we use this to extract the asymptotic scaling of various transport coefficients in physically interesting limits., Comment: 17 pages, 20 figures, published in MNRAS. Figures 9-19 and C1 and associated text have been updated relative to MNRAS version to reflect minor changes due to corrections to the numerical code. Typographic errors have been corrected in equations 8, 10 and 11 and Figure 1 caption. Figure 3 legend updated to be clearer, and colors have been coordinated between Figures 1-3

Warm Comptonization models for the soft X-ray excess in AGN do not self-consistently explain the relationship between the Comptonizing medium and the underlying accretion disc. Because of this, they cannot directly connect the fitted Comptonization temperatures and optical depths to accretion disc parameters. Since bulk velocities exceed thermal velocities in highly radiation pressure dominated discs, in these systems bulk Comptonization by turbulence may provide a physical basis in the disc itself for warm Comptonization models. We model the dependence of bulk Comptonization on fundamental accretion disc parameters, such as mass, luminosity, radius, spin, inner boundary condition, and $\alpha$. In addition to constraining warm Comptonization models, our model can help distinguish contributions from bulk Comptonization to the soft X-ray excess from those due to other physical mechanisms, such as absorption and reflection. By linking the time variability of bulk Comptonization to fluctuations in the disc vertical structure due to magnetorotational instability (MRI) turbulence, our results show that observations of the soft X-ray excess can be used to study disc turbulence in the radiation pressure dominated regime. Because our model connects bulk Comptonization to one dimensional vertical structure temperature profiles in a physically intuitive way, it will be useful for understanding this effect in future simulations run in new regimes., Comment: 35 pages, 50 figures, accepted to MNRAS