Search

While protoplanetary disks (PPDs) are generally thought to dissipate within several Myr, recent observations have revealed gas in debris disks. The origin of this gas remains uncertain, with one possibility being the unexpectedly long survival of PPDs (the primordial-origin scenario). To explore the plausibility of this scenario, we conduct 1D disk evolution simulations, varying parameters like stellar mass, disk mass, turbulent stress, and magnetohydrodynamic winds, while incorporating stellar evolution to account for time-varying photoevaporation rates. Our focus is on disks where small grains are depleted, as these are potentially long-lived due to reduced far-ultraviolet photoevaporation. Our results show that gas in these disks can survive beyond 10 Myr regardless of the stellar mass, provided they are initially massive ($M_{\mathrm{disk}}\approx 0.1M_*$) with relatively weak turbulent stress ($\alpha \ll 10^{-2}$). The longest lifetimes are consistently found for $M_* = 2 M_{\odot}$ across a wide parameter space, with gas typically persisting at $\sim 10$--$10^3 \mathrm{au}$. Roughly estimated CO masses for these disks fall within the observed range for the most massive gas-rich debris disks around early A stars. These alignments support the plausibility of the primordial-origin scenario. Additionally, our model predicts that accretion persists for as long as the disk survives, which could explain the accretion signatures detected in old disks hosted by low-mass stars, including Peter Pan disks. Our finding also suggests that ongoing accretion may exist in gas-rich debris disks. Thus, searching for accretion signatures could be a key factor to identifying the origin of gas in debris disks., Comment: 20 pages, 16 figures

Ultraviolet flare emission can drive photochemistry in exoplanet atmospheres and even serve as the primary source of uncertainty in atmospheric retrievals. Additionally, flare energy budgets are not well-understood due to a paucity of simultaneous observations. We present new near-UV (NUV) and optical observations of flares from three M dwarfs obtained at 20 s cadence with Swift and TESS, along with a re-analysis of flares from two M dwarfs in order to explore the energy budget and timing of flares at NUV--optical wavelengths. We find a 9000 K blackbody underestimates the NUV flux by $\geq$2$\times$ for 54$\pm$14% of flares and 14.8$\times$ for one flare. We report time lags between the bands of 0.5--6.6 min and develop a method to predict the qualitative flare shape and time lag to 36$\pm$30% accuracy. The scatter present in optical-NUV relations is reduced by a factor of 2.0$\pm$0.6 when comparing the total NUV energy with the TESS energy during the FWHM duration due to the exclusion of the $T_\mathrm{eff}\approx$5000 K tail. We show the NUV light curve can be used to remove flares from the optical light curve and consistently detect planets with 20% smaller transits than is possible without flare detrending. Finally, we demonstrate a 10$\times$ increase in the literature number of multi-wavelength flares with the Early eVolution Explorer (EVE), an astrophysics Small Explorer concept to observe young clusters with simultaneous NUV and optical bands in order to detect young planets, assess their photochemical radiation environments, and observe accretion., Comment: 27 pages, 13 figures, 4 tables, accepted to The Astronomical Journal

Molecular clouds (MCs) are the birthplaces of new stars in galaxies. A key component of MCs are photodissociation regions (PDRs), where far-ultraviolet radiation plays a crucial role in determining the gas's physical and chemical state. Traditional PDR models assume chemical steady state (CSS), where the rates of H$_2$ formation and photodissociation are balanced. However, real MCs are dynamic and can be out of CSS. In this study, we demonstrate that combining H$_2$ emission lines observed in the far-ultraviolet or infrared with column density observations can be used to derive the rates of H$_2$ formation and photodissociation. We derive analytical formulae that relate these rates to observable quantities, which we validate using synthetic H$_2$ line emission maps derived from the SILCC-Zoom hydrodynamical simulation. Our method estimates integrated H$_2$ formation and dissociation rates to within 29\% accuracy. Our simulations cover a wide dynamic range in H$_2$ formation and photodissociation rates, showing significant deviations from CSS, with 74\% of the MC's mass deviating from CSS by a factor greater than 2. Our analytical formulae can effectively distinguish between regions in and out of CSS. When applied to actual H$_2$ line observations, our method can assess the chemical state of MCs, providing insights into their evolutionary stages and lifetimes., Comment: Submitted to ApJ. **Comments Welcome!** This paper is the first ("Paper-I") in 2-paper series. Paper II, entitled "The Molecular Cloud Lifecycle II: Formation and Destruction of Molecular Clouds Diagnosed via H 2 Fluorescent Emission Emission" by Burkhart et al has been recently accepted for publication in ApJ, and may be found here: arXiv:2402.01587

Photoevaporation driven by hydrogen-ionizing radiation, also known as extreme-ultraviolet (EUV), profoundly shapes the lives of diverse astrophysical objects. Focusing here mainly on the dispersal of protoplanetary disks, we construct an analytical model accounting for the finite timescales of photoheating and photoionization. The model offers improved estimates for the ionization, temperature, and velocity structures versus distance from the central source, for a given EUV emission rate and spectral hardness. Compared to the classical picture of fully-ionized and isothermal winds with temperatures $\approx 10^4{\rm \,K}$ and speeds $\approx 10{\rm \,km\,s^{-1}}$, our model unveils broader hydrodynamical and thermochemical states of photoevaporative winds. In contrast to the classical picture, T~Tauri stars with EUV luminosities around $10^{30}{\rm \,erg\,s^{-1}}$ have non-isothermal ionized winds at lower temperatures than the classical value if the spectrum is soft, with an average deposited energy per photoionization less than about 3.7\,eV. Conversely, if the spectrum is hard, the winds tend to be atomic and isothermal at most radii in the disk. For lower EUV intensities, even with soft spectra, atomic winds can emerge beyond $\sim 10{\, \rm au}$ through advection. We demonstrate that the analytical model's predictions are in general agreement with detailed radiation-hydrodynamics calculations. The model furthermore illustrates how the energy efficiency of photoevaporation varies with the intensity and spectral hardness of the EUV illumination, as well as addressing discrepancies in the literature around the effectiveness of X-ray photoevaporation. These findings highlight the importance of considering the finite timescales of photoheating and photoionization, both in modeling and in interpreting observational data., Comment: Submitted to ApJ

Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ formation, stellar radiation, sink particles, and turbulence. We find that, at early times in the cloud evolution, H$_2$ formation rapidly outpaces dissociation and molecular clouds build their mass from the atomic reservoir in their environment. Rapid H$_2$ formation is also associated with a higher early star formation rate. For the clouds studied here, H$_2$ is strongly out of chemical equilibrium during the early stages of cloud formation but settles into a bursty chemical steady-state about 2 Myrs after the first stars form. At the latest stage of cloud evolution, dissociation outweighs formation and the clouds enter a dispersal phase. We discuss how theories for the molecular cloud lifecycle and the star formation efficiency may be distinguished with observational measurements of H$_2$ fluorescence with a space-based high-resolution FUV spectrometer, such as the proposed Hyperion and Eos NASA Explorer missions. Such missions would enable measurements of the H$_2$ dissociation and formation rates, which we demonstrate can be connected to different phases in a molecular cloud's star-forming life, including cloud building, rapidly star-forming, H$_2$ chemical equilibrium, and cloud destruction., Comment: Submitted to ApJ, comments welcome

While most debris disks consist of dust with little or no gas, a fraction has significant amounts of gas detected via emission lines of CO, ionized carbon, and/or atomic oxygen. Almost all such gaseous debris disks known are around A-type stars with ages up to 50 Myr. We show, using semi-analytic disk evolution modeling, that this can be understood if the gaseous debris disks are remnant protoplanetary disks that have become depleted of small grains compared to the interstellar medium. Photoelectric heating by the A-stars' FUV radiation is then inefficient, while the stars' EUV and X-ray emissions are weak owing to a lack of surface convective zones capable of driving magnetic activity. In this picture, stars outside the range of spectral types from A through early B are relatively hard to have such long-lived gas disks. Less-massive stars have stronger magnetic activity in the chromosphere, transition region, and corona with resulting EUV and X-ray emission, while more-massive stars have photospheres hot enough to produce strong EUV radiation. In both cases, primordial disk gas is likely to photoevaporate well before 50 Myr. These results come from 0D disk evolution models where we incorporate internal accretion stresses, MHD winds, and photoevaporation by EUV and X-ray photons with luminosities that are functions of the stellar mass and age. A key issue this work leaves open is how some disks become depleted in small dust so that FUV photoevaporation slows. Candidates include grains' growth, settling, radial drift, radiation force, and incorporation into planetary systems., Comment: Submitted to ApJL

Modeling self-gravitating gas flows is essential to answering many fundamental questions in astrophysics. This spans many topics including planet-forming disks, star-forming clouds, galaxy formation, and the development of large-scale structures in the Universe. However, the nonlinear interaction between gravity and fluid dynamics offers a formidable challenge to solving the resulting time-dependent partial differential equations (PDEs) in three dimensions (3D). By leveraging the universal approximation capabilities of a neural network within a mesh-free framework, physics informed neural networks (PINNs) offer a new way of addressing this challenge. We introduce the gravity-informed neural network (GRINN), a PINN-based code, to simulate 3D self-gravitating hydrodynamic systems. Here, we specifically study gravitational instability and wave propagation in an isothermal gas. Our results match a linear analytic solution to within 1\% in the linear regime and a conventional grid code solution to within 5\% as the disturbance grows into the nonlinear regime. We find that the computation time of the GRINN does not scale with the number of dimensions. This is in contrast to the scaling of the grid-based code for the hydrodynamic and self-gravity calculations as the number of dimensions is increased. Our results show that the GRINN computation time is longer than the grid code in one- and two- dimensional calculations but is an order of magnitude lesser than the grid code in 3D with similar accuracy. Physics-informed neural networks like GRINN thus show promise for advancing our ability to model 3D astrophysical flows.

The solar system's Zodiacal Cloud is visible to the unaided eye, yet the origin of its constituent dust particles is not well understood, with a wide range of proposed divisions between sources in the asteroid belt and Jupiter Family comets. The amount of dust contributed by Oort Cloud comets is uncertain. Knowledge of the Zodiacal Cloud's structure and origins would help with NASA's aim of characterizing potentially Earth-like planets around nearby stars, since the exo-Earths must be studied against the light scattered from extrasolar analogs of our cloud. As the only example where the parent bodies can be tracked, our own cloud is critical for learning how planetary system architecture governs the interplanetary dust's distribution. Our cloud has been relatively little-studied in the near-ultraviolet, a wavelength range that is important for identifying potentially-habitable planets since it contains the broad Hartley absorption band of ozone. We show through radiative transfer modeling that our cloud's shape and size at near-UV wavelengths can be measured from Earth orbit by mapping the zodiacal light's flux and linear polarization across the sky. We quantify how well the cloud's geometric and optical properties can be retrieved from a set of simulated disk observations, using a Markov chain Monte Carlo analysis. The results demonstrate that observations with sufficient precision, covering a set of fields distributed along the ecliptic and up to the poles, can be used to determine the division between asteroidal, Jupiter Family, and Oort Cloud dust components, primarily via their differing orbital inclination distributions. We find that the observations must be repeated over a time span of several months in order to disentangle the zodiacal light from the Galactic background using the Milky Way's rotation across the sky., Comment: submitted to PASP

We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful far-ultraviolet (FUV) diagnostic. Molecular hydrogen (H2) is the most abundant molecule in the universe and a key ingredient for star and planet formation, but is typically not observed directly because its symmetric atomic structure and lack of a dipole moment mean there are no spectral lines at visible wavelengths and few in the infrared. Hyperion uses molecular hydrogen's wealth of FUV emission lines to achieve three science objectives: (1) determining how star formation is related to molecular hydrogen formation and destruction at the boundaries of molecular clouds; (2) determining how quickly and by what process massive star feedback disperses molecular clouds; and (3) determining the mechanism driving the evolution of planet-forming disks around young solar-analog stars. Hyperion conducts this science using a straightforward, highly-efficient, single-channel instrument design. Hyperion's instrument consists of a 48 cm primary mirror, with an f/5 focal ratio. The spectrometer has two modes, both covering 138.5-161.5 nm bandpasses. A low resolution mode has a spectral resolution of R>10,000 with a slit length of 65 arcmin, while the high resolution mode has a spectral resolution of R>50,000 over a slit length of 5 armin. Hyperion occupies a 2 week long, high-earth, Lunar resonance TESS-like orbit, and conducts 2 weeks of planned observations per orbit, with time for downlinks and calibrations. Hyperion was reviewed as Category I, which is the highest rating possible, but was not selected., Comment: Accepted to JATIS, 9 Figures

We model substructure in the protoplanetary disks around DoAr 44 and HD 163296 in order to better understand the conditions under which planets may form. We match archival millimeter-wavelength thermal emission against models of the disks' structure that are in radiation balance with the starlight heating and in vertical hydrostatic equilibrium, and then compare to archival polarized scattered near-infrared images of the disks. The millimeter emission arises in the interior, while the scattered near-infrared radiation probes the disks' outer layers. Our best model of the HD 163296 disk has dust masses $81\pm 13$ $M_\oplus$ in the inner ring at 68 au and $82^{+26}_{-16}$ $M_\oplus$ in the outer ring at 102 au, both falling within the range of estimates from previous studies. Our DoAr 44 model has total dust mass $84^{+7.0}_{-3.5}$ $M_\oplus$. Unlike HD 163296, DoAr 44 as of yet has no detected planets. If the central cavity in the DoAr 44 disk is caused by a planet, the planet's mass must be at least 0.5 $M_J$ and is unlikely to be greater than 1.6 $M_J$. We demonstrate that the DoAr 44 disk's structure with a bright ring offset within a fainter skirt can be formed by dust particles drifting through a plausible distribution of gas., Comment: 18 pages, 19 figures

Understanding the thermal structure of protoplanetary disks is crucial for modeling planet formation and interpreting disk observations. We present a new two-layer radiative transfer model for computing the thermal structure of axisymmetric irradiated disks. Unlike the standard two-layer model, our model accounts for the radial as well as vertical transfer of the starlight reprocessed at the disk surface. The model thus allows us to compute the temperature below "shadowed" surfaces receiving no direct starlight. Thanks to the assumed axisymmetry, the reprocessed starlight flux is given in one-dimensional integral form that can be computed at a low cost. Furthermore, our model evolves the midplane temperature using a time-dependent energy equation and can therefore treat thermal instabilities. We apply our global two-layer model to disks with a planetary induced gap and confirm that the model reproduces the disks' temperature profiles obtained from more computationally expensive Monte Carlo radiative transfer calculations to an accuracy of less than 20%. We also apply the model to study the long-term behavior of the thermal wave instability in irradiated disks. Being simple and computationally efficient, the global two-layer model will be suitable for studying the interplay between disks' thermal evolution and dust evolution., Comment: 22 pages, 18 figures, accepted for publication in PASJ (Publications of the Astronomical Society of Japan)

Understanding the origin of accretion and dispersal of protoplanetary disks is fundamental for investigating planet formation. Recent numerical simulations show that launching winds are unavoidable when disks undergo magnetically driven accretion and/or are exposed to external UV radiation. Observations also hint that disk winds are common. We explore how the resulting wind mass loss rate can be used as a probe of both disk accretion and dispersal. As a proof-of-concept study, we focus on magnetocentrifugal winds, MRI (magnetorotational instability) turbulence, and external photoevapotaion. By developing a simple, yet physically motivated disk model and coupling it with simulation results available in the literature, we compute the mass loss rate as a function of external UV flux for each mechanism. We find that different mechanisms lead to different levels of mass loss rate, indicating that the origin of disk accretion and dispersal can be determined, by observing the wind mass loss rate resulting from each mechanism. This determination provides important implications for planet formation. This work thus shows that the ongoing and future observations of the wind mass loss rate for protoplanetary disks are paramount to reliably constrain how protoplanetary disks evolve with time and how planet formation takes place in the disks., Comment: 16 pages; 7 figures, 5 tables, accepted for publication in ApJL; no change in the conclusions; the additional Monte-Carlo based, parameter study was conducted, following the referee's comment

The recent high spatial/spectral resolution observations have enabled constraining formation mechanisms of giant planets, especially at the final stages. The current interpretation of such observations is that these planets undergo magnetospheric accretion, suggesting the importance of planetary magnetic fields. We explore the properties of accreting, magnetized giant planets surrounded by their circumplanetary disks, using the physical parameters inferred for PDS 70 b/c. We compute the magnetic field strength and the resulting spin rate of giant planets, and find that these planets may possess dipole magnetic fields of either a few 10 G or a few 100 G; the former is the natural outcome of planetary growth and radius evolution, while the resulting spin rate cannot reproduce the observations. For the latter, a consistent picture can be drawn, where strong magnetic fields induced by hot planetary interiors lead both to magnetospheric accretion and to spin-down due to disk locking. We also compute the properties of circumplanetary disks in the vicinity of these planets, taking into account planetary magnetic fields. The resulting surface density becomes very low, compared with the canonical models, implying the importance of radial movement of satellite-forming materials. Our model predicts a positive gradient of the surface density, which invokes the traps for both satellite migration and radially drifting dust particles. This work thus concludes that the final formation stages of giant planets are similar to those of low-mass stars such as brown dwarfs, as suggested by recent studies., Comment: 13 pages, 1 table, 2 figures, accepted for publication in ApJ

Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$\mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $\alpha \simeq 10^{-4}$. The value of $\alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $\alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation., Comment: ApJ accepted

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

The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems., Comment: A&A in press, MPIA press release 10.10.2019

Hydrodynamical simulations show that a pair of spiral arms can form in the disk around a rapidly-growing young star and that the arms are crucial in transporting angular momentum as the disk accretes material from the surrounding envelope. Here we report the detection of a pair of symmetric spiral structures in a protostellar disk, supporting the formation of spiral arms in the disk around a forming star. The HH 111 VLA 1 source is a young Class I source embedded in a massive infalling protostellar envelope and is actively accreting, driving the prominent HH 111 jet. Previous observations showed a ring of shock emission around the disk's outer edge, indicating accretion of the envelope material onto the disk at a high rate. Now with ALMA observations of thermal emission from dust particles, we detect a pair of spiral arms extending from the inner region to the disk's outer edge, similar to that seen in many simulations. Additionally, the disk is massive, with Toomre's Q parameter near unity in the outer parts where the spiral structures are detected, supporting the notion that envelope accretion is driving the outer disk gravitationally unstable. In our observations, another source, HH 111 VLA 2, is spatially resolved for the first time, showing a disk-like structure with a diameter of ~ 26 au and an orientation nearly orthogonal to that of the HH 111 VLA 1 disk., Comment: 24 pages, 3 figures

In spite of making a small contribution to total protoplanetary disk mass, dust affects the disk temperature by controlling absorption of starlight. As grains grow from their initial ISM-like size distribution, settling depletes the disk's upper layers of dust and decreases the optical depth, cooling the interior. Here we investigate the effect of collisional growth of dust grains and their dynamics on the thermal and optical profile of the disk, and explore the possibility that cooling induced by grain growth and settling could lead to gravitational instability. We develop a Monte Carlo dust collision model with a weighting technique and allow particles to collisionally evolve through sticking and fragmentation, along with vertical settling and turbulent mixing. We explore two disk models, the MMEN (minimum-mass extrasolar nebula), and a "heavy" disk with higher surface density than the MMEN, and perform simulations for both constant and spatially variable turbulence efficiency profile $\alpha(R,z)$. We then calculate mean wavelength-dependent opacities for the evolving disks and perform radiative transfer to calculate the temperature profile $T(R,z)$. Finally, we calculate the Toomre Q parameter, a measure of the disk's stability against self-gravity, for each disk model after it reaches a steady state dust-size distribution. We find that even weak turbulence can keep sub-micron sized particles stirred in the disk's upper layer, affecting its optical and thermal profiles, and the growth of large particles in the midplane can make a massive disk optically thick at millimeter wavelengths, making it difficult to calculate the surface density of dust available for planet formation in the inner disk. Also, for an initially massive disk, grain settling and growth can produce a drop in the Toomre Q parameter, driving the disk to $Q < 1.4$ and possibly triggering spiral instabilities., Comment: 27 pages, 21 Figures, Under Review (ApJ)

Young planets with masses approaching Jupiter's have tides strong enough to clear gaps around their orbits in the protostellar disk. Gas flow through the gaps regulates the planets' further growth and governs the disks' evolution. Magnetic forces may drive that flow if the gas is sufficiently ionized to couple to the fields. We compute the ionizing effects of the X-rays from the central young star, using Monte Carlo radiative transfer calculations to find the spectrum of Compton-scattered photons reaching the planet's vicinity. The scattered X-rays ionize the gas at rates similar to or greater than the interstellar cosmic ray rate near planets the mass of Saturn and of Jupiter, located at 5 au and at 10 au, in disks with the interstellar mass fraction of sub-micron dust and with the dust depleted a factor 100. Solving a gas-grain recombination reaction network yields charged particle populations whose ability to carry currents is sufficient to partly couple the magnetic fields to the gas around the planet. Most cases can undergo Hall shear instability, and some can launch magnetocentrifugal winds. However the material on the planet's orbit has diffusivities so large in all the cases we examine, that magneto-rotational turbulence is prevented and the non-ideal terms govern the magnetic field's evolution. Thus the flow of gas in the gaps opened by the young giant planets depends crucially on the finite conductivity., Comment: 24 pages, 6 figures. Gap depths are now chosen to match recent hydrodynamical results. Accepted for publication in ApJ

HH 212 is a Class 0 protostellar system found to host a "hamburger"-shaped dusty disk with a rotating disk atmosphere and a collimated SiO jet at a distance of ~ 400 pc. Recently, a compact rotating outflow has been detected in SO and SO2 toward the center along the jet axis at ~ 52 au (0.13") resolution. Here we resolve the compact outflow into a small-scale wide-opening rotating outflow shell and a collimated jet, with the observations in the same S-bearing molecules at ~ 16 au (0.04") resolution. The collimated jet is aligned with the SiO jet, tracing the shock interactions in the jet. The wide-opening outflow shell is seen extending out from the inner disk around the SiO jet and has a width of ~ 100 au. It is not only expanding away from the center, but also rotating around the jet axis. The specific angular momentum of the outflow shell is ~ 40 au km/s. Simple modeling of the observed kinematics suggests that the rotating outflow shell can trace either a disk wind or disk material pushed away by an unseen wind from the inner disk or protostar. We also resolve the disk atmosphere in the same S-bearing molecules, confirming the Keplerian rotation there., Comment: 29 pages, 11 figures

Planets are born in protostellar disks, which are now observed with enough resolution to address questions about internal gas flows. Candidates for driving the flows include magnetic forces, but ionization state estimates suggest much of the gas mass decouples from magnetic fields. Thus, hydrodynamical instabilities could play a major role. We investigate disk dynamics under conditions typical for a T Tauri system, using global 3D radiation hydrodynamics simulations with embedded particles and a resolution of 70 cells per scale height. Stellar irradiation heating is included with realistic dust opacities. The disk starts in joint radiative balance and hydrostatic equilibrium. The vertical shear instability (VSI) develops into turbulence that persists up to at least 1600 inner orbits (143 outer orbits). Turbulent speeds are a few percent of the local sound speed at the midplane, increasing to 20%, or 100 m/s, in the corona. These are consistent with recent upper limits on turbulent speeds from optically thin and thick molecular line observations of TW Hya and HD 163296. The predominantly vertical motions induced by the VSI efficiently lift particles upwards. Grains 0.1 and 1 mm in size achieve scale heights greater than expected in isotropic turbulence. We conclude that while kinematic constraints from molecular line emission do not directly discriminate between magnetic and nonmagnetic disk models, the small dust scale heights measured in HL Tau and HD 163296 favor turbulent magnetic models, which reach lower ratios of the vertical kinetic energy density to the accretion stress., Comment: Accepted for ApJ

Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s., Comment: 12 pages, including 10 figures. Accepted for publication in Astrophysical Journal. Revised following comments from community

The mechanism of angular momentum transport in protoplanetary disks is fundamental to understand the distributions of gas and dust in the disks. The unprecedented, high spatial resolution ALMA observations taken toward HL Tau and subsequent radiative transfer modeling reveal that a high degree of dust settling is currently achieved at the outer part of the HL Tau disk. Previous observations however suggest a high disk accretion rate onto the central star. This configuration is not necessarily intuitive in the framework of the conventional viscous disk model, since efficient accretion generally requires a high level of turbulence, which can suppress dust settling considerably. We develop a simplified, semi-analytical disk model to examine under what condition these two properties can be realized in a single model. Recent, non-ideal MHD simulations are utilized to realistically model the angular momentum transport both radially via MHD turbulence and vertically via magnetically induced disk winds. We find that the HL Tau disk configuration can be reproduced well when disk winds are properly taken into account. While the resulting disk properties are likely consistent with other observational results, such an ideal situation can be established only if the plasma $\beta$ at the disk midplane is $\beta_0 \simeq 2 \times 10^4$ under the assumption of steady accretion. Equivalently, the vertical magnetic flux at 100 au is about 0.2 mG. More detailed modeling is needed to fully identify the origin of the disk accretion and quantitatively examine plausible mechanisms behind the observed gap structures in the HL Tau disk., Comment: 14 pages, 9 figures; submitted to ApJ; the revision has been made, following the referee report

Gas in protostellar disks provides the raw material for giant planet formation and controls the dynamics of the planetesimal-building dust grains. Accurate gas mass measurements help map the observed properties of planet-forming disks onto the formation environments of known exoplanets. Rare isotopologues of carbon monoxide (CO) have been used as gas mass tracers for disks in the Lupus star-forming region, with an assumed interstellar CO/H$_2$ abundance ratio. Unfortunately, observations of T-Tauri disks show that CO abundance is not interstellar---a finding reproduced by models that show CO abundance decreasing both with distance from the star and as a function of time. Here we present radiative transfer simulations that assess the accuracy of CO-based disk mass measurements. We find that the combination of CO chemical depletion in the outer disk and optically thick emission from the inner disk leads observers to underestimate gas mass by more than an order of magnitude if they use the standard assumptions of interstellar CO/H$_2$ ratio and optically thin emission. Furthermore, CO abundance changes on million-year timescales, introducing an age/mass degeneracy into observations. To reach factor of a few accuracy for CO-based disk mass measurements, we suggest that observers and modelers adopt the following strategies: (1) select the low-$J$ transitions; (2) observe multiple CO isotopologues and use either intensity ratios or normalized line profiles to diagnose CO chemical depletion; and (3) use spatially resolved observations to measure the CO abundance distribution., Comment: 22 pages, 16 figures, Accepted for publication in ApJ

We perform two-dimensional axially symmetric radiation-hydrodynamic simulations to assess the impact of outflows and radiative force feedback from massive protostars by varying when the protostellar outflow starts, the ratio of ejection to accretion rates, and the strength of the wide angle disk wind component. The star formation efficiency, i.e. the ratio of final stellar mass to initial core mass, is dominated by radiative forces and the ratio of outflow to accretion rates. Increasing this ratio has three effects: First, the protostar grows slower with a lower luminosity at any given time, lowering radiative feedback. Second, bipolar cavities cleared by the outflow are larger, further diminishing radiative feedback on disk and core scales. Third, the higher momentum outflow sweeps up more material from the collapsing envelope, decreasing the protostar's potential mass reservoir via entrainment. The star formation efficiency varies with the ratio of ejection to accretion rates from 50% in the case of very weak outflows to as low as 20% for very strong outflows. At latitudes between the low density bipolar cavity and the high density accretion disk, wide angle disk winds remove some of the gas, which otherwise would be part of the accretion flow onto the disk; varying the strength of these wide angle disk winds, however, alters the final star formation efficiency by only +/-6%. For all cases, the opening angle of the bipolar outflow cavity remains below 20 degree during early protostellar accretion phases, increasing rapidly up to 65 degree at the onset of radiation pressure feedback., Comment: accepted for publication at ApJ

The gas near the midplanes of planet-forming protostellar disks remains largely unprobed by observations due to the high optical depth of commonly observed molecules such as CO and H$_2$O. However, rotational emission lines from rare molecules may have optical depths near unity in the vertical direction, so that the lines are strong enough to be detected, yet remain transparent enough to trace the disk midplane. Here we present a chemical model of an evolving T-Tauri disk and predict the optical depths of rotational transitions of $^{12}$C$^{16}$O, $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O. The MRI-active disk is primarily heated by the central star due to the formation of the dead zone. CO does not freeze out in our modeled region within $70$AU around a sunlike star. However, the abundance of CO decreases because of the formation of complex organic molecules (COM), producing an effect that can be misinterpreted as the "snow line". These results are robust to variations in our assumptions about the evolution of the gas to dust ratio. The optical depths of low-order rotational lines of C$^{17}$O are around unity, making it possible to see into the disk midplane using C$^{17}$O. Combining observations with modeled C$^{17}$O$/$H$_2$ ratios, like those we provide, can yield estimates of protoplanetary disks' gas masses., Comment: Accepted for publication in ApJ

Chondritic meteorites provide valuable opportunities to investigate the origins of the solar system. We explore impact jetting as a mechanism of chondrule formation and subsequent pebble accretion as a mechanism of accreting chondrules onto parent bodies of chondrites, and investigate how these two processes can account for the currently available meteoritic data. We find that when the solar nebula is $\le 5$ times more massive than the minimum-mass solar nebula at $a \simeq 2-3$ AU and parent bodies of chondrites are $\le 10^{24}$ g ($\le$ 500 km in radius) in the solar nebula, impact jetting and subsequent pebble accretion can reproduce a number of properties of the meteoritic data. The properties include the present asteroid belt mass, the formation timescale of chondrules, and the magnetic field strength of the nebula derived from chondrules in Semarkona. Since this scenario requires a first generation of planetesimals that trigger impact jetting and serve as parent bodies to accrete chondrules, the upper limit of parent bodies' masses leads to the following implications: primordial asteroids that were originally $\ge 10^{24}$ g in mass were unlikely to contain chondrules, while less massive primordial asteroids likely had a chondrule-rich surface layer. The scenario developed from impact jetting and pebble accretion can therefore provide new insights into the origins of the solar system., Comment: Title are modified and typographical errors are corrected, following the published version; 6 pages, 2 figures, 1 table, accepted for publication in ApJ Letters

We provide CoRoT and Spitzer light curves, as well as broad-band multi-wavelength photometry and high resolution, multi- and single-epoch spectroscopy for 17 classical T Tauris in NGC 2264 whose CoRoT light curves (LCs) exemplify the "stochastic" LC class as defined in Cody et al. (2014). The most probable physical mechanism to explain the optical variability in this LC class is time-dependent mass accretion onto the stellar photosphere, producing transient hot spots. As evidence in favor of this hypothesis, multi-epoch high resolution spectra for a subset of these stars shows that their veiling levels also vary in time and that this veiling variability is consistent in both amplitude and timescale with the optical LC morphology. Furthermore, the veiling variability is well-correlated with the strength of the HeI 6678A emission line, a feature predicted by models to arise in accretion shocks on or near the stellar photosphere. Stars with accretion burst LC morphology (Stauffer et al. 2014) are also attributed to variable mass accretion. Both the stochastic and accretion burst LCs can be explained by a simple model of randomly occurring flux bursts, with the stochastic LC class having a higher frequency of lower amplitude events. Based on their UV excesses, veiling, and mean Ha equivalent widths, members of the stochastic LC class have only moderate time-averaged mass accretion rates. The most common feature of their Ha profiles is for them to exhibit blue-shifted absorption features, most likely originating in a disk wind. The lack of periodic signatures in the LCs suggests that little of the variability is due to long-lived hot spots rotating into or out of our line of sight; instead, the primary driver of the observed photometric variability is likely to be instabilities in the inner disk that lead to variable mass accretion., Comment: 36 pages, 19 figures. Accepted to AJ

This is a joint summary of the reports from the three Astrophysics Program Analysis Groups (PAGs) in response to the "Planning for the 2020 Decadal Survey" charge given by the Astrophysics Division Director Paul Hertz. This joint executive summary contains points of consensus across all three PAGs. Additional findings specific to the individual PAGs are reported separately in the individual reports. The PAGs concur that all four large mission concepts identified in the white paper as candidates for maturation prior to the 2020 Decadal Survey should be studied in detail. These include the Far-IR Surveyor, the Habitable-Exoplanet Imaging Mission, the UV/Optical/IR Surveyor, and the X-ray Surveyor. This finding is predicated upon assumptions outlined in the white paper and subsequent charge, namely that 1) major development of future large flagship missions under consideration are to follow the implementation phases of JWST and WFIRST; 2) NASA will partner with the European Space Agency on its L3 Gravitational Wave Surveyor; 3) the Inflation Probe be classified as a probe-class mission to be developed according to the 2010 Decadal Survey report. If these key assumptions were to change, this PAG finding would need to be re-evaluated. The PAGs find that there is strong community support for the second phase of this activity - maturation of the four proposed mission concepts via Science and Technology Definition Teams (STDTs). The PAGs find that there is strong consensus that all of the STDTs contain broad and interdisciplinary representation of the science community. Finally, the PAGs find that there is community support for a line of Probe-class missions within the Astrophysics mission portfolio (condensed)., Comment: 22 pages, 2 tables, no figures

We calculate the abundances of electrons and ions in the hot (> 500 K), dusty parts of protoplanetary disks, treating for the first time the effects of thermionic and ion emission from the dust grains. High-temperature ionization modeling has involved simply assuming that alkali elements such as potassium occur as gas-phase atoms and are collisionally ionized following the Saha equation. We show that the Saha equation often does not hold, because free charges are produced by thermionic and ion emission and destroyed when they stick to grain surfaces. This means the ionization state depends not on the first ionization potential of the alkali atoms, but rather on the grains' work functions. The charged species' abundances typically rise abruptly above about 800 K, with little qualitative dependence on the work function, gas density, or dust-to-gas mass ratio. Applying our results, we find that protoplanetary disks' dead zone, where high diffusivities stifle magnetorotational turbulence, has its inner edge located where the temperature exceeds a threshold value ~1000 K. The threshold is set by ambipolar diffusion except at the highest densities, where it is set by Ohmic resistivity. We find that the disk gas can be diffusively loaded onto the stellar magnetosphere at temperatures below a similar threshold. We investigate whether the "short-circuit" instability of current sheets can operate in disks and find that it cannot, or works only in a narrow range of conditions; it appears not to be the chondrule formation mechanism. We also suggest that thermionic emission is important for determining the rate of Ohmic heating in hot Jupiters., Comment: 68 pages, 20 figures

AA Tau, a classical T Tauri star in the Taurus cloud, has been the subject of intensive photometric monitoring for more than two decades due to its quasi-cyclic variation in optical brightness. Beginning in 2011, AA Tau showed another peculiar variation -- its median optical though near-IR flux dimmed significantly, a drop consistent with a 4-mag increase in visual extinction. It has stayed in the faint state since.Here we present 4.7um CO rovibrational spectra of AA Tau over eight epochs, covering an eleven-year time span, that reveal enhanced 12CO and 13CO absorption features in the $J_{\rm low}\leqslant$13 transitions after the dimming. These newly appeared absorptions require molecular gas along the line of sight with T~500 K and a column density of log (N12CO)~18.5 cm^{-2}, with line centers that show a constant 6 km s$^{-1}$ redshift. The properties of the molecular gas confirm an origin in the circumstellar material. We suggest that the dimming and absorption are caused by gas and dust lifted to large heights by a magnetic buoyancy instability. This material is now propagating inward, and on reaching the star within a few years will be observed as an accretion outburst., Comment: Accepted for publication in ApJ

The classical T Tauri star AA Tau presented photometric variability attributed to an inner disk warp, caused by the interaction between the inner disk and an inclined magnetosphere. Previous studies of NGC 2264 have shown that similar photometric behavior is common among CTTS. The goal of this work is to investigate the main causes of the observed photometric variability of CTTS in NGC 2264 that present AA Tau-like light curves, and verify if an inner disk warp could be responsible for their variability. We investigate veiling variability in their spectra and u-r color variations and estimate parameters of the inner disk warp using an occultation model proposed for AA Tau. We compare infrared and optical light curves to analyze the dust responsible for the occultations. AA Tau-like variability is transient on a timescale of a few years. We ascribe it to stable accretion regimes and aperiodic variability to unstable accretion regimes and show that a transition, and even coexistence, between the two is common. We find evidence of hot spots associated with occultations, indicating that the occulting structures could be located at the base of accretion columns. We find average values of warp maximum height of 0.23 times its radial location, consistent with AA Tau, with variations of on average 11% between rotation cycles. We show that extinction laws in the inner disk indicate the presence of grains larger than interstellar grains. The inner disk warp scenario is consistent with observations for all but one periodic star in our sample. AA Tau-like systems comprise 14% of CTTS observed in NGC 2264, though this increases to 35% among systems of mass 0.7M_sun

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3

Protoplanetary disks are believed to accrete onto their central T Tauri star because of magnetic stresses. Recently published shearing box simulations indicate that Ohmic resistivity, ambipolar diffusion and the Hall effect all play important roles in disk evolution. In the presence of a vertical magnetic field, the disk remains laminar between 1-5au, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. Questions remain, however, about the establishment of a true physical wind solution in the shearing box simulations because of the symmetries inherent in the local approximation. We present global MHD simulations of protoplanetary disks that include Ohmic resistivity and ambipolar diffusion, where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry. Our results show that the disk remains laminar, and that a physical wind solution arises naturally in global disk models. The wind is sufficiently efficient to explain the observed accretion rates. Furthermore, the ionization fraction at intermediate disk heights is large enough for magneto-rotational channel modes to grow and subsequently develop into belts of horizontal field. Depending on the ionization fraction, these can remain quasi-global, or break-up into discrete islands of coherent field polarity. The disk models we present here show a dramatic departure from our earlier models including Ohmic resistivity only. It will be important to examine how the Hall effect modifies the evolution, and to explore the influence this has on the observational appearance of such systems, and on planet formation and migration., Comment: 18 pages, 12 figures, accepted for publication in ApJ

We carry out radiation hydrodynamical simulations of the formation of massive stars in the super-Eddington regime including both their radiative feedback and protostellar outflows. The calculations start from a prestellar core of dusty gas and continue until the star stops growing. The accretion ends when the remnants of the core are ejected, mostly by the force of the direct stellar radiation in the polar direction and elsewhere by the reradiated thermal infrared radiation. How long the accretion persists depends on whether the protostellar outflows are present. We set the mass outflow rate to 1% of the stellar sink particle's accretion rate. The outflows open a bipolar cavity extending to the core's outer edge, through which the thermal radiation readily escapes. The radiative flux is funneled into the polar directions while the core's collapse proceeds near the equator. The outflow thus extends the "flashlight effect", or anisotropic radiation field, found in previous studies from the few hundred AU scale of the circumstellar disk up to the 0.1 parsec scale of the core. The core's flashlight effect allows core gas to accrete on the disk for longer, in the same way that the disk's flashlight effect allows disk gas to accrete on the star for longer. Thus although the protostellar outflows remove material near the core's poles, causing slower stellar growth over the first few free-fall times, they also enable accretion to go on longer in our calculations. The outflows ultimately lead to stars of somewhat higher mass., Comment: accepted for publication at ApJ

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts., Comment: Accepted for publication in AJ. 39 pages; 6 tables; 25 figures, many of which are highly degraded to meet size limits. Please download the regular resolution version at http://web.ipac.caltech.edu/staff/amc/staufferetal2014.pdf

We show that the first order (non co-orbital) corotation torques are significantly modified by entropy gradients in a non-barotropic protoplanetary disk. Such non-barotropic torques can dramatically alter the balance that, for barotropic cases, results in the net eccentricity damping for giant gap-clearing planets embedded in the disk. We demonstrate that stellar illumination can heat the gap enough for the planet's orbital eccentricity to instead be excited. We also discuss the "Eccentricity Valley" noted in the known exoplanet population, where low-metallicity stars have a deficit of eccentric planets between $\sim 0.1$ and $\sim 1$ AU compared to metal-rich systems (Dawson & Murray-Clay 2013). We show that this feature in the planet distribution may be due to the self-shadowing of the disk by a rim located at the dust sublimation radius $\sim 0.1$ AU, which is known to exist for several T Tauri systems. In the shadowed region between $\sim 0.1$ and $\sim 1$ AU lack of gap insolation allows disk interactions to damp eccentricity. Outside such shadowed regions stellar illumination can heat the planetary gaps and drive eccentricity growth for giant planets. We suggest that the self-shadowing does not arise at higher metallicity due to the increased optical depth of the gas interior to the dust sublimation radius., Comment: 9 pages, 5 Figures, accepted by ApJ

We present global hydrodynamic and magnetohydrodynamic (MHD) simulations with mesh refinement of accreting planets embedded in protoplanetary disks (PPDs). The magnetized disk includes Ohmic resistivity that depends on the overlying mass column, leading to turbulent surface layers and a dead zone near the midplane. The main results are: (i) The accretion flow in the Hill sphere is intrinsically 3D for hydrodynamic and MHD models. Net inflow toward the planet is dominated by high latitude flows. A circumplanetary disk (CPD) forms. Its midplane flows outward in a pattern whose details differ between models. (ii) Gap opening magnetically couples and ignites the dead zone near the planet, leading to stochastic accretion, a quasi-turbulent flow in the Hill sphere and a CPD whose structure displays high levels of variability. (iii) Advection of magnetized gas onto the rotating CPD generates helical fields that launch magnetocentrifugally driven outflows. During one specific epoch a highly collimated, one-sided jet is observed. (iv) The CPD's surface density $\sim30{\rm\,g\,cm^{-2}}$, small enough for significant ionization and turbulence to develop. (v) The accretion rate onto the planet in the MHD simulation reaches a steady value $8 \times 10^{-3} {\rm M_\oplus yr^{-1}}$, and is similar in the viscous hydrodynamic runs. Our results suggest that gas accretion onto a forming giant planet within a magnetized PPD with dead zone allows rapid growth from Saturnian to Jovian masses. As well as being relevant for giant planet formation, these results have important implications for the formation of regular satellites around gas giant planets., Comment: 22 pages, 20 figures, second minor revision, accepted for publication in ApJ

We examine the observational appearance of partial gaps being opened by planets in protoplanetary disks, considering the effects of the inclination relative to the line of sight. The gap's trough is darkened by both shadowing and cooling, relative to the uninterrupted disk. The gap's outer wall is brightened by direct illumination and also by heating, which puffs it up so that it intercepts more starlight. In this paper, we examine the effects of inclination on resolved images of disks with and without gaps at a wide range of wavelengths. The scattering surface's offset from the disk midplane creates a brightness asymmetry along the axis of inclination, making the disk's near side appear brighter than the far side in scattered light. Finite disk thickness also causes the projected distances of equidistant points on the disk surface to be smaller on the near side of the disk as compared to the far side. Consequently, the gap shoulder on the near side of the disk should appear brighter and closer to the star than on the far side. However, if the angular resolution of the observation is coarser than the width of the brightened gap shoulder, then the gap shoulder on the far side may appear brighter because of its larger apparent size. We present a formula to recover the scale height and inclination angle of an imaged disk using simple geometric arguments and measuring disk asymmetries. Resolved images of circumstellar disks have revealed clearings and gaps, such as the transitional disk in LkCa 15. Models created using our synthetic imaging attempting to match the morphology of observed scattered light images of LkCa 15 indicate that the H-band flux deficit in the inner $\sim0.5\arcsec$ of the disk can be explained with a planet of mass greater than 0.5 Jupiter mass., Comment: 21 pages, 20 figures. Accepted to the Astrophysical Journal. Abstract is abridged. Full resolution images available at http://physics.uwyo.edu/~hannah/preprints/incldisk_hires.pdf

Magnetorotational instability (MRI) is the most promising mechanism behind accretion in low-mass protostellar disks. Here we present the first analysis of the global structure and evolution of non-ideal MRI-driven T-Tauri disks on million-year timescales. We accomplish this in a 1+1D simulation by calculating magnetic diffusivities and utilizing turbulence activity criteria to determine thermal structure and accretion rate without resorting to a 3-D magnetohydrodynamical (MHD) simulation. Our major findings are as follows. First, even for modest surface densities of just a few times the minimum-mass solar nebula, the dead zone encompasses the giant planet-forming region, preserving any compositional gradients. Second, the surface density of the active layer is nearly constant in time at roughly 10 g/cm2, which we use to derive a simple prescription for viscous heating in MRI-active disks for those who wish to avoid detailed MHD computations. Furthermore, unlike a standard disk with constant-alpha viscosity, the disk midplane does not cool off over time, though the surface cools as the star evolves along the Hayashi track. The ice line is firmly in the terrestrial planet-forming region throughout disk evolution and can move either inward or outward with time, depending on whether pileups form near the star. Finally, steady-state mass transport is a poor description of flow through an MRI-active disk. We caution that MRI activity is sensitive to many parameters, including stellar X-ray flux, grain size, gas/small grain mass ratio and magnetic field strength, and we have not performed an exhaustive parameter study here., Comment: Accepted for publication in Astrophysical Journal. 19 pages, including 8 figures

The edges of magnetically-dead zones in protostellar disks have been proposed as locations where density bumps may arise, trapping planetesimals and helping form planets. Magneto-rotational turbulence in magnetically-active zones provides both accretion of gas on the star and transport of mass to the dead zone. We investigate the location of the magnetically-active regions in a protostellar disk around a solar-type star, varying the disk temperature, surface density profile, and dust-to-gas ratio. We also consider stellar masses between 0.4 and 2 $M_\odot$, with corresponding adjustments in the disk mass and temperature. The dead zone's size and shape are found using the Elsasser number criterion with conductivities including the contributions from ions, electrons, and charged fractal dust aggregates. The charged species' abundances are found using the approach proposed by S. Okuzumi. The dead zone is in most cases defined by the ambipolar diffusion. In our maps, the dead zone takes a variety of shapes, including a fish-tail pointing away from the star and islands located on and off the midplane. The corresponding accretion rates vary with radius, indicating locations where the surface density will increase over time, and others where it will decrease. We show that density bumps do not readily grow near the dead zone's outer edge, independently of the disk parameters and the dust properties. Instead, the accretion rate peaks at the radius where the gas-phase metals freeze out. This could lead to clearing a valley in the surface density, and to a trap for pebbles located just outside the metal freeze-out line., Comment: 58 pages, 25 figures, 2 tables, accepted to ApJ

We investigate the viability of the magnetorotational instability (MRI) in accretion disks around both solar-type stars and very low mass stars. In particular, we determine the disk regions where the MRI can be shut off either by Ohmic resistivity (the so-called Dead and Undead Zones) or by ampipolar diffusion (a region we term the Zombie Zone). We consider 2 stellar masses: Mstar = 0.7 and 0.1 Msun. In each case, we assume that: the disk surface density profile is that of a scaled Minimum Mass Solar Nebula, with Mdisk/Mstar ~ 0.01 as currently estimated; disk ionisation is driven primarily by stellar X-rays, complemented by cosmic rays and radionuclides; and the stellar X-ray luminosity scales with bolometric luminosity as Lx/Lstar ~ 10^-3.5, as observed. Ionization rates are calculated with the MOCCASIN code, and ionisation balance determined using a simplified chemical network, including well-mixed 0.1 um grains at various levels of depletion. We find that (1) ambipolar diffusion is the primary factor controlling MRI activity in disks around both solar-type and very low mass stars. Assuming that the MRI yields the maximum possible field strength at each radius, we further find that: (2) the MRI-active layer constitutes only ~ 5-10% of the total disk mass; (3) the accretion rate (Mdot) varies radially in both magnitude and sign (inward or outward), implying time-variable accretion as well as the creation of disk gaps and overdensities, with consequences for planet formation and migration; (4) achieving the empirical accretion rates in solar-type and very low mass stars requires a depletion of well-mixed small grains by a factor of 10-1000 relative to the standard dust-to-gas mass ratio of 10^-2; and (5) the current non-detection of polarized emission from field-aligned grains in the outer disk regions is consistent with active MRI at those radii., Comment: 55 pages, 16 figures. ApJ, accepted

We examine the observational consequences of partial gaps being opened by planets in protoplanetary disks. We model the disk using a static alpha-disk model with detailed radiative transfer, parametrizing the shape and size of the partially cleared gaps based on the results of hydrodynamic simulations. Shadowing and illumination by stellar irradiation at the surface of the gap leads to increased contrast as the gap trough is deepened by shadowing and cooling and the far gap wall is puffed up by illumination and heating. In calculating observables, we find that multiple scattering is important and derive an approximation to include these effects. A gap produced by a 200 M_Earth (70 M_Earth) planet at 10 AU can lower/raise the midplane temperature of the disk by up to ~-25/+29% (~-11/+19%) by shadowing in the gap trough and illumination on the far shoulder of the gap. At the distance of Taurus, this gap would be resolvable with ~0.01" angular resolution. The gap contrast is most significant in scattered light and at thermal continuum wavelengths characteristic of the surface temperature, reducing or raising the surface brightness by up to order of magnitude. Since gaps sizes are correlated to planet mass, this is a promising way of finding and determining the masses of planets embedded in protoplanetary disks., Comment: 11 pages, 9 figures. Accepted to ApJ

(Abridged) Planetesimals embedded in a protoplanetary disc are stirred by gravitational torques exerted by density fluctuations in the surrounding turbulence. In particular, planetesimals in a disc supporting fully developed magneto-rotational turbulence are readily excited to velocity dispersions above the threshold for catastrophic disruption, halting planet formation. We aim to examine the stirring of planetesimals lying instead in a magnetically-decoupled midplane dead zone, stirred only by spiral density waves propagating out of the disc's magnetically-coupled turbulent surface layers. We extend previous studies to include a wider range of disc models, and explore the effects of varying the disc column density and external magnetic field strength. [...] The strength of the stirring is found to be independent of the gas surface density, which is contrary to the increase with disc mass expected from a simple linear wave picture. The discrepancy arises from the shearing out of density waves as they propagate into the dead zone, resulting in density structures near the midplane that exert weaker stochastic torques on average. We provide a simple analytic fit to our numerically obtained torque amplitudes that accounts for this effect. The stirring on the other hand depends sensitively on the net vertical magnetic flux, up to a saturation level above which magnetic forces dominate in the turbulent layers. For the majority of our models, the equilibrium planetesimal velocity dispersions lie between the thresholds for disrupting strong and weak aggregates, suggesting that collision outcomes will depend on material properties. However, discs with relatively weak magnetic fields yield reduced stirring, and their dead zones provide safe-havens even for the weakest planetesimals against collisional destruction., Comment: 20 pages, 18 figures, 3 tables, to be published in MNRAS

We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core., Comment: 20 pages, 15 figures, accepted for publication in ApJ

(abridged) Accretion in protoplanetary discs is thought to be driven by [...] turbulence via the magnetorotational instability (MRI). Recent work has shown that a planetesimal swarm embedded in a fully turbulent disc is subject to strong excitation of the velocity dispersion, leading to collisional destruction of bodies with radii R_p < 100 km. Significant diffusion of planetesimal semimajor axes also arises, leading to large-scale spreading of the planetesimal population throughout the inner regions of the protoplanetary disc, in apparent contradiction of constraints provided by the distribution of asteroids within the asteroid belt. In this paper, we examine the dynamics of planetesimals embedded in vertically stratified turbulent discs, with and without dead zones. Our main aims are to examine the turbulent excitation of the velocity dispersion, and the radial diffusion, of planetesimals in these discs. We employ three dimensional MHD simulations [...], along with an equilibrium chemistry model [...] We find that planetesimals in fully turbulent discs develop large random velocities that will lead to collisional destruction/erosion for bodies with sizes below 100 km, and undergo radial diffusion on a scale \sim 2.5 au over a 5 Myr disc life time. But planetesimals in a dead zone experience a much reduced excitation of their random velocities, and equilibrium velocity dispersions lie between the disruption thresholds for weak and strong aggregates for sizes R_p < 100 km. We also find that radial diffusion occurs over a much reduced length scale \sim 0.25 au over the disc life time, this being consistent with solar system constraints. We conclude that planetesimal growth via mutual collisions between smaller bodies cannot occur in a fully turbulent disc. By contrast, a dead zone may provide a safe haven in which km-sized planetesimals can avoid mutual destruction through collisions., Comment: 18 pages, 13 figures, 3 tables, MNRAS in press, minor corrections to match the published version

Variability is a defining characteristic of young stellar systems, and optical variability has been heavily studied to select and characterize the photospheric properties of young stars. In recent years, multi-epoch observations sampling a wider range of wavelengths and time-scales have revealed a wealth of time-variable phenomena at work during the star formation process. This splinter session was convened to summarize recent progress in providing improved coverage and understanding of time-variable processes in young stars and circumstellar disks. We begin by summarizing results from several multi-epoch Spitzer campaigns, which have demonstrated that many young stellar objects evidence significant mid-IR variability. While some of these variations can be attributed to processes in the stellar photosphere, others appear to trace short time-scale changes in the circumstellar disk which can be successfully modeled with axisymmetric or non-axisymmetric structures. We also review recent studies probing variability at shorter wavelengths that provide evidence for high frequency pulsations associated with accretion outbursts, correlated optical/X-ray variability in Classical T Tauri stars, and magnetic reversals in young solar analogs., Comment: 4 figures, 13 pages in ASP format; to be published in the proceedings of the Cool Stars 16 Workshop

The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treating the turbulence driven by the magneto-rotational instability. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the `butterfly pattern' seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible., Comment: 17 pages, 30 *.ps files for figures. Accepted 16 November 2009 in A&A