Engineering Papers⌕ Search

Engineering topics

Bae, Jaehan

Publications and source records attributed to Bae, Jaehan.

Disk Evolution Study through Imaging of Nearby Young Stars (DESTINYS): A Panchromatic View of DO Tau’s Complex Kilo-astronomical-unit 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 Very Large Telescope (VLT)/SPHERE polarimetric observations of DO Tau and combine them with archival Hubble Space Telescope (HST) scattered-light images and Atacama Large Millimeter/submillimeter Array (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 astronomical units. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand astronomical units 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.

79 ASTRONOMY AND ASTROPHYSICS↗

The Chemical Inventory of the Planet-hosting Disk PDS 70

As host to two accreting planets, PDS 70 provides a unique opportunity to probe the chemical complexity of atmosphere-forming material. We present ALMA Band 6 observations of the PDS 70 disk and report the first chemical inventory of the system. With a spatial resolution of 0."4–0."5 (~50 au), 12 species are detected, including CO isotopologs and formaldehyde, small hydrocarbons, HCN and HCO{sup +} isotopologs, and S-bearing molecules. SO and CH{sub 3}OH are not detected. All lines show a large cavity at the center of the disk, indicative of the deep gap carved by the massive planets. The radial profiles of the line emission are compared to the (sub)millimeter continuum and infrared scattered light intensity profiles. Different molecular transitions peak at different radii, revealing the complex interplay between density, temperature, and chemistry in setting molecular abundances. Column densities and optical depth profiles are derived for all detected molecules, and upper limits obtained for the nondetections. Excitation temperature is obtained for H{sub 2}CO. Deuteration and nitrogen fractionation profiles from the hydrocyanide lines show radially increasing fractionation levels. Comparison of the disk chemical inventory to grids of chemical models from the literature strongly suggests a disk molecular layer hosting a carbon-to-oxygen ratio C/O > 1, thus providing for the first time compelling evidence of planets actively accreting high C/O ratio gas at present time.

74 ATOMIC AND MOLECULAR PHYSICS↗

How Do Stars Gain Their Mass? A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star-forming Regions

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion—prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In A James Clerk Maxwell Telescope/SCUBA-2 Transient Survey of Protostars in Nearby Star-forming Regions, we are monitoring monthly with SCUBA-2 the submillimeter emission in eight fields within nearby (<500 pc) star-forming regions to measure the accretion variability of protostars. The total survey area of ∼1.6 deg{sup 2} includes ∼105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125–0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of submillimeter observations to reach a precision of ∼2%–3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any submillimeter variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for submillimeter variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

79 ASTRONOMY AND ASTROPHYSICS↗

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{sub p}⩽R⩽0.7R{sub p}, where R {sub p} denotes the semimajor axis of the planetary orbit (assumed to be 5 au), the estimated vertical diffusion rate associated with the turbulence is characterized by α{sub diff}∼ (0.2--1.2)×10{sup −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 timescales 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 that shows 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 α ∼ 5 × 10{sup −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 ν ≤ 10{sup −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.

79 ASTRONOMY AND ASTROPHYSICS↗

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{sub 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 the 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 ∼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 ≲10{sup −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 ∼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 be that protoplanetary disks have favorable conditions for vortex formation during the protostellar infall phase, which might enhance early planetary core formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Accretion outbursts in self-gravitating protoplanetary disks

We improve on our previous treatments of the 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., PdV 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 the gravitationally unstable region before they trigger outbursts at R ≲ 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 (α{sub 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 α{sub rd} triggers thermally driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when α{sub 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 M-dot {sub acc}∼10{sup −8}--10{sup −7} M{sub ⊙} yr{sup −1} over 60% of the time and spend less than 15% of the infall phase in accretion outbursts. While our models indicate that episodic mass accretion during protostellar evolution can qualitatively help explain the low accretion luminosities seen in most low-mass protostars, detailed tests of the mechanism will require model calculations for a range of protostellar masses with some constraint on the initial core angular momentum, which affects the length of time spent in a quasi-steady disk accretion phase.

79 ASTRONOMY AND ASTROPHYSICS↗

THE FIRST VERY LONG BASELINE INTERFEROMETRY IMAGE OF A 44 GHz METHANOL MASER WITH THE KVN AND VERA ARRAY (KaVA)

We have carried out the first very long baseline interferometry (VLBI) imaging of a 44 GHz class I methanol maser (7{sub 0}-6{sub 1} A {sup +}) 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 × 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than ≈ 650 km corresponding to 100 Mλ 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 × 2 mas, which corresponds to the linear size of ∼15 AU × 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ∼3.5 × 10{sup 8} to 1.0 × 10{sup 10} K, which are higher than the 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.

79 ASTRONOMY AND ASTROPHYSICS↗

THE 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 for {approx}< 20% of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that the 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 M-dot-R{sub 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 the order of 10 AU.

79 ASTRONOMY AND ASTROPHYSICS↗

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 {sup d}ead zone{sup )}. 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 {approx}< 0.5 AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

79 ASTRONOMY AND ASTROPHYSICS↗