A&A 444, L57-L60 (2005)
DOI: 10.1051/0004-6361:200500211
E. Dartois 1 - K. Pontoppidan 2 - W.-F. Thi 3 - G. M. Muñoz Caro 4
1 - Institut d'Astrophysique Spatiale, UMR 8617, Université Paris-Sud, bâtiment 121,
91405 Orsay, France
2 -
Leiden Observatory, PO Box 9513, 2300 RA Leiden, The Netherlands
3 -
ESTEC, ESA, Research Support Science Department, PO Box 2200AZ, Noordwijk, The Netherlands
4 -
Centro de Astrobiologia, Spain
Received 15 August 2005 / Accepted 3 November 2005
Abstract
The L723 low mass class 0 protostellar object has been observed
with the IRS low resolution spectrometer onboard Spitzer Space
Telescope. A huge CO2 ice column density of
6.3
1
1018 cm-2 is detected toward this
extremely red object. We constructed the spectral energy
distribution (SED) of L723 by combining the Spitzer spectral data
with Infrared Space Observatory (ISO) camera extracted fluxes, IRAS points, ISO long wavelength spectrometer (LWS) observations and millimeter ground based observations. A self consistent 1D model
was developed to fit the SED. We show that, besides the high CO2 column density, the abundance ratio of CO2 ice to silicates
grain cores is similar to other lines-of-sight.
Key words: astrochemistry - line: identification - ISM: dust, extinction - ISM: molecules - ISM: lines and bands - ISM: individual: L723
L723 is an isolated dark cloud located at a distance of about 300 pc
(Goldsmith et al. 1984). It displays a quadrupolar outflow of several
arcminutes on the sky, associated with the class 0 source
IRAS
19156+1906 (Lee et al. 2002, and reference therein). VLA observations revealed a double source (Anglada et al. 1991) named VLA1
and VLA2. The second one is associated with the observed dust
millimeter emission (Cabrit & Andre 1991). The bolometric luminosity
associated with the object is estimated between about 1.9
(Reipurth et al. 1993) to 3.3
(Shirley et al. 2000). The highest extinction is confirmed to be
associated to VLA2 by CS observations delineating the dense core
condensation (Hirano et al. 1998) of about 0.04 pc extension
(
104 AU), with about one solar mass in the circumstellar
dust (Estalella et al. 2003). In this paper we present spectroscopic
data in the 14-27
m wavelength range obtained with the Spitzer
telescope and showing the presence of large amounts of CO2 ice. These observations are described in Sect. 2 together with the existing data to build the spectral energy distribution
(SED) of this source. We then present the radiative transfer model
used to describe the SED in Sect. 3. In the last section we discuss
the results and compare the model adopted CO2 abundance with
some of the previous determinations (e.g. de Graauw et al. 1996;
Gerakines et al. 1999; Watson et al. 2004; Bergin et al. 2005).
![]() |
Figure 1:
Left: Spitzer IRS observations of IRAS
19156+1906 (L723-VLA2) in the 14-27 ![]() ![]() |
Open with DEXTER |
Low-resolution spectra (
70-120) were acquired with the
Infrared Spectrograph (IRS, Houck et al. 2004) onboard the Spitzer
Space Telescope (Werner et al. 2004) in the Long-Low wavelength
acquisition mode, covering the 14-27
m range. The integration
times of the ramps were set to 30 s. Data were reduced at the PBCD
level with in-house software extraction using the pipeline
calibration files. The distortion of the spectrum on the sky was
corrected by fitting a second order polynomial to the trace of the
PSF on the wavelength axis, to follow its evolution. The spectrum
was extracted with a Gaussian PSF profile weighting scheme centered
on this estimated trace. An estimation of the background emission
on each side of the point source was subtracted to the spectrum to
limit residual contamination. The absolute flux calibration was
estimated using the electron-to-Jy conversion polynomial given in
the appropriate Spitzer calibration file. The corresponding spectra
are displayed in Fig. 1.
Additional observations from the literature or database archive
were added to produce a spectral energy distribution (SED)
overview. The 1.3 mm flux was obtained from Motte & André (2001),
the 450 and 850
m observations from Shirley et al. (2000). In
addition to the IRAS fluxes at 25, 60 and 100
m (the 12
m
is only an upper limit), we added the Infrared Space Observatory
LWS spectrometer observations extracted from the ISO archive
and flux corrected
with the extended source flux correction given in the ISO-LWS handbook, Sect. 5.9.3. The data of each individual band were rebinned
to produce a single point, without additional gain correction. L723
has also been mapped with ISOCAM in the LW6 filter
(7.0-8.5
m). Data were retrieved from the ISO database and
reprocessed with an in house software. These additional fluxes are
summarized in Table 1.
Table 1: Summary of observations.
A 1 D ray-tracing iterative self consistent code based on the
moment method has been built to calculate the emergent spectrum of
the deeply embedded object. The iterative scheme and temperature
correction follows the ones used in e.g. Efstathiou &
Rowan-Robinson (1990). The radial density is defined by a power law
,
where p is the
index of the law, r the radial point considered (AU), and n100 is the H2 density at 100 astronomical units (AU). The dust mass equals 1/100th of the H2 one. The outer
radius
is a fixed input, whereas the inner one (
)
is self consistently determined by the refractory dust sublimation temperature.
can vary significantly
during the first iterations when the opacity is very high. We
adopted here a sublimation (destruction) temperature of 1000 K,
typical for silicate-like dust particles. Adopting a 500 K inner
cavity temperature threshold would not affect the SED output of the
model in the observed wavelengths, given the high visual
extinction. The constraints on the flux observed in the near
infrared for a pure 1D model is hidden by the presence of scattered
light probably originating in the outflow cavities walls. The model
radial grid contains 50 points equally spaced on a logaritmic
scale. 600 angles are considered to evaluate the angular dependent
intensity received by a grain at position r. The flux conservacy
is estimated by checking the constancy of the moment of second
order (
r2) and the iterations stopped when the value of
the standard deviation of the flux conservacy is below 5%, which
means a temperature determination to about 1%.
![]() |
Figure 2:
Absorption and scattering cross-section of the grains. The
real ( upper panel) and imaginary ( lower panel) part of the complex
refractive index of the core interstellar silicates (dot-dashed
line) and ice mantle (full line) are displayed. The silicates
constants are the Draine & Lee (1984) ones. The H2O:CO2 refractive index from 2 to 25 ![]() ![]() ![]() |
Open with DEXTER |
The grains determining the opacity are composed of pure silicates
in the region where the evaluated temperature is above ice
sublimation (100 K). The location of the ice sublimation front is
actualised at each iteration. Where the visual extinction is below
a given threshold, we consider bare silicates grains. Such an ice
mantle appearance threshold is known for ices since a long time
(e.g. Whittet & Duley 1991). We adopt an
threshold of 3,
the minimum expected threshold for the less volatile ices such as
H2O. In the shell delimited by the two above mentioned regions,
the grains are coated with an ice mantle whose adopted composition
is H2O:CO2 (100:15), in agreement with the correlations in
many lines of sight (e.g. de Graauw et al. 1996; Gerakines et al. 1999;
Watson et al. 2004), with an optimal ice mantle to refractory
silicates core volume ratio V=1.2, close to the one expected from
the correlation between visual extinction and water ice mantle
3
m OH stretching mode observations (Dartois 2005, and
reference therein). The aim of the study is not to explore the
possible detailed ice mantle CO2 local profile by extracting a
local continuum and fitting with pure absorption spectra of
laboratory ice films, like presented in Bergin et al. (2005), but
rather to make an overall SED ajustment to the observations in
order to gain insight into the quantities involved at cloud
scale. H2O:CO2 (2:1) was also tried but revealed to absorb
too much in the CO2 bending mode for such V. The refractive
index of the silicates (Draine & Lee 1984) and ice mantle (determined
from IAS laboratory transmittance experiments at 10 K, completed by
Warren (1984) and Trotta (1996) data, see
Fig. 2) were used to calculate the absorption
and scattering absorption coefficient, using the Dipole
Approximation DDSCAT program. We consider a distribution of
randomly oriented ellipsoids with a quadratic shape distribution
weighting (like in Fabian et al. 2001) with a mean size of 0.1
m
(i.e. in the Rayleigh limit up to the beginning of the
visible). The extinction from the UV to mm calculated for the bare
and coated grains are given in Fig. 3.
Starting from the best fits obtained in Shirley et al. (2002) and Jorgensen et al. (2002), we explored the p=1-2 parameter space. As described in Shirley et al. (2002) the flux coming from the ISRF in the outer part of the cloud plays an important role for the far infrared flux, as the temperature of the outer part of the envelopes raises or stay constant when approaching the external limit. To take this into account, the temperature in the cloud modelling cannot drop below 15 K, affecting mostly the cloud outer part and can be translated into the ISRF contribution discussed above.
![]() |
Figure 3: Mass absorption cross-sections adopted in the radiative transfer model. The dashed line is the absorption, thin one the scattering, and thick line the extinction cross-section. Upper panel: silicates ellipsoids with a MRN grain size distribution. Lower panel: silicates core/ ice mantle grains with a MRN grain size distribution. |
Open with DEXTER |
![]() |
Figure 4:
Model of the spectral energy distribution (SED) of L723. In addition to the Spitzer Long-Low spectrum, fluxes
originating from different observatories are plotted
(Table 1). The emergent SED calculated with the 1D transfer model using the parameters given in Table 2 is overplotted. The right panels are the ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The best fit model is shown in Fig. 4 (corresponding parameters in Table 2),
together with the minimization surfaces for p indexes of 1.,
1.5 and 1.9. These minimizations shown that the best fit is obtained
for p=1.5. The best constraint, obtained for the H2 density at
100 AU, is driven by the Spitzer mid infrared spectrum. The cloud
outer radius is much less defined if no additional constraints are put
on the model. The physical extension of the cloud observed by
Shirley et al. (2002) favors also the
1.5 case.
Table 2: L723 model's best fit parameters*.
The near infrared (NIR) excess observed with 2MASS, and not
reproduced in the framework of the 1D model, is due to
escaping photons scattered by the walls of the cavity digged by the
powerful quadrupolar outflow. Deep images in the K' band confirm that
the NIR source is extended (Lee et al. 2002). The NIR observations
can be reasonably reproduced with a stellar spectrum reprocessed by
dust, as shown in Fig. 4. The scattering efficiency drop
radically above 2
m, and the remaining flux can be well
fitted with a spherical envelope in which is embedded the young source.
The CO2/H2O ratio adopted to reproduce the observations is in
perfect agreement with the adopted ice to core volume ratio of one,
coupled to a CO2/H2O ratio of about 15%, numbers compatible
with the numerous observations performed with ISO
(e.g. Gerakines et al. 1999) and Spitzer (e.g. Watson et al. 2004). However,
the water ice features model should be constrained in additional
spectral regions. The ice region pencil beam column density, for the
best model, represents from 1/3.3th (500 K sublimation
temp.) to
1/7th (1000 K sublimation temp.) of the inner
"ice ree'' column density. In volume, the ice region largely
dominates over the inner "ice free'' one. The total pencil beam column
density is of the order of n(H2,
= 500 K)
4.9
1023 cm-2 or n(H2,
= 1000 K)
1.0
1024 cm-2). Taking into account only the ice region, locally,
10-5. If effectively CO2/H2O
0.15 in the ice region, at least 30% of the
cosmic available oxygen is locked into the ice mantles.
IRAS
19156+1906 in the L723 cloud display a huge quantity of ices
along the line of sight, as traced by the Spitzer observations of the
CO2 ice bending mode around 15.2 m. The estimated column
density of N(CO2) = 6.3
1
1018 cm-2 is one of the
highest ever detected. The modeling of the complete spectral
distribution reveals however that its relative abundance is not much higher
than what is expected from the correlations found in other lines of
sight. Nontheless it implies that about 30% of the cosmic oxygen is
trapped into ice mantles at this stage of the L723 cloud's evolution.
Acknowledgements
W. F. Thi acknowledges an ESA Research Fellowship at ESTEC. We thank A. Abergel for retrieving and reducing the ISOCAM data.