Source author record

Zhaohuan Zhu

Zhaohuan Zhu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR astro-ph.GA astro-ph.IM astro-ph astro-ph.HE Machine Learning physics.flu-dyn

Catalog footprint

What is connected

78works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Brightness variability in polar circumbinary disks

In binary systems with a strongly misaligned disk, the central binary stars can travel a significant vertical distance above and below the disk's orbital plane. This can cause large changes in illumination of the disk over the course of the binary orbital period. We use both analytic and radiative transfer models to examine the effect of changes in stellar illumination on the appearance of the disk, particularly in the case of the polar disk HD 98800B. We find that the observed flux from the disk can vary significantly over the binary orbital period, producing a periodically varying lightcurve which peaks twice each binary orbit. The amount of flux variation is strongly influenced by the disk geometry. We suggest that these flux variations produce several observable signatures, and that these observables may provide constraints on different properties of the disk such as its vertical structure, geometry, and cooling rate.

preprint2024arXiv

Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning

We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics of particles in the Epstein drag regime within a periodic domain of isotropic forced hydrodynamic turbulence. This setup is an idealized model relevant to the collisional growth of micron to mm-sized dust particles in early stage planet formation. The simulation data are used to train a U-Net deep learning model to predict gridded three-dimensional representations of the particle density and velocity fields, given as input the corresponding fluid fields. The trained model qualitatively captures the filamentary structure of clustered particles in a highly non-linear regime. We assess model fidelity by calculating metrics of the density field (the radial distribution function) and of the velocity field (the relative velocity and the relative radial velocity between particles). Although trained only on the spatial fields, the model predicts these statistical quantities with errors that are typically <10%. Our results suggest that, given appropriately expanded training data, deep learning could complement direct numerical simulations in predicting particle clustering within turbulent flows.

preprint2023arXiv

Structured Distributions of Gas and Solids in Protoplanetary Disks

Recent spatially-resolved observations of protoplanetary disks revealed a plethora of substructures, including concentric rings and gaps, inner cavities, misalignments, spiral arms, and azimuthal asymmetries. This is the major breakthrough in studies of protoplanetary disks since Protostars and Planets VI and is reshaping the field of planet formation. However, while the capability of imaging substructures in protoplanetary disks has been steadily improving, the origin of many substructures are still largely debated. The structured distributions of gas and solids in protoplanetary disks likely reflect the outcome of physical processes at work, including the formation of planets. Yet, the diverse properties among the observed protoplanetary disk population, for example, the number and radial location of rings and gaps in the dust distribution, suggest that the controlling process may differ between disks and/or the outcome may be sensitive to stellar or disk properties. In this review, we (1) summarize the existing observations of protoplanetary disk substructures collected from the literature; (2) provide a comprehensive theoretical review of various processes proposed to explain observed protoplanetary disk substructures; (3) compare current theoretical predictions with existing observations and highlight future research directions to distinguish between different origins; and (4) discuss implications of state-of-the-art protoplanetary disk observations to protoplanetary disk and planet formation theory.

preprint2022arXiv

ALMA Discovery of a Disk around the Planetary-mass Companion SR 12 c

We report an Atacama Large Millimeter/submillimeter Array 0.88 mm (Band 7) continuum detection of the accretion disk around SR 12 c, an $\sim$11 $M_{\rm Jup}$ planetary-mass companion (PMC) orbiting its host binary at 980 au. This is the first submillimeter detection of a circumplanetary disk around a wide PMC. The disk has a flux density of $127 \pm14~μ$Jy and is not resolved by the $\sim$0.1" beam, so the dust disk radius is likely less than 5 au and can be much smaller if the dust continuum is optically thick. If, however, the dust emission is optically thin, then the SR 12 c disk has a comparable dust mass to the circumplanetary disk around PDS 70 c but is about five times lower than that of the $\sim$12 $M_{\rm Jup}$ free-floating OTS 44. This suggests that disks around bound and unbound planetary-mass objects can span a wide range of masses. The gas mass estimated with an accretion rate of $10^{-11}~M_\odot$ yr$^{-1}$ implies a gas-to-dust ratio higher than 100. If cloud absorption is not significant, a nondetection of ${}^{12}$CO(3-2) implies a compact gas disk around SR 12 c. Future sensitive observations may detect more PMC disks at 0.88 mm flux densities of $\lesssim$100 $μ$Jy.

preprint2022arXiv

Direct Formation of Planetary Embryos in Self-Gravitating Disks

Giant planets have been discovered at large separations from the central star. Moreover, a striking number of young circumstellar disks have gas and/or dust gaps at large orbital separations, potentially driven by embedded planetary objects. To form massive planets at large orbital separations through core accretion within disk lifetime, however, an early solid body to seed pebble and gas accretion is desirable. Young protoplanetary disks are likely self-gravitating, and these gravitoturbulent disks may efficiently concentrate solid material at the midplane driven by spiral waves. We run 3D local hydrodynamical simulations of gravitoturbulent disks with Lagrangian dust particles to determine whether particle and gas self-gravity can lead to the formation of dense solid bodies, seeding later planet formation. When self-gravity between dust particles is included, solids of size $\mathrm{St} = 0.1$ to $1$ concentrate within the gravitoturbulent spiral features and collapse under their own self-gravity into dense clumps up to several $M_{\oplus}$ in mass at wide orbits. Simulations with dust that drift most efficiently, $\mathrm{St}=1$, form the most massive clouds of particles, while simulations with smaller dust particles, $\mathrm{St}=0.1$, have clumps with masses an order of magnitude lower. When the effect of dust backreaction onto the gas is included, dust clumps become smaller by a factor of a few but more numerous. The existence of large solid bodies at an early stage of the disk can accelerate the planet formation process, particularly at wide orbital separations, and potentially explain planets distant from the central stars and young protoplanetary disks with substructures.

preprint2022arXiv

Formation of Dust Rings and Gaps in Non-ideal MHD Disks Through Meridional Gas Flows

Rings and gaps are commonly observed in the dust continuum emission of young stellar disks. Previous studies have shown that substructures naturally develop in the weakly ionized gas of magnetized, non-ideal MHD disks. The gas rings are expected to trap large mm/cm-sized grains through pressure gradient-induced radial dust-gas drift. Using 2D (axisymmetric) MHD simulations that include ambipolar diffusion and dust grains of three representative sizes (1~mm, 3.3~mm, and 1~cm), we show that the grains indeed tend to drift radially relative to the gas towards the centers of the gas rings, at speeds much higher than in a smooth disk because of steeper pressure gradients. However, their spatial distribution is primarily controlled by meridional gas motions, which are typically much faster than the dust-gas drift. In particular, the grains that have settled near the midplane are carried rapidly inwards by a fast accretion stream to the inner edges of the gas rings, where they are lifted up by the gas flows diverted away from the midplane by a strong poloidal magnetic field. The flow pattern in our simulation provides an attractive explanation for the meridional flows recently inferred in HD 163296 and other disks, including both "collapsing" regions where the gas near the disk surface converges towards the midplane and a disk wind. Our study highlights the prevalence of the potentially observable meridional flows associated with the gas substructure formation in non-ideal MHD disks and their crucial role in generating rings and gaps in dust.

preprint2022arXiv

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team.

preprint2022arXiv

HST/WFC3 H$α$ Direct-Imaging Detection of a Point-like Source in the Disk Cavity of AB Aur

Accreting protoplanets enable the direct characterization of planet formation. As part of a high-contrast imaging search for accreting planets with the Hubble Space Telescope (HST) Wide Field Camera 3, we present H$α$ images of AB Aurigae (AB Aur), a Herbig Ae/Be star harboring a transition disk. The data were collected in two epochs of direct-imaging observations using the F656N narrow-band filter. After subtracting the point spread function of the primary star, we identify a point-like source located at a P.A. of $182.5^{\circ}\pm1.4^{\circ}$ and a separation of $600\pm22$~mas relative to the host star. The position is consistent with the recently identified protoplanet candidate AB Aur b. The source is visible in two individual epochs separated by ${\sim}50$ days and the H$α$ intensities in the two epochs agree. The H$α$ flux density is $F_ν=1.5\pm0.4$~mJy, $3.2\pm0.9$ times of the optical continuum determined by published HST/STIS photometry. In comparison to PDS 70 b and c, the H$α$ excess emission is weak. The central star is accreting and the stellar H$α$ emission has a similar line-to-continuum ratio as seen in AB Aur b. We conclude that both planetary accretion and scattered stellar light are possible sources of the H$α$ emission, and the H$α$ detection alone does not validate AB Aur b as an accreting protoplanet. Disentangling the origin of the emission will be crucial for probing planet formation in the AB Aur disk.

preprint2022arXiv

Investigating the Future Potential of an Upgraded ALMA to Image Planet Forming Disks at Sub-au Scales

In recent years, ALMA has been able to observe large-scale substructures within protoplanetary disks. Comparison with the predictions from models of planet-disk interaction has indicated that most of these disk substructures can be explained by the presence of planets with the mass of Neptune or larger at orbital radii of $\approx 5 - 100$ au. Better resolution is needed to observe structures closer to the star, where terrestrial planets are expected to form, as well as structures opened by planets with masses lower than Neptune. We investigate the capabilities of a possible extension to ALMA that would double the longest baseline lengths in the array to detect and resolve disk substructures opened by Earth-mass and Super Earth planets at orbital radii of $1-5$ au. By simulating observations of a family of disk models using this extended configuration in ALMA Bands 6 and 7, we show that an upgraded ALMA would detect gaps in disks formed by super-Earths as close as 1 au, as well as Earth-mass planets down to $2-3$ au from the young host stars in nearby star forming regions.

preprint2022arXiv

Maximum Temperatures in Evolving Protoplanetary Discs and Composition of Planetary Building Blocks

The maximum temperature and radial temperature profile in a protoplanetary disc are important for the condensation of different elements in the disc. We simulate the evolution of a set of protoplanetary discs from the collapse of their progenitor molecular cloud cores as well as the dust decoupling within the discs as they evolve. We show how the initial properties of the cloud cores affect the thermal history of the protoplanetary discs using a simple viscous disc model. Our results show that the maximum midplane temperature in the disc occurs within 0.5 AU. It increases with the initial cloud temperature and decreases with its angular velocity and the viscosity of the disc. From the observed properties of the molecular cloud cores we find the median value of the maximum temperature is around 1250 K, with roughly 90% of them being less than 1500 K - a value that is lower than the 50% condensation temperatures of most refractory elements. Therefore, only cloud cores with high initial temperatures or low angular velocities and/or low viscosities within the planet-forming discs will result in refractory-rich planetesimals. To reproduce the volatile depletion pattern of CM, CO, and CV chondrites and the terrestrial planets in Solar system, one must either have rare properties of the initial molecular cloud cores like high core temperature, or other sources of energy to heat the disc to sufficiently high temperatures. Alternatively, the volatile depletion observed in these chondrites may be inherited from the progenitor molecular cloud.

preprint2022arXiv

Thermal instability in radiation hydrodynamics: instability mechanisms, position-dependent S-curves, and attenuation curves

Local thermal instability can plausibly explain the formation of multiphase gas in many different astrophysical environments, but the theory is only well understood in the optically thin limit of the equations of radiation hydrodynamics (RHD). Here we lay groundwork for transitioning from this limit to a full RHD treatment assuming a gray opacity formalism. We consider a situation where the gas becomes thermally unstable due to the hardening of the radiation field when the main radiative processes are free-free cooling and Compton heating. We identify two ways in which this can happen: (i) when the Compton temperature increases with time, through a rise in either the intensity or energy of a hard X-ray component; and (ii) when attenuation reduces the flux of the thermal component so that the Compton temperature increases with depth through the slab. Both ways likely occur in the broad line region of active galactic nuclei where columns of gas can be ionization bounded. In such instances where attenuation is significant, thermal equilibrium solution curves become position-dependent and it no longer suffices to assess the stability of an irradiated column of gas at all depths using a single equilibrium curve. We demonstrate how to analyze a new equilibrium curve -- the attenuation curve -- for this purpose, and we show that by Field's instability criterion, a negative slope along this curve indicates that constant density slabs are thermally unstable whenever the gas temperature increases with depth.

preprint2021arXiv

Constraining Protoplanetary Disk Accretion and Young Planets Using ALMA Kinematic Observations

Recent ALMA molecular line observations have revealed 3-D gas velocity structure in protoplanetary disks, shedding light on mechanisms of disk accretion and structure formation. 1) By carrying out viscous simulations, we confirm that the disk's velocity structure differs dramatically using vertical stress profiles from different accretion mechanisms. Thus, kinematic observations tracing flows at different disk heights can potentially distinguish different accretion mechanisms. On the other hand, the disk surface density evolution is mostly determined by the vertically integrated stress. The sharp disk outer edge constrained by recent kinematic observations can be caused by a radially varying $α$ in the disk. 2) We also study kinematic signatures of a young planet by carrying out 3-D planet-disk simulations. The relationship between the planet mass and the "kink" velocity is derived, showing a linear relationship with little dependence on disk viscosity, but some dependence on disk height when the planet is massive, e.g. $10 M_J$. We predict the "kink" velocities for the potential planets in DSHARP disks. At the gap edge, the azimuthally-averaged velocities at different disk heights deviate from the Keplerian velocity at similar amplitudes, and its relationship with the planet mass is consistent with that in 2-D simulations. After removing the planet, the azimuthally-averaged velocity barely changes within the viscous timescale, and thus the azimuthally-averaged velocity structure at the gap edge is due to the gap itself and not directly caused to the planet. Combining both axisymmetric kinematic observations and the residual "kink" velocity is needed to probe young planets in protoplanetary disks.

preprint2021arXiv

Kozai-Lidov oscillations triggered by a tilt instability of detached circumplanetary discs

Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star-planet system. We use three-dimensional hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long term evolution. Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, of $H/r\gtrsim 0.05$, with a weak dependence on viscosity in the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs that satisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where the circumplanetary disc becomes unstable to Kozai--Lidov (KL) oscillations. Dissipation in the Kozai--Lidov unstable regime caps further tilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodic accretion events. We discuss implications of the joint tilt--KL instability for the detectability of circumplanetary discs, for the obliquity evolution of forming giant planets, and for the formation of satellite systems.

preprint2021arXiv

The ALMA Early Science view of FUor/EXor objects. I. Through the looking-glass of V2775 Ori

As part of an ALMA survey to study the origin of episodic accretion in young eruptive variables, we have observed the circumstellar environment of the star V2775 Ori. This object is a very young, pre-main sequence object which displays a large amplitude outburst characteristic of the FUor class. We present Cycle-2 band 6 observations of V2775 Ori with a continuum and CO (2-1) isotopologue resolution of 0.25\as (103 au). We report the detection of a marginally resolved circumstellar disc in the ALMA continuum with an integrated flux of $106 \pm 2$ mJy, characteristic radius of $\sim$ 30 au, inclination of $14.0^{+7.8}_{-14.5}$ deg, and is oriented nearly face-on with respect to the plane of the sky. The \co~emission is separated into distinct blue and red-shifted regions that appear to be rings or shells of expanding material from quasi-episodic outbursts. The system is oriented in such a way that the disc is seen through the outflow remnant of V2775 Ori, which has an axis along our line-of-sight. The $^{13}$CO emission displays similar structure to that of the \co, while the C$^{18}$O line emission is very weak. We calculated the expansion velocities of the low- and medium-density material with respect to the disc to be of -2.85 km s$^{-1}$ (blue), 4.4 km s$^{-1}$ (red) and -1.35 and 1.15 km s$^{-1}$ (for blue and red) and we derived the mass, momentum and kinetic energy of the expanding gas. The outflow has an hourglass shape where the cavities are not seen. We interpret the shapes that the gas traces as cavities excavated by an ancient outflow. We report a detection of line emission from the circumstellar disc and derive a lower limit of the gas mass of 3 \MJup.

preprint2020arXiv

A fast-growing tilt instability of detached circumplanetary disks

Accretion disks in binary systems can exhibit a tilt instability, arising from the interaction between components of the tidal potential and dissipation. Using a linear analysis, we show that the aspect ratios and outer radii of circumplanetary disks provide favorable conditions for tilt growth. We quantify the growth rate of the instability using particle-based ({\sc phantom}) and grid-based ({\sc athena++}) hydrodynamic simulations. For a disk with outer aspect ratio $H/r \simeq 0.1$, initially moderate tilts double on a time scale of about 15-30 binary orbits. Our results imply that detached circumplanetary disks, whose evolution is not entirely controlled by accretion from the circumstellar disk, may commonly be misaligned to the planetary orbital plane. We discuss implications for planetary spin evolution, and possible interactions between the tilt instability and Kozai-Lidov dynamics.

preprint2020arXiv

A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $α$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."

preprint2020arXiv

A Study of Millimeter Variability in FUor Objects

FU Orionis objects (FUors) are rapidly-accreting, pre-main sequence objects that are known to exhibit large outbursts at optical and near-infrared wavelengths, with post-eruption, small-scale photometric variability superimposed on longer-term trends. In contrast, little is known about the variability of FUors at longer wavelengths. To explore this further, we observed six FUor objects using the NOrthern Extended Millimeter Array (NOEMA) and for a subset of three objects we obtained coordinated observations with NOEMA and the Lowell Discovery Telescope (LDT). In combination with previously published NOEMA observations from 2014, our 2017 observations of V1735 Cyg provide the first detection of variability in an FUor object at 2.7 mm. In the absence of significant optical variability, we discount the possibility that the mm flux density changed as a result of irradiation from the central disk. In addition, a change in the dust mass due to infall is highly unlikely. A plausible explanation for the change in 2.7 mm flux density is variability in free-free emission due to changes in the object's jet/wind. Thus, it may be that free-free emission in some FUor objects is significant at $\sim$3 mm and must be considered when deriving disk masses in order to help constrain the mechanism responsible for triggering FUor outbursts.

preprint2020arXiv

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits.

preprint2020arXiv

ALMA 0.88 mm Survey of Disks around Planetary-mass Companions

Characterizing the physical properties and compositions of circumplanetary disks can provide important insights into the formation of giant planets and satellites. We report ALMA 0.88 mm (Band 7) continuum observations of six planetary-mass (10-20 $M_{\rm Jup}$) companions: CT Cha b, 1RXS 1609 b, ROXs 12 b, ROXs 42B b, DH Tau b, and FU Tau b. No continuum sources are detected at the locations of the companions down to 3$σ$ limits of 120-210 $μ$Jy. Given these non-detections, it is not clear whether disks around planetary-mass companions indeed follow the disk flux-host mass trend in the stellar regime. The faint radio brightness of these companion disks may result from a combination of fast radial drift and a lack of dust traps. Alternatively, as disks in binary systems are known to have significantly lower millimeter fluxes due to tidal interactions, these companion disks may instead follow the relationship of moderate-separation binary stars. This scenario can be tested with sensitive continuum imaging at rms levels of $\lesssim$10 $μ$Jy.

preprint2020arXiv

Dust Condensation in Evolving Discs and the Composition of Planetary Building Blocks

Partial condensation of dust from the Solar nebula is likely responsible for the diverse chemical compositions of chondrites and rocky planets/planetesimals in the inner Solar system. We present a forward physical-chemical model of a protoplanetary disc to predict the chemical compositions of planetary building blocks that may form from such a disc. Our model includes the physical evolution of the disc and the condensation, partial advection, and decoupling of the dust within it. The chemical composition of the condensate changes with time and radius. We compare the results of two dust condensation models: one where an element condenses when the midplane temperature in the disc is lower than the 50\% condensation temperature ($\rm T_{50}$) of that element and the other where the condensation of the dust is calculated by a Gibbs free energy minimization technique assuming chemical equilibrium at local disc temperature and pressure. The results of two models are generally consistent with some systematic differences of $\sim 10$\% depending upon the radial distance and an element's condensation temperature. Both models predict compositions similar to CM, CO, and CV chondrites provided that the decoupling timescale of the dust is on the order of the evolution timescale of the disc or longer. If the decoupling timescale is too short, the composition deviates significantly from the measured values. These models may contribute to our understanding of the chemical compositions of chondrites, and ultimately the terrestrial planets in the solar system, and may constrain the potential chemical compositions of rocky exoplanets.

preprint2020arXiv

Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori's Accretion Disk and Observational Signatures of Magnetic Fields

FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori's accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori's inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z$\sim$R, similar to the "surface accretion" mechanism in previous locally-isothermal MHD simulations. (2) A moderate disk wind is launched in the vertical field simulations with a terminal speed of $\sim$300-500 km/s and a mass loss rate of 1-10\% the disk accretion rate, which is consistent with observations. Disk wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disk photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disk wind, the disk photosphere temperature is lower than that predicted by the thin disk theory, and the previously inferred disk accretion rate may be lower than the real accretion rate by a factor of $\sim$2-3.

preprint2020arXiv

Large-scale CO spiral arms and complex kinematics associated with the T Tauri star RU Lup

While protoplanetary disks often appear to be compact and well-organized in millimeter continuum emission, CO spectral line observations are increasingly revealing complex behavior at large distances from the host star. We present deep ALMA maps of the $J=2-1$ transition of $^{12}$CO, $^{13}$CO, and C$^{18}$O, as well as the $J=3-2$ transition of DCO$^+$, toward the T Tauri star RU Lup at a resolution of $\sim0.3''$ ($\sim50$ au). The CO isotopologue emission traces four major components of the RU Lup system: a compact Keplerian disk with a radius of $\sim120$ au, a non-Keplerian ``envelope-like'' structure surrounding the disk and extending to $\sim260$ au from the star, at least five blueshifted spiral arms stretching up to 1000 au, and clumps outside the spiral arms located up to 1500 au in projection from RU Lup. We comment on potential explanations for RU Lup's peculiar gas morphology, including gravitational instability, accretion of material onto the disk, or perturbation by another star. RU Lup's extended non-Keplerian CO emission, elevated stellar accretion rate, and unusual photometric variability suggest that it could be a scaled-down Class II analog of the outbursting FU Ori systems.

preprint2020arXiv

Near-Infrared Imaging of a Spiral in the CQ Tau Disk

We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.

preprint2020arXiv

The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures

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

preprint2020arXiv

The effects of disk self-gravity and radiative cooling on the formation of gaps and spirals by young planets

We have carried out two-dimensional hydrodynamical simulations to study the effects of disk self-gravity and radiative cooling on the formation of gaps and spirals. (1) With disk self-gravity included, we find stronger, more tightly-wound spirals and deeper gaps in more massive disks. The deeper gaps are due to the larger Angular Momentum Flux (AMF) of the waves excited in more massive disks, as expected from the linear theory. The position of the secondary gap does not change, provided that the disk is not extremely massive ($Q \gtrsim 2$). (2) With radiative cooling included, the excited spirals become monotonically more open (less tightly-wound) as the disk's cooling timescale increases. On the other hand, the amplitude and strength of the spirals decrease when the cooling time increases from a small value to $\sim 1/Ω$, but then the amplitude starts to increase again when the cooling time continues to increase. This indicates that radiative dissipation becomes important for waves with $T_{cool}\sim$ 1. Consequently, the induced primary gap is narrower and the secondary gap becomes significantly shallower when the cooling time becomes $\sim 1/Ω$. When the secondary gap is present, the position of it moves to the inner disk from the fast cooling cases to the slow cooling cases. The dependence of gap properties on the cooling timescale (e.g. in AS 209) provides a new way to constrain the disk optical depth and thus disk surface density.

preprint2020arXiv

The Planetary Luminosity Problem: "Missing Planets" and the Observational Consequences of Episodic Accretion

The high occurrence rates of spiral arms and large central clearings in protoplanetary disks, if interpreted as signposts of giant planets, indicate that gas giants form commonly as companions to young stars ($<$ few Myr) at orbital separations of 10--300\,au. However, attempts to directly image this giant planet population as companions to more mature stars ($> 10$\, Myr) have yielded few successes. This discrepancy could be explained if most giant planets form "cold start," i.e., by radiating away much of their formation energy as they assemble their mass, rendering them faint enough to elude detection at later times. In that case, giant planets should be bright at early times, during their accretion phase, and yet forming planets are detected only rarely through direct imaging techniques. Here we explore the possibility that the low detection rate of accreting planets is the result of episodic accretion through a circumplanetary disk. We also explore the possibility that the companion orbiting the Herbig Ae star HD~142527 may be a giant planet undergoing such an accretion outburst.

preprint2019arXiv

Circumplanetary Disk Dynamics in the Isothermal and Adiabatic Limits

Circumplanetary disks (CPDs) may be essential to the formation of planets, regulating their spin and accretion evolution. We perform a series of 3D hydrodynamics simulations in both the isothermal and adiabatic limits to systematically measure the rotation rates, sizes, and masses of CPDs as functions of $q_{\rm thermal}$, the ratio of the planet mass to the disk thermal mass. Our $q_{\rm thermal}$ ranges from 0.1 to 4; for our various disk temperatures, this corresponds to planet masses between 1 Earth mass and 4 Jupiter masses. Within this parameter space, we find that isothermal CPDs are disky and bound within $\sim$10\% of the planet's Bondi radius $r_{\rm B}$, with the innermost $\sim0.05\,r_{\rm B}$ in full rotational support. Adiabatic CPDs are spherical (therefore not actually "disks"), bound within $\sim0.2\,r_{\rm B}$, and mainly pressure-supported with rotation rates scaling linearly with $q_{\rm thermal}$; extrapolation suggests full rotational support of adiabatic envelopes at $\sim10\,q_{\rm thermal}$. Fast rotation and 3D supersonic flow render isothermal CPDs significantly different in structure from --- and orders of magnitude less massive than --- their 1D isothermal hydrostatic counterparts. Inside a minimum-mass solar nebula, even a maximally cooled, isothermal CPD around a 10 Earth-mass core may have less than 1 Earth mass, suggesting that gas giant formation may hinge on angular momentum transport processes in CPDs. Our CPD sizes and masses appear consistent with the regular satellites orbiting solar system giants.

preprint2019arXiv

Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

preprint2019arXiv

Morphological signatures induced by dust back reaction in discs with an embedded planet

Recent observations have revealed a gallery of substructures in the dust component of nearby protoplanetary discs, including rings, gaps, spiral arms, and lopsided concentrations. One interpretation of these substructures is the existence of embedded planets. Not until recently, however, most of the modelling effort to interpret these observations ignored the dust back reaction to the gas. In this work, we conduct local-shearing-sheet simulations for an isothermal, inviscid, non-self-gravitating, razor-thin dusty disc with a planet on a fixed circular orbit. We systematically examine the parameter space spanned by planet mass ($0.1M_\text{th} \le M_\text{p} \le 1M_\text{th}$, where $M_\text{th}$ is the thermal mass), dimensionless stopping time ($10^{-3} \le τ_\text{s} \le 1$), and solid abundance ($0 < Z \le 1$). We find that when the dust particles are tightly coupled to the gas ($τ_\text{s} < 0.1$), the spiral arms are less open and the gap driven by the planet becomes deeper with increasing $Z$, consistent with a reduced speed of sound in the approximation of a single dust-gas mixture. By contrast, when the dust particles are marginally coupled ($0.1 \lesssim τ_\text{s} \lesssim 1$), the spiral structure is insensitive to $Z$ and the gap structure in the gas can become significantly skewed and unidentifiable. When the latter occurs, the pressure maximum radially outside of the planet is weakened or even extinguished, and hence dust filtration by a low-mass ($M_\text{p} < M_\text{th}$) planet could be reduced or eliminated. Finally, we find that the gap edges where the dust particles are accumulated as well as the lopsided large-scale vortices driven by a massive planet, if any, are unstable, and they are broken into numerous small-scale dust-gas vortices.

preprint2019arXiv

Resolving the FU Ori System with ALMA: Interacting Twin Disks?

FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts of several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypical FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system with ALMA. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of $\sim$11 au. The ratio between major and minor axes shows that the inclination of the disk is $\sim$37 deg. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between north and south components is 0.6'', i.e. $\sim$250 au. Maps of $^{12}$CO emission show a complex kinematic environment with signatures disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models (Zhu et al.), yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M$_{\odot}$ and 3.8$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, respectively.

preprint2016arXiv

A Dwarf Protoplanetary Disk around XZ Tau B

We report the discovery of a dwarf protoplanetary disk around the star XZ Tau B that shows all the features of a classical transitional disk but on a much smaller scale. The disk has been imaged with the Atacama Large Millimeter/Submillimeter Array (ALMA), revealing that its dust emission has a quite small radius of ~ 3.4 au and presents a central cavity of ~ 1.3 au in radius that we attribute to clearing by a compact system of orbiting (proto)planets. Given the very small radii involved, evolution is expected to be much faster in this disk (observable changes in a few months) than in classical disks (observable changes requiring decades) and easy to monitor with observations in the near future. From our modeling we estimate that the mass of the disk is large enough to form a compact planetary system.

preprint2016arXiv

Architecture design study and technology roadmap for the Planet Formation Imager (PFI)

The Planet Formation Imager (PFI) Project has formed a Technical Working Group (TWG) to explore possible facility architectures to meet the primary PFI science goal of imaging planet formation in situ in nearby star- forming regions. The goals of being sensitive to dust emission on solar system scales and resolving the Hill-sphere around forming giant planets can best be accomplished through sub-milliarcsecond imaging in the thermal infrared. Exploiting the 8-13 micron atmospheric window, a ground-based long-baseline interferometer with approximately 20 apertures including 10km baselines will have the necessary resolution to image structure down 0.1 milliarcseconds (0.014 AU) for T Tauri disks in Taurus. Even with large telescopes, this array will not have the sensitivity to directly track fringes in the mid-infrared for our prime targets and a fringe tracking system will be necessary in the near-infrared. While a heterodyne architecture using modern mid-IR laser comb technology remains a competitive option (especially for the intriguing 24 and 40μm atmospheric windows), the prioritization of 3-5μm observations of CO/H2O vibrotational levels by the PFI-Science Working Group (SWG) pushes the TWG to require vacuum pipe beam transport with potentially cooled optics. We present here a preliminary study of simulated L- and N-band PFI observations of a realistic 4-planet disk simulation, finding 21x2.5m PFI can easily detect the accreting protoplanets in both L and N-band but can see non-accreting planets only in L band. (abridged -- see PDF for full abstract)

preprint2016arXiv

Dust dynamics in 2D gravito-turbulent disks

The dynamics of solid bodies in protoplanetary disks are subject to the properties of any underlying gas turbulence. Turbulence driven by disk self-gravity shows features distinct from those driven by the magnetorotational instability (MRI). We study the dynamics of solids in gravito-turbulent disks with two-dimensional (in the disk plane), hybrid (particle and gas) simulations. Gravito-turbulent disks can exhibit stronger gravitational stirring than MRI-active disks, resulting in greater radial diffusion and larger eccentricities and relative speeds for large particles (those with dimensionless stopping times $t_{stop} Ω> 1$, where $Ω$ is the orbital frequency). The agglomeration of large particles into planetesimals by pairwise collisions is therefore disfavored in gravito-turbulent disks. However, the relative speeds of intermediate-size particles $t_{stop} Ω\sim 1$ are significantly reduced as such particles are collected by gas drag and gas gravity into coherent filament-like structures with densities high enough to trigger gravitational collapse. First-generation planetesimals may form via gravitational instability of dust in marginally gravitationally unstable gas disks.

preprint2016arXiv

Exoplanet Orbital Eccentricities Derived From LAMOST-Kepler Analysis

The nearly circular (mean eccentricity <e>~0.06) and coplanar (mean mutual inclination <i>~3 deg) orbits of the Solar System planets motivated Kant and Laplace to put forth the hypothesis that planets are formed in disks, which has developed into the widely accepted theory of planet formation. Surprisingly, the first several hundred extrasolar planets (mostly Jovian) discovered using the Radial Velocity (RV) technique are commonly on eccentric orbits (<e> ~ 0.3). This raises a fundamental question: Are the Solar System and its formation special? The Kepler mission has found thousands of transiting planets dominated by sub-Neptunes, but most of their orbital eccentricities remain unknown. By using the precise spectroscopic host star parameters from the LAMOST observations, we measure the eccentricity distributions for a large (698) and homogeneous Kepler planet sample with transit duration statistics. Nearly half of the planets are in systems with single transiting planets (singles), while the other half are multiple-transiting planets (multiples). We find an eccentricity dichotomy: on average, Kepler singles are on eccentric orbits with <e>~0.3, while the multiples are on nearly circular (<e> = 0.04^{+0.03}_{-0.04}) and coplanar (<i> = 1.4^{+0.8}_{-1.1} deg) orbits similar to the Solar System planets. Our results are consistent with previous studies of smaller samples and individual systems. We also show that Kepler multiples and solar system objects follow a common relation <e>~(1-2)x<i> between mean eccentricities and mutual inclinations. The prevalence of circular orbits and the common relation may imply that the solar system is not so atypical in the galaxy after all.

preprint2016arXiv

Gas and dust hydrodynamical simulations of massive lopsided transition discs - I. Gas distribution

Motivated by lopsided structures observed in some massive transition discs, we have carried out 2D numerical simulations to study vortex structure in massive discs, including the effects of disc self-gravity and the indirect force which is due to the displacement of the central star from the barycenter of the system by the lopsided structure. When only the indirect force is included, we confirm the finding by Mittal & Chiang (2015) that the vortex becomes stronger and can be more than two pressure scale heights wide, as long as the disc-to-star mass ratio is >1%. Such wide vortices can excite strong density waves in the disc and therefore migrate inwards rapidly. However, when disc self-gravity is also considered in simulations, self-gravity plays a more prominent role on the vortex structure. We confirm that when the disc Toomre Q parameter is smaller than pi/(2h), where h is the disc's aspect ratio, the vortices are significantly weakened and their inward migration slows down dramatically. Most importantly, when the disc is massive enough (e.g. Q~3), we find that the lopsided gas structure orbits around the star at a speed significantly slower than the local Keplerian speed. This sub-Keplerian pattern speed can lead to the concentration of dust particles at a radius beyond the lopsided gas structure (as shown in Paper II). Overall, disc self-gravity regulates the vortex structure in massive discs and the radial shift between the gas and dust distributions in vortices within massive discs may be probed by future observations.

preprint2016arXiv

Global MHD Simulations of Accretion Disks in Cataclysmic Variables (CVs): I. The Importance of Spiral Shocks

We present results from the first global 3D MHD simulations of accretion disks in Cataclysmic Variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared to that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component to the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states and both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. 1) Linear wave dispersion relation fits the pitch angles of spiral density waves very well. 2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. 3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate alpha_{eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. 4) Using the mass accretion rate measured in our simulations to directly measure alpha defined in standard thin-disk theory, we find 0.02 < alpha_{eff} < 0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae (DNe). In this regime the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime, and therefore future models with more realistic thermodynamics and non-ideal MHD are warranted.

preprint2016arXiv

Imaging the water-snow line during a protostellar outburst

A snow-line is the region of a protoplanetary disk at which a major volatile, such as water or carbon monoxide, reaches its condensation temperature. Snow-lines play a crucial role in disk evolution by promoting the rapid growth of ice-covered grains. Signatures of the carbon monoxide snow-line (at temperatures of around 20 kelvin) have recently been imaged in the disks surrounding the pre-main-sequence stars TW Hydra and HD163296, at distances of about 30 astronomical units (au) from the star. But the water snow-line of a protoplanetary disk (at temperatures of more than 100 kelvin) has not hitherto been seen, as it generally lies very close to the star (less than 5 au away for solar-type stars). Water-ice is important because it regulates the efficiency of dust and planetesimal coagulation, and the formation of comets, ice giants and the cores of gas giants. Here we report ALMA images at 0.03-arcsec resolution (12 au) of the protoplanetary disk around V883 Ori, a protostar of 1.3 solar masses that is undergoing an outburst in luminosity arising from a temporary increase in the accretion rate. We find an intensity break corresponding to an abrupt change in the optical depth at about 42 au, where the elevated disk temperature approaches the condensation point of water, from which we conclude that the outburst has moved the water snow-line. The spectral behaviour across the snow-line confirms recent model predictions: dust fragmentation and the inhibition of grain growth at higher temperatures results in soaring grain number densities and optical depths. As most planetary systems are expected to experience outbursts caused by accretion during their formation our results imply that highly dynamical water snow-lines must be considered when developing models of disk evolution and planet formation.

preprint2016arXiv

Planet Formation Imager (PFI): science vision and key requirements

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

preprint2016arXiv

Planetary Signatures in the SAO 206462 (HD 135344B) Disk: A Spiral Arm Passing Through Vortex?

The disk surrounding SAO 206462, an 8 Myr-old Herbig Ae star, has recently been reported to exhibit spiral arms, an asymmetric dust continuum, and a dust-depleted inner cavity. By carrying out two-dimensional, two-fluid hydrodynamic calculations, we find that a planetary-mass companion located at the outer disk could be responsible for these observed structures. In this model, the planet excites primary and secondary arms interior to its orbit. It also carves a gap and generates a local pressure bump at the inner gap edge where a vortex forms through Rossby wave instability. The vortex traps radially drifting dust particles, forming a dust-depleted cavity in the inner disk. We propose that the vortex is responsible for the brightest southwestern peak seen in infrared scattered light and sub-millimeter dust continuum emission. In particular, it is possible that the scattered light is boosted as one of the spiral arms passes through the high density vortex region, although the vortex alone may be able to explain the peak. We suggest that a planetary companion with a mass of 10-15 M_J is orbiting SAO 206462 at 100-120 au. Monitoring of the brightest peak over the next few years will help reveal its origin because the spiral arms and vortex will show distinguishable displacement.

preprint2016arXiv

Ringed Substructure and a Gap at 1 AU in the Nearest Protoplanetary Disk

We present long-baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870 micron continuum emission from the nearest gas-rich protoplanetary disk, around TW Hya, that trace millimeter-sized particles down to spatial scales as small as 1 AU (20 mas). These data reveal a series of concentric ring-shaped substructures in the form of bright zones and narrow dark annuli (1-6 AU) with modest contrasts (5-30%). We associate these features with concentrations of solids that have had their inward radial drift slowed or stopped, presumably at local gas pressure maxima. No significant non-axisymmetric structures are detected. Some of the observed features occur near temperatures that may be associated with the condensation fronts of major volatile species, but the relatively small brightness contrasts may also be a consequence of magnetized disk evolution (the so-called zonal flows). Other features, particularly a narrow dark annulus located only 1 AU from the star, could indicate interactions between the disk and young planets. These data signal that ordered substructures on ~AU scales can be common, fundamental factors in disk evolution, and that high resolution microwave imaging can help characterize them during the epoch of planet formation.

preprint2016arXiv

Shock-driven Accretion in Circumplanetary Disks: Observables and Satellite Formation

Circumplanetary disks (CPDs) control the growth of planets, supply material for satellites to form, and provide observational signatures of young forming planets. We have carried out two dimensional hydrodynamical simulations with radiative cooling to study CPDs, and suggested a new mechanism to drive the disk accretion. Two spiral shocks are present in CPDs, excited by the central star. We find that spiral shocks can at least contribute to, if not dominate the angular momentum transport and energy dissipation in CPDs. Meanwhile, dissipation and heating by spiral shocks have a positive feedback on shock-driven accretion itself. As the disk is heated up by spiral shocks, the shocks become more open, leading to more efficient angular momentum transport. This shock driven accretion is, on the other hand, unsteady on a timescale of months/years due to production and destruction of vortices in disks. After being averaged over time, a quasi-steady accretion is reached from the planet's Hill radius all the way to the planet surface, and the disk $α$-coefficient characterizing angular momentum transport due to spiral shocks is $\sim$0.001-0.02. The disk surface density ranges from 10 to 1000 g cm$^{-2}$ in our simulations, which is at least 3 orders of magnitude smaller than the "minimum mass sub-nebula" model used to study satellite formation; instead it is more consistent with the "gas-starved" satellite formation model. Finally, we calculate the millimeter flux emitted by CPDs at ALMA and EVLA wavelength bands and predict the flux for several recently discovered CPD candidates, which suggests that ALMA is capable of discovering these accreting CPDs.

preprint2016arXiv

The Eruption of the Candidate Young Star ASASSN-15qi

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

preprint2016arXiv

The VLA view of the HL Tau Disk - Disk Mass, Grain Evolution, and Early Planet Formation

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings.

preprint2016arXiv

Very Low-Mass Stellar and Substellar Companions to Solar-like Stars From MARVELS VI: A Giant Planet and a Brown Dwarf Candidate in a Close Binary System HD 87646

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

preprint2015arXiv

An M Dwarf Companion and Its Induced Spiral Arms in the HD 100453 Protoplanetary Disk

Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to examine the scattered light morphology of the HD 100453 disk as perturbed by the companion on a circular and coplanar orbit. We find that the companion truncates the disk at ~45 AU in scattered light images, and excites two spiral arms in the remaining (circumprimary) disk with a near m=2 rotational symmetry. Both the truncated disk size and the morphology of the spirals are in excellent agreement with the SPHERE observations at Y , J, H, and K1-bands, suggesting that the M dwarf companion is indeed responsible for the observed double-spiral-arm pattern. Our model suggests that the disk is close to face on (inclination angle ~5 degree), and that the entire disk-companion system rotates counterclockwise on the sky. The HD 100453 observations, along with our modeling work, demonstrate that double spiral arm patterns in near-infrared scattered light images can be generically produced by companions, and support future observations to identify the companions responsible for the arms observed in the MWC 758 and SAO 206462 systems.

preprint2015arXiv

An Overview of Inside-Out Planet Formation

The Kepler-discovered Systems with Tightly-packed Inner Planets (STIPs), typically with several planets of Earth to super-Earth masses on well-aligned, sub-AU orbits may host the most common type of planets, including habitable planets, in the Galaxy. They pose a great challenge for planet formation theories, which fall into two broad classes: (1) formation further out followed by inward migration; (2) formation in situ, in the very inner regions of the protoplanetary disk. We review the pros and cons of these classes, before focusing on a new theory of sequential in situ formation from the inside-out via creation of successive gravitationally unstable rings fed from a continuous stream of small (~cm-m size) "pebbles," drifting inward via gas drag. Pebbles first collect at the pressure trap associated with the transition from a magnetorotational instability (MRI)-inactive ("dead zone") region to an inner, MRI-active zone. A pebble ring builds up that begins to dominate the local mass surface density of the disk and spawns a planet. The planet continues to grow, most likely by pebble accretion, until it becomes massive enough to isolate itself from the accretion flow via gap opening. This reduces the local gas density near the planet, leading to enhanced ionization and a retreat of the dead zone inner boundary. The process repeats with a new pebble ring gathering at the new pressure maximum associated with this boundary. We discuss the theory's predictions for planetary masses, relative mass scalings with orbital radius, and minimum orbital separations, and their comparison with observed systems. Finally, we discuss open questions, including potential causes of diversity of planetary system architectures, i.e., STIPs versus Solar System analogs.

preprint2015arXiv

Are protoplanetary disks born with vortices? -- Rossby wave instability driven by protostellar infall

We carry out two-fluid, two-dimensional global hydrodynamic simulations to test whether protostellar infall can trigger Rossby wave instability (RWI) in protoplanetry disks. Our results show that infall can trigger the RWI and generate vortices near the outer edge of the mass landing on the disk (i.e. centrifugal radius). We find that the RWI is triggered under a variety of conditions, although the details depend on the disk parameters and the infall pattern. The common key feature of triggering the RWI is the steep radial gradient of the azimuthal velocity induced by the local increase in density at the outer edge of the infall region. Vortices form when the instability enters the nonlinear regime. In our standard model where self-gravity is neglected, vortices merge together to a single vortex within $\sim 20$ local orbital times, and the merged vortex survives for the remaining duration of the calculation ($> 170$ local orbital times). The vortex takes part in outward angular momentum transport, with a Reynolds stress of $\lesssim10^{-2}$. Our two-fluid calculations show that vortices efficiently trap dust particles with stopping times of the order of the orbital time, locally enhancing the dust to gas ratio for particles of the appropriate size by a factor of $\sim 40$ in our standard model. When self-gravity is considered, however, vortices tend to be impeded from merging and may eventually dissipate. We conclude it may well have that protoplanetary disks have favorable conditions for vortex formation during the protostellar infall phase, which might enhance early planetary core formation.

preprint2015arXiv

Gas and dust hydrodynamical simulations of massive lopsided transition discs - II. Dust concentration

We investigate the dynamics of large dust grains in massive lopsided transition discs via 2D hydrodynamical simulations including both gas and dust. Our simulations adopt a ring-like gas density profile that becomes unstable against the Rossby-wave instability and forms a large crescent-shaped vortex. When gas self-gravity is discarded, but the indirect force from the displacement of the star by the vortex is included, we confirm that dust grains with stopping times of order the orbital time, which should be typically a few centimetres in size, are trapped ahead of the vortex in the azimuthal direction, while the smallest and largest grains concentrate towards the vortex centre. We obtain maximum shift angles of about 25 degrees. Gas self-gravity accentuates the concentration differences between small and large grains. At low to moderate disc masses, the larger the grains, the farther they are trapped ahead of the vortex. Shift angles up to 90 degrees are reached for 10 cm-sized grains, and we show that such large offsets can produce a double-peaked continuum emission observable at mm/cm wavelengths. This behaviour comes about because the large grains undergo horseshoe U-turns relative to the vortex due to the vortex's gravity. At large disc masses, since the vortex's pattern frequency becomes increasingly slower than Keplerian, small grains concentrate slightly beyond the vortex and large grains form generally non-axisymmetric ring-like structures around the vortex's radial location. Gas self-gravity therefore imparts distinct trapping locations for small and large dust grains which may be probed by current and future observations, and which suggest that the formation of planetesimals in vortices might be more difficult than previously thought.

preprint2015arXiv

Inside-Out Planet Formation. III. Planet-disk interaction at the dead zone inner boundary

The Kepler mission has discovered more than 4000 exoplanet candidates. Many are in systems with tightly packed inner planets. Inside-Out Planet Formation (IOPF) has been proposed to explain these systems. It involves sequential in situ planet formation at the local pressure maximum of a retreating dead zone inner boundary (DZIB). Pebbles accumulate at this pressure trap, which builds up a ring, and then a planet. The planet is expected to grow until it opens a gap, which helps to both truncate pebble accretion and induce DZIB retreat that sets the location of formation of the next planet. This simple scenario may be modified if the planet migrates significantly from its formation location. Thus planet-disk interactions play a crucial role in the IOPF scenario. We present numerical simulations that first assess migration of planets of various masses that are forming at the DZIB of an active accretion disk, where the effective viscosity rapidly increases in the radially inward direction. We find that the disk's torques on the planet tend to trap the planet at a location very close to the initial pressure maximum where it formed. We then study gap opening by these planets to assess at what mass a significant gap is created. Finally we present a simple model for DZIB retreat due to penetration of X-rays from the star to the disk midplane. Overall, these simulations help to quantify both the mass scale of first,"Vulcan," planet formation and the orbital separation to the location of second planet formation.

preprint2015arXiv

Observational Signatures of Planets in Protoplanetary Disks I: Gaps Opened by Single and Multiple Young Planets in Disks

It has been suggested that the gaps and cavities recently discovered in transitional disks are opened by planets. To explore this scenario, we combine two-dimensional two fluid (gas + particle) hydrodynamical calculations with three-dimensional Monte Carlo Radiative Transfer simulations, and study the observational signatures of gaps opened by one or several planets, making qualitative comparisons with observations. We find that a single planet as small as 0.2 MJ can produce a deep gap at millimeter (mm) wavelengths and almost no features at near-infrared (NIR) wavelengths, while multiple planets can open up a few *10 AU wide common gap at both wavelengths. Both the contrast ratio of the gaps and the wavelength dependence of the gap sizes are broadly consistent with data. We also confirm previous results that NIR gap sizes may be smaller than mm gap sizes due to dust-gas coupling and radiative transfer effects. When viewed at a moderate inclination angle, a physically circular on-centered gap could appear to be off-centered from the star due to shadowing. Planet-induced spiral arms are more apparent at NIR than at mm wavelengths. Overall, our results suggest that the planet-opening-gap scenario is a promising way to explain the origin of the transitional disks. Finally, inspired by the recent ALMA release of the image of the HL Tau disk, we show that multiple narrow gaps, well separated by bright rings, can be opened by 0.2 MJ planets soon after their formation in a relatively massive disk.

preprint2015arXiv

Observational Signatures of Planets in Protoplanetary Disks II: Spiral Arms Observed in Scattered Light Imaging Can be Induced by Planets

Using 3D global hydro simulations coupled with radiative transfer calculations, we study the appearance of density waves induced by giant planets in direct imaging observations at near infrared wavelengths. We find that a 6 MJ planet in a typical disk around a 1 M_sun star can produce prominent and detectable spiral arms both interior and exterior to its orbit. The inner arms have (1) two well separated arms in roughly m=2 symmetry, (2) exhibit ~10-15 degrees pitch angles, (3) ~180-270 degrees extension in the azimuthal direction, and (4) ~150% surface brightness enhancement, all broadly consistent with observed spiral arms in the SAO 206462 and MWC 758 systems. The outer arms cannot explain observations as they are too tightly wound given typical disk scale height. We confirm previous results that the outer density waves excited by a 1 MJ planet exhibit low contrast in the IR and are practically not detectable. We also find that 3D effects of the waves are important. Compared to isothermal models, density waves in adiabatic disks exhibit weaker contrast in surface density but stronger contrast in scattered light images, due to a more pronounced vertical structure in the former caused by shock heating. To drive observed pairs of arms with an external companion on a circular orbit, a massive planet, possibly a brown dwarf, is needed at around [r~0.7", PA~10 degrees] (position angle PA from north to east) in SAO 206462 and [r~0.6", PA~10 degrees] in MWC 758. Their existence may be confirmed by direct imaging planet searches.

preprint2015arXiv

Survival and Structure of Dusty Vortices in Protoplanetary Discs

We have studied the impact of dust feedback on the survival and structure of vortices in protoplanetary discs using 2-D shearing box simulations with Lagrangian dust particles. We consider dust with a variety of sizes (stopping time $t_s = 10^{-2}Ω^{-1} - 10^{2}Ω^{-1}$, from fully coupled with the gas to the decoupling limit. We find that a vortex is destroyed by dust feedback when the total dust-to-gas mass ratio within the vortex is larger than 30-50%, independent of the dust size. The dust distribution can still be asymmetric in some cases after the vortex has been destroyed. With smaller amounts of dust, a vortex can survive for at least 100 orbits, and the maximum dust surface density within the vortex can be more than 100 times larger than the gas surface density, potentially facilitating planetesimal formation. On the other hand, in these stable vortices, small ($t_s < Ω^{-1}$) and large ($t_s > Ω^{-1}$) dust grains concentrate differently and affect the gas dynamics in different ways. The distribution of large dust is more elongated than that of small dust. Large dust ($t_s > Ω^{-1}$) concentrates in the centre of the vortex and feedback leads to turn-over in vorticity towards the centre, forming a quiescent region within an anticyclonic vortex. Such a turn-over is absent if the vortex is loaded with small grains. We demonstrate that, in protoplanetary discs where both large and small dust grains are present and under the right condition, the concentration of large dust towards the vortex centre can lead to a quiescent centre, repelling the small dust and forming a small dust ring around the vortex centre. Such anticorrelations between small and large dust within vortices may explain the discrepancy between ALMA and near-IR scattered light observations in the asymmetric region of transitional discs.

preprint2015arXiv

The Structure of Spiral Shocks Excited by Planetary-mass Companions

Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3-D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.

preprint2014arXiv

Accreting Circumplanetary Disks: Observational Signatures

I calculate the spectral energy distributions (SEDs) of accreting circumplanetary disks using atmospheric radiative transfer models. Circumplanetary disks only accreting at $10^{-10} M_{\odot} yr^{-1}$ around a 1 M$_{J}$ planet can be brighter than the planet itself. A moderately accreting circumplanetary disk ($\dot{M}\sim 10^{-8}M_{\odot} yr^{-1}$; enough to form a 10 M$_{J}$ planet within 1 Myr) around a 1 M$_{J}$ planet has a maximum temperature of $\sim$2000 K, and at near-infrared wavelengths ($J$, $H$, $K$ bands), this disk is as bright as a late M-type brown dwarf or a 10 M$_{J}$ planet with a "hot start". To use direct imaging to find the accretion disks around low mass planets (e.g., 1 M$_{J}$) and distinguish them from brown dwarfs or hot high mass planets, it is crucial to obtain photometry at mid-infrared bands ($L'$, $M$, $N$ bands) because the emission from circumplanetary disks falls off more slowly towards longer wavelengths than those of brown dwarfs or planets. If young planets have strong magnetic fields ($\gtrsim$100 G), fields may truncate slowly accreting circumplanetary disks ($\dot{M}\lesssim10^{-9} M_{\odot} yr^{-1}$) and lead to magnetospheric accretion, which can provide additional accretion signatures, such as UV/optical excess from the accretion shock and line emission.

preprint2014arXiv

Accretion Outbursts in Self-gravitating Protoplanetary Disks

We improve on our previous treatments of long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e. $P\rm{d}V$ work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead-zone. The GI-induced spiral waves propagate well inside of gravitationally unstable region before they trigger outbursts at $R \lesssim 1$ AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local $α$ treatments for GI. In our standard model where zero dead-zone residual viscosity ($α_{\rm rd}$) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as $0.01$ in terms of the generic $α$ parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite $α_{\rm rd}$ triggers thermally-driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when $α_{\rm rd}=0$. The inclusion of non-zero residual viscosity in the dead-zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a quiescent accretion phase with $\dot{M}_{acc}\sim10^{-8}-10^{-7}~M_{\odot}~{\rm yr^{-1}}$ over $60~\%$ of the time and spend less than $15~\%$ of the infall phase in accretion outbursts.

preprint2014arXiv

An Analytic Model for Buoyancy Resonances in Protoplanetary Disks

Zhu, Stone, and Rafikov (2012) found in 3D shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with results of Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber k_y > 1/h (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.

preprint2014arXiv

An Observational Perspective of Transitional Disks

Transitional disks are objects whose inner disk regions have undergone substantial clearing. The Spitzer Space Telescope produced detailed spectral energy distributions (SEDs) of transitional disks that allowed us to infer their radial dust disk structure in some detail, revealing the diversity of this class of disks. The growing sample of transitional disks also opened up the possibility of demographic studies, which provided unique insights. There now exist (sub)millimeter and infrared images that confirm the presence of large clearings of dust in transitional disks. In addition, protoplanet candidates have been detected within some of these clearings. Transitional disks are thought to be a strong link to planet formation around young stars and are a key area to study if further progress is to be made on understanding the initial stages of planet formation. Here we provide a review and synthesis of transitional disk observations to date with the aim of providing timely direction to the field, which is about to undergo its next burst of growth as ALMA reaches its full potential. We discuss what we have learned about transitional disks from SEDs, color-color diagrams, and imaging in the (sub)mm and infrared. We then distill the observations into constraints for the main disk clearing mechanisms proposed to date (i.e., photoevaporation, grain growth, and companions) and explore how the expected observational signatures from these mechanisms, particularly planet-induced disk clearing, compare to actual observations. Lastly, we discuss future avenues of inquiry to be pursued with ALMA, JWST, and next generation of ground-based telescopes.

preprint2014arXiv

Dust Transport in MRI Turbulent Disks: Ideal and Non-ideal MHD with Ambipolar Diffusion

We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. The turbulence is driven by the magnetorotational instability (MRI) with either ideal or non-ideal MHD that includes ambipolar diffusion (AD). In ideal MHD simulations, the surface density evolution (except for dust that drifts fastest), turbulent diffusion, and vertical scale height of dust can all be reproduced by simple one-dimensoinal and/or analytical models. However, in AD dominated simulations which simulate protoplanetary disks beyond 10s of AU, the vertical scale height of dust is larger than previously predicted. To understand this anomaly in more detail, we carry out both unstratified and stratified local shearing box simulations with Lagrangian particles, and find that turbulence in AD dominated disks has very different properties (e.g., temporal autocorrelation functions and power spectra) than turbulence in ideal MHD disks, which leads to quite different particle diffusion efficiency. For example, MRI turbulence with AD has a longer correlation time for the vertical velocity, which causes significant vertical particle diffusion and large dust scale height. In ideal MHD the Schmidt numbers ($Sc$) for radial and vertical turbulent diffusion are $Sc_{r}\sim 1$ and $Sc_{z}\gtrsim 3$, but in the AD dominated regime both $Sc_{r}$ and $Sc_{z}$ are $\lesssim 1$. Particle concentration in pressure bumps induced by MRI turbulence has also been studied. Since non-ideal MHD effects dominate most regions in protoplanetary disks, our study suggests that modeling dust transport in turbulence driven by MRI with non-ideal MHD effects is important for understanding dust transport in realistic protoplanetary disks.

preprint2014arXiv

Dust Trapping by Vortices in Transitional Disks: Evidence for Non-ideal MHD Effects in Protoplanetary Disks

We study particle trapping at the edge of a gap opened by a planet in a protoplanetary disk. In particular, we explore the effects of turbulence driven by the magnetorotational instability on particle trapping, using global three-dimensional magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. We study disks either in the ideal MHD limit or dominated by ambipolar diffusion (AD) that plays an essential role at the outer regions of a protoplanetary disk. With ideal MHD, strong turbulence (the equivalent viscosity parameter $α\sim 10^{-2}$) in disks prevents vortex formation at the edge of the gap opened by a 9 $M_{J}$ planet, and most particles (except the particles that drift fastest) pile up at the outer gap edge almost axisymmetrically. When AD is considered, turbulence is significantly suppressed ($α\lesssim 10^{-3}$), and a large vortex forms at the edge of the planet induced gap, which survives $\sim$ 1000 orbits. The vortex can efficiently trap dust particles that span 3 orders of magnitude in size within 100 planetary orbits. We have also carried out two-dimensional hydrodynamical simulations using viscosity as an approximation to MHD turbulence. These hydrodynamical simulations can reproduce vortex generation at the gap edge as seen in MHD simulations. Finally, we use our simulation results to generate synthetic images for ALMA dust continuum observations on Oph IRS 48 and HD 142527, which show good agreement with existing observations. Predictions for future ALMA cycle 2 observations have been made. We conclude that the asymmetry in ALMA observations can be explained by dust trapping vortices and the existence of vortices could be the evidence that the outer protoplanetary disks are dominated by ambipolar diffusion with $α<10^{-3}$ at the disk midplane.

preprint2014arXiv

Episodic Accretion in Young Stars

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

preprint2014arXiv

On the Metallicities of Kepler Stars

We use 12000 stars from Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectroscopic data to show that the metallicities of Kepler field stars as given in the Kepler Input Catalog (KIC) systematically underestimate both the true metallicity and the dynamic range of the Kepler sample. Specifically, to the first order approximation, we find [Fe/H]_KIC = -0.20 + 0.43 [Fe/H]_LAMOST, with a scatter of ~0.25 dex, due almost entirely to errors in KIC. This relation is most secure for -0.3<[Fe/H]_LAMOST<+0.4 where we have >200 comparison stars per 0.1 dex bin and good consistency is shown between metallicities determined by LAMOST and high-resolution spectra. It remains approximately valid in a slightly broader range. When the relation is inverted, the error in true metallicity as derived from KIC is (0.25 dex)/0.43~0.6 dex. We thereby quantitatively confirm the cautionary note by Brown et al. (2011) that KIC estimates of [Fe/H] should not be used by "anyone with a particular interest in stellar metallicities". Fortunately, many more LAMOST spectroscopic metallicities will be available in the near future.

preprint2014arXiv

Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks

We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low mass planet ($8 M_{\oplus}$), two narrow gaps start to open in the gas on each side of the planet where the density waves shock. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet ($>0.1 M_{J}$) carves out a deeper gap with sharp edges, which are unstable to the vortex formation. Particles with a wide range of sizes ($0.02<Ωt_{s}<20$) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with $Ωt_{s}\sim 1$. The dust concentration is highly elongated in the $ϕ$ direction, and can be as wide as 4 disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10$^2$ in a very non-axisymmetric fashion. For very big particles ($Ωt_{s}\gg 1$) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap opening there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g. the "dead zone") dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with $M_{p}<M_{J}$, potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and EVLA.

preprint2014arXiv

The Science Case for the Planet Formation Imager (PFI)

Among the most fascinating and hotly-debated areas in contemporary astrophysics are the means by which planetary systems are assembled from the large rotating disks of gas and dust which attend a stellar birth. Although important work has already been, and is still being done both in theory and observation, a full understanding of the physics of planet formation can only be achieved by opening observational windows able to directly witness the process in action. The key requirement is then to probe planet-forming systems at the natural spatial scales over which material is being assembled. By definition, this is the so-called Hill Sphere which delineates the region of influence of a gravitating body within its surrounding environment. The Planet Formation Imager project (PFI) has crystallized around this challenging goal: to deliver resolved images of Hill-Sphere-sized structures within candidate planet-hosting disks in the nearest star-forming regions. In this contribution we outline the primary science case of PFI. For this purpose, we briefly review our knowledge about the planet-formation process and discuss recent observational results that have been obtained on the class of transition disks. Spectro-photometric and multi-wavelength interferometric studies of these systems revealed the presence of extended gaps and complex density inhomogeneities that might be triggered by orbiting planets. We present detailed 3-D radiation-hydrodynamic simulations of disks with single and multiple embedded planets, from which we compute synthetic images at near-infrared, mid-infrared, far-infrared, and sub-millimeter wavelengths, enabling a direct comparison of the signatures that are detectable with PFI and complementary facilities such as ALMA. From these simulations, we derive some preliminary specifications that will guide the array design and technology roadmap of the facility.

preprint2013arXiv

Fast Rise of "Neptune-Size" Planets ($4-8 R_{\rm Earth}$) from $P\sim10$ to $\sim250$ days -- Statistics of Kepler Planet Candidates Up to $\sim 0.75 {\rm AU}$

We infer the period ($P$) and size ($R_p$) distribution of Kepler transiting planet candidates with $R_p\ge 1 R_{\rm Earth}$ and $P < 250$ days hosted by solar-type stars. The planet detection efficiency is computed by using measured noise and the observed timespans of the light curves for $\sim 120,000$ Kepler target stars. We focus on deriving the shape of planet period and radius distribution functions. We find that for orbital period $P>10$ days, the planet frequency d$N_p$/d$\log$P for "Neptune-size" planets ($R_p = 4-8 R_{\rm Earth}$) increases with period as $\propto P^{0.7\pm0.1}$. In contrast, d$N_p$/d$\log$P for "super-Earth-size" ($2-4 R_{\rm Earth}$) as well as "Earth-size" ($1-2 R_{\rm Earth}$) planets are consistent with a nearly flat distribution as a function of period ($\propto P^{0.11\pm0.05}$ and $\propto P^{-0.10\pm0.12}$, respectively), and the normalizations are remarkably similar (within a factor of $\sim 1.5$ at $50 $ days). Planet size distribution evolves with period, and generally the relative fractions for big planets ($\sim 3-10 R_{\rm Earth}$) increase with period. The shape of the distribution function is not sensitive to changes in selection criteria of the sample. The implied nearly flat or rising planet frequency at long period appears to be in tension with the sharp decline at $\sim 100$ days in planet frequency for low mass planets (planet mass $m_p < 30 M_{\rm Earth}$) recently suggested by HARPS survey. Within $250$ days, the cumulative frequencies for Earth-size and super-Earth-size planets are remarkably similar ($\sim 28 %$ and $25%$), while Neptune-size and Jupiter-size planets are $\sim 7%$, and $\sim 3%$, respectively. A major potential uncertainty arises from the unphysical impact parameter distribution of the candidates.

preprint2013arXiv

Long-term Evolution of Photoevaporating Protoplanetary Disks

We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically-thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting of $\lesssim20 %$ of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly-rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the $\dot{M}-R_{\rm wall}$ plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of order 10 AU.

preprint2013arXiv

Low mass planets in protoplanetary disks with net vertical magnetic fields: the Planetary Wake and Gap Opening

We study wakes and gap opening by low mass planets in gaseous protoplanetary disks threaded by net vertical magnetic fields which drive magnetohydrodynamical (MHD) turbulence through the magnetorotational instabilty (MRI), using three dimensional simulations in the unstratified local shearing box approximation. The wakes, which are excited by the planets, are damped by shocks similar to the wake damping in inviscid hydrodynamic (HD) disks. Angular momentum deposition by shock damping opens gaps in both MHD turbulent disks and inviscid HD disks even for low mass planets, in contradiction to the "thermal criterion" for gap opening. To test the "viscous criterion", we compared gap properties in MRI-turbulent disks to those in viscous HD disks having the same stress, and found that the same mass planet opens a significantly deeper and wider gap in net vertical flux MHD disks than in viscous HD disks. This difference arises due to the efficient magnetic field transport into the gap region in MRI disks, leading to a larger effective αwithin the gap. Thus, across the gap, the Maxwell stress profile is smoother than the gap density profile, and a deeper gap is needed for the Maxwell stress gradient to balance the planetary torque density. We also confirmed the large excess torque close to the planet in MHD disks, and found that long-lived density features (termed zonal flows) produced by the MRI can affect planet migration. The comparison with previous results from net toroidal flux/zero flux MHD simulations indicates that the magnetic field geometry plays an important role in the gap opening process. Overall, our results suggest that gaps can be commonly produced by low mass planets in realistic protoplanetary disks, and caution the use of a constant α-viscosity to model gaps in protoplanetary disks.

preprint2012arXiv

ALMA and VLA observations of the outflows in IRAS 16293-2422

We present ALMA and VLA observations of the molecular and ionized gas at 0.1-0.3 arcsec resolution in the Class 0 protostellar system IRAS 16293-2422. These data clarify the origins of the protostellar outflows from the deeply embedded sources in this complex region. Source A2 is confirmed to be at the origin of the well known large scale north-east--south-west flow. The most recent VLA observations reveal a new ejection from that protostar, demonstrating that it drives an episodic jet. The central compact part of the other known large scale flow in the system, oriented roughly east-west, is well delineated by the CO(6-5) emission imaged with ALMA and is confirmed to be driven from within component A. Finally, a one-sided blueshifted bubble-like outflow structure is detected here for the first time from source B to the north-west of the system. Its very short dynamical timescale (~ 200 yr), low velocity, and moderate collimation support the idea that source B is the youngest object in the system, and possibly one of the youngest protostars known.

preprint2012arXiv

Dust Filtration by Planet-Induced Gap Edges: Implications for Transitional Disks

By carrying out two-dimensional two-fluid global simulations, we have studied the response of dust to gap formation by a single planet in the gaseous component of a protoplanetary disk - the so-called "dust filtration" mechanism. We have found that a gap opened by a giant planet at 20 AU in a α=0.01, \dot{M}=10^{-8} Msun/yr disk can effectively stop dust particles larger than 0.1 mm drifting inwards, leaving a sub-millimeter dust cavity/hole. However, smaller particles are difficult to filter by a planet-induced gap due to 1) dust diffusion, and 2) a high gas accretion velocity at the gap edge. An analytic model is also derived to understand what size particles can be filtered by the gap edge. Finally, with our updated understanding of dust filtration, we have computed Monte-Carlo radiative transfer models with variable dust size distributions to generate the spectral energy distributions (SEDs) of disks with gaps. By comparing with transitional disk observations (e.g. GM Aur), we have found that dust filtration alone has difficulties to deplete small particles sufficiently to explain the near-IR deficit of transitional disks, except under some extreme circumstances. The scenario of gap opening by multiple planets studied previously suffers the same difficulty. One possible solution is by invoking both dust filtration and dust growth in the inner disk. In this scenario, a planet induced gap filters large dust particles in the disk, and the remaining small dust particles passing to the inner disk can grow efficiently without replenishment from fragmentation of large grains. Predictions for ALMA have also been made based on all these scenarios. We conclude that dust filtration with planet(s) in the disk is a promising mechanism to explain submm observations of transitional disks but it may need to be combined with other processes (e.g. dust growth) to explain the near-IR deficit.

preprint2012arXiv

Planet-Disk interaction in 3D: the importance of buoyancy waves

We carry out local three dimensional (3D) hydrodynamic simulations of planet-disk interaction in stratified disks with varied thermodynamic properties. We find that whenever the Brunt-Vaisala frequency (N) in the disk is nonzero, the planet exerts a strong torque on the disk in the vicinity of the planet, with a reduction in the traditional "torque cutoff". In particular, this is true for adiabatic perturbations in disks with isothermal density structure, as should be typical for centrally irradiated protoplanetary disks. We identify this torque with buoyancy waves, which are excited (when N is non-zero) close to the planet, within one disk scale height from its orbit. These waves give rise to density perturbations with a characteristic 3D spatial pattern which is in close agreement with the linear dispersion relation for buoyancy waves. The torque due to these waves can amount to as much as several tens of per cent of the total planetary torque, which is not expected based on analytical calculations limited to axisymmetric or low-m modes. Buoyancy waves should be ubiquitous around planets in the inner, dense regions of protoplanetary disks, where they might possibly affect planet migration.

preprint2012arXiv

The Structure of Pre-transitional Protoplanetary Disks I: Radiative Transfer Modeling of the Disk+Cavity in the PDS 70 system

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

preprint2012arXiv

Variable Accretion Outbursts in Protostellar Evolution

We extend the one-dimensional, two-zone models of long-term protostellar disk evolution with infall of Zhu et al. to consider the potential effects of a finite viscosity in regions where the ionization is too low for the magnetorotational instability (MRI) to operate (the "dead zone"). We find that the presence of a small but finite dead zone viscosity, as suggested by simulations of stratified disks with MRI-active outer layers, can trigger inside-out bursts of accretion, starting at or near the inner edge of the disk, instead of the previously-found outside-in bursts with zero dead zone viscosity, which originate at a few AU in radius. These inside-out bursts of accretion bear a qualitative resemblance to the outburst behavior of one FU Ori object, V1515 Cyg, in contrast to the outside-in burst models which more closely resemble the accretion events in FU Ori and V1057 Cyg. Our results suggest that the type and frequency of outbursts are potentially a probe of transport efficiency in the dead zone. Simulations must treat the inner disk regions, $R \lesssim 0.5$ AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

preprint2011arXiv

Challenges in Forming Planets by Gravitational Instability: Disk Irradiation and Clump Migration, Accretion & Tidal Destructio

We present two-dimensional hydrodynamic simulations of self-gravitating protostellar disks subject to axisymmetric infall from envelopes and irradiation from the central star, to explore disk fragmentation due to gravitational instability (GI), and the fragmented clump evolution. We assume that the disk is built gradually and smoothly by the infall, resulting in good numerical convergence. We confirm that for disks around solar-mass stars, infall at high rates at radii beyond ~50 AU leads to disk fragmentation. At lower infall rates <1e-5 Msun/yr, however, irradiation suppresses fragmentation. We find that, once formed, the fragments or clumps migrate inward on typical type-I time scales of ~2e3 yr initially, but later migration deviates from the type-I time scale when the clump becomes more massive than the local disk mass, and/or when they starts to open gaps. As they migrate, the clumps accrete from the disk at a rate 1e-3 to 1e-1 MJupiter/yr, consistent with analytic estimates that assume a 1-2 Hill radii cross section. Most clumps can grow to their isolation masses >0.1 Msun quickly. The eventual fates of these clumps, however, diverges depending on the migration speed: 3 out of 13 clumps become massive enough (brown dwarf mass) to open gaps in the disk and essentially stop migrating; 4 are tidally destroyed during inward migration; 6 migrate across the inner simulated disk boundary. A simple analytic model for clump evolution is derived to explain these different fates. Overall, our results indicate that fast migration, accretion, and tidal destruction of the clumps pose challenges to the scenario of giant planet formation by GI in situ, but may provide a formation mechanism for close binary systems.

preprint2011arXiv

On Rapid Disk Accretion and Initial Conditions in Protostellar Evolution

Low-mass protostars may accrete most of their material through short-lived episodes of rapid disk accretion; yet until recently evolutionary tracks for these protostars assumed only constant or slowly-varying accretion. Important initial steps toward examining the potential effects of rapid accretion were recently made by Baraffe, Chabrier, & Gallardo, who showed that in the limit of low-temperature ("cold") accretion, protostars may have much smaller radii than found in previous treatments. Such small radii at the end of protostellar accretion would have the effect of making some young stars appear much older - perhaps as much as 10 Myr - than they really are. However, we argue that very rapid disk accretion is unlikely to be cold, because observations of the best-studied pre-main sequence disks with rapid disk accretion outbursts - the FU Ori objects - have spectral energy distributions which imply large, not small, protostellar radii. In addition, theory indicates models that at high accretion rates, protostellar disks become internally hot and geometrically thick, making it much more likely that hot material is added to the star. In addition, the very large luminosity of the accretion disk is likely to irradiate the central star strongly, heating up the outer layers and potentially expanding them. Nevertheless, the Baraffe et al. calculations emphasize the importance of initial protostellar radii for subsequent evolution.

preprint2010arXiv

Long-term Evolution of Protostellar and Protoplanetary Disks. I. Outbursts

As an initial investigation into the long-term evolution of protostellar disks, we explore the conditions required to explain the large outbursts of disk accretion seen in some young stellar objects. We use one-dimensional time-dependent disk models with a phenomenological treatment of the magnetorotational instability (MRI) and gravitational torques to follow disk evolution over long timescales. Comparison with our previous two-dimensional disk model calculations (Zhu et al. 2009b, Z2009b) indicates that the neglect of radial effects and two-dimensional disk structure in the one-dimensional case makes only modest differences in the results; this allows us to use the simpler models to explore parameter space efficiently. We find that the mass infall rates typically estimated for low-mass protostars generally result in AU-scale disk accretion outbursts, as predicted by our previous analysis (Zhu et al. 2009a,Z2009a). We also confirm quasi-steady accretion behavior for high mass infall rates if the values of $α$-parameter for the magnetorotational instability is small, while at this high accretion rate convection from the thermal instability may lead to some variations. We further constrain the combinations of the $α$-parameter and the MRI critical temperature, which can reproduce observed outburst behavior. Our results suggest that dust sublimation may be connected with full activation of the MRI. This is consistent with the idea that small dust captures ions and electrons to suppress the MRI. In a later paper we will explore both long-term outburst and disk evolution with this model, allowing for infall from protostellar envelopes with differing angular momenta.

preprint2010arXiv

Long-term Evolution of Protostellar and Protoplanetary Disks. II. Layered Accretion with Infall

We use one-dimensional two-zone time-dependent accretion disk models to study the long-term evolution of protostellar disks subject to mass addition from the collapse of a rotating cloud core. Our model consists of a constant surface density magnetically coupled active layer, with transport and dissipation in inactive regions only via gravitational instability. We start our simulations after a central protostar has formed, containing ~ 10% of the mass of the protostellar cloud. Subsequent evolution depends on the angular momentum of the accreting envelope. We find that disk accretion matches the infall rate early in the disk evolution because much of the inner disk is hot enough to couple to the magnetic field. Later infall reaches the disk beyond ~10 AU, and the disk undergoes outbursts of accretion in FU Ori-like events as described in Zhu et al. 2009c. If the initial cloud core is moderately rotating most of the central star's mass is built up by these outburst events. Our results suggest that the protostellar "luminosity problem" is eased by accretion during these FU Ori-like outbursts. After infall stops the disk enters the T Tauri phase. An outer, viscously evolving disk has structure that is in reasonable agreement with recent submillimeter studies and its surface density evolves from $Σ\propto R^{-1}$ to $R^{-1.5}$. An inner, massive belt of material-- the "dead zone" -- would not have been observed yet but should be seen in future high angular resolution observations by EVLA and ALMA. This high surface density belt is a generic consequence of low angular momentum transport efficiency at radii where the disk is magnetically decoupled, and would strongly affect planet formation and migration.

preprint2010arXiv

Transitional and Pre-Transitional disks: Gap Opening by Multiple Planets?

We use 2D hydrodynamic simulations of viscous disks to examine whether dynamically-interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence of inner disk dust depletion, we find that it requires three to four giant planets to open up large enough gaps to be consistent with inferences from spectral energy distributions, because the gap width is limited by the tendency of the planets to be driven together into 2:1 resonances. With very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In these cases, however, the reduction of surface density produces a corresponding reduction in the inner disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of the pre/transitional disks. We find that even with four planets in disks, additional substantial dust depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, and an earlier history of dust growth is consistent with the hypothesis that pre/transitional disks are explained by the presence of giant planets. We conclude that the depths and widths of gaps, and disk accretion rates in pre/transitional disks cannot be reproduced by a planet-induced gap opening scenario alone. Significant dust depletion is also required within the gaps/holes. Order of magnitude estimates suggest the mass of small dust particles (<1 micron) relative to the gas must be depleted to 1e-5 -- 1e-2 of the interstellar medium value, implying a very efficient mechanism of small dust removal or dust growth.

preprint2009arXiv

2-D simulations of FU Orionis disk outbursts

We have developed time-dependent models of FU Ori accretion outbursts to explore the physical properties of protostellar disks. Our two-dimensional, axisymmetric models incorporate full vertical structure with a new treatment of the radiative boundary condition for the disk photosphere. We find that FU Ori-type outbursts can be explained by a slow accumulation of matter due to gravitational instability. Eventually this triggers the magnetorotational instability, which leads to rapid accretion. The thermal instability is triggered in the inner disk but this instability is not necessary for the outburst. An accurate disk vertical structure, including convection, is important for understanding the outburst behavior. Large convective eddies develop during the high state in the inner disk. The models are in agreement with Spitzer IRS spectra and also with peak accretion rates and decay timescales of observed outbursts, though some objects show faster rise timescale. We also propose that convection may account for the observed mild-supersonic turbulence and the short-timescale variations of FU Orionis objects.

preprint2008arXiv

Non-steady Accretion in Protostars

Observations indicate that mass accretion rates onto low-mass protostars are generally lower than the rates of infall to their disks; this suggests that much of the protostellar mass must be accreted during rare, short outbursts of rapid accretion. We explore when protostellar disk accretion is likely to be highly variable. While constant $α$ disks can in principle adjust their accretion rates to match infall rates, protostellar disks are unlikely to have constant $α$. In particular we show that neither models with angular momentum ransport due solely to the magnetorotational instability (MRI) nor ravitational instability (GI) are likely to transport disk mass at rotostellar infall rates over the large range of radii needed to move infalling envelope material down to the central protostar. We show that the MRI and GI are likely to combine to produce outbursts of rapid accretion starting at a few AU. Our analysis is consistent with the time-dependent models of Armitage, Livio, & Pringle (2001) and agrees with our observational study of the outbursting object FU Ori.

Zhaohuan Zhu

What is connected

Connect this record

See the researcher in context

Building this map preview

78 published item(s)

Brightness variability in polar circumbinary disks

Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning

Structured Distributions of Gas and Solids in Protoplanetary Disks

ALMA Discovery of a Disk around the Planetary-mass Companion SR 12 c

Direct Formation of Planetary Embryos in Self-Gravitating Disks

Formation of Dust Rings and Gaps in Non-ideal MHD Disks Through Meridional Gas Flows

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

HST/WFC3 H$α$ Direct-Imaging Detection of a Point-like Source in the Disk Cavity of AB Aur

Investigating the Future Potential of an Upgraded ALMA to Image Planet Forming Disks at Sub-au Scales

Maximum Temperatures in Evolving Protoplanetary Discs and Composition of Planetary Building Blocks

Thermal instability in radiation hydrodynamics: instability mechanisms, position-dependent S-curves, and attenuation curves

Constraining Protoplanetary Disk Accretion and Young Planets Using ALMA Kinematic Observations

Kozai-Lidov oscillations triggered by a tilt instability of detached circumplanetary discs

The ALMA Early Science view of FUor/EXor objects. I. Through the looking-glass of V2775 Ori

A fast-growing tilt instability of detached circumplanetary disks

A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk

A Study of Millimeter Variability in FUor Objects

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

ALMA 0.88 mm Survey of Disks around Planetary-mass Companions

Dust Condensation in Evolving Discs and the Composition of Planetary Building Blocks

Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori's Accretion Disk and Observational Signatures of Magnetic Fields

Large-scale CO spiral arms and complex kinematics associated with the T Tauri star RU Lup

Near-Infrared Imaging of a Spiral in the CQ Tau Disk

The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures

The effects of disk self-gravity and radiative cooling on the formation of gaps and spirals by young planets

The Planetary Luminosity Problem: "Missing Planets" and the Observational Consequences of Episodic Accretion

Circumplanetary Disk Dynamics in the Isothermal and Adiabatic Limits

Irregular dust features around intermediate-mass young stars with GPI: signs of youth or misaligned disks?

Morphological signatures induced by dust back reaction in discs with an embedded planet

Resolving the FU Ori System with ALMA: Interacting Twin Disks?

A Dwarf Protoplanetary Disk around XZ Tau B

Architecture design study and technology roadmap for the Planet Formation Imager (PFI)

Dust dynamics in 2D gravito-turbulent disks

Exoplanet Orbital Eccentricities Derived From LAMOST-Kepler Analysis

Gas and dust hydrodynamical simulations of massive lopsided transition discs - I. Gas distribution

Global MHD Simulations of Accretion Disks in Cataclysmic Variables (CVs): I. The Importance of Spiral Shocks

Imaging the water-snow line during a protostellar outburst

Planet Formation Imager (PFI): science vision and key requirements

Planetary Signatures in the SAO 206462 (HD 135344B) Disk: A Spiral Arm Passing Through Vortex?

Ringed Substructure and a Gap at 1 AU in the Nearest Protoplanetary Disk

Shock-driven Accretion in Circumplanetary Disks: Observables and Satellite Formation

The Eruption of the Candidate Young Star ASASSN-15qi

The VLA view of the HL Tau Disk - Disk Mass, Grain Evolution, and Early Planet Formation

Very Low-Mass Stellar and Substellar Companions to Solar-like Stars From MARVELS VI: A Giant Planet and a Brown Dwarf Candidate in a Close Binary System HD 87646

An M Dwarf Companion and Its Induced Spiral Arms in the HD 100453 Protoplanetary Disk

An Overview of Inside-Out Planet Formation

Are protoplanetary disks born with vortices? -- Rossby wave instability driven by protostellar infall

Gas and dust hydrodynamical simulations of massive lopsided transition discs - II. Dust concentration

Inside-Out Planet Formation. III. Planet-disk interaction at the dead zone inner boundary

Observational Signatures of Planets in Protoplanetary Disks I: Gaps Opened by Single and Multiple Young Planets in Disks

Observational Signatures of Planets in Protoplanetary Disks II: Spiral Arms Observed in Scattered Light Imaging Can be Induced by Planets

Survival and Structure of Dusty Vortices in Protoplanetary Discs

The Structure of Spiral Shocks Excited by Planetary-mass Companions

Accreting Circumplanetary Disks: Observational Signatures

Accretion Outbursts in Self-gravitating Protoplanetary Disks

An Analytic Model for Buoyancy Resonances in Protoplanetary Disks

An Observational Perspective of Transitional Disks

Dust Transport in MRI Turbulent Disks: Ideal and Non-ideal MHD with Ambipolar Diffusion

Dust Trapping by Vortices in Transitional Disks: Evidence for Non-ideal MHD Effects in Protoplanetary Disks

Episodic Accretion in Young Stars

On the Metallicities of Kepler Stars

Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks

The Science Case for the Planet Formation Imager (PFI)

Fast Rise of "Neptune-Size" Planets ($4-8 R_{\rm Earth}$) from $P\sim10$ to $\sim250$ days -- Statistics of Kepler Planet Candidates Up to $\sim 0.75 {\rm AU}$

Long-term Evolution of Photoevaporating Protoplanetary Disks

Low mass planets in protoplanetary disks with net vertical magnetic fields: the Planetary Wake and Gap Opening

ALMA and VLA observations of the outflows in IRAS 16293-2422

Dust Filtration by Planet-Induced Gap Edges: Implications for Transitional Disks

Planet-Disk interaction in 3D: the importance of buoyancy waves

The Structure of Pre-transitional Protoplanetary Disks I: Radiative Transfer Modeling of the Disk+Cavity in the PDS 70 system

Variable Accretion Outbursts in Protostellar Evolution

Challenges in Forming Planets by Gravitational Instability: Disk Irradiation and Clump Migration, Accretion & Tidal Destructio

On Rapid Disk Accretion and Initial Conditions in Protostellar Evolution

Long-term Evolution of Protostellar and Protoplanetary Disks. I. Outbursts