Effects of disc midplane evolution on CO snowline location

Temperature changes in the planet forming disc midplanes carry important physico-chemical consequences, such as the effect on the locations of the condensation fronts of molecules - the snowlines. Snowlines impose major chemical gradients and possibly foster grain growth. The aim of this paper is to understand how disc midplane temperature changes with gas and dust evolution, and identify trends that may influence planet formation or allow to constrain disc evolution observationally. We calculate disc temperature, hydrostatic equilibrium and dust settling in a mutually consistent way from a grid of disc models at different stages of gas loss, grain growth and hole opening. We find that the CO snowline location depends very strongly on disc properties. The CO snowline location migrates closer to the star for increasing degrees of gas dispersal and dust growth. Around a typical A type star, the snowline can be anywhere between several tens and a few hundred au, depending on the disc properties such as gas mass and grain size. In fact, gas loss is as efficient as dust evolution in settling discs, and flat discs may be gas-poor counterparts of flared discs. Our results, in the context of different pre-main sequence evolution of the luminosity in low- and intermediate-mass stars suggests very different thermal (and hence chemical) histories in these two types of discs. Discs of T Tauri stars settle and cool down while discs of Herbig Ae stars may remain rather warm throughout the pre-main sequence

High angular resolution studies of protoplanetary discs

The planets, comets, asteroids... all objects in the Solar system have become from a single disc of matter at the time Sun was a young star. In the recent years it has become possible to take images of such discs elsewhere in our Galaxy, but it soon became clear that each different 'camera' we use tells a different story about the disc. In this thesis, the state-of-the-art telescopes around the globe are used to image the distribution of matter in such discs around young distant stars. We combine different pieces of information together, in a single model, and test this model against images that probe the basic disc properties: its size, mass and shape. These hi-tech snapshots of the childhood of a planetary system show that a disc may extend much further and not be as smooth as previously thought. This thesis concludes that the structure of the disc as a whole cannot be constrained without the interferometric images, and demonstrates that it is essential to combine them with the existing techniques and theory.Leiden University, Leiden Observatory, Leids Kerkhoven Boscha Fonds, Marie Curie FP6 programmeUBL - phd migration 201

Probing the radial temperature structure of protoplanetary disks with Herschel/HIFI

Herschel/HIFI spectroscopic observations of CO J=10-9, CO J=16-15 and [CII] towards HD 100546 are presented. The objective is to resolve the velocity profile of the lines to address the emitting region of the transitions and directly probe the distribution of warm gas in the disk. The spectra reveal double-peaked CO line profiles centered on the systemic velocity, consistent with a disk origin. The J=16-15 line profile is broader than that of the J=10-9 line, which in turn is broader than those of lower J transitions (6-5, 3-2, observed with APEX), thus showing a clear temperature gradient of the gas with radius. A power-law flat disk model is used to fit the CO line profiles and the CO rotational ladder simultaneously, yielding a temperature of T_0=1100 \pm 350 K (at r_0 = 13 AU) and an index of q=0.85 \pm 0.1 for the temperature radial gradient. This indicates that the gas has a steeper radial temperature gradient than the dust (mean q_{dust} ~ 0.5), providing further proof of the thermal decoupling of gas and dust at the disk heights where the CO lines form. The [CII] line profile shows a strong single-peaked profile red-shifted by 0.5 km s-1 compared to the systemic velocity. We conclude that the bulk of the [CII] emission has a non-disk origin (e.g., remnant envelope or diffuse cloud).Comment: Accepted for publication in ApJ

Resolving HD 100546 disc in the mid-infrared: Small inner disc and asymmetry near the gap

A region of roughly half of the solar system scale around the star HD 100546 is largely cleared of gas and dust, in contrast to the bright outer disc. However, some material is observed in the immediate vicinity of the star. We investigate how the dust is distributed within and outside the gap, and constrain the disc geometry with mid-infrared interferometric observations using VLTI/MIDI. With baseline lengths of 40m, our long baseline observations are sensitive to the inner few AU from the star, and we combined them with observations at shorter, 15m baselines, to probe emission beyond the gap at up to 20AU from the star. We modelled the mid-infrared emission using radial temperature profiles. Our model is composed of infinitesimal concentric annuli emitting as black bodies, and it has distinct inner and outer disc components. We derived an upper limit of 0.7AU for the radial size of the inner disc, from our longest baseline data. This small dusty disc is separated from the edge of the outer disc by a large, roughly 10AU wide gap. Our short baseline data place a bright ring of emission at 11+-1AU, consistent with prior observations of the transition region between the gap and the outer disc, known as the disc wall. The inclination and position angle are constrained by our data to i=53+-8deg and PA=145+-5deg. Compared to the rim and outer disc geometry this suggests co-planarity. Brightness asymmetry is evident in both short and long baseline data, and it is unequivocally discernible from any atmospheric or instrumental effects. The origin of the asymmetry is consistent with the bright disc wall, which we find to be 1-2AU wide. The gap is cleared of micron-sized dust, but we cannot rule out the presence of larger particles and/or perturbing bodies.Comment: 12 pages, 9 figures, accepted for publication in A&

Warm molecular gas and kinematics in the disc around HD 100546

The disc around the Herbig Ae/Be star HD 100546 is one of the most extensively studied discs in the southern sky. Although there is a wealth of information about its dust content and composition, not much is known about its gas and large scale kinematics. We detect and study the molecular gas in the disc at spatial resolution from 7.7" to 18.9" using the APEX telescope. The lines 12CO J=7-6, J=6-5, J=3-2, 13CO J=3-2 and [C I] 3P2-3P1 are observed, diagnostic of disc temperature, size, chemistry, and kinematics. We use parametric disc models that reproduce the low-J 12CO emission from Herbig~Ae stars and vary the basic disc parameters - temperature, mass and size. Using the molecular excitation and radiative transfer code RATRAN we fit the observed spectral line profiles. Our observations are consistent with more than 0.001 Msun of molecular gas in a disc of approximately 400 AU radius in Keplerian rotation around a 2.5 Msun star, seen at an inclination of 50 degrees. The detected 12CO lines are dominated by gas at 30-70~K. The non-detection of the [C I] line indicates excess ultraviolet emission above that of a B9 type model stellar atmosphere. Asymmetry in the 12CO line emission suggests that one side of the outer disc is colder by 10-20~K than the other, possibly due to a shadow by a warped geometry of the inner disc. Pointing offsets, foreground cloud absorption and asymmetry in the disc extent are excluded scenarios. Efficient heating of the outer disc ensures that low- and high-J 12CO lines are dominated by the outermost disc regions, indicating a 400 AU radius. The 12CO J=6--5 line arises from a disc layer higher above disc midplane, and warmer by 15-20~K than the layer emitting the J=3--2 line. The existing models of discs around Herbig Ae stars, assuming a B9.5 type model stellar atmosphere overproduce the [CI] 3P2--3P1 line intensity from HD 100546 by an order of magnitude.Comment: 9pages, 3figures, Accepted for publication in Astronomy & Astrophysic

TW Hya: an old protoplanetary disc revived by its planet

Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of a ∼40M⊕ Core Accretion planet is disrupted, e.g. as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way – by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes

Asymmetric mid-plane gas in ALMA images of HD 100546

In this paper we present new ALMA observations towards the proto-planet hosting transitional disc of Herbig Ae/Be star HD 100546. This includes resolved 1.3 mm continuum, $^{13}$CO and the first detection of C$^{18}$O in this disc, which displays azimuthal asymmetry in regions spatially coincident with structures previously identified in HST images related to spiral arms. The lower limit on the mass of the dust disc is calculated to be 9.6x10$^{-4}$M$_\odot$. A firm lower-limit on the total gas mass calculated from optically thin, mid-plane tracing C$^{18}$O (2-1) emission is 0.018M$_\odot$ assuming ISM abundances. These mass estimates provide an estimate of gas-to-dust ratio in the disc of 19, the ratio will increase if C$^{18}$O is relatively under-abundant in the disc compared to CO and H2. Through deprojection and azimuthal averaging of the image plane we detect 1.3 mm continuum emission out to 290+/-10 au,$^{13}$CO to 390+/-10 au and C$^{18}$O to 300+/-10au. We measure a radially increasing millimetre spectral index between wavelengths of 867$\mu$m and 1.3 mm, which shows that grain sizes increase towards the star, with solid particles growing to cm scales in the inner disc

Constraining the Gap Size in the Disk around HD 100546 in the Mid-infrared

We refine the gap size measurements of the disk surrounding the Herbig Ae star HD 100546 in the N band. Our new mid-infrared interferometric (MIDI) data have been taken with the UT baselines and span the full range of orientations. The correlated fluxes show a wavy pattern in which the minima separation links to a geometrical structure in the disk. We fit each correlated flux measurement with a spline function, deriving the corresponding spatial scale, while assuming that the pattern arises interferometrically due to the bright emission from the inner disk and the opposing sides of the wall of the outer disk. We then fit an ellipse to the derived separations at their corresponding position angles, thereby using the observations to constrain the disk inclination to i = 47° ± 1° and the disk position angle to PA = 135⁰0 ± 2⁰5 east of north, both of which are consistent with the estimated values in previous studies. We also derive the radius of the ellipse to 15.7 ± 0.8 au. To confirm that the minima separations translate to a geometrical structure in the disk, we model the disk of HD 100546 using a semianalytical approach taking into account the temperature and optical depth gradients. Using this model, we simultaneously reproduce the level and the minima of the correlated fluxes and constrain the gap size of the disk for each observation. The values obtained for the projected gap size in different orientations are consistent with the separation found by the geometrical model

Detection of the Water Reservoir in a Forming Planetary System

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extra-solar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on-board the Herschel Space Observatory has detected emission from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface hinting at a water ice reservoir equivalent to several thousand Earth Oceans in mass. The water's ortho-to-para ratio falls well below that of Solar System comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.Comment: 18 pages, 2 figures. Corrected typo in reported mass (in g) of detected water vapor reservoir. All conclusions are unchange