Observational signature of planet formation: The ALMA view

Protoplanetary disks are the most probable sites where planet formation takes place. According to theory, planet formation in protoplanetary disks should show remarkable signatures, such as a gap/hole or a spiral structure. In fact, recent high-angular and high-sensitivity observations in millimeter and submillimeter wavelengths, as well as optical/near-IR wavelengths, have shown such structures in protoplanetary disks. Two particular examples of such disks around AB Aurigae and HD 142527 are discussed here, with an emphasis on results obtained using the Submillimeter Array. These disks—and their probable planet formation—will be very important future targets for ALMA to study the physical process of planet formation in detail.     

A submillimeter view of protoplanetary dust disks

We present some results from our submillimeter single-dish and aperture synthesis imaging surveys of protoplanetary disks using the JCMT, CSO, and Submillimeter Array (SMA) on Mauna Kea, Hawaii. Employing a simple disk model, we simultaneously fit the spectral energy distributions and spatially resolved submillimeter continuum emission from our SMA survey to constrain disk structure properties, including surface density profiles and sizes. The typical disk structure we infer is consistent with a fiducial accretion disk model with a viscosity parameter α≈0.01. Combined with a large, multiwavelength single-dish survey of similar disks, we show how these observations provide evidence for significant grain growth and rapid evolution in the outer regions of disks, perhaps due to an internal photoevaporation process. In addition, we discuss SMA observations of the disks in the Orion Trapezium (proplyds) in the context of disk evolution in a more extreme environment.     

Star and planet-formation with ALMA: an overview

Submillimeter observations with ALMA will be the essential next step in our understanding of how stars and planets form. Key projects range from detailed imaging of the collapse of pre-stellar cores and measuring the accretion rate of matter onto deeply embedded protostars, to unravelling the chemistry and dynamics of high-mass star-forming clusters and high-spatial resolution studies of protoplanetary disks down to the 1 AU scale.     

High-mass star forming regions: An ALMA view

The advent of ALMA is bound to improve our knowledge of OB star formation dramatically. Here, we present an overview of this topic outlining how high angular resolution and sensitivity may contribute to shed light on the structure of high-mass star forming regions and hence on the process itself of massive star formation. The impact of this new generation instrument will range from establishing the mass function of pre-stellar cores inside IR-dark clouds, to investigating the kinematics of the gas from which OB stars are built up, to assessing or ruling out the existence of circumstellar accretion disks in these objects.     

Observations of molecular clouds in nearby galaxies with ALMA

We present recent results of the observations of giant molecular clouds in nearby galaxies with the Nobeyama 45 m telescope and Millimeter Array. We give some brief comments about observations of GMCs in nearby galaxies with ALMA.     

Observations of circumstellar disks with infrared interferometry

The study of circumstellar disks around young stellar objects is arguably the area of astrophysics on which the technique of infrared interferometry has had the biggest impact. Here I will review the existing set of observations in this field, concentrating on disks but also including jets/winds and stellar properties. At the end, there is a brief discussion of how ongoing technical developments and observational improvements will expand the impact of infrared interferometry on the study of star formation.     

The study of evolved stars with ALMA

Intense mass loss occurs for low- and intermediate-mass stars on the asymptotic giant branch (AGB), and for the higher mass (≳8 M_⊙) stars during their red supergiant evolution. These winds affect the evolution of the stars profoundly, creates circumstellar envelopes of gas and dust, as well as enrich the interstellar medium with heavy elements and grain particles. The mass loss characteristics are well-studied, but the basic processes are still not understood in detail, and the mass-loss rate of an individual star cannot be derived from first principles. These objects also provide us with fascinating systems, in which intricate interplays between various physical and chemical processes take place, and their relative simplicity in terms of geometry, density distribution, and kinematics makes them excellent astrophysical laboratories. The review concentrates on the aspects of AGB stars and their mass loss which are of particular interest in connection with ALMA.     

On oligarchic growth of planets in protoplanetary disks

In this paper we present a new semianalytical model of oligarchic growth of planets considering a distribution of planetesimal sizes, fragmentation of planetesimals in mutual collisions, sublimation of ices through the snow line, random velocities out of equilibrium and merging of planetary embryos. We show that the presence of several planetary embryos growing simultaneously at different locations in the protoplanetary disk affects the whole accretion history, specially for the innermost planets. The results presented here clearly indicate the relevance of considering a distribution of planetesimal sizes. Fragmentation occurring during planetesimal-planetesimal collisions represent only a marginal effect in shaping the surface density of solid material in the protoplanetary disc.     

Particle stirring in turbulent gas disks: Including orbital oscillations

We describe the diffusion and random velocities of solid particles due to stochastic forcing by turbulent gas. We include the orbital dynamics of Keplerian disks, both in-plane epicycles and vertical oscillations. We obtain a new result for the diffusion of solids. The Schmidt number (ratio of gas to particle diffusivity) is Sc≈1+(Ωt_stop)², in terms of the particle stopping time t_stop and the orbital frequency Ω. The standard result, Sc=1+t_stop/t_eddy, in terms of the eddy turnover time, t_eddy, is shown to be incorrect. The main difference is that Sc rises quadratically, not linearly, with stopping time. Consequently, particles larger than 10 cm in protoplanetary disks will suffer less radial diffusion and will settle closer to the midplane. Such a layer of boulders would be more prone to gravitational collapse. Our predictions of RMS speeds, vertical scale height and diffusion coefficients will help interpret numerical simulations. We confirm previous results for the vertical stirring of particles (scale heights and random velocities), and add a correction for arbitrary ratios of eddy to orbital times. The particle layer becomes thinner for t_eddy>1/Ω with the strength of turbulent diffusion held fixed. We use two analytic techniques—the Hinze–Tchen formalism and the Fokker–Planck equation with velocity diffusion—with identical results when the regimes of validity overlap. We include simple physical arguments for the scaling of our results.     

Two-dimensional transport of solids in viscous protoplanetary disks

Large-scale radial transport of solids appears to be a fundamental consequence of protoplanetary disk evolution based on the presence of high temperature minerals in comets and the outer regions of protoplanetary disks around other stars. Further, inward transport of solids from the outer regions of the solar nebula has been postulated to be the manner in which short-lived radionuclides were introduced to the terrestrial planet region and the cause of the variations in oxygen isotope ratios in the primitive materials. Here, both outward and inward transport of solids are investigated in the context of a two-dimensional, viscously evolving protoplanetary disk. The dynamics of solids are investigated to determine how they depend on particle size and the particular stage of protoplanetary disk evolution, corresponding to different rates of mass transport. It is found that the outward flows that arise around the disk midplane of a protoplanetary disk aid in the outward transport of solids up to the size of CAIs s and can increase the crystallinity fraction of silicate dust at 10 AU around a solar mass star to as much as ∼40% in the case of rapidly evolving disks, decreasing as the accretion rate onto the star slows. High velocity, inward flows along the disk surface aid in the rapid transport of solids from the outer disk to the inner disk, particularly for small dust. Despite the diffusion that occurs throughout the disk, the large-scale, meridonal flows associated with mass transport prevent complete homogenization of the disk, allowing compositional gradients to develop that vary in intensity for a timescale of one million of years. The variations in the rates and the preferred direction of radial transport with height above the disk midplane thus have important implications for the dynamics and chemical evolution of primitive materials.     

Studying the first galaxies with ALMA

We discuss observations of the first galaxies, within cosmic reionization, at centimeter and millimeter wavelengths. We present a summary of current observations of the host galaxies of the most distant QSOs (z∼6). These observations reveal the gas, dust, and star formation in the host galaxies on kpc-scales. These data imply an enriched ISM in the QSO host galaxies within 1 Gyr of the big bang, and are consistent with models of coeval supermassive black hole and spheroidal galaxy formation in major mergers at high redshift. Current instruments are limited to studying truly pathologic objects at these redshifts, meaning hyper-luminous infrared galaxies (L _FIR ∼10¹³ L _⊙). ALMA will provide the one to two orders of magnitude improvement in millimeter astronomy required to study normal star forming galaxies (i.e. Ly-α emitters) at z∼6. ALMA will reveal, at sub-kpc spatial resolution, the thermal gas and dust—the fundamental fuel for star formation—in galaxies into cosmic reionization.     

Studies of dense cores with ALMA

Dense cores are the simplest star-forming sites that we know, but despite their simplicity, they still hold a number of mysteries that limit our understanding of how solar-type stars form. ALMA promises to revolutionize our knowledge of every stage in the life of a core, from the pre-stellar phase to the final disruption by the newly born star. This contribution presents a brief review of the evolution of dense cores and illustrates particular questions that will greatly benefit from the increase in resolution and sensitivity expected from ALMA.     

The study of young substellar objects with ALMA

We discuss the potential of ALMA for studying the formation of substellar objects. We first review briefly the various formation mechanisms proposed so far and stress the unique capability of ALMA to detect and study pre-brown dwarf cores and to confirm the core-collapse scenario to the lowest possible masses. We then discuss the properties of disks around substellar objects. We show how it will be possible to detect with ALMA most disks around objects with mass as low as few Jupiter masses, and to resolve spatially their emission in the more favorable cases.     

Formation of Massive Stars with High Accretion

We make models of the formation of massive stars by an accretion which is growing with the mass already accreted. These models are supported by several observations. They have implications for the IMF and the maximum stellar mass. This revised version was published online in September 2006 with corrections to the Cover Date.     

Gemini observations of disks and jets in young stellar objects and in active galaxies

We present first results from the Near-infrared Integral Field Spectrograph (NIFS) located at Gemini North. For the active galaxies Cygnus A and Perseus A we observe rotationally-supported accretion disks and adduce the existence of massive central black holes and estimate their masses. In Cygnus A we also see remarkable high-excitation ionization cones dominated by photoionization from the central engine. In the T-Tauri stars HV Tau C and DG Tau we see highly-collimated bipolar outflows in the [Fe II] λ 1.644 micron line, surrounded by a slower molecular bipolar outflow seen in the H₂ lines, in accordance with the model advocated by Pyo, T.-S., et al., Astrophys. J. 570, 724 (2002).     

The interstellar chemistry of protostellar disks

A study has been undertaken of the gas-grain chemistry of protostellar disks which are sufficiently cool that in the outer regions, where the gas density is less than 10¹³ cm^–3 and the ionization rate highest, a bimolecular chemistry resembling that of dark clouds can occur. Since the gas-grain collision rate is so high, outgassing mantle molecules effectively determine the gas phase composition at any position in the disk. In contrast to previous work, a detailed gas phase chemistry is considered along with the accretion and desorption of mantle species which is controlled locally by the dust temperature.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: I. Effects of interaction with planetesimals and other protoplanets

We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (10⁶-10⁷ years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 10⁵ years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: II. Solid surface density evolution with type-I migration

This paper investigates the surface density evolution of a planetesimal disk due to the effect of type-I migration by carrying out N-body simulation and through analytical method, focusing on terrestrial planet formation. The coagulation and the growth of the planetesimals take place in the abundant gas disk except for a final stage. A protoplanet excites density waves in the gas disk, which causes the torque on the protoplanet. The torque imbalance makes the protoplanet suffer radial migration, which is known as type-I migration. Type-I migration time scale derived by the linear theory may be too short for the terrestrial planets to survive, which is one of the major problems in the planet formation scenario. Although the linear theory assumes a protoplanet being in a gas disk alone, Kominami et al. [Kominami, J., Tanaka, H., Ida, S., 2005. Icarus 167, 231-243] showed that the effect of the interaction with the planetesimal disk and the neighboring protoplanets on type-I migration is negligible. The migration becomes pronounced before the planet's mass reaches the isolation mass, and decreases the solid component in the disk. Runaway protoplanets form again in the planetesimal disk with decreased surface density. In this paper, we present the analytical formulas that describe the evolution of the solid surface density of the disk as a function of gas-to-dust ratio, gas depletion time scale and semimajor axis, which agree well with our results of N-body simulations. In general, significant depletion of solid material is likely to take place in inner regions of disks. This might be responsible for the fact that there is no planet inside Mercury's orbit in our Solar System. Our most important result is that the final surface density of solid components (Σd) and mass of surviving planets depend on gas surface density (Σg) and its depletion time scale (τdep) but not on initial Σd; they decrease with increase in Σg and τdep. For a fixed gas-to-dust ratio and τdep, larger initial Σd results in smaller final Σd and smaller surviving planets, because of larger Σg. To retain a specific amount of Σd, the efficient disk condition is not an initially large Σd but the initial Σd as small as the specified final one and a smaller gas-to-dust ratio. To retain Σd comparable to that of the minimum mass solar nebula (MMSN), a disk must have the same Σd and a gas-to-dust ratio that is smaller than that of MMSN by a factor of 1.3×(τdep/1 Myr) at ∼1 AU. (Equivalently, type-I migration speed is slower than that predicted by the linear theory by the same factor.) The surviving planets are Mars-sized ones in this case; in order to form Earth-sized planets, their eccentricities must be pumped up to start orbit crossing and coagulation among them. At ∼5 AU, Σd of MMSN is retained under the same condition, but to form a core massive enough to start runaway gas accretion, a gas-to-dust ratio must be smaller than that of MMSN by a factor of 3×τdep/1 Myr.