A&A 402, 183-187 (2003)
DOI: 10.1051/0004-6361:20030194
R. Liseau1 - A. Brandeker1 - M. Fridlund2 - G. Olofsson1 - T. Takeuchi3 - P. Artymowicz1
1 - Stockholm Observatory, SCFAB, Roslagstullsbacken 21, 106 91 Stockholm, Sweden
2 -
ESTEC/ESA, PO Box 299, 2200AG Noordwijk, The Netherlands
3 -
Lick Observatory, University of California Santa Cruz, CA 95064, USA
Received 26 July 2002 / Accepted 5 February 2003
Abstract
We present millimeter imaging observations in the 1200 m continuum
of the disk around
.
With the 25
beam, the
disk is
unresolved perpendicularly to the disk plane (
10
), but slightly resolved
in the northeast-southwest direction (26
). Peak emission is observed at the stellar
position. A secondary maximum is found 1000 AU along the disk plane in the southwest,
which does not positionally coincide with a similar feature reported earlier at 850
m.
Arguments are presented which could be seen in support of the reality of these features.
The observed submm/mm emission is consistent with thermal emission from dust grains, which are
significantly larger than those generally found in the interstellar medium, including
mm-size particles, and thus more reminiscent of the dust observed in protostellar disks.
Modelling the observed scattered light in the visible and the emission in the submm/mm
provides evidence for the particles dominating the scattering in the visible/NIR and
those primarily responsible for the thermal emission at longer wavelengths belonging
to different populations.
Key words: stars: individual: Pictoris - circumstellar matter -
planetary systems: formation - protoplanetary disks - ISM: dust, extinction
Since its discovery by IRAS (e.g., Aumann 1984), the
toris system has presented
the prime example of a dusty disk around a main sequence star, partly
because of its high degree of "dustiness'' (
,
e.g. Lagrange et al. 2000) and partly because of its relatively close distance to the Earth
(19.3 pc, Crifo et al. 1997), which makes it possible to obtain high quality data over
the entire spectral range. Recent papers reviewing the physics of
the disk around
include those of Artymowicz (2000), Lagrange et al. (2000) and Zuckerman (2001).
Considerable uncertainty existed regarding the age of the system, but most recent estimates place
the stellar age close to only ten million years (
12+8-4 Myr, Zuckerman et al. 2001). This
could open up the possibility that planet formation (nearing its final phases?) might actually become observable.
Since these reviews were written, new relevant information has been added to our knowledge of
the
system: Olofsson et al. (2001) reported the discovery of widespread atomic gas in the
disk, a result recently confirmed and extended by Brandeker et al. (2002). These observations revealed
the sense of disk rotation and that the northeast (NE) part
of the gaseous disk is extending to the limit of the observations by Brandeker et al., viz. to at
least 17
(330 AU) from the star. Several difficulties were encountered with these discoveries,
such as the observed fact that the gas stays on (quasi-)Keplerian orbits, although
radiation pressure forces in the resonance lines should accelerate the gas to high velocities
and remove it on time scales comparable to the orbital period.
This needs to be addressed in the context of
the origin and evolution of the gas and dust, i.e. whether one or both components
are presently produced in situ in the disk or whether they (at least to some degree)
constitute "left-overs'' from the star formation process, being of primordial origin.
Takeuchi & Artymowicz (2001) have considered the interaction of gas and dust in a
circumstellar disk, the dynamical evolution of which is critically dependent on the relative
abundance of these species (see also Lecavelier des Etangs et al. 1998). Possible observational consequences,
even relatively far from the central star, may become assessible with modern mm/submm cameras.
With the aim to compare and to extend the results obtained at 850 m with SCUBA by
Holland et al. (1998), we performed imaging observations at longer wavelengths and, in this paper,
we present the image of the
disk at 1200
m (1.2 mm).
The SIMBA observations and the reduction of the data are presented in Sect. 2. Our basic result,
i.e. the 1.2 mm image of
and its circumstellar disk, is found in Sect. 3, and in Sect. 4
the possible implications of these observations are discussed, where also other data are consulted.
Finally, in Sect. 5, we briefly summarise our main conclusions from this work.
![]() |
Figure 1:
The normalised SIMBA 1.2 mm image with 8
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 2:
The
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The observations were performed with the SEST Imaging Bolometer Array (SIMBA) at the 15 m
Swedish-ESO Submillimeter Telescope (SEST), La Silla, Chile, during the periods November 17-23
and 27-29, 2001. SIMBA has been developed by the Max-Planck-Institut für Radioastronomie, Bonn,
in collaboration with the Astronomisches Institut der Ruhr-Universität Bochum.
The 37 liquid helium cooled semiconductor elements are n-doped silicon chips, mounted on a
saphire substrate. Projected on the sky, the 37 horn antennae have each an HPBW of 23
and
are spaced by 44
in a hexagonal arrangement, covering about 4
.
The spectral band
pass is centered at 250 GHz (1200
m) and has a full width at half maximum (FWHM) of 90 GHz.
As map centre, the equatorial coordinates of
were used, viz.
05
47
17
1,
-51
03
59
(J 2000). Generally,
the mapping was done in fast mode by scanning 600
in azimuth and 400
in elevation
at the rate of 80
per second and with a step size in elevation of 8
,
oversampling
the beam by about a factor of three. The zenith opacities were obtained by means of frequent
sky dips and were on the average 0.22 at the beginning of our observing run, but improved
to about 0.15 after a couple of days. The pointing of the telescope was regularly checked using
the extragalactic radio source
0537-441 and/or the planet Saturn and was found to be within a
third of the beamwidth.
Flux calibrations were based on observations of the planet Uranus. The brightness temperature of
94.8 K at 250 GHz implies a Uranus flux density of 38 Jy/beam, with an
uncertainty of 5% (Griffin & Orton 1993). As is evident from Fig. 1,
the spatial flux distribution of this point source
(diameter
)
is consistent with a circular Gaussian telescope beam, with the measured
FWHM(
(24
9
0
3,
0
).
The deviation from a circularly symmetric Gaussian beam pattern occurs at the 1% level of the peak value
(-20 db).
In total, we recorded 171 maps with 21 hours of integration.
These data were reduced by making use of the MOPSI reduction software
package
developed by Robert Zylka. This involved correcting
for atmospheric extinction, cosmic rays and the variable sky background,
as well as producing maps from the fast scanning mode. The sky noise
was greatly reduced by correlating and removing the simultaneous flux
level variations of the different bolometer channels.
The 1200 m image of the
region is displayed in Fig. 2, where the
displayed contours are chosen in compliance with Holland et al. (1998). The flux maximum,
designated as feature A in the figure, is centered on the position of the star
(0
, 0
).
At position angle 237
relative to A, an elongation of the emission in the NE-SW direction is
discernable. For an assumed Gaussian source flux distribution, the deconvolved minor
and major axes are
10
and 26
,
respectively, the latter corresponding to 500 AU.
The disk emission extends to at least 55
(1050 AU) in the SW direction. In addition,
faint emission protrudes south toward a "blob'' at (-27
,
-44
). This feature C,
at position angle 211
5, is 52
(1000 AU) distant from the star and is thus not
(whether real or not) positionally coincident with the 850
m-SW blob of Holland et al. 1998 at
(-21
, -26
)
and identified as B in Fig. 2, but we note the following important facts:
(1) this faint emission, extending straight south, is also apparent in the 850
m image,
(2) as is evident in Fig. 2, the position angle to blob C coincides with that of
the midplane of the
disk inferred from optical data (= 31
5 + 180
;
Kalas & Jewitt 1995, Heap et al. 2000), and (3) faint dust scattered light in this
direction has been observed far from the star (1450 AU, Larwood & Kalas 2001). Toward the NE,
the scattering disk has been claimed to extend even further from the star, to 1835 AU.
An about
feature at (+48
,
+32
)
and PA = 57
(1100 AU)
is discernable in our image, but which would again be significantly offset from a similar blob
at (+28
,
+25
)
in the 850
m data. In order to assess the reality of these faint features
we have divided the raw data into different portions and then applied the same reduction procedures
to these subsamples. The result of this exercise indicates that the (+48
,
+32
)
feature
is an artefact introduced by the noise, as it is not seen in all frames, whereas blob C is persistently
present in the final sub-maps.
The flux densities of both the
disk (A) and the SW features (B and C)
are presented in Table 1. The result for the peak flux obtained by Chini et al. (1991)
with a single element bolometer at the SEST at 1300
m, viz.
mJy/beam, is
in agreement with our array value of
mJy/beam at 1200
m (cf. Fig. 3).
The ratio of the peak flux A to that of blob C is
.
For widely adopted parameters, the stellar disk subtends an angle in the sky of less than 1 mas and the
photosphere generates a flux density at the Earth of less than 1 mJy at 1200 m. The stellar contribution
to our SEST measurements can therefore be safely ignored. Also, any line emission in this band pass is
likely to be totally negligible (Liseau & Artymowicz 1998; Liseau 1999 and references therein).
Assuming an opacity law of the form
,
the flux density
at long wavelenghts from an optically thin source of a certain dust population can be expressed as
![]() |
(1) |
![]() |
(2) |
Feature | Relative |
![]() ![]() |
Remarks |
Offset (
![]() |
(mJy/beam) | ||
A | (0, 0) |
![]() |
![]() |
![]() |
integrated over a radius of 40
![]() |
||
B | (-21, -26) | ![]() |
SW blob (Holland et al. 1998): |
contaminated by
![]() |
|||
C | (-27, -44) |
![]() |
SW blob (this paper) |
The asymmetric flux distribution displayed in Fig. 2 may be surprising. In agreement with
our observations, Chini et al. (1991) and Dent et al. (2000)
too were unable to detect any emission in the northeast part of the disk,
where we place a -upper limit on the mass of 0.2
(see below; the dust temperature at 500 AU
K, when
K, see Liseau & Artymowicz 1998).
Given the low signal-to-noise ratio (S/N), the reality of this lopsidedness is difficult to assess, but
asymmetries in the
disk have been noticed also at other wavelengths. For instance, in scattered
light, the receding NE side of the disk extends much further and is much brighter than the SW disk.
In contrast, the shorter, approaching SW disk seems much thicker (Kalas & Jewitt 1995; Larwood & Kalas 2001). The
situation is reversed in the thermal infrared (albeit on smaller spatial scales), where the SW disk
appears significantly brighter and more extended than the NE side (Lagage & Pantin 1994; Wahhaj et al. 2002;
Weinberger et al. 2002). This could be due to a "Janus-effect'', i.e. the NE being dominated
by "bright'' dust particles (high albedo, silicates), whereas in the SW, the majority of dust grains is "dark''
(high absorptivity, carbonaceous?). What would accomplish such uneven distribution in the disk
is not clear, but large differences in albedo, by more than one order of magnitude, are not uncommon,
for instance, in solar system material. Also, in order to understand the nature of feature C
(and B) velocity information would be valuable.
Blob C would be situated in the disk midplane and on the second contour in the scattered light image of Larwood & Kalas (2001) ( 22 < R < 25 mag/arcsec2, see also Kalas & Jewitt 1995). No obvious distinct feature is seen at its position. However, to be detectable with SIMBA at the SEST, any point source at mm-wavelengths would not be point-like at visual wavelengths, and its optical surface brightness could be very low.
To gain some
quantitative insight we ran numerical models, exploiting Mie-theory, for a variety of plausible dust
mixtures regarding the chemical composition and grain size distributions (for details, see Pantin et al. 1997).
The equilibrium temperatures were found to be in the interval 16 K to 58 K
and depending on the dust albedo and scattering phase function,
the predicted integrated scattered light, to be consistent with our SEST observations, spans 8 magnitudes.
The two most extreme cases considered were (1) bright cold dust (albedo = 1 at visual wavelengths,
K) which is
scattering isotropically and has albedo = 0.2 at thermal wavelengths, which results in a spatially integrated
R-magnitude of 17.3 and (2) dark warm dust (albedo = 0.02 in the visual,
K) which
gives rise to "comet scattering'', i.e. 14% of isotropic at 90
,
and has zero albedo at thermal wavelengths,
resulting in an integrated R-magnitude of 25.6. These extreme cases are felt to be either too optimistic
or too conservative and an intermediate case might be more appropriate. Our adopted model includes
isotropically scattering dust at
K with albedo = 0.2 at visual and zero
albedo at thermal wavelengths, yielding an integrated R-magnitude of 20. An about 10
source
(R=25 mag/arcsec2) would thus be consistent with both the optical data and our SIMBA measurement
and a deep R-band search might become successful. Similarly, the integrated 850
m flux density is
predicted to be slightly less than 18 mJy and should become readily detectable. However, blob C is
situated outside the figure of Holland et al. 1998.
Feature C is perfectly aligned with the optical disk plane, and its (hypothetical) mass can
be estimated from
![]() |
(3) |
![]() |
Figure 3:
The SED of the
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
For radiation mechanisms generating power law spectra and/or thermal dust emission from blob B,
the spectral slope is given by
,
where
,
as before, refers to A = (0
,
0
), i.e. the
disk, and where
and
m or 1300
m. Because of the relatively low S/N of the SCUBA and SIMBA data, the actual flux ratios are
highly uncertain and, furthermore, calibration uncertainties and telescope beam effects could
potentially introduce large errors. The combined observations of blob B
suggest
,
yielding
,
i.e. consistent
with the spectrum of A.
Based on 1200 m imaging observations of the circumstellar dust of
we conclude the following:
Acknowledgements
We are grateful to Dr. José Afonso, who made available to us the Uranus calibration data, and to the SEST staff for providing additional observations. The critical comments by the anonymous referee are highly appreciated.