Context. Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their counterparts in the solar system are the asteroid and Edgeworth-Kuiper belts. Aims. The aim of this paper is to provide robust numbers for the incidence of debris discs around FGK stars in the solar neighbourhood. Methods. The full sample of 177 FGK stars with d <= 20 pc proposed for the DUst around Nearby Stars (DUNES) survey is presented. Herschel/PACS observations at 100 and 160 mu m were obtained, and were complemented in some cases with data at 70 mu m and at 250, 350, and 500 mu m SPIRE photometry. The 123 objects observed by the DUNES collaboration were presented in a previous paper. The remaining 54 stars, shared with the Disc Emission via a Bias-free Reconnaissance in IR and Sub-mm (DEBRIS) consortium and observed by them, and the combined full sample are studied in this paper. The incidence of debris discs per spectral type is analysed and put into context together with other parameters of the sample, like metallicity, rotation and activity, and age. Results. The subsample of 105 stars with d <= 15 pc containing 23 F, 33 G, and 49 K stars is complete for F stars, almost complete for G stars, and contains a substantial number of K stars from which we draw solid conclusions on objects of this spectral type. The incidence rates of debris discs per spectral type are 0.26(-0.14)(+0.21) (6 objects with excesses out of 23 F stars), 0.21(-0.11)(+0.17) (7 out of 33 G stars), and 0.20(-0.09)(+0.14) (10 out of 49 K stars); the fraction for all three spectral types together is 0.22(-0.07)(+0.08) (23 out of 105 stars). The uncertainties correspond to a 95% confidence level. The medians of the upper limits of L-dust/L-* for each spectral type are 7.8 x 10(-7) (F), 1.4 x 10(-6) (G), and 2.2 x 10(-6) (K); the lowest values are around 4.0 x 10(-7). The incidence of debris discs is similar for active (young) and inactive (old) stars. The fractional luminosity tends to drop with increasing age, as expected from collisional erosion of the debris belts.

Using Keck/HIRES spectra (Delta v similar to 7 km s(-1)) we analyze forbidden lines of [O I] 6300 angstrom, [O I] 5577 angstrom. and [S II] 6731 angstrom. from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high-velocity component (HVC) associated with microjets, we study the properties of the low-velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (magnetohydrodynamic) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material. LVC emission is seen in all 30 stars with detected [O. I] but only in two out of eight with detected [S. II], so our analysis is largely based on the properties of the [O. I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks. The full width at half maximum of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 au and 0.5 to 5 au, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC, however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive.

We develop a simple model to predict the radial distribution of planetesimal formation. The model is based on the observed growth of dust to millimeter-sized particles, which drift radially, pile-up, and form planetesimals where the stopping time and dust-to-gas ratio intersect the allowed region for streaming instability-induced gravitational collapse. Using an approximate analytic treatment, we first show that drifting particles define a track in metallicity-stopping time space whose only substantial dependence is on the disk's angular momentum transport efficiency. Prompt planetesimal formation is feasible for high particle accretion rates (relative to the gas, (M) over dot(p)/(M) over dot greater than or similar to 3 x 10(-2) for alpha = 10(-2)), which could only be sustained for a limited period of time. If it is possible, it would lead to the deposition of a broad and massive belt of planetesimals with a sharp outer edge. Numerically including turbulent diffusion and vapor condensation processes, we find that a modest enhancement of solids near the snow line occurs for centimeter-sized particles, but that this is largely immaterial for planetesimal formation. We note that radial drift couples planetesimal formation across radii in the disk, and suggest that considerations of planetesimal formation favor a model in which the initial deposition of material for giant planet cores occurs well beyond the snow line.

Context. High-contrast scattered light observations have revealed the surface morphology of several dozen protoplanetary disks at optical and near-infrared wavelengths. Inclined disks offer the opportunity to measure part of the phase function of the dust grains that reside in the disk surface which is essential for our understanding of protoplanetary dust properties and the early stages of planet formation. Aims. We aim to construct a method which takes into account how the flaring shape of the scattering surface of an optically thick protoplanetary disk projects onto the image plane of the observer. This allows us to map physical quantities (e.g., scattering radius and scattering angle) onto scattered light images and retrieve stellar irradiation corrected images (r(2)-scaled) and dust phase functions. Methods. The scattered light mapping method projects a power law shaped disk surface onto the detector plane after which the observed scattered light image is interpolated backward onto the disk surface. We apply the method on archival polarized intensity images of the protoplanetary disk around HD 100546 that were obtained with VLT/SPHERE in the R' band and VLT/NACO in the H and K-s bands. \Results. The brightest side of the r(2)-scaled R-0 band polarized intensity image of HD 100546 changes from the far to the near side of the disk when a flaring instead of a geometrically flat disk surface is used for the r(2)-scaling. The decrease in polarized surface brightness in the scattering angle range of similar to 40 degrees-70 degrees is likely a result of the dust phase function and degree of polarization which peak in different scattering angle regimes. The derived phase functions show part of a forward scattering peak, which indicates that large, aggregate dust grains dominate the scattering opacity in the disk surface. Conclusions. Projection effects of a protoplanetary disk surface need to be taken into account to correctly interpret scattered light images. Applying the correct scaling for the correction of stellar irradiation is crucial for the interpretation of the images and the derivation of the dust properties in the disk surface layer.

Thesis (Ph. D.)--University of Washington, 2005. / Vita. Includes bibliographical references (p. 166-174).

Aims. We numerically studied close encounters between a young stellar system hosting a massive, gravitationally fragmenting disk and an intruder diskless star with the aim of determining the evolution of fragments that have formed in the disk prior to the encounter. Methods. Numerical hydrodynamics simulations in the non-inertial frame of reference of the host star were employed to simulate the prograde and retrograde co-planar encounters. The initial configuration of the target system (star plus disk) was obtained via a separate numerical simulation featuring the gravitational collapse of a solar-mass pre-stellar core. Results. We found that close encounters can lead to the ejection of fragments that have formed in the disk of the target prior to collision. In particular, prograde encounters are more efficient in ejecting the fragments than the retrograde encounters. The masses of ejected fragments are in the brown-dwarf mass regime. They also carry away an appreciable amount of gas in their gravitational radius of influence, implying that these objects may possess extended disks or envelopes, as also previously suggested. Close encounters can also lead to the ejection of entire spiral arms, followed by fragmentation and formation of freely-floating objects straddling the planetary mass limit. However, numerical simulations with a higher resolution are needed to confirm this finding.

We present a high angular resolution (similar to 0 ''.2), high-sensitivity (sigma similar to 0.2 mJy) survey of the 870 mu m continuum emission from the circumstellar material around 49 pre-main-sequence stars in the rho Ophiuchus molecular cloud. Because most millimeter instruments have resided in the northern hemisphere, this represents the largest high-resolution, millimeter-wave survey of the circumstellar disk content of this cloud. Our survey of 49 systems comprises 63 stars; we detect disks associated with 29 single sources, 11 binaries, 3 triple systems, and 4 transition disks. We present flux and radius distributions for these systems; in particular, this is the first presentation of a reasonably complete probability distribution of disk radii at millimeter wavelengths. We also compare the flux distribution of these protoplanetary disks with that of the disk population of the Taurus-Auriga molecular cloud. We find that disks in binaries are both significantly smaller and have much less flux than their counterparts around isolated stars. We compute truncation calculations on our binary sources and find that these disks are too small to have been affected by tidal truncation and posit some explanations for this. Lastly, our survey found three candidate gapped disks, one of which is a newly identified transition disk with no signature of a dip in infrared excess in extant observations.

Context. Whether or not a self-gravitating accretion disk fragments is still an open issue. There are many different physical and numerical explanations for fragmentation, but simulations often show a non-convergent behavior for ever better resolution. Aims. We aim to investigate the influence of different numerical limiters in Godunov type schemes on the fragmentation boundary in self-gravitating disks. Methods. We have compared the linear and non-linear outcomes in two-dimensional shearingsheet simulations using the VANLE ER and the SUPERBEE limiter. Results. We show that choosing inappropriate limiting functions to handle shock-capturing in Godunov type schemes can lead to an overestimation of the surface density in regions with shallow density gradients. The effect amplifies itself on timescales comparable to the dynamical timescale even at high resolutions. This is exactly the environment in which clumps are expected to form. The effect is present without, but scaled up by, self-gravity and also does not depend on cooling. Moreover it can be backtracked to a well known effect called oversteepening. If the effect is also observed in the linear case, the fragmentation limit is shifted to larger values of the critical cooling timescale.

Disks around brown dwarfs (BDs) are excellent laboratories to study the first steps of planet formation in cold and low-mass disk conditions. The radial-drift velocities of dust particles in BD disks higher than in disks around more massive stars. Therefore, BD disks are expected to be more depleted in millimeter-sized grains compared to disks around T Tauri or Herbig Ae/Be stars. However, recent millimeter observations of BD disks revealed low millimeter spectral indices, indicating the presence of large grains in these disks and challenging models of dust evolution. We present 3 mm photometric observations carried out with the IRAM/Plateau de Bure Interferometer (PdBI) of three BD disks in the Taurus star-forming region, which have been observed with ALMA at 0.89 mm. The disks were not resolved and only one was detected with enough confidence (similar to 3.5 sigma) with PdBI. Based on these observations, we obtain the values and lower limits of the spectral index and find low values (alpha(mm) less than or similar to 3.0). We compare these observations in the context of particle trapping by an embedded planet, a promising mechanism to explain the observational signatures in more massive and warmer disks. We find, however, that this model cannot reproduce the current millimeter observations for BD disks, and multiple-strong pressure bumps globally distributed in the disk remain as a favorable scenario to explain observations. Alternative possibilities are that the gas masses in the BD disk are very low (similar to 2 x 10(-3) M-Jup) such that the millimeter grains are decoupled and do not drift, or fast growth of fluffy aggregates.

The mass of a protoplanetary disk limits the formation and future growth of any planet. Masses of protoplanetary disks are usually calculated from measurements of the dust continuum emission by assuming an interstellar gas-to-dust ratio. To investigate the utility of CO as an alternate probe of disk mass, we use ALMA to survey (CO)-C-13 and (CO)-O-18 J = 3-2 line emission from a sample of 93 protoplanetary disks around stars and brown dwarfs with masses from 0.03 to 2 M-circle dot in the nearby Chamaeleon I star-forming region. We detect (CO)-C-13 emission from 17 sources and (CO)-O-18 from only one source. Gas masses for disks are then estimated by comparing the CO line luminosities to results from published disk models that include CO freeze-out and isotope-selective photodissociation. Under the assumption of a typical interstellar medium CO-to-H-2 ratio of 10(-4), the resulting gas masses are implausibly low, with an average gas mass of similar to 0.05M(Jup) as inferred from the average flux of stacked (CO)-C-13 lines. The low gas masses and gas-to-dust ratios for Cha I disks are both consistent with similar results from disks in the Lupus star-forming region. The faint CO line emission may instead be explained if disks have much higher gas masses, but freeze-out of CO or complex C-bearing molecules is underestimated in disk models. The conversion of CO flux to CO gas mass also suffers from uncertainties in disk structures, which could affect gas temperatures. CO emission lines will only be a good tracer of the disk mass when models for C and CO depletion are confirmed to be accurate.