Issue |
A&A
Volume 532, August 2011
|
|
---|---|---|
Article Number | A7 | |
Number of page(s) | 8 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201117078 | |
Published online | 12 July 2011 |
The cool magnetic DAZ white dwarf NLTT 10480⋆
Astronomický ústav, Akademie věd České republiky, Fričova 298, 251 65 Ondřejov, Czech Republic
e-mail: kawka@sunstel.asu.cas.cz; vennes@sunstel.asu.cas.cz
Received: 14 April 2011
Accepted: 1 June 2011
We have identified a new cool magnetic white dwarf in the New Luyten Two-Tenths (NLTT) catalogue. The high proper-motion star NLTT 10480 (μ = 0.5″ yr-1) shows weak Zeeman-split lines of calcium as well as characteristic Hα and β Zeeman triplets. Using VLT X-shooter spectra, we measured a surface-averaged magnetic field BS ~ 0.5 MG. The relative intensity of the π and σ components of the calcium and hydrogen lines imply a high inclination (i ≳ 60°). The optical-to-infrared V − J colour index and the Ca I/Ca II ionization balance indicate a temperature between 4900 and 5200 K, while the Balmer line profiles favour a higher temperature of 5400 K. The discrepancy is potentially resolved by increasing the metallicity to 0.03 × solar, hence increasing the electron pressure. However, the measured calcium abundance and abundance upper limits for other elements (Na, Al, Si, and Fe) imply a low photospheric metallicity ≲ 10-4 solar. Assuming diffusion steady-state, a calcium accretion rate of log Ṁ(g s-1) = 5.6 ± 0.3 is required to sustain a calcium abundance of log n(Ca)/n(H) = −10.30 ± 0.30 in the white dwarf atmosphere. We examine the implications of this discovery for the incidence of planetary debris and weak magnetic fields in cool white dwarf stars.
Key words: white dwarfs / stars: individual: NLTT 10480 / stars: atmospheres / stars: abundances
© ESO, 2011
1. Introduction
NLTT 10480 (LHS 5070, LP 887-66) is a high proper-motion star (Luyten 1979, 1980) that was also listed as a white dwarf candidate in Luyten’s white dwarf catalogue (Luyten 1977). Using an optical-infrared reduced proper-motion diagram (Salim & Gould 2002) and additional colourimetric criteria (Kawka et al. 2004), we selected NLTT 10480 for spectroscopic observations to investigate stellar properties such as the effective temperature (Teff), surface gravity (log g), chemical composition, and magnetic field strength, and to constrain its cooling age and mass.
The chemical composition of cool white dwarfs shows great diversity. Heavy elements, particularly calcium, are detected in close to a quarter of cool, hydrogen-rich (DA) white dwarfs, but with abundances well below solar (Zuckerman et al. 2003). The abundance of heavy elements decreases with white dwarf cooling ages mainly because of the increasing depth of the mixed convective layers in aging white dwarfs (Paquette et al. 1986; Koester 2009). Few examples of very cool (Teff ≲ 5000 K) polluted white dwarfs (DAZ) are known, such as G 77-50 (WD 0322 − 019, Hintzen & Strittmatter 1974; Sion et al. 1990), which is also harbouring a weak magnetic field (Farihi et al. 2011). New high-dispersion and high signal-to-noise ratio spectroscopic observations of faint high proper-motion stars are likely to contribute new objects to the current sample.
Heavy elements in the atmosphere of cool white dwarfs are almost certainly accreted from their immediate environment. Kilic et al. (2006) and Farihi et al. (2009) reported infrared observations of a sample of cool white dwarfs contaminated with heavy elements, and the authors noted an infrared-excess incidence of ≈ 10 − 20%. This excess was attributed to debris discs of temperatures ranging from a few hundred degrees to over 1000 K. Some cool DAZ white dwarfs, such as G 174 − 74 (WD 0245+541), do not show an infrared excess (Debes et al. 2007), and the presence of heavy elements in cool white dwarfs with ages in excess of several billion years suggests, instead, the effect of episodic accretion from small asteroids rather than from a stable debris disc (Jura 2008). Therefore, the identification of new cool DAZ white dwarfs is of interest to constrain the phenomenon.
We present a first report on a programme aimed at identifying and characterizing new DAZ white dwarfs in the New Luyten Two-Tenths (NLTT) catalogue. Section 2 describes observations obtained at the European Southern Observatories (ESO) using the New Technology Telescope (3.6-m) and the Very Large Telescopes (VLTs). Section 3 presents our model atmosphere analysis including details of the model structures (Sect. 3.1), heavy element line opacities, and Zeeman effect on line profiles (Sect. 3.2), and measurements (Sect. 3.3). We summarize and discuss some implications of our results in Sect. 4.
2. Observations
We first observed NLTT 10480 with the focal reducer and low-dispersion spectrograph (FORS1) attached to the 8 m UT2 (Kueyen) at Paranal Observatory as part of our spectro-polarimetric survey of white dwarfs. The purpose of the survey was to search for white dwarfs with weak magnetic fields. We used the 600B grism combined with a slit-width of 1 arcsec that provided a resolution of 6.0 Å. The spectra covered the range between 3780 and 6180 Å. The observations were conducted on UT 2007 Nov. 1 and consisted of a sequence of two consecutive exposures with an exposure time of 1360 s each. In the first exposure the Wollaston prism is rotated to −45° and it is followed by a second exposure with the Wollaston prism rotated to +45° from which we extracted the flux and circular polarization spectra.
Since the FORS1 spectra did not include Hα, we obtained two additional low-dispersion spectra with the ESO Faint Object Spectrograph and Camera (EFOSC2) attached to the New Technology Telescope (NTT) at La Silla Observatory on UT 2008 Oct. 23. We used Grism 11 which has 300 lines per mm and a blaze wavelength of 4000 Å. The slit-width was set to 1 arcsec, which resulted in a spectral resolution of ~14 Å. The exposure time of each spectrum was 1500 s. Observations were carried out at the parallactic angle and were flux calibrated with the flux standard Feige 110.
Figure 1 compares the low-resolution spectrum of Ca II lines obtained with EFOSC2 to the higher resolution spectrum of FORS1. The flux and circular polarization spectra clearly show the σ components, proving NLTT 10480 to be magnetic.
![]() |
Fig. 1 Flux spectra of the Ca II H&K lines in NLTT 10480 obtained with EFOSC2 (top) and FORS1 (middle), and circular polarization (bottom) spectrum obtained with FORS1. The dotted lines mark the position of the σ components and the dashed lines mark the π components (see Sect. 3.2.3). |
Finally, we obtained a set of echelle spectra of NLTT 10480 with the X-shooter spectrograph attached to the UT2 (Kueyen) at Paranal Observatory. The spectra were obtained at four epochs (UT 2010 Dec. 10, 2011 Jan. 9 and 26, and March 3). The 2011 Mar. 3 spectrum was of much poorer quality and was excluded from the analysis. The spectra were obtained with the slit-width set to 0.5, 0.9, and 0.6 arcsec for the UVB, VIS, and NIR arms, respectively. This set-up delivered a resolving power of 9100 for UVB, 8800 for VIS, and 6200 for NIR. The exposure time for each exposure was 2400 s.
Table 1 summarizes our spectroscopic observations. The signal-to-noise ratio in the summed X-shooter spectrum is S/N ≈ 53 per binned-pixel at 4000 Å, while it is ≈ 120 in the summed FORS1 spectrum and ≈65 in the summed EFOSC2 spectrum at the same wavelength.
Log of spectroscopic observations.
We searched for photometric measurements using VizieR. The Two Micron All Sky Survey (2MASS) listed infrared JHK magnitudes, but only the J magnitude was of acceptable quality. We also used the acquisition images from the X-shooter observations to estimate a V magnitude for NLTT 10480. We used 11 acquisition images of Feige 110 obtained between UT 2010 Dec. 10 and 2011 Jan. 1, and set the zero point. Next, we determined an average V magnitude for NLTT 10480 using our five acquisition images. We employed the atmospheric extinction table of Patat et al. (2011). Table 2 lists the photometric measurements.
Photometry and astrometry.
The SPM Catalog provides a photographic B magnitude (18.08 ± 0.16). However, due to large uncertainties in the SPM and other photographic magnitudes, these are not as useful as our V magnitude and the 2MASS J magnitudes. The 2MASS K magnitude is unreliable and the H magnitude is uncertain. The index B − V = 0.59 ± 0.17 loosely constrains the temperature in the range 5100 ≲ Teff ≲ 6700 K, but the more precise index V − J = 1.49 ± 0.09 implies that NLTT 10480 is a cool white dwarf with Teff ≈ 5050 ± 150 K (see Table 3 and Sect. 3.1).
NLTT 10480 is also characterized by a large proper-motion of ~0.5″ yr-1 (Table 2). The kinematical properties were determined using the proper-motion and the X-shooter radial velocity measurement (see Sect. 4).
3. Analysis
3.1. Model atmospheres and spectral syntheses
We calculated a grid of model atmospheres for cool hydrogen-rich white dwarfs. The models are in convective and radiative equilibrium with the total flux converged to better than 0.01% in all layers. The model grid covers the effective temperature range 4900 ≤ Teff ≤ 6000 K in 100 K steps and the surface gravity range 7.5 ≤ log g ≤ 8.5 in steps of 0.25 dex. The adopted treatment of the convective energy transport and pressure ionization effects are described in Kawka & Vennes (2006), but further improvements to the models will be described in a forthcoming publication. All relevant species (H, H+, H2, H, H
) are included in the statistical equilibrium equation. We employed the H
partition function of Neale & Tennyson (1995). Electrons contributed by identifiable trace elements (e.g., calcium) are also included in the charge conservation equation, although the ionization of hydrogen atoms and molecules dominate the electron budget.
The model atmospheres include opacities caused by H bound-bound, bound-free and free-free transitions, H− bound-free and free-free transitions, H2 − H2 collision-induced absorption (CIA, Borysow et al. 1997, 2001), and the H2-H and H-H collision-induced absorptions in the far Lyα wing (see Kowalski & Saumon 2006) using opacity tables from Rohrmann et al. (2011). Finally, the H2 and H Rayleigh scattering are included along with electron scattering.
Synthetic colours as well as detailed hydrogen and heavy element line profiles are computed using the model structures. Table 3 lists some photometric properties of the cool models. The colour indices at shorter wavelengths are effected by the Lyα collision-induced absorptions. Allard et al. (2008) calculated new Hα line profiles including self-broadening effects and found a half-width at half maximum (HWHM) ≈ 40% larger than estimated by Ali & Griem (1965, 1966) and comparable to the calculations of Barklem et al. (2000a). We adopted the Hα HWHM from Allard et al. (2008) and the Hβ, and Hγ cross-sections and velocity parameters, converted into HWHM, from Barklem et al. (2000a). We describe the heavy element line profiles in the following section.
Selected synthetic colours.
3.2. Neutral and ionized line profiles
3.2.1. Line broadening
The dominant broadening mechanism is collision with hydrogen atoms. We employed the coefficients of Barklem et al. (2000b), where the full-width at half-maximum (FWHM) of the Lorentzian profiles is given by (1)where log Γ = −7.562 for Ca Iλ4226 and − 7.76 for Ca II H&K at T = 10 000 K, and α = 0.238 for Ca Iλ4226 and 0.223 for Ca II H&K. Although the adopted broadening parameters do not include the effect of hydrogen molecules, these provide > 50% of the gas pressure in some layers. Following the approximate treatment of Kurucz & Avrett (1981), the broadening parameter Γ ∝ nH + 0.85 n(H2), and, therefore, hydrogen molecules may contribute to the total line width. We found that neither Ca I or Ca II equivalent widths are significantly affected in models at Teff = 5400 K, but abundances inferred from Ca Iλ4226 may be underestimated by a factor of ~2 in models at 4900 K.
3.2.2. Zeeman effect
Fundamentals of stellar line formation in the presence of a magnetic field are described by Unno (1956) and Martin & Wickramasinghe (1981), who also describe significant magneto-optical effects. In particular, Unno (1956) showed that the effect of field inclination with respect to the line-of-sight on the relative intensity of σ and π components reaches a maximum at an angle of 55°. Moreover, Martin & Wickramasinghe (1981) showed that taking into account magneto-optical effects may enhance the depth of the π components upon certain conditions.
Kemic (1975) studied the quadratic Zeeman effect for the Ca II H and K lines and showed that the linear Zeeman effect at fields of 15 MG is still dominant.
The Ca II H and K lines are the result of transitions between the ground state with , and the excited states with
and
, respectively, where J is the total angular momentum, L is the orbital angular momentum and S is the spin angular momentum. The Ca Iλ4226 line is the result of transitions between the ground state with J,L,S = 0,0,0 and the excited state with J,L,S = 1,1,0.
The levels are split by a magnetic field into 2J + 1 components defined by the magnetic quantum number m = −J,...,J: (2)where λ is the wavelength in Å, B is the magnetic field in MG, e is the electron charge, me is the electron rest mass and c is the speed of light. The Landé factor and the magnetic quantum number of the upper and lower levels are given by gu,mu and gl,ml, respectively. Landé factors for calcium and other elements except iron were calculated assuming LS coupling:
(3)For the 4S1/2 state of Ca II, the experimentally determined g factor of 2.00225664 (Tommaseo et al. 2003) agrees with the theoretically calculated factor based on the LS coupling scheme. Similarly, for the 3p 2P3/2 state of AlI, the experimentally determined g factor of 1.33474 ± 0.00005 (Martin et al. 1968) agrees with the LS coupling scheme calculated factor. For iron, the Landé factors were obtained from the Vienna Atomic Line Database (VALD)1.
The permitted transitions are defined by Δm = 0, ± 1, where Δm = 0 defines the π components and Δm = ± 1 the σ components.
3.2.3. Opacity and line intensities
The relative intensities of the Zeeman components are computed following Condon & Shortley (1963). They found for the transitions where ΔJ = 0 and Δm = 0 (π components) (4)or where Δm = ± 1 (σ components)
(5)Similarly, they found for transitions where J → J + 1 and for the π components
(6)and for the σ components
(7)Finally, for transitions where J → J − 1, the relative intensities are for the π components
(8)and for the σ components
(9)Table 4 lists the total angular momenta, Landé factors, magnetic quantum numbers for the various transitions of calcium in a magnetic field of 0.513 MG. The relative intensities for the different components are also provided.
Zeeman splitting at B = 0.513 MG ( Ca I and Ca II).
![]() |
Fig. 2 Linear Zeeman splitting of Ca II (H and K) (top), Ca I 4226 Å (bottom left), and Hα (bottom right) lines at B = 0.519 MG compared to the observed Zeeman splitting in the spectrum of NLTT 10480. The dotted lines indicate the σ components and dashed lines indicate the π components. |
The observed relative intensities of the π and σ components will also vary as a function of the angle between the magnetic field lines and the line of sight (Unno 1956; Martin & Wickramasinghe 1981). The monochromatic opacity of the π components depends on the angle ψ between the magnetic field axis and the line of sight: (10)Similarly, the monochromatic opacity of the σ components vary as
(11)where the opacities χπ(90°) and χσ(90°) are given at ψ = 90° and, as shown earlier, are calculated following Condon & Shortley (1963). For example, for hydrogen Balmer lines or Ca Iλ4226 we have that
(12)Figure 2 shows the observed linear Zeeman splitting of Ca II H&K, Ca Iλ4226 Å, and Hα lines obtained with the X-shooter spectrograph and the predicted line positions as a function of the magnetic field strength. Interestingly, the observed π components for both Ca II lines appear weaker than the σ components showing the effect of inclination of the field with respect to the line-of-sight. No significant variation in the line positions or intensities were noted between the three usable exposures.
![]() |
Fig. 3 Balmer lines (Hα to Hγ) and best-fitting model (Teff = 5410 K, log g = 8.0) showing prominent Zeeman-split Hα line (left panel), and (right panel) χ2 contours at 66, 90, and 99% confidence (dashed lines) in the (Teff, log g) plane. The measured calcium abundance from Ca II lines (log n(Ca)/n(H)) is drawn (full lines) as a function of the assumed temperature and surface gravity and labelled with the logarithm of the abundance. The thick line shows the range of temperature and surface gravity where the abundance measurements based on Ca I or Ca II are consistent. |
We fitted the ratio of the Ca II π and σ line equivalent widths varying the inclination and determined ψ = 60 ± 3°. However, including magneto-optical effects would likely decrease the measured inclination.
3.3. Magnetic field and atmospheric parameters
We determined the averaged surface magnetic field of NLTT 10480 using Hα and the calcium lines. We first measured the centroids of the Zeeman components of Ca II, Ca I and Hα, where the spectrum was already adjusted to the solar system barycentre. We then fitted these lines to the predicted line positions by varying the magnetic field strength (assumed constant over the surface) and the velocity and by minimizing the χ2. For calcium, we determined a surface-averaged magnetic field of Bs = 0.513 ± 0.005 MG with a velocity of 3 = 50.8 ± 3.0 km s-1. Using Hα, we obtained Bs = 0.525 ± 0.005 MG and 3 = 51.8 ± 3.0 km s-1. The measurements are consistent within 2σ and taking the average we found Bs = 0.519 ± 0.004 and 3 = 51.3 ± 2 km s-1. Subtracting the gravitational redshift ( km s-1) of the white dwarf from the measured velocity results in the actual line of sight velocity of the star
km s-1.
In principle, the effective temperature and surface gravity may be constrained simultaneously by fitting the Balmer line profiles. The strong Hα and weaker β show the effect of Zeeman line splitting, while Hγ is extremely weak. Figure 3 shows the best-fitting model to the Balmer lines that includes an approximate treatment of line opacities in a weak magnetic field. The 99% confidence contour corresponds to an uncertainty of 200 K in the temperature, and 0.5 dex in the logarithm of the surface gravity. The errors are statistical only and do not take into account possible systematic errors in the modelling of the line profiles.
3.3.1. Abundance of heavy elements
We determined the calcium abundance by fitting the Ca II line profiles using χ2 minimization techniques where we varied the calcium abundance at each grid point in the (Teff,log g) plane while fixing the magnetic field strength at 0.519 MG (Sect. 3.3). Figure 3 shows the resulting abundance map: the calcium abundance and its uncertainty depend on the adopted atmospheric parameters Teff and log g.
Figure 4 shows best-fitting models to the calcium lines at Teff = 5400 (log g = 8) and at a lower temperature of 4900 K (log g = 8). The calcium abundance varies between log n(Ca)/n(H) = −10.3 (highest temperature) and − 10.45 (lowest temperature), but the calcium ionization balance favours a lower temperature than estimated with the Balmer lines (Fig. 3). The discrepancy is partially resolved by increasing the metallicity of the atmosphere. Using a set of heavy elements with low first-ionization potentials and higher solar abundances including C, N, and O, then Na to Si, and, finally, K to Cu, we found that a high metallicity of 0.03 × solar would restore the calcium ionization balance with an abundance of log n(Ca)/n(H) = −10.2. However, a much lower abundance of heavy elements is present in the atmosphere of NLTT 10480.
Figure 5 shows the predicted location of the strongest lines of Fe I, Si I, Al I and Na I in the X-shooter spectrum. We calculated the position of the Zeeman-split lines for these elements assuming a magnetic field of 0.519 MG. For aluminium, some weak lines appear to match the predicted positions. However, the putative Al Iλ3967.8532 component should be accompanied by stronger λ3957.7353 and λ3960.2458 components that are not clearly identified. Another possible identification for this feature is an interstellar Ca H line, although Ca K is not detected in the spectrum. Spectra with higher signal-to-noise ratios and better resolution are needed to clarify this identification.
We estimated sodium, aluminium, silicon and iron abundance upper limits using the spectral ranges covered in Fig. 5. We found log n(Na)/n(H), log n(Al)/n(H), and log n(Fe)/n(H) ≲ − 9.3, while log n(Si)/n(H) ≲ − 8.7. The abundances relative to solar range from 2 × 10-5 to 3 × 10-4 times solar, or a few orders of magnitude below the level required to significantly increase the electron density and alter the calcium ionization balance.
![]() |
Fig. 4 X-shooter spectrum of the Ca I (bottom) and Ca II (H & K: middle and top, respectively) lines and best-fitting models at (log g = 8) Teff = 5400 K and log n(Ca)/n(H) = −10.3 (thick lines), and Teff = 4900 K and log n(Ca)/n(H) = −10.45 (thin lines). The wavelength scale is normalized to Δλ1 ≡ 4.67 × 10-7λ2B, where B = 0.513 MG. Therefore, the line positions match the factors glml − gumu listed in Table 4 (vertical dashed lines). All models assume ψ = 60°. |
![]() |
Fig. 5 Spectra in the vicinity of the strongest Fe I (3734.8638, 3749.4851, and 3758.2330 Å), Si I (3905.523 Å), Al I (3944.0060 and 3961.5200 Å), and Na I (5889.951 and 5895.924 Å) lines compared to the predicted Zeeman-split line positions assuming a magnetic field of 0.519 MG. |
4. Summary and discussion
We found that the high proper-motion star NLTT 10480 is a rare example of cool white dwarfs with trace heavy elements and a weak magnetic field revealed in both the Hα circular polarization spectrum and Zeeman line-splitting. Other examples of this class of white dwarfs are the DZ white dwarfs LHS 2534 (Reid et al. 2001) and G 165-7 (Dufour et al. 2006), the DAZ white dwarfs G 77-50 (Farihi et al. 2011) and possibly LTT 8381 (Koester et al. 2009). Based on independent diagnostics (Table 5), we estimated a temperature of Teff = 5200 ± 200 K and a surface gravity close to log g = 8. However, we noted systematic differences in temperature measurements based on the calcium ionization ratio, the colour index, and the Balmer line profiles amounting to ~400 K. The weaker Ca I lines favour a lower temperature than estimated using Balmer lines alone. The temperature measured with the V − J colour index also favours a lower temperature. We found that increasing the heavy-element contribution to the electron density helps restore the calcium ionization balance, but we also found that the required abundance exceeds upper limits on the abundance of Na, Al, Si, and Fe by a few orders of magnitude. We are left with the possibility that subtle effects on line formation (broadening parameters, magnetic-optical, ...) caused by the magnetic field may influence temperature measurements based on Balmer line profiles.
Although our modelling of the hydrogen line profiles takes into account the effect of inclination, we neglected the magneto-optical effects and only approximated the full solution of the radiative transfer equations that ought to include all Stokes parameters (see Martin & Wickramasinghe 1981). The effect of this approximation on the determination of the stellar parameters using Balmer line profile fitting may well amount to a few hundred degrees.
Fortunately, because ionized calcium is the dominant species, the abundance of calcium based on Ca II lines is not sensitive to temperature. On the other hand, as we have demonstrated, it does show a mild dependence on surface gravity because a higher electron pressure favours neutral calcium.
Properties of NLTT 10480.
Adopting conservative error bars for the temperature and surface gravity we calculated an absolute V magnitude using the mass-radius relations of Benvenuto & Althaus (1999). The distance modulus implies a photometric distance
pc. The object is relatively old with a cooling age tcool = 3.5 − 8.0 Gyr, but with an uncertain mass (0.35 − 0.91 M⊙). We determined the Galactic velocity vector UVW (Table 5) using our distance estimate and radial velocity measurement (Sect. 3.3), and published proper-motion. We employed the algorithm of Johnson & Soderblom (1987). The kinematics imply membership to the old thin disk (Sion et al. 1988) consistent with the upper range of our age estimate.
The acquisition of broadband UBV and JHK photometry and of a parallax measurement should help determine the atmospheric parameters more precisely. The stellar radius, hence surface gravity measurement would be improved with a parallax measurement. Accordingly, the error on the calcium abundance measurement would be reduced. Accurate JHK photometry would also allow us to investigate possible infrared excess and the presence of a debris disc.
The DAZ white dwarfs G 77 − 50 and G 174 − 74 are part of a survey including the coolest known DAZ white dwarfs (Zuckerman et al. 2003). In cool convective white dwarfs, heavy elements diffuse below the mixed convective layers, and their presence in white dwarf atmospheres is transitory. Koester & Wilken (2006) estimated the diffusion time-scale for various heavy elements. The accretion rate required to sustain a given mass fraction X in the atmosphere is given by (13)where Mcvz is the mass of the convection zone and τ the diffusion time-scale at the bottom of the convection zone, where diffusion is allowed to take place. The ratio X/Xacc is the ratio of the measured mass-fraction to the accreted mass fraction. Therefore, the mass accretion rate of any particular element is
(14)Adopting, in the appropriate temperature range, a value for the slow-varying ratio Mcvz to diffusion time scale of log (Mcvz/τCa) ≈ − 11.5 in units of M⊙ yr-1 or = 14.3 in units of g s-1 (Koester & Wilken 2006), we estimated the mass accretion rate of calcium (in g s-1) to be
(15)where XCa ≈ [ACa n(Ca) ] / [AH n(H) ] = 2 × 10-9, and ACa and AH are the atomic weights. Assuming calcium is accreted as part of a solar-composition flow, the total mass accretion rate (including hydrogen) is 5 × 109 g s-1 or 8 × 10-17 M⊙ yr-1.
Our measured error on the calcium abundance of ± 0.3 dex translates into a similar error on the calculated accretion rate onto the white dwarf surface. The true error may well be much larger. Koester (2009) considers that the application of the mixing-length theory to the structure of convection zones may underestimate the mass of the mixed layers by orders of magnitude. The effect of “under”-shooting below the convection zone may affect diffusion time-scale estimates. Therefore, the precision claimed in measuring abundance of parent bodies may be over-estimated.
Farihi et al. (2011) propose that the model of Potter & Tout (2010)2 for the presence of a magnetic field in post-common envelope (CE) binaries could also be applied to CE episodes with planetary rather than stellar secondary components. Whether the magnetic field is acquired during such a process, or whether it is a fossil field cannot be ascertained for individual objects but rather from population studies (see, e.g., Kawka & Vennes 2004; Kawka et al. 2007; Wickramasinghe & Ferrario 2005). In the case of NLTT 10480, which is old (> 3.5 Gyr) with relatively short diffusion time-scales, the present-day metallicity is not linked to the CE event that potentially generated the magnetic field, but more likely to a recent accretion event.
A low incidence of planetary systems would imply a low incidence of weak magnetic fields (B ≲ 1 MG). Current data indicate a low incidence of weak magnetic fields, and Kawka et al. (2007) found that 6 out of 53 local white dwarfs (d ≤ 20 pc) observed with sufficient accuracy to unveil fields weaker than 1 MG were found to harbour such a low field. Moreover, Kawka & Vennes (2004) found that low-field white dwarfs lack progenitors, a gap that could be filled with the CE-mechanism.
Early results from the Kepler survey also indicate a low incidence of very large planets in short-period orbits (P ≲ 50 days, or a ≲ 0.2 AU) that are likely to trigger the field-generating CE events postulated by Farihi et al. (2011) based on the model of Potter & Tout (2010). Planets with sizes ranging from 8 to 32 R ⊕ (= 0.7 − 2.9 RJ) may occur within 0.25 AU of solar-type stars with a frequency of 1.3% (Howard et al. 2011). Borucki et al. (2011) quote a similar fraction for Jupiter-sized or larger planets within 0.2 AU based on the first data set from Kepler. Both studies noted a declining occurrence with increasing separation. Overall, large planets that are likely to participate in a CE phase may surround a few percents of white dwarf progenitors and generate, as observed, a similar percentage of low-field white dwarfs.
See also Nordhaus et al. (2011).
Acknowledgments
S.V. and A.K. are supported by GA AV grant numbers IAA300030908 and IAA301630901, respectively, and by GA ČR grant number P209/10/0967. A.K. also acknowledges support from the Centre for Theoretical Astrophysics (LC06014). We thank the anonymous referee and R. Napiwotzki for several comments that improved the paper. This research has made use of the VizieR catalogue access tool, CDS, Strasbourg, France. This publication makes use of data products 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.
References
- Ali, A. W., & Griem, H. R. 1965, Phys. Rev., 140, 1044 [NASA ADS] [CrossRef] [Google Scholar]
- Ali, A. W., & Griem, H. R. 1966, Phys. Rev., 144, 366 [NASA ADS] [CrossRef] [Google Scholar]
- Allard, N. F., Kielkopf, J. F., Cayrel, R., & van’t Veer-Menneret, C. 2008, A&A, 480, 581 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barklem, P. S., Piskunov, N., & O’Mara, B. J. 2000a, A&A, 363, 1091 [NASA ADS] [Google Scholar]
- Barklem, P. S., Piskunov, N., & O’Mara, B. J. 2000b, A&AS, 142, 467 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
- Benvenuto, O. G., & Althaus, L. G. 1999, MNRAS, 303, 30 [NASA ADS] [CrossRef] [Google Scholar]
- Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 728, 117 [NASA ADS] [CrossRef] [Google Scholar]
- Borysow, A., Jørgensen, U. G., & Zheng, C. 1997, A&A, 324, 185 [NASA ADS] [Google Scholar]
- Borysow, A., Jørgensen, U. G., & Fu, Y. 2001, JQSRT, 68, 235 [NASA ADS] [CrossRef] [Google Scholar]
- Condon, E. U., & Shortley, G. H. 1963 (Cambridge University Press) [Google Scholar]
- Debes, J. H., Sigurdsson, S., & Hansen, B. 2007, AJ, 134, 1662 [NASA ADS] [CrossRef] [Google Scholar]
- Dufour, P., Bergeron, P., Schmidt, G. D., et al. 2006, ApJ, 651, 1112 [NASA ADS] [CrossRef] [Google Scholar]
- Farihi, J., Jura, M., & Zuckerman, B. 2009, ApJ, 694, 805 [NASA ADS] [CrossRef] [Google Scholar]
- Farihi, J., Dufour, P., Napiwotzki, R., & Koester, D. 2011, MNRAS, 413, 2559 [NASA ADS] [CrossRef] [Google Scholar]
- Hintzen, P., & Strittmatter, P. A. 1974, ApJ, 193, L111 [NASA ADS] [CrossRef] [Google Scholar]
- Howard, A. W., Marcy, G. W., Bryson, S. T., et al. 2011, ApJ, submitted [arXiv:1103.2541] [Google Scholar]
- Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [NASA ADS] [CrossRef] [Google Scholar]
- Jura, M. 2008, AJ, 135, 1785 [NASA ADS] [CrossRef] [Google Scholar]
- Kawka, A., & Vennes, S. 2004, in A-Star Puzzle (Cambridge University Press), ed. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, Proc. IAU Symp., 224, 879 [Google Scholar]
- Kawka, A., & Vennes, S. 2006, ApJ, 643, 402 [NASA ADS] [CrossRef] [Google Scholar]
- Kawka, A., Vennes, S., & Thorstensen, J. R. 2004, AJ, 127, 1702 [NASA ADS] [CrossRef] [Google Scholar]
- Kawka, A., Vennes, S., Schmidt, G. D., Wickramasinghe, D. T., & Koch, R. 2007, ApJ, 654, 499 [NASA ADS] [CrossRef] [Google Scholar]
- Kemic, S. B. 1975, Ap&SS, 36, 459 [NASA ADS] [CrossRef] [Google Scholar]
- Kilic, M., von Hippel, T., Leggett, S. K., & Winget, D. E. 2006, ApJ, 646, 474 [NASA ADS] [CrossRef] [Google Scholar]
- Koester, D. 2009, A&A, 498, 517 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Koester, D., & Wilken, D. 2006, A&A, 453, 1051 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Koester, D., Voss, B., Napiwotzki, R., et al. 2009, A&A, 505, 441 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kowalski, P. M., & Saumon, D. 2006, ApJ, 651, L137 [Google Scholar]
- Kupka, F., Ryabchikova T. A., Piskunov N. E., Stempels H. C., & Weiss, W. W., 2000, Baltic Astron., 9, 590 [Google Scholar]
- Kurucz, R. L., & Avrett, E. H. 1981, SAO Spec. Rep., 391 [Google Scholar]
- Luyten, W. J. 1977, White Dwarfs II (Minneapolis: Univ. Minnesota Press) [Google Scholar]
-
Luyten, W. J. 1979, A catalogue of stars with proper motions exceeding
annually (Minneapolis: Univ. Minnesota Press) [Google Scholar]
- Luyten, W. J. 1980, NLTT Catalogue, Volume III, 0 to − 30° (Minneapolis: Univ. Minnesota Press) [Google Scholar]
- Martin, B., & Wickramasinghe, D. T. 1981, MNRAS, 196, 23 [NASA ADS] [Google Scholar]
- Martin, N. J., Sandars, P. G. H., & Woodgate, G. K. 1968, Roy. Soc. London Proc. Ser. A, 305, 139 [NASA ADS] [CrossRef] [Google Scholar]
- Neale, L., & Tennyson, J. 1995, ApJ, 454, L169 [NASA ADS] [CrossRef] [Google Scholar]
- Nordhaus, J., Wellons, S., Spiegel, D. S., Metzger, B. D., & Blackman, E. G. 2011, Proc. National Academy of Science, 108, 3135 [Google Scholar]
- Paquette, C., Pelletier, C., Fontaine, G., & Michaud, G. 1986, ApJS, 61, 197 [NASA ADS] [CrossRef] [Google Scholar]
- Patat, F., Moehler, S., O’Brien, K., et al. 2011, A&A, 527, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Platais, I., Girard, T. M., Kozhurina-Platais, V., et al. 1998, AJ, 116, 2556 [NASA ADS] [CrossRef] [Google Scholar]
- Pokorny, R. S., Jones, H. R. A., Hambly, N. C., & Pinfield, D. J. 2004, A&A, 421, 763 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Potter, A. T., & Tout, C. A. 2010, MNRAS, 402, 1072 [NASA ADS] [CrossRef] [Google Scholar]
- Reid, I. N., Liebert, J., & Schmidt, G. D. 2001, ApJ, 550, L61 [NASA ADS] [CrossRef] [Google Scholar]
- Rohrmann, R. D., Althaus, L. G., & Kepler, S. O. 2011, MNRAS, 411, 781 [NASA ADS] [CrossRef] [Google Scholar]
- Salim, S., & Gould, A. 2002, ApJ, 575, L83 [NASA ADS] [CrossRef] [Google Scholar]
- Salim, S., & Gould, A. 2003, ApJ, 582, 1011 [NASA ADS] [CrossRef] [Google Scholar]
- Sion, E. M., Fritz, M. L., McMullin, J. P., & Lallo, M. D. 1988, AJ, 96, 251 [NASA ADS] [CrossRef] [Google Scholar]
- Sion, E. M., Kenyon, S. J., & Aannestad, P. A. 1990, ApJS, 72, 707 [NASA ADS] [CrossRef] [Google Scholar]
- Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
- Tommaseo, G., Pfeil, T., Revalde, G., et al. 2003, Eur. Phys. J. D, 25, 113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Unno, W. 1956, PASJ, 8, 108 [NASA ADS] [Google Scholar]
- Wickramasinghe, D. T., & Ferrario, L. 2005, MNRAS, 356, 1576 [NASA ADS] [CrossRef] [Google Scholar]
- Zuckerman, B., Koester, D., Reid, I. N., & Hünsch, M. 2003, ApJ, 596, 477 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
![]() |
Fig. 1 Flux spectra of the Ca II H&K lines in NLTT 10480 obtained with EFOSC2 (top) and FORS1 (middle), and circular polarization (bottom) spectrum obtained with FORS1. The dotted lines mark the position of the σ components and the dashed lines mark the π components (see Sect. 3.2.3). |
In the text |
![]() |
Fig. 2 Linear Zeeman splitting of Ca II (H and K) (top), Ca I 4226 Å (bottom left), and Hα (bottom right) lines at B = 0.519 MG compared to the observed Zeeman splitting in the spectrum of NLTT 10480. The dotted lines indicate the σ components and dashed lines indicate the π components. |
In the text |
![]() |
Fig. 3 Balmer lines (Hα to Hγ) and best-fitting model (Teff = 5410 K, log g = 8.0) showing prominent Zeeman-split Hα line (left panel), and (right panel) χ2 contours at 66, 90, and 99% confidence (dashed lines) in the (Teff, log g) plane. The measured calcium abundance from Ca II lines (log n(Ca)/n(H)) is drawn (full lines) as a function of the assumed temperature and surface gravity and labelled with the logarithm of the abundance. The thick line shows the range of temperature and surface gravity where the abundance measurements based on Ca I or Ca II are consistent. |
In the text |
![]() |
Fig. 4 X-shooter spectrum of the Ca I (bottom) and Ca II (H & K: middle and top, respectively) lines and best-fitting models at (log g = 8) Teff = 5400 K and log n(Ca)/n(H) = −10.3 (thick lines), and Teff = 4900 K and log n(Ca)/n(H) = −10.45 (thin lines). The wavelength scale is normalized to Δλ1 ≡ 4.67 × 10-7λ2B, where B = 0.513 MG. Therefore, the line positions match the factors glml − gumu listed in Table 4 (vertical dashed lines). All models assume ψ = 60°. |
In the text |
![]() |
Fig. 5 Spectra in the vicinity of the strongest Fe I (3734.8638, 3749.4851, and 3758.2330 Å), Si I (3905.523 Å), Al I (3944.0060 and 3961.5200 Å), and Na I (5889.951 and 5895.924 Å) lines compared to the predicted Zeeman-split line positions assuming a magnetic field of 0.519 MG. |
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.