Search

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., Comment: 13 pages, 10 figures, published on ApJ

ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks., Comment: 15 pages, 11 figures, Accepted for publication in ApJL

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

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

The magneto-rotational instability (MRI) is one of the most important processes in sufficiently ionized astrophysical disks. Grid-based simulations, especially those using the local shearing box approximation, provide a powerful tool to study the ensuing nonlinear turbulence. On the other hand, while meshless methods have been widely used in both cosmology, galactic dynamics, and planet formation they have not been fully deployed on the MRI problem. We present local unstratified and vertically stratified MRI simulations with two meshless MHD schemes: a recent implementation of SPH MHD (Price2012), and a MFM MHD scheme with a constrained gradient divergence cleaning scheme, as implemented in the GIZMO code \citep{Hopkins2017}. Concerning variants of the SPH hydro force formulation we consider both the "vanilla" SPH and the PSPH variant included in GIZMO. We find, as expected, that the numerical noise inherent in these schemes affects turbulence significantly. A high order kernel, free of the pairing instability, is necessary. Both schemes can adequately simulate MRI turbulence in unstratified shearing boxes with net vertical flux. The turbulence, however, dies out in zero-net-flux unstratified boxes, probably due to excessive and numerical dissipation. In zero-net-flux vertically stratified simulations, MFM can reproduce the MRI dynamo and its characteristic butterfly diagram for several tens of orbits before ultimately decaying. In contrast, extremely strong toroidal fields, as opposed to sustained turbulence, develop in equivalent simulations using SPH MHD. This unphysical state in SPH MHD is likely caused by a combination of excessive artificial viscosity, numerical resistivity, and the relatively large residual errors in the divergence of the magnetic field remaining even after cleaning procedures are applied., Comment: This version has been accepted by ApJS

Protoplanetary disks (PPDs) are characterized by different kinds of gas dynamics and chemistry, which are coupled via ionization, heating and cooling processes, as well as advective and turbulent transport. However, directly coupling gas dynamics with time-dependent chemistry is prohibitively computationally expensive when using comprehensive chemical reaction networks. In this paper, we evaluate the utility of a species-based network reduction method in different disk environments to produce small chemical networks that reproduce the abundances of major species found in large gas-phase chemistry networks. We find that the method works very well in disk midplane and surfaces regions, where approximately 20-30 gas phase species, connected by $\sim$50-60 gas phase reactions, are sufficient to reproduce the targeted ionization fraction and chemical abundances. Most species of the reduced networks, including major carriers of oxygen, carbon and nitrogen, also have similar abundances in the reduced and complete network models. Our results may serve as an initial effort for future hydrodynamic/magnetohydrodynamic simulations of PPDs incorporating time-dependent chemistry in appropriate regions. Accurately modeling the abundances of major species at intermediate disk heights, however, will require much more extended network incorporating gas-grain chemistry and are left for future studies., Comment: accepted for publication in ApJ

To characterize the substructures induced in protoplanetary disks by the interaction between stars in multiple systems, we study the $1.25\,$mm continuum and the $^{12}$CO$(J=2-1)$ spectral line emission of the triple systems HT Lup and AS 205, at scales of $\approx 5\,$au, as part of the 'Disk Substructures at High Angular Resolution Project' (DSHARP). In the continuum emission, we find two symmetric spiral arms in the disk around AS 205 N, with pitch angle of $14^\circ$, while the southern component AS 205 S, itself a spectroscopic binary, is surrounded by a compact inner disk and a bright ring at a radius of $34\,$au. The $^{12}$CO line exhibits clear signatures of tidal interactions, with spiral arms, extended arc-like emission, and high velocity gas, possible evidence of a recent close encounter between the disks in the AS 205 system, as these features are predicted by hydrodynamic simulations of fly-by encounters. In the HT Lup system, we detect continuum emission from all three components. The primary disk, HT Lup A, also shows two-armed symmetric spiral structure with a pitch angle of $4^\circ$, while HT Lup B and C, located at $25$ and $434\,$au in projected separation from HT Lup A, are barely resolved with $\sim5$ and $\sim10\,$au in diameter, respectively. The gas kinematics for the closest pair indicates a different sense of rotation for each disk, which could be explained by either a counter rotation of the two disks in different, close to parallel, planes, or by a projection effect of these disks with a close to $90^\circ$ misalignment between them., Comment: 15 pages, 10 figures, accepted to ApJ Letters

We present an analysis of ALMA 1.25 millimeter continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of $\sim40-60$ mas ($\sim6-7$ au). All three disks feature $m=2$ spiral patterns in conjunction with annular substructures. Gas kinematics established by $^{12}$CO $J=2-1$ observations indicate that the continuum spiral arms are trailing. The arm-interarm intensity contrasts are modest, typically less than 3. The Elias 27 spiral pattern extends throughout much of the disk, and the arms intersect the gap at $R\sim69$ au. The spiral pattern in the IM Lup disk is particularly complex-it extends about halfway radially through the disk, exhibiting pitch angle variations with radius and interarm features that may be part of ring substructures or spiral arm branches. Spiral arms also extend most of the way through the WaOph 6 disk, but the source overall is much more compact than the other two disks. We discuss possible origins for the spiral structures, including gravitational instability and density waves induced by a stellar or planetary companion. Unlike the millimeter continuum counterparts of many of the disks with spiral arms detected in scattered light, these three sources do not feature high-contrast crescent-like asymmetries or large ($R>20$ au) emission cavities. This difference may point to multiple spiral formation mechanisms operating in disks., Comment: 19 pages, 9 figures, minor typos corrected

ALMA observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project (DSHARP) resolve the dust and gas emission on angular scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this article, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm dust continuum and $^{12}$CO J=2-1 rotational line emission at a spatial resolution of 4 and 10 au, respectively. The continuum observations resolve and allow us to characterize the previously discovered dust rings at radii of 67 and 100 au. They also reveal new small scale structures, such as a dark gap at 10 au, a bright ring at 15 au, a dust crescent at a radius of 55 au, and several fainter azimuthal asymmetries. The observations of the CO and dust emission inform about the vertical structure of the disk and allow us to directly constrain the dust extinction optical depth at the dust rings. Furthermore, the observed asymmetries in the dust continuum emission corroborate to the hypothesis that the complex structure of the HD 163296 disk is the result of the gravitational interaction with yet unseen planets., Comment: ApJ Letters, in press. 19 pages, 11 figures

A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this paper we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bumps, and if confirmed, put constraints on the physics of the dust trapping mechanism. We model this process analytically in 1D, assuming axisymmetry. By comparing this model to the data, we find that all rings are consistent with dust trapping. Based on a plausible model of the dust temperature we find that several rings are narrower than the pressure scale height, providing strong evidence for dust trapping. The rings have peak absorption optical depth in the range between 0.2 and 0.5. The dust masses stored in each of these rings is of the order of tens of Earth masses, though much ambiguity remains due to the uncertainty of the dust opacities. The dust rings are dense enough to potentially trigger the streaming instability, but our analysis cannot give proof of this mechanism actually operating. Our results show, however, that the combination of very low alpha_turb << 5e-4 and very large grains a_grain >> 0.1 cm can be excluded by the data for all the rings studied in this paper., Comment: Accepted for publication in ApJ Letters

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk., Comment: 33 pages, 12 figures, captions corrected

We present a detailed analysis of new ALMA observations of the disk around the T-Tauri star HD 143006, which at 46 mas (7.6 au) resolution reveal new substructures in the 1.25 mm continuum emission. The disk resolves into a series of concentric rings and gaps together with a bright arc exterior to the rings that resembles hydrodynamics simulations of a vortex, and a bridge-like feature connecting the two innermost rings. Although our $^{12}$CO observations at similar spatial resolution do not show obvious substructure, they reveal an inner disk depleted of CO emission. From the continuum emission and the CO velocity field we find that the innermost ring has a higher inclination than the outermost rings and the arc. This is evidence for either a small ($\sim8^{\circ}$) or moderate ($\sim41^{\circ}$) misalignment between the inner and outer disk, depending on the specific orientation of the near/far sides of the inner/outer disk. We compare the observed substructures in the ALMA observations with recent scattered light data from VLT/SPHERE of this object. In particular, the location of narrow shadow lanes in the SPHERE image combined with pressure scale height estimates, favor a large misalignment of about $41^{\circ}$. We discuss our findings in the context of a dust-trapping vortex, planet-carved gaps, and a misaligned inner disk due to the presence of an inclined companion to HD 143006., Comment: 15 pages, 8 figures, accepted for publication at ApJ Letters

The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass and 1) the gaseous gap depth/width, and 2) the sub/super-Keplerian motion across the gap. We then compute dust continuum intensity maps at the frequency of the DSHARP observations. We provide the relationship between the planet mass and 1) the depth/width of the gaps at millimeter intensity maps, 2) the gap edge ellipticity and asymmetry, and 3) the position of secondary gaps induced by the planet. With these relationships, we lay out the procedure to constrain the planet mass using gap properties, and study the potential planets in the DSHARP disks. We highlight the excellent agreement between observations and simulations for AS 209 and the detectability of the young Solar System analog. Finally, under the assumption that the detected gaps are induced by young planets, we characterize the young planet population in the planet mass-semimajor axis diagram. We find that the occurrence rate for $>$ 5 $M_J$ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques., Comment: 36 pages, 22 figures, ApJ Letters, in press

We introduce the Disk Substructures at High Angular Resolution Project (DSHARP), one of the initial Large Programs conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal of DSHARP is to find and characterize substructures in the spatial distributions of solid particles for a sample of 20 nearby protoplanetary disks, using very high resolution (0.035 arcsec, or 5 au FWHM) observations of their 240 GHz (1.25 mm) continuum emission. These data provide a first homogeneous look at the small-scale features in disks that are directly relevant to the planet formation process, quantifying their prevalence, morphologies, spatial scales, spacings, symmetry, and amplitudes, for targets with a variety of disk and stellar host properties. We find that these substructures are ubiquitous in this sample of large, bright disks. They are most frequently manifested as concentric, narrow emission rings and depleted gaps, although large-scale spiral patterns and small arc-shaped azimuthal asymmetries are also present in some cases. These substructures are found at a wide range of disk radii (from a few au to more than 100 au), are usually compact ($<$10 au), and show a wide range of amplitudes (brightness contrasts). Here we discuss the motivation for the project, describe the survey design and the sample properties, detail the observations and data calibration, highlight some basic results, and provide a general overview of the key conclusions that are presented in more detail in a series of accompanying articles. The DSHARP data -- including visibilities, images, calibration scripts, and more -- are released for community use at https://almascience.org/alma-data/lp/DSHARP., Comment: ApJ Letters, in press; some figure resolutions reduced, full-res version at https://www.cfa.harvard.edu/~sandrews/DSHARP_I.pdf ; Data Release at https://almascience.org/alma-data/lp/DSHARP

We present a detailed analysis of the high-angular resolution (0.''037, corresponding to 5 au) observations of the 1.25 mm continuum and $^{12}$CO $2-1$ emission from the disk around the T Tauri star AS 209. AS 209 hosts one of the most unusual disks from the DSHARP sample, the first high angular resolution ALMA survey of disks (Andrews et al. 2018), as nearly all of the emission can be explained with concentric Gaussian rings. In particular, the dust emission consists of a series of narrow and closely spaced rings in the inner $\sim$60 au, two well-separated bright rings in the outer disk, centered at 74 and 120 au, and at least two fainter emission features at 90 and 130 au. We model the visibilities with a parametric representation of the radial surface brightness profile, consisting of a central core and 7 concentric Gaussian rings. Recent hydro-dynamical simulations of low viscosity disks show that super-Earth planets can produce the multiple gaps seen in AS 209 millimeter continuum emission. The $^{12}$CO line emission is centrally peaked and extends out to $\sim$300 au, much farther than the millimeter dust emission. We find axisymmetric, localized deficits of CO emission around four distinct radii, near 45, 75, 120 and 210 au. The outermost gap is located well beyond the edge of the millimeter dust emission, and therefore cannot be due to dust opacity and must be caused by a genuine CO surface density reduction, due either to chemical effects or depletion of the overall gas content., Comment: ApJ Letters, in press

The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very non-trivial task. This is in part due to the uncertainty in the dust opacities, an uncertainty which is not likely to be resolved any time soon. It is also partly due to the fact that, as the DSHARP survey has shown, these disk often contain regions of intermediate to high optical depth, even at millimeter wavelengths and at relatively large radius in the disk. This makes the interpretation challenging, in particular if the grains are large and have a large albedo. On the other hand, the highly structured features seen in the DSHARP survey, of which strong indications were already seen in earlier observations, provide a unique opportunity to study the dust growth and dynamics. To provide continuity within the DSHARP project, its follow-up projects, and projects by other teams interested in these data, we present here the methods and opacity choices used within the DSHARP collaboration to link the measured intensity $I_\nu$ to dust surface density $\Sigma_d$., Comment: ApJ Letters, in press

Magnetized winds may be important in dispersing protoplanetary disks and influencing planet formation. We carry out global full magnetohydrodynamic simulations in axisymmetry, coupled with ray-tracing radiative transfer, consistent thermochemistry, and non-ideal MHD diffusivities. Magnetized models lacking EUV photons ($h\nu>13.6\ \mathrm{eV}$) feature warm molecular outflows that have typical poloidal speeds $\gtrsim 4\ \mathrm{km\ s}^{-1}$. When the magnetization is sufficient to drive accretion rates $\sim 10^{-8}\ M_\odot\ \mathrm{yr}^{-1}$, the wind mass-loss rate is comparable. Such outflows are driven not centrifugally but by the pressure of toroidal magnetic fields produced by bending the poloidal field. Both the accretion and outflow rates increase with the poloidal field energy density, the former almost linearly. The mass-loss rate is also strongly affected by ionization due to UV and X-ray radiation near the wind base. Adding EUV irradiation to the system heats, ionizes, and accelerates the part of the outflow nearest the symmetry axis, but reduces the overall mass-loss rate by exerting pressure on the wind base. Most of our models are non-turbulent, but some with reduced dust abundance and therefore higher ionization fractions exhibit magnetorotational instabilities near the base of the wind., Comment: 25 pages, 16 figures; submitted to ApJ

The mechanism behind angular momentum transport in protoplanetary disks, and whether this transport is turbulent in nature, is a fundamental issue in planet formation studies. Recent ALMA observations have suggested that turbulent velocities in the outer regions of these disks are less than ~5-10% of the sound speed, contradicting theoretical predictions of turbulence driven by the magnetorotational instability (MRI). These observations have generally been interpreted to be consistent with a large-scale laminar magnetic wind driving accretion. Here, we carry out local, shearing box simulations with varying ionization levels and background magnetic field strengths in order to determine which parameters produce results consistent with observations. We find that even when the background magnetic field launches a strong largely laminar wind, significant turbulence persists and is driven by localized regions of vertical magnetic field (the result of zonal flows) that are unstable to the MRI. The only conditions for which we find turbulent velocities below the observational limits are weak background magnetic fields and ionization levels well below that usually assumed in theoretical studies. We interpret these findings within the context of a preliminary model in which a large scale magnetic field, confined to the inner disk, hinders ionizing sources from reaching large radial distances, e.g., through a sufficiently dense wind. Thus, in addition to such a wind, this model predicts that for disks with weakly turbulent outer regions, the outer disk will have significantly reduced ionization levels compared to standard models and will harbor only a weak vertical magnetic field., Comment: Accepted to ApJ, 13 pages, 6 figures, modified from original version to include more extensive discussion and analysis and a title change

It has been realized in recent years that the accretion of pebble-sized dust particles onto planetary cores is an important mode of core growth, which enables the formation of giant planets at large distances and assists planet formation in general. The pebble accretion theory is built upon the orbit theory of dust particles in a laminar protoplanetary disk (PPD). For sufficiently large core mass (in the "Hill regime"), essentially all particles of appropriate sizes entering the Hill sphere can be captured. However, the outer regions of PPDs are expected to be weakly turbulent due to the magnetorotational instability (MRI), where turbulent stirring of particle orbits may affect the efficiency of pebble accretion. We conduct shearing-box simulations of pebble accretion with different levels of MRI turbulence (strongly turbulent assuming ideal magnetohydrodynamics, weakly turbulent in the presence of ambipolar diffusion, and laminar) and different core masses to test the efficiency of pebble accretion at a microphysical level. We find that accretion remains efficient for marginally coupled particles (dimensionless stopping time tau_s ~ 0.1 - 1) even in the presence of strong MRI turbulence. Though more dust particles are brought toward the core by the turbulence, this effect is largely canceled by a reduction in accretion probability. As a result, the overall effect of turbulence on the accretion rate is mainly reflected in the changes in the thickness of the dust layer. On the other hand, we find that the efficiency of pebble accretion for strongly coupled particles (down to tau_s ~ 0.01) can be modestly reduced by strong turbulence for low-mass cores., Comment: 15 pages, 7 figures, accepted for publication in ApJ

The gas dynamics of weakly ionized protoplanetary disks (PPDs) is largely governed by the coupling between gas and magnetic fields, described by three non-ideal magnetohydrodynamical (MHD) effects (Ohmic, Hall, ambipolar). Previous local simulations incorporating these processes have revealed that the inner regions of PPDs are largely laminar accompanied by wind-driven accretion. We conduct 2D axisymmetric, fully global MHD simulations of these regions ($\sim1-20$ AU), taking into account all non-ideal MHD effects, with tabulated diffusion coefficients and approximate treatment of external ionization and heating. With net vertical field aligned with disk rotation, the Hall-shear instability strongly amplifies horizontal magnetic field, making the overall dynamics dependent on initial field configuration. Following disk formation, the disk likely relaxes into an inner zone characterized by asymmetric field configuration across the midplane that smoothly transitions to a more symmetric outer zone. Angular momentum transport is driven by both MHD winds and laminar Maxwell stress, with both accretion and decretion flows present at different heights, and modestly asymmetric winds from the two disk sides. With anti-aligned field polarity, weakly magnetized disks settle into an asymmetric field configuration with supersonic accretion flow concentrated at one side of disk surface, and highly asymmetric winds between the two disk sides. In all cases, the wind is magneto-thermal in nature characterized by mass loss rate exceeding the accretion rate. More strongly magnetized disks give more symmetric field configuration and flow structures. Deeper far-UV penetration leads to stronger and less stable outflows. Implications for observations and planet formation are also discussed., Comment: 32 pages, 22 figures, accepted for publication in ApJ

Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combined to map the turbulence throughout the vertical extent of a protoplanetary disk. Here we present ALMA observations of three molecules (DCO$^+$(3-2), C$^{18}$O(2-1) and CO(2-1)) from the disk around HD 163296. We are able to place stringent upper limits ($v_{\rm turb}<$0.06c$_s$, $<$0.05c$_s$ and $<$0.04c$_s$ for CO(2-1), C$^{18}$O(2-1) and DCO$^+$(3-2) respectively), corresponding to $\alpha\lesssim$3$\times$10$^{-3}$, similar to our prior limit derived from CO(3-2). This indicates that there is little turbulence throughout the vertical extent of the disk, contrary to theoretical predictions based on the magneto-rotational instability and gravito-turbulence. In modeling the DCO$^+$ emission we also find that it is confined to three concentric rings at 65.7$\pm$0.9 au, 149.9$^{+0.5}_{-0.7}$ au and 259$\pm$1 au, indicative of a complex chemical environment., Comment: Accepted to ApJ

The global evolution of protoplanetary disks (PPDs) has recently been shown to be largely controlled by the amount of poloidal magnetic flux threading the disk, which is further controlled by the poorly understood process of magnetic flux transport. In weakly ionized gas as in PPDs, magnetic flux is largely frozen in the electron fluid, except when resistivity is large. When the disk is largely laminar, we show that the relative drift between the electrons and ions (the Hall-drift), and the ions and neutral fluids (ambipolar-drift) can play a dominant role on the transport of magnetic flux. Using two-dimensional simulations that incorporate the Hall effect and ambipolar diffusion (AD) with prescribed diffusivities, we show that when large-scale poloidal field is aligned with disk rotation, the Hall effect rapidly drags magnetic flux inward at the midplane region, while it slowly pushes flux outward above/below the midplane. This leads to a highly radially elongated field configuration as a global manifestation of the Hall-shear instability. This field configuration further promotes rapid outward flux transport by AD at the midplane, leading to instability saturation. In quasi-steady state, magnetic flux is transported outward at approximately the same rate at all heights, and the rate is comparable to the Hall-free case. For anti-aligned field polarity, the Hall effect consistently transports magnetic flux outward, leading to a largely vertical field configuration in the midplane region. The field lines in the upper layer first bend radially inward and then outward to launch a disk wind. Overall, the net rate of outward flux transport is about twice faster than the aligned case. In addition, the rate of flux transport increases with increasing disk magnetization. The absolute rate of transport is sensitive to disk microphysics which remains to be explored in future studies., Comment: 19 pages, 10 figures, accepted for publication in ApJ