Issue |
A&A
Volume 507, Number 2, November IV 2009
|
|
---|---|---|
Page(s) | 929 - 937 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/200912885 | |
Published online | 15 September 2009 |
A&A 507, 929-937 (2009)
Orbital periods of cataclysmic variables identified by the SDSS
VI. The 4.5-h period
eclipsing system SDSS J100658.40+233724.4![[*]](/icons/foot_motif.png)
J. Southworth1 - R. D. G. Hickman1 - T. R. Marsh1 - A. Rebassa-Mansergas1,2 - B. T. Gänsicke1 - C. M. Copperwheat1 - P. Rodríguez-Gil3,4
1 - Department of Physics, University of Warwick, Coventry, CV4 7AL, UK
2 - Departamento de Física y Astronomía, Universidad de Valparaíso, Avenida Gran Bretana 1111,
Valparaíso, Chile
3 - Isaac Newton Group of Telescopes, Apdo. de Correos 321, 38700 Santa Cruz de La Palma, Spain
4 - Instituto de Astrofísica de Canarias, vía Líctea, s/n, La Laguna, 38205 Tenerife, Spain
Received 14 July 2009 / Accepted 13 September 2009
Abstract
We present time-resolved spectroscopy and photometry of
SDSS J100658.40+233724.4, which we have discovered to be an
eclipsing cataclysmic variable with an orbital period of
0.18591324 days (267.71507 min). The observed velocity
amplitude of the secondary star is km s-1, which an irradiation correction reduces to
km s-1.
Doppler tomography of emission lines from the infrared calcium triplet
supports this measurement. We have modelled the light curve using the LCURVE code and Markov Chain Monte Carlo simulations, finding a mass ratio of
.
From the velocity amplitude and the light curve analysis we find the mass of the white dwarf to be
and the masses and radii of the secondary star to be
and
,
respectively. The secondary component is less dense than a normal main
sequence star but its properties are in good agreement with the
expected values for a CV of this orbital period. By modelling the
spectral energy distribution of the system we find a distance of
pc and estimate a white dwarf effective temperature of
K.
Key words: stars: novae, cataclysmic variables - stars: dwarf novae - stars: binaries: eclipsing - stars: binaries: spectroscopic - stars: white dwarfs - stars: individual: SDSS J100658.40+233724
1 Introduction
Cataclysmic variables (CVs) are interacting binary systems containing a low-mass secondary star losing material to a white dwarf primary star. The Sloan Digital Sky Survey (SDSS) has spectroscopically identified 252 of these objects, 204 of which are new discoveries (see Szkody et al. 2009, and references therein). This sample of SDSS CVs is valuable because of its large size and homogeneity (Gänsicke et al. 2009) and we are undertaking a project to characterise its constituent objects (see Dillon et al. 2008; Southworth et al. 2008a,2006; Gänsicke et al. 2006; Southworth et al. 2008b, and references therein). In the course of this work we have discovered that SDSS J100658.40+233724.4 (hereafter SDSS J1006) shows deep eclipses which are identifiable both spectroscopically and photometrically. The presence of eclipses allows us to determine the basic physical properties of the system, information which is difficult or impossible to obtain for the great majority of CVs (Littlefair et al. 2006; Smith & Dhillon 1998; Knigge 2006).
SDSS J1006 was discovered to be a CV by Szkody et al. (2007) on the basis of an SDSS spectrum which showed strong and wide Balmer emission lines. The continuum is blue at bluer wavelengths but clearly shows the spectral features of the M-type secondary star at redder wavelengths. SDSS J1006 is one of a select bunch of long-period CVs in which the eclipse of the white dwarf (WD) is directly visible in the light curve. In this work we present and analyse time-resolved spectroscopy and photometry of SDSS J1006, from which we measure the masses and radii of the WD and secondary star.
2 Observations and data reduction
Table 1: Log of the observations presented in this work. The mean magnitudes are calculated excluding observations taken during eclipse.
![]() |
Figure 1: Flux-calibrated average spectrum of SDSS J1006. Data from the blue arm of ISIS are shown in the left panel, and from the red arm in the right panel. The most prominent emission and absorption features are labelled. |
Open with DEXTER |
2.1 Spectroscopy
Spectroscopic observations were obtained in 2008 January, using the ISIS double-beam spectrograph on the William Herschel Telescope (WHT) at La Palma (Table 1). For the red arm we used the R316R grating and Marconi CCD binned by factors of 2 (spectral) and 3 (spatial), giving a wavelength range of 6115-8840 Å at a reciprocal dispersion of 1.85 Å per binned pixel. For the blue arm we used the R600B grating and EEV12 CCD with the same binning factors as for the Marconi CCD, giving a wavelength coverage of 3575-5155 Å at 0.88 Å per binned pixel. From measurements of the full widths at half maximum of arc-lamp and night-sky spectral emission lines, we find that this gave resolutions of 3.5 Å (red arm) and 1.8 Å (blue arm).
Data reduction was undertaken by optimal extraction (Horne 1986) as implemented in the PAMELA code (Marsh 1989), which also uses the STARLINK
packages FIGARO and KAPPA; further details can be found in Southworth et al. (2007b,a).
Copper-neon and copper-argon arc lamp exposures were taken every hour
during our observations and the wavelength calibration for each science
exposure was linearly interpolated from the two arc observations
bracketing it. We removed the telluric lines and flux-calibrated the
target spectra using observations of Feige 110, treating each
night separately.
The averaged WHT spectra are shown in Fig. 1. Trailed greyscale plots of the phase-binned spectra are shown in Fig. 2, for the H,
He I 6678 Å and Ca II 8662 Å emission lines, and Na I 8183 and 8194 Å absorption lines, and will be discussed in Sect. 3.
![]() |
Figure 2:
Greyscale plot of the continuum-normalised
and phase-binned trailed spectra of SDSS J1006.
From left to right the panels show H |
Open with DEXTER |
2.2 Photometry
Time-series photometry of SDSS J1006 was obtained in service mode using two telescopes equipped with imaging spectrographs: the Nordic Optical Telescope (NOT) with ALFOSC, and the Calar Alto (CAHA) 2.2 m telescope with CAFOS. For the NOT observations we used the No. 92 filter, which has a wide-V passband with points of half transmission at approximately 4400 and 7000 Å. The CAHA observations were made unfiltered to maximise throughput. The CCDs were mostly binned and windowed to reduce readout time, and short exposure times were used to maximise the cadence of the observations.
The 2008 December observations obtained with the NOT were reduced with the pipeline described by Southworth et al. (2009b,a), which uses an IDL implementation of DAOPHOT to perform aperture photometry. The remaining photometric data were reduced using the pipeline described by Gänsicke et al. (2004), which performs bias and flat-field corrections within MIDAS and aperture photometry with the
SE XTRACTOR package (Bertin & Arnouts 1996). Instrumental differential magnitudes were converted into V-band apparent magnitudes, using the SDSS ugriz apparent magnitudes of several comparison stars and the transformation equations given by Lupton
.
The light curves are plotted in Fig. 3 with the measured eclipse midpoints indicated. It is apparent from this plot that the depth of the eclipse is dependent on the wavelength of observation: the Calar Alto data were unfiltered, so are more affected by the light of the secondary star and thus show shallower eclipses.
![]() |
Figure 3: Plot of the four light curves obtained covering eclipses of SDSS J1006. The eclipse midpoints have been aligned on the panels. |
Open with DEXTER |
3 Analysis
3.1 Orbital ephemeris
Our first observations of SDSS J1006 were spectroscopic. Radial velocities (RVs) measured from the emission lines (see below) clearly showed anomalies due to three eclipses, on an unambiguous period of 267.9 min. The resulting preliminary ephemeris was sufficiently accurate for us to photometrically observe eclipses, on which more precise period measurements could be based.
For each of the four eclipses, a mirror-image of the light
curve was shifted until the two representations of the central parts of
the eclipse were in the best possible agreement. The time defining the
axis of reflection was taken as the midpoint and uncertainties were
estimated based on how far this could be shifted before the agreement
was clearly poorer. We have fitted a linear ephemeris to these times of
minimum light, finding

where E is the cycle number and the parenthesised quantities indicate the uncertainty in the last digit of the preceding number. This corresponds to an orbital period of

Table 2: Times of eclipse for SDSS J1006 and the residuals with respect to the linear ephemeris given in Sect. 3.1.
3.2 Emission-line radial velocities
![]() |
Figure 4:
Diagnostic diagram showing the
variation of the best-fitting spectroscopic orbital parameters
for RVs measured with a range of separations using the double
Gaussian function.
|
Open with DEXTER |
The spectrum of SDSS J1006 (Fig. 1) shows strong emission at the wavelengths of the hydrogen Balmer lines and some helium lines. These emission lines are produced by the accretion disc which surrounds the WD, so variations in their velocity hold information on the motion of the WD itself. However, spectroscopic studies of CVs often show a phase difference between the RV variation of emission lines and the orbital phases measured using other methods (Thorstensen 2000; Unda-Sanzana et al. 2006; Thoroughgood et al. 2005; Steeghs et al. 2007). This casts doubt on whether emission lines are good indicators of the motion of WDs in CVs, and due to this we did not use emission-line RVs in calculating the physical properties of SDSS J1006.
We measured RVs from the H
emission, which is the strongest emission line, using the double-Gaussian method (Schneider & Young 1980) as implemented in MOLLY. The full width half maximum of the Gaussian functions was set to 300 km s-1,
which is a good compromise between resolving emission-line features and
minimising the random noise in the RV measurements. The separation
of the two Gaussians,
,
was varied from 800 to 3000 km s-1 in jumps of 100 km s-1. For each value of
a spectroscopic orbit was fitted to the measured RVs using the SBOP
code, which we find gives reliable error estimates for the optimised parameters (Southworth et al. 2005). The orbital period was fixed at the ephemeris value (Sect. 3.1),
a circular orbit was assumed, and the phase zeropoint was included as a
fitted parameter. RVs between phases 0.9 and 0.1 were
rejected as they are affected by the eclipse of the accretion disc by
the secondary star.
We have constructed a diagnostic diagram (Shafter 1983; Shafter et al. 1986) for SDSS J1006 (Fig. 4), which shows that the properties of the spectroscopic orbit change only slowly for
-2200 km s-1, and that the lowest scatters in the residuals (
)
occurs for
-2000 km s-1.
The offset between the orbital phase and the phase of greatest negative
change in the RVs is only about 0.04 for these separations, which
indicates that the emission-line RVs might trace the motion of the WD
with reasonable accurately. We have adopted the spectroscopic orbit for
km s-1, which gives the lowest
,
and these quantities are given in Table 3. The RVs and best fit are shown in Fig. 5. Our error estimates include the standard errors given by SBOP, plus a contribution to account for variations between the solutions for
-2200 km s-1 (where the
values are the lowest). We also calculated a diagnostic diagram for the H
emission line, which yielded similar results but a greater scatter due to the weaker emission-line flux.
Table 3: Best-fitting spectroscopic orbits found using SBOP.
![]() |
Figure 5:
The measured RVs (circles) and the spectroscopic
orbits fitted to them (solid curves). The emission-line RVs (filled circles)
were calculated using
|
Open with DEXTER |
3.3 Absorption-line radial velocities
The secondary component of SDSS J1006 is clearly visible in our red-arm WHT/ISIS spectra, but very few features can be seen by the naked eye to vary in velocity, due to the modest signal-to-noise ratio of individual spectra. This velocity variation plays a vital role in constraining the properties of the system, so we have used two methods to tease out the absorption-line velocity amplitude.
Firstly, the observed spectra were augmented with a set of
template M dwarf spectra from the SDSS, then velocity-binned and
subjected to a cross-correlation analysis. The orbital ephemeris was
fixed to the numbers in Sect. 3.1
after verifying that this does not cause a significant change in the
results. The cross-correlation functions were examined interactively
and measured for velocity if they contained a clear peak from the
secondary star, and the resulting RVs were fitted with a circular orbit
using SBOP. We did this for many different spectral
regions and template spectra, finding that the resulting velocity
amplitudes were always in the interval 270-282 km s-1. For illustration, in Fig. 5 we plot the absorption-line RVs found using an M4 spectral template and the full red-arm wavelength interval (with the H
and helium emission lines masked out).
The second method is designed to cope well with spectra of a low
signal-to-noise ratio, and to target the strongest spectral features
observed to come from the secondary star. Using the MGFIT routine in MOLLY
we fitted a double Gaussian function plus spectroscopic orbit to the
sodium doublet at 8183.3 and 8194.8 Å. All 79 spectra were
fitted simultaneously, yielding a direct measurement of the velocity
amplitude:
km s-1. Figure 2
shows the phase-binned and trailed spectra of SDSS J1006 in the
region of the Na doublet. We have been unable to completely remove the
effects of telluric absorption from our spectra, so have also performed
fits with extra Gaussians added to account for the residual absorption.
We find that our K2 measurement is not significantly affected.
Given the good agreement between the two methods, we adopt a value of
km s-1,
where the error estimate accounts for both the random errors and the
variation in results from different analysis techniques (Table 3). The two methods agree well on the value of K2
but produce slightly discrepant systemic velocity measurements. This is
likely due to difficulties in placing the continuum, due to the complex
spectrum of the secondary star. We adopt a value of
km s-1, which encompasses most of the systemic velocities found during our analysis. A better measurement of
will require further data.
3.3.1 K-correction for the absorption-line velocities
![]() |
Figure 6: Variation of the equivalent widths of the Na doublet and TiO molecular band with orbital phase. The spectra were combined into 20 phase bins prior to measurement. The black filled circles represent measurements outside primary and secondary eclipses, and the grey open circles those within eclipse. The wavelength intervals over which the equivalent widths were measured were 8180-8200 Å for Na and 7090-7345 Å for TiO. |
Open with DEXTER |
Our measured K2 cannot be assumed to represent the motion of the centre of mass of the secondary star,
,
due to irradiation of the inner hemisphere by the WD and accretion
disc. The irradiated surface has a lower vertical temperature gradient
and thus weaker absorption lines. RV measurements from these lines are
therefore skewed towards the motion of the outward-facing part of the
star, causing K2 to overestimate
(Wade & Horne 1988; Billington et al. 1996; Hessman et al. 1984).
To estimate the correction
we have measured the equivalent widths of absorption features arising
from the secondary star, as a function of orbital phase. The wavelength
scales of the spectra were moved to shift out the motion of the star,
and the spectra were then rectified to a continuum level of 1 and
binned into twenty phase intervals. The resulting plots (Fig. 6)
show that the equivalent widths vary by approximately 30% outside
eclipse. Extrapolating to the phase of secondary mid-eclipse and
considering the errors on this approach and on the equivalent width
measurements, we find a total variation in equivalent width (and thus
in the light from the secondary star) of
%.
Wade & Horne (1988) obtained the expression

where f is the size of the displacement as a fraction of








3.4 Doppler tomography
![]() |
Figure 7:
Doppler maps of H |
Open with DEXTER |
To investigate the properties of the accretion disc of SDSS J1006
we have constructed Doppler maps of several of the emission lines using
the maximum entropy method (Marsh & Horne 1988). The maps are shown in Fig. 7 and phase-binned and trailed plots of the emission lines are shown in Fig. 2. The value
for the Doppler maps were chosen to be marginally larger than the
values for which noise features start to be visible, and the
orientation of the maps was specified using the eclipse ephemeris.
Overlaid on the Doppler maps are interpretations of the system
properties, adopting
km s-1 and
km s-1.
The Doppler maps for the Balmer emission lines (see the H
map in Fig. 7) have an unusual wide double-lobed structure. The H
map also shows emission attributable to the secondary star, although
this is oddly offset from the line of centres of the system. The shape
of the accretion disc and the offset of the secondary star emission may
be artefacts of the breakdown of an important assumption of Doppler
tomography: that emitting regions are optically thin.
Doppler maps of the He I emission lines show weak and diffuse emission in the region of the bright spot, which is where the accretion stream from secondary star encounters the edge of the accretion disc. The bright spot is not a major source of He I emission, but very little else is seen in the He I maps.
We have also constructed a Doppler map of the Ca II 8662 Å emission, which is the line of the calcium triplet which is least affected by night-sky emission and telluric lines. The map (Fig. 7) shows a circular accretion disc feature and clear emission arising from the secondary star. The latter emission can be seen describing a S-wave in trailed spectra (Fig. 2), and its position in the map supports our measurement of the velocity amplitude of this star.
3.5 Light curve modelling
Our light curves show deep eclipses due to obscuration of the WD and accretion disc by the secondary star. To obtain constraints on the system properties, the best dataset (2008 December) was compared to synthetic light curves created using the LCURVE code written by TRM (see Pyrzas et al. 2009). This uses grids of points to model the WD as a sphere, the secondary star using Roche geometry, a flat circular accretion disc, and an exponentially decreasing bright spot. The best fit to the observed data was obtained with a combination of the downhill simplex and Levenberg-Marquardt algorithms (Press et al. 1992). The outside-eclipse data show strong stochastic variation (termed flickering; see Bruch 1992 and Bruch 2000) arising from the mass-transfer process in SDSS J1006. We have down-weighted data outside the phase interval [ -0.07,0.08] by a factor of three, to limit their influence on the fit.
The small gap in the light curve at HJD 2454821.69 is unfortunate, as
the egress of the bright spot occurs somewhere during this time. We
find two main families of good fits to the light curve corresponding to
different WD radii and mass ratios: the first family is in the region
of
and q=0.51 and is our preferred solution. The second centres on
,
which is unphysically large, and q=0.60.
After extensive exploration of the parameter space we adopt the first
solution but increase the errorbars of the light curve parameters to
include the full range of plausible solutions which we found.
Internal parameter errors were determined by 105 Markov Chain Monte Carlo (MCMC) simulations. For these simulations we accepted a certain fraction of random jumps in parameter values and evaluated how they changed the quality of the fit. After the simulations showed reasonable convergence, the errors and covariances could be computed from examining the parameters from the accepted jumps. We rejected typically the first 10% of values to avoid a dependence on the initial parameter values. This gives a more realistic view of the parameter uncertainties compared to the values computed simply from the analytic errors alone, an aspect which is particularly important given the correlated noise due to flickering.
The results of the light curve modelling process are given in Table 4 and the best-fitting model is compared to the data in Fig. 8. The radii of the stars and accretion disc are given in units of the orbital semimajor axis, a. The uncertainties yielded by the MCMC analysis still do not fully take into account the flickering or the range of plausible solutions we found. We have increased the uncertainties to include the full range of reasonable trial solutions we found (Table 4), and regard the results as conservative. A substantial improvement will require high-speed photometry of several eclipses of the SDSS J1006 system.
Table 4: Results of the light curve modelling process.
![]() |
Figure 8: Light curve of SDSS J1006 from the 2008 December observations (points) compared to the best fit found using LCURVE (solid curve). The residuals of the fit are plotted at the base of the figure. |
Open with DEXTER |
4 The physical properties of SDSS J1006
The spectroscopic and eclipsing characteristics of SDSS J1006
allow the determination of the masses and radii of the WD and its
low-mass companion. From measurements of the times of mid-eclipse we
have obtained an accurate orbital period of 0.18591324(41) d. From
the infrared sodium doublet we have measured the velocity amplitude
km s-1. A correction for irradiation effects leads to the secondary star velocity amplitude
km s-1. From modelling the eclipse morphology of SDSS J1006 we have found an orbital inclination of
and a mass ratio of
.
Combining these results yields the masses and radii of the WD and secondary star in SDSS J1006 (Table 5). The mass of the former, 0.78
,
is higher than the average for single WDs, in agreement with previous results for CVs (Smith & Dhillon 1998; Littlefair et al. 2008). Its radius is consistent (within its large uncertainty) with the theoretical mass-radius relationship for a
K WD (Bergeron et al. 1995).
The secondary star has a mass of 0.40
and a radius of 0.47
,
which is distended compared to a normal object - a mass-radius relation based on detached eclipsing binary star systems (Southworth 2009) predicts a radius of 0.41
- but is in excellent agreement with the semi-empirical sequence for CV secondary stars constructed by Knigge (2006).
This is expected because the mass transfer timescale becomes similar to
the thermal timescale for the secondary components of CVs above the
period gap, allowing continued mass transfer to drive the star out of
thermal equilibrium.
Table 5: Physical properties of the stellar components of SDSS J1006.
4.1 Distance and white dwarf temperature
![]() |
Figure 9: Black line: the SDSS spectrum of SDSS J1006. Red line: an M3.2 template spectrum scaled to match the strengths of the spectral features of the companion star in SDSS J1006. Blue line: the residual spectrum obtained after subtracting the M-dwarf template from the spectrum of SDSS J1006. |
Open with DEXTER |
The secondary star dominates the red end of the optical spectrum of
SDSS J1006, exhibiting the strong TiO bandheads that are
characteristic of mid-to-late M dwarfs. We have obtained the
star's spectral type using the M dwarf template library that Rebassa-Mansergas et al. (2007)
assembled from SDSS spectroscopy, interpolated onto a finer grid
spanning types M0 to M9 in steps of 0.2 subtypes. These
templates were scaled and subtracted from the SDSS spectrum of
SDSS J1006 until the smoothest residual spectrum was obtained,
resulting in a spectral type of M3.2 0.2 (Fig. 9). We find a good agreement with the relationship between spectral type and mass presented by Rebassa-Mansergas et al. (2007). From the best-fit template, we calculated
,
which is the flux difference between the bands 7450-7550 Å and 7140-7190 Å, as defined by Beuermann (2006). Using the polynomial expressions of Beuermann (2006), we obtained the
surface fluxes for secondary stars with spectral types in the range M3.0 to M3.4. Taking R2=0.466
,
we find a distance of

where the uncertainty in d is dominated by that in R2.
SDSS J1006 has also been detected in the ultraviolet (UV) all-sky survey carried out by GALEX (Martin et al. 2005), from which we can estimate the WD's effective temperature,
.
The GALEX observations were obtained during the phase interval
0.117-0.124, which is outside the WD eclipse to high confidence. We
adopt
for the WD (Table 5), and d
= 676 pc. Under the assumption that the UV flux is entirely due to
the unobscured photospheric emission of the WD, the only free parameter
to reproduce the observed GALEX fluxes is then
.
The observed far-UV flux implies
K, or
K if a maximum reddening of
(Schlegel et al. 1998) is assumed (Fig. 10), with an uncertainty of
1500 K. We therefore adopt
K.
This value is only an estimate, because it is possible that the WD is
partially veiled by the accretion disc, and that the disc, bright spot
and boundary layer contribute to the observed UV fluxes. A more
reliable
measurement could be obtained from UV spectroscopy.
This
is unusually low for a dwarf nova with
h (e.g. Townsley & Gänsicke 2009; Urban & Sion 2006)
- by comparison the well-studied dwarf nova U Gem has overall
properties which are very similar to SDSS J1006 but harbours a WD
with
K (Long et al. 2006; Sion et al. 1998).
Given that SDSS J1006 is a high-inclination system, it might be
possible that the WD is partially veiled by extended structures above
the accretion disc, similar to those seen in OY Car (Horne et al. 1994). If veiling is not the cause of the low observed UV flux (no veiling is observed in U Gem), the low
implies a secular mean accretion rate of a few 10-10
yr-1 (Townsley & Gänsicke 2009), which is a factor of about three lower than in U Gem.
![]() |
Figure 10:
Spectral energy distribution of SDSS J1006,
showing the observed (solid error bars) and dereddened (
|
Open with DEXTER |
4.2 Outburst characteristics
SDSS J1006 has been observed by the Catalina Sky Survey (Drake et al. 2009), who obtained 209 unfiltered magnitude measurements between 2005 April and 2008 April.
Most of these observations found the object in the magnitude interval
17-18, but two dwarf nova outbursts have also been observed (at
JDs 2453678 and 2454259). The durations of these outbursts are not
known, but are constrained to be more than one day in the first case.
We investigated the dwarf nova outburst frequency of SDSS J1006, using Monte Carlo simulations and the times of the Catalina observations. Assuming that each outburst is observable for a 10 d period (Ak et al. 2002), we obtained a detection efficiency of 39% over the full time span of the Catalina observations. If we further assume that the outbursts occur periodically, we can obtain the detection efficiency (and thus the probable number of outbursts observed) as a function of outburst frequency. The results of this calculation are shown in Fig. 11 and favour an outburst interval in the region of 400 d. This is a very long interval for a 4-5 h period CV: Ak et al. (2002) find a mean outburst interval of 62.0 d for U Gem-type systems.
Based on its photometric and spectroscopic properties, SDSS J1006 can be classified as a dwarf nova of U Gem type, for which two dwarf nova outbursts have been detected and the outburst interval is long. Continued observations of this object would be very useful in refining its outburst frequency.
![]() |
Figure 11: Outburst detection efficiency for the Catalina Sky Survey observations of SDSS J1006 (blue lower line) and the number of outbursts which would have been detected for SDSS J1006 (black upper line), as a function of the outburst period. |
Open with DEXTER |
5 Summary and discussion
Table 6: Eclipsing cataclysmic variables for which masses and radii of one or both components has been measured accurately and precisely.
From the observations presented in this work we have discovered that
SDSS J1006 is an eclipsing CV, measured the orbital period, and
calculated the masses and radii of its component stars. This was
achieved through a parameteric model of its eclipses, combined with a
spectroscopic velocity amplitude for the secondary star corrected for
the effects of irradiation. Doppler maps of the infrared calcium
triplet reveal emission from the secondary star and are in agreement
with this
.
From the spectral characteristics of the system we have also found a WD effective temperature of
K, a secondary star spectral type of M3.2
0.2, and a distance of
pc.
A dwarf nova outburst interval of roughly 400 d agrees with the
available photometric observations of SDSS J1006.
We measured radial velocities for the WD from the H
and H
emission lines, finding a velocity amplitude
km s-1.
Despite using the diagnostic diagram approach, our RVs still have a
phase offset of 0.04 from the eclipse ephemeris. We therefore could not
assume that they represent the motion of the WD, so did not use them in
obtaining the physical properties of SDSS J1006. Notwithstanding
this, the K1 we measured turns out to be in excellent agreement with the expected WD velocity amplitude of 131.5 km s-1.
The mass of the WD is
,
and its radius is consistent with theoretical expectations. The secondary star has a mass and radius of
and
,
respectively, which is in excellent agreement with the semi-empirical sequence for CV secondary stars constructed by Knigge (2006).
The uncertainties in the system parameters are dominated by the
moderate quality of the light curve, and an improved photometric study
of this object is warranted.
In Table 6 we
have assembled a list of the component masses and radii of eclipsing
CVs with long orbital periods (greater than 3.0 h). We discount
systems with uncertain properties or whose analysis rests on
emission-line RVs (not always reliable) or mass-radius relations for
the secondary star (to avoid circular arguments). The list is
worryingly short: only 10 systems (including SDSS J1006)
satisfy our criteria, of which one is magnetic (DQ Her). The
weighted mean and standard deviation of the WD masses is
.
The masses and radii of the secondary stars display a clear negative
correlation with orbital period, as expected by our current
understanding of the evolution of CVs. The WD masses display no
correlation with orbital period or with the secondary star masses, in
agreement with studies which show that WDs in CVs do not undergo large
overall changes in mass (Knigge 2006; Prialnik & Kovetz 1995).
The reduced observational data presented in this work will be made available at the CDS (http://cdsweb.u-strasbg.fr/) and at http://www.astro.keele.ac.uk/ jkt/. We are grateful to Andrew Drake for providing the Catalina Sky Survey observations of SDSS J1006, and to the anonymous referee for a positive report. JS, TRM, BTG and CMC acknowledge financial support from STFC in the form of grant number ST/F002599/1. A.R.M. acknowledges financial support from ESO, and Gemini/Conicyt in the form of grant number 32080023. Based on observations made with the William Herschel Telescope, operated by the Isaac Newton Group, and the Nordic Optical Telescope, operated jointly by Denmark, Finland, Iceland, Norway, and Sweden, both on the island of La Palma in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. Based on observations collected at the Centro Astronómico Hispano Alemán (CAHA) at Calar Alto, Spain, operated jointly by the Max-Planck Institut für Astronomie and the Instituto de Astrofísica de Andalucia (CSIC). The following internet-based resources were used in research for this paper: the ESO Digitized Sky Survey; the NASA Astrophysics Data System; the SIMBAD database operated at CDS, Strasbourg, France; and the ar
iv scientific paper preprint service operated by Cornell University.
References
- Ak, T., Ozkan, M. T., & Mattei, J. A. 2002, A&A, 389, 478 [EDP Sciences] [CrossRef] [NASA ADS]
- Baptista, R., & Catalán, M. S. 2001, MNRAS, 324, 599 [CrossRef] [NASA ADS]
- Baptista, R., Catalán, M. S., & Costa, L. 2000, MNRAS, 316, 529 [CrossRef] [NASA ADS]
- Bergeron, P., Saumon, D., & Wesemael, F. 1995, ApJ, 443, 764 [CrossRef] [NASA ADS]
- Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [EDP Sciences] [CrossRef] [NASA ADS]
- Beuermann, K. 2006, A&A, 460, 783 [EDP Sciences] [CrossRef] [NASA ADS]
- Billington, I., Marsh, T. R., & Dhillon, V. S. 1996, MNRAS, 278, 673 [NASA ADS]
- Bruch, A. 1992, A&A, 266, 237 [NASA ADS]
- Bruch, A. 2000, A&A, 359, 998 [NASA ADS]
- Dillon, M., Gänsicke, B. T., Aungwerojwit, A., et al. 2008, MNRAS, 386, 1568 [CrossRef] [NASA ADS]
- Drake, A. J., Djorgovski, S. G., Mahabal, A., et al. 2009, ApJ, 696, 870 [CrossRef] [NASA ADS]
- Echevarría, J., de la Fuente, E., & Costero, R. 2007, AJ, 134, 262 [CrossRef] [NASA ADS]
- Feline, W. J., Dhillon, V. S., Marsh, T. R., Watson, C. A., & Littlefair, S. P. 2005, MNRAS, 364, 1158 [CrossRef] [NASA ADS]
- Fiedler, H., Barwig, H., & Mantel, K. H. 1997, A&A, 327, 173 [NASA ADS]
- Gänsicke, B. T., Fried, R. E., Hagen, H.-J., et al. 2000, A&A, 356, L79 [NASA ADS]
- Gänsicke, B. T., Araujo-Betancor, S., Hagen, H.-J., et al. 2004, A&A, 418, 265 [EDP Sciences] [CrossRef] [NASA ADS]
- Gänsicke, B. T., Rodríguez-Gil, P., Marsh, T. R., et al. 2006, MNRAS, 365, 969 [CrossRef] [NASA ADS]
- Gänsicke, B. T., Dillon, M., Southworth, J., et al. 2009, MNRAS, 397, 2170 [CrossRef] [NASA ADS]
- Hessman, F. V., Robinson, E. L., Nather, R. E., & Zhang, E.-H. 1984, ApJ, 286, 747 [CrossRef] [NASA ADS]
- Horne, K. 1986, PASP, 98, 609 [CrossRef] [NASA ADS]
- Horne, K., Welsh, W. F., & Wade, R. A. 1993, ApJ, 410, 357 [CrossRef] [NASA ADS]
- Horne, K., Marsh, T. R., Cheng, F. H., Hubeny, I., & Lanz, T. 1994, ApJ, 426, 294 [CrossRef] [NASA ADS]
- Knigge, C. 2006, MNRAS, 373, 484 [CrossRef] [NASA ADS]
- Littlefair, S. P., Dhillon, V. S., Marsh, T. R., et al. 2006, Science, 314, 1578 [CrossRef] [NASA ADS]
- Littlefair, S. P., Dhillon, V. S., Marsh, T. R., et al. 2008, MNRAS, 388, 1582 [CrossRef] [NASA ADS]
- Long, K. S., & Gilliland, R. L. 1999, ApJ, 511, 916 [CrossRef] [NASA ADS]
- Long, K. S., Brammer, G., & Froning, C. S. 2006, ApJ, 648, 541 [CrossRef] [NASA ADS]
- Marsh, T. R. 1989, PASP, 101, 1032 [CrossRef] [NASA ADS]
- Marsh, T. R., & Horne, K. 1988, MNRAS, 235, 269 [NASA ADS]
- Marsh, T. R., Horne, K., Schlegel, E. M., Honeycutt, R. K., & Kaitchuck, R. H. 1990, ApJ, 364, 637 [CrossRef] [NASA ADS]
- Martin, D. C., Fanson, J., Schiminovich, D., et al. 2005, ApJ, 619, L1 [CrossRef] [NASA ADS]
- Naylor, T., Allan, A., & Long, K. S. 2005, MNRAS, 361, 1091 [CrossRef] [NASA ADS]
- Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN 77, The art of scientific computing, 2nd edn. (Cambridge: University Press)
- Prialnik, D., & Kovetz, A. 1995, ApJ, 445, 789 [CrossRef] [NASA ADS]
- Pyrzas, S., Gänsicke, B. T., Marsh, T. R., et al. 2009, MNRAS, 394, 978 [CrossRef] [NASA ADS]
- Rebassa-Mansergas, A., Gänsicke, B. T., Rodríguez-Gil, P., Schreiber, M. R., & Koester, D. 2007, MNRAS, 382, 1377 [NASA ADS]
- Ribeiro, T., Baptista, R., Harlaftis, E. T., Dhillon, V. S., & Rutten, R. G. M. 2007, A&A, 474, 213 [EDP Sciences] [CrossRef] [NASA ADS]
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [CrossRef] [NASA ADS]
- Schneider, D. P., & Young, P. 1980, ApJ, 238, 946 [CrossRef] [NASA ADS]
- Shafter, A. W. 1983, ApJ, 267, 222 [CrossRef] [NASA ADS]
- Shafter, A. W., Szkody, P., & Thorstensen, J. R. 1986, ApJ, 308, 765 [CrossRef] [NASA ADS]
- Sion, E. M., Cheng, F. H., Szkody, P., et al. 1998, ApJ, 496, 449 [CrossRef] [NASA ADS]
- Smith, D. A., & Dhillon, V. S. 1998, MNRAS, 301, 767 [CrossRef] [NASA ADS]
- Southworth, J. 2009, MNRAS, 394, 272 [CrossRef] [NASA ADS]
- Southworth, J., Smalley, B., Maxted, P. F. L., Claret, A., & Etzel, P. B. 2005, MNRAS, 363, 529 [CrossRef] [NASA ADS]
- Southworth, J., Gänsicke, B. T., Marsh, T. R., et al. 2006, MNRAS, 373, 687 [CrossRef] [NASA ADS]
- Southworth, J., Gänsicke, B. T., Marsh, T. R., de Martino, D., & Aungwerojwit, A. 2007a, MNRAS, 378, 635 [CrossRef] [NASA ADS]
- Southworth, J., Marsh, T. R., Gänsicke, B. T., et al. 2007b, MNRAS, 382, 1145 [NASA ADS]
- Southworth, J., Gänsicke, B. T., Marsh, T. R., et al. 2008a, MNRAS, 391, 591 [CrossRef] [NASA ADS]
- Southworth, J., Townsley, D. M., & Gänsicke, B. T. 2008b, MNRAS, 388, 709 [CrossRef] [NASA ADS]
- Southworth, J., Hinse, T. C. Burgdorf, M. J., et al. 2009a, MNRAS, in press [arXiv:0907:3356]
- Southworth, J., Hinse, T. C., Jørgensen, U. G., et al. 2009b, MNRAS, 396, 1023 [CrossRef] [NASA ADS]
- Steeghs, D., Howell, S. B., Knigge, C., et al. 2007, ApJ, 667, 442 [CrossRef] [NASA ADS]
- Szkody, P., Henden, A., Mannikko, L., et al. 2007, AJ, 134, 185 [CrossRef] [NASA ADS]
- Szkody, P., Anderson, S. F., Hayden, M., et al. 2009, AJ, 137, 4011 [CrossRef] [NASA ADS]
- Thoroughgood, T. D., Dhillon, V. S., Watson, C. A., et al. 2004, MNRAS, 353, 1135 [CrossRef] [NASA ADS]
- Thoroughgood, T. D., Dhillon, V. S., Steeghs, D., et al. 2005, MNRAS, 357, 881 [CrossRef] [NASA ADS]
- Thorstensen, J. R. 2000, PASP, 112, 1269 [CrossRef] [NASA ADS]
- Townsley, D. M., & Gänsicke, B. T. 2009, ApJ, 693, 1007 [CrossRef] [NASA ADS]
- Unda-Sanzana, E., Marsh, T. R., & Morales-Rueda, L. 2006, MNRAS, 369, 805 [CrossRef] [NASA ADS]
- Urban, J. A., & Sion, E. M. 2006, ApJ, 642, 1029 [CrossRef] [NASA ADS]
- Wade, R. A., & Horne, K. 1988, ApJ, 324, 411 [CrossRef] [NASA ADS]
- Wood, M. A., Robertson, J. R., Simpson, J. C., et al. 2005, ApJ, 634, 570 [CrossRef] [NASA ADS]
Footnotes
- ...
SDSS J100658.40+233724.4
- Spectra are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/507/929
- ... PAMELA
- PAMELA and MOLLY were written by TRM and can be obtained from http://www.warwick.ac.uk/go/trmarsh
- ... STARLINK
- The Starlink software and documentation can be obtained from http://starlink.jach.hawaii.edu/
- ... MIDAS
- http://www.eso.org/projects/esomidas/
- ... Lupton
- The ugriz-BVRI transformation equations are attributed to ``Lupton (2005)'' but appear to be unpublished. They can be found at http://www.sdss.org/dr6/algorithms/sdssUBVRITransform.html
- ... SBOP
- Spectroscopic Binary Orbit Program, written by Etzel, http://mintaka.sdsu.edu/faculty/etzel/
- ... April
- See http://nesssi.cacr.caltech.edu/catalina/20050301/SDSSCV.html#table77
All Tables
Table 1: Log of the observations presented in this work. The mean magnitudes are calculated excluding observations taken during eclipse.
Table 2: Times of eclipse for SDSS J1006 and the residuals with respect to the linear ephemeris given in Sect. 3.1.
Table 3: Best-fitting spectroscopic orbits found using SBOP.
Table 4: Results of the light curve modelling process.
Table 5: Physical properties of the stellar components of SDSS J1006.
Table 6: Eclipsing cataclysmic variables for which masses and radii of one or both components has been measured accurately and precisely.
All Figures
![]() |
Figure 1: Flux-calibrated average spectrum of SDSS J1006. Data from the blue arm of ISIS are shown in the left panel, and from the red arm in the right panel. The most prominent emission and absorption features are labelled. |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Greyscale plot of the continuum-normalised
and phase-binned trailed spectra of SDSS J1006.
From left to right the panels show H |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Plot of the four light curves obtained covering eclipses of SDSS J1006. The eclipse midpoints have been aligned on the panels. |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Diagnostic diagram showing the
variation of the best-fitting spectroscopic orbital parameters
for RVs measured with a range of separations using the double
Gaussian function.
|
Open with DEXTER | |
In the text |
![]() |
Figure 5:
The measured RVs (circles) and the spectroscopic
orbits fitted to them (solid curves). The emission-line RVs (filled circles)
were calculated using
|
Open with DEXTER | |
In the text |
![]() |
Figure 6: Variation of the equivalent widths of the Na doublet and TiO molecular band with orbital phase. The spectra were combined into 20 phase bins prior to measurement. The black filled circles represent measurements outside primary and secondary eclipses, and the grey open circles those within eclipse. The wavelength intervals over which the equivalent widths were measured were 8180-8200 Å for Na and 7090-7345 Å for TiO. |
Open with DEXTER | |
In the text |
![]() |
Figure 7:
Doppler maps of H |
Open with DEXTER | |
In the text |
![]() |
Figure 8: Light curve of SDSS J1006 from the 2008 December observations (points) compared to the best fit found using LCURVE (solid curve). The residuals of the fit are plotted at the base of the figure. |
Open with DEXTER | |
In the text |
![]() |
Figure 9: Black line: the SDSS spectrum of SDSS J1006. Red line: an M3.2 template spectrum scaled to match the strengths of the spectral features of the companion star in SDSS J1006. Blue line: the residual spectrum obtained after subtracting the M-dwarf template from the spectrum of SDSS J1006. |
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Spectral energy distribution of SDSS J1006,
showing the observed (solid error bars) and dereddened (
|
Open with DEXTER | |
In the text |
![]() |
Figure 11: Outburst detection efficiency for the Catalina Sky Survey observations of SDSS J1006 (blue lower line) and the number of outbursts which would have been detected for SDSS J1006 (black upper line), as a function of the outburst period. |
Open with DEXTER | |
In the text |
Copyright ESO 2009
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.