Disk Midplane Research Articles

UV processing of icy pebbles in the outer parts of VSI-turbulent disks

Icy dust particles emerge in star-forming clouds and are subsequently incorporated in protoplanetary disks, where they coagulate into larger pebbles up to millimeter in size. In the disk midplane, ices are shielded from UV radiation, but moderate levels of disk turbulence can lift small particles to the disk surface, where they can be altered, or destroyed. Nevertheless, studies of comets and meteorites generally find that ices at least partly retain their interstellar medium (ISM) composition before being accreted onto these minor bodies. We modeled this process through hydrodynamical simulations with vertical shear instability (VSI) driven turbulence in the outer protoplanetary disk. We used the PLUTO code in a 2.5 D global accretion setup and included Lagrangian dust particles of 0.1 and 1 mm sizes. In a post-processing step, we used the RADMC3D code to generate the local UV radiation field to assess the level of ice processing of pebbles. We find that a small fraction (∼17%) of 100 µm size particles are frequently lifted up to Z/R = 0.2, which can result in the loss of their pristine composition as their residence time in this layer allow effective CO and water photodissociation. The larger 1 mm size particles remain UV-shielded in the disk midplane throughout the dynamical evolution of the disk. Our results indicate that the assembly of icy bodies via the accretion of drifting millimeter-sized icy pebbles can explain the presence of pristine ice from the ISM, even in VSI-turbulent disks. Nevertheless, particles ≤100 µm experience efficient UV processing and may mix with unaltered icy pebbles, resulting in a less ISM-like composition in the midplane.

3D Gap opening in non-ideal MHD protoplanetary disks: Asymmetric accretion, meridional vortices, and observational signatures

Abstract Recent high-angular resolution ALMA observations have revealed rich information about protoplanetary disks, including ubiquitous substructures and three-dimensional gas kinematics at different emission layers. One interpretation of these observations is embedded planets. Previous 3-D planet-disk interaction studies are either based on viscous simulations, or non-ideal magnetohydrodynamics (MHD) simulations with simple prescribed magnetic diffusivities. This study investigates the dynamics of gap formation in 3-D non-ideal MHD disks using non-ideal MHD coefficients from the look-up table that is self-consistently calculated based on the thermo-chemical code. We find a concentration of the poloidal magnetic flux in the planet-opened gap (in agreement with previous work) and enhanced field-matter coupling due to gas depletion, which together enable efficient magnetic braking of the gap material, driving a fast accretion layer significantly displaced from the disk midplane. The fast accretion helps deplete the gap further and is expected to negatively impact the planet growth. It also affects the corotation torque by shrinking the region of horseshoe orbits on the trailing side of the planet. Together with the magnetically driven disk wind, the fast accretion layer generates a large, persistent meridional vortex in the gap, which breaks the mirror symmetry of gas kinematics between the top and bottom disk surfaces. Finally, by studying the kinematics at the emission surfaces, we discuss the implications of planets in realistic non-ideal MHD disks on kinematics observations.

Dust ring and gap formation by gas flow induced by low-mass planets embedded in protoplanetary disks

The observed dust rings and gaps in protoplanetary disks could be imprints of forming planets. Even low-mass planets in the 1-10 Earth-mass regime, which have not yet carved deep gas gaps, can generate such dust rings and gaps by driving a radially-outward gas flow, as shown in previous work. However, understanding the creation and evolution of these dust structures is challenging due to dust drift and diffusion, requiring an approach beyond previous steady state models. Here we investigate the time evolution of the dust surface density influenced by the planet-induced gas flow, based on post-processing three-dimensional hydrodynamical simulations. We find that planets larger than a dimensionless thermal mass of m = 0.05, corresponding to 0.3 Earth mass at 1 au or 1.7 Earth masses at 10 au, generate dust rings and gaps, provided that solids have small Stokes numbers (St ≲ 10−2) and that the disk midplane is weakly turbulent (αdiff ≲10−4). As dust particles pile up outside the orbit of the planet, the interior gap expands with time when the advective flux dominates over diffusion. Dust gap depths range from a factor of a few to several orders of magnitude, depending on planet mass and the level of midplane particle diffusion. We constructed a semi-analytic model describing the width of the dust ring and gap, and then compared it with the observational data. We find that up to 65% of the observed wide-orbit gaps could be explained as resulting from the presence of a low-mass planet, assuming αdiff = 10−5 and St = 10−3. However, it is more challenging to explain the observed wide rings, which in our model would require the presence of a population of small particles (St = 10−4). Further work is needed to explore the role of pebble fragmentation, planet migration, and the effect of multiple planets.

Evidence for a Sharp CO Snow Line Transition in a Protoplanetary Disk and Implications for Millimeter-wave Observations of CO Isotopologues

Observations of CO isotopologue emission from protoplanetary disks at millimeter wavelengths are a powerful tool for probing the CO snow line, an important marker for disk chemistry, and also for estimating total disk gas mass, a key quantity for planet formation. We use simple models to demonstrate that the vertical thickness of an isothermal layer around the disk midplane has important effects on the CO column density radial profile, with a thick layer producing a sharp CO snow line transition. We simulate next-generation Very Large Array (ngVLA) and Atacama Large Millimeter/submillimeter Array (ALMA) images to show that this sharp change in the CO column density can be detected in the derivative of the radial profile of emission from optically thin CO isotopologue lines. We apply this method to archival ALMA observations of the disk around the Herbig Ae star HD 163296 in the C17O and C18O J = 1–0 and J = 2–1 lines to identify a sharp CO snow line transition near ∼80 au (0.″8 at 101 pc), and show the CO column density decreases by more than a factor of 20. This finding is consistent with previous inferences from the steep rise of N2H+ emission, which marks the location where CO depletes. We also demonstrate that the disk’s thermal structure introduces significant systematic uncertainty to estimates of total disk gas mass derived from these lines. The substantial improvement in sensitivity envisioned for the ngVLA over ALMA for observations of ground-state lines of CO isotopologues has the potential to extend this approach to a much larger population of disks.

ALMA Observations of Proper Motions of the Dust Clumps in the Protoplanetary Disk MWC 758

To study the dust dynamics in the dust-trapping vortices in the protoplanetary disk around MWC 758, we analyzed the 1.3 mm continuum images of the MWC 758 disk obtained with the Atacama Large Millimeter/submillimeter Array in 2017 and 2021. We detect proper motions of 22 mas and 24 mas in the two dust clumps at radii of 0.″32 and 0.″54 in the disk on the plane of the sky, respectively. On the assumption that the dust clumps are located in the disk midplane, the velocities of the observed proper motions along the azimuthal direction of the inner and outer dust clumps are sub- and super-Keplerian, respectively, and both have angular velocities corresponding to the Keplerian angular velocity at a radius of 0.″46 ± 0.″04. This deviation from the Keplerian motion is not expected in the conventional theory of vortices formed by the Rossby wave instability. The observed non-Keplerian proper motions of the dust clumps are unlikely due to the disk warp and eccentricity, nor are they associated with any predicted planets. The two dust clumps are likely spatially coincident with the infrared spirals. In addition, we detect the changes in the intensity profiles of the dust clumps over the 4 yr span. Therefore, we suggest that the observed proper motions are possibly due to changes in the density distributions in the dust clumps caused by their interaction with the spirals in the disk.

Effect of Time-varying X-Ray Emission from Stellar Flares on the Ionization of Protoplanetary Disks

X-rays have significant impacts on cold, weakly ionized protoplanetary disks by increasing the ionization rate and driving chemical reactions. Stellar flares are explosions that emit intense X-rays and are the unique source of hard X-rays with an energy of ≳10 keV in the protoplanetary disk systems. Hard X-rays should be carefully taken into account in models as they can reach the disk midplane as a result of scattering in the disk atmospheres. However, previous models are insufficient to predict the hard X-ray spectra because of simplifications in flare models. We develop a model of X-ray spectra of stellar flares based on observations and flare theories. The flare temperature and nonthermal electron emissions are modeled as functions of flare energy, which allows us to better predict the hard X-ray photon flux than before. Using our X-ray model, we conduct radiative transfer calculations to investigate the impact of flare hard X-rays on disk ionization, with a particular focus on the protoplanetary disk around a T Tauri star. We demonstrate that for a flare with an energy of 1035 erg, X-ray photons with ≳5 keV increase the ionization rates more than galactic cosmic rays down to z ≈ 0.1R. The contribution of flare X-rays to the ionization at the midplane depends on the disk parameters such as disk mass and dust settling. We also find that the 10 yr averaged X-rays from multiple flares could certainly contribute to the ionization. These results emphasize the importance of stellar flares on the disk evolution.

A Search for Collisions and Planet–Disk Interactions in the Beta Pictoris Disk with 26 Years of High-precision HST/STIS Imaging

β Pictoris's well-studied debris disk and two known giant planets, in combination with the stability of the Hubble Space Telescope’s Space Telescope Imaging Spectrograph (HST/STIS) (and now also JWST), offers a unique opportunity to test planet–disk interaction models and to observe recent planetesimal collisions. We present HST/STIS coronagraphic imaging from two new epochs of data taken between 2021 and 2023, complementing earlier data taken in 1997 and 2012. This data set enables a temporal comparison with the longest baseline and highest precision of any debris disk to date, with sensitivity to variations in temporal surface brightness of sub-percent levels in the midplane of the disk. While no localized changes in surface brightness are detected, which would be indicative of a recent planetesimal collision, there is a tentative brightening of the southeast side of the disk over the past decade. We link the constraints on surface brightness variations to dynamical models of the planetary system’s evolution and to the collisional history of planetesimals. Using a coupled collisional model and injection/recovery framework, we estimate sensitivity to expanding collisional debris down to a Ceres mass per progenitor in the most sensitive regions of the disk midplane. These results demonstrate the capabilities of long-baseline, temporal studies with HST (and also soon with JWST) for constraining the physical processes occurring within debris disks.

Vertical Shear Instability with Partially Reflecting Boundary Conditions

Abstract The vertical shear instability (VSI) is widely believed to be effective in driving turbulence in protoplanetary disks. Prior studies on VSI exclusively exploit the reflecting boundary conditions (BCs) at the disk surfaces. VSI depends critically on the boundary behaviors of waves at the disk surfaces. We extend earlier studies by performing a comprehensive numerical analysis of VSI with partially reflecting BCs for both the axisymmetric and non-axisymmetric unstable VSI modes. We find that the growth rates of the unstable modes diminish when the outgoing component of the flow is greater than the incoming one for high-order body modes. When the outgoing wave component dominates, the growth of VSI is notably suppressed. We find that the non-axisymmetric modes are unstable and they grow at a rate that decreases with the azimuthal wavenumber. The different BCs at the lower and upper disk surfaces naturally lead to non-symmetric modes relative to the disk midplane. The potential implications of our studies for further understanding planetary formation and evolution in protoplanetary disks (PPDs) are also briefly discussed.

Nonideal Magnetohydrodynamic Instabilities in Protoplanetary Disks: Vertical Modes and Reflection Asymmetry

Magnetized disk winds and wind-driven accretion are an essential and intensively studied dispersion mechanism of protoplanetary disks (PPDs). However, the stability of these mechanisms has yet to be adequately examined. This paper employs semi-analytic linear perturbation theories under nonideal magnetohydrodynamics, focusing on disk models whose magnetic diffusivities vary by a few orders of magnitude from the disk midplane to its surface. Linear modes are distinguished by their symmetry with respect to the midplane. These modes have qualitatively different growth rates: symmetric modes almost always decay, while at least one antisymmetric mode always has a positive growth rate. This growth rate decreases faster than the Keplerian angular velocity with cylindrical radius R in the disk and scales steeper than R −5/2 in the fiducial disk model. The growth of antisymmetric modes breaks the reflection symmetry across the disk equatorial plane, and may occur even in the absence of the Hall effect. In the disk regions where fully developed antisymmetric modes occur, accretion flows appear only on one side of the disk, while disk winds occur only on the other. This may explain the asymmetry of some observed PPD outflows.

Magnetically Driven Turbulence in the Inner Regions of Protoplanetary Disks

Given the important role turbulence plays in the settling and growth of dust grains in protoplanetary disks, it is crucial that we determine whether these disks are turbulent and to what extent. Protoplanetary disks are weakly ionized near the midplane, which has led to a paradigm in which largely laminar magnetic field structures prevail deeper in the disk, with angular momentum being transported via magnetically launched winds. Yet, there has been little exploration of the precise behavior of the gas within the bulk of the disk. We carry out 3D, local shearing box simulations that include all three low-ionization effects (ohmic diffusion, ambipolar diffusion, and the Hall effect) to probe the nature of magnetically driven gas dynamics 1–30 au from the central star. We find that gas turbulence can persist with a generous yet physically motivated ionization prescription (order unity Elsässer numbers). The gas velocity fluctuations range from 0.03 to 0.09 of the sound speed c s at the disk midplane to ∼c s near the disk surface, and are dependent on the initial magnetic field strength. However, the turbulent velocities do not appear to be strongly dependent on the field polarity, and thus appear to be insensitive to the Hall effect. The midplane turbulence has the potential to drive dust grains to collision velocities exceeding their fragmentation limit, and likely reduces the efficacy of particle clumping in the midplane, though it remains to be seen if this level of turbulence persists in disks with lower ionization levels.

JWST-TST High Contrast: JWST/NIRCam Observations of the Young Giant Planet β Pic b

We present the first JWST/NIRCam observations of the directly imaged gas giant exoplanet β Pic b. Observations in six filters using NIRCam's round coronagraphic masks provide a high-signal-to-noise-ratio detection of β Pic b and the archetypal debris disk around β Pic over a wavelength range of ∼1.7–5 μm. This paper focuses on the detection of β Pic b and other potential point sources in the NIRCam data, following a paper by Rebollido et al. that presented the NIRCam and MIRI view of the debris disk around β Pic. We develop and validate approaches to obtaining accurate photometry of planets in the presence of bright, complex circumstellar backgrounds. By simultaneously fitting the planet’s point-spread function and a geometric model for the disk, we obtain planet photometry that is in good agreement with previous measurements from the ground. The NIRCam data support the cloudy nature of β Pic b’s atmosphere and the discrepancy between its mass as inferred from evolutionary models and the dynamical mass reported in the literature. We further identify five additional localized sources in the data, but all of them are found to be background stars or galaxies based on their color or spatial extent. We can rule out additional planets in the disk midplane above 1 M Jup outward of 2″ (∼40 au) and away from the disk midplane above 0.05 M Jup outward of 4″ (∼80 au). The inner giant planet β Pic c remains undetected behind the coronagraphic masks of NIRCam in our observations.

Uncovering an Excess of X-Ray Point Sources in the Halos of Virgo Late-type Galaxies

We present a systematic search for extraplanar X-ray point sources around 19 late-type, highly inclined disk galaxies residing in the Virgo cluster, based on archival Chandra observations reaching a source detection sensitivity of L(0.5–8 keV) ∼ 1038 erg s−1. Based on the cumulative source surface density distribution as a function of projected vertical distance from the disk midplane, we identify a statistically significant (∼3.3σ) excess of ∼20 X-ray sources within a projected vertical off-disk distance of 0.′92–2.′5 (∼4.4–12 kpc), the presence of which cannot be explained by the bulk stellar content of the individual galaxies, nor by the cosmic X-ray background. On the other hand, there is no significant evidence for an excess of extraplanar X-ray sources in a comparison sample of field late-type edge-on galaxies, for which Chandra observations reaching a similar source detection sensitivity are available. We discuss the possible origins for the observed excess, including low-mass X-ray binaries (LMXBs) associated with globular clusters, supernova-kicked LMXBs, and high-mass X-ray binaries born in recent star formation induced by ram pressure stripping of the disk gas, as well as a class of intracluster X-ray sources previously identified around early-type member galaxies of Virgo. We find that none of these X-ray populations can naturally dominate the observed extraplanar excess, although supernova-kicked LMXBs and the effect of ram pressure are most likely relevant.

Asymmetric Drift Map of the Milky Way Disk Populations between 8 -16 kpc with LAMOST and Gaia datasets

The application of asymmetric drift (AD) tomography across different populations provides valuable insights into the kinematics, dynamics, and rotation curves of the Galactic disk. By leveraging common stars identified in both the LAMOST and Gaia surveys, alongside Gaia DR3’s circular velocity curve, we conducted a qualitative exploration of asymmetric drift distributions within the Galactic disk spanning distances from 8 to 16 kpc. In the R-Z plane, we observed that the asymmetric drift is minimal near the mid-plane of the Galactic disk and gradually increases with vertical distance, resulting in a distinctive “horn” shape. Additionally, our analysis revealed that populations with higher [ α/Fe] ratios exhibit greater asymmetric drift compared to those with lower [ α/Fe] ratios. Specifically, we found the asymmetric drift around the solar location to be approximately 6 km s −1, with a median value of 16 km s −1 across the entire sample. Notably, the median asymmetric drift in the northern region of the Galactic disk (20 km s −1) surpasses that in the southern region (13 km s −1), with errors remaining within 2 km s −1. Furthermore, our investigation into mono-age stellar populations unveiled that older stellar populations tend to exhibit larger asymmetric drift and velocity dispersion, aligning closely with predictions from previous numerical models. Finally, based on chemical compositions, we observed that the median asymmetric drift of the thick disk significantly exceeds that of the thin disk and found that star formation within the thick disk primarily occurred 8-10 billion years ago, whereas the thin disk’s predominant star formation period spanned 6-8 billion years ago.

Protoplanetary Disk Polarization at Multiple Wavelengths: Are Dust Populations Diverse?

Millimeter and submillimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 μm and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can exhibit different correlations between observing wavelength and polarization morphology and fraction. We explore possible origins for these differences in polarization, including differences in dust populations and protostar properties. For RY Tau and MWC 480, which are consistent with scattering at both wavelengths, we present models of the scattering polarization from several dust grain size distributions. These models aim to reproduce two features of the observational results for these disks: (1) both disks have an observable degree of polarization at both wavelengths; and (2) the polarization fraction is higher at 3 mm than at 870 μm in the centers of the disks. For both disks, these features can be reproduced by a power-law distribution of spherical dust grains with a maximum radius of 200 μm and high optical depth. In MWC 480, we can also reproduce features (1) and (2) with a model containing large grains (a max = 490 μm) near the disk midplane and small grains (a max = 140 μm) above and below the midplane.

The First Spatially Resolved Detection of 13CN in a Protoplanetary Disk and Evidence for Complex Carbon Isotope Fractionation

Recent measurements of carbon isotope ratios in both protoplanetary disks and exoplanet atmospheres have suggested a possible transfer of significant carbon isotope fractionation from disks to planets. For a clearer understanding of the isotopic link between disks and planets, it is important to measure the carbon isotope ratios in various species. In this paper, we present a detection of the 13CN N = 2 − 1 hyperfine lines in the TW Hya disk with the Atacama Large Millimeter/submillimeter Array. This is the first spatially resolved detection of 13CN in disks, which enables us to measure the spatially resolved 12CN/13CN ratio for the first time. We conducted nonlocal thermal equilibrium modeling of the 13CN lines in conjunction with previously observed 12CN lines to derive the kinetic temperature, H2 volume density, and column densities of 12CN and 13CN. The H2 volume density is found to range between (4 − 10) × 107 cm−3, suggesting that CN molecules mainly reside in the disk's upper layer. The 12CN/13CN ratio is measured to be 70−6+9 at 30 < r < 80 au from the central star, which is similar to the 12C/13C ratio in the interstellar medium. However, this value differs from the previously reported values found for other carbon-bearing molecules (CO and HCN) in the TW Hya disk. This could be self-consistently explained by different emission layer heights for different molecules combined with preferential sequestration of 12C into the solid phase toward the disk midplane. This study reveals the complexity of the carbon isotope fractionation operating in disks.

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of h I radio ∼ h PI radio = 0.23 ± 0.03 kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane (h I FIR ∼ h PI FIR = 0.11 ± 0.01 kpc). Radio emission traces magnetic fields at scales of ≳300 pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution (<300 pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism.

Dust-gas coupling in turbulence- and MHD wind-driven protoplanetary disks: Implications for rocky planet formation

The degree of coupling between dust particles and their surrounding gas in protoplanetary disks is quantified by the dimensionless Stokes number. The Stokes number (St) governs particle size and spatial distributions, in turn establishing the dominant mode of planetary accretion in different disk regions. In this paper, we model the characteristic St of particles across time in disks evolving under both turbulent viscosity and magnetohydrodynamic (MHD) disk winds. In both turbulence- and wind-dominated disks, we find that collisional fragmentation is the limiting mechanism of particle growth. The water-ice sublimation line constitutes a critical transition point in dust settling, drift, and size regimes. For a fiducial disk evolution parameter α̃≃10−3, silicate particles inteior to the ice-line are characterized by low St (≲10−2) and sizes in the sub-mm- to 1cm-scale. Icy particles/boulders beyond the ice-line are characterized by high St (≳10−2) and sizes in the cm to dm size range. Hence, icy particles settle into a thin layer at the outer disk midplane and drift inward at velocities exceeding the gaseous accretionary flow due to substantial headwind drag. Silicate particles in the inner disk remain relatively well dispersed and are to a large extent advected inward with their surrounding gas.The St dichotomy across the ice-line translates to distinct planet formation pathways between the inner and outer disk. While pebble accretion proceeds slowly for rocky embryos within the ice-line (across most of parameter space), it does so rapidly for volatile-rich embryos beyond it, allowing for the growth of giant planet cores before disk dissipation. Through simulations of rocky planet growth, we evaluate the competition between pebble accretion and classical pairwise collisions between planetesimals. We conclude that the dominance of pebble accretion can only be realized in disks that are driven by MHD winds, slow-evolving (α̃≲10−3.5), and devoid of pressure maxima that may concentrate solids and give rise of planetesimal rings in which classical growth is enhanced. Such disks are extremely quiescent, with Shakura-Sunyaev turbulence parameters αν≲10−4. We conclude that for most of parameter space corresponding to values of αν reflected in observations of protoplanetary disks (≳10−4), pairwise planetesimal collisions constitute the dominant pathway of rocky planet accretion. Our results are discussed in the context of super-Earth origins, and lend support to the emerging view that they formed in planetesimal rings. Moreover, these results argue against a significant contribution (≳ 10%) of outer disk, carbonaceous material to the proto-Earth in the form of pebbles, in agreement with chemical and isotopic investigations of Earth’s accretion history.

MRI turbulence in vertically stratified accretion discs at large magnetic Prandtl numbers

ABSTRACT The discovery of the first binary neutron star merger, GW170817, has spawned a plethora of global numerical relativity simulations. These simulations are often ideal (with dissipation determined by the grid) and/or axisymmetric (invoking ad hoc mean-field dynamos). However, binary neutron star mergers (similar to X-ray binaries and active galactic nuclei inner discs) are characterized by large magnetic Prandtl numbers, $\rm Pm$, (the ratio of viscosity to resistivity). $\rm Pm$ is a key parameter determining dynamo action and dissipation but it is ill-defined (and likely of order unity) in ideal simulations. To bridge this gap, we investigate the magnetorotational instability (MRI) and associated dynamo at large magnetic Prandtl numbers using fully compressible, three-dimensional, vertically stratified, isothermal simulations of a local patch of a disc. We find that, within the bulk of the disc (z ≲ 2H, where H is the scale-height), the turbulent intensity (parametrized by the stress-to-thermal-pressure ratio α), and the saturated magnetic field energy density, Emag, produced by the MRI dynamo, both scale as a power with Pm at moderate Pm (4 ≲ Pm ≲ 32): Emag ∼ Pm0.74 and α ∼ Pm0.71, respectively. At larger Pm (≳ 32), we find deviations from power-law scaling and the onset of a plateau. Compared to our recent unstratified study, this scaling with Pm becomes weaker further away from the disc mid-plane, where the Parker instability dominates. We perform a thorough spectral analysis to understand the underlying dynamics of small-scale MRI-driven turbulence in the mid-plane and of large-scale Parker-unstable structures in the atmosphere.

JWST MIRI MRS Images of Disk Winds, Water, and CO in an Edge-on Protoplanetary Disk

We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young subsolar-mass star Tau 042021, acquired as part of the Cycle 1 GO program “Mapping Inclined Disk Astrochemical Signatures.” These data resolve the mid-IR spatial distributions of H2, revealing X-shaped emission extending to ∼200 au above the disk midplane with a semiopening angle of 35° ± 5°. We do not velocity-resolve the gas in the spectral images, but the measured semiopening angle of the H2 is consistent with a magnetohydrodynamic wind origin. A collimated, bipolar jet is seen in forbidden emission lines from [Ne ii], [Ne iii], [Ni ii], [Fe ii], [Ar ii], and [S iii]. Extended H2O and CO emission lines are also detected, reaching diameters of ∼90 and 190 au, respectively. Hot molecular emission is not expected at such radii, and we interpret its extended spatial distribution as scattering of inner disk molecular emission by dust grains in the outer disk surface. H i recombination lines, characteristic of inner disk accretion shocks, are similarly extended and are likely also scattered light from the innermost star–disk interface. Finally, we detect extended polycyclic aromatic hydrocarbon (PAH) emission at 11.3 μm cospatial with the scattered-light continuum, making this the first low-mass T Tauri star around which extended PAHs have been confirmed, to our knowledge. MIRI MRS line images of edge-on disks provide an unprecedented window into the outflow, accretion, and scattering processes within protoplanetary disks, allowing us to constrain the disk lifetimes and accretion and mass-loss mechanisms.

The XMM-Newton Line Emission Analysis Program (X-LEAP). II. The Multiscale Temperature Structures in the Milky Way Hot Gas

This paper presents the multiscale temperature structures in the Milky Way (MW) hot gas, as part of the XMM-Newton Line Emission Analysis Program, surveying the O vii, O viii, and Fe-L band emission features in the XMM-Newton archive. In particular, we define two temperature tracers, I OVIII/I OVII (O87) and I FeL/(I OVII + I OVIII) (FeO). These two ratios cannot be explained simultaneously using single-temperature collisional ionization models, which indicates the need for multitemperature structures in hot gas. In addition, we show three large-scale features in the hot gas: the eROSITA bubbles around the Galactic center (GC), the disk, and the halo. In the eROSITA bubbles, the observed line ratios can be explained by a log-normal temperature distribution with a median of logT/K≈6.4 and a scatter of σ T ≈ 0.2 dex. Beyond the bubbles, the line ratio dependence on the Galactic latitude suggests higher temperatures around the midplane of the MW disk. The scale height of the temperature variation is estimated to be ≈2 kpc assuming an average distance of 5 kpc for the hot gas. The halo component is characterized by the dependence on the distance to the GC, showing a temperature decline from logT/K≈6.3 to 5.8. Furthermore, we extract the autocorrelation and cross-correlation functions to investigate the small-scale structures. O87 and FeO ratios show a consistent autocorrelation scale of ≈5° (i.e., ≈400 pc at 5 kpc), which is consistent with the expected physical sizes of X-ray bubbles associated with star-forming regions or supernova remnants. Finally, we examine the cross-correlation between the hot and UV-detected warm gas and show an intriguing anticorrelation.