Source author record

Jaehan Bae

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

30 published item(s)

preprint2026arXiv

Circumplanetary Disk Candidate in the Disk of HD 163296 Traced by Localized Emission from Simple Organics

Atacama Large Millimeter/submillimeter Array observations suggest that the disc of HD 163296 is being actively shaped by embedded, yet unseen protoplanets, as indicated by numerous gas and dust substructures consistent with planet-disc interaction models. We report the first detection of simple organic molecules, HCN and C2H, tracing a candidate circumplanetary disc (CPD) in the HD 163296 system, located at an orbital radius of $R=88\pm7$ au and azimuth $ϕ=46\pm3^\circ$ (or $R=0.75''$, $\rm{PA}=350^\circ$ in projected sky coordinates), and originating near the midplane of the circumstellar disc. The signature is localised but spectrally resolved, and it overlaps with a previously reported planet candidate, P94, identified through kinematic perturbations traced by CO lines. We propose a scenario in which the observed chemical anomalies arise from increased heating driven by the forming planet and ongoing accretion through its CPD, facilitating the thermal desorption of species that would otherwise remain frozen out in the disc midplane, and potentially triggering the activation barriers of chemical reactions that lead to enhanced molecular production. Based on a first-order dynamical analysis of the HCN spectrum from the CPD--isolated with a 7$σ$ significance--we infer an upper limit on the planet mass of 1.8 $M_{\rm Jup}$, consistent with predictions from CO kinematics and constraints from direct imaging studies. By comparing the CPD sizes derived from our models with theoretical expectations where the CPD radius corresponds to roughly one-third of the planet's Hill radius, we favor CPD gas temperatures $T > 150$ K, planet masses $M_{\rm p} < 1.0$ $M_{\rm Jup}$, and CPD radii $R_{\rm CPD} < 2$ au.

preprint2026arXiv

Tracing Pebble Drift History in Two Protoplanetary Disks with CO Enhancement

Pebble drift is an important mechanism for supplying the materials needed to build planets in the inner region of protoplanetary disks. Thus, constraining pebble drift's timescales and mass flux is essential to understanding planet formation history. Current pebble drift models suggest pebble fluxes can be constrained from the enhancement of gaseous volatile abundances when icy pebbles sublimate after drifting across key snowlines. In this work, we present ALMA observations of spatially resolved $^{13}$C$^{18}$O J=2-1 line emission inside the midplane CO snowline of the HD 163296 and MWC 480 protoplanetary disks. We use radiative transfer and thermochemical models to constrain the spatial distribution of CO gas column density. We find that both disks display centrally peaked CO abundance enhancement of up to ten times of ISM abundance levels. For HD 163296 and MWC 480, the inferred enhancements require 250-350 and 480-660 Earth Masses of pebbles to have drifted across their CO snowlines, respectively. These ranges fall within cumulative pebble mass flux ranges to grow gas giants in the interior to the CO snowline. The centrally peaked CO enhancement is unexpected in current pebble drift models, which predict CO enhancement peaks at the CO snowline or is uniform inside the snowline. We propose two hypotheses to explain the centrally-peaked CO enhancement, including a large CO desorption distance and CO trapped in water ice. By testing both hypotheses with the 1D gas and dust evolution code chemcomp, we find that volatile trapping (about 30\%) best reproduces the centrally peaked CO enhancement observed.

preprint2024arXiv

Mapping the Vertical Gas Structure of the Planet-hosting PDS 70 Disk

PDS 70 hosts two massive, still-accreting planets and the inclined orientation of its protoplanetary disk presents a unique opportunity to directly probe the vertical gas structure of a planet-hosting disk. Here, we use high-spatial-resolution (${\approx}$0."1;10 au) observations in a set of CO isotopologue lines and HCO$^+$ J=4-3 to map the full 2D $(r,z)$ disk structure from the disk atmosphere, as traced by $^{12}$CO, to closer to the midplane, as probed by less abundant isotopologues and HCO$^+$. In the PDS 70 disk, $^{12}$CO traces a height of $z/r\approx0.3$, $^{13}$CO is found at $z/r\approx0.1$, and C$^{18}$O originates at, or near, the midplane. The HCO$^+$ surface arises from $z/r\approx0.2$ and is one of the few non-CO emission surfaces constrained with high fidelity in disks to date. In the $^{12}$CO J=3-2 line, we resolve a vertical dip and steep rise in height at the cavity wall, making PDS 70 the first transition disk where this effect is directly seen in line emitting heights. In the outer disk, the CO emission heights of PDS 70 appear typical for its stellar mass and disk size and are not substantially altered by the two inner embedded planets. By combining CO isotopologue and HCO$^+$ lines, we derive the 2D gas temperature structure and estimate a midplane CO snowline of ${\approx}$56-85 au. This implies that both PDS 70b and 70c are located interior to the CO snowline and are likely accreting gas with a high C/O ratio of ${\approx}$1.0, which provides context for future planetary atmospheric measurements from, e.g., JWST, and for properly modeling their formation histories.

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

Anisotropic Infall and Substructure formation in Embedded Disks

The filamentary nature of accretion streams found around embedded sources suggest that protostellar disks experience heterogenous infall from the star-forming environment, consistent with the accretion behavior onto star-forming cores in top-down star-cluster formation simulations. This may produce disk substructures in the form of rings, gaps, and spirals continuing to be identified by high-resolution imaging surveys in both embedded Class 0/I and later Class II sources. We present a parameter study of anisotropic infall, informed by the properties of accretion flows onto protostellar cores in numerical simulations, and varying the relative specific angular momentum of incoming flows as well as their flow geometry. Our results show that anisotropic infall perturbs the disk and readily launches the Rossby wave instability (RWI). It forms vortices at the inner and outer edge of the infall zone where material is deposited. These vortices drive spiral waves and angular momentum transport, with some models able to drive stresses corresponding to a viscosity parameter on the order of $α\sim 10^{-2}$. The resulting azimuthal shear forms robust pressure bumps that act as barriers to radial drift of dust grains, as demonstrated by post-processing calculations of drift-dominated dust evolution. We discuss how a self-consistent model of anisotropic infall can account for the formation of millimeter rings in the outer disk as well as producing compact dust disks, consistent with observations of embedded sources.

preprint2022arXiv

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$.

preprint2022arXiv

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): A Panchromatic View of DO Tau's Complex Kilo-au Environment

While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales.

preprint2022arXiv

Gas and Dust Shadows in the TW Hydrae Disk

We present new observations of CO J=2-1 emission from the protoplanetary disk around TW Hya. Emission is detected out to 240 au (4") and found to exhibit azimuthal variations up to 20% beyond 180 au (3"), with the west side of the disk brighter than the east. This asymmetry is interpreted as tracing the shadow previously seen in scattered light. A reanalysis of the multi-epoch observations of the dust shadow in scattered light from Debes et al. (2017) suggests that an oscillatory motion would provide a better model of the temporal evolution of the dust shadow rather than orbital motion. Both models predict an angular offset between the dust shadow and the gas shadow of up to ~100 deg. We attribute this offset to the finite rate at which dust grains and gas molecules can exchange heat, dominated by the collisional rate between gas molecules and dust grains, $t_{\rm coll}$. The angular offsets derived are equivalent to collisional timescales that range from the near instantaneous up to $t_{\rm coll}$ ~ 10 years, depending on whether a straight or a curved dust shadow, as suggested by HST observations reported by Debes et al. (2017), is adopted. The inferred range of $t_{\rm coll}$ are consistent with those predictions based on representative gas densities, temperatures, gas-to-dust ratios and grain sizes. These results represent the first time empirical constraints can be placed on $t_{\rm coll}$.

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

Molecules with ALMA at Planet-forming Scales (MAPS) III: Characteristics of Radial Chemical Substructures

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS). A Circumplanetary Disk Candidate in Molecular Line Emission in the AS 209 Disk

We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion.

preprint2022arXiv

Scattering and sublimation: a multi-scale view of $μ$m-sized dust in the inclined disc of HD 145718

We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angle, PA\,$=-1.0-0.6^{\circ}$, consistent with previous estimates. The NIR scattering surface extends out to $\sim75\,$au and we infer an aspect ratio, $h_{\rm{scat}}(r)/r\sim0.24$ in $J$-band; $\sim0.22$ in $H$-band. Our GPI images and VLTI+CHARA NIR interferometry suggest that the disc surface layers are populated by grains $\gtrsim λ/2π$ in size, indicating these grains are aerodynamically supported against settling and/or the density of smaller grains is relatively low. We demonstrate that our geometric analysis provides a reasonable assessment of the height of the NIR scattering surface at the outer edge of the disc and, if the inclination can be independently constrained, has the potential to probe the flaring exponent of the scattering surface in similarly inclined ($i\gtrsim70^{\circ}$) discs. In re-evaluating HD~145718's stellar properties, we found that the object's dimming events - previously characterised as UX~Or and dipper variability - are consistent with dust occultation by grains larger, on average, than found in the ISM. This occulting dust likely originates close to the inferred dust sublimation radius at $0.17\,$au.

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

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

Original use of MUSE's laser tomography adaptive optics to directly image young accreting exoplanets

We present recent results obtained with the VLT/MUSE Integral Field Spectrograph fed by the 4LGSF and its laser tomography adaptive optics module GALACSI. While this so-called narrow-field mode of MUSE was not designed to perform directly imaging of exoplanets and outflows, we show that it can be a game changer to detect and characterize young exoplanets with a prominent emission lines (i.e Hα, tracer of accretion), at moderate contrasts. These performances are achieved thanks to the combo of a near-diffraction limited PSF and a medium resolution spectrograph and a cross-correlation approach in post-processing . We discuss this in the context of ground and space, infrared and visible wavelengths, preparing for missions like JWST and WFIRST in great synergy and as pathfinder for future ELT/GSMT (Extremely Large and/or Giant Segmented Mirror Telescopes) instruments.

preprint2020arXiv

SCExAO/CHARIS High-Contrast Imaging of Spirals and Darkening Features in the HD 34700 A Protoplanetary Disk

We present Subaru/SCExAO+CHARIS broadband ($JHK$-band) integral field spectroscopy of HD 34700 A. CHARIS data recover HD 34700 A's disk ring and confirm multiple spirals discovered in Monnier et al. (2019). We set limits on substellar companions of $\sim12\ M_{\rm Jup}$ at $0\farcs3$ (in the ring gap) and $\sim5\ M_{\rm Jup}$ at $0\farcs75$ (outside the ring). The data reveal darkening effects on the ring and spiral, although we do not identify the origin of each feature such as shadows or physical features related to the outer spirals. Geometric albedoes converted from the surface brightness suggests a higher scale height and/or prominently abundant sub-micron dust at position angle between $\sim45^\circ$ and $90^\circ$. Spiral fitting resulted in very large pitch angles ($\sim30-50^\circ$) and a stellar flyby of HD 34700 B or infall from a possible envelope is perhaps a reasonable scenario to explain the large pitch angles.

preprint2020arXiv

Variable Accretion onto Protoplanet Host Star PDS 70

The PDS 70 system has been subject to many studies in the past year following the discovery of two accreting planets in the gap of its circumstellar disk. Nevertheless, the mass accretion rate onto the star is still not well known. Here we determined the stellar mass accretion rate and its variability based on TESS and HARPS observations. The stellar light curve shows a strong signal with a $3.03\pm0.06$ days period, which we attribute to stellar rotation. Our analysis of the HARPS spectra shows a rotational velocity of $v\sin\,i=16.0\pm0.5\,{\rm km\,s^{-1}}$, indicating that the inclination of the rotation axis is $50\pm8$ degrees. This implies that the rotation axes of the star and its circumstellar disk are parallel within the measurement error. We apply magnetospheric accretion models to fit the profiles of the H$α$ line and derive mass accretion rates onto the star in the range of $0.6-2.2\times10^{-10}\,{\rm M_{\odot}yr^{-1}}$, varying over the rotation phase. The measured accretion rates are in agreement with those estimated from NUV fluxes using accretion shock models. The derived accretion rates are higher than expected from the disk mass and planets properties for the low values of the viscous parameter $α$ suggested by recent studies, potentially pointing to an additional mass reservoir in the inner disk to feed the accretion, such as a dead zone. We find that the He I $λ$10830 line shows a blueshifted absorption feature, indicative of a wind. The mass-loss rate estimated from the line depth is consistent with an accretion-driven inner disk MHD wind.

preprint2020arXiv

Visualizing the Kinematics of Planet Formation

A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on the dynamical structure of its parental disk has progressed significantly, revealing a host of potentially observable features which would betray the presence of a deeply embedded planet. In concert, recent observations have shown that subtle perturbations in the kinematic structure of protoplanetary disks are found in multiple sources, potentially the characteristic disturbances associated with embedded planets. In this work, we review the theoretical background of planet-disk interactions, focusing on the kinematical features, and the current methodologies used to observe these interactions in spatially and spectrally resolved observations. We discuss the potential pit falls of such kinematical detections of planets, providing best-practices for imaging and analysing interferometric data, along with a set of criteria to use as a benchmark for any claimed detection of embedded planets. We finish with a discussion on the current state of simulations in regard to planet-disk interactions, highlighting areas of particular interest and future directions which will provide the most significant impact in our search for embedded planets. This work is the culmination of the 'Visualizing the Kinematics of Planet Formation' workshop, held in October 2019 at the Center for Computational Astrophysics at the Flatiron Institute in New York City.

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

Magnetospheric Accretion as a Source of H$α$ Emission from Proto-planets around PDS 70

Advances in high-resolution imaging have revealed H$α$ emission separated from the host star. It is generally believed that the emission is associated with forming planets in protoplanetary disks. However, the nature of this emission is still not fully understood. Here we report a modeling effort of H$α$ emission from the planets around the young star PDS 70. Using standard magnetospheric accretion models previously applied to accreting young stars, we find that the observed line fluxes can be reproduced using a range of parameters relevant to PDS 70b and c, with the mean mass accretion rate of log(${\rm \dot{M}}$) = $-8.0\pm0.6$ M$_{\rm Jup}$ yr$^{-1}$ and $-8.1\pm0.6$ M$_{\rm Jup}$ yr$^{-1}$ for PDS 70b and PDS 70c, respectively. Our results suggest that H$α$ emission from young planets can originate in the magnetospheric accretion of mass from the circumplanetary disk. We find that empirical relationships between mass accretion rate and H$α$ line properties frequently used in T Tauri stars are not applicable in the planetary mass regime. In particular, the correlations between line flux and mass accretion rate underpredict the accretion rate by about an order of magnitude, and the width at the 10% height of the line is insensitive to the accretion rate at ${\rm \dot{M}}$ $< 10^{-8}$ M$_{\rm Jup}$ yr$^{-1}$.

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

Self-Destructing Spiral Waves: Global Simulations of a Spiral Wave Instability in Accretion Disks

We present results from a suite of three-dimensional global hydrodynamic simulations which show that spiral density waves propagating in circumstellar disks are unstable to the growth of a parametric instability that leads to break-down of the flow into turbulence. This spiral wave instability (SWI) arises from a resonant interaction between pairs of inertial waves, or inertial-gravity waves, and the background spiral wave. The development of the instability in the linear regime involves the growth of a broad spectrum of inertial modes, with growth rates on the order of the orbital time, and results in a nonlinear saturated state in which turbulent velocity perturbations are of a similar magnitude to those induced by the spiral wave. The turbulence induces angular momentum transport, and vertical mixing, at a rate that depends locally on the amplitude of the spiral wave (we obtain a stress parameter $α\sim 5 \times 10^{-4}$ in our reference model). The instability is found to operate in a wide-range of disk models, including those with isothermal or adiabatic equations of state, and in viscous disks where the dimensionless kinematic viscosity $ν\le 10^{-5}$. This robustness suggests that the instability will have applications to a broad range of astrophysical disk-related phenomena, including those in close binary systems, planets embedded in protoplanetary disks (including Jupiter in our own Solar System) and FU Orionis outburst models. Further work is required to determine the nature of the instability, and to evaluate its observational consequences, in physically more complete disk models than we have considered in this paper.

preprint2016arXiv

The Spiral Wave Instability Induced by a Giant Planet: I. Particle Stirring in the Inner Regions of Protoplanetary Disks

We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially-alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval $0.3~R_p \leq R \leq 0.7~R_p$, where $R_p$ denotes the semi-major axis of the planetary orbit (assumed to be 5~au), the estimated vertical diffusion rate associated with the turbulence is characterized by $α_{\rm diff} \sim (0.2-1.2) \times 10^{-2}$. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr time scales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.

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.

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

The First Very Long Baseline Interferometry Image of 44 GHz Methanol Maser with the KVN and VERA Array (KaVA)

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

preprint2014arXiv

VLBI observations of bright AGN jets with KVN and VERA Array (KaVA): Evaluation of Imaging Capability

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

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.

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.