Search

We implement a local model for a spherical collapsing/expanding gas cloud into the Athena++ magnetohydrodynamic code. This local model consists of a Cartesian periodic box with time-dependent geometry. We present a series of benchmark test problems, including non-linear solutions and linear perturbations of the local model, confirming the code's desired performance. During a spherical collapse, a horizontal shear flow is amplified, corresponding to angular momentum conservation of zonal flows in the global problem; wave speed and amplitude of sound waves increase in the local frame, due to the reduction in the characteristic length scale of the box, which can lead to an anisotropic effective sound speed in the local box. Our code conserves both mass and momentum to machine precision. This numerical implementation of the local model has potential applications to the study of local physics and hydrodynamic instabilities during protostellar collapse, providing a powerful framework for better understanding the earliest stages of star and planet formation., Comment: 14 pages, 9 figures, plus appendix (3 pages, 1 figure). Accepted for publication in Astronomy & Astrophysics

Disk winds are thought to play a critical role in the evolution and dispersal of protoplanetary disks. A primary diagnostic of this physics is emission from the wind, especially in the low-velocity component of the [O I] $\lambda6300$ line. However, the interpretation of the line is usually based on spectroscopy alone, which leads to confusion between magnetohydrodynamic winds and photoevaporative winds. Here we report that in high-resolution spectral mapping of TW Hya by the multi-unit spectroscopic explorer at the Very Large Telescope, 80% of the [O I] emission is confined to within 1 AU radially from the star. A generic model of a magnetothermal wind produces [O I] emission at the base of the wind that broadly matches the flux and the observed spatial and spectral profiles. The emission at large radii is much fainter that predicted from models of photoevaporation, perhaps because the magnetothermal wind partially shields the outer disk from energetic radiation from the central star. This result calls into question the previously assessed importance of photoevaporation in disk dispersal predicted by models., Comment: Published in Nature Astronomy

Mammalian brain is a complex organ that contains billions of neurons. These neurons form various neural circuits that control the perception, cognition, emotion and behavior. Developing in vivo neuronal labeling and imaging techniques is crucial for studying the structure and function of neural circuits. In vivo techniques can provide true physiological information that cannot be provided by ex vivo methods. In this study, we describe a new strategy for in vivo neuronal labeling and quantification using MRI. To demonstrate the ability of this new method, we used neurotropic virus to deliver oatp1a1 gene to the target neural circuit. OATP1A1 protein is expressed on the neuronal membrane and can increase the uptake of a specific MRI contrast agent (Gd-EOB-DTPA). By using T1-weighted images for observation, labeled neurons "light up" on MRI. We further use a dynamic-contrast-enhancement based method to obtain measures that provide quantitative information of labeled neurons in vivo.

Recent observations of protoplanetary disks (PPDs) in the sub-mm have revealed the ubiquity of annular substructures, indicative of pebble-sized dust particles trapped in turbulent ring-like gas pressure bumps. This major paradigm shift also challenges the leading theory of planetesimal formation from such pebbles by the streaming instability, which operates in a pressure gradient and can be suppressed by turbulence. Here we conduct three-dimensional local shearing-box, non-ideal magnetohydrodynamic (MHD) simulations of dust trapping in enforced gas pressure bumps including dust backreaction. Under the moderate level of turbulence generated by the magnetorotational instability (MRI) with ambipolar diffusion that is suitable for outer disk conditions, we achieve quasi-steady states of dust trapping balanced by turbulent diffusion. We find strong dust clumping in all simulations near the gas pressure maxima, reaching a maximum density well above the threshold of triggering gravitational collapse to form planetesimals. A strong pressure bump concentrates dust particles towards bump center. With a weak pressure bump, dust can also concentrate in secondary filaments off the bump center due to dust backreaction, but strong clumping still occurs mainly in the primary ring around the bump center. Our results reveal dust-trapping rings as robust locations for planetesimal formation in outer PPDs, while they may possess diverse observational properties., Comment: 13 pages, 4 figures, plus appendix (4 pages, 4 figures). Accepted for publication in ApJ Letter

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from ALMA and JWST, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, NIR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry., Comment: accepted to ApJ

Planetesimal formation is a crucial yet poorly understood process in planet formation. It is widely believed that planetesimal formation is the outcome of dust clumping by the streaming instability (SI). However, recent analytical and numerical studies have shown that the SI can be damped or suppressed by external turbulence, and at least the outer regions of protoplanetary disks are likely weakly turbulent due to magneto-rotational instability (MRI). We conduct high-resolution local shearing-box simulations of hybrid particle-gas magnetohydrodynamics (MHD), incorporating ambipolar diffusion as the dominant non-ideal MHD effect, applicable to outer disk regions. We first show that dust backreaction enhances dust settling towards the midplane by reducing turbulence correlation time. Under modest level of MRI turbulence, we find that dust clumping is in fact easier than the conventional SI case, in the sense that the threshold of solid abundance for clumping is lower. The key to dust clumping includes dust backreaction and the presence of local pressure maxima, which in our work is formed by the MRI zonal flows overcoming background pressure gradient. Overall, our results support planetesimal formation in the MRI-turbulent outer protoplanetary disks, especially in ring-like substructures., Comment: 21 pages, 9 figures, accepted for publication in ApJ

We present an analysis of wind absorption in the C II ${\lambda}1335$ doublet towards 40 classical T Tauri stars with archival far-ultraviolet (FUV) spectra obtained by the Hubble Space Telescope. Absorption features produced by fast or slow winds are commonly detected (36 out of 40 targets) in our sample. The wind velocity of the fast wind decreases with disk inclination, consistent with expectations for a collimated jet. Slow wind absorption is detected mostly in disks with intermediate or high inclination, without a significant dependence of wind velocity on disk inclination. Both the fast and slow wind absorption are preferentially detected in FUV lines of neutral or singly ionized atoms. The Mg II ${\lambda}{\lambda}2796,2804$ lines show wind absorption consistent with the absorption in the C II lines. We develop simplified semi-analytical disk/wind models to interpret the observational disk wind absorption. Both fast and slow winds are consistent with expectations from a thermal-magnetized disk wind model and are generally inconsistent with a purely thermal wind. Both the models and the observational analysis indicate that wind absorption occurs preferentially from the inner disk, offering a wind diagnostic in complement to optical forbidden line emission that traces the wind in larger volumes., Comment: 27 pages, 16 figures, plus appendix. Accepted for publication in ApJ

We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring., Comment: Accepted for publication in the Astrophysical Journal

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