Evidence for CO depletion in the inner regions of gas-rich protoplanetary disks

G. van der Plas; M. E. van den Ancker; B. Acke; A. Carmona; C. Dominik; D. Fedele; L. B. F. M. Waters

doi:10.1051/0004-6361/200811148

Free Access

Issue		A&A Volume 500, Number 3, June IV 2009


Page(s)		1137 - 1141
Section		Interstellar and circumstellar matter
DOI		https://doi.org/10.1051/0004-6361/200811148
Published online		29 April 2009

Evidence for CO depletion in the inner regions of gas-rich protoplanetary disks

G. van der Plas^1,2 - M. E. van den Ancker¹ - B. Acke^3, - A. Carmona^4,5 - C. Dominik^2,6 - D. Fedele⁷ - L. B. F. M. Waters^2,3

1 - European Southern Observatory, Karl-Schwarzschild-Str.2, 85748 Garching bei München, Germany
2 - Astronomical Institute Anton Pannekoek, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands
3 - Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
4 - ISDC Data Centre for Astrophysics, Ch. d'Ecogia 16, 1290 Versoix, Switzerland
5 - Observatoire de Genève, Ch. des Maillettes 51, 1290 Sauverny, Switzerland
6 - Institute for Astrophysics, Radbout University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
7 - Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany

Received 13 October 2008 / Accepted 17 April 2009

Abstract
Aims. We investigate the physical properties and spatial distribution of Carbon Monoxide (CO) gas in the disks around the Herbig Ae/Be stars HD 97048 and HD 100546.
Methods. Using high-spectral-resolution 4.588-4.715 $\mu$ m spectra containing fundamental CO emission taken with CRIRES on the VLT, we probe the circumstellar gas and model the kinematics of the emission lines. By using spectro-astrometry on the spatially resolved targets, we constrain the physical size of the emitting regions in the disks.
Results. We resolve, spectrally and spatially, the emission of the ¹³CO v(1-0) vibrational band and the ¹²CO v=1-0, v=2-1, v=3-2 and v=4-3 vibrational bands in both targets, as well as the ¹²CO v=5-4 band in HD 100546. Modeling of the CO emission with a homogeneous disk in Keplerian motion, yields a best fit with an inner and outer radius of the CO emitting region of 11 and $\geq$ 100 AU for HD 97048. HD 100546 is not fit well with our model, but we derive a lower limit on the inner radius of 8 AU. The fact that gaseous [OI] emission was previously detected in both targets at significantly smaller radii suggests that CO may be effectively destroyed at small radii in the surface layers of these disks.

Key words: circumstellar matter - stars: pre-main sequence - stars: individual: HD 97048 - stars: individual: HD 100546 - planetary systems: protoplanetary disks

1 Introduction

We present a study of warm carbon monoxide (CO) gas in two circumstellar disks surrounding young intermediate-mass stars. CO is commonly detected in disks (Najita et al. 2003), and often used as a tracer of the physical conditions in the disk and the kinematics of the gas. Ro-vibrational transitions of warm (T = 120-1010 K) CO gas have already been found by Mitchell et al. (1990) in 8 out of 9 surveyed embedded infrared sources. Brittain et al. (2007) had a 100 % detection rate for the ro-vibrational CO transitions in the disks around 9 intermediate mass stars with optically thick inner disks, while their detection rate is only 1 out of 5 for disks with an optically thin inner disk, a case they explain by UV fluorescence. Blake & Boogert (2004) have detected CO gas in 5/5 Herbig Ae stars and recently, Pontoppidan et al. (2008) have spatially and spectrally resolved the 4.7 $\mu$ m CO ro vibrational lines in the disks around 3 young solar mass stars with known dust gaps or inner holes, and detected CO well inside the dust gaps in all disks.

$\begin{figure} \par\includegraphics[height=6.7cm,width=17.7cm]{1148fig1.eps} \end{figure}$

Figure 1:

Reduced CRIRES spectra of HD 97048 ( top) and HD 100546 ( bottom). The black spectrum has been used for the line measurements, whereas the red line shows the areas that have been ignored due to large errors, mostly caused by the correction for saturated telluric lines. The line identifications plotted above show detections of the isotope (¹²CO and ¹³CO), the vibrational transition (e.g. v=1-0), and the rotational transition (R(xx) or P(xx)). The broad hydrogen recombination lines Pfund $\beta$ at 4654 nm and Humphreys $\epsilon$ at 4673 nm are also shown.

Open with DEXTER

HD 97048 (Sp. type A0pshe, van den Ancker et al. 1998) is a well-studied nearby Herbig Ae star, whose spectral energy distribution (SED) displays a large infrared excess above the emission from the star, commonly explained by emission from warm dust in a circumstellar disk. The disk is a flaring disk according to the classification proposed by Meeus et al. (2001). Lagage et al. (2006) resolved the disk around HD 97048 in the mid-IR and determined the flaring angle. PAH emission from radii up to 200-300 AU has been detected by van Boekel et al. (2004). Habart et al. (2004) have resolved emission from the 3.43 $\mu$ m and 3.54 $\mu$ m diamond bands from within the inner 15 AU in the disk. Circumstellar gas has also been detected, e.g. by Acke & van den Ancker (2006) who detect spatially and spectrally resolved [OI] emission from a few tenths to several tens of AU. HD 97048 is one of the only two Herbig stars where the H₂ S(1) line at 17.035 $\mu$ m has been detected (Martin-Zaïdi et al. 2007), coming from warm ( $T_{\rm ex}$ $\leq$ 570 K) gas within 35 AU.

HD 100546 (Sp. type B9Vne, van den Ancker et al. 1998) also has a SED characteristic for a flared disk. Its IR spectrum indicates a high crystalline silicate dust fraction in the disk surface layers (Bouwman et al. 2003). There is observational evidence for a disk gap, or ``second wall'', at $\approx$ 10 AU from the central star. This was derived based on the SED (Bouwman et al. 2003), and confirmed by nulling interferometry (Liu et al. 2007). The [OI] 6300 Å line was studied by Acke & van den Ancker (2006), who also find evidence for a disk gap at that distance, possibly induced by a very low mass stellar companion or planet of ${\approx}20~M_{\rm jup}$ , or a planetary system. The latter authors spatially and spectrally resolve the [OI] emission, and determine it to originate from a few tenths to several tens of AU. Geers et al. (2007) have detected PAH emission from 12 $\pm$ 3 AU, and the circumstellar material has been spatially resolved many times, revealing rich structures in the disk.

In this paper, we present high spectral resolution observations with milli-arcsecond spatial precision of CO emission originating from the disks around HD 97048 and HD 100546.

2 Observations and data reduction

High-spectral-resolution (R $\approx$ 94 000, determined from telluric lines) spectra of HD 97048 and HD 100546 were obtained on the early morning of June 16 ${\rm th}$ 2007 with the VLT cryogenic high-resolution infrared echelle spectrograph (CRIRES, Käufl et al. (2004)). Adaptive Optics were used to optimize the signal-to-noise and the spatial resolution of the observations. These cover the wavelength range of 4.588-4.711 $\mu$ m, and contain many fundamental ro-vibrational transitions of the CO molecule. The observations were made with a slit width of 0.2'', and with the slit rotated along the parallactic angle. HD 97048 and HD 100546 were observed with 2 wavelength settings, each for 10 and 3 min respectively, and the slit Position Angle (PA) during the observations spans a range of 41.1-50.3 degrees for HD 97048 and 74.8-78.2 degrees for HD 100546. This spread in PA, close to the minor axis of the disks, is not optimal for our astrometry purpose, but was enforced by the observational circumstances. The data were reduced using the CRIRES pipeline V1.5.0, which performs wavelength calibration, background subtraction and flatfield correction.

To correct for telluric absorption we have observed a telluric standard star directly after each science target: HIP 052419 (Sp. type B0V) for HD 97048 and HIP 061585 (Sp. type B2IV) for HD 100546. Each standard was chosen to be as close as possible on the sky so that their spectra are affected by similar atmospheric conditions. We compared their spectra to an appropriate Kurucz (1991) stellar atmosphere model to determine the instrumental response. The optical depth of the telluric lines of the standard was scaled to that of the science target. The science spectra were derived by division with the instrumental response. A sub-pixel wavelength shift was manually added to reduce telluric residuals (i.e. spikes due to minor wavelength mismatches). Some telluric absorption lines are fully saturated, causing problems with the division of the spectra (i.e. division by 0). In the further data reduction, we will disregard these areas. We have shifted the spectrum to center the v=2-1 lines on their rest wavelength. The spectra are shown in Fig. 1.

3 Observational analysis

To check whether the targets are spatially resolved, we perform spectro-astrometry on the reduced 2D frames. We determine the location of the photocenter of the PSF and the spatial FWHM as function of wavelength (or equivalently velocity for each line), by fitting a Gaussian to the spatial profile in each spectral column. The position of the photocenter relative to the continuum, the Spatial Peak Position (SPP), is a measure for the spatial offset of the emitting region compared to the continuum. Our observations are made with the slit positioned over a small range in position angles. We note that Pontoppidan et al. (2008) have observed three nearby T Tauri stars at 6 different PA's with CRIRES and the same observational set-up as employed by us, and find, after pairwise subtracting parallel and anti parallel spectro-astrometric signals, that no telluric feature produces false spectro-astrometric signals in excess of 0.2-0.5 mas. Since we do not possess anti-parallel spectra, we show the telluric absorption in Figs. 2 and 3, and note, by comparing with CO lines un-affected by telluric absorption, that their effect on the astrometric signal is small. We compare the SPP in the CO lines to the SPP of the underlying continuum (SPP = 0 by definition). If the target is resolved, and the photocenter of the CO emission region is offset to that of the dust continuum emission from the disk, the SPP is expected to be non-zero. The photocenter displacement can be derived from the measured SPP if we correct for dilution by the continuum flux, by multiplying the SPP with a factor $1 + F_{\rm c}(v)/F_l(v)$ , with $F_{\rm c}(v)$ and $F_{\rm l}(v)$ the continuum and line flux at velocity v (Pontoppidan et al. 2008).

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig2.eps} \end{figure}$

Figure 2:

In all three windows the thick black line represents the spectrum of HD 97048 (orange thin line) smoothed over 2 pixels, and the red line denotes the same data for the telluric standard. Shown with the vertical dotted lines are the locations of the ( from left to right) ¹²CO v=4-3 R(23), v=3-2 R(14), v=2-1 R(06), v=5-4 R(33) and ¹³CO v=1-0 R(12) emission lines. Top: close up of the spectrum of HD 97048 (as shown in Fig. 1). Middle: the SPP of HD 97048 and of the telluric standard, shifted down by 10 mas. Bottom: the spatial FWHM of the PSF for each spectral row for HD 97048 and the telluric standard.

Open with DEXTER

The Gaussian fit to the PSF also yields a spatial FWHM. If the continuum is unresolved, this quantity is a measure for the achieved effective spatial resolution in the direction parallel to the slit. When resolved, the FWHM allows us to determine the spatial extent of the emission. Outside the spectral lines it then is a direct probe for the size of the continuum emitting disk, whilst inside the lines it is a function of both the line- and continuum emission, and its value a lower limit to the CO emission. Thus, measuring the SPP and the FWHM of the PSF in the continuum and the CO lines allows us to determine [1] the size of the continuum emitting region; and [2] the spatial offset of each CO line with respect to the continuum. This procedure is shown in Figs. 2 and 3.

4 Results

We detect in the spectrum of HD 97048 the ¹³CO v=1-0, as well as in the ¹²CO v=1-0 to v=4-3 vibrational bands. All lines are spectrally resolved with a remarkably constant FWHM of 15 $\rm km~s^{-1}$ . The Half Width at Zero Intensity (HWZI) is more difficult to estimate due to the continuum S/N, but we make a conservative estimate of 15 $\rm km~s^{-1}$ corresponding to an inner radius of 4.5 AU, assuming the stellar parameters and inclination as given in Table 1. Average line profiles and astrometric signals of the ¹²CO v=1-0 to v=3-2 are shown in Fig. 5 and Table 3. The 4.6 $\mu$ m continuum emission of the disk around HD 97048 is resolved. Its FWHM of 0.231'' $\pm$ 0.004'' corresponds to 23 $\pm$ 3 AU, after correcting for the FWHM of the PSF. The latter (0.192'' $\pm$ 0.001'') was estimated from 10 unresolved telluric standards, all sufficiently bright for the AO loop to close at maximum frequency, observed with the same instrumental settings, and within 2 days of our observations. The red- and blue- shifted wings of the SPP form a sinusoidal signal, typical for an inclined circumstellar disk.

In the spectrum of HD 100546 we detect the ¹³CO v=1-0, as well as in the ¹²CO v=1-0 to v=5-4 vibrational bands. All lines are spectrally resolved and have a FWHM of 15 $\rm km~s^{-1}$ , irrespective of the transition. A conservative estimate for the HWZI is 17.5 $\rm km~s^{-1}$ corresponding to an inner radius of 3.3 AU. The average line profiles and astrometric signals of the ¹²CO v=1-0 to v=4-3 transitions are shown in Fig. 7 and their values in Table 3. The 4.6 $\mu$ m continuum emission of the disk around HD 100546 is resolved. Its FWHM of 0.218'' $\pm$ 0.002'' corresponds to 11 $\pm$ 1 AU, after correcting for the FWHM of the PSF. The SPP is very asymmetric, i.e the negative shift of the blue wing in the SPP is larger than the positive SPP of the red wing. This asymmetry is also seen in the FWHM and is more pronounced for the lower vibrational transitions.

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig3.eps} \end{figure}$	Figure 3: Same as Fig. 2 for HD 100546.
Open with DEXTER

Table 1: Astrophysical parameters of the programme stars.

Table 2: Log of spectroscopic observations at June $16{\rm th}$ 2007.

Table 3: Spatial information for the averaged vibrational transitions shown in Figs. 5 and 7, and the radii used to model the kinematics.

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig4.eps} \end{figure}$	Figure 4: Projected images of the simulated CO disk with the best fit parameters for HD 97048 and HD 100546, together with the slit of the spectrograph at the time of observation. North is up, East is left.
Open with DEXTER

5 A Simple model

A powerful method to determine the extent of the emitting gas is to model the kinematics of the line- and astrometric profiles. We construct the average profiles for the v=1-0 to v=4-3 vibrational bands in Figs. 5 and 7, and fit these to the model shown in Fig. 4. In this model, the gas is in Keplerian orbit in a flat disk with known inclination and PA around a star with known stellar mass and distance. Further, the intensity of the emission decreases as $I(R) = I_{\rm in}\left(\frac{R}{R_{\rm in}}\right)^{-\alpha}$ , with $I_{\rm in}$ the intensity at the inner radius $R_{\rm in}$ , and R the radial distance from the star (Carmona et al. 2007). The astrometric signal is consequently determined by fitting a Gaussian to the simulated PSF. Because the stellar parameters are all well constrained, the free parameters are $\alpha$ , $R_{\rm in}$ , and $R_{\rm out}$ . $\alpha$ describes the change in intensity of the CO emission as function of radius, where $\alpha$ = 2 reflects an intensity directly proportional to the stellar radiation field assuming the surface density has no radial dependence. We compare the line profile, the SPP and the FWHM with their simulated counterparts for the parameter-space $1~{\rm AU} \leq R_{\rm in} \leq 29~{\rm AU}$ , $25~{\rm AU} \leq R_{\rm out} \leq 200~{\rm AU}$ and $0 \leq \alpha \leq 5.5$ , and quantify the best fit by minimizing the reduced $\chi^2$ . We estimate the variance on the data from the continuum adjacent to the line, and determine the best fit from the v=2-1 line profile, because the v=1-0 data is affected by telluric absorption, and because this transitions has a higher S/N compared to the higher vibrational transitions. Since the continuum corrected SPP is similar for each vibrational transition within the separate stars, we assume that all vibrational transitions have the same spatial origin. For HD 100546, and in lesser extent HD 97048, our model produces a more distinctly double peaked line profile than the observed profiles. This phenomenon is also observed for 9 out of 60 protoplanetary disks (Bast et al., in prep.), and suggests there is an extra low projected velocity gas component present in those disks. We have repeated the fitting procedure for only the part of the line profile sensitive to the inner part of the disk, the line wings ( $\left\vert V_{\rm proj}\right\vert \geq 4.5~{\rm km~s}^{-1}$ ), and find no difference for the best fit.

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig5.eps} \end{figure}$

Figure 5:

Top: the average line profiles of the HD 97048 v=1-0, v=2-1 and v=3-2 transitions, created with 5 (5 < J < 7), 12 ( 5 < J < 14) and 7 ( 14 < J < 21) lines respectively, are shown with the black histogram. The blue wing of the v=1-0 transition is contaminated by telluric CO lines, shown with dots. Over plotted is the 6300 Å [OI] emission (red dotted line), and the model fit to the CO emission with the solid red line. Middle: the corresponding, continuum corrected (Sect. 3, the uncorrected SPP is shown with the dotted line), average astrometric signal. The model is over plotted in red. Bottom: the average spatial FWHM in black, with the model over plotted in red.

Open with DEXTER

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig6.eps} \end{figure}$

Figure 6:

$\chi ^2_{\rm red}$ surfaces for the best fitting disk model for, from left to right, the line profile, SPP and FWHM of the v=2-1 transition of HD 97048. The minimal $\chi ^2_{\rm red}$ for each fit is given in the panels, and the 1-3 $\sigma$ lines represent $\geq$ 68 $\%$ , $\geq$ 90 $\%$ and $\geq$ 99 $\%$ confidence levels. The best fit parameters for each panel are shown by a red circle. To better show the $\chi ^2_{\rm red}$ surface landscape we have overplotted the $\chi ^2_{\rm red}$ = 5, 10 ( left two panels, $\chi ^2_{\rm red}$ = 1.5, 1.7 in the right panel) contours.

Open with DEXTER

The best fit to the v=1-0, v=2-1, and v=3-2 vibrational transitions of HD 97048, are $R_{\rm in}$ = 11 AU, $R_{\rm out}$ $\geq$ 100 AU and $\alpha$ = 3.00. We show the $\chi ^2_{\rm red}$ surfaces for the v=2-1 fit to the line profile, the SPP and the FWHM in Fig. 6. These values have been derived using a PA of 180 $^{\circ}$ , estimated from direct imaging of the circumstellar disk (Fig. S1 of Lagage et al. 2006), compared to 160 $^{\pm 19 \circ}$ , found by Acke & van den Ancker (2006) using astrometry. Using the latter value makes the simulated line profile slightly more double peaked, and thus increases the minimal $\chi ^2_{\rm red}$ , but does not change the best fit parameters.

The spatial distribution of gas in HD 100546 is more complex than the ``simple'' homogeneous disk in Keplerian motion we model, and the best fit values for each transition vary, with a lowest value of $R_{\rm in}$ $\geq$ 4 AU and $\alpha$ $\geq$ 1.5. If we assume that all vibrational states have the same spatial distribution and origin, and we only fit the line wings, values of $R_{\rm in}$ = 8 AU and $\alpha$ = 1.5 are found. The poor quality of the fit is possibly due to complex disk geometry and gas kinematics, related to the presence of a disk gap and consequent ``second wall'' at 10 AU (Bouwman et al. 2003). The average profiles of all vibrational bands have a FWHM of 15 $\rm km~s^{-1}$ , hinting that they share a common spatial origin. The [OI] emission has much wider wings (FWHM = 24 and 25 $\rm km~s^{-1}$ for HD 97048 and HD 100546), corresponding to a smaller $R_{\rm in}$ (80% of the [OI] emission emanates from the region between 0.8 and 20 AU, Fig. 7 in Acke & van den Ancker 2006).

$\begin{figure} \par\includegraphics[width=8.9cm,clip]{1148fig7.eps} \end{figure}$	Figure 7: Same as Fig. 5 for HD 100546, but the v=1-0, v=2-1, v=3-2 and v=4-3 profiles are made from 11 (0 < J < 7), 12 ( 5 < J < 14), 9 ( 11 < J < 26) and 12 ( 19 < J < 34) lines respectively.
Open with DEXTER

6 Discussion

Fundamental CO emission from Herbig stars, and gaps therein, have been detected in more circumstellar disks, e.g. in SR 21 with an $R_{\rm in}$ of 7.6 $\pm$ 0.4 AU for the ¹²CO v(1-0) line (Pontoppidan et al. 2008), while there is an inner dust component at $\leq$ 0.5 AU, followed by a dust gap up to 18 AU (Brown et al. 2007). Goto et al. (2006) find an inner clearing in the CO emission at 11 $\pm$ 2 AU around HD 141569, and suggest the cavity most likely is cleared by photo-evaporation working together with viscous accretion. Brittain et al. (2007) observe the same star, find CO emission starting between 6 and 15 AU, and explain the emission as due to UV fluorescence. Blake & Boogert (2004) study a sample of 5 Herbig stars, not known to have disk gaps or inner holes, and find that the CO gas traces both collisional driven emission in the inner (<0.5-1 AU) disk, as well as CO gas excited by resonance fluorescence (R < 50-100 AU). These observations suggest that CO depletion in the inner disk is associated with large disk gaps or inner holes. Finally, in all cases discussed above (except SR 21 for which no data exists), the [OI] emission, when detected, arises at similar or higher velocities than the CO fundamental ro-vibrational emission.

Contrary to the above, the lack of CO emission at small radii in HD 97048 and HD 100546 can not be caused by a gas-free inner disk. [OI] 6300 Å emission, a tracer of OH molecular gas (Acke et al. 2005) in the surface layers of flared disks, and induced by the stellar UV field, has been observed for both stars by Acke & van den Ancker (2006), with $80 \%$ of the emission originating between 0.8 < R < 20 AU. The SEDs of our targets also display large NIR excesses, typical for hot dust close to the star (Acke & van den Ancker 2004; Bouwman et al. 2003). In the cases of HD 97048 and HD 100546, the absence of CO emission at Small radii can thus not be attributed to the paucity of matter in the inner disk, nor to the absence of UV flux. Disk holes and/or large gaps therefore are not the only cause of CO depletion in the inner disk. Here, we suggest that CO is efficiently destroyed in the inner disk. A possible mechanism is the photo dissociation of CO and the subsequent formation of simple organic species such as C₂H₂, HCN and CH₄ (Agúndez et al. 2008). Further analysis (e.g. observations of hot transitions of carbon-bearing molecules) is needed to determine the reason why the CO gas is not detected in the inner regions, and may offer us a glimpse of the molecular chemistry of these inner disks.

The similar line widths and shapes for all vibrational states in HD 97048, as well as the similar continuum corrected SPP, suggest a common excitation mechanism and location. The best fit parameter for the radial dependency of the CO emission intensity $\alpha$ = 3.0 can guide models that aim to explain the highly excited vibrational line emission, such as UV fluorescence. The similar line shape, continuum corrected SPP and FWHM of the higher (v=3-2 and v=4-3) vibrational transitions in HD 100546 also suggest a common excitation mechanism and location, but the v=2-1 and especially v=1-0 vibrational transitions deviate significantly. The asymmetry in the data on HD 100546, and the poor fit to the model, point to a more complicated disk geometry, in agreement with finding of Bouwman et al. (2003) and Acke & van den Ancker (2006).

Acknowledgements

The authors of this paper wish to thank the Paranal nighttime astronomer E. Valenti for her cheerful assistance making the observations, and also the anonymous referee for helping improving this paper.

References

Acke, B., & van den Ancker, M. E. 2004, A&A, 426, 151 [NASA ADS] [CrossRef] [EDP Sciences]
Acke, B., & van den Ancker, M. E. 2006, A&A, 449, 267 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Acke, B., van den Ancker, M. E., Dullemond, C. P., van Boekel, R., & Waters, L. B. F. M. 2004, A&A, 422, 621 [NASA ADS] [CrossRef] [EDP Sciences]
Acke, B., van den Ancker, M. E., & Dullemond, C. P. 2005, A&A, 436, 209 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Ardila, D. R., Golimowski, D. A., Krist, J. E., et al. 2007, ApJ, 665, 512 [NASA ADS] [CrossRef]
Agúndez, M., Cernicharo, J., & Goicoechea, J. R. 2008, A&A, 483, 831 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Boogert, A. C. A., Hogerheijde, M. R., & Blake, G. A. 2002, ApJ, 568, 761 [NASA ADS] [CrossRef]
Blake, G. A., & Boogert, A. C. A. 2004, ApJ, 606, L73 [NASA ADS] [CrossRef] (In the text)
Bouwman, J., de Koter, A., Dominik, C., & Waters, L. B. F. M. 2003, A&A, 401, 577 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Brannigan, E., Takami, M., Chrysostomou, A., & Bailey, J. 2006, MNRAS, 367, 315 [NASA ADS] [CrossRef]
Brittain, S. D., Simon, T., Najita, J. R., & Rettig, T. W. 2007, ApJ, 659, 685 [NASA ADS] [CrossRef] (In the text)
Brown, J. M., Blake, G. A., Dullemond, C. P., et al. 2007, ApJ, 664, L107 [NASA ADS] [CrossRef] (In the text)
Carmona, A., van den Ancker, M. E., Henning, T., et al. 2007, A&A, 476, 853 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Dullemond, C. P., & Dominik, C. 2004, A&A, 417, 159 [NASA ADS] [CrossRef] [EDP Sciences]
Geers, V. C., van Dishoeck, E. F., Visser, R., et al. 2007, A&A, 476, 279 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Goto, M., Usuda, T., Dullemond, C. P., et al. 2006, ApJ, 652, 758 [NASA ADS] [CrossRef] (In the text)
Habart, E., Testi, L., Natta, A., & Carbillet, M. 2004, ApJ, 614, L129 [NASA ADS] [CrossRef] (In the text)
Käufl, H.-U., Ballester, P., Biereichel, P., et al. 2004, Proc. SPIE, 5492, 1218 [NASA ADS] (In the text)
Kurucz, R. L. 1991, Precision Photometry: Astrophysics of the Galaxy, 27 (In the text)
Lagage, P.-O., Doucet, C., Pantin, E., et al. 2006, Science, 314, 621 [NASA ADS] [CrossRef] (In the text)
Liu, W. M., Hinz, P. M., Meyer, M. R., et al. 2007, ApJ, 658, 1164 [NASA ADS] [CrossRef] (In the text)
Martin-Zaïdi, C., Lagage, P.-O., Pantin, E., & Habart, E. 2007, ApJ, 666, L117 [NASA ADS] [CrossRef] (In the text)
Martin-Zaidi, C., Habart, E., Augereau, J.-C., et al. 2009 [arXiv:0902.1614]
Meeus, G., Waters, L. B. F. M., Bouwman, J., et al. 2001, A&A, 365, 476 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
Mitchell, G. F., Maillard, J.-P., Allen, M., Beer, R., & Belcourt, K. 1990, ApJ, 363, 554 [NASA ADS] [CrossRef] (In the text)
Najita, J., Carr, J. S., & Mathieu, R. D. 2003, ApJ, 589, 931 [NASA ADS] [CrossRef] (In the text)
Pontoppidan, K. M., Blake, G. A., van Dishoeck, E. F., et al. 2008, ApJ, 684, 1323 [NASA ADS] [CrossRef] (In the text)
Salyk, C., Blake, G. A., Boogert, A. C. A., & Brown, J. M. 2007, ApJ, 655, L105 [NASA ADS] [CrossRef]
van Boekel, R., Waters, L. B. F. M., Dominik, C., et al. 2004, A&A, 418, 177 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 1998, A&A, 330, 145 [NASA ADS] (In the text)

Footnotes

... disks: Based on observations collected at the European Southern Observatory, Paranal, Chile. (Program ID 079.C-0349A).
...: Postdoctoral Fellow of the Fund for Scientific Research, Flanders.
... (CRIRES: http://www.eso.org/sci/facilities/paranal/instruments/crires/
... pipeline V1.5.0: http://www.eso.org/sci/data-processing/software/pipelines/index.html

All Tables

Table 1: Astrophysical parameters of the programme stars.

Table 2: Log of spectroscopic observations at June $16{\rm th}$ 2007.

Table 3: Spatial information for the averaged vibrational transitions shown in Figs. 5 and 7, and the radii used to model the kinematics.

All Figures

$\begin{figure} \par\includegraphics[height=6.7cm,width=17.7cm]{1148fig1.eps} \end{figure}$	Figure 1: Reduced CRIRES spectra of HD 97048 ( top) and HD 100546 ( bottom). The black spectrum has been used for the line measurements, whereas the red line shows the areas that have been ignored due to large errors, mostly caused by the correction for saturated telluric lines. The line identifications plotted above show detections of the isotope (¹²CO and ¹³CO), the vibrational transition (e.g. v=1-0), and the rotational transition (R(xx) or P(xx)). The broad hydrogen recombination lines Pfund $\beta$ at 4654 nm and Humphreys $\epsilon$ at 4673 nm are also shown.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig2.eps} \end{figure}$	Figure 2: In all three windows the thick black line represents the spectrum of HD 97048 (orange thin line) smoothed over 2 pixels, and the red line denotes the same data for the telluric standard. Shown with the vertical dotted lines are the locations of the ( from left to right) ¹²CO v=4-3 R(23), v=3-2 R(14), v=2-1 R(06), v=5-4 R(33) and ¹³CO v=1-0 R(12) emission lines. Top: close up of the spectrum of HD 97048 (as shown in Fig. 1). Middle: the SPP of HD 97048 and of the telluric standard, shifted down by 10 mas. Bottom: the spatial FWHM of the PSF for each spectral row for HD 97048 and the telluric standard.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig3.eps} \end{figure}$	Figure 3: Same as Fig. 2 for HD 100546.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig4.eps} \end{figure}$	Figure 4: Projected images of the simulated CO disk with the best fit parameters for HD 97048 and HD 100546, together with the slit of the spectrograph at the time of observation. North is up, East is left.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig5.eps} \end{figure}$	Figure 5: Top: the average line profiles of the HD 97048 v=1-0, v=2-1 and v=3-2 transitions, created with 5 (5 < J < 7), 12 ( 5 < J < 14) and 7 ( 14 < J < 21) lines respectively, are shown with the black histogram. The blue wing of the v=1-0 transition is contaminated by telluric CO lines, shown with dots. Over plotted is the 6300 Å [OI] emission (red dotted line), and the model fit to the CO emission with the solid red line. Middle: the corresponding, continuum corrected (Sect. 3, the uncorrected SPP is shown with the dotted line), average astrometric signal. The model is over plotted in red. Bottom: the average spatial FWHM in black, with the model over plotted in red.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=9cm,clip]{1148fig6.eps} \end{figure}$	Figure 6: $\chi ^2_{\rm red}$ surfaces for the best fitting disk model for, from left to right, the line profile, SPP and FWHM of the v=2-1 transition of HD 97048. The minimal $\chi ^2_{\rm red}$ for each fit is given in the panels, and the 1-3 $\sigma$ lines represent $\geq$ 68 $\%$ , $\geq$ 90 $\%$ and $\geq$ 99 $\%$ confidence levels. The best fit parameters for each panel are shown by a red circle. To better show the $\chi ^2_{\rm red}$ surface landscape we have overplotted the $\chi ^2_{\rm red}$ = 5, 10 ( left two panels, $\chi ^2_{\rm red}$ = 1.5, 1.7 in the right panel) contours.
Open with DEXTER
In the text

$\begin{figure} \par\includegraphics[width=8.9cm,clip]{1148fig7.eps} \end{figure}$	Figure 7: Same as Fig. 5 for HD 100546, but the v=1-0, v=2-1, v=3-2 and v=4-3 profiles are made from 11 (0 < J < 7), 12 ( 5 < J < 14), 9 ( 11 < J < 26) and 12 ( 19 < J < 34) lines respectively.
Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.