Search

In the core accretion hypothesis, giant planets form by gas accretion onto solid protoplanetary cores. The minimum (or critical) core mass to form a gas giant is typically quoted as 10 M⊕. The actual value depends on several factors: the location in the protoplanetary disk, atmospheric opacity, and the accretion rate of solids. Motivated by ongoing direct imaging searches for giant planets, this study investigates core mass requirements in the outer disk. To determine the fastest allowed rates of gas accretion, we consider solid cores that no longer accrete planetesimals, as this would heat the gaseous envelope. Our spherical, two-layer atmospheric cooling model includes an inner convective region and an outer radiative zone that matches onto the disk. We determine the minimum core mass for a giant planet to form within a typical disk lifetime of 3 Myr. The minimum core mass declines with disk radius, from ∼8.5 M⊕ at 5 AU to ∼3.5 M⊕ at 100 AU, with standard interstellar grain opacities. Lower temperatures in the outer disk explain this trend, while variations in disk density are less influential. At all distances, a lower dust opacity or higher mean molecular weight reduces the critical core mass. Our non-self-gravitating, analytic cooling model reveals that self-gravity significantly affects early atmospheric evolution, starting when the atmosphere is only ∼10% as massive as the core. [ABSTRACT FROM AUTHOR]

Some regions in protoplanetary disks are turbulent, while some regions are quiescent (e.g. the dead zone). In order to study how planets open gaps in both inviscid hydrodynamic disk (e.g. the dead zone) and the disk subject to magnetorotational instability (MRI), we carried out both shearing box two-dimensional inviscid hydrodynamical simulations and three-dimensional unstratified magnetohydrodynamical (MHD) simulations (having net vertical magnetic fields) with a planet at the box center. We found that, due to the nonlinear wave steepening, even a low mass planet can open gaps in both cases, in contradiction to the "thermal criterion" for gap opening. In order to understand if we can represent the MRI turbulent stress with the viscous α prescription for studying gap opening, we compare 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. 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. We also found that long-lived density features (termed zonal flows) produced by the MRI can affect planet migration. 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. [ABSTRACT FROM AUTHOR]

Carbon monoxide is a commonly used IR/submillimeter tracer of gas in protoplanetary disks. We present an analysis of ultraviolet CO emission in Hubble Space Telescope Cosmic Origins Spectrograph spectra for 12 Classical T Tauri stars (CTTSs). Several ro-vibrational bands of the CO A¹Π-X¹Σ+ (Fourth Positive) electronic transition system are spectrally resolved from emission of other atoms and H2. The CO A¹Π v' = 14 state is populated by absorption of Lyα photons, created at the accretion column on the stellar surface. For targets with strong CO emission, we model the Lyα radiation field as an input for a simple fluorescence model to estimate CO rotational excitation temperatures and column densities. Typical column densities range from NCO = 1018 to 1019 cm-2. Our measured excitation temperatures are mostly below TCO= 600 K, cooler than typical M-band CO emission. These temperatures and the emission line widths imply that the UV emission originates in a different population of CO than that which is IR-emitting. We also find a significant correlation between CO emission and the disk accretion rate Ṁ and age. Our analysis shows that ultraviolet CO emission can be a useful diagnostic of CTTS disk gas. [ABSTRACT FROM AUTHOR]

We present high-resolution (0".3 = 16 AU), high signal-to-noise ratio Submillimeter Array observations of the 870 µm (345 GHz) continuum and CO J = 3 - 2 line emission from the protoplanetary disk around TW Hya. Using continuum and line radiative transfer calculations, these data and the multiwavelength spectral energy distribution are analyzed together in the context of simple two-dimensional parametric disk structure models. Under the assumptions of a radially invariant dust population and gas-to-dust mass ratio, we are unable to simultaneously reproduce the CO and dust observations with model structures that employ either a single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared scattered light) that extends to a much larger radius of at least 215 AU. We discuss some possible explanations for the observed radial distribution of millimeter-sized dust grains and the apparent CO-dust size discrepancy, and suggest that they may be hallmarks of substructure in the dust disk or natural signatures of the growth and radial drift of solids that might be expected for disks around older pre-main-sequence stars like TW Hya. [ABSTRACT FROM AUTHOR]

We study the turbulence induced in the dust layer of a protoplanetary disk based on the energetics of dust accretion due to gas drag. We estimate turbulence strength from the energy supplied by dust accretion, using the radial drift velocity of the dust particles in a laminar disk. Our estimate of the turbulence strength agrees with previous analytical and numerical research on the turbulence induced by Kelvin-Helmholtz and/or streaming instabilities for particles with stopping time less than the Keplerian time. For such small particles, the strongest turbulence is expected to occur when the dust-to-gas ratio of the disk is ~Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed.(hg/r) ~ 10-2, where Ceff ≈ 0.2 represents the energy supply efficiency to turbulence and hg/r ~ 5 x 10-2 is the aspect ratio of the gas disk. The maximum viscosity parameter is αmax ~ CeffTs(hg/r)² ~ 10-4Ts, where Ts(< 1) is the non-dimensional stopping time of the dust particles. Modification to the dust-to-gas ratio from the standard value, 10-2, by any process, results in weaker turbulence and a thinner dust layer, and consequently may accelerate the growth process of the dust particles. [ABSTRACT FROM AUTHOR]

The direct study of molecular gas in inner protoplanetary disks is complicated by uncertainties in the spatial distribution of the gas, the time variability of the source, and the comparison of observations across a wide range of wavelengths. Some of these challenges can be mitigated with far-ultraviolet spectroscopy. Using new observations obtained with the Hubble Space Telescope Cosmic Origins Spectrograph, we measure column densities and rovibrational temperatures for CO and H2 observed on the line of sight through the AA Tauri circumstellar disk. CO A - X absorption bands are observed against the far-UV continuum. The CO absorption is characterized by log10 (N(12CO)) = 17.5 ± 0.5 cm-2 and Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed K, although this rotational temperature may underestimate the local kinetic temperature of the CO-bearing gas. We also detect 13CO in absorption with an isotopic ratio of ~20. We do not observe H2 absorption against the continuum; however, hot H2 (v > 0) is detected in absorption against the Lyα emission line. We measure the column densities in eight individual rovibrational states, determining a total log10 (N(H2)) = Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed with a thermal temperature of Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed K. The high temperature of the molecules, the relatively small H2 column density, and the high inclination of the AA Tauri disk suggest that the absorbing gas resides in an inner disk atmosphere. If the H2 and CO are cospatial within a molecular layer ~0.6 AU thick, this region is characterized by Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed with an observed (CO/H2) ratio of ~0.4. We also find evidence for a departure from a purely thermal H2 distribution, suggesting that excitation by continuum photons and H2 formation may be altering the level populations in the molecular gas. [ABSTRACT FROM AUTHOR]

Tiny grains such as polycyclic aromatic hydrocarbons (PAHs) have been thought to dramatically reduce the coupling between the gas and magnetic fields in weakly ionized gas such as in protoplanetary disks (PPDs) because they provide a tremendous surface area to recombine free electrons. The presence of tiny grains in PPDs thus raises the question of whether the magnetorotational instability (MRI) is able to drive rapid accretion consistent with observations. Charged tiny grains have similar conduction properties as ions, whose presence leads to qualitatively new behaviors in the conductivity tensor, characterized by n̄/ne where n̄e and n̄ denote the number densities of free electrons and all other charged species, respectively. In particular, Ohmic conductivity becomes dominated by charged grains rather than by electrons when n̄ /ne exceeds about 10³, and Hall and ambipolar diffusion (AD) coefficients are reduced by a factor of (n̄/ne)² in the AD-dominated regime relative to that in the Ohmic regime. Applying the methodology of Bai, we find that in PPDs, when PAHs are sufficiently abundant (≳10-9 per H2 molecule), there exists a transition radius rtrans of about 10-20 AU, beyond which the MRI active layer extends to the disk midplane. At r < rtrans, the optimistically predicted MRI-driven accretion rate M is one to two orders of magnitude smaller than that in the grain-free case, which is too small compared with the observed rates, but is in general no smaller than the predicted M with solar-abundance 0.1 µm grains. At r > rtrans, we find that, remarkably, the predicted M exceeds the grain-free case due to a net reduction of AD by charged tiny grains and reaches a few times 10-1 M⊙ yr-1. This is sufficient to account for the observed M in transitional disks. Larger grains (≳0.1 µm) are too massive to reach such high abundance as tiny grains and to facilitate the accretion process. [ABSTRACT FROM AUTHOR]

Non-ideal MHD effects play an important role in the gas dynamics in protoplanetary disks (PPDs). This paper addresses the influence of non-ideal MHD effects on the magnetorotational instability (MRI) and angular momentum transport in PPDs using the most up-to-date results from numerical simulations. We perform chemistry calculations using a complex reaction network with standard prescriptions for X-ray and cosmic-ray ionizations. We first show that whether or not grains are included, the recombination time is at least one order of magnitude less than the orbital time within five disk scale heights, justifying the validity of local ionization equilibrium and strong couplin limit in PPDs. The full conductivity tensor at different disk radii and heights is evaluated, with the MRI active region determined by requiring that (1) the Ohmic Elsasser number A be greater than 1 and (2) the ratio of gas to magnetic pressure β be greater than βmin (Am) as identified in the recent study by Bai & Stone, where Am is the Elsasser number for ambipolar diffusion. With full flexibility as to the magnetic field strength, we provide a general framework for estimating the MRI-driven accretion rate M and the magnetic field strength in the MRI active layer. We find that the MRI active layer always exists.at any disk radius as long as the magnetic field in PPDs is sufficiently weak. However, the optimistically predicted M in the inner disk (r = 1-10 AU) appears insufficient to account for the observed range of accretion rates in PPDs (around 10-8 M⊙ yr-1) even in the grain-free calculation, and the presence of solar abundance sub-micron grains further reduces M by one to two orders of magnitude. Moreover, we find that the predicted M increases with radius in the inner disk where accretion is layered, which would lead to runaway mass accumulation if disk accretion is solely driven by the MRI. Our results suggest that stronger sources of ionization and/or additional mechanisms such as magnetized wind are needed to explain the observed accretion rates in PPDs. In contrast, our predicted M is on the order of 10-9 M⊙ yr-1 in the outer disk, consistent with the observed accretion rates in transitional disks. [ABSTRACT FROM AUTHOR]

In the outer part of a protoplanetary disk, dust grains consist of silicate core covered by an ice mantle. A temporal heating event in the disk results in sublimation of the ice mantle. After the end of the heating event, as the temperature decreases, H2O molecules recondense on the surface of the dust grain. Ultimately, the dust grain is covered by an ice mantle. Because the equilibrium vapor pressure on the grain surface decreases with the grain size, a large grain grows faster than a small grain. As a result, the size of an icy dust grain changes as a result of the heating event. The change in size also affects the mechanical properties of the dust aggregates formed by the icy grains. In this paper, we investigated the evolution of the size distribution of icy dust grains during sublimation and condensation. We found that the size evolution of icy grains can be divided into two stages. In the first stage, the icy grains grow through condensation of H2O molecules. In the second stage, the size of grains changes further as H2O molecules are transferred between icy grains while the surrounding gas condenses. The size distribution of the icy dust grains becomes bimodal, with a small number of relatively large grains and many small grains without an icy mantle. Possible effects of the size change on the evolution of icy dust aggregates are discussed. [ABSTRACT FROM AUTHOR]

We consider the effects of self-gravity on the hydrostatic balance in the vertical direction of a gaseous disk and discuss the possible signature of the self-gravity that may be captured by direct imaging observations of protoplanetary disks in the future. In this paper, we consider a vertically isothermal disk in order to isolate the effects of self-gravity. The specific disk model we consider in this paper is the one with a radial surface density gap, at which the Toomre's Q-parameter of the disk varies rapidly in the radial direction. We calculate the vertical structure of the disk including the effects of self-gravity. We then calculate the scattered light and the dust thermal emission. We find that if the disk is massive enough and the effects of self-gravity come into play, a weak bump-like structure at the gap edge appears in the near-infrared (NIR) scattered light, while no such bump-like structure is seen in the submillimeter (sub-mm) dust continuum image. The appearance of the bump is caused by the variation of the height of the surface in the NIR wavelength. If such a bump-like feature is detected in future direct imaging observations, combined with sub-mm observations, it will give us useful information about the physical states of the disk. [ABSTRACT FROM AUTHOR]

Evolution of a snow line in an optically thick protoplanetary disk is investigated with numerical simulations. The ice-condensing region in the disk is obtained by calculating the temperature and the density with the 1+1D approach. The snow line migrates as the mass accretion rate (M) in the disk decreases with time. Calculations are carried out from an early phase with high disk accretion rates (M ~ 10-7 M⊙ yr-l) to a later phase with low disk accretion rates (M ~ 10-12 M⊙ yr-1) using the same numerical method. It is found that the snow line moves inward for M ≤ 10-10 M⊙ yr-1, while it gradually moves outward in the later evolution phase with M ≤ 10-10 Mo yr-1. In addition to the silicate opacity, the ice opacity is taken into consideration. In the inward migration phase, the additional ice opacity increases the distance of the snow line from the central star by a factor of 1.3 for dust grains ≤ 10µm in size and of 1.6 for ≥ 100 µm. It is inevitable that the snow line comes inside Earth's orbit in the course of the disk evolution if the viscosity parameter α is in the range 0.001-0.1, the dust-to-gas mass ratio is higher than a tenth of the solar abundance value, and the dust grains are smaller than 1 mm. The formation of water-devoid planetesimals in the terrestrial planet region seems to be difficult throughout the disk evolution, which imposes a new challenge to planet formation theory. [ABSTRACT FROM AUTHOR]

Although there has yet been no undisputed discovery of a still-forming planet embedded in a gaseous protoplanetary disk, the cleared inner holes of transitional disks may be signposts of young planets. Here, we show that the subset of accreting transitional disks with wide, optically thin inner holes of 15 AU or more can only be sculpted by multiple planets orbiting inside each hole. Multiplanet systems provide two key ingredients for explaining the origins of transitional disks. First, multiple planets can clear wide inner holes where single planets open only narrow gaps. Second, the confined, non-axisymmetric accretion flows produced by multiple planets provide a way for an arbitrary amount of mass transfer to occur through an apparently optically thin hole without overproducing infrared excess flux. Rather than assuming that the gas and dust in the hole are evenly and axisymmetrically distributed, one can construct an inner hole with apparently optically thin infrared fluxes by covering a macroscopic fraction of the hole's surface area with locally optically thick tidal tails. We also establish that other clearing mechanisms, such as photoevaporation, cannot explain our subset of accreting transitional disks with wide holes. Transitional disks are therefore high-value targets for observational searches for young planetary systems. [ABSTRACT FROM AUTHOR]

Collisional growth of dust aggregates is a plausible root of planetesimals forming in protoplanetary disks. However, a rebound of colliding dust aggregates prevents dust from growing into planetesimals. In fact, rebounding aggregates are observed in laboratory experiments but not in previous numerical simulations. Therefore, the condition of rebound between dust aggregates should be clarified to better understand the processes of dust growth and planetesimal formation. We have carried out numerical simulations of aggregate collisions for various types of aggregates and succeeded in reproducing a rebound of colliding aggregates under specific conditions. Our finding is that in the rebound process, the key factor of the aggregate structure is the coordination number, namely, the number of particles in contact with a particle. A rebound is governed by the energy dissipation along with restructuring of the aggregates and a large coordination number inhibits the restructuring at collisions. Results of our numerical simulation for various aggregates indicate that they stick to each other when the mean coordination number is less than 6, regardless of their materials and structures, as long as their collision velocity is less than the critical velocity for fragmentation. This criterion of the coordination number would correspond to a filling factor of ~0.3, which is somewhat larger than that reported in laboratory experiments. In protoplanetary disks, dust aggregates are expected to have low bulk densities (<0.1 g cm-3) during their growth, which would prevent dust aggregates from rebounding. This result supports the formation of planetesimals with direct dust growth in protoplanetary disks. [ABSTRACT FROM AUTHOR]

Protoplanetary disk simulations show that a single planet can excite more than one spiral arm, possibly explaining the recent observations of multiple spiral arms in some systems. In this paper, we explain the mechanism by which a planet excites multiple spiral arms in a protoplanetary disk. Contrary to previous speculations, the formation of both primary and additional arms can be understood as a linear process when the planet mass is sufficiently small. A planet resonantly interacts with epicyclic oscillations in the disk, launching spiral wave modes around the Lindblad resonances. When a set of wave modes is in phase, they can constructively interfere with each other and create a spiral arm. More than one spiral arm can form because such constructive interference can occur for different sets of wave modes, with the exact number and launching position of the spiral arms being dependent on the planet mass as well as the disk temperature profile. Nonlinear effects become increasingly important as the planet mass increases, resulting in spiral arms with stronger shocks and thus larger pitch angles. This is found to be common for both primary and additional arms. When a planet has a sufficiently large mass (≳3 thermal masses for (h/r)p = 0.1), only two spiral arms form interior to its orbit. The wave modes that would form a tertiary arm for smaller mass planets merge with the primary arm. Improvements in our understanding of the formation of spiral arms can provide crucial insights into the origin of observed spiral arms in protoplanetary disks. [ABSTRACT FROM AUTHOR]

We revisit calculations of nebular hydrogen Lyα and He ii λ1640 line strengths for Population III (Pop III) galaxies, undergoing continuous, and bursts of, star formation. We focus on initial mass functions (IMFs) motivated by recent theoretical studies, which generally span a lower range of stellar masses than earlier works. We also account for case-B departures and the stochastic sampling of the IMF. In agreement with previous work, we find that departures from case-B can enhance the Lyα flux by a factor of a few, but we argue that this enhancement is driven mainly by collisional excitation and ionization, and not due to photoionization from the n = 2 state of atomic hydrogen. The increased sensitivity of the Lyα flux to the high-energy end of the galaxy spectrum makes it more subject to stochastic sampling of the IMF. The latter introduces a dispersion in the predicted nebular line fluxes around the deterministic value by as much as a factor of ∼4. In contrast, the stochastic sampling of the IMF has less impact on the emerging Lyman Werner photon flux. When case-B departures and stochasticity effects are combined, nebular line emission from Pop III galaxies can be up to one order of magnitude brighter than predicted by “standard” calculations that do not include these effects. This enhances the prospects for detection with future facilities such as the James Webb Space Telescope and large, ground-based telescopes. [ABSTRACT FROM AUTHOR]

We present the first high-resolution sub-millimeter survey of both dust and gas for a large population of protoplanetary disks. Characterizing fundamental properties of protoplanetary disks on a statistical level is critical to understanding how disks evolve into the diverse exoplanet population. We use the Atacama Large Millimeter/Submillimeter Array (ALMA) to survey 89 protoplanetary disks around stars with in the young (1–3 Myr), nearby (150–200 pc) Lupus complex. Our observations cover the 890 μm continuum and the 13CO and C18O 3–2 lines. We use the sub-millimeter continuum to constrain to a few Martian masses (0.2–0.4 M⊕) and the CO isotopologue lines to constrain to roughly a Jupiter mass (assuming an interstellar medium (ISM)-like abundance). Of 89 sources, we detect 62 in continuum, 36 in 13CO, and 11 in C18O at significance. Stacking individually undetected sources limits their average dust mass to Lunar masses (0.03 M⊕), indicating rapid evolution once disk clearing begins. We find a positive correlation between and M*, and present the first evidence for a positive correlation between and M*, which may explain the dependence of giant planet frequency on host star mass. The mean dust mass in Lupus is 3× higher than in Upper Sco, while the dust mass distributions in Lupus and Taurus are statistically indistinguishable. Most detected disks have and gas-to-dust ratios , assuming an ISM-like abundance; unless CO is very depleted, the inferred gas depletion indicates that planet formation is well underway by a few Myr and may explain the unexpected prevalence of super-Earths in the exoplanet population. [ABSTRACT FROM AUTHOR]

To what degree the cluster environment influences the sizes of protoplanetary disks surrounding young stars is still an open question. This is particularly true for the short-lived clusters typical for the solar neighborhood, in which the stellar density and therefore the influence of the cluster environment change considerably over the first 10 Myr. In previous studies, the effect of the gas on the cluster dynamics has often been neglected; this is remedied here. Using the code NBody6++, we study the stellar dynamics in different developmental phases—embedded, expulsion, and expansion—including the gas, and quantify the effect of fly-bys on the disk size. We concentrate on massive clusters (Mcl ≥ 103–6 ∗ 104MSun), which are representative for clusters like the Orion Nebula Cluster (ONC) or NGC 6611. We find that not only the stellar density but also the duration of the embedded phase matters. The densest clusters react fastest to the gas expulsion and drop quickly in density, here 98% of relevant encounters happen before gas expulsion. By contrast, disks in sparser clusters are initially less affected, but because these clusters expand more slowly, 13% of disks are truncated after gas expulsion. For ONC-like clusters, we find that disks larger than 500 au are usually affected by the environment, which corresponds to the observation that 200 au-sized disks are common. For NGC 6611-like clusters, disk sizes are cut-down on average to roughly 100 au. A testable hypothesis would be that the disks in the center of NGC 6611 should be on average ≈20 au and therefore considerably smaller than those in the ONC. [ABSTRACT FROM AUTHOR]

The presence of excesses of short-lived radionuclides in the early solar system evidenced in meteorites has been taken as testament to close encounters with exotic nucleosynthetic sources, including supernovae or AGB stars. An analysis of the likelihoods associated with different sources of these extinct nuclides in the early solar system indicates that, rather than being exotic, their abundances were typical of star-forming regions like those observed today in the Galaxy. The radiochemistry of the early solar system is therefore unexceptional, being the consequence of extensive averaging of solids from molecular clouds. [ABSTRACT FROM AUTHOR]

Protoplanetary disks with non-axisymmetric structures have been observed. The Rossby wave instability (RWI) is considered as one of the origins of the non-axisymmetric structures. We perform linear stability analyses of the RWI in barotropic flow using four representative types of the background flow on a wide parameter space. We find that the co-rotation radius is located at the background vortensity minimum with large concavity if the system is marginally stable to the RWI, and this allows us to easily check the stability against the RWI. We newly derive the necessary and sufficient condition for the onset of the RWI in semi-analytic form. We discuss the applicability of the new condition in realistic systems and the physical nature of the RWI. [ABSTRACT FROM AUTHOR]

In protoplanetary disks, the distribution and abundance of small (sub)micron grains are important for a range of physical and chemical processes. For example, they dominate the optical depth at short wavelengths and their surfaces are the sites of many important chemical reactions, such as the formation of water. Based on their aerodynamical properties (i.e., their strong dynamical coupling with the surrounding gas) it is often assumed that these small grains are well-mixed with the gas. Our goal is to study the vertical (re)distribution of grains taking into account settling, turbulent diffusion, and collisions with other dust grains. Assuming a fragmentation-limited background dust population, we developed a Monte Carlo approach that follows single monomers as they move through a vertical column of gas and become incorporated in different aggregates as they undergo sticking and fragmenting collisions. We find that (sub)micron grains are not necessarily well-mixed vertically, but can become trapped in a thin layer with a scale height that is significantly smaller than that of the gas. This collisional trapping occurs when the timescale for diffusion is comparable to or longer than the collision timescale in the midplane and its effect is strongest when the most massive particles in the size distribution show significant settling. Based on simulations and analytical considerations, we conclude that for typical dust-to-gas ratios and turbulence levels, the collisional trapping of small grains should be a relatively common phenomenon. The absence of trapping could then indicate a low dust-to-gas ratio, possibly because a large portion of the dust mass has been removed through radial drift or is locked up in planetesimals. [ABSTRACT FROM AUTHOR]

The latest observation of HL Tau by ALMA revealed spectacular concentric dust rings in its circumstellar disk. We attempt to explain the multiple ring structure as a consequence of aggregate sintering. Sintering is known to reduce the sticking efficiency of dust aggregates and occurs at temperatures slightly below the sublimation point of the constituent material. We present a dust growth model that incorporates sintering and use it to simulate global dust evolution due to sintering, coagulation, fragmentation, and radial inward drift in a modeled HL Tau disk. We show that aggregates consisting of multiple species of volatile ices experience sintering, collisionally disrupt, and pile up at multiple locations slightly outside the snow lines of the volatiles. At wavelengths of 0.87–1.3 mm, these sintering zones appear as bright, optically thick rings with a spectral slope of , whereas the non-sintering zones appear as darker, optically thinner rings of a spectral slope of –2.5. The observational features of the sintering and non-sintering zones are consistent with those of the major bright and dark rings found in the HL Tau disk, respectively. Radial pileup and vertical settling occur simultaneously if disk turbulence is weak and if monomers constituting the aggregates are in radius. For the radial gas temperature profile of , our model perfectly reproduces the brightness temperatures of the optically thick bright rings and reproduces their orbital distances to an accuracy of . [ABSTRACT FROM AUTHOR]

A global picture of the evolution of protoplanetary disks (PPDs) is key to understanding almost every aspect of planet formation, where standard α-disk models have been continually employed for their simplicity. In the meantime, disk mass loss has been conventionally attributed to photoevaporation, which controls disk dispersal. However, a paradigm shift toward accretion driven by magnetized disk winds has taken place in recent years, thanks to studies of non-ideal magnetohydrodynamic effects in PPDs. I present a framework of global PPD evolution aiming to incorporate these advances, highlighting the role of wind-driven accretion and wind mass loss. Disk evolution is found to be largely dominated by wind-driven processes, and viscous spreading is suppressed. The timescale of disk evolution is controlled primarily by the amount of external magnetic flux threading the disks, and how rapidly the disk loses the flux. Rapid disk dispersal can be achieved if the disk is able to hold most of its magnetic flux during the evolution. In addition, because wind launching requires a sufficient level of ionization at the disk surface (mainly via external far-UV (FUV) radiation), wind kinematics is also affected by the FUV penetration depth and disk geometry. For a typical disk lifetime of a few million years, the disk loses approximately the same amount of mass through the wind as through accretion onto the protostar, and most of the wind mass loss proceeds from the outer disk via a slow wind. Fractional wind mass loss increases with increasing disk lifetime. Significant wind mass loss likely substantially enhances the dust-to-gas mass ratio and promotes planet formation. [ABSTRACT FROM AUTHOR]

We irradiated a (4:1) mixture of water and carbon monoxide with soft X-rays of energies up to 1.2 keV. The experiments were performed using the spherical grating monochromator beamline at National Synchrotron Radiation Research Center in Taiwan. Both monochromatic (300 and 900 eV) and broader energy fluxes (250–1200 eV) were employed. During the irradiation, the H2O + CO mixture was ionized, excited, and fragmented, producing a number of reactive species. The composition of the ice has been monitored throughout both the irradiation and warm-up phases. We identified several products, which can be related through a plausible chemical reaction scheme. Such chemistry is initiated by the injection of energetic photoelectrons that produce multiple ionization events generating a secondary electron cascade. The results have been discussed in light of a model for protoplanetary disks around young solar-type stars. [ABSTRACT FROM AUTHOR]

Protoplanetary disks originate in the collapse of molecular cloud cores. The formation and evolution of disks are influenced by the properties of molecular cloud cores. In this paper we investigate the dependence of disk lifetimes and accretion rates on cloud core properties. We find that the lifetime increases as the angular velocities and the mass of cloud cores increase and that the lifetime decreases as the core temperature increases. We have calculated the distribution of disk lifetimes and disk fractions with stellar age. Our calculations show that the lifetime is in the range of 1–15 Myr and that the typical lifetime is 1–3 Myr. There are a few disks with lifetimes greater than 10 Myr and of the disks have lifetimes less than 1 Myr. We also fit the disk fraction by an exponential decay curve with characteristic time ∼3.7 Myr. Our results explain the observations of disk lifetimes. We also find that the accretion rate does not change significantly with ω and generally decreases with . At the early evolution of the disks, the relation is about . Since the effects of the photoevaporation are weak at this stage, this relation is the consequence of the cloud core properties. At the late evolution of the disks, the relation is about . For low accretion rates at this stage, the relation results from the effects of X-ray photoevaporation. The calculated relations are consistent with the observations. [ABSTRACT FROM AUTHOR]

We investigated whether outward planetesimal-driven migration (PDM) takes place or not in simulations when the self-gravity of planetesimals is included. We performed N-body simulations of planetesimal disks with a large width (0.7–4 au) that ranges over the ice line. The simulations consisted of two stages. The first-stage simulations were carried out to see the runaway growth phase using the planetesimals of initially the same mass. The runaway growth took place both at the inner edge of the disk and at the region just outside the ice line. This result was utilized for the initial setup of the second-stage simulations, in which the runaway bodies just outside the ice line were replaced by the protoplanets with about the isolation mass. In the second-stage simulations, the outward migration of the protoplanet was followed by the stopping of the migration due to the increase of the random velocity of the planetesimals. Owing to this increase of random velocities, one of the PDM criteria derived in Minton & Levison was broken. In the current simulations, the effect of the gas disk is not considered. It is likely that the gas disk plays an important role in PDM, and we plan to study its effect in future papers. [ABSTRACT FROM AUTHOR]

Weakly ionized protoplanetary disks (PPDs) are subject to nonideal magnetohydrodynamic (MHD) effects, including ohmic resistivity, the Hall effect, and ambipolar diffusion (AD), and the resulting magnetic diffusivities (, and ) largely control the disk gas dynamics. The presence of grains not only strongly reduces the disk ionization fraction, but also modifies the scalings of and with magnetic field strength. We analytically derive asymptotic expressions of and in both the strong and weak field limits and show that toward a strong field, can change sign (at a threshold field strength ), mimicking a flip of field polarity, and AD is substantially reduced. Applied to PPDs, we find that when small ∼0.1 (0.01)μm grains are sufficiently abundant (mass ratio ∼0.01 (10−4)), can change sign up to ∼2–3 scale heights above the midplane at a modest field strength (plasma ) over a wide range of disk radii. The reduction of AD is also substantial toward the AD-dominated outer disk and may activate the magnetorotational instability. We further perform local nonideal MHD simulations of the inner disk (within 10 au) and show that, with sufficiently abundant small grains, the magnetic field amplification due to the Hall-shear instability saturates at a very low level near the threshold field strength . Together with previous studies, we conclude by discussing the grain-abundance-dependent phenomenology of PPD gas dynamics. [ABSTRACT FROM AUTHOR]

The global evolution and dispersal of protoplanetary disks (PPDs) are governed by disk angular-momentum transport and mass-loss processes. Recent numerical studies suggest that angular-momentum transport in the inner region of PPDs is largely driven by magnetized disk wind, yet the wind mass-loss rate remains unconstrained. On the other hand, disk mass loss has conventionally been attributed to photoevaporation, where external heating on the disk surface drives a thermal wind. We unify the two scenarios by developing a one-dimensional model of magnetized disk winds with a simple treatment of thermodynamics as a proxy for external heating. The wind properties largely depend on (1) the magnetic field strength at the wind base, characterized by the poloidal Alfvén speed vAp, (2) the sound speed cs near the wind base, and (3) how rapidly poloidal field lines diverge (achieve scaling). When , corotation is enforced near the wind base, resulting in centrifugal acceleration. Otherwise, the wind is accelerated mainly by the pressure of the toroidal magnetic field. In both cases, the dominant role played by magnetic forces likely yields wind outflow rates that exceed purely hydrodynamical mechanisms. For typical PPD accretion-rate and wind-launching conditions, we expect vAp to be comparable to cs at the wind base. The resulting wind is heavily loaded, with a total wind mass-loss rate likely reaching a considerable fraction of the wind-driven accretion rate. Implications for modeling global disk evolution and planet formation are also discussed. [ABSTRACT FROM AUTHOR]

Recent discoveries of circumbinary planets by the Kepler mission provide motivation for understanding their birthplaces—protoplanetary disks around stellar binaries with separations . We explore properties and evolution of such circumbinary disks focusing on modification of their structure caused by tidal coupling to the binary. We develop a set of analytical scaling relations describing viscous evolution of the disk properties, which are verified and calibrated using 1D numerical calculations with realistic inputs. Injection of angular momentum by the central binary suppresses mass accretion onto the binary and causes radial distribution of the viscous angular momentum flux to be different from that in a standard accretion disk around a single star with no torque at the center. Disks with no mass accretion at the center develop an profile that is flat in radius. Radial profiles of temperature and surface density are also quite different from those in disks around single stars. Damping of the density waves driven by the binary and viscous dissipation dominates heating of the inner disk (within 1–2 AU), pushing the ice line beyond 3–5 AU, depending on disk mass and age. Irradiation by the binary governs disk thermodynamics beyond ∼10 AU. However, self-shadowing by the hot inner disk may render central illumination irrelevant out to ∼20 AU. Spectral energy distribution of a circumbinary disk exhibits a distinctive bump around 10 μm, which may facilitate identification of such disks around unresolved binaries. Efficient tidal coupling to the disk drives orbital inspiral of the binary and may cause low-mass and relatively compact binaries to merge into a single star within the disk lifetime. We generally find that circumbinary disks present favorable sites for planet formation (despite their wider zone of volatile depletion), in agreement with the statistics of Kepler circumbinary planets. [ABSTRACT FROM AUTHOR]

The Kepler mission has discovered more than 4000 exoplanet candidates. Many of them are in systems with tightly packed inner planets. Inside-out planet formation (IOPF) has been proposed as a scenario 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 pebble ring and then a planet. The planet is expected to grow in mass until it opens a gap, which helps to both truncate pebble accretion and also induce DZIB retreat that sets the location of formation of the next planet. This simple scenario may be modified if the planet undergoes significant migration from its formation location. Thus, planet–disk interactions play a crucial role in the IOPF scenario. Here we present numerical simulations that first assess the degree of migration for planets of various masses that are forming at the DZIB of an active accretion disk, where the effective viscosity is undergoing a rapid increase in the radially inward direction. We find that torques exerted on the planet by the disk 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. [ABSTRACT FROM AUTHOR]

Excitation and propagation of waves in a thermally stratified disk with an arbitrary vertical temperature profile are studied. Previous analytical studies of three-dimensional waves had been focused on either isothermal or polytropic vertical disk structures. However, at the location in a protoplanetary disk where the dominant heating source is stellar irradiation, the temperature gradient may become positive in the vertical direction. We extend the analysis to study the effects of the vertical temperature structure on the waves that are excited at the Lindblad resonances. For a hotter disk atmosphere, the f-mode contributes less to the torque and remains confined near the midplane as it propagates away from the resonances. On the other hand, the excitation of the g-modes is stronger. As they propagate, they channel to the top of disk atmosphere and their group velocities decrease. The differences compared to previous studies may have implications in understanding the wave dynamics in a realistic disk structure. [ABSTRACT FROM AUTHOR]

The growth of small planetesimals into large planetary embryos occurs far before the dispersal of the gas from the protoplanetary disk. The planetesimal–gaseous disk interactions give rise to migration and orbital evolution of the planetesimals/planets. Small planetesimals are dominated by aerodynamic gas drag. Large protoplanets, are dominated by type I migration differential torque. There is an additional mass range of intermediate-mass planetesimals (IMPs), where gravitational interactions with the disk dominate over aerodynamic gas drag, but for which such interactions were typically neglected. Here, we model these interactions using the gas dynamical friction (GDF) approach, previously used to study the disk–planet interactions at the planetary mass range. We find the critical size where GDF dominates over gas drag, and then we study the implications of GDF on single IMPs. We find that planetesimals with small inclinations rapidly become co-planar. Eccentric orbits circularize within a few Myr, provided the the planetesimal mass is large, , and that the initial eccentricity is low, Planetesimals of higher masses, , inspiral on a timescale of a few Myr, leading to an embryonic migration to the inner disk. This can lead to an overabundance of rocky material (in the form of IMPs) in the inner protoplanetary disk () and induce rapid planetary growth. This can explain the origin of super-Earth planets. In addition, GDF damps the velocities of IMPs, thereby cooling the planetesimal disk and affecting its collisional evolution through quenching the effects of viscous stirring by the large bodies. [ABSTRACT FROM AUTHOR]

The vertical shear instability (VSI) offers a potential hydrodynamic mechanism for angular momentum transport in protoplanetary disks (PPDs). The VSI is driven by a weak vertical gradient in the disk’s orbital motion, but must overcome vertical buoyancy, a strongly stabilizing influence in cold disks, where heating is dominated by external irradiation. Rapid radiative cooling reduces the effective buoyancy and allows the VSI to operate. We quantify the cooling timescale tc needed for efficient VSI growth, through a linear analysis of the VSI with cooling in vertically global, radially local disk models. We find the VSI is most vigorous for rapid cooling with in terms of the Keplerian orbital frequency, , the disk’s aspect-ratio, , the radial power-law temperature gradient, q, and the adiabatic index, γ. For longer tc, the VSI is much less effective because growth slows and shifts to smaller length scales, which are more prone to viscous or turbulent decay. We apply our results to PPD models where tc is determined by the opacity of dust grains. We find that the VSI is most effective at intermediate radii, from ∼5 to ∼50 AU with a characteristic growth time of ∼30 local orbital periods. Growth is suppressed by long cooling times both in the opaque inner disk and the optically thin outer disk. Reducing the dust opacity by a factor of 10 increases cooling times enough to quench the VSI at all disk radii. Thus the formation of solid protoplanets, a sink for dust grains, can impede the VSI. [ABSTRACT FROM AUTHOR]

The intermediate phases of planet formation are not directly observable due to lack of emission from planetesimals. Planet formation is, however, a dynamically active process resulting in collisions between the evolving planetesimals and the production of dust. Thus, indirect observation of planet formation may indeed be possible in the near future. In this paper we present synthetic observations based on numerical N-body simulations of the intermediate phase of planet formation including a state-of-the-art collision model, EDACM, which allows multiple collision outcomes, such as accretion, erosion, and bouncing events. We show that the formation of planetary embryos may be indirectly observable by a fully functioning ALMA telescope if the surface area involved in planetesimal evolution is sufficiently large and/or the amount of dust produced in the collisions is sufficiently high in mass. [ABSTRACT FROM AUTHOR]

A recent Atacama Large Millimeter/Submillimeter Array image revealed several concentric gaps in the protoplanetary disk surrounding the young star HL Tau. We consider the hypothesis that these gaps are carved by planets, and present a general framework for understanding the dynamical stability of such systems over typical disk lifetimes, providing estimates for the maximum planetary masses. We collect these easily evaluated constraints into a workflow that can help guide the design and interpretation of new observational campaigns and numerical simulations of gap opening in such systems. We argue that the locations of resonances should be significantly shifted in massive disks like HL Tau, and that theoretical uncertainties in the exact offset, together with observational errors, imply a large uncertainty in the dynamical state and stability in such disks. This presents an important barrier to using systems like HL Tau as a proxy for the initial conditions following planet formation. An important observational avenue to breaking this degeneracy is to search for eccentric gaps, which could implicate resonantly interacting planets. Unfortunately, massive disks like HL Tau should induce swift pericenter precession that would smear out any such eccentric features of planetary origin. This motivates pushing toward more typical, less massive disks. For a nominal non-resonant model of the HL Tau system with five planets, we find a maximum mass for the outer three bodies of approximately 2 Neptune masses. In a resonant configuration, these planets can reach at least the mass of Saturn. The inner two planets’ masses are unconstrained by dynamical stability arguments. [ABSTRACT FROM AUTHOR]

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet–disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1–10 for a 1 star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets () produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet–disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory. [ABSTRACT FROM AUTHOR]

We analyze the concentration of solid particles in vortices created and sustained by radial buoyancy in protoplanetary disks, e.g., baroclinic vortex growth. Besides the gas drag acting on particles, we also allow for back-reaction from dust onto the gas. This becomes important when the local dust-to-gas ratio approaches unity. In our two-dimensional, local, shearing sheet simulations, we see high concentrations of grains inside the vortices for a broad range of Stokes numbers, St. An initial dust-to-gas ratio of 1:100 can easily be reversed to 100:1 for St = 1.0. The increased dust-to-gas ratio triggers the streaming instability, thus counter-intuitively limiting the maximal achievable overdensities. We find that particle trapping inside vortices opens the possibility for gravity assisted planetesimal formation even for small particles (St = 0.01) and a low initial dust-to-gas ratio of 1:104, e.g., much smaller than in the previously studied magnetohydrodynamic zonal flow case. [ABSTRACT FROM AUTHOR]

The vacuum ultraviolet (VUV) synthesis of CO and CO2 on ice-coated graphite and isotopic labeled 13C graphite has been examined for temperatures between 40 and 120 K. The results show that CO and CO2 can be formed at the buried ice:graphite interface with Lyα photon irradiation via the reaction of radicals (O and OH) produced by direct photodissociation and the dissociative electron attachment of the interfacial water molecules. The synthesized CO and CO2 molecules can desorb in hot photon-dominated regions and are lost to space when ice coated carbonaceous dust grains cycle within the protoplanetary disks. Thus, the nonthermal formation of CO and CO2 at the buried ice:grain interface by VUV photons may help regulate the carbon inventory during the early stage of planet formation. This may contribute to the carbon deficits in our solar system and suggests that a universal carbon deficit gradient may be expected within astrophysical bodies surrounding center stars. [ABSTRACT FROM AUTHOR]

We study the evolutionary behavior of the snow line with a time-dependent disk model that includes the mass influx onto the disk from the collapse of the progenitor cloud core. We are able to investigate the evolution of the snow line during the disk formation and early evolution. We find that the radius of the snow line in a disk first increases with time, reaches a maximum, , and then decreases with time. The mass influx in our model inherits the progenitor core properties and we quantitatively illustrate the relationship between the behavior of the snow line and the core properties. We find that the mass influx onto the disk becomes dilute for high core angular momentum and that surface density near the snow line decreases. Thus viscosity heating decreases and becomes small. When the angular momentum is low, the amount of material falling onto the disk is small and thus the surface density is small. This limits the value of . We find that the maximum value of is for the core mass . We suggest that diffusion of water vapor to makes a preferred location for gas giant planet formation. Water inside is substantially depleted due to diffusion. The observed lack of cold water vapor in protoplanetary disks by Herschel might be related to the depletion of water inside . [ABSTRACT FROM AUTHOR]

We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified 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 and , but in the AD dominated regime both Scr and Scz are . 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. [ABSTRACT FROM AUTHOR]

We report on recent Ks-band interferometric observations of the young pre-main-sequence star AB Aurigae obtained with the Palomar Fiber Nuller (PFN). Reaching a contrast of a few 10–4 inside a field of view extending from 35 to 275 mas (5-40 AU at AB Aur's distance), the PFN is able to explore angular scales that are intermediate between those accessed by coronagraphic imaging and long baseline interferometry. This intermediate region is of special interest given that many young stellar objects are believed to harbor extended halos at such angular scales. Using destructive interference (nulling) between two sub-apertures of the Palomar 200 inch telescope and rotating the telescope pupil, we measured a resolved circumstellar excess at all probed azimuth angles. The astrophysical null measured over the full rotation is fairly constant, with a mean value of 1.52%, and a slight additional azimuthal modulation of ±0.2%. The isotropic astrophysical null is indicative of circumstellar emission dominated by an azimuthally extended source, possibly a halo, or one or more rings of dust, accounting for several percent of the total Ks-band flux. The modest azimuthal variation may be explained by some skewness or anisotropy of the spatially extended source, e.g., an elliptical or spiral geometry, or clumping, but it could also be due to the presence of a point source located at a separation of ∼120 mas (17 AU) with ∼6 × 10–3 of the stellar flux. We combine our results with previous Infrared Optical Telescope Array observations of AB Aur at H band, and demonstrate that a dust ring located at ∼30 mas (4.3 AU) represents the best-fitting model to explain both sets of visibilities. We are also able to test a few previously hypothesized models of the incoherent component evident at longer interferometric baselines. [ABSTRACT FROM AUTHOR]

Dust grains in neutral gas behave as aerodynamic particles, so they can develop large local density fluctuations entirely independent of gas density fluctuations. Specifically, gas turbulence can drive order-of-magnitude “resonant” fluctuations in the dust density on scales where the gas stopping/drag timescale is comparable to the turbulent eddy turnover time. Here we show that for large grains (size ≳ 0.1 μm, containing most grain mass) in sufficiently large molecular clouds (radii ≳ 1-10 pc, masses ≳ 104M☼), this scale becomes larger than the characteristic sizes of prestellar cores (the sonic length), so large fluctuations in the dust-to-gas ratio are imprinted on cores. As a result, star clusters and protostellar disks formed in large clouds should exhibit significant abundance spreads in the elements preferentially found in large grains (C, O). This naturally predicts populations of carbon-enhanced stars, certain highly unusual stellar populations observed in nearby open clusters, and may explain the “UV upturn” in early-type galaxies. It will also dramatically change planet formation in the resulting protostellar disks, by preferentially “seeding” disks with an enhancement in large carbonaceous or silicate grains. The relevant threshold for this behavior scales simply with cloud densities and temperatures, making straightforward predictions for clusters in starbursts and high-redshift galaxies. Because of the selective sorting by size, this process is not necessarily visible in extinction mapping. We also predict the shape of the abundance distribution—when these fluctuations occur, a small fraction of the cores may actually be seeded with abundances Z ∼ 100 〈 Z〉 such that they are almost “totally metal” (Z ∼ 1)! Assuming the cores collapse, these totally metal stars would be rare (1 in ∼104 in clusters where this occurs), but represent a fundamentally new stellar evolution channel. [ABSTRACT FROM AUTHOR]

We present Atacama Large Millimeter/Submillimeter Array (ALMA) observations of a wide binary system in Orion, with projected separation 440 AU, in which we detect submillimeter emission from the protoplanetary disks around each star. Both disks appear moderately massive and have strong line emission in CO 3-2, HCO+ 4-3, and HCN 3-2. In addition, CS 7-6 is detected in one disk. The line-to-continuum ratios are similar for the two disks in each of the lines. From the resolved velocity gradients across each disk, we constrain the masses of the central stars, and show consistency with optical-infrared spectroscopy, both indicative of a high mass ratio ∼9. The small difference between the systemic velocities indicates that the binary orbital plane is close to face-on. The angle between the projected disk rotation axes is very high, ∼72°, showing that the system did not form from a single massive disk or a rigidly rotating cloud core. This finding, which adds to related evidence from disk geometries in other systems, protostellar outflows, stellar rotation, and similar recent ALMA results, demonstrates that turbulence or dynamical interactions act on small scales well below that of molecular cores during the early stages of star formation. [ABSTRACT FROM AUTHOR]

The collision dynamics of dusty bodies are crucial for planetesimal formation. Decimeter agglomerates are especially important in the different formation models. Therefore, in continuation of our experiments on mutual decimeter collisions, we investigate collisions of centimeter onto decimeter dust agglomerates in a small drop tower under vacuum conditions (p ≲ 5 × 10–1 mbar) at a mean collision velocity of 6.68 ± 0.67 m s–1. We use quartz dust with irregularly shaped micrometer grains. Centimeter projectiles with different diameters, masses, and heights are used, their typical volume filling factor is Φp, m = 0.466 ± 0.02. The decimeter agglomerates have a mass of about 1.5 kg, a diameter and height of 12 cm, and a mean filling factor of Φt, m = 0.44 ± 0.004. At lower collision energies, only the projectile gets destroyed and mass is transferred to the target. The accretion efficiency decreases with increasing obliquity and increasing difference in filling factor, if the projectile is more compact than the target. The accretion efficiency increases with increasing collision energy for collision energies under a certain threshold. Beyond this threshold at 298 ± 25 mJ, catastrophic disruption of the target can be observed. This corresponds to a critical fragmentation strength Q* = 190 ± 16 mJ kg–1, which is a factor of four larger than expected. Analyses of the projectile fragments show a power-law size distribution with an average exponent of –3.8 ± 0.3. The mass distributions suggest that the fraction of smallest fragments increases for higher collision energies. This is interesting for impacts of small particles on large target bodies within protoplanetary disks, as smaller fragments couple better to the surrounding gas and re-accretion by gas drag is more likely. [ABSTRACT FROM AUTHOR]

We describe an open source GPU implementation of a hybrid symplectic N-body integrator, GENGA (Gravitational ENcounters with Gpu Acceleration), designed to integrate planet and planetesimal dynamics in the late stage of planet formation and stability analyses of planetary systems. GENGA uses a hybrid symplectic integrator to handle close encounters with very good energy conservation, which is essential in long-term planetary system integration. We extended the second-order hybrid integration scheme to higher orders. The GENGA code supports three simulation modes: integration of up to 2048 massive bodies, integration with up to a million test particles, or parallel integration of a large number of individual planetary systems. We compare the results of GENGA to Mercury and pkdgrav2 in terms of energy conservation and performance and find that the energy conservation of GENGA is comparable to Mercury and around two orders of magnitude better than pkdgrav2. GENGA runs up to 30 times faster than Mercury and up to 8 times faster than pkdgrav2. GENGA is written in CUDA C and runs on all NVIDIA GPUs with a computing capability of at least 2.0. [ABSTRACT FROM AUTHOR]

Super-Earths with orbital periods less than 100 days are extremely abundant around Sun-like stars. It is unlikely that these planets formed at their current locations. Rather, they likely formed at large distances from the star and subsequently migrated inward. Here we use N-body simulations to study the effect of super-Earths on the accretion of rocky planets. In our simulations, one or more super-Earths migrate inward through a disk of planetary embryos and planetesimals embedded in a gaseous disk. We tested a wide range of migration speeds and configurations. Fast-migrating super-Earths (τmig ∼ 0.01-0.1 Myr) only have a modest effect on the protoplanetary embryos and planetesimals. Sufficient material survives to form rocky, Earth-like planets on orbits exterior to the super-Earths'. In contrast, slowly migrating super-Earths shepherd rocky material interior to their orbits and strongly deplete the terrestrial planet-forming zone. In this situation any Earth-sized planets in the habitable zone are extremely volatile-rich and are therefore probably not Earth-like. [ABSTRACT FROM AUTHOR]

The instability in protoplanetary disks due to gas-dust friction and self-gravity of gas and dust is investigated using linear analysis. In the case where the dust-to-gas ratio is enhanced and turbulence is weak, the instability grows, even in gravitationally stable disks, on a timescale of order 104-5 yr at a radius of order 100 AU. If we ignore the dynamical feedback from dust grains in the gas equation of motion, the instability reduces to the so-called “secular gravitational instability,” which was investigated previously to be an instability of dust in a fixed background gas flow. In this work, we solve the equations of motion for both gas and dust consistently and find that long-wavelength perturbations are stable, in contrast to the secular gravitational instability in the simplified treatment. This may indicate that we should not neglect small terms in the equation of motion if the growth rate is small. The instability is expected to form ring structures in protoplanetary disks. The width of the ring formed at a radius of 100 AU is a few tens of AU. Therefore, the instability is a candidate for the formation mechanism of observed ring-like structures in disks. Another aspect of the instability is the accumulation of dust grains, and hence the instability may play an important role in the formation of planetesimals, rocky protoplanets, and cores of gas giants located at radii ∼100 AU. If these objects survive the dispersal of the gaseous component of the disk, they may be the origin of debris disks. [ABSTRACT FROM AUTHOR]

We study the influence of the turbulent transport on ice chemistry in protoplanetary disks, focusing on carbon- and nitrogen-bearing molecules. Chemical rate equations are solved with the diffusion term, mimicking the turbulent mixing in the vertical direction. Turbulence can bring ice-coated dust grains from the midplane to the warm irradiated disk surface, and the ice mantles are reprocessed by photoreactions, thermal desorption, and surface reactions. The upward transport decreases the abundance of methanol and ammonia ices at r ≲ 30 AU because warm dust temperature prohibits their reformation on grain surfaces. This reprocessing could explain the smaller abundances of carbon and nitrogen bearing molecules in cometary coma than those in low-mass protostellar envelopes. We also show the effect of mixing on the synthesis of complex organic molecules (COMs) in two ways: (1) transport of ices from the midplane to the disk surface and (2) transport of atomic hydrogen from the surface to the midplane. The former enhances the COMs formation in the disk surface, while the latter suppresses it in the midplane. Then, when mixing is strong, COMs are predominantly formed in the disk surface, while their parent molecules are (re)formed in the midplane. This cycle expands the COMs distribution both vertically and radially outward compared with that in the non-turbulent model. We derive the timescale of the sink mechanism by which CO and N2 are converted to less volatile molecules to be depleted from the gas phase and find that the vertical mixing suppresses this mechanism in the inner disks. [ABSTRACT FROM AUTHOR]

Recent simulations have shown that the formation of planets in circumbinary configurations (such as those recently discovered by Kepler) is dramatically hindered at the planetesimal accretion stage. The combined action of the binary and the protoplanetary disk acts to raise impact velocities between kilometer-sized planetesimals beyond their destruction threshold, halting planet formation within at least 10 AU from the binary. It has been proposed that a primordial population of “large” planetesimals (100 km or more in size), as produced by turbulent concentration mechanisms, would be able to bypass this bottleneck; however, it is not clear whether these processes are viable in the highly perturbed circumbinary environments. We perform two-dimensional hydrodynamical and N-body simulations to show that kilometer-sized planetesimals and collisional debris can drift and be trapped in a belt close to the central binary. Within this belt, planetesimals could initially grow by accreting debris, ultimately becoming “indestructible” seeds that can accrete other planetesimals in situ despite the large impact speeds. We find that large, indestructible planetesimals can be formed close to the central binary within 105 yr, therefore showing that even a primordial population of “small” planetesimals can feasibly form a planet. [ABSTRACT FROM AUTHOR]

We present precision radial velocity (RV) data sets from Keck-HIRES and from Lick Observatory's new Automated Planet Finder Telescope and Levy Spectrometer on Mt. Hamilton that reveal a multiple-planet system orbiting the nearby, slightly evolved, K-type star HD 141399. Our 91 observations over 10.5 yr suggest the presence of four planets with orbital periods of 94.35, 202.08, 1070.35, and 3717.35 days and minimum masses of 0.46, 1.36, 1.22, and 0.69 MJ, respectively. The orbital eccentricities of the three inner planets are small, and the phase curves are well sampled. The inner two planets lie just outside the 2:1 resonance, suggesting that the system may have experienced dissipative evolution during the protoplanetary disk phase. The fourth companion is a Jupiter-like planet with a Jupiter-like orbital period. Its orbital eccentricity is consistent with zero, but more data will be required for an accurate eccentricity determination. [ABSTRACT FROM AUTHOR]

A correction to the article "Long-Term Evoluaton of Photoevaporating Protoplanetary Disks," by Jaehan Bae and colleagues published in the October 7, 2013 issue of the journal.