Source author record

C. P. Dullemond

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

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR astro-ph astro-ph.GA astro-ph.IM

Catalog footprint

What is connected

55works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2019arXiv

Spiral arms in the proto-planetary disc HD100453 detected with ALMA: evidence for binary-disc interaction and a vertical temperature gradient

Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.

preprint2016arXiv

Investigating dust trapping in transition disks with millimeter-wave polarization

Spatially resolved polarized (sub-)mm emission has been observed for example in the protoplanetary disk around HL Tau. Magnetically aligned grains are commonly interpreted as the source of polarization. However, self-scattering by large dust grains with a high enough albedo is another polarization mechanism, becoming a compelling method independent of the spectral index to constrain the dust grain size in protoplanetary disks. We study the dust polarization at mm wavelength in the dust trapping scenario proposed for transition disks, when a giant planet opens a gap in the disk. We investigate the characteristic polarization patterns and their dependence on disk inclination, dust size evolution, planet position, and observing wavelength. We combine two-dimensional hydrodynamical simulations of planet-disk interactions with self-consistent dust growth models. These size-dependent dust density distributions are used for follow-up three-dimensional radiative transfer calculations to predict the polarization degree at ALMA bands due to scattered thermal emission. We find that the polarization pattern of a disk with a planetary gap after 1 Myr of dust evolution shows a distinctive three ring structure. Two narrow inner rings are located at the planet gap edges. For increasing observing wavelengths all three rings slightly change their position, where the innermost and outermost rings move inward. This distance is detectable comparing the results at ALMA bands 3, 6 and 7. Within the highest polarized intensity regions the polarization vectors are oriented in the azimuthal direction. For an inclined disk there is an interplay between polarization originating from a flux gradient and inclination-induced quadrupole polarization. For intermediate inclined transition disks the polarization degree is as high as ~ 2% at band 3, which is well above the detection limit of future ALMA observations.

preprint2016arXiv

Shadows and spirals in the protoplanetary disk HD 100453

Understanding the diversity of planets requires to study the morphology and the physical conditions in the protoplanetary disks in which they form. We observed and spatially resolved the disk around the ~10 Myr old protoplanetary disk HD 100453 in polarized scattered light with SPHERE/VLT at optical and near-infrared wavelengths, reaching an angular resolution of ~0.02", and an inner working angle of ~0.09". We detect polarized scattered light up to ~0.42" (~48 au) and detect a cavity, a rim with azimuthal brightness variations at an inclination of 38 degrees, two shadows and two symmetric spiral arms. The spiral arms originate near the location of the shadows, close to the semi major axis. We detect a faint spiral-like feature in the SW that can be interpreted as the scattering surface of the bottom side of the disk, if the disk is tidally truncated by the M-dwarf companion currently seen at a projected distance of ~119 au. We construct a radiative transfer model that accounts for the main characteristics of the features with an inner and outer disk misaligned by ~72 degrees. The azimuthal brightness variations along the rim are well reproduced with the scattering phase function of the model. While spirals can be triggered by the tidal interaction with the companion, the close proximity of the spirals to the shadows suggests that the shadows could also play a role. The change in stellar illumination along the rim, induces an azimuthal variation of the scale height that can contribute to the brightness variations. Dark regions in polarized images of transition disks are now detected in a handful of disks and often interpreted as shadows due to a misaligned inner disk. The origin of such a misalignment in HD 100453, and of the spirals, is unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit.

preprint2015arXiv

A multi-wavelength analysis for interferometric (sub-)mm observations of protoplanetary disks: radial constraints on the dust properties and the disk structure

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

preprint2015arXiv

Scattered light images of spiral arms in marginally gravitationally unstable discs with an embedded planet

Scattered light images of transition discs in the near-infrared often show non-axisymmetric structures in the form of wide-open spiral arms in addition to their characteristic low-opacity inner gap region. We study self-gravitating discs and investigate the influence of gravitational instability on the shape and contrast of spiral arms induced by planet-disc interactions. Two-dimensional non-isothermal hydrodynamical simulations including viscous heating and a cooling prescription are combined with three-dimensional dust continuum radiative transfer models for direct comparison to observations. We find that the resulting contrast between the spirals and the surrounding disc in scattered light is by far higher for pressure scale height variations, i.e. thermal perturbations, than for pure surface density variations. Self-gravity effects suppress any vortex modes and tend to reduce the opening angle of planet-induced spirals, making them more tightly wound. If the disc is only marginally gravitationally stable with a Toomre parameter around unity, an embedded massive planet (planet-to-star mass ratio of $10^{-2}$) can trigger gravitational instability in the outer disc. The spirals created by this instability and the density waves launched by the planet can overlap resulting in large-scale, more open spiral arms in the outer disc. The contrast of these spirals is well above the detection limit of current telescopes.

preprint2014arXiv

Interferometer predictions with triangulated images: solving the multi-scale problem

Interferometers play an increasingly important role for spatially resolved observations. If employed at full potential, interferometry can probe an enormous dynamic range in spatial scale. Interpretation of the observed visibilities requires the numerical compu- tation of Fourier integrals over the synthetic model images. To get the correct values of these integrals, the model images must have the right size and resolution. Insufficient care in these choices can lead to wrong results. We present a new general-purpose scheme for the computation of visibilities of radiative transfer images. Our method requires a model image that is a list of intensities at arbitrarily placed positions on the image-plane. It creates a triangulated grid from these vertices, and assumes that the intensity inside each triangle of the grid is a linear function. The Fourier integral over each triangle is then evaluated with an analytic expression and the complex visibility of the entire image is then the sum of all triangles. The result is a robust Fourier trans- form that does not suffer from aliasing effects due to grid regularities. The method automatically ensures that all structure contained in the model gets reflected in the Fourier transform.

preprint2014arXiv

Mean gas opacity for circumstellar environments and equilibrium temperature degeneracy

In a molecular cloud dust opacity typically dominates over gas opacity, yet in the vicinities of forming stars dust is depleted, and gas is the sole provider of opacity. In the optically thin circumstellar environments the radiation temperature cannot be assumed to be equal to the gas temperature, hence the two-temperature Planck means are necessary to calculate the radiative equilibrium. By using the two-temperature mean opacity one does obtain the proper equilibrium gas temperature in a circumstellar environment, which is in a chemical equilibrium. A careful consideration of a radiative transfer problem reveals that the equilibrium temperature solution can be degenerate in an optically thin gaseous environment. We compute mean gas opacities based on the publicly available code DFSYNTHE by Kurucz and Castelli. We performed the calculations assuming local thermodynamic equilibrium and an ideal gas equation of state. The values were derived by direct integration of the high-resolution opacity spectrum. We produced two sets of gas opacity tables: Rosseland means and two-temperature Planck means (the tables available via http://cdsweb.u-strasbg.fr/ as well as via http://www.mpia-hd.mpg.de/homes/malygin). For three metallicities [Me/H] = 0.0,+/-0.3 we covered the parameter range 3.48 <= log T_rad[K] <= 4.48 in radiation temperature, 2.8 <= log T_gas[K]} <= 6.0 in gas temperature, and -10 <= log P[dyn cm^-2] <= 6 in gas pressure. We show that in the optically thin circumstellar environment for a given stellar radiation field and local gas density there are several equilibrium gas temperatures possible. We conclude that, in general, equilibrium gas temperature cannot be determined without treating the temperature evolution.

preprint2014arXiv

Millimetre spectral indices of transition disks and their relation to the cavity radius

Transition disks are protoplanetary disks with inner depleted dust cavities and excellent candidates to investigate the dust evolution under the existence of a pressure bump. A pressure bump at the outer edge of the cavity allows dust grains from the outer regions to stop their rapid inward migration towards the star and efficiently grow to millimetre sizes. Dynamical interactions with planet(s) have been one of the most exciting theories to explain the clearing of the inner disk. We look for evidence of the presence of millimetre dust particles in transition disks by measuring their spectral index with new and available photometric data. We investigate the influence of the size of the dust depleted cavity on the disk integrated millimetre spectral index. We present the 3mm photometric observations carried out with PdBI of four transition disks: LkHa330, UXTauA, LRLL31, and LRLL67. We use available values of their fluxes at 345GHz to calculate their spectral index, as well as the spectral index for a sample of twenty transition disks. We compare the observations with two kind of models. In the first set of models, we consider coagulation and fragmentation of dust in a disk in which a cavity is formed by a massive planet located at different positions. The second set of models assumes disks with truncated inner parts at different radius and with power-law dust size distributions, where the maximum size of grains is calculated considering turbulence as the source of destructive collisions. We show that the integrated spectral index is higher for transition disks than for regular protoplanetary disks. For transition disks, the probability that the measured spectral index is positively correlated with the cavity radius is 95%. High angular resolution imaging of transition disks is needed to distinguish between the dust trapping scenario and the truncated disk case.

preprint2014arXiv

Planet-vortex interaction:How a vortex can shepherd a planetary embryo

Context: Anticyclonic vortices are considered as a favourable places for trapping dust and forming planetary embryos. On the other hand, they are massive blobs that can interact gravitationally with the planets in the disc. Aims: We aim to study how a vortex interacts gravitationally with a planet which migrates toward it or a planet which is created inside the vortex. Methods: We performed hydrodynamical simulations of a viscous locally isothermal disc using GFARGO and FARGO-ADSG. We set a stationary Gaussian pressure bump in the disc in a way that RWI is triggered. After a large vortex is established, we implanted a low mass planet in the outer disc or inside the vortex and allowed it to migrate. We also examined the effect of vortex strength on the planet migration and checked the validity of the final result in the presence of self-gravity. Results: We noticed regardless of the planet's initial position, the planet is finally locked to the vortex or its migration is stopped in a farther orbital distance in case of a stronger vortex. For the model with the weaker vortex, we studied the effect of different parameters such as background viscosity, background surface density, mass of the planet and different planet positions. In these models, while the trapping time and locking angle of the planet vary for different parameters, the main result, which is the planet-vortex locking, remains valid. We discovered that even a planet with a mass less than 5 * 10^{-7} M_{\star} comes out from the vortex and is locked to it at the same orbital distance. For a stronger vortex, both in non-self-gravitated and self-gravitating models, the planet migration is stopped far away from the radial position of the vortex. This effect can make the vortices a suitable place for continual planet formation under the condition that they save their shape during the planetary growth.

preprint2014arXiv

Radiation hydrodynamics including irradiation and adaptive mesh refinement with AZEuS. I. Methods

Aims. The importance of radiation to the physical structure of protoplanetary disks cannot be understated. However, protoplanetary disks evolve with time, and so to understand disk evolution and by association, disk structure, one should solve the combined and time-dependent equations of radiation hydrodynamics. Methods. We implement a new implicit radiation solver in the AZEuS adaptive mesh refinement magnetohydrodynamics fluid code. Based on a hybrid approach that combines frequency-dependent ray-tracing for stellar irradiation with non-equilibrium flux limited diffusion, we solve the equations of radiation hydrodynamics while preserving the directionality of the stellar irradiation. The implementation permits simulations in Cartesian, cylindrical, and spherical coordinates, on both uniform and adaptive grids. Results. We present several hydrostatic and hydrodynamic radiation tests which validate our implementation on uniform and adaptive grids as appropriate, including benchmarks specifically designed for protoplanetary disks. Our results demonstrate that the combination of a hybrid radiation algorithm with AZEuS is an effective tool for radiation hydrodynamics studies, and produces results which are competitive with other astrophysical radiation hydrodynamics codes.

preprint2013arXiv

A quantification of hydrodynamical effects on protoplanetary dust growth

Context. The growth process of dust particles in protoplanetary disks can be modeled via numerical dust coagulation codes. In this approach, physical effects that dominate the dust growth process often must be implemented in a parameterized form. Due to a lack of these parameterizations, existing studies of dust coagulation have ignored the effects a hydrodynamical gas flow can have on grain growth, even though it is often argued that the flow could significantly contribute either positively or negatively to the growth process. Aims. We intend to provide a quantification of hydrodynamical effects on the growth of dust particles, such that these effects can be parameterized and implemented in a dust coagulation code. Methods. We numerically integrate the trajectories of small dust particles in the flow of disk gas around a proto-planetesimal, sampling a large parameter space in proto-planetesimal radii, headwind velocities, and dust stopping times. Results. The gas flow deflects most particles away from the proto-planetesimal, such that its effective collisional cross section, and therefore the mass accretion rate, is reduced. The gas flow however also reduces the impact velocity of small dust particles onto a proto-planetesimal. This can be beneficial for its growth, since large impact velocities are known to lead to erosion. We also demonstrate why such a gas flow does not return collisional debris to the surface of a proto-planetesimal. Conclusions. We predict that a laminar hydrodynamical flow around a proto-planetesimal will have a significant effect on its growth. However, we cannot easily predict which result, the reduction of the impact velocity or the sweep-up cross section, will be more important. Therefore, we provide parameterizations ready for implementation into a dust coagulation code.

preprint2013arXiv

Asymmetric transition disks: Vorticity or eccentricity?

Context. Transition disks typically appear in resolved millimeter observations as giant dust rings surrounding their young host stars. More accurate observations with ALMA have shown several of these rings to be in fact asymmetric: they have lopsided shapes. It has been speculated that these rings act as dust traps, which would make them important laboratories for studying planet formation. It has been shown that an elongated giant vortex produced in a disk with a strong viscosity jump strikingly resembles the observed asymmetric rings. Aims. We aim to study a similar behavior for a disk in which a giant planet is embedded. However, a giant planet can induce two kinds of asymmetries: (1) a giant vortex, and (2) an eccentric disk. We studied under which conditions each of these can appear, and how one can observationally distinguish between them. This is important because only a vortex can trap particles both radially and azimuthally, while the eccentric ring can only trap particles in radial direction. Methods. We used the FARGO code to conduct the hydro-simulations. We set up a disk with an embedded giant planet and took a radial grid spanning from 0.1 to 7 times the planet semi-major axis. We ran the simulations with various viscosity values and planet masses for 1000 planet orbits to allow a fully developed vortex or disk eccentricity. Afterwards, we compared the dust distribution in a vortex-holding disk with an eccentric disk using dust simulations. Results. We find that vorticity and eccentricity are distinguishable by looking at the azimuthal contrast of the dust density. While vortices, as particle traps, produce very pronounced azimuthal asymmetries, eccentric features are not able to accumulate millimeter dust particles in azimuthal direction, and therefore the asymmetries are expected to be modest.

preprint2013arXiv

Explaining millimeter-sized particles in brown dwarf disks

Planets have been detected around a variety of stars, including low-mass objects, such as brown dwarfs. However, such extreme cases are challenging for planet formation models. Recent sub-millimeter observations of disks around brown dwarf measured low spectral indices of the continuum emission that suggest that dust grains grow to mm-sizes even in these very low mass environments. To understand the first steps of planet formation in scaled-down versions of T-Tauri disks, we investigate the physical conditions that can theoretically explain the growth from interstellar dust to millimeter-sized grains in disks around brown dwarf. We modeled the evolution of dust particles under conditions of low-mass disks around brown dwarfs. We used coagulation, fragmentation and disk-structure models to simulate the evolution of dust, with zero and non-zero radial drift. For the non-zero radial drift, we considered strong inhomogeneities in the gas surface density profile that mimic long-lived pressure bumps in the disk. We studied different scenarios that could lead to an agreement between theoretical models and the spectral slope found by millimeter observations. We find that fragmentation is less likely and rapid inward drift is more significant for particles in brown dwarf disks than in T-Tauri disks. We present different scenarios that can nevertheless explain millimeter-sized grains. As an example, a model that combines the following parameters can fit the millimeter fluxes measured for brown dwarf disks: strong pressure inhomogeneities of $\sim$ 40% of amplitude, a small radial extent $\sim$ 15 AU, a moderate turbulence strength $α_{\mathrm{turb}}= 10^{-3}$, and average fragmentation velocities for ices $v_f = 10 m s^{-1}$.

preprint2013arXiv

Herbig Stars' Near-Infrared Excess: An Origin in the Protostellar Disk's Magnetically-Supported Atmosphere

Young stars with masses 2-8 Suns, called the Herbig Ae and Be stars, often show a near-infrared excess too large to explain with a hydrostatically-supported circumstellar disk of gas and dust. At the same time the accretion flow carrying the circumstellar gas to the star is thought to be driven by magneto-rotational turbulence, which according to numerical MHD modeling yields an extended low-density atmosphere supported by the magnetic fields. We demonstrate that the base of the atmosphere can be optically-thick to the starlight and that the parts lying near 1 AU are tall enough to double the fraction of the stellar luminosity reprocessed into the near-infrared. We generate synthetic spectral energy distributions (SEDs) using Monte Carlo radiative transfer calculations with opacities for sub-micron silicate and carbonaceous grains. The synthetic SEDs closely follow the median Herbig SED constructed recently by Mulders and Dominik, and in particular match the large near-infrared flux, provided the grains have a mass fraction close to interstellar near the disk's inner rim.

preprint2013arXiv

Lopsided dust rings in transition disks

Context. Particle trapping in local or global pressure maxima in protoplanetary disks is one of the new paradigms in the theory of the first stages of planet formation. However, finding observational evidence for this effect is not easy. Recent work suggests that the large ring-shaped outer disks observed in transition disk sources may in fact be lopsided and constitute large banana-shaped vortices. Aims. We wish to investigate how effective dust can accumulate along the azimuthal direction. We also want to find out if the size- sorting resulting from this can produce a detectable signatures at millimeter wavelengths. Methods. To keep the numerical cost under control we develop a 1+1D method in which the azimuthal variations are treated sepa- rately from the radial ones. The azimuthal structure is calculated analytically for a steady-state between mixing and azimuthal drift. We derive equilibration time scales and compare the analytical solutions to time-dependent numerical simulations. Results. We find that weak, but long-lived azimuthal density gradients in the gas can induce very strong azimuthal accumulations of dust. The strength of the accumulations depends on the Péclet number, which is the relative importance of advection and diffusion. We apply our model to transition disks and our simulated observations show that this effect would be easily observable with ALMA and in principle allows to put constraints on the strength of turbulence and the local gas density.

preprint2013arXiv

Planetesimal formation via sweep-up growth at the inner edge of dead zones

The early stages of planet formation are still not well understood. Coagulation models have revealed numerous obstacles to the dust growth, such as the bouncing, fragmentation and radial drift barriers. We study the interplay between dust coagulation and drift in order to determine the conditions in protoplanetary disk that support the formation of planetesimals. We focus on planetesimal formation via sweep-up and investigate whether it can take place in a realistic protoplanetary disk. We have developed a new numerical model that resolves spatial distribution of dust in the radial and vertical dimension. The model uses representative particles approach to follow the dust evolution in protoplanetary disk. The coagulation and fragmentation of solids is taken into account using Monte Carlo method. A collision model adopting the mass transfer effect, that can occur for different-sized dust aggregate collisions, is implemented. We focus on a protoplanetary disk including a pressure bump caused by a steep decline of turbulent viscosity around the snow line. Our results show that sufficient resolution of the vertical disk structure in dust coagulation codes is necessary to obtain adequately short growth timescales, especially in the case of a low turbulence region. We find that a sharp radial variation of the turbulence strength at the inner edge of dead zone promotes planetesimal formation in several ways. It provides a pressure bump that efficiently prevents the dust from drifting inwards. It also causes a radial variation in the size of aggregates at which growth barriers occur, favoring the growth of large aggregates via sweeping up of small particles. In our model, by employing an ad hoc alpha viscosity change near the snow line, it is possible to grow planetesimals by incremental growth on timescales of approximately 10^5 years.

preprint2012arXiv

Kinematics of the CO Gas in the Inner Regions of the TW Hya Disk

We present a detailed analysis of the spatially and spectrally resolved 12CO J=2-1 and J=3-2 emission lines from the TW Hya circumstellar disk, based on science verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These lines exhibit substantial emission in their high-velocity wings (with projected velocities out to 2.1 km/s, corresponding to intrinsic orbital velocities >20 km/s) that trace molecular gas as close as 2 AU from the central star. However, we are not able to reproduce the intensity of these wings and the general spatio-kinematic pattern of the lines with simple models for the disk structure and kinematics. Using three-dimensional non-local thermodynamic equilibrium molecular excitation and radiative transfer calculations, we construct some alternative models that successfully account for these features by modifying either (1) the temperature structure of the inner disk (inside the dust-depleted disk cavity; r < 4 AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) the distribution of disk inclination angles (a warp). The latter approach is particularly compelling because a representative warped disk model qualitatively reproduces the observed azimuthal modulation of optical light scattered off the disk surface. In any model scenario, the ALMA data clearly require a substantial molecular gas reservoir located inside the region where dust optical depths are known to be substantially diminished in the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. The results from these updated model prescriptions are discussed in terms of their potential physical origins, which might include dynamical perturbations from a low-mass companion with an orbital separation of a few AU.

preprint2012arXiv

Mid-Infrared Spectral Variability Atlas of Young Stellar Objects

Optical and near-infrared variability is a well-known property of young stellar objects. However, a growing number of recent studies claim that a considerable fraction of them also exhibit mid-infrared flux changes. With the aim of studying and interpreting variability on a decadal timescale, here we present a mid-infrared spectral atlas containing observations of 68 low- and intermediate mass young stellar objects. The atlas consists of 2.5-11.6 um low-resolution spectra obtained with the ISOPHOT-S instrument on-board the Infrared Space Observatory (ISO) between 1996 and 1998, as well as 5.2-14.5 um low-resolution spectra obtained with the IRS instrument on-board the Spitzer Space Telescope between 2004 and 2007. The observations were retrieved from the ISO and Spitzer archives and were post-processed interactively by our own routines. For those 47 objects where multi-epoch spectra were available, we analyze mid-infrared spectral variability on annual and/or decadal timescales. We identify 37 variable candidate sources. Many stars show wavelength-independent flux changes, possibly due to variable accretion rate. In several systems, all exhibiting 10 um silicate emission, the variability of the 6-8 um continuum and the silicate feature exhibit different amplitudes. A possible explanation is variable shadowing of the silicate emitting region by an inner disk structure of changing height or extra silicate emission from dust clouds in the disk atmosphere. Our results suggest that mid-infrared variability, in particular the wavelength-dependent changes, are more ubiquitous than was known before. Interpreting this variability is a new possibility to explore the structure of the disk and its dynamical processes.

preprint2011arXiv

Accretion through the inner hole of transitional disks: What happens to the dust?

We study the effect of radiation pressure on the dust in the inner rim of transitional disks with large inner holes. In particular, we evaluate whether radiation pressure can be responsible for keeping the inner holes dust-free, while allowing gas accretion to proceed. This has been proposed in a paper by Chiang and Murray-Clay (2007, Nature Physics 3, p. 604) who explain the formation of these holes as an inside-out evacuation due to X- ray-triggered accretion of the innermost layer of the disk rim outside of the hole. We show that radiation pressure is clearly incapable of stopping dust from flowing into the hole because of dust pile-up and optical depth effects, and also because of viscous mixing. Other mechanisms need to be found to explain the persistence of the opacity hole in the presence of accretion, and we speculate on possible solutions.

preprint2011arXiv

Near-infrared spectroscopy of EX Lupi in outburst

EX Lup is the prototype of the EXor class of young eruptive stars: objects showing repetitive brightenings due to increased accretion from the circumstellar disk to the star. In this paper, we report on medium-resolution near-infrared spectroscopy of EX\,Lup taken during its extreme outburst in 2008, as well as numerical modeling with the aim of determining the physical conditions around the star. We detect emission lines from atomic hydrogen, helium, and metals, as well as first overtone bandhead emission from carbon monoxide. Our results indicate that the emission lines are originating from gas located in a dust-free region within ~ 0.2 AU of the star. The profile of the CO bandhead indicates that the CO gas has a temperature of 2500 K, and is located in the inner edge of the disk or in the outer parts of funnel flows. The atomic metals are probably co-located with the CO. Some metallic lines are fluorescently excited, suggesting direct exposure to ultraviolet photons. The Brackett series indicates emission from hot (10000 K) and optically thin gas. The hydrogen lines display a strong spectro-astrometric signal, suggesting that the hydrogen emission is probably not coming from an equatorial boundary layer; a funnel flow or disk wind origin is more likely. This picture is broadly consistent with the standard magnetospheric accretion model usually assumed for normally accreting T Tauri stars. Our results also set constraints on the eruption mechanism, supporting a model where material piles up around the corotation radius and episodically falls onto the star.

preprint2011arXiv

Possible planet-forming regions on submillimetre images

Submillimetre images of transition discs are expected to reflect the distribution of the optically thin dust. Former observation of three transition discs LkHa330, SR21N, and HD1353444B at submillimetre wavelengths revealed images which cannot be modelled by a simple axisymmetric disc. We show that a large-scale anticyclonic vortex that develops where the viscosity has a large gradient (e.g., at the edge of the disc dead zone), might be accountable for these large-scale asymmetries. We modelled the long-term evolution of vortices being triggered by the Rossby wave instability. We found that a horseshoe-shaped (azimuthal wavenumber m=1) large-scale vortex forms by coalescing of smaller vortices within 5x10^4 yr, and can survive on the disc life-time (~5x10^6 yr), depending on the magnitude of global viscosity and the thickness of the viscosity gradient. The two-dimensional grid-based global disc simulations with local isothermal approximation and compressible-gas model have been done by the GPU version of hydrodynamic code FARGO (GFARGO). To calculate the dust continuum image at submillimetre wavelengths, we combined our hydrodynamical results with a 3D radiative transfer code. By the striking similarities of the calculated and observed submillimetre images, we suggest that the three transition discs can be modelled by a disc possessing a large-scale vortex formed near the disc dead zone edge. Since the larger dust grains (larger than mm in size) are collected in these vortices, the non-axisymmetric submillimetre images of the above transition discs might be interpreted as active planet and planetesimal forming regions situated far (> 50 AU) from the central stars.

preprint2011arXiv

Resolved Images of Large Cavities in Protoplanetary Transition Disks

Circumstellar disks are thought to experience a rapid "transition" phase in their evolution that can have a considerable impact on the formation and early development of planetary systems. We present new and archival high angular resolution (0.3" = 40-75 AU) Submillimeter Array (SMA) observations of the 880 micron dust continuum emission from 12 such transition disks in nearby star-forming regions. In each case, we directly resolve a dust-depleted disk cavity around the central star. Using radiative transfer calculations, we interpret these dust disk structures in a homogeneous, parametric model framework by reproducing their SMA visibilities and SEDs. The cavities in these disks are large (R_cav = 15-73 AU) and substantially depleted of small (~um-sized) dust grains, although their mass contents are still uncertain. The structures of the remnant material at larger radii are comparable to normal disks. We demonstrate that these large cavities are common among the millimeter-bright disk population, comprising at least 20% of the disks in the bright half of the millimeter luminosity (disk mass) distribution. Utilizing these results, we assess some of the physical mechanisms proposed to account for transition disk structures. As has been shown before, photoevaporation models do not produce the large cavity sizes, accretion rates, and disk masses representative of this sample. It would be difficult to achieve a sufficient decrease of the dust optical depths in these cavities by particle growth alone: substantial growth (to meter sizes or beyond) must occur in large (tens of AU) regions of low turbulence without also producing an abundance of small particles. Given those challenges, we suggest instead that the observations are most commensurate with dynamical clearing due to tidal interactions with low-mass companions --young brown dwarfs or giant planets on long-period orbits.

preprint2011arXiv

Spectral signatures of disk eccentricity in young binary systems: I. Circumprimary case

Star formation occurs via fragmentation of molecular clouds, which means that the majority of stars born are a members of binaries. There is growing evidence that planets might form in circumprimary disks of medium-separation binaries. The tidal forces caused by the secondary generally act to distort the originally circular disk to an eccentric one. To infer the disk eccentricity from high-res NIR spectroscopy, we calculate the fundamental band emission lines of the CO molecule emerging from the atmosphere of the disk. We model circumprimary disk evolution under the gravitational perturbation of the orbiting secondary using a 2D grid-based hydrodynamical code, assuming alpha-type viscosity. The hydrodynamical results are combined with our spectral code based on the double-layer disk model to calculate the CO molecular line profiles. We find that the orbital velocity distribution of the gas parcels differs significantly from the circular Keplerian fashion, thus the line profiles are asymmetric in shape. The magnitude of asymmetry is insensitive to the binary mass ratio, the magnitude of viscosity, and the disk mass. In contrast, the disk eccentricity, thus the level of the line profile asymmetry, is influenced significantly by the binary eccentricity and the disk geometrical thickness. We demonstrate that the disk eccentricity profile in the planet-forming region can be determined by fitting the high-resolution CO line profile asymmetry using a simple 2D spectral model that accounts for the velocity distortions caused by the disk eccentricity. Thus, with our novel approach the disk eccentricity can be inferred with high-resolution near-IR spectroscopy prior to the era of high angular resolution optical or radio direct-imaging. By determining the disk eccentricity in medium-separation young binaries, we might be able to constrain the planet formation theories.

preprint2011arXiv

The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? III. Sedimentation driven coagulation inside the snow-line

We investigate dust growth due to settling in a 1D vertical column of a protoplanetary disk. It is known from the observed 10 micron feature in disk SEDs, that small micron-sized grains are present at the disk atmosphere throughout the lifetime of the disk. We hope to explain such questions as what process can keep the disk atmospheres dusty for the lifetime of the disk and how does the particle properties change as a function of height above the midplane. We use a Monte Carlo code to follow the mass and porosity evolution of the particles in time. The used collision model is based on laboratory experiments performed on dust aggregates. As the experiments cannot cover all possible collision scenarios, the largest uncertainty of our model is the necessary extrapolations we had to perform. We simultaneously solve for the particle growth and motion. Particles can move vertically due to settling and turbulent mixing. We assume that the vertical profile of the gas density is fixed in time and only the solid component evolves. We find that the used collision model strongly influences the masses and sizes of the particles. The laboratory experiment based collision model greatly reduces the particle sizes compared to models that assume sticking at all collision velocities. We find that a turbulence parameter of alpha = 10^-2 is needed to keep the dust atmospheres dusty, but such strong turbulence can produce only small particles at the midplane which is not favorable for planetesimal formation models. We also see that the particles are larger at the midplane and smaller at the upper layers of the disk. At 3-4 pressure scale heights micron-sized particles are produced. These particle sizes are needed to explain the 10 micron feature of disk SEDs. Turbulence may therefore help to keep small dust particles in the disk atmosphere.

preprint2011arXiv

Trapping dust particles in the outer regions of protoplanetary disks

In order to explain grain growth to mm sized particles and their retention in outer regions of protoplanetary disks, as it is observed at sub-mm and mm wavelengths, we investigate if strong inhomogeneities in the gas density profiles can slow down excessive radial drift and can help dust particles to grow. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million years time scales. In addition, we introduced a comparison between the pressure inhomogeneities considered in this work and the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study if there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. We present the favorable conditions to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of the Taurus Auriga or Ophiucus star forming regions.

preprint2011arXiv

Warm gas at 50 AU in the disk around Herbig Be star HD 100546

The disk atmosphere is one of the fundamental elements of theoretical models of a protoplanetary disk. However, the direct observation of the warm gas (>> 100 K) at large radius of a disk (>> 10 AU) is challenging, because the line emission from warm gas in a disk is usually dominated by the emission from an inner disk. Our goal is to detect the warm gas in the disk atmosphere well beyond 10 AU from a central star in a nearby disk system of the Herbig Be star HD 100546. We measured the excitation temperature of the vibrational transition of CO at incremental radii of the disk from the central star up to 50 AU, using an adaptive optics system combined with the high-resolution infrared spectrograph CRIRES at the VLT. The observation successfully resolved the line emission with 0".1 angular resolution, which is 10 AU at the distance of HD 100546. Population diagrams were constructed at each location of the disk, and compared with the models calculated taking into account the optical depth effect in LTE condition. The excitation temperature of CO is 400-500 K or higher at 50 AU away from the star, where the blackbody temperature in equilibrium with the stellar radiation drops as low as 90 K. This is unambiguous evidence of a warm disk atmosphere far away from the central star.

preprint2010arXiv

A new condition for the transition from runaway to oligarchic growth

Accretion among macroscopic bodies of ~km size or larger is enhanced significantly due to gravitational focusing. Two regimes can be distinguished. Initially, the system experiences runaway growth, in which the gravitational focusing factors increase, and bodies at the high-mass tail of the distribution grow fastest. However, at some point the runaway body dynamically heats its environment, gravitational focusing factors decrease, and runaway growth passes into oligarchic growth. Based on the results of recent simulations, we reconsider the runaway growth-oligarchy transition. In contrast to oligarchy, we find that runaway growth cannot be approximated with a two component model (of small and large bodies) and that the criterion of Ida & Makino (1993), which is frequently adopted as the start of oligarchy, is not a sufficient condition to signify the transition. Instead, we propose a new criterion based on timescale arguments. We then find a larger value for the runaway growth-oligarchy transition: from several hundreds of km in the inner disk regions up to ~10^3 km for the outer disk. These findings are consistent with the view that runaway growth has been responsible for the size distribution of the present day Kuiper belt objects. Our finding furthermore outlines the proper initial conditions at the start of the oligarchy stage.

preprint2010arXiv

Accretion among preplanetary bodies: the many faces of runaway growth

(abridged) When preplanetary bodies reach proportions of ~1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R_0. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion.

preprint2010arXiv

Dust size distributions in coagulation/fragmentation equilibrium: Numerical solutions and analytical fits

Context. Grains in circumstellar disks are believed to grow by mutual collisions and subsequent sticking due to surface forces. Results of many fields of research involving circumstellar disks, such as radiative transfer calculations, disk chemistry, magneto-hydrodynamic simulations largely depend on the unknown grain size distribution. Aims. As detailed calculations of grain growth and fragmentation are both numerically challenging and computationally expensive, we aim to find simple recipes and analytical solutions for the grain size distribution in circumstellar disks for a scenario in which grain growth is limited by fragmentation and radial drift can be neglected. Methods. We generalize previous analytical work on self-similar steady-state grain distributions. Numerical simulations are carried out to identify under which conditions the grain size distributions can be understood in terms of a combination of power-law distributions. A physically motivated fitting formula for grain size distributions is derived using our analytical predictions and numerical simulations. Results. We find good agreement between analytical results and numerical solutions of the Smoluchowski equation for simple shapes of the kernel function. The results for more complicated and realistic cases can be fitted with a physically motivated "black box" recipe presented in this paper. Our results show that the shape of the dust distribution is mostly dominated by the gas surface density (not the dust-to-gas ratio), the turbulence strength and the temperature and does not obey an MRN type distribution.

preprint2010arXiv

Dust, Ice and Gas in Time (DIGIT) Herschel program first results: A full PACS-SED scan of the gas line emission in protostar DK Cha

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

preprint2010arXiv

Fundamental Vibrational Transition of CO during the Outburst of EX Lupi in 2008

We report monitoring observations of the T Tauri star EX Lupi during its outburst in 2008 in the CO fundamental band at 4.6-5.0 um. The observations were carried out at the VLT and the Subaru Telescope at six epochs from April to August 2008, covering the plateau of the outburst and the fading phase to a quiescent state. The line flux of CO emission declines with the visual brightness of the star and the continuum flux at 5 um, but composed of two subcomponents that decay with different rates. The narrow line emission (50 km s-1 in FWHM) is near the systemic velocity of EX Lupi. These emission lines appear exclusively in v=1-0. The line widths translate to a characteristic orbiting radius of 0.4 AU. The broad line component (FWZI ~ 150 km s-1) is highly excited upto v<=6. The line flux of the component decreases faster than the narrow line emission. Simple modeling of the line profiles implies that the broad-line emitting gas is orbiting around the star at 0.04-0.4 AU. The excitation state, the decay speed of the line flux, and the line profile, indicate that the broad-line emission component is physically distinct from the narrow-line emission component, and more tightly related to the outburst event.

preprint2010arXiv

Gas- and dust evolution in protoplanetary disks

Context. Current models of the size- and radial evolution of dust in protoplanetary disks generally oversimplify either the radial evolution of the disk (by focussing at one single radius or by using steady state disk models) or they assume particle growth to proceed monodispersely or without fragmentation. Further studies of protoplanetary disks - such as observations, disk chemistry and structure calculations or planet population synthesis models - depend on the distribution of dust as a function of grain size and radial position in the disk. Aims. We attempt to improve upon current models to be able to investigate how the initial conditions, the build-up phase, and the evolution of the protoplanetary disk influence growth and transport of dust. Methods. We introduce a new version of the model of Brauer et al. (2008) in which we now include the time-dependent viscous evolution of the gas disk, and in which more advanced input physics and numerical integration methods are implemented. Results. We show that grain properties, the gas pressure gradient, and the amount of turbulence are much more influencing the evolution of dust than the initial conditions or the build-up phase of the protoplanetary disk. We quantify which conditions or environments are favorable for growth beyond the meter size barrier. High gas surface densities or zonal flows may help to overcome the problem of radial drift, however already a small amount of turbulence poses a much stronger obstacle for grain growth.

preprint2010arXiv

High-resolution spectroscopic view of planet formation sites

Theories of planet formation predict the birth of giant planets in the inner, dense, and gas-rich regions of the circumstellar disks around young stars. These are the regions from which strong CO emission is expected. Observations have so far been unable to confirm the presence of planets caught in formation. We have developed a novel method to detect a giant planet still embedded in a circumstellar disk by the distortions of the CO molecular line profiles emerging from the protoplanetary disk's surface. The method is based on the fact that a giant planet significantly perturbs the gas velocity flow in addition to distorting the disk surface density. We have calculated the emerging molecular line profiles by combining hydrodynamical models with semianalytic radiative transfer calculations. Our results have shown that a giant Jupiter-like planet can be detected using contemporary or future high-resolution near-IR spectrographs such as VLT/CRIRES or ELT/METIS. We have also studied the effects of binarity on disk perturbations. The most interesting results have been found for eccentric circumprimary disks in mid-separation binaries, for which the disk eccentricity - detectable from the asymmetric line profiles - arises from the gravitational effects of the companion star. Our detailed simulations shed new light on how to constrain the disk kinematical state as well as its eccentricity profile. Recent findings by independent groups have shown that core-accretion is severely affected by disk eccentricity, hence detection of an eccentric protoplanetary disk in a young binary system would further constrain planet formation theories.

preprint2010arXiv

Impact of grain evolution on the chemical structure of protoplanetary disks

We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensive set of gas phase and grain surface chemical reactions is used to simulate the chemical structure of the disk. The main effect of the grain evolution on the disk chemical composition comes from sedimentation, and, to a lesser degree, from the reduction of the total grain surface area. The net effect of grain growth is suppressed by the fragmentation process which maintains a population of small grains, dominating the total grain surface area. We consider three models of dust properties. In model GS both growth and sedimentation are taken into account. In models A5 and A4 all grains are assumed to have the same size (10(-5) cm and 10(-4) cm, respectively) with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three-layer" pattern (midplane, molecular layer, hot atmosphere) in the disk chemical structure is preserved in all models, but shifted closer to the midplane in models with increased grain size (GS and A4). Unlike other similar studies, we find that in models GS and A4 column densities of most gas-phase species are enhanced by 1-3 orders of magnitude relative to those in a model with pristine dust (A5), while column densities of their surface counterparts are decreased. We show that column densities of certain species, like C2H, HC(2n+1)N (n=0-3), H2O and some other molecules, as well as the C2H2/HCN abundance ratio which are accessible with Herschel and ALMA can be used as observational tracers of early stages of the grain evolution process in protoplanetary disks.

preprint2010arXiv

Protoplanetary Disk Structures in Ophiuchus II: Extension to Fainter Sources

We present new results from a significant extension of our previous high angular resolution (0.3" = 40 AU) Submillimeter Array survey of the 880 um continuum emission from dusty circumstellar disks in the ~1 Myr-old Ophiuchus star-forming region. An expanded sample is constructed to probe disk structures that emit significantly lower millimeter luminosities (hence dust masses), down to the median value for T Tauri stars. Using a Monte Carlo radiative transfer code, the millimeter visibilities and broadband spectral energy distribution for each disk are simultaneously reproduced with a two-dimensional parametric model for a viscous accretion disk. We find wide ranges of characteristic radii (14-198 AU) and disk masses (0.004-0.143 M_sun), but a narrow distribution of surface density gradients (0.4-1.1) that is consistent with a uniform value $γ$ = 0.9 +/- 0.2 and independent of mass (or millimeter luminosity). In this sample, we find a correlation between the disk luminosity/mass and characteristic radius, such that fainter disks are both smaller and less massive. We suggest that this relationship is an imprint of the initial conditions inherited by the disks at their formation epoch, compare their angular momenta with those of molecular cloud cores, and speculate on how future observations can help constrain the distribution of viscous evolution timescales. No other correlations between disk and star properties are found. The inferred disk structures are briefly compared with theoretical models for giant planet formation, although resolution limitations do not permit us to directly comment on material inside R = 20 AU. However, there is some compelling evidence for dust evolution in the planet formation region: 4/17 disks in the sample show resolved regions of significantly reduced optical depths within ~20-40 AU of their central stars.

preprint2010arXiv

Sub-Keplerian accretion onto circumstellar disks

Models of the formation, evolution and photoevaporation of circumstellar disks are an essential ingredient in many theories of the formation of planetary systems. The ratio of disk mass over stellar mass in the circumstellar phase of a disk is largely determined by the angular momentum of the original cloud core from which the system was formed. While full 3D or 2D axisymmetric hydrodynamical models of accretion onto the disk automatically treat all aspects of angular momentum, this is not so trivial for 1D and semi-2D viscous disk models. Since 1D and semi-2D disk models are still very useful for long-term evolutionary modelling of disks with relatively little numerical effort, we investigate how the 2D nature of accretion affects the formation and evolution of the disk in such models. A proper treatment of this problem requires a correction for the sub-Keplerian velocity at which accretion takes place. We develop an update of our semi-2D time-dependent disk evolution model to properly treat the effects of sub-Keplerian accretion. The new model also accounts for the effects of the vertical extent of the disk on the accretion streamlines from the envelope. The disks produced with the new method are smaller than those obtained previously, but their mass is mostly unchanged. The new disks are a few degrees warmer in the outer parts, so they contain less solid CO. Otherwise, the results for ices are unaffected. The 2D treatment of the accretion results in material accreting at larger radii, so a smaller fraction comes close enough to the star for amorphous silicates to be thermally annealed into crystalline form. The lower crystalline abundances thus predicted correspond more closely to observed abundances than did earlier model predictions. We argue that thermal annealing followed by radial mixing must be responsible for at least part of the observed crystalline material.

preprint2010arXiv

Testing the theory of grain growth and fragmentation by millimeter observations of protoplanetary disks

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.

preprint2010arXiv

The inner regions of protoplanetary disks

To understand how planetary systems form in the dusty disks around pre-main-sequence stars a detailed knowledge of the structure and evolution of these disks is required. While this is reasonably well understood for the regions of the disk beyond about 1 AU, the structure of these disks inward of 1 AU remains a puzzle. This is partly because it is very difficult to spatially resolve these regions with current telescopes. But it is also because the physics of this region, where the disk becomes so hot that the dust starts to evaporate, is poorly understood. With infrared interferometry it has become possible in recent years to directly spatially resolve the inner AU of protoplanetary disks, albeit in a somewhat limited way. These observations have partly confirmed current models of these regions, but also posed new questions and puzzles. Moreover, it has turned out that the numerical modeling of these regions is extremely challenging. In this review we give a rough overview of the history and recent developments in this exciting field of astrophysics.

preprint2010arXiv

The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier

The sticking of micron sized dust particles due to surface forces in circumstellar disks is the first stage in the production of asteroids and planets. The key ingredients that drive this process are the relative velocity between the dust particles in this environment and the complex physics of dust aggregate collisions. Here we present the results of a collision model, which is based on laboratory experiments of these aggregates. We investigate the maximum aggregate size and mass that can be reached by coagulation in protoplanetary disks. We model the growth of dust aggregates at 1 AU at the midplane at three different gas densities. We find that the evolution of the dust does not follow the previously assumed growth-fragmentation cycles. Catastrophic fragmentation hardly occurs in the three disk models. Furthermore we see long lived, quasi-steady states in the distribution function of the aggregates due to bouncing. We explore how the mass and the porosity change upon varying the turbulence parameter and by varying the critical mass ratio of dust particles. Particles reach Stokes numbers of roughly 10^-4 during the simulations. The particle growth is stopped by bouncing rather than fragmentation in these models. The final Stokes number of the aggregates is rather insensitive to the variations of the gas density and the strength of turbulence. The maximum mass of the particles is limited to approximately 1 gram (chondrule-sized particles). Planetesimal formation can proceed via the turbulent concentration of these aerodynamically size-sorted chondrule-sized particles.

preprint2010arXiv

The thermal structure and the location of the snow line in the protosolar nebula: axisymmetric models with full 3-D radiative transfer

The precise location of the water ice condensation front ('snow line') in the protosolar nebula has been a debate for a long time. Its importance stems from the expected substantial jump in the abundance of solids beyond the snow line, which is conducive to planet formation, and from the higher stickiness in collisions of ice-coated dust grains, which may help the process of coagulation of dust and the formation of planetesimals. In an optically thin nebula, the location of the snow line is easily calculated to be around 3 AU. However, in its first 5 to 10 million years, the solar nebula was optically thick, implying a smaller snow line radius due to shielding from direct sunlight, but also a larger radius because of viscous heating. Several models have attempted to treat these opposing effects. However, until recently treatments beyond an approximate 1+1D radiative transfer were unfeasible. We revisit the problem with a fully self-consistent 3D treatment in an axisymmetric disk model, including a density-dependent treatment of the dust and ice sublimation. We find that the location of the snow line is very sensitive to the opacities of the dust grains and the mass accretion rate of the disk. We show that previous approximate treatments are quite efficient at determining the location of the snow line if the energy budget is locally dominated by viscous accretion. We derive an analytic estimate of the location of the snow line that compares very well with results from this and previous studies. Using solar abundances of the elements we compute the abundance of dust and ice and find that the expected jump in solid surface density at the snow line is smaller than previously assumed. We further show that in the inner few AU the refractory species are partly evaporated, leading to a significantly smaller solid state surface density in the regions where the rocky planets were formed.

preprint2009arXiv

Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation

Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration. We test the accuracy of seven independently developed radiative transfer codes by comparing the temperature structures, spectral energy distributions, scattered light images, and linear polarisation maps that each model predicts for a variety of disc opacities and viewing angles. The test cases have been chosen to be numerically challenging, with midplane optical depths up 10^6, a sharp density transition at the inner edge and complex scattering matrices. We also review recent progress in the implementation of the Monte Carlo method that allow an efficient solution to these kinds of problems and discuss the advantages and limitations of Monte Carlo codes compared to those of discrete ordinate codes. For each of the test cases, the predicted results from the radiative transfer codes are within good agreement. The results indicate that these codes can be confidently used to interpret present and future observations of protoplanetary discs.

preprint2009arXiv

C2D Spitzer-IRS spectra of disks around T Tauri stars. IV. Crystalline silicates

Dust grains in the planet forming regions around young stars are expected to be heavily processed due to coagulation, fragmentation and crystallization. This paper focuses on the crystalline silicate dust grains in protoplanetary disks. As part of the Cores to Disks Legacy Program, we obtained more than a hundred Spitzer/IRS spectra of TTauri stars. More than 3/4 of our objects show at least one crystalline silicate emission feature that can be essentially attributed to Mg-rich silicates. Observational properties of the crystalline features seen at lambda > 20 mu correlate with each other, while they are largely uncorrelated with the properties of the amorphous silicate 10 mu feature. This supports the idea that the IRS spectra essentially probe two independent disk regions: a warm zone (< 1 AU) emitting at lambda ~ 10 mu and a much colder region emitting at lambda > 20 mu (< 10 AU). We identify a crystallinity paradox, as the long-wavelength crystalline silicate features are 3.5 times more frequently detected (~55 % vs. ~15%) than the crystalline features arising from much warmer disk regions. This suggests that the disk has an inhomogeneous dust composition within ~10 AU. The abundant crystalline silicates found far from their presumed formation regions suggests efficient outward radial transport mechanisms in the disks. The analysis of the shape and strength of both the amorphous 10 mu feature and the crystalline feature around 23 mu provides evidence for the prevalence of micron-sized grains in upper layers of disks. Their presence in disk atmospheres suggests efficient vertical diffusion, likely accompanied by grain-grain fragmentation to balance the efficient growth expected. Finally, the depletion of submicron-sized grains points toward removal mechanisms such as stellar winds or radiation pressure.

preprint2009arXiv

Episodic formation of cometary material in the outburst of a solar-like young star

Our Solar System originated in interstellar gas and dust; the latter is in the form of amorphous silicate particles and carbonaceous dust. The composition of cometary material shows that a significant fraction of the amorphous silicates was transformed into crystalline form during the early evolution of the protosolar nebula. How and when this transformation happened has been controversial, with the main options being heating by the young Sun or shock heating. Here we report mid-infrared features in the outburst spectrum of the young solar-like star EX Lupi that were not present in quiescence. We attribute them to crystalline forsterite; the crystals were produced via thermal annealing in the surface layer of the inner disk by heat from the outburst, a process that has hitherto not been considered. The observed lack of cold crystals excludes shock heating at larger radii.

preprint2009arXiv

Evidence for dust clearing through resolved submillimeter imaging

Mid-infrared spectrophotometric observations have revealed a small sub-class of circumstellar disks with spectral energy distributions (SEDs) suggestive of large inner gaps with low dust content. However, such data provide only an indirect and model-dependent method of finding central holes. Imaging of protoplanetry disks provides an independent check of SED modeling. We present here the direct characterization of three 33-47 AU radii inner gaps, in the disks around LkHa 330, SR 21N and HD 135344B, via 340 GHz (880 micron) dust continuum aperture synthesis observations obtained with the Submillimeter Array (SMA). The large gaps are fully resolved at ~0\farcs3 by the SMA observations and mostly empty of dust, with less than 1 - 7.5 x 10^-6 Msolar of fine grained solids inside the holes. Gas (as traced by atomic accretion markers and CO 4.7 micron rovibrational emission) is still present in the inner regions of all three disks. For each, the inner hole exhibits a relatively steep rise in dust emission to the outer disk, a feature more likely to originate from the gravitational influence of a companion body than from a process expected to show a more shallow gradient like grain growth. Importantly, the good agreement of the spatially resolved data and spectrophotometry-based models lends confidence to current interpretations of SEDs, wherein the significant dust emission deficits arise from disks with inner gaps or holes. Further SED-based searches can therefore be expected to yield numerous additional candidates that can be examined at high spatial resolution.

preprint2009arXiv

Radiative transfer models of mid-infrared H2O lines in the Planet-forming Region of Circumstellar Disks

The study of warm molecular gas in the inner regions of protoplanetary disks is of key importance for the study of planet formation and especially for the transport of H2O and organic molecules to the surfaces of rocky planets/satellites. Recent Spitzer observations have shown that the mid-infrared spectra of protoplanetary disks are covered in emission lines due to water and other molecules. Here, we present a non-LTE 2D radiative transfer model of water lines in the 10-36 mum range that can be used to constrain the abundance structure of water vapor, given an observed spectrum, and show that an assumption of local thermodynamic equilibrium (LTE) does not accurately estimate the physical conditions of the water vapor emission zones. By applying the model to published Spitzer spectra we find that: 1) most water lines are subthermally excited, 2) the gas-to-dust ratio must be one to two orders of magnitude higher than the canonical interstellar medium ratio of 100-200, and 3) the gas temperature must be higher than the dust temperature, and 4) the water vapor abundance in the disk surface must be truncated beyond ~ 1 AU. A low efficiency of water formation below ~ 300 K may naturally result in a lower water abundance beyond a certain radius. However, we find that chemistry, may not be sufficient to produce an abundance drop of many orders of magnitude and speculate that the depletion may also be caused by vertical turbulent diffusion of water vapor from the superheated surface to regions below the snow line, where the water can freeze out and be transported to the midplane as part of the general dust settling. Such a vertical cold finger effect is likely to be efficient due to the lack of a replenishment mechanism of large, water-ice coated dust grains to the disk surface.

preprint2009arXiv

Truncated Disks in TW Hya Association Multiple Star Systems

We present high angular resolution (down to 0.3" = 13 AU in diameter) Submillimeter Array observations of the 880 micron (340 GHz) thermal continuum emission from circumstellar dust disks in the nearby HD 98800 and Hen 3-600 multiple star systems. In both cases, the dust emission is resolved and localized around one stellar component - the HD 98800 B and Hen 3-600 A spectroscopic binaries - with no evidence for circum-system material. Using two-dimensional Monte Carlo radiative transfer calculations, we compare the SMA visibilities and broadband spectral energy distributions with truncated disk models to empirically locate the inner and outer edges of both disks. The HD 98800 B disk appears to be aligned with the spectroscopic binary orbit, is internally truncated at a radius of 3.5 AU, and extends to only 10-15 AU from the central stars. The Hen 3-600 A disk is slightly larger, with an inner edge at 1 AU and an outer radius of 15-25 AU. These inferred disk structures compare favorably with theoretical predictions of their truncation due to tidal interactions with the stellar companions.

preprint2007arXiv

Abundant crystalline silicates in the disk of a very low mass star

We announce the discovery of SST-Lup3-1, a very low mass star close to the brown dwarf boundary in Lupus III with a circum(sub)stellar disk, discovered by the `Cores to Disks' Spitzer Legacy Program from mid-, near-infrared and optical data, with very conspicuous crystalline silicate features in its spectrum. It is the first of such objects with a full 5 to 35 micron spectrum taken with the IRS and it shows strong 10 and 20 micron silicate features with high feature to continuum ratios and clear crystalline features out to 33 micron. The dust in the disk upper layer has a crystalline silicate grain fraction between 15% and 33%, depending on the assumed dust continuum. The availability of the full Spitzer infrared spectrum allows an analysis of the dust composition as a function of temperature and position in the disk. The hot (~ 300 K) dust responsible for the 10 micron feature consists of a roughly equal mix of small (~ 0.1 micron) and large (~ 1.5 micron) grains, whereas the cold (~ 70 K) dust responsible for the longer wavelength silicate features contains primarily large grains (> 1 micron). Since the cold dust emission arises from deeper layers in the inner (< 3 AU) disk as well as from the surface layers of the outer (3-5 AU) disk, this provides direct evidence for combined grain growth and settling in the disk. The inferred crystalline mass fractions in the two components are comparable. Since only the inner 0.02 AU of the disk is warm enough to anneal the amorphous silicate grains, even the lowest fraction of 15% of crystalline material requires either very efficient mixing or other formation mechanisms.

preprint2007arXiv

Cold Disks: Spitzer Spectroscopy of Disks around Young Stars with Large Gaps

We have identified four circumstellar disks with a deficit of dust emission from their inner 15-50 AU. All four stars have F-G spectral type, and were uncovered as part of the Spitzer Space Telescope ``Cores to Disks'' Legacy Program Infrared Spectrograph (IRS) first look survey of ~100 pre-main sequence stars. Modeling of the spectral energy distributions indicates a reduction in dust density by factors of 100-1000 from disk radii between ~0.4 and 15-50 AU, but with massive gas-rich disks at larger radii. This large contrast between the inner and outer disk has led us to use the term `cold disks' to distinguish these unusual systems. However, hot dust [0.02-0.2 Mmoon] is still present close to the central star (R ~0.8 AU). We introduce the 30/13 micron, flux density ratio as a new diagnostic for identifying cold disks. The mechanisms for dust clearing over such large gaps are discussed. Though rare, cold disks are likely in transition from an optically thick to an optically thin state, and so offer excellent laboratories for the study of planet formation.

preprint2006arXiv

Inner Rim of A Molecular Disk Spatially Resolved in Infrared CO Emission Lines

We present high-resolution infrared spectroscopy of the Herbig Ae star HD 141569 A in the CO v=2-1 transition. With the angular resolution attained by the adaptive optics system, the gas disk around HD 141569 A is spatially resolved down to its inner-rim truncation. The size of the inner clearing is 11+-2 AU in radius, close to the gravitational radius of the star. The rough coincidence to the gravitational radius indicates that the viscous accretion working together with the photoevaporation by the stellar radiation has cleared the inner part of the disk.

preprint2006arXiv

Mid-InfraRed imaging of the circumstellar dust around three Herbig Ae stars : HD135344, CQTau, HD163296

Planet formation has been known for many years to be tied to the spatial distribution of gas and dust in disks around young stars. To constrain planet formation models, imaging observations of protoplanetary disks are required. In this framework, we have undertaken a mid-infrared imaging survey of Herbig Ae stars, which are pre-main sequence stars of intermediate mass still surrounded by a large amount of circumstellar material. The observations were made at a wavelength of 20.5 $μ$m with the CAMIRAS camera mounted at the Cassegrain focus of the Canada France Hawaii Telescope. We report the observations of three stars, HD135344, CQTau and HD163296. The circumstellar material around the three objects is spatially resolved. The extensions feature a disk like shape. The images provide direct information on two key parameters of the disk : its inclination and its outer radius. The outer radius is found to be quite different from the one deduced from disk models only constrained by fitting the Spectral Energy Distribution of the object. Other parameters of the disk, such as flaring, dust mass have been deduced from fitting both the observed extension and the spectral energy distribution with sophisticated disk models. Our results show how important imaging data are to tighten constraints on the disk model parameters.

preprint2006arXiv

Probing Protoplanetary Disks with Silicate Emission: Where is the silicate emission zone?

Recent results indicate that the grain size and crystallinity inferred from observations of silicate features may be correlated with spectral type of the central star and/or disk geometry. In this paper, we show that grain size, as probed by the 10 um silicate feature peak-to-continuum and 11.3-to-9.8 um flux ratios, is inversely proportional to log L_star. These trends can be understood using a simple two-layer disk model for passive irradiated flaring disks, CGPLUS. We find that the radius, R_10, of the 10 um silicate emission zone in the disk goes as (L_star/L_sun)^0.56, with slight variations depending on disk geometry (flaring angle, inner disk radius). The observed correlations, combined with simulated emission spectra of olivine and pyroxene mixtures, imply a grain size dependence on luminosity. Combined with the fact that R_10 is smaller for less luminous stars, this implies that the apparent grain size of the emitting dust is larger for low-luminositysources. In contrast, our models suggest that the crystallinity is only marginally affected, because for increasing luminosity, the zone for thermal annealing (assumed to be at T>800 K) is enlarged by roughly the same factor as the silicate emission zone. The observed crystallinity is affected by disk geometry, however, with increased crystallinity in flat disks. The apparent crystallinity may also increase with grain growth due to a corresponding increase in contrast between crystalline and amorphous silicate emission bands.

preprint2005arXiv

Protostellar holes: Spitzer Space Telescope observations of the protostellar binary IRAS16293-2422

Mid-infrared (23-35 micron) emission from the deeply embedded "Class 0" protostar IRAS16293-2422 is detected with the Spitzer Space Telescope infrared spectrograph. A detailed radiative transfer model reproducing the full spectral energy distribution (SED) from 23 micron to 1.3 mm requires a large inner cavity of radius 600 AU in the envelope to avoid quenching the emission from the central sources. This is consistent with a previous suggestion based on high angular resolution millimeter interferometric data. An alternative interpretation using a 2D model of the envelope with an outflow cavity can reproduce the SED but not the interferometer visibilities. The cavity size is comparable to the centrifugal radius of the envelope and therefore appears to be a natural consequence of the rotation of the protostellar core, which has also caused the fragmentation leading to the central protostellar binary. With a large cavity such as required by the data, the average temperature at a given radius does not increase above 60-80 K and although hot spots with higher temperatures may be present close to each protostar, these constitute a small fraction of the material in the inner envelope. The proposed cavity will also have consequences for the interpretation of molecular line data, especially of complex species probing high temperatures in the inner regions of the envelope.

preprint2004arXiv

Dust coagulation in protoplanetary disks: a rapid depletion of small grains

We model the process of dust coagulation in protoplanetary disks and calculate how it affects their observational appearance. Our model involves the detailed solution of the coagulation equation at every location in the disk. At regular time intervals we feed the resulting 3-D dust distribution functions into a continuum radiative transfer code to obtain spectral energy distributions. We find that, even if only the very basic -- and well understood -- coagulation mechanisms are included, the process of grain growth is much too quick to be consistent with infrared observations of T Tauri disks. Small grains are removed so efficiently that, long before the disk reaches an age of 1E6 years typical of T Tauri stars, the SED shows only very weak infrared excess. This is inconsistent with observed SEDs of most classical T Tauri stars. Small grains somehow need to be replenished, for instance by aggregate fragmentation through high-speed collisions. A very simplified calculation shows that when aggregate fragmentation is included, a quasi-stationary grain size distribution is obtained in which growth and fragmentation are in equilibrium. This quasi-stationary state may last 1E6 years or even longer, dependent on the circumstances in the disk, and may bring the time scales into the right regime. We conclude that a simple evolutionary scenario in which grains slowly grow from pristine 0.1 micron grains to larger grains over a period of a few Myr is most likely incorrect.

preprint2004arXiv

Flaring vs. self-shadowed disks: the SEDs of Herbig Ae/Be stars

Isolated Herbig Ae stars can be divided into two groups (Meeus et al. 2001): those with an almost flat spectral energy distribution in the mid-infrared (`group I'), and those with a strong decline towards the far-infrared (`group II'). In this paper we show that the group I vs. II distinction can be understood as arising from flaring vs. self-shadowed disks. We show that these two types of disks are natural solutions of the 2-D radiation-hydrostatic structure equations. Disks with high optical depth turn out to be flaring and have a strong far-IR emission, while disks with an optical depth below a certain threshold drop into the shadow of their own puffed-up inner rim and are weak in the far-IR. In spite of not having a directly irradiated surface layer, self-shadowed disks still display dust features in emission, in agreement with observations of group II sources. We propose an evolutionary scenario in which a disk starts out with a flaring shape (group I source), and then goes through the process of grain growth, causing the optical depth of the disk to drop and the disk to become self-shadowed (group II source). We show that this scenario predicts that the (sub-)millimeter slope of the disk changes from steep (small grains) to Rayleigh-Jeans-like (large grains) in the early stages of evolution, so that all group II sources are expected to have Rayleigh-Jeans-like slopes, while some group I sources may still have steep (sub-)millimeter slopes.

preprint2003arXiv

The effect of scattering on the structure and SED of protoplanetary disks

In this paper we investigate how the inclusion of scattering of the stellar radiation into a passive flaring disk model affects its structure and spectral energy distribution, and whether neglecting it could significantly decrease the model reliability. In order to address these questions we construct a detailed 1+1D vertical structure model in which the scattering properties of the dust can be varied. Models are presented with and without dust scattering, and for different albedos and phase functions. It is found that scattering has the effect of reducing the disk temperature at all heights, so that the disk "shrinks", i.e., the the density at all intermediate heights decreases. However, this effect in most cases is more than compensated by the increase of the total extinction (absorption + scattering) cross section, so that the surface scale height increases, and images in scattered light will see a slightly thicker disk. The integrated infrared emission decreases as the albedo increases, because an increasing part of the flux captured by the disk is reflected away instead of absorbed and reprocessed. The reduction of the infrared thermal emission of the disk is stronger at short wavelengths (near infrared) and practically negligible at millimeter wavelengths. For relatively low albedo (alb <~ 0.5), or for strongly forward-peaked scattering (g roughly >0.8), the infrared flux reduction is relatively small.

C. P. Dullemond

What is connected

Connect this record

See the researcher in context

Building this map preview

55 published item(s)

Spiral arms in the proto-planetary disc HD100453 detected with ALMA: evidence for binary-disc interaction and a vertical temperature gradient

Investigating dust trapping in transition disks with millimeter-wave polarization

Shadows and spirals in the protoplanetary disk HD 100453

A multi-wavelength analysis for interferometric (sub-)mm observations of protoplanetary disks: radial constraints on the dust properties and the disk structure

Scattered light images of spiral arms in marginally gravitationally unstable discs with an embedded planet

Interferometer predictions with triangulated images: solving the multi-scale problem

Mean gas opacity for circumstellar environments and equilibrium temperature degeneracy

Millimetre spectral indices of transition disks and their relation to the cavity radius

Planet-vortex interaction:How a vortex can shepherd a planetary embryo

Radiation hydrodynamics including irradiation and adaptive mesh refinement with AZEuS. I. Methods

A quantification of hydrodynamical effects on protoplanetary dust growth

Asymmetric transition disks: Vorticity or eccentricity?

Explaining millimeter-sized particles in brown dwarf disks

Herbig Stars' Near-Infrared Excess: An Origin in the Protostellar Disk's Magnetically-Supported Atmosphere

Lopsided dust rings in transition disks

Planetesimal formation via sweep-up growth at the inner edge of dead zones

Kinematics of the CO Gas in the Inner Regions of the TW Hya Disk

Mid-Infrared Spectral Variability Atlas of Young Stellar Objects

Accretion through the inner hole of transitional disks: What happens to the dust?

Near-infrared spectroscopy of EX Lupi in outburst

Possible planet-forming regions on submillimetre images

Resolved Images of Large Cavities in Protoplanetary Transition Disks

Spectral signatures of disk eccentricity in young binary systems: I. Circumprimary case

The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? III. Sedimentation driven coagulation inside the snow-line

Trapping dust particles in the outer regions of protoplanetary disks

Warm gas at 50 AU in the disk around Herbig Be star HD 100546

A new condition for the transition from runaway to oligarchic growth

Accretion among preplanetary bodies: the many faces of runaway growth

Dust size distributions in coagulation/fragmentation equilibrium: Numerical solutions and analytical fits

Dust, Ice and Gas in Time (DIGIT) Herschel program first results: A full PACS-SED scan of the gas line emission in protostar DK Cha

Fundamental Vibrational Transition of CO during the Outburst of EX Lupi in 2008

Gas- and dust evolution in protoplanetary disks

High-resolution spectroscopic view of planet formation sites

Impact of grain evolution on the chemical structure of protoplanetary disks

Protoplanetary Disk Structures in Ophiuchus II: Extension to Fainter Sources

Sub-Keplerian accretion onto circumstellar disks

Testing the theory of grain growth and fragmentation by millimeter observations of protoplanetary disks

The inner regions of protoplanetary disks

The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier

The thermal structure and the location of the snow line in the protosolar nebula: axisymmetric models with full 3-D radiative transfer

Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation

C2D Spitzer-IRS spectra of disks around T Tauri stars. IV. Crystalline silicates

Episodic formation of cometary material in the outburst of a solar-like young star

Evidence for dust clearing through resolved submillimeter imaging

Radiative transfer models of mid-infrared H2O lines in the Planet-forming Region of Circumstellar Disks

Truncated Disks in TW Hya Association Multiple Star Systems

Abundant crystalline silicates in the disk of a very low mass star

Cold Disks: Spitzer Spectroscopy of Disks around Young Stars with Large Gaps

Inner Rim of A Molecular Disk Spatially Resolved in Infrared CO Emission Lines

Mid-InfraRed imaging of the circumstellar dust around three Herbig Ae stars : HD135344, CQTau, HD163296

Probing Protoplanetary Disks with Silicate Emission: Where is the silicate emission zone?

Protostellar holes: Spitzer Space Telescope observations of the protostellar binary IRAS16293-2422

Dust coagulation in protoplanetary disks: a rapid depletion of small grains

Flaring vs. self-shadowed disks: the SEDs of Herbig Ae/Be stars

The effect of scattering on the structure and SED of protoplanetary disks