Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We investigate the proto-atmospheres on low-mass planets with eccentric orbits with radiation-hydrodynamics simulations. A 2D radiative model of the proto-atmosphere is established with tabulated opacities for the gas and dust. The solutions reveal large-scale gas recycling inside a bow shock structure. The atmospheres on eccentric planets are typically three to four orders of magnitude less massive than those of planets with circular orbits. Overall, however, a supersonic environment is favorable for planets to keep an early stable atmosphere, rather than harmful, due to the steady gas supply through the recycling flow. We also quantitatively explore how such atmospheres are affected by the relative velocity of the planet to the gas, the planet mass, and the background gas density. Our time-dependent simulations track the orbital evolution of the proto-atmosphere with the planet--disk parameters changing throughout the orbit. Atmospheric properties show oscillatory patterns as the planet travels on an eccentric orbit, with a lag in phase. To sum up, low-mass eccentric planets can retain small proto-atmospheres despite the stripping effects of bow shocks. The atmospheres are always connected to and interacting with the disk gas. These findings provide important insights into the impacts of migration and scattering on planetary proto-atmospheres.

When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $\alpha/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($\alpha/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $\alpha$ turbulent parameter ($\alpha \sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.

Recent ALMA observations revealed concentric annular structures in several young class-II objects. In an attempt to produce the rings and gaps in some of these systems, they have been modeled numerically with a single embedded planet assuming a locally isothermal equation of state. This is often justified by observations targeting the irradiation-dominated outer regions of disks (approximately 100 au). We test this assumption by conducting hydrodynamics simulations of embedded planets in thin locally isothermal and radiative disks that mimic the systems HD 163296 and AS 209 in order to examine the effect of including the energy equation in a seemingly locally isothermal environment as far as planet-disk interaction is concerned. We find that modeling such disks with an ideal equation of state makes a difference in terms of the number of produced rings and the spiral arm contrast in the disk. Locally isothermal disks produce sharper annular or azimuthal features and overestimate a single planet's gap-opening capabilities by producing multiple gaps. In contrast, planets in radiative disks carve a single gap for typical disk parameters. Consequently, for accurate modeling of planets with semimajor axes up to about 100 au, radiative effects should be taken into account even in seemingly locally isothermal disks. In addition, for the case of AS 209, we find that the primary gap is significantly different between locally isothermal and radiative models. Our results suggest that multiple planets are required to explain the ring-rich structures in such systems.

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $\alpha$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."

Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present the phenomena emerging in the framework of two-dimensional axisymmetric, non-ideal magnetohydrodynamic simulations including EUV-/ X-ray driven photoevaporation. Of particular interest are the examination of the transition region between photoevaporation and magnetically driven wind, the possibility of emerging magneto-centrifugal wind effects, as well as the morphology of the wind itself depending on the strength of the magnetic field. We use the PLUTO code in a 2.5D axisymmetric configuration with additional treatment of EUV-/ X-ray heating and dynamic ohmic diffusion based on a semi-analytical chemical model. We identify the transition between both outflow types to occur for values of the initial plasma beta $\beta \geq 10^7$, while magnetically driven winds generally outperform photoevaporation for stronger fields. In our simulations we observe irregular and asymmetric outflows for stronger magnetic fields. In the weak field regime the photoevaporation rates are slightly lowered by perturbations of the gas density in the inner regions of the disk. Overall, our results predict a wind with a lever arm smaller than 1.5, consistent with a hot magneto-thermal wind. Stronger accretion flows are present for values of $\beta < 10^7$.

Knowledge of the midplane temperature of protoplanetary disks is one of the key ingredients in theories of dust growth and planet formation. However, direct measurement of this quantity is complicated, and often depends on the fitting of complex models to the data. In this paper we demonstrate a method to directly measure the midplane gas temperature from an optically thick molecular line, if the disk is moderately inclined. The only model assumption that enters is that the line is very optically thick, also in the midplane region where we wish to measure the temperature. Freeze-out of the molecule onto dust grains could thwart this. However, in regions that are expected to be warm enough to avoid freeze-out, this method should work. We apply the method to the CO 2-1 line channel maps of the disk around HD 163296. We find that the midplane temperature between 100 and 400 au drops only mildly from 25 K down to 18 K. While we see no direct evidence of the midplane being optically thin due to strong CO depletion by freeze-out, we cannot rule it out either. The fact that the inferred temperatures are close to the expected CO freeze-out temperature could be an indication of this. Incidently, for the disk around HD 163296 we also find dynamic evidence for a rather abrupt outer edge of the disk, suggestive of outside-in photoevaporation or truncation by an unseen companion.

Both observations of arc-like structures and luminosity bursts of stars > 1 Myr in age indicate that at least some stars undergo late infall events. We investigate scenarios of replenishing the mass reservoir around a star via capturing and infalling events of cloudlets. We carry out altogether 24 three-dimensional hydrodynamical simulations of cloudlet encounters with a Herbig star of mass 2.5 solar mass using the moving mesh code AREPO. To account for the two possibilities of a star or a cloudlet traveling through the interstellar medium (ISM), we put either the star or the cloudlet at rest with respect to the background gas. For absent cooling in the adiabatic runs, almost none of the cloudlet gas is captured due to high thermal pressure. However, second-generation disks easily form when accounting for cooling of the gas. The disk radii range from several 100 au to about 1000 au and associated arc-like structures up to 10 000 au in length form around the star for runs with and without stellar irradiation. Consistent with angular momentum conservation, the arcs and disks are larger for larger impact parameters. Accounting for turbulence in the cloudlet only mildly changes the model outcome. In the case of the star being at rest with the background gas, the disk formation and mass replenishment process is more pronounced and the associated arc-shaped streamers are longer-lived. The results of our models confirm that late encounter events lead to the formation of transitional disks associated with arc-shaped structures such as observed for AB Aurigae or HD 100546. In addition, we find that second-generation disks and their associated filamentary arms are longer lived (>100 000 yrs) in infall events, when the star is at rest with the background gas.

After its formation, a young star spends some time traversing the molecular cloud complex in which it was born. It is therefore not unlikely that, well after the initial cloud collapse event which produced the star, it will encounter one or more low mass cloud fragments, which we call "cloudlets" to distinguish them from full-fledged molecular clouds. Some of this cloudlet material may accrete onto the star+disk system, while other material may fly by in a hyperbolic orbit. In contrast to the original cloud collapse event, this process will be a "cloudlet flyby" and/or "cloudlet capture" event: A Bondi-Hoyle-Lyttleton type accretion event, driven by the relative velocity between the star and the cloudlet. As we will show in this paper, if the cloudlet is small enough and has an impact parameter similar or less than $GM_{*}/v^2_\infty$ (with $v_\infty$ being the approach velocity), such a flyby and/or capture event would lead to arc-shaped or tail-shaped reflection nebulosity near the star. Those shapes of reflection nebulosity can be seen around several transitional disks and FU Orionis stars. Although the masses in the those arcs appears to be much less than the disk masses in these sources, we speculate that higher-mass cloudlet capture events may also happen occasionally. If so, they may lead to the tilting of the outer disk, because the newly infalling matter will have an angular momentum orientation entirely unrelated to that of the disk. This may be one possible explanation for the highly warped/tilted inner/outer disk geometries found in several transitional disks. We also speculate that such events, if massive enough, may lead to FU Orionis outbursts.

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.

Planetary migration is a key link between planet formation models and observed exoplanet statistics. So far the theory of migration has focused on the interaction of planets with an inviscid or viscously evolving disk. Turbulent viscosity is thought to be the main driver of disk evolution and is known to affect the migration process. Recently, the topic of wind-driven accretion is experiencing a renaissance, as evidence is mounting that PPDs may be less turbulent than previously thought, and 3-D non-ideal MHD modeling of the wind-launching process is maturing. Aim: We wish to investigate how wind-driven accretion may affect migration. We aim for a qualitative exploration of the main effects, rather than a quantitative prediction. Methods: We perform 2-D hydrodynamic simulations with the FARGO3D code in the $(r,\phi)$-plane. The vertical coordinate and the launching of the wind are not treated explicitly. Instead, the torque of the wind onto the disk is treated using a simple 2-parameter formula treating the wind mass loss rate and the lever arm. Results: We find that the wind-driven accretion process has a different way of replenishing the co-orbital region than the viscous accretion. The former always injects mass from the outer edge of the co-orbital region and removes mass from the inner edge, while the latter injects or removes mass from the co-orbital region depending on the radial density gradients in the disk. The migration behavior can differ very much and under certain conditions it can drive rapid type-III-like outward migration. We derive an analytic expression for the parameters under which this outward migration occurs. Conclusion: If wind-driven accretion plays a role in the secular evolution of PPDs, migration studies have to include this process as well, because it can strongly affect the resulting migration rate and direction.

Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure. Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric nature of icelines in systems with massive planets. Methods: We conduct a parameter study using numerical hydrodynamics simulations where we treat viscous heating, thermal cooling and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters. Results: Both gap opening and shock heating can displace the iceline, with the effects being amplified for massive planets in optically thick disks. The gap region can split an initially hot (T>170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric iceline into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context. Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and in general move and/or reshape the water iceline. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones, "islands" of vapor and ice around spirals which could affect the accretion or growth of icy aggregates.

While numerical simulations have been playing a key role in the studies of planet-disk interaction, testing numerical results against observations has been limited so far. With the two directly imaged protoplanets embedded in its circumstellar disk, PDS 70 offers an ideal testbed for planet-disk interaction studies. Using two-dimensional hydrodynamic simulations we show that the observed features can be well explained with the two planets in formation, providing strong evidence that previously proposed theories of planet-disk interaction are in action, including resonant migration, particle trapping, size segregation, and filtration. Our simulations suggest that the two planets are likely in 2:1 mean motion resonance and can remain dynamically stable over million-year timescales. The growth of the planets at $10^{-8}-10^{-7}~M_{\rm Jup}~{\rm yr}^{-1}$, rates comparable to the estimates from H$\alpha$ observations, does not destabilize the resonant configuration. Large grains are filtered at the gap edge and only small, (sub-)$\mu$m grains can flow to the circumplanetary disks and the inner circumstellar disk. With the sub-millimeter continuum ring observed outward of the two directly imaged planets, PDS 70 provides the first observational evidence of particle filtration by gap-opening planets. The observed sub-millimeter continuum emission at the vicinity of the planets can be reproduced when (sub-)$\mu$m grains survive over multiple circumplanetary disk gas viscous timescales and accumulate therein. One such possibility is if (sub-)$\mu$m grains grow in size and remain trapped in pressure bumps, similar to what we find happening in circumstellar disks. We discuss potential implications to planet formation in the solar system and mature extrasolar planetary systems.

Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of 1.25 mm. This surprising result could either be coincidental or indicate that the optical depth in all of the rings is regulated by the same process. We investigated if ongoing planetesimal formation could explain the "fine-tuned" optical depths in the DSHARP rings by removing dust and transforming it into "invisible" planetesimals. We performed a one-dimensional simulation of dust evolution in the second dust ring of the protoplanetary disk around HD 163296, including radial transport of gas and dust, dust growth and fragmentation, and planetesimal formation via gravitational collapse of sufficiently dense pebble concentrations. We show that planetesimal formation can naturally explain the observed optical depths if streaming instability regulates the midplane dust-to-gas ratio to unity. Furthermore, our simple monodisperse analytical model supports the hypothesis that planetesimal formation in dust rings should universally limit their optical depth to the observed range.

Detection of forming planets means detection of the circumplanetary disk (CPD) in reality, since the planet is still surrounded by a disk at this evolutionary stage. Yet, no comprehensive CPD modeling was done in near-infrared wavelengths, where high contrast imaging is a powerful tool to detect these objects. We combined 3D radiative hydrodynamic simulations of various embedded planets with RADMC-3D radiative transfer post-processing that includes scattering of photons on dust particles. We made synthetic images for VLT NaCo/ERIS in the Ks, L', M' bands as well as examined the spectral energy distributions (SEDs) of disks between 1 $\mu m$ and 10 cm. We found that the observed magnitudes from the planet's vicinity will mostly depend on the CPD parameters, not on the planet's. The CPD is 20-100x brighter than the embedded planet in near-IR. We also show how the CPD parameters, e.g. the dust-to-gas ratio will affect the resulting CPD magnitudes. According to the SEDs, the best contrast ratio between the CPD and circumstellar disks is in sub-mm/radio wavelengths and between 8-33 microns in case if the planet opened a resolvable, deep gap ($\ge 5 \rm{M_{Jup}}$), while the contrast is particularly poor in the near-IR. Hence, to detect the forming planet and its CPD, the best chance today is targeting the sub-mm/radio wavelengths and the 10-micron silicate feature vicinity. In order to estimate the forming planet's mass from the observed brightness, it is necessary to run system specific disk modeling.

ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $\alpha$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $\alpha<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks.

We present an analysis of ALMA 1.25 millimeter continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of $\sim40-60$ mas ($\sim6-7$ au). All three disks feature $m=2$ spiral patterns in conjunction with annular substructures. Gas kinematics established by $^{12}$CO $J=2-1$ observations indicate that the continuum spiral arms are trailing. The arm-interarm intensity contrasts are modest, typically less than 3. The Elias 27 spiral pattern extends throughout much of the disk, and the arms intersect the gap at $R\sim69$ au. The spiral pattern in the IM Lup disk is particularly complex-it extends about halfway radially through the disk, exhibiting pitch angle variations with radius and interarm features that may be part of ring substructures or spiral arm branches. Spiral arms also extend most of the way through the WaOph 6 disk, but the source overall is much more compact than the other two disks. We discuss possible origins for the spiral structures, including gravitational instability and density waves induced by a stellar or planetary companion. Unlike the millimeter continuum counterparts of many of the disks with spiral arms detected in scattered light, these three sources do not feature high-contrast crescent-like asymmetries or large ($R>20$ au) emission cavities. This difference may point to multiple spiral formation mechanisms operating in disks.

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.

To characterize the substructures induced in protoplanetary disks by the interaction between stars in multiple systems, we study the $1.25\,$mm continuum and the $^{12}$CO$(J=2-1)$ spectral line emission of the triple systems HT Lup and AS 205, at scales of $\approx 5\,$au, as part of the 'Disk Substructures at High Angular Resolution Project' (DSHARP). In the continuum emission, we find two symmetric spiral arms in the disk around AS 205 N, with pitch angle of $14^\circ$, while the southern component AS 205 S, itself a spectroscopic binary, is surrounded by a compact inner disk and a bright ring at a radius of $34\,$au. The $^{12}$CO line exhibits clear signatures of tidal interactions, with spiral arms, extended arc-like emission, and high velocity gas, possible evidence of a recent close encounter between the disks in the AS 205 system, as these features are predicted by hydrodynamic simulations of fly-by encounters. In the HT Lup system, we detect continuum emission from all three components. The primary disk, HT Lup A, also shows two-armed symmetric spiral structure with a pitch angle of $4^\circ$, while HT Lup B and C, located at $25$ and $434\,$au in projected separation from HT Lup A, are barely resolved with $\sim5$ and $\sim10\,$au in diameter, respectively. The gas kinematics for the closest pair indicates a different sense of rotation for each disk, which could be explained by either a counter rotation of the two disks in different, close to parallel, planes, or by a projection effect of these disks with a close to $90^\circ$ misalignment between them.

A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this paper we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bumps, and if confirmed, put constraints on the physics of the dust trapping mechanism. We model this process analytically in 1D, assuming axisymmetry. By comparing this model to the data, we find that all rings are consistent with dust trapping. Based on a plausible model of the dust temperature we find that several rings are narrower than the pressure scale height, providing strong evidence for dust trapping. The rings have peak absorption optical depth in the range between 0.2 and 0.5. The dust masses stored in each of these rings is of the order of tens of Earth masses, though much ambiguity remains due to the uncertainty of the dust opacities. The dust rings are dense enough to potentially trigger the streaming instability, but our analysis cannot give proof of this mechanism actually operating. Our results show, however, that the combination of very low alpha_turb << 5e-4 and very large grains a_grain >> 0.1 cm can be excluded by the data for all the rings studied in this paper.

The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very non-trivial task. This is in part due to the uncertainty in the dust opacities, an uncertainty which is not likely to be resolved any time soon. It is also partly due to the fact that, as the DSHARP survey has shown, these disk often contain regions of intermediate to high optical depth, even at millimeter wavelengths and at relatively large radius in the disk. This makes the interpretation challenging, in particular if the grains are large and have a large albedo. On the other hand, the highly structured features seen in the DSHARP survey, of which strong indications were already seen in earlier observations, provide a unique opportunity to study the dust growth and dynamics. To provide continuity within the DSHARP project, its follow-up projects, and projects by other teams interested in these data, we present here the methods and opacity choices used within the DSHARP collaboration to link the measured intensity $I_\nu$ to dust surface density $\Sigma_d$.

We introduce the Disk Substructures at High Angular Resolution Project (DSHARP), one of the initial Large Programs conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal of DSHARP is to find and characterize substructures in the spatial distributions of solid particles for a sample of 20 nearby protoplanetary disks, using very high resolution (0.035 arcsec, or 5 au FWHM) observations of their 240 GHz (1.25 mm) continuum emission. These data provide a first homogeneous look at the small-scale features in disks that are directly relevant to the planet formation process, quantifying their prevalence, morphologies, spatial scales, spacings, symmetry, and amplitudes, for targets with a variety of disk and stellar host properties. We find that these substructures are ubiquitous in this sample of large, bright disks. They are most frequently manifested as concentric, narrow emission rings and depleted gaps, although large-scale spiral patterns and small arc-shaped azimuthal asymmetries are also present in some cases. These substructures are found at a wide range of disk radii (from a few au to more than 100 au), are usually compact ($<$10 au), and show a wide range of amplitudes (brightness contrasts). Here we discuss the motivation for the project, describe the survey design and the sample properties, detail the observations and data calibration, highlight some basic results, and provide a general overview of the key conclusions that are presented in more detail in a series of accompanying articles. The DSHARP data -- including visibilities, images, calibration scripts, and more -- are released for community use at https://almascience.org/alma-data/lp/DSHARP.

We present a detailed analysis of the high-angular resolution (0.''037, corresponding to 5 au) observations of the 1.25 mm continuum and $^{12}$CO $2-1$ emission from the disk around the T Tauri star AS 209. AS 209 hosts one of the most unusual disks from the DSHARP sample, the first high angular resolution ALMA survey of disks (Andrews et al. 2018), as nearly all of the emission can be explained with concentric Gaussian rings. In particular, the dust emission consists of a series of narrow and closely spaced rings in the inner $\sim$60 au, two well-separated bright rings in the outer disk, centered at 74 and 120 au, and at least two fainter emission features at 90 and 130 au. We model the visibilities with a parametric representation of the radial surface brightness profile, consisting of a central core and 7 concentric Gaussian rings. Recent hydro-dynamical simulations of low viscosity disks show that super-Earth planets can produce the multiple gaps seen in AS 209 millimeter continuum emission. The $^{12}$CO line emission is centrally peaked and extends out to $\sim$300 au, much farther than the millimeter dust emission. We find axisymmetric, localized deficits of CO emission around four distinct radii, near 45, 75, 120 and 210 au. The outermost gap is located well beyond the edge of the millimeter dust emission, and therefore cannot be due to dust opacity and must be caused by a genuine CO surface density reduction, due either to chemical effects or depletion of the overall gas content.

ALMA observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project (DSHARP) resolve the dust and gas emission on angular scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this article, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm dust continuum and $^{12}$CO J=2-1 rotational line emission at a spatial resolution of 4 and 10 au, respectively. The continuum observations resolve and allow us to characterize the previously discovered dust rings at radii of 67 and 100 au. They also reveal new small scale structures, such as a dark gap at 10 au, a bright ring at 15 au, a dust crescent at a radius of 55 au, and several fainter azimuthal asymmetries. The observations of the CO and dust emission inform about the vertical structure of the disk and allow us to directly constrain the dust extinction optical depth at the dust rings. Furthermore, the observed asymmetries in the dust continuum emission corroborate to the hypothesis that the complex structure of the HD 163296 disk is the result of the gravitational interaction with yet unseen planets.

We present a detailed analysis of new ALMA observations of the disk around the T-Tauri star HD 143006, which at 46 mas (7.6 au) resolution reveal new substructures in the 1.25 mm continuum emission. The disk resolves into a series of concentric rings and gaps together with a bright arc exterior to the rings that resembles hydrodynamics simulations of a vortex, and a bridge-like feature connecting the two innermost rings. Although our $^{12}$CO observations at similar spatial resolution do not show obvious substructure, they reveal an inner disk depleted of CO emission. From the continuum emission and the CO velocity field we find that the innermost ring has a higher inclination than the outermost rings and the arc. This is evidence for either a small ($\sim8^{\circ}$) or moderate ($\sim41^{\circ}$) misalignment between the inner and outer disk, depending on the specific orientation of the near/far sides of the inner/outer disk. We compare the observed substructures in the ALMA observations with recent scattered light data from VLT/SPHERE of this object. In particular, the location of narrow shadow lanes in the SPHERE image combined with pressure scale height estimates, favor a large misalignment of about $41^{\circ}$. We discuss our findings in the context of a dust-trapping vortex, planet-carved gaps, and a misaligned inner disk due to the presence of an inclined companion to HD 143006.

The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass and 1) the gaseous gap depth/width, and 2) the sub/super-Keplerian motion across the gap. We then compute dust continuum intensity maps at the frequency of the DSHARP observations. We provide the relationship between the planet mass and 1) the depth/width of the gaps at millimeter intensity maps, 2) the gap edge ellipticity and asymmetry, and 3) the position of secondary gaps induced by the planet. With these relationships, we lay out the procedure to constrain the planet mass using gap properties, and study the potential planets in the DSHARP disks. We highlight the excellent agreement between observations and simulations for AS 209 and the detectability of the young Solar System analog. Finally, under the assumption that the detected gaps are induced by young planets, we characterize the young planet population in the planet mass-semimajor axis diagram. We find that the occurrence rate for $>$ 5 $M_J$ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques.

Context: Vertically hydrostatic protoplanetary disk models are based on the assumption that the main heating source, stellar irradiation, does not vary much with time. However, it is known that accreting young stars are variable sources of radiation. This is particularly evident for outbursting sources such as EX Lupi and FU Orionis stars. Aim: We investigate how such outbursts affect the vertical structure of the outer regions of the protoplanetary disk, in particular their appearance in scattered light at optical and near-infrared wavelengths. Methods: We employ the 3D FARGOCA radiation-hydrodynamics code, in polar coordinates, to compute the time-dependent behavior of the axisymmetric disk structure. The outbursting inner disk region is not included explicitly. Instead, its luminosity is added to the stellar luminosity and is thus included in the irradiation of the outer disk regions. For time snapshots of interest we insert the density structure into the RADMC-3D radiative transfer code and compute the appearance of the disk at optical/near-infrared wavelengths. Results: We find that, depending on the amplitude of the outbursts, the vertical structure of the disk can become highly dynamic, featuring circular surface waves of considerable amplitude. These "hills" and "valleys" on the disk's surface show up in the scattered light images as bright and dark concentric rings. Initially these rings are small and act as standing waves, but they subsequently lead to outward propagating waves, like the waves produced by a stone thrown into a pond. These waves continue long after the actual outburst has died out. Conclusions: We propose that some of the multi-ringed structures seen in optical/infrared images of several protoplanetary disks may have their origin in outbursts that occurred decades or centuries ago.

While planet formation is thought to occur early in the history of a protoplanetary disk, the presence of planets embedded in disks, or of other processes driving disk evolution, might be traced from their imprints on the disk structure. We observed the T Tauri star HD 143006, located in the 5-11 Myr-old Upper Sco region, in polarized scattered light with VLT/SPHERE at near-infrared wavelengths, reaching an angular resolution of ~0.037" (~6 au). We obtained two datasets, one with a 145 mas diameter coronagraph, and the other without, enabling us to probe the disk structure down to an angular separation of ~0.06" (~10 au). In our observations, the disk of HD 143006 is clearly resolved up to ~0.5" and shows a clear large-scale asymmetry with the eastern side brighter than the western side. We detect a number of additional features, including two gaps and a ring. The ring shows an overbrightness at a position angle (PA) of ~140 deg, extending over a range in position angle of ~60 deg, and two narrow dark regions. The two narrow dark lanes and the overall large-scale asymmetry are indicative of shadowing effects, likely due to a misaligned inner disk. We demonstrate the remarkable resemblance between the scattered light image of HD 143006 and a model prediction of a warped disk due to an inclined binary companion. The warped disk model, based on the hydrodynamic simulations combined with 3D radiative transfer calculations, reproduces all major morphological features. However, it does not account for the observed overbrightness at PA~140 deg. Shadows have been detected in several protoplanetary disks, suggesting that misalignment in disks is not uncommon. However, the origin of the misalignment is not clear. As-yet-undetected stellar or massive planetary companions could be responsible for them, and naturally account for the presence of depleted inner cavities.

We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.

Models of dust coagulation and subsequent planetesimal formation are usually computed on the backdrop of an already fully formed protoplanetary disk model. At the same time, observational studies suggest that planetesimal formation should start early, possibly even before the protoplanetary disk is fully formed. In this paper, we investigate under which conditions planetesimals already form during the disk buildup stage, in which gas and dust fall onto the disk from its parent molecular cloud. We couple our earlier planetesimal formation model at the water snow line to a simple model of disk formation and evolution. We find that under most conditions planetesimals only form after the buildup stage when the disk becomes less massive and less hot. However, there are parameters for which planetesimals already form during the disk buildup. This occurs when the viscosity driving the disk evolution is intermediate ($\alpha_v \sim 10^{-3}-10^{-2}$) while the turbulent mixing of the dust is reduced compared to that ($\alpha_t \lesssim 10^{-4}$), and with the assumption that water vapor is vertically well-mixed with the gas. Such $\alpha_t \ll \alpha_v$ scenario could be expected for layered accretion, where the gas flow is mostly driven by the active surface layers, while the midplane layers, where most of the dust resides, are quiescent.

Protoplanetary disks often appear as multiple concentric rings in dust continuum emission maps and scattered light images. These features are often associated with possible young planets in these disks. Many non-planetary explanations have also been suggested, including snow lines, dead zones and secular gravitational instabilities in the dust. In this paper we suggest another potential origin. The presence of copious amounts of dust tends to strongly reduce the conductivity of the gas, thereby inhibiting the magneto-rotational instability, and thus reducing the turbulence in the disk. From viscous disk theory it is known that a disk tends to increase its surface density in regions where the viscosity (i.e. turbulence) is low. Local maxima in the gas pressure tend to attract dust through radial drift, increasing the dust content even more. We investigate mathematically if this could potentially lead to a feedback loop in which a perturbation in the dust surface density could perturb the gas surface density, leading to increased dust drift and thus amplification of the dust perturbation and, as a consequence, the gas perturbation. We find that this is indeed possible, even for moderately small dust grain sizes, which drift less efficiently, but which are more likely to affect the gas ionization degree. We speculate that this instability could be triggered by the small dust population initially, and when the local pressure maxima are strong enough, the larger dust grains get trapped and lead to the familiar ring-like shapes. We also discuss the many uncertainties and limitations of this model.

Chondrules are the dominant bulk silicate constituent of chondritic meteorites and originate from highly energetic, local processes during the first million years after the birth of the Sun. So far, an astrophysically consistent chondrule formation scenario, explaining major chemical, isotopic and textural features, remains elusive. Here, we examine the prospect of forming chondrules from planetesimal collisions. We show that intensely melted bodies with interior magma oceans became rapidly chemically equilibrated and physically differentiated. Therefore, collisional interactions among such bodies would have resulted in chondrule-like but basaltic spherules, which are not observed in the meteoritic record. This inconsistency with the expected dynamical interactions hints at an incomplete understanding of the planetary growth regime during the protoplanetary disk phase. To resolve this conundrum, we examine how the observed chemical and isotopic features of chondrules constrain the dynamical environment of accreting chondrite parent bodies by interpreting the meteoritic record as an impact-generated proxy of planetesimals that underwent repeated collision and reaccretion cycles. Using a coupled evolution-collision model we demonstrate that the vast majority of collisional debris feeding the asteroid main belt must be derived from planetesimals which were partially molten at maximum. Therefore, the precursors of chondrite parent bodies either formed primarily small, from sub-canonical aluminum-26 reservoirs, or collisional destruction mechanisms were efficient enough to shatter planetesimals before they reached the magma ocean phase. Finally, we outline the window in parameter space for which chondrule formation from planetesimal collisions can be reconciled with the meteoritic record and how our results can be used to further constrain early solar system dynamics.

Hydrodynamic, non-magnetic instabilities can provide turbulent stress in the regions of protoplanetary discs, where the MRI can not develop. The induced motions influence the grain growth, from which formation of planetesimals begins. Thermal relaxation of the gas constrains origins of the identified hydrodynamic sources of turbulence in discs. We estimate the radiative relaxation timescale of temperature perturbations and study the dependence of this timescale on the perturbation wavelength, the location within the disc, the disc mass, and the dust-to-gas mass ratio. We then apply thermal relaxation criteria to localise modes of the convective overstability, the vertical shear instability, and the zombie vortex instability. Our calculations employed the latest tabulated dust and gas mean opacities and we account for the collisional coupling to the emitting species. The relaxation criterion defines the bulk of a typical T Tauri disc as unstable to the development of linear hydrodynamic instabilities. The midplane is unstable to the convective overstability from at most $2\mbox{ au}$ and up to $40\mbox{ au}$, as well as beyond $140\mbox{ au}$. The vertical shear instability can develop between $15\mbox{ au}$ and $180\mbox{ au}$. The successive generation of (zombie) vortices from a seeded noise can work within the inner $0{.}8\mbox{ au}$. Dynamic disc modelling with the evolution of dust and gas opacities is required to clearly localise the hydrodynamic turbulence, and especially its non-linear phase.

Context. DG Tau is a low-mass pre-main sequence star, whose strongly accreting protoplanetary disk exhibits a so-far enigmatic behavior: its mid-infrared thermal emission is strongly time-variable, even turning the 10 $\mu$m silicate feature from emission to absorption temporarily. Aims. We look for the reason for the spectral variability at high spatial resolution and at multiple epochs. Methods. We study the temporal variability of the mid-infrared interferometric signal, observed with the VLTI/MIDI instrument at six epochs between 2011 and 2014. We fit a geometric disk model to the observed interferometric signal to obtain spatial information about the disk. We also model the mid-infrared spectra by template fitting to characterize the profile and time dependence of the silicate emission. We use physically motivated radiative transfer modeling to interpret the mid-infrared interferometric spectra. Results. The inner disk (r<1-3 au) spectra exhibit a 10 $\mu$m absorption feature related to amorphous silicate grains. The outer disk (r>1-3 au) spectra show a crystalline silicate feature in emission, similar to the spectra of comet Hale-Bopp. The striking difference between the inner and outer disk spectral feature is highly unusual among T Tauri stars. The mid-infrared variability is dominated by the outer disk. The strength of the silicate feature changed by more than a factor of two. Between 2011 and 2014 the half-light radius of the mid-infrared-emitting region decreased from 1.15 to 0.7 au. Conclusions. For the origin of the absorption we discuss four possible explanations: a cold obscuring envelope, an accretion heated inner disk, a temperature inversion on the disk surface and a misaligned inner geometry. The silicate emission in the outer disk can be explained by dusty material high above the disk plane, whose mass can change with time, possibly due to turbulence in the disk.

Context. Ice lines are suggested to play a significant role in grain growth and planetesimal formation in protoplanetary disks. Evaporation fronts directly influence the gas and ice abundances of volatile species in the disk and therefore the coagulation physics and efficiency and the chemical composition of the resulting planetesimals. Aims. In this work we investigate the influence of the existence of the CO ice line on the particle growth and on the distribution of CO in the disk. Methods. We include the possibility of tracking the CO content and/or other volatiles in particles and in the gas in our existing dust coagulation and disk evolution model and developed a method for evaporation and condensation of CO using the Hertz-Knudsen equation. Our model does not include fragmentation, yet, which will be part of further investigations. Results. We find no enhanced grain growth just outside the ice line where the particle size is limited by radial drift. Instead we find a depletion of solid material inside the ice line which is solely due to evaporation of the CO. Such a depression inside the ice line may be observable and may help to quantify the processes described in this work. Furthermore, we find that the viscosity and diffusivity of the gas heavily influence the re-distribution of vaporized CO at the ice line and can lead to an increase in the CO abundance by up to a factors of a few in the region just inside the ice line. Depending on the strength of the gaseous transport mechanisms the position of the ice line in our model can change by up to 10 AU and consequently, the temperature at that location can range from 21 K to 23 K.

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.

We present the polarization observations toward the circumstellar disk around HD 142527 by using Atacama Large Millimeter/submillimeter Array (ALMA) at the frequency of 343 GHz. The beam size is $0.51 " \times 0.44 "$, which corresponds to the spatial resolution of $\sim$ 71 $\times$ 62 AU. The polarized intensity displays a ring-like structure with a peak located on the east side with a polarization fraction of $P= 3.26 \pm 0.02$ %, which is different from the peak of the continuum emission from the northeast region. The polarized intensity is significantly weaker at the peak of the continuum where $P= 0.220 \pm 0.010$ %. The polarization vectors are in the radial direction in the main ring of the polarized intensity, while there are two regions outside at the northwest and northeast areas where the vectors are in the azimuthal direction. If the polarization vectors represent the magnetic field morphology, the polarization vectors indicate the toroidal magnetic field configuration on the main ring and the poloidal fields outside. On the other hand, the flip of the polarization vectors is predicted by the self-scattering of thermal dust emission due to the change of the direction of thermal radiation flux. Therefore, we conclude that self-scattering of thermal dust emission plays a major role in producing polarization at millimeter wavelengths in this protoplanetary disk. Also, this puts a constraint on the maximum grain size to be approximately 150 ${\rm \mu m}$ if we assume compact spherical dust grains.

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

Solving the mystery of the origin of chondrules is one of the most elusive goals in the field of meteoritics. Recently the idea of planet(esimal) collisions releasing splashes of lava droplets, long considered out of favor, has been reconsidered as a possible origin of chondrules by several papers. One of the main problems with this idea is the lack of quantitative and simple models that can be used to test this scenario by directly comparing to the many known observables of chondrules. In Paper I of this series we presented a simple thermal evolution model of a spherically symmetric expanding cloud of molten lava droplets that is assumed to emerge from a collision between two planetesimals. In the present paper, number II of this series, we use this model to calculate whether or not volatile elements such as Na and K will remain abundant in these droplets or whether they will get depleted due to evaporation. The high density of the droplet cloud (e.g. small distance between adjacent droplets) causes the vapor to quickly reach saturation pressure and thus shutting down further evaporation. We show to which extent, and under which conditions, this keeps the abundances of these elements high, as is seen in chondrules. We find that for most parameters of our model (cloud mass, expansion velocity, initial temperature) the volatile elements Mg, Si and Fe remain entirely in the chondrules. The Na and K abundances inside the droplets will initially stay mostly at their initial values due to the saturation of the vapor pressure, but at some point start to drop due to the cloud expansion. However, as soon as the temperature starts to decrease, most or all of the vapor recondenses again. At the end the Na and K elements retain most of their initial abundances, albeit occasionally somewhat reduced, depending on the parameters of the expanding cloud model.

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.

The millimeter-wave polarization of the protoplanetary disk around HL Tau has been interpreted as the emission from elongated dust grains aligned with the magnetic field in the disk. However, the self-scattering of thermal dust emission may also explain the observed millimeter-wave polarization. In this paper, we report a modeling of the millimeter-wave polarization of the HL Tau disk with the self-polarization. Dust grains are assumed to be spherical and to have a power-law size distribution. We change the maximum grain size with a fixed dust composition in a fixed disk model to find the grain size to reproduce the observed signature. We find that the direction of the polarization vectors and the polarization degree can be explained with the self-scattering. Moreover, the polarization degree can be explained only if the maximum grain size is $\sim 150 {\rm~\mu m}$. The obtained grain size from the polarization is different from a size of millimeter or larger that has been previously expected from the spectral index of the dust opacity coefficient if the emission is optically thin. We discuss that porous dust aggregates may solve the inconsistency of the maximum grain size between the two constraints.

To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.

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.

We combined hydrodynamical simulations of planet-disk interactions with dust evolution models that include coagulation and fragmentation of dust grains over a large range of radii and derived observational properties using radiative transfer calculations. We studied the role of the snow line in the survival of the inner disk of transition disks. Inside the snow line, the lack of ice mantles in dust particles decreases the sticking efficiency between grains. As a consequence, particles fragment at lower collision velocities than in regions beyond the snow line. This effect allows small particles to be maintained for up to a few Myrs within the first astronomical unit. These particles are closely coupled to the gas and do not drift significantly with respect to the gas. For lower mass planets (1$M_{\rm{Jup}}$), the pre-transition appearance can be maintained even longer because dust still trickles through the gap created by the planet, moves invisibly and quickly in the form of relatively large grains through the gap, and becomes visible again as it fragments and gets slowed down inside of the snow line. The global study of dust evolution of a disk with an embedded planet, including the changes of the dust aerodynamics near the snow line, can explain the concentration of millimetre-sized particles in the outer disk and the survival of the dust in the inner disk if a large dust trap is present in the outer disk. This behaviour solves the conundrum of the combination of both near-infrared excess and ring-like millimetre emission observed in several transition disks.

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen.

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.

We present a new method to constrain the grain size in protoplanetary disks with polarization observations at millimeter wavelengths. If dust grains are grown to the size comparable to the wavelengths, the dust grains are expected to have a large scattering opacity and thus the continuum emission is expected to be polarized due to self-scattering. We perform 3D radiative transfer calculations to estimate the polarization degree for the protoplanetary disks having radial Gaussian-like dust surface density distributions, which have been recently discovered. The maximum grain size is set to be $100 {\rm~\mu m}$ and the observing wavelength to be 870 ${\rm \mu m}$. We find that the polarization degree is as high as 2.5 % with a subarcsec spatial resolution, which is likely to be detected with near-future ALMA observations. The emission is polarized due to scattering of anisotropic continuum emission. The map of the polarization degree shows a double peaked distribution and the polarization vectors are in the radial direction in the inner ring and in the azimuthal direction in the outer ring. We also find the wavelength dependence of the polarization degree: the polarization degree is the highest if dust grains have a maximum size of $a_{\rm max}\sim\lambda/2\pi$, where $\lambda$ is the observing wavelength. Hence, multi-wave and spatially resolved polarization observations toward protoplanetary disks enable us to put a constraint on the grain size. The constraint on the grain size from polarization observations is independent of or may be even stronger than that from the opacity index.

One of the striking discoveries of protoplanetary disc research in recent years are the spiral arms seen in several transitional discs in polarised scattered light. An interesting interpretation of the observed spiral features is that they are density waves launched by one or more embedded (proto-)planets in the disc. In this paper we investigate whether planets can be held responsible for the excitation mechanism of the observed spirals. We use locally isothermal hydrodynamic simulations as well as analytic formulae to model the spiral waves launched by planets. Then H-band scattered light images are calculated using a 3D continuum radiative transfer code to study the effect of surface density and pressure scale height perturbation on the detectability of the spirals. We find that a relative change of about 3.5 in the surface density is required for the spirals to be detected with current telescopes in the near-infrared for sources at the distance of typical star-forming regions (140pc). This value is a factor of eight higher than what is seen in hydrodynamic simulations. We also find that a relative change of only 0.2 in pressure scale height is sufficient to create detectable signatures under the same conditions. Therefore, we suggest that the spiral arms observed to date in protoplanetary discs are the results of changes in the vertical structure of the disc (e.g. pressure scale height perturbation) instead of surface density perturbations.

Protoplanetary disks are dispersed by viscous evolution and photoevaporation in a few million years; in the interim small, sub-micron sized dust grains must grow and form planets. The time-varying abundance of small grains in an evolving disk directly affects gas heating by far-ultraviolet photons, while dust evolution affects photoevaporation by changing the disk opacity and resulting penetration of FUV photons in the disk. Photoevaporative flows, in turn, selectively carry small dust grains leaving the larger particles---which decouple from the gas---behind in the disk. We study these effects by investigating the evolution of a disk subject to viscosity, photoevaporation by EUV, FUV and X-rays, dust evolution, and radial drift using a 1-D multi-fluid approach (gas + different dust grain sizes) to solve for the evolving surface density distributions. The 1-D evolution is augmented by 1+1D models constructed at each epoch to obtain the instantaneous disk structure and determine photoevaporation rates. The implementation of a dust coagulation/fragmentation model results in a marginal decrease in disk lifetimes when compared to models with no dust evolution; the disk lifetime is thus found to be relatively insensitive to the evolving dust opacity. We find that photoevaporation can cause significant reductions in the gas/dust mass ratio in the planet-forming regions of the disk as it evolves, and may result in a corresponding increase in heavy element abundances relative to hydrogen. We discuss implications for theories of planetesimal formation and giant planet formation, including the formation of gas-poor giants. After gas disk dispersal, $\sim 3\times 10^{-4}$ \ms\ of mass in solids typically remain, comparable to the solids inventory of our solar system.

The formation of chondrules is one of the oldest unsolved mysteries in meteoritics and planet formation. Recently an old idea has been revived: the idea that chondrules form as a result of collisions between planetesimals in which the ejected molten material forms small droplets which solidify to become chondrules. Pre-melting of the planetesimals by radioactive decay of 26Al would help producing sprays of melt even at relatively low impact velocity. In this paper we study the radiative cooling of a ballistically expanding spherical cloud of chondrule droplets ejected from the impact site. We present results from a numerical radiative transfer models as well as analytic approximate solutions. We find that the temperature after the start of the expansion of the cloud remains constant for a time t_cool and then drops with time t approximately as T ~ T_0[(3/5)t/t_cool+ 2/5]^(-5/3) for t>t_cool. The time at which this temperature drop starts t_cool depends via an analytical formula on the mass of the cloud, the expansion velocity and the size of the chondrule. During the early isothermal expansion phase the density is still so high that we expect the vapor of volatile elements to saturate so that no large volatile losses are expected.

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.

Streaming instability can be a very efficient way of overcoming growth and drift barriers to planetesimal formation. However, it was shown that strong clumping, which leads to planetesimal formation, requires a considerable number of large grains. State-of-the-art streaming instability models do not take into account realistic size distributions resulting from the collisional evolution of dust. We investigate whether a sufficient quantity of large aggregates can be produced by sticking and what the interplay of dust coagulation and planetesimal formation is. We develop a semi-analytical prescription of planetesimal formation by streaming instability and implement it in our dust coagulation code based on the Monte Carlo algorithm with the representative particles approach. We find that planetesimal formation by streaming instability may preferentially work outside the snow line, where sticky icy aggregates are present. The efficiency of the process depends strongly on local dust abundance and radial pressure gradient, and requires a super-solar metallicity. If planetesimal formation is possible, the dust coagulation and settling typically need ~100 orbits to produce sufficiently large and settled grains and planetesimal formation lasts another ~1000 orbits. We present a simple analytical model that computes the amount of dust that can be turned into planetesimals given the parameters of the disk model.

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.

In recent years many models of chondrule formation have been proposed. One of those models is the processing of dust in shock waves in protoplanetary disks. In this model, the dust and the chondrule precursors are overrun by shock waves, which heat them up by frictional heating and thermal exchange with the gas. In this paper we reanalyze the nebular shock model of chondrule formation and focus on the downstream boundary condition. We show that for large-scale plane-parallel chondrule-melting shocks the postshock equilibrium temperature is too high to avoid volatile loss. Even if we include radiative cooling in lateral directions out of the disk plane into our model (thereby breaking strict plane-parallel geometry) we find that for a realistic vertical extent of the solar nebula disk the temperature decline is not fast enough. On the other hand, if we assume that the shock is entirely optically thin so that particles can radiate freely, the cooling rates are too high to produce the observed chondrules textures. Global nebular shocks are therefore problematic as the primary sources of chondrules.

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.

Simple toy models are often not sufficient to cover the complexity of the dust coagulation process, and a number of numerical approaches are therefore used, among which integration of the Smoluchowski equation and various versions of Monte Carlo algorithm are the most popular. In this paper, we directly compare the Smoluchowski and Monte Carlo approaches and we find a general agreement for most of the coagulation problems. However, for the sweep-up growth driven by the "lucky" breakthrough mechanism, the methods exhibit very different resolution dependencies. With too few mass bins, the Smoluchowski algorithm tends to overestimate the growth rate and the probability of breakthrough. The Monte Carlo method is less resolution dependent in the growth timescale aspect but it tends to underestimate the breakthrough chance due to its limited dynamic mass range. We discuss the features and drawbacks of both the approaches, which may limit their astrophysical applications.

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

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.

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.

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.

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.

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.

The statistics of discovered exoplanets suggest that planets form efficiently. However, there are fundamental unsolved problems, such as excessive inward drift of particles in protoplanetary disks during planet formation. Recent theories invoke dust traps to overcome this problem. We report the detection of a dust trap in the disk around the star Oph IRS 48 using observations from the Atacama Large Millimeter/submillimeter Array (ALMA). The 0.44-millimeter-wavelength continuum map shows high-contrast crescent-shaped emission on one side of the star originating from millimeter-sized grains, whereas both the mid-infrared image (micrometer-sized dust) and the gas traced by the carbon monoxide 6-5 rotational line suggest rings centered on the star. The difference in distribution of big grains versus small grains/gas can be modeled with a vortex-shaped dust trap triggered by a companion.

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 $\alpha_{\mathrm{turb}}= 10^{-3}$, and average fragmentation velocities for ices $v_f = 10 m s^{-1}$.

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.

We present a possible identification strategy for first hydrostatic core (FHSC) candidates and make predictions of ALMA dust continuum emission maps from these objects. We analyze the results given by the different bands and array configurations and identify which combinations of the two represent our best chance of solving the fragmentation issue in these objects. If the magnetic field is playing a role, the emission pattern will show evidence of a pseudo-disk and even of a magnetically driven outflow, which pure hydrodynamical calculations cannot reproduce.

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\'eclet 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.

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

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.

It is unknown how far dust growth can proceed by coagulation. Obstacles to collisional growth are the fragmentation and bouncing barriers. However, in all previous simulations of the dust-size evolution in protoplanetary disks, only the mean collision velocity has been considered, neglecting that a small but possibly important fraction of the collisions will occur at both much lower and higher velocities. We study the effect of the probability distribution of impact velocities on the collisional dust growth barriers. Assuming a Maxwellian velocity distribution for colliding particles to determine the fraction of sticking, bouncing, and fragmentation, we implement this in a dust-size evolution code. We also calculate the probability of growing through the barriers and the growth timescale in these regimes. We find that the collisional growth barriers are not as sharp as previously thought. With the existence of low-velocity collisions, a small fraction of the particles manage to grow to masses orders of magnitude above the main population. A particle velocity distribution softens the fragmentation barrier and removes the bouncing barrier. It broadens the size distribution in a natural way, allowing the largest particles to become the first seeds that initiate sweep-up growth towards planetesimal sizes.

The low-mass star formation evolutionary sequence is relatively well-defined both from observations and theoretical considerations. The first hydrostatic core is the first protostellar equilibrium object that is formed during the star formation process. Using state-of-the-art radiation-magneto-hydrodynamic 3D adaptive mesh refinement calculations, we aim to provide predictions for the dust continuum emission from first hydrostatic cores. We investigate the collapse and the fragmentation of magnetized one solar mass prestellar dense cores and the formation and evolution of first hydrostatic cores using the RAMSES code. We use three different magnetization levels for the initial conditions, which cover a large variety of early evolutionary morphology, e.g., the formation of a disk or a pseudo-disk, outflow launching, and fragmentation. We post-process the dynamical calculations using the 3D radiative transfer code RADMC-3D. We compute spectral energy distributions and usual evolutionary stage indicators such as bolometric luminosity and temperature. We find that the first hydrostatic core lifetimes depend strongly on the initial magnetization level of the parent dense core. We derive, for the first time, spectral energy distribution evolutionary sequences from high-resolution radiation-magneto-hydrodynamic calculations. We show that under certain conditions, first hydrostatic cores can be identified from dust continuum emission at 24 microns and 70 microns. We also show that single spectral energy distributions cannot help to distinguish between the formation scenarios of the first hydrostatic core, i.e., between the magnetized and non-magnetized models. Spectral energy distributions are a first useful and direct way to target first hydrostatic core candidates but high-resolution interferometry is definitively needed to determine the evolutionary stage of the observed sources.

Hydrogen recombination lines are one of the major diagnostics of H II region physical properties and kinematics. In the near future, the Expanded Very Large Array (EVLA) and the Atacama Large Millimeter Array (ALMA) will allow observers to study recombination lines in the radio and sub-mm regime in unprecedented detail. In this paper, we study the properties of recombination lines, in particular at ALMA wavelengths. We find that such lines will lie in almost every wideband ALMA setup and that the line emission will be equally detectable in all bands. Furthermore, we present our implementation of hydrogen recombination lines in the adaptive-mesh radiative transfer code RADMC-3D. We particularly emphasize the importance of non-LTE (local thermodynamical equilibrium) modeling since non-LTE effects can drastically affect the line shapes and produce asymmetric line profiles from radially symmetric H II regions. We demonstrate how these non-LTE effects can be used as a probe of systematic motions (infall & outflow) in the gas. We use RADMC-3D to produce synthetic observations of model H II regions and study the necessary conditions for observing such asymmetric line profiles with ALMA and EVLA.

We explore the accretion mechanisms in EX Lupi, prototype of EXor variables, during its quiescence and outburst phases. We analyse high-resolution optical spectra taken before, during, and after its 2008 outburst. In quiescence and outburst, the star presents many permitted emission lines, including typical CTTS lines and numerous neutral and ionized metallic lines. During the outburst, the number of emission lines increases to over a thousand, with narrow plus broad component structure (NC+BC). The BC profile is highly variable on short timescales (24-72h). An active chromosphere can explain the metallic lines in quiescence and the outburst NC. The dynamics of the BC line profiles suggest an origin in a hot, dense, non-axisymmetric, and non-uniform accretion column that suffers velocity variations along the line-of-sight on timescales of days. Assuming Keplerian rotation, the emitting region would be located at ~0.1-0.2 AU, consistent with the inner disk rim, but the velocity profiles of the lines reveal a combination of rotation and infall. Line ratios of ions and neutrals can be reproduced with a temperature of T~6500 K for electron densities of a few times 10$^{12}$cm$^{-3}$ in the line-emitting region. The data confirm that the 2008 outburst was an episode of increased accretion, albeit much stronger than previous EX Lupi and typical EXors outbursts. The line profiles are consistent with the infall/rotation of a non-axisymmetric structure that could be produced by clumpy accretion during the outburst phase. A strong inner disk wind appears in the epochs of higher accretion. The rapid recovery of the system after the outburst and the similarity between the pre-outburst and post-outburst states suggest that the accretion channels are similar during the whole period, and only the accretion rate varies, providing a superb environment for studying the accretion processes.

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.

We present the analysis of an IRS 5-38 {\mu}m spectrum and MIPS photometric measurements of an infrared echo near the Cassiopeia A supernova remnant observed with the Spitzer Space Telescope. We have modeled the recorded echo accounting for PAHs, quantum-heated carbon and silicate grains, as well as thermal carbon and silicate particles. Using the fact that optical light echo spectroscopy has established that Cas A originated from a type IIb supernova explosion showing an optical spectrum remarkably similar to the prototypical type IIb SN 1993J, we use the latter to construct template data input for our simulations. We are then able to reproduce the recorded infrared echo spectrum by combining the emission of dust heated by the UV burst produced at the shock breakout after the core-collapse and dust heated by optical light emitted near the visual maximum of the supernova light curve, where the UV burst and optical light curve characteristics are based on SN 1993J. We find a mean density of \sim680 H cm^{-3} for the echo region, with a size of a few light years across. We also find evidence of dust processing in the form of a lack of small PAHs with less than \sim300 carbon atoms, consistent with a scenario of PAHs destruction by the UV burst via photodissociation at the estimated distance of the echo region from Cas A. Furthermore, our simulations suggest that the weak 11 {\mu}m features of our recorded infrared echo spectrum are consistent with a strong dehydrogenated state of the PAHs. This exploratory study highlights the potential of investigating dust processing in the interstellar medium through infrared echoes.

The formation of planetesimals is often accredited to collisional sticking of dust grains. The exact process is unknown, as collisions between larger aggregates tend to lead to fragmentation or bouncing rather than sticking. Recent laboratory experiments have however made great progress in the understanding and mapping of the complex physics involved in dust collisions. We want to study the possibility of planetesimal formation using the results from the latest laboratory experiments, particularly by including the fragmentation with mass transfer effect, which might lead to growth even at high impact velocities. We present a new experimentally and physically motivated dust collision model capable of predicting the outcome of a collision between two particles of arbitrary masses and velocities. It is used together with a continuum dust-size evolution code that is both fast in terms of execution time and able to resolve the dust well at all sizes, allowing for all types of interactions to be studied without biases. We find that for the general dust population, bouncing collisions prevent the growth above millimeter-sizes. However, if a small number of cm-sized particles are introduced, for example due to vertical mixing or radial drift, they can act as a catalyst and start to sweep up the smaller particles. At a distance of 3 AU, 100-meter-sized bodies are formed on a timescale of 1 Myr. We conclude that direct growth of planetesimals might be a possibility thanks to a combination of the existence of a bouncing barrier and the fragmentation with mass transfer effect. The bouncing barrier is here even beneficial, as it prevents the growth of too many large particles that would otherwise only fragment among each other, and creates a reservoir of small particles that can be swept up by larger bodies. However, for this process to work, a few seeds of cm in size or larger have to be introduced.

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.

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.

EX Lup is the prototype of the EXor class of eruptive young stars. These objects show optical outbursts which are thought to be related to runaway accretion onto the star. In a previous study we observed in-situ crystal formation in the disk of EX Lup during its latest outburst in 2008, making the object an ideal laboratory to investigate circumstellar crystal formation and transport. This outburst was monitored by a campaign of ground-based and Spitzer Space Telescope observations. Here we modeled the spectral energy distribution of EX Lup in the outburst from optical to millimeter wavelengths with a 2D radiative transfer code. Our results showed that the shape of the SED at optical wavelengths was more consistent with a single temperature blackbody than a temperature distribution. We also found that this single temperature component emitted 80-100 % of the total accretion luminosity. We concluded that a thermal instability, the most widely accepted model of EXor outbursts, was likely not the triggering mechanism of the 2008 outburst of EX Lup. Our mid-infrared Spitzer spectra revealed that the strength of all crystalline bands between 8 and 30 um increased right after the end of the outburst. Six months later, however, the crystallinity in the 10 um silicate feature complex decreased. Our modeling of the mid-infrared spectral evolution of EXLup showed that, although vertical mixing should be stronger during the outburst than in the quiescent phase, fast radial transport of crystals (e.g., by stellar/disk wind) was required to reproduce the observed mid-infrared spectra.

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.

We present Spitzer/IRS spectra of 31 TTS and IRAM/1.3mm observations for 34 low- and intermediate-mass stars in the Cep OB2 region. Including our previously published data, we analyze 56 TTS and the 3 intermediate-mass stars with silicate features in Tr 37 (~4 Myr) and NGC 7160 (~12 Myr). The silicate emission features are well reproduced with a mixture of amorphous (with olivine, forsterite, and silica stoichiometry) and crystalline grains (forsterite, enstatite). We explore grain size and disk structure using radiative transfer disk models, finding that most objects have suffered substantial evolution (grain growth, settling). About half of the disks show inside-out evolution, with either dust-cleared inner holes or a radially-dependent dust distribution, typically with larger grains and more settling in the innermost disk. The typical strong silicate features require nevertheless the presence of small dust grains, and could be explained by differential settling according to grain size, anomalous dust distributions, and/or optically thin dust populations within disk gaps. M-type stars tend to have weaker silicate emission and steeper SEDs than K-type objects. The inferred low dust masses are in a strong contrast with the relatively high gas accretion rates, suggesting global grain growth and/or an anomalous gas to dust ratio. Transition disks (TD) in the Cep OB2 region display strongly processed grains, suggesting that they are dominated by dust evolution and settling. Finally, the presence of rare but remarkable disks with strong accretion at old ages reveals that some very massive disks may still survive to grain growth, gravitational instabilities, and planet formation.

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.

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.

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.

We investigate how the X factor, the ratio of H_2 column density (NH2) to velocity-integrated CO intensity (W), is determined by the physical properties of gas in model molecular clouds (MCs). We perform radiative transfer calculations on chemical-MHD models to compute X. Using integrated NH2 and W reproduces the limited range in X found in observations, resulting in a mean value X=2\times10^20 s/cm^2/K^1/km^1 from the Galactic MC model. However, in limited velocity intervals, X can take on a much larger range due to CO line saturation. Thus, X strongly depends on both the range in gas velocities and volume densities. The temperature (T) variations within individual MCs do not strongly affect X, as dense gas contributes most to setting X. For fixed velocity and density structure, gas with higher T has higher W, yielding X ~ T^-1/2 for T~20-100 K. We demonstrate that the linewidth-size scaling relation does not influence the X factor - only the range in velocities is important. Clouds with larger linewidths, regardless of the linewidth-size relation, have a higher W, corresponding to a lower value of X, scaling roughly as X ~ sigma^-1/2. The "mist" model, consisting of optically thick cloudlets with well-separated velocities, does not accurately reflect the conditions in a turbulent MC. We propose that the observed cloud-average values of X ~ XGal is simply a result of the limited range in NH2, temperatures, and velocities found in Galactic MCs - a ~constant value of X therefore does not require any linewidth-size relation, or that MCs are virialized objects. Since gas properties likely differ (slightly) between clouds, masses derived through a standard X should only be considered as a rough first estimate. For temperatures T~10-20 K, velocity dispersions ~1-6 km/s, and NH2~2-20\times10^21 cm^-2, we find cloud-averaged X ~ 2-4\times10^20 s/cm^2/K^1/km^1 for Solar-metallicity models.

Massive stars form deeply embedded in dense molecular gas, which they stir and heat up and ionize. During an early phase, the ionization is confined to hypercompact HII regions, and the stellar radiation is entirely absorbed by dust, giving rise to a hot molecular core. To investigate the innermost structure of such high-mass star-forming regions, we observed vibrationally excited HCN (via the direct $\ell$-type transition of v2=1, $\Delta J$=0, J=13, which lies 1400 K above ground) toward the massive hot molecular cores G10.47+0.03, SgrB2-N, and SgrB2-M with the Very Large Array (VLA) at 7 mm, reaching a resolution of about 1000 AU (0.1"). We detect the line both in emission and in absorption against HII regions. The latter allows to derive lower limits on the column densities of hot HCN, which are several times $10^{19}$ cm$^{-2}$. We see indication of expansion motions in G10.47+0.03 and detect velocity components in SgrB2-M at 50, 60, and 70 km/s relative to the Local Standard of Rest. The emission originates in regions of less than 0.1 pc diameter around the hypercompact HII regions G10.47+0.03 B1 and SgrB2-N K2, and reaches brightness temperatures of more than 200 K. Using the three-dimensional radiative transfer code RADMC-3D, we model the sources as dense dust cores heated by stars in the HII regions, and derive masses of hot (>300 K) molecular gas of more than 100 solar masses (for an HCN fractional abundance of 10$^{-5}$), challenging current simulations of massive star formation. Heating only by the stars in the HII regions is sufficient to produce such large quantities of hot molecular gas, provided that dust is optically thick to its own radiation, leading to high temperatures through diffusion of radiation.

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.

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.

One of the longstanding unsolved problems of planet formation is how solid bodies of a few decimeters in size can "stick" to form large planetesimals. This is known as the "meter size barrier". In recent years it has become increasingly clear that some form of "particle trapping" must have played a role in overcoming the meter size barrier. Particles can be trapped in long-lived local pressure maxima, such as those in anticyclonic vortices, zonal flows or those believed to occur near ice lines or at dead zone boundaries. Such pressure traps are the ideal sites for the formation of planetesimals and small planetary embryos. Moreover, they likely produce large quantities of such bodies in a small region, making it likely that subsequent N-body evolution may lead to even larger planetary embryos. The goal of this Letter is to show that this indeed happens, and to study how efficient it is. In particular, we wish to find out if rocky/icy bodies as large as 10 M_earth can form within 1 Myr, since such bodies are the precursors of gas giant planets in the core accretion scenario.

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.

Theoretical and observational investigations have indicated that the abundance of carbon monoxide (CO) is very sensitive to intrinsic properties of the gaseous medium, such as density, metallicity, and the background UV field. In order to accurately interpret CO observations, it is thus important to understand how well CO traces the gas, which in molecular clouds (MCs) is predominantly molecular hydrogen (H2). Recent hydrodynamic simulations by Glover & Mac Low have explicitly followed the formation and destruction of molecules in model MCs under varying conditions, confirming that CO formation strongly depends on the cloud properties. Conversely, the H2 formation is primarily determined by the age of the MC. We apply radiative transfer calculations to these MC models in order to investigate the properties of CO line emission. We focus on integrated CO (J=1-0) intensities emerging from individual clouds, including its relationship to the total, H2, and CO column densities, as well as the "X factor," the ratio of H2 column density to CO intensity. Models with high CO abundances have a threshold CO intensity ~65 K km/s at sufficiently large extinctions. Clouds with low CO abundances show no such intensity thresholds. The distribution of H2 column densities are well described as log-normal functions, though the distributions of CO intensities and column densities are usually not log-normal. In general, the PDFs of the integrated intensity do not follow the distribution functions of CO column densities. In the model with Milky Way-like conditions, the X factor is in agreement with the near constant value determined from observations. In clouds with lower CO abundances the X factor can vary appreciably - sometimes by > 4 orders of magnitude. In models with high densities, the CO line is fully saturated, so that the X factor is directly proportional to the molecular column density.

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.

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.

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.

We examine the first phase of the core accretion model, namely the dust growth/fragmentation in binary systems. In our model, a gas and dust disk is present around the primary star and is perturbed by the secondary. We study the effects of a secondary with/without eccentricity on the dust population to determine what sizes the aggregates can reach and how that compares to the dust population in disks around single stars. We find that the secondary star has two effects on the dust population. 1.) The disk is truncated due to the presence of the secondary star and the maximum mass of the particles is decreased in the lowered gas densities. This effect is dominant in the outer disk. 2.) The perturbation of the secondary pumps up the eccentricity of the gas disk, which in turn increases the relative velocity between the dust and the gas. Therefore the maximum particle sizes are further decreased. The second effect of the secondary influences the entire disk. Coagulation is efficiently reduced even at the very inner parts of the disk. The average mass of the particles is reduced by four orders of magnitude (as a consequence, the stopping time is reduced by one order of magnitude) in disks around binary systems compared to dust in disks around single stars.

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.

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 $\gamma$ = 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.

Previous studies suggest that the planet-forming disks around very-low-mass stars/brown dwarfs may be flatter than those around more massive stars, in contrast to model predictions of larger scale heights for gas-disks around lower-mass stars. We conducted a statistically robust study to determine whether there is evidence for stellar-mass-dependent disk structure in planet-forming disks. We find a statistically significant difference in the Spitzer/IRAC color distributions of disks around very-low-mass and low-mass stars all belonging to the same star-forming region, the Chamaeleon I star-forming region. We show that self consistently calculated disk models cannot fit the median spectral energy distributions (SEDs) of the two groups. These SEDs can be only explained by flatter disk models, consistent with the effect of dust settling in disks. We find that relative to the disk structure predicted for flared disks the required reduction in disk scale height is anti-correlated with the stellar mass, i.e. disks around lower-mass stars are flatter. Our results show that the initial and boundary conditions of planet formation are stellar-mass-dependent, an important finding that must be considered in planet formation models.

In this paper we intend to provide an indirect method to detect Jovian planets by studying near infrared emission spectra originating in the protoplanetary disks around T Tauri stars. Our idea is to investigate whether a massive planet could induce any observable effect on the spectral lines emerging in the disks atmosphere. As a tracer molecule we propose CO, which is excited in the ro-vibrational fundamental band in the disk atmosphere to a distance of ~2-3 AU (depending on the stellar mass) where terrestrial planets are thought to form. The synthetic molecular spectral line profiles were calculated by an own developed semi-analytical double layer disk model. 2D gas dynamics were incorporated in the calculation of synthetic spectral lines. We demonstrate that a massive planet embedded in a protoplanetary disk strongly influences the originally circular Keplerian gas dynamics. The perturbed motion of the gas can be detected by comparing the CO line profiles in emission, which is emerge from planet-bearing to those of planet-free disk models. The planet signal has two major characteristics: a permanent line profile asymmetry, and short timescale variability correlated with the orbital phase of the giant planet. We have found that the strength of the asymmetry depends on the physical parameters of the star-planet-disk system, such as the disk inclination angle, the planetary and stellar masses, the orbital distance, and the size of the disk inner cavity. The permanent line profile asymmetry is caused by a disk in an eccentric state in the gap opened by the giant planet. However, the variable component is a consequence of the local dynamical perturbation by the orbiting giant planet. We show that a forming giant planet, still embedded in the protoplanetary disk, can be detected using contemporary or future high-resolution near-IR spectrographs like VLT/CRIRES and ELT/METIS.

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.

(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.

The spatial morphology, spectral characteristics, and time variability of ultracompact H II regions provide strong constraints on the process of massive star formation. We have performed simulations of the gravitational collapse of rotating molecular cloud cores, including treatments of the propagation of ionizing and non-ionizing radiation. We here present synthetic radio continuum observations of H II regions from our collapse simulations, to investigate how well they agree with observation, and what we can learn about how massive star formation proceeds. We find that intermittent shielding by dense filaments in the gravitationally unstable accretion flow around the massive star leads to highly variable H II regions that do not grow monotonically, but rather flicker, growing and shrinking repeatedly. This behavior appears able to resolve the well-known lifetime problem. We find that multiple ionizing sources generally form, resulting in groups of ultracompact H II regions, consistent with observations. We confirm that our model reproduces the qualitative H II region morphologies found in surveys, with generally consistent relative frequencies. We also find that simulated spectral energy distributions (SEDs) from our model are consistent with the range of observed H II region SEDs, including both regions showing a normal transition from optically thick to optically thin emission, and those with intermediate spectral slopes. In our models, anomalous slopes are solely produced by inhomogeneities in the H II region, with no contribution from dust emission at millimeter or submillimeter wavelengths. We conclude that many observed characteristics of ultracompact H II regions appear consistent with massive star formation in fast, gravitationally unstable, accretion flows.

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.

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.

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.

We present a new systematic way of setting up galactic gas disks based on the assumption of detailed hydrodynamic equilibrium. To do this, we need to specify the density distribution and the velocity field which supports the disk. We first show that the required circular velocity has no dependence on the height above or below the midplane so long as the gas pressure is a function of density only. The assumption of disks being very thin enables us to decouple the vertical structure from the radial direction. Based on that, the equation of hydrostatic equilibrium together with the reduced Poisson equation leads to two sets of second-order non-linear differential equation, which are easily integrated to set-up a stable disk. We call one approach `density method' and the other one `potential method'. Gas disks in detailed balance are especially suitable for investigating the onset of the gravitational instability. We revisit the question of global, axisymmetric instability using fully three-dimensional disk simulations. The impact of disk thickness on the disk instability and the formation of spontaneously induced spirals is studied systematically with or without the presence of the stellar potential. In our models, the numerical results show that the threshold value for disk instability is shifted from unity to 0.69 for self-gravitating thick disks and to 0.75 for combined stellar and gas thick disks. The simulations also show that self-induced spirals occur in the correct regions and with the right numbers as predicted by the analytic theory.

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.

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.

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.

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.

A long-standing problem in low-mass star formation is the "luminosity problem," whereby protostars are underluminous compared to the accretion luminosity expected both from theoretical collapse calculations and arguments based on the minimum accretion rate necessary to form a star within the embedded phase duration. Motivated by this luminosity problem, we present a set of evolutionary models describing the collapse of low-mass, dense cores into protostars, using the Young & Evans (2005) model as our starting point. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse in two dimensions. From the resulting spectral energy distributions, we calculate standard observational signatures to directly compare to observations. We incorporate several modifications and additions to the original Young & Evans model in an effort to better match observations with model predictions. We find that scattering, 2-D geometry, mass-loss, and outflow cavities all affect the model predictions, as expected, but none resolve the luminosity problem. A cycle of episodic mass accretion, however, can resolve this problem and bring the model predictions into better agreement with observations. Standard assumptions about the interplay between mass accretion and mass loss in our model give star formation efficiencies consistent with recent observations that compare the core mass function (CMF) and stellar initial mass function (IMF). The combination of outflow cavities and episodic mass accretion reduce the connection between observational Class and physical Stage to the point where neither of the two common observational signatures (bolometric temperature and ratio of bolometric to submillimeter luminosity) can be considered reliable indicators of physical Stage.

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.

The growth processes from protoplanetary dust to planetesimals are not fully understood. Laboratory experiments and theoretical models have shown that collisions among the dust aggregates can lead to sticking, bouncing, and fragmentation. However, no systematic study on the collisional outcome of protoplanetary dust has been performed so far so that a physical model of the dust evolution in protoplanetary disks is still missing. We intend to map the parameter space for the collisional interaction of arbitrarily porous dust aggregates. This parameter space encompasses the dust-aggregate masses, their porosities and the collision velocity. With such a complete mapping of the collisional outcomes of protoplanetary dust aggregates, it will be possible to follow the collisional evolution of dust in a protoplanetary disk environment. We use literature data, perform own laboratory experiments, and apply simple physical models to get a complete picture of the collisional interaction of protoplanetary dust aggregates. In our study, we found four different types of sticking, two types of bouncing, and three types of fragmentation as possible outcomes in collisions among protoplanetary dust aggregates. We distinguish between eight combinations of porosity and mass ratio. For each of these cases, we present a complete collision model for dust-aggregate masses between 10^-12 and 10^2 g and collision velocities in the range 10^-4 to 10^4 cm/s for arbitrary porosities. This model comprises the collisional outcome, the mass(es) of the resulting aggregate(s) and their porosities. We present the first complete collision model for protoplanetary dust. This collision model can be used for the determination of the dust-growth rate in protoplanetary disks.

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.

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.

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.

The material that formed the present-day Solar System originated in feeding zones in the inner Solar Nebula located at distances within ~20 AU from the Sun, known as the planet-forming zone. Meteoritic and cometary material contain abundant evidence for the presence of a rich and active chemistry in the planet-forming zone during the gas-rich phase of Solar System formation. It is a natural conjecture that analogs can be found amoung the zoo of protoplanetary disks around nearby young stars. The study of the chemistry and dynamics of planet formation requires: 1) tracers of dense gas at 100-1000 K and 2) imaging capabilities of such tracers with 5-100 (0.5-20 AU) milli-arcsec resolution, corresponding to the planet-forming zone at the distance of the closest star-forming regions. Recognizing that the rich infrared (2-200 micron) molecular spectrum recently discovered to be common in protoplanetary disks represents such a tracer, we present a new general raytracing code, RADLite, that is optimized for producing infrared line spectra and images from axisymmetric structures. RADLite can consistently deal with a wide range of velocity gradients, such as those typical for the inner regions of protoplanetary disks. The code is intended as a backend for chemical and excitation codes, and can rapidly produce spectra of thousands of lines for grids of models for comparison with observations. Such radiative transfer tools will be crucial for constraining both the structure and chemistry of planet-forming regions, including data from current infrared imaging spectrometers and extending to the Atacama Large Millimeter Array and the next generation of Extremely Large Telescopes, the James Webb Space Telescope and beyond.

Context: Protoplanetary disks are observed to remain dust-rich for up to several million years. Theoretical modeling, on the other hand, raises several questions. Firstly, dust coagulation occurs so rapidly, that if the small dust grains are not replenished by collisional fragmentation of dust aggregates, most disks should be observed to be dust poor, which is not the case. Secondly, if dust aggregates grow to sizes of the order of centimeters to meters, they drift so fast inwards, that they are quickly lost. Aims: We attempt to verify if collisional fragmentation of dust aggregates is effective enough to keep disks 'dusty' by replenishing the population of small grains and by preventing excessive radial drift. Methods: With a new and sophisticated implicitly integrated coagulation and fragmentation modeling code, we solve the combined problem of coagulation, fragmentation, turbulent mixing and radial drift and at the same time solve for the 1-D viscous gas disk evolution. Results: We find that for a critical collision velocity of 1 m/s, as suggested by laboratory experiments, the fragmentation is so effective, that at all times the dust is in the form of relatively small particles. This means that radial drift is small and that large amounts of small dust particles remain present for a few million years, as observed. For a critical velocity of 10 m/s, we find that particles grow about two orders of magnitude larger, which leads again to significant dust loss since larger particles are more strongly affected by radial drift.

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.

We present a high angular resolution (0.3" = 40 AU) SMA survey of the 870 micron thermal continuum emission from 9 of the brightest, and therefore most massive, circumstellar disks in the ~1 Myr-old Ophiuchus star-forming region. Using 2-D radiative transfer calculations, we simultaneously fit the observed continuum visibilities and broadband spectral energy distribution for each disk with a parametric structure model. Compared to previous millimeter studies, this survey includes significant upgrades in modeling, data quality, and angular resolution that provide improved constraints on key structure parameters, particularly those that characterize the spatial distribution of mass in the disks. In the context of a surface density profile motivated by similarity solutions for viscous accretion disks, the best-fit models for the sample disks have characteristic radii R_c = 20-200 AU, high disk masses M_d = 0.005-0.14 M_sun, and a narrow range of radial surface density gradients around a median $\gamma$ = 0.9. These density structures are used in conjunction with accretion rate estimates from the literature to help characterize the viscous evolution of the disk material. Using the standard prescription for disk viscosities, those combined constraints indicate that $\alpha$ = 0.0005-0.08. Three of the sample disks show large (R = 20-40 AU) central cavities in their continuum emission morphologies, marking extensive zones where dust has been physically removed and/or has significantly diminished opacities. Based on the current requirements of planet formation models, these emission cavities and the structure constraints for the sample as a whole suggest that these young disks may eventually produce planetary systems, and have perhaps already started. (abridged)

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.

In this paper we present two efficient implementations of the diffusion approximation to be employed in Monte Carlo computations of radiative transfer in dusty media of massive circumstellar disks. The aim is to improve the accuracy of the computed temperature structure and to decrease the computation time. The accuracy, efficiency and applicability of the methods in various corners of parameter space are investigated. The effects of using these methods on the vertical structure of the circumstellar disk as obtained from hydrostatic equilibrium computations are also addressed. Two methods are presented. First, an energy diffusion approximation is used to improve the accuracy of the temperature structure in highly obscured regions of the disk, where photon counts are low. Second, a modified random walk approximation is employed to decrease the computation time. This modified random walk ensures that the photons that end up in the high-density regions can quickly escape to the lower density regions, while the energy deposited by these photons in the disk is still computed accurately. A new radiative transfer code, MCMax, is presented in which both these diffusion approximations are implemented. These can be used simultaneously to increase both computational speed and decrease statistical noise. We conclude that the diffusion approximations allow for fast and accurate computations of the temperature structure, vertical disk structure and observables of very optically thick circumstellar disks.

The spectral region around 10 micrometer, showing prominent dust emission bands, is commonly used to derive the chemical composition of protoplanetary dust. Different analysis methods have been proposed for this purpose, but so far, no comparative test has been performed to test their validity. We calculated model spectra of disk models with different geometries and central sources, using a 2D radiative transfer code. These spectra were then fitted in a blind test using four different spectral decomposition methods. We studied the effect of disk structure (flared vs. flat), inclination angle, size of the inner disk hole and stellar luminosity on the fitted chemical composition. Our results show that the derived dust compositions by all methods deviate systematically from the real chemical composition. Out of the four tested spectral decomposition methods, our new two-layer temperature distribution method, differs the least from the input dust composition and the results show the weakest systematic effects. The reason for the deviations of the results given by other methods lies in their simplifying assumptions (e.g. single average grain temperature or one component continuum). We also tested the influence of different noise levels on the results of the spectral decomposition methods. We find that, for ground-based observations (8-13 micrometer) the expected uncertainty in the value of the crystallinity is about 11% for a signal-to-noise ratio of 100, while for space-based observations (7-17 micrometer) the uncertainty is about 5%. On the basis of our results, we propose a recipe for the analysis and interpretation of dust spectroscopy data in the mid-infrared which should be especially valuable for analysing Spitzer data and ground-based infrared spectroscopy data in the 10 micrometer window.

(Abridged) Aims & Methods. A two-dimensional, semi-analytical model is presented that follows, for the first time, the chemical evolution from a collapsing molecular cloud (a pre-stellar core) to a protostar and circumstellar disk. The model computes infall trajectories from any point in the cloud and tracks the radial and vertical motion of material in the viscously evolving disk. It includes a full time-dependent radiative transfer treatment of the dust temperature, which controls much of the chemistry. A small parameter grid is explored to understand the effects of the sound speed and the mass and rotation of the cloud. The freeze-out and evaporation of carbon monoxide (CO) and water (H2O), as well as the potential for forming complex organic molecules in ices, are considered as important first steps to illustrate the full chemistry. Results. Both species freeze out towards the centre before the collapse begins. Pure CO ice evaporates during the infall phase and re-adsorbs in those parts of the disk that cool below the CO desorption temperature of ~18 K. H2O remains solid almost everywhere during the infall and disk formation phases and evaporates within ~10 AU of the star. Mixed CO-H2O ices are important in keeping some solid CO above 18 K and in explaining the presence of CO in comets. Material that ends up in the planet- and comet-forming zones of the disk is predicted to spend enough time in a warm zone during the collapse to form first-generation complex organic species on the grains. The dynamical timescales in the hot inner envelope (hot core or hot corino) are too short for abundant formation of second-generation molecules by high-temperature gas-phase chemistry.

PAHs have been detected toward molecular clouds and some young stars with disks, but have not yet been associated with embedded young stars. We present a sensitive mid-IR spectroscopic survey of PAH features toward a sample of low-mass embedded YSOs. The aim is to put constraints on the PAH abundance in the embedded phase of star formation using radiative transfer modeling. VLT-ISAAC L-band spectra for 39 sources and Spitzer IRS spectra for 53 sources are presented. Line intensities are compared to recent surveys of Herbig Ae/Be and T Tauri stars. The radiative transfer codes RADMC and RADICAL are used to model the PAH emission from embedded YSOs consisting of a PMS star with a circumstellar disk embedded in an envelope. The dependence of the PAH feature on PAH abundance, stellar radiation field, inclination and the extinction by the surrounding envelope is studied. The 3.3 micron PAH feature is undetected for the majority of the sample (97%), with typical upper limits of 5E-16 W/m^2. Compact 11.2 micron PAH emission is seen directly towards 1 out of the 53 Spitzer Short-High spectra, for a source that is borderline embedded. For all 12 sources with both VLT and Spitzer spectra, no PAH features are detected in either. In total, PAH features are detected toward at most 1 out of 63 (candidate) embedded protostars (<~ 2%), even lower than observed for class II T Tauri stars with disks (11-14%). Assuming typical class I stellar and envelope parameters, the absence of PAHs emission is most likely explained by the absence of emitting carriers through a PAH abundance at least an order of magnitude lower than in molecular clouds but similar to that found in disks. Thus, most PAHs likely enter the protoplanetary disks frozen out in icy layers on dust grains and/or in coagulated form.

We present a parameter study of self-consistent models of protoplanetary disks around Herbig AeBe stars. We use the code developed by Dullemond and Dominik, which solves the 2D radiative transfer problem including an iteration for the vertical hydrostatic structure of the disk. This grid of models will be used for several studies on disk emission and mineralogy in followup papers. In this paper we take a first look on the new models, compare them with previous modeling attempts and focus on the effects of various parameters on the overall structure of the SED that leads to the classification of Herbig AeBe stars into two groups, with a flaring (group I) or self-shadowed (group II) SED. We find that the parameter of overriding importance to the SED is the total mass in grains smaller than 25um, confirming the earlier results by Dullemond and Dominik. All other parameters studied have only minor influences, and will alter the SED type only in borderline cases. We find that there is no natural dichotomy between group I and II. From a modeling point of view, the transition is a continuous function of the small dust mass. We also show that moderate grain growth produces spectra with weak or no 1um feature, both for flaring (Group I) and non-flaring (Group II) sources. The fact that sources with weak features have been found mostly in Group I sources is therefore surprising and must be due to observational biases or evolutionary effects.

Turbulent coagulation in protoplanetary disks is known to operate on timescale far shorter than the lifetime of the disk. In the absence of mechanisms that replenish the small dust grain population, protoplanetary disks would rapidly lose their continuum opacity-bearing dust. This is inconsistent with infrared observations of disks around T Tauri stars and Herbig Ae/Be stars, which are usually optically thick at visual wavelengths and show signatures of small (a<~ 3um) grains. A plausible replenishing mechanism of small grains is collisional fragmentation or erosion of large dust aggregates, which model calculations predict to play an important role in protoplanetary disks. If optically thick disks are to be seen as proof for ongoing fragmentation or erosion, then alternative explanations for the existence of optically thick disks must be studied carefully. In this study we explore two scenarios. First, we study the effect of residual, low-level infall of matter onto the disk surface. We find that infall rates as low as 10^{-11} Msun/yr can, in principle, replenish the small grain population to a level that keeps the disk marginally optically thick. However, it remains to be seen if the assumption of such inflow is realistic for star+disk systems at the age of several Myrs, at which winds and jets are expected to have removed any residual envelope. In summary, fragmentation or erosion still appear to be the most promising processes to explain the abundant presence of small grains in old disks.

Context: There is increasing need for good algorithms for modeling the aggregation and fragmentation of solid particles (dust grains, dust aggregates, boulders) in various astrophysical settings, including protoplanetary disks, planetary- and sub-stellar atmospheres and dense molecular cloud cores. Here we describe a new algorithm that combines advantages of various standard methods into one. Aims: The aim is to develop a method that 1) can solve for aggregation and fragmentation, 2) can easily include the effect and evolution of grain properties such as compactness, composition, etc., and 3) can be built as a coagulation/fragmentation module into a hydrodynamics simulations. Methods: We develop a Monte-Carlo method in which we follow the 'life' of a limited number of representative particles. Each of these particles is associated with a certain fraction of the total dust mass and thereby represents a large number of true particles which all are assumed to have the same properties as their representative particle. Under the assumption that the total number of true particles vastly exceeds the number of representative particles, the chance of a representative particle colliding with another representative particle is negligibly small, and we therefore ignore this possibility. This now makes it possible to employ a statistical approach to the evolution of the representative particles. Results: The method reproduces the known analytic solutions of simplified coagulation kernels, and compares well to numerical results for Brownian motion using other methods. For reasonably well-behaved kernels it produces good results even for moderate number of swarms.

Recent observations of two black hole candidates (GX 339-4 and J1753.5-0127) in the low-hard state (L_X/L_Edd ~ 0.003-0.05) suggest the presence of a cool accretion disk very close to the innermost stable orbit of the black hole. This runs counter to models of the low-hard state in which the cool disk is truncated at a much larger radius. We study the interaction between a moderately truncated disk and a hot inner flow. Ion-bombardment heats the surface of the disk in the overlap region between a two-temperature advection-dominated accretion flow and standard accretion disk, producing a hot (kT_e ~70 keV) layer on the surface of the cool disk. The hard X-ray flux from this layer heats the inner parts of the underlying cool disk, producing a soft X-ray excess. Together with interstellar absorption these effects mimic the thermal spectrum from a disk extending to the last stable orbit. The results show that soft excesses in the low-hard state are a natural feature of truncated disk models.

(Abridged) Classification schemes for YSOs are based on evaluating the degree of dissipation of the surrounding envelope, whose main effects are the extinction of the optical radiation from the central YSO and re-emission in the far-infrared. Since extinction is a property of column density along the line of sight, the presence of a protoplanetary disk may lead to a misclassification when the system is viewed edge-on. We performed radiative transfer calculations, using the axysimmetric 3D radiative transfer codes RADMC and RADICAL, to show the effects of different geometries on the main indicators of YSO evolutionary stage, like the slope of the flux between 2 and 24mum, the bolometric temperature and the optical depth of silicates and ices. We show that for systems viewed at intermediate angles the 'classical' indicators of evolution accurately trace the envelope column density, and they all agree with each other. On the other hand, edge-on system are misclassified for inclinations larger than ~65deg. In particular, silicate emission, typical of pre-main sequence stars with disks, turns into absorption when the disk column density reaches 1e22cm-2, corresponding e.g. to a 5e-3 Msun flaring disk viewed at 64deg. A similar effect is noticed in all the other classification indicators studied alpha, Tbol, and the H2O and CO2 ices absorption strengths. This misclassification has a big impact on the nature of the flat-spectrum sources (alpha ~0), whose number can be explained by simple geometrical arguments without invoking evolution. A reliable classification scheme using a minimal number of observations is constituted by observations of the mm-flux with both a single dish and an interferometer.

The formation of planetesimals in protoplanetary disks due to collisional sticking of smaller dust aggregates has to face at least two severe obstacles, namely the rapid loss of material due to radial inward drift and particle fragmentation due to destructive collisions. We present a scenario to circumvent these two hurdles. Our dust evolution model involves two main mechanisms. First, we consider a disk with a dead zone. In an almost laminar region close to the midplane, the relative velocities of the turbulent particles are comparatively small, which decreases the probability of destructive particle collisions. Second, turbulence is not the only source of violent relative particle velocities, because high radial drift speeds can also lead to boulder fragmentation. For this reason, we focus additionally on the snow line. Evaporation fronts can be associated with gas pressure maxima in which radial drift basically vanishes. This implies that particle fragmentation becomes even less likely. Our simulation results suggest that particles can overcome the fragmentation barrier. We find that boulders of several 100 m can form within only a few thousand years.

Aims: We wish to investigate what the effect of dust sedimentation is on the observed 10 mum feature of protoplanetary disks and how this may affect the interpretation of the observations. Methods: Using a combination of modeling tools, we simulated the sedimentation of a dust grain size distribution in an axisymmetric 2-D model of a turbulent protoplanetary disk, and we used a radiative transfer program to compute the resulting spectra. Results: We find that the sedimentation can turn a flat feature into a pointy one, but only to a limited degree and for a very limited set of particle size distributions. Only if we have a bimodal size distribution, i.e. a very small grain population and a bigger grain population, do we find that the transformation from a flat to a pointy feature upon dust sedimentation is strong. However, our model shows that, if sedimentation is the sole reason for the variety of silicate feature strengths observed in protoplanetary disks, then we would expect to find a correlation such that disks with weak mid- to far-infrared excess have a stronger 10 mum silicate feature than disks with a strong mid- to far-infrared excess. If this is contrary to what is observed, then this would indicate that sedimentation cannot be the main reason for the variety of 10 mum silicate features observed in protoplanetary disks.

We observed molecular hydrogen around a sample of pre-main sequence stars in order to better characterize their gaseous CS environments. We analyzed the FUSE (Far Ultraviolet Spectroscopic Explorer) spectra of a sample of Herbig Ae/Be stars (HAeBes) covering a broad spectral range, including the main-sequence A5 star Beta-Pictoris. To better diagnose the origin of the detected H2 and its excitation conditions, we used a model of a photodissociation region. Our analysis demonstrates that the excitation of H2 is clearly different around most of the HAeBes compared to the interstellar medium. Moreover, the characteristics of H2 around Herbig Ae and Be stars give evidence for different excitation mechanisms. For the most massive stars of our sample (B8 to B2 type), the excitation diagrams are reproduced well by a model of photodissociation regions (PDR). Our results favor an interpretation in terms of large CS envelopes, remnants of the molecular clouds in which the stars were formed. On the other hand, the group of Ae stars (later than B9 type) known to possess disks is more inhomogeneous. In most cases, when CS H2 is detected, the lines of sight do not pass through the disks. The excitation conditions of H2 around Ae stars cannot be reproduced by PDR models and correspond to warm and/or hot excited media very close to the stars. In addition, no clear correlation has been found between the ages of the stars and the amount of circumstellar H2. Our results suggest structural differences between Herbig Ae and Be star environments. Herbig Be stars do evolve faster than Ae stars, and consequently, most Herbig Be stars are younger than Ae ones at the time we observe them. It is thus more likely to find remnants of their parent cloud around them.

To treat the problem of growing protoplanetary disc solids across the meter barrier, we consider a very simplified two-component coagulation-fragmentation model that consists of macroscopic boulders and smaller dust grains, the latter being the result of catastrophic collisions between the boulders. Boulders in turn increase their radii by sweeping up the dust fragments. An analytical solution to the dynamical equations predicts that growth by coagulation-fragmentation can be efficient and allow agglomeration of 10-meter-sized objects within the time-scale of the radial drift. These results are supported by computer simulations of the motion of boulders and fragments in 3-D time-dependent magnetorotational turbulence. Allowing however the fragments to diffuse freely out of the sedimentary layer of boulders reduces the density of both boulders and fragments in the mid-plane, and thus also the growth of the boulder radius, drastically. The reason is that the turbulent diffusion time-scale is so much shorter than the collisional time-scale that dust fragments leak out of the mid-plane layer before they can be swept up by the boulders there. Our conclusion that coagulation-fragmentation is not an efficient way to grow across the meter barrier in fully turbulent protoplanetary discs confirms recent results by Brauer, Dullemond, & Henning who solved the coagulation equation in a parameterised turbulence model with collisional fragmentation, cratering, radial drift, and a range of particle sizes. We find that a relatively small population of boulders in a sedimentary mid-plane layer can populate the entire vertical extent of the disc with small grains and that these grains are not first generation dust, but have been through several agglomeration-destruction cycles during the simulations.

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. We present here the direct characterization of a 40 AU radius inner gap in the disk around LkHa 330 through 340 GHz (880 micron) dust continuum imaging with the Submillimeter Array (SMA). This large gap is fully resolved by the SMA observations and mostly empty of dust with less than 1.3 x 10^-6 M_solar of solid particles inside of 40 AU. Gas (as traced by accretion markers and CO M-band emission) is still present in the inner disk and the outer edge of the gap rises steeply -- features in better agreement with the underlying cause being gravitational perturbation than a more gradual process such as grain growth. Importantly, the good agreement of the spatially resolved data and spectrophometry-based model lends confidence to current interpretations of SEDs with significant dust emission deficits as arising 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.

Recent work has suggested that the fraction of obscured AGN declines with increasing luminosity, but it has been difficult to quantify this trend. Here, we attempt to measure this fraction as a function of luminosity by studying the ratio of mid-infrared to intrinsic nuclear bolometric luminosity in unobscured AGN. Because the mid-infrared is created by dust reprocessing of shorter wavelength nuclear light, this ratio is a diagnostic of f_obsc, the fraction of solid angle around the nucleus covered by obscuring matter. In order to eliminate possible redshift-dependences while also achieving a large dynamic range in luminosity, we have collected archival 24 micron MIPS photometry from objects with z~1 in the Sloan Digital Sky Survey (SDSS), the Great Observatories Origins Deep Survey (GOODS) and the Cosmic Evolution Survey (COSMOS). To measure the bolometric luminosity for each object, we used archival optical data supplemented by GALEX data. We find that the mean ratio of 24 microns to bolometric luminosity decreases by a factor of ~3 in the L_bol=10^44-3x10^47 ergs s^-1 range, but there is also a large scatter at constant L_bol. Using radiation transfer solutions for model geometries, we show how the IR/bolometric ratio relates to f_obsc and compare these values with those obtained obtained from samples of X-ray selected AGN. Although we find approximate agreement, our method indicates somewhat higher values of f_obsc, particularly in the middle range of luminosities, suggesting that there may be a significant number of heavily obscured AGN missed by X-ray surveys.

We observed the Herbig Ae/Be stars UX Ori, HD 34282, HD 100453, HD 101412, HD 104237 and HD 142666, and the T Tauri star HD 319139 and searched for H2 0-0 S(2) emission at 12.278 micron and H2 0-0 S(1) emission at 17.035 micron with VISIR, ESO-VLT's high-resolution MIR spectrograph. None of the sources present evidence for H2 emission. Stringent 3sigma upper limits to the integrated line fluxes and the mass of optically thin warm gas in the disks are derived. The disks contain less than a few tenths of Jupiter mass of optically thin H2 gas at 150 K at most, and less than a few Earth masses of optically thin H2 gas at 300 K and higher temperatures. We compare our results to a Chiang and Goldreich (1997, CG97) two-layer disk model. The upper limits to the disk's optically thin warm gas mass are smaller than the amount of warm gas in the interior layer of the disk, but they are much larger than the amount of molecular gas in the surface layer. We present a calculation of the expected thermal H2 emission from optically thick disks, assuming a CG97 disk structure, a gas-to-dust ratio of 100 and Tgas = Tdust. The expected H2 thermal emission fluxes from typical disks around Herbig Ae/Be stars (10^-16 to 10^-17 erg/s/cm2 at 140 pc) are much lower than the detection limits of our observations (5*10^-15 erg/s/cm2). H2 emission levels are very sensitive to departures from the thermal coupling between the molecular gas and dust. Additional sources of heating of gas in the disk's surface layer could have a major impact on the expected H2 disk emission. In the observed sources the molecular gas and dust in the surface layer have not significantly departed from thermal coupling (Tgas/Tdust< 2) and that the gas-to-dust ratio in the surface layer is very likely lower than 1000.

The growth of solid particles towards meter sizes in protoplanetary disks has to circumvent at least two hurdles, namely the rapid loss of material due to radial drift and particle fragmentation due to destructive collisions. In this paper, we present the results of numerical simulations with more and more realistic physics involved. Step by step, we include various effects, such as particle growth, radial/vertical particle motion and dust particle fragmentation in our simulations. We demonstrate that the initial dust-to-gas ratio is essential for the particles to overcome the radial drift barrier. If this value is increased by a factor of 2 compared with the canonical value for the interstellar medium, km-sized bodies can form in the inner disk <2 AU within 10 thousand years. However, we find that solid particles get destroyed through collisional fragmentation. Only with the unrealistically high-threshold velocities needed for fragmentation to occur (>30 m/s), particles are able to grow to larger sizes in low turbulent disks. We also find that less than 5% of the small dust grains remain in the disk after 1 Myrs due to radial drift, no matter whether fragmentation is included in the simulations or not. In this paper, we also present considerable improvements to existing algorithms for dust-particle coagulation, which speed up the coagulation scheme by a factor of 10 thousand.

We want to improve knowledge of the structure of the inner few AU of the circumstellar disk around the nearby T Tauri star TW Hya. Earlier studies have suggested the existence of a large inner hole, possibly caused by interactions with a growing protoplanet. We used interferometric observations in the N-band obtained with the MIDI instrument on the Very Large Telescope Interferometer, together with 10 micron spectra recorded by the infrared satellite Spitzer. The fact that we were able to determine N-band correlated fluxes and visibilities for this comparatively faint source shows that MIR interferometry can be applied to a large number of low-mass young stellar objects. The MIR spectra obtained with Spitzer reveal emission lines from HI (6-5), HI (7-6), and [Ne II] and show that over 90% of the dust we see in this wavelength regime is amorphous. According to the correlated flux measured with MIDI, most of the crystalline material is in the inner, unresolved part of the disk, about 1 AU in radius. The visibilities exclude the existence of a very large (3-4 AU radius) inner hole in the circumstellar disk of TW Hya, which was required in earlier models. We propose instead a geometry of the inner disk where an inner hole still exists, but at a much reduced radius, with the transition from zero to full disk height between 0.5 and 0.8 AU, and with an optically thin distribution of dust inside. Such a model can comply with SED and MIR visibilities, as well as with visibility and extended emission observed in the NIR at 2 micron. If a massive planet was the reason for this inner hole, as has been speculated, its orbit would have to be closer to the star than 0.3 AU. Alternatively, we may be witnessing the end of the accretion phase and an early phase of an inward-out dispersal of the circumstellar disk.

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.

The process of turbulent radial mixing in protoplanetary disks has strong relevance to the analysis of the spatial distribution of crystalline dust species in disks around young stars and to studies of the composition of meteorites and comets in our own solar system. A debate has gone on in the recent literature on the ratio of the effective viscosity coefficient $\nu$ (responsible for accretion) to the turbulent diffusion coefficient $D$ (responsible for mixing). Numerical magneto-hydrodynamic simulations have yielded values between $\nu/D\simeq 10$ (Carballido, Stone & Pringle, 2005) and $\nu/D\simeq 0.85$ (Johansen & Klahr, 2005}). Here we present two analytic arguments for the ratio $\nu/D=1/3$ which are based on elegant, though strongly simplified assumptions. We argue that whichever of these numbers comes closest to reality may be determined {\em observationally} by using spatially resolved mid-infrared measurements of protoplanetary disks around Herbig stars. If meridional flows are present in the disk, then we expect less abundance of crystalline dust in the surface layers, a prediction which can likewise be observationally tested with mid-infrared interferometers.

We present spatially resolved mid-infrared images of the disk surrounding the young star IRS 48 in the Ophiuchus cloud complex. The disk exhibits a ring-like structure at 18.7 micron, and is dominated by very strong emission from polycyclic aromatic hydrocarbons at shorter wavelengths. This allows a detailed study of the relative distributions of small and large dust grains. Images of IRS 48 in 5 mid-infrared bands from 8.6 to 18.7 micron as well as a low resolution N-band spectrum are obtained with VLT-VISIR. Optical spectroscopy is used to determine the spectral type of the central star and to measure the strength of the Halpha line. The 18.7 micron ring peaks at a diameter of 110 AU, with a gap of ~ 60 AU. The shape of the ring is consistent with an inclination of i = 48 +- 8 degrees. In contrast, the 7.5-13 micron PAH emission bands are centered on the source and appear to fill the gap within the ring. The measured PAH line strengths are 10-100x stronger than those typically measured for young M0 stars and can only be explained with a high PAH abundance and/or strong excess optical/UV emission. The morphology of the images, combined with the absence of a silicate emission feature, imply that the inner disk has been cleared of micron-sized dust but with a significant population of PAHs remaining. We argue that the gap can be due to grain growth and settling or to clearing by an unseen planetary or low-mass companion. IRS 48 may represent a short-lived transitional phase from a classical to a weak-line T Tauri star.

Millimeter interferometry provides evidence for the presence of mm to cm size "pebbles" in the outer parts of disks around pre-main-sequence stars. The observations suggest that large grains are produced relatively early in disk evolution (< 1 Myr) and remain at large radii for longer periods of time (5 to 10 Myr). Simple theoretical estimates of the radial drift time of solid particles, however, imply that they would drift inward over a time scale of less than 0.1 Myr. In this paper, we address this conflict between theory and observation, using more detailed theoretical models, including the effects of sedimentation, collective drag forces and turbulent viscosity. We find that, although these effects slow down the radial drift of the dust particles, this reduction is not sufficient to explain the observationally determined long survival time of mm/cm-sized grains in protoplanetary disks. However, if for some reason the gas to dust ratio in the disk is reduced by at least a factor of 20 from the canonical value of 100 (for instance through photoevaporation of the gas), then the radial drift time scales become sufficiently large to be in agreement with observations.

We present deep Spitzer-IRS low-resolution (lambda/Delta lambda ~ 100) 5-35 micron spectroscopy of the edge-on disk ``the Flying Saucer'' (2MASS J16281370-2431391) in the Ophiuchus molecular cloud. The spectral energy distribution exhibits the characteristic two-peak shape predicted for a circumstellar disk viewed very close to edge-on. The short-wavelength peak is entirely due to photons scattered off the surface of the disk, while the long-wavelength peak is due to thermal emission from the disk itself. The Spitzer spectrum represents the first spectroscopic detection of scattered light out to 15 micron from a bona-fide, isolated edge-on disk around a T Tauri star. The depth and the wavelength of the mid-infrared "valley" of the SED give direct constraints on the size distribution of large grains in the disk. Using a 2D continuum radiative transfer model, we find that a significant amount of 5-10 micron-sized grains is required in the surface layers of the disk at radii of 50-300 AU. The detection of relatively large grains in the upper layers implies that vertical mixing is effective, since grain growth models predict the grains would otherwise settle deep in the disk on short time scales. Additionally, we tentatively detect the 9.66 micron S(3) line of H2 and the 11.2 micron emission feature due to PAHs.

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.

Aims. The chemistry of, and infrared (IR) emission from, polycyclic aromatic hydrocarbons (PAHs) in disks around Herbig Ae/Be and T Tauri stars are investigated. The equilibrium distribution of the PAHs over all accessible charge/hydrogenation states depends on the size and shape of the PAHs and on the physical properties of the star and surrounding disk. Methods. A chemistry model is created to calculate this equilibrium distribution. Destruction of PAHs by ultraviolet (UV) photons, possibly in multi-photon absorption events, is taken into account. The chemistry model is coupled to a radiative transfer code to provide the physical parameters and to combine the PAH emission with the spectral energy distribution (SED) from the star+disk system. Results. Normally hydrogenated PAHs in Herbig Ae/Be disks account for most of the observed PAH emission, with neutral and positively ionized species contributing in roughly equal amounts. Close to the midplane, the PAHs are more strongly hydrogenated and negatively ionized, but these species do not contribute to the overall emission because of the low UV/optical flux deep inside the disk. PAHs of 50 carbon atoms are destroyed out to 100 AU in the disk's surface layer, and the resulting spatial extent of the emission does not agree well with observations. Rather, PAHs of about 100 carbon atoms or more are predicted to cause most of the observed emission. The emission is extended on a scale similar to that of the size of the disk. Furthermore, the emission from T Tauri disks is much weaker and concentrated more towards the central star than that from Herbig Ae/Be disks. Positively ionized PAHs are predicted to be largely absent in T Tauri disks because of the weaker radiation field.

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.

We investigate the long-term optical-infrared variability of SV Cep, and explain it in the context of an existing UX Ori (UXOR) model. A 25-month monitoring programme was completed with the Infrared Space Observatory in the 3.3-100 um wavelength range. Following a careful data reduction, the infrared light curves were correlated with the variations of SV Cep in the V-band. A remarkable correlation was found between the optical and the far-infrared light curves. In the mid-infrared regime the amplitude of variations is lower, with a hint for a weak anti-correlation with the optical changes. In order to interpret the observations, we modelled the spectral energy distribution of SV Cep assuming a self-shadowed disc with a puffed-up inner rim, using a 2-dimensional radiative transfer code. We found that modifying the height of the inner rim, the wavelength-dependence of the long-term optical-infrared variations is well reproduced, except the mid-infrared domain. The origin of variation of the rim height might be fluctuation in the accretion rate in the outer disc. In order to model the mid-infrared behaviour we tested to add an optically thin envelope to the system, but this model failed to explain the far-infrared variability. Infrared variability is a powerful tool to discriminate between models of the circumstellar environment. The proposed mechanism of variable rim height may not be restricted to UXOR stars, but might be a general characteristic of intermediate-mass young stars.

Aims: To calculate chemistry and gas temperature of evolving protoplanetary disks with decreasing mass or dust settling, and to explore the sensitivity of gas-phase tracers. Methods: The density and dust temperature profiles for a range of models of flaring and self-shadowed disks around a typical Herbig Ae star are used together with 2-dimensional ultraviolet (UV) radiative transfer to calculate the chemistry and gas temperature. In each model the line profiles and intensities for the fine structure lines of [O I], [C II] and [C I] and the pure rotational lines of CO, CN, HCN and HCO+ are determined. Results: The chemistry shows a strong correlation with disk mass. Molecules that are easily dissociated, like HCN, require high densities and large extinctions before they can become abundant. The products of photodissociation, like CN and C2H, become abundant in models with lower masses. Dust settling mainly affects the gas temperature, and thus high temperature tracers like the O and C+ fine structure lines. The carbon chemistry is found to be very sensitive to the adopted PAH abundance. The line ratios CO/13CO, CO/HCO+ and [O I] 63 um/146 um can be used to distinguish between disks where dust growth and settling takes place, and disks that undergo overall mass loss.

We present mid-infrared Spitzer IRAC and MIPS images of the UX Orionis star VV Ser and the surrounding cloud. The 5.6--70 micron images show bright, localized and nebulous emission extended over 4 arcmin centered on VV Ser. We interpret the nebulosity as being due to transiently heated grains excited by UV photons emitted by VV Ser. A companion paper describes how the physical structure of the VV Ser disk has been constrained using a wide range of observational data modeled by an axisymmetric Monte Carlo radiative transfer code. In this paper we employ the model to study the nebulosity surrounding VV Ser using quantum-heated PAH molecules and Very Small Grains (VSGs) consisting of amorphous carbon in the thermal cooling approximation. Imprinted on the nebulosity is a wedge-shaped dark band, centered on the star. We interpret this dark wedge as the shadow cast by the inner regions of a near-edge-on disk in UV light, allowing the PAHs to be excited only outside of this shadow. The presence of a disk shadow strongly constrains the inclination as well as the position angle of the disk. Although depending on the adopted PAH opacity, the abundance of PAHs in the surrounding cloud is constrained to 5+/-2% of the total dust mass, given the opacity. The extent of the nebulosity constrains the density of the gas surrounding the VV Ser disk to 500+/-200 cm-3 for a gas-to-dust ratio of 100. This low density suggests that the quantum heated material is not part of the original envelope of VV Ser and that it is rather a quiescent part of the Serpens molecular cloud that the star has passed into after being formed. [Abstract abridged]

We present mid-infrared Spitzer-IRS spectra of the well-known UX Orionis star VV Ser. We combine the Spitzer data with interferometric and spectroscopic data from the literature covering UV to submillimeter wavelengths. The full set of data are modeled by a two-dimensional axisymmetric Monte Carlo radiative transfer code. The model is used to test the prediction of (Dullemond et al. 2003) that disks around UX Orionis stars must have a self-shadowed shape, and that these disks are seen nearly edge-on, looking just over the edge of a puffed-up inner rim, formed roughly at the dust sublimation radius. We find that a single, relatively simple model is consistent with all the available observational constraints spanning 4 orders of magnitude in wavelength and spatial scales, providing strong support for this interpretation of UX Orionis stars. The grains in the upper layers of the puffed-up inner rim must be small (0.01-0.4 micron) to reproduce the colors (R_V ~ 3.6) of the extinction events, while the shape and strength of the mid-infrared silicate emission features indicate that grains in the outer disk (> 1-2 AU) are somewhat larger (0.3-3.0 micron). From the model fit, the location of the puffed-up inner rim is estimated to be at a dust temperature of 1500 K or at 0.7-0.8 AU for small grains. This is almost twice the rim radius estimated from near-infrared interferometry. A best fitting model for the inner rim in which large grains in the disk mid-plane reach to within 0.25 AU of the star, while small grains in the disk surface create a puffed-up inner rim at ~0.7-0.8 AU, is able to reproduce all the data, including the near-infrared visibilities. [Abstract abridged]

Aims: We search for PAH features towards T Tauri stars and compare them with surveys of Herbig Ae/Be stars. The presence and strength of the PAH features are interpreted with disk radiative transfer models exploring the PAH feature dependence on the incident UV radiation, PAH abundance and disk parameters. Methods: Spitzer Space Telescope 5-35 micron spectra of 54 pre-main sequence stars with disks were obtained, consisting of 38 T Tauri, 7 Herbig Ae/Be and 9 stars with unknown spectral type. Results: Compact PAH emission is detected towards at least 8 sources of which 5 are Herbig Ae/Be stars. The 11.2 micron PAH feature is detected in all of these sources, as is the 6.2 micron PAH feature where short wavelength data are available. However, the 7.7 and 8.6 micron features appear strongly in only 1 of these 4 sources. PAH emission is observed towards at least 3 T Tauri stars (8% detection rate). The lowest mass source with PAHs in our sample is T Cha (G8). All 4 sources in our sample with evidence for dust holes in their inner disk show PAH emission, increasing the feature/continuum ratio. Typical 11.2 micron line intensities are an order of magnitude lower than those observed for the more massive Herbig Ae/Be stars. Measured line fluxes indicate PAH abundances that are factors of 10-100 lower than standard interstellar values. Conversely, PAH features from disks exposed to stars with Teff<=4200K without enhanced UV are predicted to be below the current detection limit, even for high PAH abundances. Disk modeling shows that the 6.2 and 11.2 micron features are the best PAH tracers for T Tauri stars, whereas the 7.7 and 8.6 micron bands have low feature over continuum ratios due to the strongly rising silicate emission.

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 $\mu$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.

We report the first detections of the Class 0 protostellar source IRAM 04191+1522 at wavelengths shortward of 60 microns with the Spitzer Space Telescope. We see extended emission in the Spitzer images that suggests the presence of an outflow cavity in the circumstellar envelope. We combine the Spitzer observations with existing data to form a complete dataset ranging from 3.6 to 1300 microns and use these data to construct radiative transfer models of the source. We conclude that the internal luminosity of IRAM 04191+1522, defined to be the sum of the luminosity from the internal sources (a star and a disk), is L_int = 0.08 +/- 0.04 L_sun, placing it among the lowest luminosity protostars known. Though it was discovered before the launch of the Spitzer Space Telescope, IRAM 04191+1522 falls within a new class of Very Low Luminosity Objects being discovered by Spitzer. Unlike the two other well-studied objects in this class, which are associated either with weak, compact outflows or no outflows at all, IRAM 04191+1522 has a well-defined molecular outflow with properties consistent with those expected based on relations derived from higher luminosity (L_int > 1 L_sun) protostars. We discuss the difficulties in understanding IRAM 04191+1522 in the context of the standard model of star formation, and suggest a possible explanation for the very low luminosity of this source.

Using the Hectochelle multifiber spectrograph, we have obtained high-resolution (R~34,000) spectra in the Halpha region for a large number of stars in the 4 Myr-old cluster Tr 37, containing 146 previously known members and 26 newly identified ones. We present the Halpha line profiles of all members, compare them to our IR observations of dusty disks (2MASS/JHK + IRAC + MIPS 24 micron), use the radial velocities as a membership criterion, and calculate the rotational velocities. We find a good correlation between the accretion-broadened profiles and the presence of protoplanetary disks, noting that a small fraction of the accreting stars presents broad profiles with Halpha equivalent widths smaller than the canonical limit separating CTTS and WTTS. The number of strong accretors appears to be lower than in younger regions, and a large number of CTTS have very small accretion rates (dM/dt<10^{-9} Msun/yr). Taking into account that the spectral energy distributions are consistent with dust evolution (grain growth/settling) in the innermost disk, this suggests a parallel evolution of the dusty and gaseous components. We also observe that about half of the "transition objects" (stars with no IR excesses at wavelengths shorter than ~6 micron) do not show any signs of active accretion, whereas the other half is accreting with accretion rates <10^{-9} Msun/yr. These zero or very low accretion rates reveal important gas evolution and/or gas depletion in the innermost disk, which could be related to grain growth up to planetesimal or even planet sizes. Finally, we examine the rotational velocities of accreting and non accreting stars, finding no significant differences that could indicate disk locking at these ages.

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.

In this Letter we present a theoretical scenario to explain the steep correlation between disk accretion rates and stellar masses observed in pre-main sequence stars. We show that the correlations and spread observed in the two best studied regions, rho-Ophiuchus and Taurus, can be reproduced by a simple model of single star formation from a rotating collapsing core and the viscous evolution of the circumstellar disk. In this model, the rate of rotation of the parent core sets the value of the `centrifugal radius' within which the infalling matter is loaded onto the surface of the disk. As a consequence, the disk accretion rate measured long after the dispersal of the parental core bears the imprint of the initial conditions of star formation. The observed trend results naturally if, at the onset of the collapse, cores of all masses rotate with the same distribution of angular velocities, measured in units of the break-up rotation rate.

We review advances in the modeling of protoplanetary disks. This review will focus on the regions of the disk beyond the dust sublimation radius, i.e. beyond 0.1 - 1 AU, depending on the stellar luminosity. We will be mostly concerned with models that aim to fit spectra of the dust continuum or gas lines, and derive physical parameters from these fits. For optically thick disks, these parameters include the accretion rate through the disk onto the star, the geometry of the disk, the dust properties, the surface chemistry and the thermal balance of the gas. For the latter we are mostly concerned with the upper layers of the disk, where the gas and dust temperature decouple and a photoevaporative flow may originate. We also briefly discuss optically thin disks, focusing mainly on the gas, not the dust. The evolution of these disks is dominated by accretion, viscous spreading, photoevaporation, and dust settling and coagulation. The density and temperature structure arising from the surface layer models provide input to models of photoevaporation, which occurs largely in the outer disk. We discuss the consequences of photoevaporation on disk evolution and planet formation.

We present a new perspective on the crystallinity of dust in protoplanetary disks. The dominant crystallization by thermal annealing happens in the very early phases of disk formation and evolution. Both the disk properties and the level of crystallinity are thereby directly linked to the properties of the molecular cloud core from which the star+disk system was formed. We show that, under the assumption of single star formation, rapidly rotating clouds produce disks which, after the main infall phase (i.e. in the optically revealed class II phase), are rather massive and have a high accretion rate but low crystallinity. Slowly rotating clouds, on the other hand, produce less massive disks with lower accretion rate, but high levels of crystallinity. Cloud fragmentation and the formation of multiple stars complicates the problem and necessitates further study. The underlying physics of the model is insufficiently understood to provide the precise relationship between crystallinity, disk mass and accretion rate. But the fact that with `standard' input physics the model produces disks which, in comparison to observations, appear to have either too high levels of crystallinity or too high disk masses, demonstrates that the comparison of these models to observations can place strong contraints on the disk physics. The question to ask is not why some sources are so crystalline, but why some other sources have such a low level of crystallinity.

Spitzer Space Telescope spectra of the low mass young stellar object (YSO) IRS 46 (L_bol ~ 0.6 L_sun) in Ophiuchus reveal strong vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2. This is the only source out of a sample of ~100 YSO's that shows these features and the first time they are seen in the spectrum of a solar-mass YSO. Analysis of the Spitzer data combined with Keck L- and M-band spectra gives excitation temperatures of > 350 K and abundances of 10(-6)-10(-5) with respect to H2, orders of magnitude higher than those found in cold clouds. In spite of this high abundance, the HCN J=4-3 line is barely detected with the James Clerk Maxwell Telescope, indicating a source diameter less than 13 AU. The (sub)millimeter continuum emission and the absence of scattered light in near-infrared images limits the mass and temperature of any remnant collapse envelope to less than 0.01 M_sun and 100 K, respectively. This excludes a hot-core type region as found in high-mass YSO's. The most plausible origin of this hot gas rich in organic molecules is in the inner (<6 AU radius) region of the disk around IRS 46, either the disk itself or a disk wind. A nearly edge-on 2-D disk model fits the spectral energy distribution (SED) and gives a column of dense warm gas along the line of sight that is consistent with the absorption data. These data illustrate the unique potential of high-resolution infrared spectroscopy to probe organic chemistry, gas temperatures and kinematics in the planet-forming zones close to a young star.

The onset of planet formation in protoplanetary disks is marked by the growth and crystallization of sub-micron-sized dust grains accompanied by dust settling toward the disk mid-plane. Here we present infrared spectra of disks around brown dwarfs and brown dwarf candidates. We show that all three processes occur in such cool disks in a way similar or identical to that in disks around low- and intermediate-mass stars. These results indicate that the onset of planet formation extends to disks around brown dwarfs, suggesting that planet formation is a robust process occurring in most young circumstellar disks.

Infrared ~5--35 um spectra for 40 solar-mass T Tauri stars and 7 intermediate-mass Herbig Ae stars with circumstellar disks were obtained using the Spitzer Space Telescope as part of the c2d IRS survey. This work complements prior spectroscopic studies of silicate infrared emission from disks, which were focused on intermediate-mass stars, with observations of solar-mass stars limited primarily to the 10 um region. The observed 10 and 20 um silicate feature strengths/shapes are consistent with source-to-source variations in grain size. A large fraction of the features are weak and flat, consistent with um-sized grains indicating fast grain growth (from 0.1--1.0 um in radius). In addition, approximately half of the T Tauri star spectra show crystalline silicate features near 28 and 33 um indicating significant processing when compared to interstellar grains. A few sources show large 10-to-20 um ratios and require even larger grains emitting at 20 um than at 10 um. This size difference may arise from the difference in the depth into the disk probed by the two silicate emission bands in disks where dust settling has occurred. The 10 um feature strength vs. shape trend is not correlated with age or Halpha equivalent width, suggesting that some amount of turbulent mixing and regeneration of small grains is occurring. The strength vs. shape trend is related to spectral type, however, with M stars showing significantly flatter 10 um features (larger grain sizes) than A/B stars. The connection between spectral type and grain size is interpreted in terms of the variation in the silicate emission radius as a function of stellar luminosity, but could also be indicative of other spectral-type dependent factors (e.g, X-rays, UV radiation, stellar/disk winds, etc.).

The structure and chemistry of protoplanetary disks depends strongly on the nature of the central star around which it has formed. The dust temperature is mainly set by the stellar luminosity, while the chemistry of the upper disk layers depends on the amount of intercepted UV and X-ray flux. We will study the differences in chemistry, thermal sturcture and line emission around Herbig Ae/Be, T Tauri stars and low mass M dwarfs. Predictions will be made for future observations with SOFIA and Herschel.

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.

We present simulations of the interferometric visibilities of Herbig Ae star disks. We investigate whether interferometric measurements in the 10 micrometer atmospheric window are sensitive to the presence of an increased scale height at the inner disk edge, predicted by recent models. Furthermore, we investigate whether such measurements can discriminate between disks with a ``flaring'' geometry and disks with a ``flat'' geometry. We show that both these questions can be addressed, using measurements at a small number of appropriately chosen baselines. The classification of Herbig Ae stars in two groups, based on the appearance of the spectral energy distribution (SED), has been attributed to a difference in disk geometry. Sources with a group I SED would have a flaring outer disk geometry, whereas the disk of group II sources is proposed to be flat (or ``self-shadowed''). We show that this hypothesis can be tested using long-baseline interferometric measurements in the micrometer atmospheric window.

We present high spectral-resolution optical spectra of 49 Herbig Ae/Be stars in a search for the [O I] 6300A line. The vast majority of the stars in our sample show narrow (FWHM < 100 km/s) emission lines, centered on the stellar radial velocity. In only three sources is the feature much broader (~400 km/s), and strongly blueshifted (-200 km/s) compared to the stellar radial velocity. Some stars in our sample show double-peaked lines profiles, with peak-to-peak separations of ~10 km/s. The presence and strength of the [O I] line emission appears to be correlated with the far-infrared energy distribution of each source: stars with a strong excess at 60 micron have in general stronger [O I] emission than stars with weaker 60 micron excesses. We interpret these narrow [O I] 6300A line profiles as arising in the surface layers of the protoplanetary disks surrounding Herbig Ae/Be stars. A simple model for [O I] 6300A line emission due to the photodissociation of OH molecules shows that our results are in quantitative agreement with that expected from the emission of a flared disk if the fractional OH abundance is ~5 x 10^{-7}.

We use a 3D Monte Carlo radiative transfer code to study the projection of large shadows by circumstellar disks around young stellar objects on surrounding reflection nebulosity. It is shown that for a wide range of parameters a small (10-100 AU) circumstellar disk can project a large (1 000-10 000 AU) dark band in the near-infrared that often resembles a massive edge-on disk. The disk shadows are divided into two basic types, depending on the distribution of the reflecting material and the resulting morphology of the shadows in the near-infrared. Two YSOs associated with bipolar nebulosity, CK 3/EC 82 illuminating the Serpens Reflection Nebula (SRN) and Ced 110 IRS 4 in the Chamaeleon I molecular cloud, are modelled in detail as disk shadows. Spectral energy distributions of the two sources are collected using both archival ISO data and new Spitzer-IRS data. An axisymmetric model consisting of a small disk and a spherically symmetric envelope can reproduce the near-infrared images and full spectral energy distributions of the two disk shadow candidates. It is shown that the model fits can be used to constrain the geometry of the central disks due to the magnifying effect of the projection. We find that a disk unresolved in near-infrared images, but casting a large disk shadow, can be modelled at a level of sophistication approaching that of an edge-on disk with resolved near-infrared images. It is found that the most obvious observable difference between a disk shadow and a large optically thick disk is that the disk shadows have a compact near-infrared source near the center of the dark band. High resolution imaging and/or polarimetry should reveal the compact source in the center of a disk shadow as an edge-on disk. [Abstract abridged]

We discuss the question whether the matter in dusty tori around active galactic nuclei has a smooth or a clumpy structure. Nenkova, Ivezic & Elitzur (2002) have argued that the lack of emission feature in the SEDs of type 1 AGN galaxies combined with a clear absorption feature in type 2 AGN can be explained if the circumnuclear dust is distributed in discrete clumps. Our aim is to verify this. We use multi-dimensional radiative transfer models of smooth and clumpy tori, and compare the SEDs of equivalent smooth and clumpy models. We find that the 10 micron emission feature of the clumpy models, when seen almost face-on, is not appreciably reduced compared to the equivalent smooth models. Some of the clumpy models have a weak or even absent 10 micron feature, but so do some of the smooth models. On the whole the SEDs of clumpy and smooth tori are similar, but some details are different. The absorption feature seen at edge-on inclinations appears to be less deep in the clumpy models than in the smooth models, and the average flux in the near-infrared regime is stronger in the clumpy models. Moreover, at these inclinations the clumpy models have a slightly wider SED. Whether these differences are unique enough to be used as a diagnostic for clumpiness of AGN tori is not yet clear.

We calculate the evaporation of a cool accretion disk around a black hole due to the ion-bombardment by an ion supported accretion flow (here ISAF, or optically thin ADAF). As first suggested by Spruit & Deufel (2002), this evaporation takes place in two stages: ion bombardment of the cool disk (Shakura-Sunyaev disk: SSD) produces an intermediate-temperature layer on top of the disk (`warm layer') which constitutes an independent accretion flow on both sides of the SSD. As this warm material accretes inward of the inner radius of the SSD, it becomes thermally unstable by lack of cooling of photons, and evaporates into the ISAF, thereby feeding the latter. Angular momentum conservation forces a certain fraction of the ISAF material to move outward, where it can bombard the SSD with its hot ions. The flow geometry is derived by computing stationary solutions of the continuity- and angular momentum equations for the three components (ISAF, warm flow and SSD). The overall radiative output is dominated by hard X-rays. They are produced mostly from the warm component, rather than the ISAF. The expected time dependence and stability of the flow, not computed here, is discussed briefly.

We have analyzed the total and nuclear SED for NGC 4261 and find that the dominant process for the mid- and far-infrared emission in this object is non-thermal emission from the active nucleus. Modeling the emission from the optically detected 300 pc dust disk yields no significant disk contribution at any wavelength. To explain the observations, either the disk has an inflated inner region which partly absorbs the core, or the intrinsic core spectrum is curved. The inner 10 pc of the disk can potentially be conceived as an obscuring torus, albeit with optical depth around unity.

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.

(Abriged) We present the results of millimeter observations and a suitable chemical and radiative transfer model of the AB Aur (HD 31293) circumstellar disk and surrounding envelope. The integral molecular content of this system is studied by observing CO, C$^{18}$O, CS, HCO$^+$, DCO$^+$, H$_2$CO, HCN, HNC, and SiO rotational lines with the IRAM 30-m antenna, while the disk is mapped in the HCO$^+$(1-0) transition with the Plateau de Bure interferometer. Using a flared disk model with a vertical temperature gradient and an isothermal spherical envelope model with a shadowed midplane and two unshielded cones together with a gas-grain chemical network, time-dependent abundances of observationally important molecules are calculated. Then a 2D non-LTE line radiative transfer code is applied to compute excitation temperatures of several rotational transitions of HCO$^+$, CO, C$^{18}$O, and CS molecules. We synthesize the HCO$^+$(1-0) interferometric map along with single-dish CO(2-1), C$^{18}$O(2-1), HCO$^+$(1-0), HCO$^+$(3-2), CS(2-1), and CS(5-4) spectra and compared them with the observations. Our disk model successfully reproduces observed interferometric HCO$^+$(1-0) data, thereby constraining the following disk properties: (1) the inclination angle $\iota=17^{+6}_{-3}\degr$, (2) the position angle $\phi=80\pm30\degr$, (3) the size $R_\mathrm{out}=400\pm200$ AU, (4) the mass $M_\mathrm{disk}=1.3\cdot10^{-2} M_{\sun}$ (with a factor of $\sim7$ uncertainty), and (5) that the disk is in Keplerian rotation. Furthermore, indirect evidence for a local inhomogeneity of the envelope at $\ga600$ AU is found...

We present 5.2-37.2 micron spectroscopy of the edge-on circumstellar disk CRBR 2422.8-3423 obtained using the InfraRed Spectrograph (IRS) of the Spitzer Space Telescope. The IRS spectrum is combined with ground-based 3-5 micron spectroscopy to obtain a complete inventory of solid state material present along the line of sight toward the source. We model the object with a 2D axisymmetric (effectively 3D) Monte Carlo radiative transfer code. It is found that the model disk, assuming a standard flaring structure, is too warm to contain the very large observed column density of pure CO ice, but is possibly responsible for up to 50% of the water, CO2 and minor ice species. In particular the 6.85 micron band, tentatively due to NH4+, exhibits a prominent red wing, indicating a significant contribution from warm ice in the disk. It is argued that the pure CO ice is located in the dense core Oph-F in front of the source seen in the submillimeter imaging, with the CO gas in the core highly depleted. The model is used to predict which circumstances are most favourable for direct observations of ices in edge-on circumstellar disks. Ice bands will in general be deepest for inclinations similar to the disk opening angle, i.e. ~70 degrees. Due to the high optical depths of typical disk mid-planes, ice absorption bands will often probe warmer ice located in the upper layers of nearly edge-on disks. The ratios between different ice bands are found to vary by up to an order of magnitude depending on disk inclination due to radiative transfer effects caused by the 2D structure of the disk. Ratios between ice bands of the same species can therefore be used to constrain the location of the ices in a circumstellar disk. [Abstract abridged]

A 850 micron map of the interacting spiral galaxy M51 shows well-defined spiral arms, closely resembling the structures seen in CO and HI emission. However, most of the 850 micron emission originates in an underlying exponential disk, a component that has not been observed before in a face-on galaxy at these wavelengths. The scale-length of this disk is 5.45 kpc, which is somewhat larger than the scale-length of the stellar disk, but somewhat smaller than that of atomic hydrogen. Its profile can not be explained solely by a radial disk temperature gradient but requires the underlying dust to have an exponential distribution as well. This reinforces the view that the submm emission from spiral galaxy disks traces total hydrogen column density, i.e.the sum of H2 and HI. A canonical gas-to-dust ratio of 100+/-26 is obtained for kappa(850)=1.2 g**-1 cm**2, where kappa(850) is the dust opacity at 850 micron.

We present accurate mid-infrared observations of the disk around the young, bona-fide brown dwarf CFHT-BD-Tau 4. We report GEMINI/T-ReCS measurements in the 7.9, 10.4 and 12.3 micron filters, from which we infer the presence of a prominent, broad silicate emission feature. The shape of the silicate feature is dominated by emission from 2 micron amorphous olivine grains. Such grains, being an order of magnitude larger than those in the interstellar medium, are a first proof of dust processing and grain growth in disks around brown dwarfs. The object's spectral energy distribution is below the prediction of the classical flared disk model but higher than that of the two-layer flat disk. A good match can be achieved by using an intermediate disk model with strongly reduced but non-zero flaring. Grain growth and dust settling processes provide a natural explanation for this disk geometry and we argue that such intermediate flaring might explain the observations of several other brown dwarf disks as well.

We have imaged two bona-fide brown dwarfs with TReCS/GEMINI-S and find mid-infrared excess emission that can be explained by optically thick dust disk models. In the case of the young ($\approx$2Myr) Cha H$\alpha$1 we measure fluxes at 10.4$\mu$m and 12.3$\mu$m that are fully consistent with a standard flared disk model and prominent silicate emission. For the $\approx$ 10Myr old brown dwarf 2MASS1207-3932 located in the TW Hydrae association we find excess emission at 8.7$\mu$m and 10.4$\mu$m with respect to its photosphere, and confirm disk accretion as likely cause of its strong activity. Disks around brown dwarfs likely last at least as long as their low-mass stellar counterparts in the T-Tauri phase. Grain growth, dust settling, and evolution of the geometry of brown dwarfs disks may appear on a timescale of 10Myr and can be witnessed by observations in the mid-infrared.

Models for the structure of protoplanetary disks have so far been based on the assumption that the gas and the dust temperature are equal. The gas temperature, an essential ingredient in the equations of hydrostatic equilibrium of the disk, is then determined from a continuum radiative transfer calculation, in which the continuum opacity is provided by the dust. It has been long debated whether this assumption still holds in the surface layers of the disk, where the dust infrared emission features are produced. In this paper we compute the temperature of the gas in the surface layers of the disk in a self-consistent manner. The gas temperature is determined from a heating-cooling balance equation in which processes such as photoelectric heating, dissociative heating, dust-gas thermal heat exchange and line cooling are included. The abundances of the dominant cooling species such as CO, C, C+ and O are determined from a chemical network based on the atomic species H, He, C, O, S, Mg, Si, Fe (Kamp & Bertoldi 2000). The underlying disk models to our calculations are the models of Dullemond, van Zadelhoff & Natta (2002). We find that in general the dust and gas temperature are equal to withing 10% for A_V >~ 0.1, which is above the location of the `super-heated surface layer' in which the dust emission features are produced (e.g. Chiang & Goldreich 1997). High above the disk surface the gas temperature exceeds the dust temperature and can can become -- in the presence of polycyclic aromatic hydrocarbons -- as high as 600 K at a radius of 100 AU. This is a region where CO has fully dissociated, but a significant fraction of hydrogen is still in molecular form. The densities are still high enough for non-negligible H_2 emission to be produced.....(see paper for full abstract)

We present high spectral-resolution optical spectra of 49 Herbig Ae/Be stars in a search for the [O I] 6300.2 A line. The vast majority of the stars in our sample show narrow (FWHM < 100 km/s) emission lines, centered on the stellar radial velocity. Some stars in our sample show double-peaked lines profiles, with peak-to-peak separations of about 10 km/s. The presence and strength of the [O I] line emission appears to be well correlated with the far-infrared energy distribution of each source: stars with a strong excess at 60 micron have in general stronger [O I] emission than stars with weaker 60 micron excesses. We interpret the observed [O I] 6300.2 A line profiles as arising in the surface layers of the protoplanetary disks surrounding Herbig Ae/Be stars.

We present observations of L1014, a dense core in the Cygnus region previously thought to be starless, but data from the Spitzer Space Telescope shows the presence of an embedded source. We propose a model for this source that includes a cold core, heated by the interstellar radiation field, and a low-luminosity internal source. The low luminosity of the internal source suggests a substellar object. If L1014 is representative, other "starless" cores may turn out to harbor central sources.

Sensitive 5-38 um Spitzer Space Telescope (SST) and ground based 3-5 um spectra of the embedded low mass protostars B5 IRS1 and HH46 IRS show deep ice absorption bands superposed on steeply rising mid-infrared continua. The ices likely originate in the circumstellar envelopes. The CO2 bending mode at 15 um is a particularly powerful tracer of the ice composition and processing history. Toward these protostars, this band shows little evidence for thermal processing at temperatures above 50 K. Signatures of lower temperature processing are present in the CO and OCN- bands, however. The observed CO2 profile indicates an intimate mixture with H2O, but not necessarily with CH3OH, in contrast to some high mass protostars. This is consistent with the low CH3OH abundance derived from the ground based L band spectra. The CO2/H2O column density ratios are high in both B5 IRS1 and HH46 IRS (~35%). Clearly, the SST spectra are essential to study ice evolution in low mass protostellar environments, and to eventually determine the relation between interstellar and solar system ices.

We analyze how the process of dust settling affects the spectral energy distribution and optical appearance of protoplanetary disks. Using simple analytic estimates on the one hand, and detailed 1+1-D models on the other hand, we show that, while the time scale for settling down to the equator may exceed the life time of the disk, it takes much less time for even small grains of 0.1 micron to settle down to a few pressure scale heights. This is often well below the original location of the disk's photosphere, and the disk therefore becomes effectively 'flatter'. If turbulent stirring is included, a steady state solution can be found, which is typically reached after a few times 1E5 years. Dependent on the strength of the turbulence, the shape of the disk in such a steady state can be either fully flaring, or flaring only up to a certain radius and self-shadowed beyond that radius. We show that these partly self-shadowed disks have a much weaker mid- to far-infrared flux than the fully flaring ones. We also show that these self-shadowed regions of the disk are very weak in resolved images of scattered light, in contrast to the fully flaring disks. From the calculations with compact grains it follows that, after about 1E6 years, most disks should be self-shadowed. The fact that some older disks are still observed with the characteristics of flaring disks therefore seems somewhat inconsistent with the time scales predicted by the settling model based on compact grains. This suggests that perhaps even the small grains (lesssim 0.1 micron) have a porous or fractal structure, or that the different geometries of observed disks is merely a reflection of the turbulent state of these disks.

We have calculated the (sub-)mm spectral indices of 26 Herbig Ae/Be stars, for which we can determine the infrared spectral energy distribution (SED). We find a clear correlation between the strength of the ratio of the near- to mid-infrared excess of these sources, and the slope of the (sub-)mm energy distribution. Based on earlier multi-dimensional modeling of disks around Herbig Ae stars, we interpret this as a correlation between the geometry of the disk (flared or self-shadowed) and the size of the grains: self-shadowed disks have, on average, larger grains than their flared counterparts. These data suggest that the geometry of a young stellar disk evolves from flared to self-shadowed.

We present benchmark problems and solutions for the continuum radiative transfer (RT) in a 2D disk configuration. The reliability of three Monte-Carlo and two grid-based codes is tested by comparing their results for a set of well-defined cases which differ for optical depth and viewing angle. For all the configurations, the overall shape of the resulting temperature and spectral energy distribution is well reproduced. The solutions we provide can be used for the verification of other RT codes.We also point out the advantages and disadvantages of the various numerical techniques applied to solve the RT problem.

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.

In this Letter we propose a new view on UX Orionis type variability. The idea is based on the earlier proposal by various authors that UXORs are nearly-edge-on disks in which hydrodynamic fluctuations could cause clumps of dust and gas to cross the line of sight. However, because the standard disk models have a flaring geometry, it is mostly the outer regions of the disk that obscure the star. The time scales for such obscuration events would be too long to match the observed time scales of weeks to months. Recent 2-D self-consistent models of Herbig Ae/Be protoplanetary disks (Dullemond et al. 2002,2003 henceforth D02/DD03), however, have indicated that for Herbig Ae/Be star disks there exists, in addition to the usual flared disks, also a new class of disks: those that are fully self-shadowed. Only their puffed-up inner rim (at the dust evaporation radius) is directly irradiated by the star, while the disk at larger radius resides in the shadow of the rim. For these disks there exist inclinations at which the line of sight towards the star skims the upper parts of the puffed-up inner rim, while passing high over the surface of outer disk regions. Small hydrodynamic fluctuations in the puffed-up inner rim could then be held responsible for the extinction events seen in UXORs. If this idea is correct, it makes a prediction for the shape of the SEDs of these stars. It was shown by D02/DD03 that flared disks have a strong far-IR excess and can be classified as `group I' (in the classification of Meeus et al. 2001), while self-shadowed disks have a relatively weak far-IR excess and are classified as `group II'. Our model therefore predicts that UXORs belong to the `group II' sources. We show that this correlation is indeed found within a sample of 86 Herbig Ae/Be stars.

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.

The combination of mid-infrared and recent submm/mm measurements allows us to set up the first comprehensive spectral energy distribution (SED) of the circumstellar material around a young Brown Dwarf. Simple arguments suggest that the dust is distributed in the form of a disk. We compare basic models to explore the disk parameters. The modeling shows that a flat disk geometry fits well the observations. A flared disk explains the SED only if it has a puffed-up inner rim and an inner gap much larger than the dust sublimation radius. Similarities and differences with disks around T Tauri stars are discussed.

The two-layer disk models of Chiang & Goldreich (1997, henceforth CG) and its derivatives are popular among astronomers because of their simplicity and the clear predictions they make for the SEDs of T Tauri stars and Herbig Ae/Be stars. Moreover, they can be computed quickly, which is a great advantage when fitting observations using automated procedures. In this paper we wish to assess the accuracy and reliability of 2-layer models, by comparing them to detailed vertical structure models with accurate 1+1D radiative transfer. We focus on the shape of the SED, and the predicted height and "flaring index" of the disk. We first consider models where scattering is set to zero. We find that 2-layer models overestimate significantly the near-infrared flux, and we suggest a simple way of correcting this effect, at least in part. At longer wavelengths, the SED of two-layer models often show a two-bump structure, which is absent in 1+1D models. Nevertheless, overall agreement is reasonably good, and the differences are in most cases within 30%. At (sub)-mm wavelengths the differences may even be less. The shape of the disk, as measured by its pressure and surface scale height and by the flaring angle are also well reproduced by two-layer models. When scattering is included in the 1+1D models, the differences become larger, especially in the near-infrared. We suggest simple ways to include scattering in two-layer models and discuss their reliability. We do not compare the two-layer models to full 2-D/3-D models, so the conclusions remain valid only within the annulus-by-annulus approximation.

Two-dimensional radiative transfer is employed to obtain the broad-band infrared spectrum of active galaxies. In the models we vary the geometry and size of the obscuring medium, the surface density, the opacity and the grain size distribution. Resulting spectral energy distributions are constructed for different orientations of the toroid. Colour-colour comparisons with observational data are consistent with previous observations that the emission longward of 60 micron is produced by star-formation and unrelated to the presence of an obscuring torus. We also find that the toroid cannot be flat, but is rather conical or flaring. The density is most likely constant with radius, and the size is relatively large with an inner radius around 10 pc. A direct comparison with radio galaxy Cygnus A yields a best fit for a conical disk with constant surface density, and a size from 10 to 30 pc assuming the far-infrared emission is due to star-formation in the host galaxy.

We discuss spectral energy distributions of a sample of Herbig Ae/Be stars in the context of a passive irradiated disk model. The data have been presented earlier by Meeus et al (2001), and preliminary interpretations of these data were given in that paper. While the spectra of Herbig Ae stars all show similarities, there is significant variation between the spectra, in particular in the shape of the mid-IR rise and in the presence or absence of a silicate feature. We explore the hypothesis that all these different spectra can be interpreted as pure disk spectra without additional components. Using the model of Dullemond, Dominik and Natta (2001) we deduce the disk parameters of a number of the sources, and find that for a large fraction of investigated sources, satisfactory fits can be obtained. The derived model parameters show that some group Ia sources can only be fit with radially increasing surface densities, indicating the presence of depleted inner disk regions. The steep-sloped SEDs of group IIa sources can be fit with very compact disks, probably representing disks with collapsed outer regions. The largest difficulties arise from sources that do not show significant silicate emission features. Our attempts to explain these objects with a pure geometric effect are only partially successful. It seems that these stars indeed require a strong depletion of small silicate grains.

We present the Spectral Energy Distribution of HR 4049 based on literature data and new continuum measurements at 850 micron. The SED shows variable absorption in the UV, and a large IR excess, both caused by circumstellar dust. The shape of the IR excess from 1 micron all the way down to 850 micron can be nearly perfectly fitted with a single blackbody function at T ~ 1150 K or alternatively with a sum of blackbodies in a narrow temperature range. The energy emitted in this IR continuum radiation is about one-third of the stellar luminosity. We show that this blackbody radiation must be due to the presence of a circumbinary disk with a large height. This disk must also be gas-rich, in agreement with the observations of molecular bands in the ISO-SWS spectrum. We present two possible scenario's for explaining the shape and the intensity of the IR excess. The first scenario involves large grains (a >1 mm) that each radiate like a blackbody. The second scenario argues that the blackbody radiation is due to a very optically thick circumbinary disk. We investigate if such a disk would indeed produce blackbody radiation by presenting results from radiative transfer calculations. We further quantify the properties of such a disk and its stability in the framework of (hydro)dynamics, grain settling, radiation pressure and grain drift. The virtues and shortcomings of both models for the origin of the IR blackbody are discussed by contrasting them with other observations and assessing them in the framework of (binary) (post-)AGB evolution.

In this paper the two-dimensional structure of protoplanetary disks around Herbig Ae/Be stars is studied. This is done by constructing a self-consistent model based on 2-D radiative transfer coupled to the equation of vertical hydrostatics. As a simplifying assumption a grey opacity is used. It is found that the disk can adopt four different structures, dependent on the surface density distribution Sigma(R) as a function of radius, i.e. on radial- and vertical optical depth of the disk. For the case of high to intermediate vertical optical depth, the temperature and density structures are in agreement with the simple ``disk with inner hole'' model of Dullemond, Dominik & Natta (2001, henceforth DDN01). At large radii the disk adopts a flaring shape as expected, and near the dust destruction radius (located at about 0.5 AU for most Herbig Ae stars) the disk is superheated and puffs up. The region directly behind this ``puffed-up inner dust wall'' is shadowed, as predicted by DDN01. For the case of intermediate to low vertical optical depth, but still high radial optical depth, the 2-D models show that the shadow can cover the entire disk. The disk is then no longer flaring. When the disk is optically thin even in radial direction, it becomes unshadowed again because the inner rim can no longer block the stellar light. Finally, for disks that flare at intermediate radii, but become too optically thin at large radii, the outer parts once again become shadowed. But this time the shadowing is caused by the flaring part of the disk instead of the inner rim. Observational consequences of these four geometries are discussed.

Comparison is made between a number of independent computer programs for radiative transfer in molecular rotational lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO+. For the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (tau=4800). The problem of the collapsing cloud in HCO+ has a larger spread in results, ranging up to 12% for the J=4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO+ J=4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.

In this paper we present detailed models of the vertical structure (temperature and density) of passive irradiated circumstellar disks around T Tauri and Herbig Ae/Be stars. In contrast to earlier work, we use full frequency- and angle-dependent radiative transfer instead of the usual moment equations. We find that this improvement of the radiative transfer has strong influence on the resulting vertical structure of the disk, with differences in temperature as large as 70 %. However, the spectral energy distribution (SED) is only mildly affected by this change. In fact, the SED compares reasonably well with that of improved versions of the Chiang & Goldreich (CG) model. This shows that the latter is a reasonable model for the SED, in spite of its simplicity. It also shows that from the SED alone, little can be learned about the vertical structure of a passive circumstellar disk. The molecular line emission from these disks is more sensitive to the vertical temperature and density structure, and we show as an example how the intensity and profiles of various CO lines depend on the adopted disk model. The models presented in this paper can also serve as the basis of theoretical studies of e.g. dust coagulation and settling in disks.

The X-ray spectrum from a cool accretion disk heated by virialized protons is computed. The cool disk is either embedded in a magnetically heated accretion disk corona or partly extends into an ion supported torus (or ADAF). We calculate the stationary equilibrium between proton heating, electron thermal conduction and the radiative losses by bremsstrahlung and Compton scattering. A heated surface layer on top of the accretion disk is produced with temperatures between 60--90 keV above a cool layer with temperatures of 0.01 keV (AGN) and 1keV (galactic black hole candidates). The spectra produced by the surface layer are reminiscent of hard state spectra, but a bit too steep, especially for AGN's. Near the inner edge of the disk, where the optical depth of the disk $\tau lesssim 1$, we find that the cool component of the disk disappears. Instead, the hot protons from the corona/ADAF heat the disk, on a dynamical time-scale, to temperatures of several 100 keV, limited by pair production. This region, here called a `warm disk', could contribute significantly to the hard X-ray spectra and could be important for feeding material into an ADAF.

We consider the interaction of a slowly rotating unmagnetized neutron star with a hot (ion supported, ADAF) accretion flow. The virialized protons of the ADAF penetrate into the neutron star atmosphere, heating a surface layer. Detailed calculations are presented of the equilibrium between heating by the protons, electron thermal conduction, bremsstrahlung and multiple Compton scattering in this layer. Its temperature is of the order 40-70 keV. Its optical depth increases with the incident proton energy flux, and is of the order unity for accretion at $10^{-2}$--$10^{-1}$ of the Eddington rate. At these rates, the X-ray spectrum produced by the layer has a hard tail extending to 100 keV, and is similar to the observed spectra of accreting neutron stars in their hard states. The steep gradient at the base of the heated layer gives rise to an excess of photons at the soft end of the spectrum (compared to a blackbody) through an `inverse photosphere effect'. The differences with respect to previous studies of similar problems are discussed, they are due mostly to a more accurate treatment of the proton penetration process and the vertical structure of the heated layer.

A model for irradiated dust disks around Herbig Ae stars is proposed. The model is based on the flaring disk model of Chiang & Goldreich (1997), but with the central regions of the disk removed. The inner rim of the disk is puffed up and is much hotter than the rest of the disk, because it is directly exposed to the stellar flux. If located at the dust evaporation radius, its reemitted flux produces a conspicuous bump in the SED which peaks at 2-3 micron. We propose that this emission is the explanation for the near-infrared bump observed in the SEDs of Herbig Ae stars. We study for which stellar parameters this bump would be observable, and find that it is the case for Herbig Ae stellar parameters but not for T-Tauri stars, confirming what is found from the observations. We also study the effects of the shadow cast by the inner rim over the rest of the flaring disk. The shadowed region can be quite large, and under some circumstances the entire disk may lie in the shadow. This shadowed region will be much cooler than an unshadowed flaring disk, since its only heating sources are radial radiative diffusion and possible indirect sources of irradiation. Under certain special circumstances the shadowing effect can suppress, or even completely eliminate, the 10 micron emission feature from the spectrum, which might explain the anomalous SEDs of some isolated Herbig Ae stars in the sample of Meeus et al. (2001). At much larger radii the disk emerges from the shadow, and continues as a flaring disk towards the outer edge. The complete model, including structure of the inner edge, shadowed region and the flared outer part, is described in detail in this paper, and we show examples of the general behavior of the model for varying parameters.

The uniqueness and stability of irradiated flaring passive protostellar disks is investigated in the context of a simplified set of equations for the vertical height H as a function of radius R. It is found that the well-known flaring disk solution with H ~ R^{9/7} is not unique. Diverging solutions and asymptotically conical (H ~ R) solutions are also found. Moreover, using time-dependent linear perturbation analysis, it is found that the flaring disk solution may become unstable to self-shadowing. A local enhancement in the vertical height alters the functional form of irradiation grazing angle, and causes the 'sunny side' of the enhancement to grow and the 'shadow side' to collapse in a run-away fashion. This instability operates in regions of the disk in which the cooling time is much shorter than the vertical sound crossing time, which may occur in the outer regions of the passive irradiated disk if dust and gas are sufficiently strongly thermally coupled. Processes that may stabilize the disk, which include active accretion, irradiation from above (e.g. a scattering corona) and low disk optical depth, are likely to operate only at small or at large radius. The simple analysis of this Letter therefore suggests that the instability may alter the flaring disk structure at intermediate radii (between the actively accreting and fast rotating inner regions and the optically thin outer regions).

We present an algorithm for two-dimensional radiative transfer in axisymmetric, circumstellar media. The formal integration of the transfer equation is performed by a generalization of the short characteristics (SC) method to spherical coordinates. Accelerated Lambda Iteration (ALI) and Ng's algorithm are used to converge towards a solution. By taking a logarithmically spaced radial coordinate grid, the method has the natural capability of treating problems that span several decades in radius, in the most extreme case from the stellar radius up to parsec scale. Flux conservation is guaranteed in spherical coordinates by a particular choice of discrete photon directions and a special treatment of nearly-radially outward propagating radiation. The algorithm works well from zero up to very high optical depth, and can be used for a wide variety of transfer problems, including non-LTE line formation, dust continuum transfer and high temperature processes such as compton scattering. In this paper we focus on multiple scattering off dust grains and on non-LTE transfer in molecular and atomic lines. Line transfer is treated according to an ALI scheme for multi-level atoms/molecules, and includes both random and systematic velocity fields. The algorithms are implemented in a multi-purpose user-friendly radiative transfer program named RADICAL. We present two example computations: one of dust scattering in the Egg Nebula, and one of non-LTE line formation in rotational transitions of HCO$^{+}$ in a flattened protostellar collapsing cloud.

In this Letter we investigate the properties of advection-dominated accretion flows (ADAFs) fed by the evaporation of a Shakura-Sunyaev accretion disk (SSD). In our picture the ADAF fills the central cavity evacuated by the SSD and extends beyond the transition radius into a coronal region. We find that, because of global angular momentum conservation, a significant fraction of the hot gas flows away from the black hole forming a transsonic wind, unless the injection rate depends only weakly on radius (if $r^2\dot\sigma\propto r^{-\xi}$, $\xi< 1/2$). The Bernoulli number of the inflowing gas is negative if the transition radius is $\lesssim 100$ Schwarzschild radii, so matter falling into the hole is gravitationally bound. The ratio of inflowing to outflowing mass is $\approx 1/2$, so in these solutions the accretion rate is of the same order as in standard ADAFs and much larger than in advection-dominated inflow/outflow models (ADIOS). The possible relevance of evaporation-fed solutions to accretion flows in black hole X-ray binaries is briefly discussed.

The conditions for disk evaporation by electron thermal conduction are examined, using a simplified semi--analytical 1-D model. The model is based on the mechanism proposed by Meyer & Meyer-Hofmeister (1994) in which an advection dominated accretion flow evaporates the top layers from the underlying disk by thermal conduction. The evaporation rate is calculated as a function of the density of the advective flow, and an analysis is made of the time scales and length scales of the dynamics of the advective flow. It is shown that evaporation can only completely destroy the disk if the conductive length scale is of the order of the radius. This implies that radial conduction is an essential factor in the evaporation process. The heat required for evaporation is in fact produced at small radii and transported radially towards the evaporation region.

The existence of bimodal disks is investigated. Following a simple argument based on energetic considerations we show that stationary, bimodal accretion disk models in which a Shakura--Sunyaev disk (SSD) at large radii matches an advection dominated accretion flow (ADAF) at smaller radii are never possible using the standard slim disk approach, unless some extra energy flux is present. The same argument, however, predicts the possibility of a transition from an outer Shapiro--Lightman--Eardley (SLE) disk to an ADAF, and from a SLE disk to a SSD. Both types of solutions have been found.