A&A 435, 161-166 (2005)
DOI: 10.1051/0004-6361:20041989
S. De Ruyter1 - H. Van Winckel2 - C. Dominik3 - L. B. F. M Waters2,3 - H. Dejonghe1
1 - Sterrenkundig Observatorium, Universiteit Gent, Krijgslaan
281 S9, 9000 Gent, Belgium
2 - Instituut voor Sterrenkunde,
KU Leuven, Celestijnenlaan 200B, 3001 Leuven, Belgium
3 -
Sterrenkundig Instituut "Anton Pannekoek'', Kruislaan 403, 1098
Amsterdam, The Netherlands
Received 10 September 2004 / Accepted 20 January 2005
Abstract
We present extended Spectral Energy Distributions (SEDs) of
seven classical RV Tauri stars, using newly obtained submillimetre
continuum measurements and Geneva optical photometry supplemented with
literature data. The broad-band SEDs show a large IR excess with a
black-body slope at long wavelengths in six of the seven stars, R Sct
being the noticeable exception. This long wavelength slope is best
explained assuming the presence of a dust component of large grains in
the circumstellar material. We show that the most likely distribution
of the circumstellar dust around the six systems is that the dust
resides in a disc. Moreover, very small outflow velocities are
needed to explain the presence of dust near the sublimation
temperature and we speculate that the discs are Keplerian. The
structure and evolution of these compact discs are as yet not
understood but a likely prerequisite for their formation is that the
dusty RV Tauri stars are binaries.
Key words: stars: AGB and post-AGB - stars: binaries: general - stars: circumstellar matter
RV Tauri stars form a distinct but loosely defined class of classical pulsating objects. They are located in the brightest part of the Population II instability strip and are defined as luminous (I-II) mid-F to K supergiants that show alternating deep and shallow minima in their light curves. They have formal periods (defined as the time between two successive deep minima) between 30 and 150 days, but often show strong cycle-to-cycle variability. RV Tauri stars are subdivided into two photometric classes on the basis of their light curves. The RVa class contains the stars with constant mean magnitude, while the RVb stars show variation of the mean magnitude on a timescale in the range of 600 to 1500 days.
Table 1:
The name, the HD number, the equatorial coordinates and
(J2000), the Galactic coordinates l and b, the
formal period P, the amplitude in the V-magnitude
,
the
photometric classification, the effective temperature
,
the surface gravity
and the
metallicity [ Fe/ H] of the seven classical
RV Tauri stars from our sample.
After Gehrz (1972) and Lloyd Evans (1985) pointed out that many RV Tauri stars present near-infrared excesses, Jura (1986) identified the low-mass RV Tauri variables as young post-AGB objects. IRAS detected in many RV Tauri stars considerable amounts of cool, circumstellar dust, which is interpreted as being a relic of the strong dusty mass-loss on the AGB.
Chemically, however, RV Tauri stars do not show the expected post-AGB abundances: no high C-abundances or s-process overabundances are observed. Instead, RV Tauri photospheres often show chemical anomalies similar to the depletion patterns seen in the gas phase of the ISM (e.g. Giridhar et al. 2000). These patterns are a result of a poorly understood chemical process in which the separation of the circumstellar dust from the circumstellar gas is followed by a selective reaccretion of only the gas, which is then rich in non-refractory elements. Waters et al. (1992) proposed that the most likely circumstance for the process to occur is when the dust is trapped in a circumstellar disc. Note, however, that there is to date no detailed quantitative model of this reaccretion scenario.
Originally, depletion patterns were found in
Boo stars and
only in five peculiar post-AGB stars, which were indeed proven
to be binaries likely surrounded by a stable circumbinary disc
(Van Winckel et al. 1995). The depletion patterns observed in many
RV Tauri stars might call for the presence of such a disc as well,
but since there is no relation found between the observed IR excess
and the chemical composition of the stellar photosphere, there is a
priori no direct observational evidence for this claim
(e.g. Giridhar et al. 2000, and references therein).
However, for individual objects like AC Her, the presence of a stable circumstellar structure has been proposed to explain the ISO SWS spectrum (Van Winckel et al. 1998) and the CO data (Jura & Kahane 1999). The disc was even claimed to be resolved in the N and Q-band (Jura et al. 2000) but this was not confirmed by higher resolution imaging (Close et al. 2003). Indirect observational evidence for the existence of a disc in more objects includes the presence of a near-IR excess in the SED (Lloyd Evans 1999). Also CO (J=1-0) rotational line emission is remarkably weak and narrow for the few objects detected (Bujarrabal et al. 1988), indicating either very slow winds or Keplerian kinematics.
The dusty particles in a disc may grow by coagulation to sizes much
larger than found in AGB outflows, and one way to identify
those large particles is to observe at long wavelengths. We have
therefore obtained submillimetre continuum observations at m
(Sect. 2) with the Submillimetre Common-User
Bolometer Array (SCUBA) at the James Clerk Maxwell Telescope (JCMT) of
7 classical RV Tauri stars (Table 1).
In Sect. 3 the SEDs of the 7 classical RV Tauri stars are presented from the UV up to submillimetre wavelengths. Although the objects are well-known RV Tauri stars, we realized that the total line-of-sight extinction of the sample was very poorly constrained in the literature. We therefore quantify in a systematic way the total extinction of all objects based on broad-band photometry. As a first-order approximation, an optically thin dust model is used to fit the observed IR excess in Sect. 4. The results are evaluated with respect to the dust geometry in Sect. 5 and in Sect. 6 conclusions are formulated.
Table 2:
Observations with SCUBA: m fluxes.
![]() |
Figure 1:
The SEDs of 6 RV Tauri stars: the dereddened fluxes
are given together with the scaled photospheric Kurucz model
representing the unattenuated stellar photosphere (solid line). An
optically thin dust fit was used to model the IR excess
(dotted line). Data found in the literature together with our 7 band
Geneva photometry (only the maxima) are plotted as triangles. The
minimal data points (squares) were not used for the determination of
E(B-V). They are only plotted to give an idea of the amplitude of
the pulsations. Crosses represent our ![]() ![]() ![]() |
Open with DEXTER |
The main difficulty in constructing the SEDs of pulsating stars with large amplitudes and in many cases severe cycle-to-cycle variability is the acquisition of equally phased data over a wide spectral domain. Since these data are not available, we limited our study of the broad-band energetics to the phase of minimal and maximal brightness.
We acquired Geneva optical photometry (Table A.1) at
random epochs for RV Tau (53 measurements), SU Gem (29),
U Mon (119), AC Her (62) and R Sct (36) with the Flemish
Mercator Telescope at La Palma, using the refurbished Geneva
photometer P7 (Raskin et al. 2004). We added additional Geneva optical
photometry from the Geneva database
(http://obswww.unige.ch/gcpd/gcpd.html). Johnson and Cousins
broad-band photometry was found in the literature
(Table A.2). For TW Cam and UY CMa no Geneva
photometry is available, so we used only literature data. Near-IR data
were taken from the 2 MASS and DENIS projects complemented with data
from the literature (Table A.3). Far-IR data come
mainly from IRAS (Table A.4) while for SU Gem, U Mon and
R Sct, MSX data are also available (Table A.5). At the
short wavelength side we used the IUE (m-0.320
m)
data from the newly extracted spectral data release (INES). The
resulting SEDs are given in Figs. 1 and 2.
![]() |
Figure 2: The SED of R Sct is clearly different from these of the other programme stars. Symbols are the same as in Fig. 1. |
Open with DEXTER |
Table 3:
The five parameters of the dust model: the
normalization temperature T0, the inner radius of the dust shell
,
the outer radius of the dust shell
,
the spectral index p and the density parameter m. The last column lists the energy ratio
.
Because of the circumstellar material, the total extinction must contain both a circumstellar and an interstellar component. Whether the circumstellar extinction affects the line-of-sight depends on the geometry. The estimate for E(B-V) is found by minimizing the difference between the dereddened observed fluxes in the UV-optical and the appropriate Kurucz model (Kurucz 1979). We scale our data to the J filter, the reddest filter for which no dust excess is expected. For the dereddening of the observed maximal fluxes we used the average interstellar extinction law given by Savage & Mathis (1979). This implies that we also adopt the interstellar extinction law for the circumstellar component of the total extinction.
The total colour excess for all objects lies between 0.0 and 0.8(Table 4). The intrinsic error on this E(B-V)-determination is typically 0.1. Together with the error of 0.2 induced by the uncertainty of the temperature of the underlying photosphere, we have an uncertainty of 0.3 on the total extinction during maximal light.
Despite the significant IR excess, the total line-of-sight extinction is very small for UY CMa, U Mon and R Sct, even within the error. While for R Sct the dust excess is marginal and the large amplitude pulsations make the determination of the total reddening uncertain, the very large IR excess of UY CMa and U Mon is clearly in contradiction with the lack of line-of-sight reddening. This implies either a significant grey extinction or a non-spherically symmetric circumstellar dust distribution.
Table 4:
The estimated total colour excess E(B-V) and the optical
depth estimates in the V-band.
is determined from
the SED minimalization and gives the total line-of-sight
extinction. The contribution from the ISM is given in Col. 3
assuming the distance range given in
Table 5.
was determined
by integrating the spherically symmetric dust model and in the last
column the mean optical depth is given, determined by the energy ratio
from the dust excess over the optical fluxes.
Table 5:
The bolometric corrections BCV (from Bessell et al. 1998),
the luminosity L (in
), the distance D (in kpc)
and the dust mass (
)
are given, together with the
estimated errors.
An important conclusion of our dust fitting is also that the
inner radii
are very small: in all cases less than
14 AU from the central object. The dust is located very close to the
star and dust excesses start near the sublimation temperature. The
distance of the inner boundary to the stellar surface would be covered
by a normal AGB wind velocity of 10 km s-1 in only about 7 years. If indeed the dust we observe is a relic of the AGB mass-loss,
the expansion velocity of the circumstellar material must be extremely
small! Note that there is no evidence of current day dusty
mass-loss. Moreover, RV Tauri stars with very similar pulsational
amplitudes and spectral types exist without detected infrared
excesses.
Although the outer radius estimates are not well constrained and are
very dependent on the exact choice of the spectral index, it is clear
that the outer radii are not very large either. In most cases, the
dust can be fitted well, even with a single temperature
black-body (e.g. Dominik et al. 2003). This is also clear from visual
inspection of the SED plots, which show a lack of cool
(
K) dust. If the circumstellar
geometry is optically thick, the extent of the dust distribution would
be even more compact than in the optically thin case. Given the
distances to the objects, interferometric observations will be needed
to resolve the circumstellar dust geometry.
Nevertheless, the different optical depths in Table 4 yield significant information on individual objects. At first sight, only the total reddening of AC Her is the sum of the ISM and the circumstellar component. UY CMa, U Mon and RV Tau are examples for which the large IR excess is inconsistent with the small total line-of-sight reddening, pointing to a strongly non-spherical geometry of the circumstellar material. For SU Gem the total reddening is larger than the sum of the likely ISM contribution and the circumstellar extinction expected from a spherical dust shell. For TW Cam, the large expected ISM contribution indicates that our distance estimates are too large. The empirical extinction model from Neckel et al. (1980) shows a steep rise up to about 1 kpc in the direction of TW Cam. Again also in this star the circumstellar line-of-sight extinction must be very small.
R Sct is a case for which the period-luminosity relation does not
hold. The Galactic coordinates imply a significant mean
interstellar component of the extinction using the distance obtained
with the P-L relation (Table 4). Even with the upper limit
of the Hipparcos parallax, one expects the ISM extinction to be
significant (
using the mean exinction
estimates from Hakkila et al. (1997) with a distance of 320 pc).
R Sct is known to be a very irregular pulsator with a large
amplitude, but strict simultaneous photometric studies over a wide
wavelength domain show a small total extinction (e.g.
E(B-V)=0.3-0.5, Shenton et al. 1994) for which the interstellar
component is low (Cardelli 1985,
). We
conclude that R Sct is much less luminous than given in
Table 5. This RV Tauri object is, in line with other
studies (e.g. Matsuura et al. 2002), considered as an exception.
Our energy balance calculations indicate a non-spherical distribution of the circumstellar dust in 5 sources, with the possible exception of SU Gem, where we are likely looking into the disc.
The amount of energy reprocessed by the dust grains
(
)
is larger than 34% for all objects,
except for U Mon, AC Her and R Sct (Table 3). The
mean optical depth of the circumstellar dust (Table 4) is
defined as
.
Assuming
that the IR emission is produced by an infinitely optically
thick disc, the opening angle of the disc as seen from the star must
be
for objects with 40% flux
conversion. The significant scaleheight is likely to be
sustained by gas pressure in the gas-rich discs
(e.g. Dominik et al. 2003).
Despite the large scaleheights, we still prefer to use the word "disc'' instead of torus. Dust tori are often resolved around Proto Planetary Nebulae, but these have very different SED characteristics since they are much colder and are likely expanding (e.g. Sánchez Contreras & Sahai 2004). Moreover, the physical sizes of the resolved tori are much larger than we expect the circumstellar material to be in the RV Tauri stars.
In each case, the dust mass lies within the range of
and
,
with the noticeable exception
of R Sct for which
is found
(Table 5).
The most obvious way to explain the formation of a disc around evolved RV Tauri stars is to assume that the disc was created during binary interaction, when the primary star was at giant dimensions. Binarity in RV Tauri stars, however, is hard to prove with radial velocity monitoring since the photospheric pulsations have large amplitudes in radial velocity. Only in a few examples was orbital motion indeed found: AC Her (Van Winckel et al. 1998), U Mon (Pollard & Cottrell 1995) and RU Cen and SX Cen (Maas et al. 2002). In this analysis we confirm the proposition made earlier (Van Winckel et al. 1999) that the six well-studied dusty RV Tauri stars are likely all binaries in which the binary interaction has played a fundamental role in creating a stable dusty disc.
The actual structure of the disc, as well as its formation, stability and evolution are not well understood at this stage. We argue that realistic models of circumstellar material of RV Tauri stars should be constructed in the framework of stable discs, similar to the post-AGB star HR 4049 (Dominik et al. 2003).
Acknowledgements
The staff and service observers of the Mercator Observatory at La Palma are acknowledged for the photometric Geneva data. This publication makes also use of data from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. We also used data from the DENIS project, which is partly funded by the European Commission through SCIENCE and Human Capital and Mobility grants. It is also supported in France byINSU, the Education Ministry and CNRS, in Germany by the Land of Baden-Würtenberg, in Spain by DGICYT, in Italy by CNR, in Austria by the Fonds zur Förderung der Wissenschaftlichen Forschung and Bundesministerium für Wissenschaft und Forschung. We also like to thank the anonymous referee for his contribution to this paper.
Table A.1: Geneva data were acquired with the Flemish Mercator Telescope at La Palma, Spain. Our total dataset was scanned for the maximum and minimum magnitudes. Observation dates, number of measurements and total timebases of these maxima and minima are given as well. For AC Her, our observational monitoring data are divided in two major blocks of 8 and 446 days with 4612 days inbetween. Additional data were found in the Geneva database (http://obswww.unige.ch/gcpd/gcpd.html).
Table A.2: Ground-based optical data, acquired over a long period, found in the literature.
Table A.3: Ground-based near-IR photometry found in the literature. If there was more than one measurement we used both the maximum and minimum data points. For the data from 2 MASS and DENIS we made use of the catalogues found in VIZIER (http://vizier.u-strasbg.fr/viz-bin/VizieR).
Table A.4:
IRAS photometry points at 12, 25, 60 and m. Note however that in some cases the
m is an upper
limit (L); these observations are probably contaminated by
interstellar cirrus clouds. For RV Tau the
m data point is
a lower limit (:).
Table A.5:
Data of the Midcourse Space eXperiment (MSX). The instrument
on board MSX is the SPIRIT III (Spatial Infrared Imaging Telescope
III). The approximate effective wavelengths of the 6 MSX filters are in m.