On the Dynamics of Dust Grains in a Hierarchical Environment

Most main sequence stars are binaries or higher multiplicity Systems and it appears that at birth most stars have circumstellar disks. It is commonly accepted that planetary systems arise from the material of these disks; consequently, binary and multiple systems may have a main role in planet formation. In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimetre-to-metre sized primordial aggregates. We investigate the response of the particles, in a protoplanetary disk with radius R_D = 100 AU around a solar-like star, to the gravitational field of bound perturbing companions in a moderately wide (300–1600 AU) orbit. For this purpose, we have carried out a series of simulations of coplanar hierarchical configurations using a direct integration code that models gravitational and viscous forces. The massive protoplanetary disk is around one of the components of the binary. The evolution in time of the dust sub-disk depends mainly on the nature (prograde or retrograde) of the relative revolution of the stellar companion, and on the temperature and mass of the circumstellar disk. Our results show that for binary companions near the limit of tidal truncation of the disk, the perturbation leads to an enhanced accretion rate onto the primary, decreasing the lifetime of the particles in the protoplanetary disk with respect to the case of a single star. As a consequence of an enhanced accretion rate the mass of the disk decreases faster, which leads to a longer resultant lifetime for particles in the disk. On the other hand, binary companions may induce tidal arms in the dust phase of protoplanetary disks. Spiral perturbations with m = 1 may increase in a factor 10 or more the dust surface density in the neighbourhood of the arm, facilitating the growth of the particles. Moreover, in a massive disk (0.01M⊙) the survival time of particles is significantly shorter than in a less massive nebula (0.001M⊙) and the temperature of the disk severely influences the spiral-in time of particles. The rapid evolution of the dust component found in post T Tauri stars can be explained as a result of their binary nature. Binarity may also influence the evolution of circumpulsar disks. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the nature of the proposed blob in the inner preplanetary disk of the Herbig Ae/Be star HD141569: evidence for a giant vortex?

Recently, Brittain and Rettig, using the cryogenic echelle spectrograph at the Infrared Telescope Facility to study the infrared emission from the inner preplanetary disk of the Herbig Ae/Be star HD141569, detected CO and H₃⁺ ion emission. This emission has been tentatively interpreted as due to the existence of a forming gas giant planet. The suggested protoplanetary blob appears to be orbiting its host star at about 7 AU being perhaps 2 AU across and roughly five times the mass of Jupiter. Based on numerical modeling of the evolution of the dust disk we show that their observational results are compatible with the presence of an evolved giant vortex in the disk. Our calculations suggest that vortices formed in disks similar to the one found around HD141569 are more effective at capturing solid material than equivalent structures around solar-like stars. On the other hand, we investigate the possibility to find evidence for large-scale vortices in preplanetary disks by submillimeter interferometry. Disks around Herbig Ae/Be stars may be primary targets for giant vortex detection using this technique.     

Two-dimensional models of proto planetary disk chemistry

We have developed a two-dimensional model of a flared protoplanetary disk (PPD) incorporating a self-consistent treatment of gas and dust temperature, and a detailed treatment of the gas-phase chemistry as well as the freeze-out and desorption of material from dust grains. The results show that, in the inner 10 AU of the disk, the gas-phase abundances are dominated by material evaporated from dust grains. The surface layer of the disk shows many of the characteristics of photon-dominated regions. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Study on the Characteristics of the Structure of Vega's Debris Disk

The clumpy structure in the Vega's debris disk was reported at millimeter wavelengths previously, and attributed to the concentration of dust grains trapped in resonances with a potential high-eccentricity planet. However, current imaging at multi-wavelengths with higher sensitivity indicates that the Vega's debris disk has a smooth structure. But a planet orbiting Vega could not be neglected, and the present-day observations may place a severe constraint on the orbital parameters for the potential planet. Herein, we utilize the modi- fied MERCURY codes to numerically simulate the Vega system, which consists of a debris disk and a planet. In our simulations, the initial inner and outer boundaries of the debris disk are assumed to be 80 AU and 120 AU, respectively. The dust grains in the disk have the sizes from 10 μm to 100 μm, and the nearly coplanar orbits. From the outcomes, we show that the evolution of debris disk is consistent with recent observations, if there is no planet orbiting Vega. However, if Vega owns a planet with a high eccentricity (e.g., e = 0.6), the planet's semi- major axis cannot be larger than 60 AU, otherwise, an aggregation phenomenon will occur in the debris disk due to the existence of the postulated planet. In addition, the 2:1 mean motion resonances may play a significant role in forming the structure of debris disk.     

Magnetic fields in protoplanetary disks

Magnetic fields likely play a key role in the dynamics and evolution of protoplanetary disks. They have the potential to efficiently transport angular momentum by MHD turbulence or via the magnetocentrifugal acceleration of outflows from the disk surface. Magnetically-driven mixing has implications for disk chemistry and evolution of the grain population, and the effective viscous response of the disk determines whether planets migrate inwards or outwards. However, the weak ionisation of protoplanetary disks means that magnetic fields may not be able to effectively couple to the matter. I examine the magnetic diffusivity in a minimum solar nebula model and present calculations of the ionisation equilibrium and magnetic diffusivity as a function of height from the disk midplane at radii of 1 and 5 AU. Dust grains tend to suppress magnetic coupling by soaking up electrons and ions from the gas phase and reducing the conductivity of the gas by many orders of magnitude. However, once grains have grown to a few microns in size their effect starts to wane and magnetic fields can begin to couple to the gas even at the disk midplane. Because ions are generally decoupled from the magnetic field by neutral collisions while electrons are not, the Hall effect tends to dominate the diffusion of the magnetic field when it is able to partially couple to the gas, except at the disk surfaces where the low density of neutrals permits the ions to remain attached to the field lines. For a standard population of 0.1 μm grains the active surface layers have a combined column Σ_active≈2 g cm⁻² at 1 AU; by the time grains have aggregated to 3 μm, Σ_active≈80 g cm⁻². Ionisation in the active layers is dominated by stellar X-rays. In the absence of grains, X-rays maintain magnetic coupling to 10% of the disk material at 1 AU (i.e. Σ_active≈150 g cm⁻²). At 5 AU the Σ_active≈Σ_total once grains have aggregated to 1 μm in size.     

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.     

Clumps in the outer disk by disk instability: Why they are initially gas giants and the legacy of disruption

We explore the initial conditions for fragments in the extended regions of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H₂ dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution.     

Mid-infrared spectra of the shocked Murchison CM chondrite: Comparison with astronomical observations of dust in debris disks

We present laboratory mid-infrared transmission/absorption spectra obtained from matrix of the hydrated Murchison CM meteorite experimentally shocked at peak pressures of 10-49 GPa, and compare them to astronomical observations of circumstellar dust in different stages of the formation of planetary systems. The laboratory spectra of the Murchison samples exhibit characteristic changes in the infrared features. A weakly shocked sample (shocked at 10 GPa) shows almost no changes from the unshocked sample dominated by hydrous silicate (serpentine). Moderately shocked samples (21-34 GPa) have typical serpentine features gradually replaced by bands of amorphous material and olivine with increasing shock pressure. A strongly shocked sample (36 GPa) shows major changes due to decomposition of the serpentine and due to devolatilization. A shock melted sample (49 GPa) shows features of olivine recrystallized from melted material.The infrared spectra of the shocked Murchison samples show similarities to astronomical spectra of dust in various young stellar objects and debris disks. The spectra of highly shocked Murchison samples (36 and 49 GPa) are similar to those of dust in the debris disks of HD113766 and HD69830, and the transitional disk of HD100546. The moderately shocked samples (21-34 GPa) exhibit spectra similar to those of dust in the debris disks of Beta Pictoris and BD+20307, and the transitional disk of GM Aur. An average of the spectra of all Murchison samples (0-49 GPa) has a similarity to the spectrum of the older protoplanetary disk of SU Auriga. In the gas-rich transitional and protoplanetary disks, the abundances of amorphous silicates and gases have widely been considered to be a primary property. However, our study suggests that impact processing may play a significant role in generating secondary amorphous silicates and gases in those disks. Infrared spectra of the shocked Murchison samples also show similarities to the dust from comets (C/2002 V1, C/2001 RX14, 9P/Tempel 1, and Hale Bopp), suggesting that the comets also contain shocked Murchison-like material.     

The Exo-Zodiacal Disk Mapper: A Space Interferometer to Detect and Map Zodiacal Dust Disks Around Nearby Stars

We present a concept for a space mission designed to make a mid-IR survey of potential zodiacal dust disks around nearby stars. We show that a two-aperture (0.6 m diameter), 10-m baseline, nulling interferometer located in a 1 × 4 AU, 4-yr solar orbit would allow for the survey of 400 stars in the solar neighborhood and permit a first-order determination of the disk inclination and radial dependences of density and temperature. The high dynamic range of the instrument may also be used to study additional astrophysical phenomena. Beyond its own scientific merit, such a mission would serve as a technological precursor to a larger interferometer of the type being considered for the detection of earth-like planets.     

Gravitational Instability in the Dust Layer of a Protoplanetary Disk with Interaction between the Layer and the Surrounding Gas

We consider gravitational instability of the dust layer in the midplane of a protoplanetary disk with turbulence and shear stresses between the gas in the disk and that in the dust layer. We solve a linearized system of hydrodynamic equations for perturbations of dust (monodisperse) and gas phases in the incompressible gas approximation. We take into account the gas drag of solid particles (dust aggregates), turbulent diffusion and the velocity dispersion of particles, and the perturbation of the azimuthal velocity of gas in the layer upon the transfer of angular momentum from solid particles to it and from this gas to the surrounding gas in the disk. We obtain and solve the dispersion equation for the layer with the ratio of surface densities of the dust phase and gas being well above unity. The following parameters of gravitational instability in the dust layer are calculated: the critical surface density of solid matter and the Stokes number of particles corresponding to the onset of instability, the wavelength range in which instability occurs, and the rate of its growth as a function of the perturbation wavelength in the circumsolar disk at radial distances of 1 and 10 AU. We show that at 10 AU, the maximum instability growth rate increases due to the transfer of angular momentum of gas in the layer to gas outside it, a new maximum emerges at a longer wavelength, a long-wavelength instability “tail” forms, and the critical surface density initiating instability decreases relative to that determined without the transfer of angular momentum to gas outside the layer. None of these effects are observed at 1 AU, since instability in this region probably develops faster than the transfer of angular momentum to the surrounding gаs of a protoplanetary disk occurs.     

A Submillimeter HCN Laser in IRC +10216

We present subarcsecond thermal infrared imaging of HD 98800, a young quadruple system composed of a pair of low-mass spectroscopic binaries separated by 0&farcs;8 (38 AU), each with a K-dwarf primary. Images at wavelengths ranging from 5 to 24.5 μm show unequivocally that the optically fainter binary, HD 98800B, is the sole source of a comparatively large infrared excess on which a silicate emission feature is superposed. The excess is detected only at wavelengths of 7.9 μm and longer, peaks at 25 μm, and has a best-fit blackbody temperature of 150 K, indicating that most of the dust lies at distances greater than the orbital separation of the spectroscopic binary. We estimate the radial extent of the dust with a disk model that approximates radiation from the spectroscopic binary as a single source of equivalent luminosity. Given the data, the most likely values of disk properties in the ranges considered are Rin=5.0+/-2.5 AU, DeltaR=13+/-8 AU, lambda0=2+4-1.5 μm, gamma=0+/-2.5, and sigmatotal=16+/-3 AU2, where Rin is the inner radius, DeltaR is the radial extent of the disk, lambda0 is the effective grain size, gamma is the radial power-law exponent of the optical depth tau, and sigmatotal is the total cross section of the grains. The range of implied disk masses is 0.001-0.1 times that of the Moon. These results show that, for a wide range of possible disk properties, a circumbinary disk is far more likely than a narrow ring.     

Dust ejection from planetary bodies by temperature gradients: Laboratory experiments

Laboratory experiments show that dusty bodies in a gaseous environment eject dust particles if they are illuminated. We find that even more intense dust eruptions occur when the light source is turned off. We attribute this to a compression of gas by thermal creep in response to the changing temperature gradients in the top dust layers. The effect is studied at a light flux of 13 kW/m² and 1 mbar ambient pressure. The effect is applicable to protoplanetary disks and Mars. In the inner part of protoplanetary disks, planetesimals can be eroded especially at the terminator of a rotating body. This leads to the production of dust which can then be transported towards the disk edge or the outer disk regions. The generated dust might constitute a significant fraction of the warm dust observed in extrasolar protoplanetary disks. We estimate erosion rates of about 1 kg s^?1 for 100 m parent bodies. The dust might also contribute to subsequent planetary growth in different locations or on existing protoplanets which are large enough not to be susceptible to particle loss by light induced ejection. Due to the ejections, planetesimals and smaller bodies will be accelerated or decelerated and drift outward or inward, respectively. The effect might also explain the entrainment of dust in dust devils on Mars, especially at high altitudes where gas drag alone might not be sufficient.     

Gravitational instability in the dust layer of a protoplanetary disk: Interaction of solid particles with turbulent gas in the layer

Gravitational instability of the dust layer formed after the aggregates of dust particles settle toward the midplane of a protoplanetary disk under turbulence is considered. A linearized system of hydrodynamic equations for perturbations of dust (monodisperse) and gas phases in the incompressible gas approximation is solved. Turbulent diffusion and the velocity dispersion of solid particles and the perturbation of gas azimuthal velocity in the layer upon the transfer of angular momentum from the dust phase due to gas drag are taken into account. Such an interaction of the particles and the gas establishes upper and lower bounds on the perturbation wavelength that renders the instability possible. The dispersion equation for the layer in the case when the ratio of surface densities of the dust phase and the gas in the layer is well above unity is obtained and solved. An approximate gravitational instability criterion, which takes the size-dependent stopping time of a particle (aggregate) in the gas into account, is derived. The following parameters of the layer instability are calculated: the wavelength range of its subsistence and the dependence of the perturbation growth rate on the perturbation wavelength in the circumsolar disk at a radial distance of 1 and 10 AU. It is demonstrated that at a distance of 1 AU, the gas–dust disk should be enriched with solids by a factor of 5–10 relative to the initial abundance as well as the particle aggregates should grow to the sizes higher than about 0.3 m in order for the instability to emerge in the layer in the available turbulence models. Such high disk enrichment and aggregate growth is not needed at a distance of 10 AU. The conditions under which this gravitational instability in the layer may be examined with no allowance made for the transfer of angular momentum from the gas in the layer to the gas in a protoplanetary disk outside the layer are discussed.     

A new view of proto-planetary disks with ALMA

The dynamical, physical and chemical processes which lead to planet formation constitute an astrophysical domain which will strongly benefit from ALMA in terms of frequency coverage, sensitivity and angular resolution. Recent results from current mm/submm interferometers obtained on molecules and dust in proto-planetary disks are presented. The observational coupling between gas and dust is discussed and it is shown that dust disks must be analyzed with the knowledge provided by gas disks, and respectively, both from the chemical and physical points. For these purposes, the methods of analysis of mm/submm interferometric data specific to disks are summarized. Emphasis is given on recent, unexpected, findings obtained in the highest sensitivity and resolution observations obtained so far, as they provide a hint of what ALMA could discover. A comparison with the expected sensitivities for ALMA illustrates how ALMA can enhance our knowledge of the disk physics, either by providing statistics or by allowing much more detailed studies of representative objects.     

Temperature distribution in the solar nebula at successive stages of its evolution

We have constructed a model of the solar nebula that allows for the temperature and pressure distributions at various stages of its evolution to be calculated. The mass flux from the accretion envelope to the disk and from the disk to the Sun, the turbulent viscosity parameter α, the opacity of the disk material, and the initial angular momentum of the protosun are the input model parameters that are varied. We also take into account the changes in the luminosity and radius of the young Sun. The input model parameters are based mostly on data obtained from observations of young solar-type stars with disks. To correct the input parameters, we use the mass and chemical composition of Jupiter, as well as models of its internal structure and formation that allow constraints to be imposed on the temperature and surface density of the protoplanetary disk in Jupiter’s formation zone. Given the derived constraints on the input parameters, we have calculated models of the solar nebula at successive stages of its evolution: the formation inside the accretion envelope, the evolution around the young Sun going through the T Tauri stage, and the formation and compaction of a thin dust layer (subdisk) in the disk midplane. We have found the following evolutionary trend: an increase in the temperature of the disk at the stage of its formation, cooling at the T Tauri stage, and the subsequent internal heating of the dust subdisk by turbulence dissipation that causes a temperature rise in the formation zone of the terrestrial planets at the high subdisk density and the opacity in this zone. We have obtained the probable ranges of temperatures in the disk midplane, i.e., the temperatures of the protoplanetary material in the formation region of the terrestrial planets at the initial stage of their formation.     

Astronomical evidence for the rapid growth of millimeter‐sized particles in protoplanetary disks

Abstract– I summarize recent surveys of protoplanetary disks at millimeter wavelengths and show that the distribution of luminosity, equivalent to the mass in small dust grains, declines rapidly. This contrasts with statistics on the lifetime of disks from infrared observations and the high occurrence of planets from radial velocity and transit surveys. I suggest that these disparate results can be reconciled if most of the dust in a disk is locked up in millimeter and larger‐sized particles within about 2 Myr. This statistical result on disk evolution agrees with detailed modeling of a small number of individual disks and with cosmochemical measurements of rapid planetesimal formation.     

Dust in the disk winds from young stars as a source of the circumstellar extinction

We consider the problem of dust grain survival in the disk winds from T Tauri and Herbig Ae stars. For our analysis, we have chosen a disk wind model in which the gas component of the wind is heated through ambipolar diffusion to a temperature of ～10⁴ K. We show that the heating of dust grains through their collisions with gas atoms is inefficient compared to their heating by stellar radiation and, hence, the grains survive even in the hot wind component. As a result, the disk wind can be opaque to the ultraviolet and optical stellar radiation and is capable of absorbing an appreciable fraction of it. Calculations show that the fraction of the wind-absorbed radiation for T Tauri stars can be from 20 to 40% of the total stellar luminosity at an accretion rate ? _a = 10^?8-10^?6 M _⊙yr^?1. This means that the disk winds from T Tauri stars can play the same role as the puffed-up inner rim in current accretion disk models. In Herbig Ae stars, the inner layers of the disk wind (r ≤ 0.5 AU) are dust-free, since the dust in this region sublimates under the effect of stellar radiation. Therefore, the fraction of the radiation absorbed by the disk wind in this case is considerably smaller and can be comparable to the effect from the puffed-up inner rim only at an accretion rate of the order of or higher than 10^?6 M _⊙yr^?1. Since the disk wind is structurally inhomogeneous, its optical depth toward the observer can be variable, which should be reflected in the photometric activity of young stars. For the same reason, moving shadows from gas and dust streams with a spiral-like shape can be observed in high-angular-resolution circumstellar disk images.     

eta Carinae: Binarity Confirmed

Gas giant planets have been detected in orbit around an increasing number of nearby stars. Two theories have been advanced for the formation of such planets: core accretion and disk instability. Core accretion, the generally accepted mechanism, requires several million years or more to form a gas giant planet in a protoplanetary disk like the solar nebula. Disk instability, on the other hand, can form a gas giant protoplanet in a few hundred years. However, disk instability has previously been thought to be important only in relatively massive disks. New three-dimensional, "locally isothermal," hydrodynamical models without velocity damping show that a disk instability can form Jupiter-mass clumps, even in a disk with a mass (0.091 M middle dot in circle within 20 AU) low enough to be in the range inferred for the solar nebula. The clumps form with initially eccentric orbits, and their survival will depend on their ability to contract to higher densities before they can be tidally disrupted at successive periastrons. Because the disk mass in these models is comparable to that apparently required for the core accretion mechanism to operate, the models imply that disk instability could obviate the core accretion mechanism in the solar nebula and elsewhere.     

High angular resolution millimeter observations of circumstellar disks

In this lecture, we review the properties of protoplanetary disks as derived from high angular resolution observations at millimeter wavelengths. We discuss how the combination of several different high angular resolution techniques allow us to probe different regions of the disk around young stellar objects and to derive the properties of the dust when combined with sophisticated disk models. The picture that emerges is that the dust in circumstellar disks surrounding pre-main sequence stars is in many cases significantly evolved compared to the dust in molecular clouds and the interstellar medium. It is however still difficult to derive a consistent picture and timeline for dust evolution in disks as the observations are still limited to small samples of objects.We also review the evidence for and properties of disks around high-mass young stellar objects and the implications on their formation mechanisms. The study of massive YSOs is complicated by their short lifetimes and larger average distances. In most cases high angular resolution data at millimeter wavelengths are the only method to probe the structure of disks in these objects.We provide a summary of the characteristics of available high angular resolution millimeter and submillimeter observatories. We also describe the characteristics of the ALMA observatory being constructed in the Chilean Andes. ALMA is going to be the world leading observatory at millimeter wavelengths in the coming decades, the project is now in its main construction phase with early science activities envisaged for 2010 and full science operations for 2012.