Free Access
Volume 585, January 2016
Article Number A115
Number of page(s) 17
Section Stellar atmospheres
Published online 05 January 2016

© ESO, 2016

1. Introduction

HD 188112 is one of the brightest known (V = 10.2 mag) hot subdwarf B (sdB) stars, and it is a peculiar one. Subdwarf B stars are low-mass (~0.5 M) helium-core burning objects with a hydrogen envelope that is not massive enough to sustain significant hydrogen shell burning. These stars form the extreme horizontal branch (EHB), which is a hot extension (Teff 21 000 K) of the horizontal branch. The majority of sdB stars are thus found between the zero-age EHB and the terminal-age EHB, which corresponds to the central helium burning phase (Heber 2009). The subsequent, more rapid evolution during helium shell burning brings the star above the EHB.

Given its fundamental parameters (Teff = 21 500 K, log g = 5.66), HD 188112 was classified as an sdB; however, it lies at an odd position in the log gTeff diagram, i.e., significantly below the zero-age EHB (Heber et al. 2003). This is an indication that the star has too little mass to be He-core burning. The low-mass nature of HD 188112 was confirmed in Heber et al. (2003) by comparing the trigonometric distance of the star (from the Hipparcos parallax) with the spectroscopic distance; a mass of ~0.24 M is required for both values to match. Comparing the fundamental parameters of HD 188112 with models for post red giant branch evolution (Driebe et al. 1998) led to a similar mass of ~0.23 M (see Figure 4 in Heber et al. 2003). Such a star is likely to be formed via unstable mass transfer during its red giant phase, which leads to the formation of a common envelope. The subsequent spiralling-in of both components, and the eventual ejection of the envelope will leave behind a short-period binary system. If the envelope is removed before the core helium burning begins, the star evolves with an inert helium core and becomes a helium core white dwarf (WD; Han et al. 2002, 2003). Such low-mass evolved stars need to be formed via close binary evolution because the universe is not old enough to produce such a star through single-star evolution (Marsh et al. 1995). It is thus no surprise that HD 188112 was found to be radial velocity (RV) variable with an orbital period P = 0.606585 d, a semi-amplitude K1 = 188.3 km s-1, a systemic velocity γ0 = 26.6 km s-1, and a miminum companion mass of 0.73 M (Heber et al. 2003). This indicates that the companion is a rather massive WD, as a main-sequence (MS) star of such a mass would be detectable in the spectral energy distribution of HD 188112, which does not show any deviation from that of a single star.

At the time of the discovery of HD 188112, very few low-mass He-core objects were known. The first was discovered as the companion of a neutron star in a millisecond pulsar (van Kerkwijk et al. 1996); another serendipitous discovery followed a few years later, among the first data release of the Sloan Digital Sky Survey (SDSS; Liebert et al. 2004). No more than a handful of such objects had been spectroscopically confirmed (Eisenstein et al. 2006; Kawka et al. 2006; Kilic et al. 2007) when a dedicated search for what are now called extremely low-mass (ELM) WDs was undertaken in the SDSS data base (Gianninas et al. 2015; Brown et al. 2013, and references therein). Because these stars are the product of close binary evolution, they are monitored for radial velocity variations, and the majority of them have indeed been found to have an orbital period shorter than a day. By determining the mass of the star and the orbital parameters of the system, it is possible to derive a minimum mass for the companion, which allows to study the future outcome of these systems. Some of them will merge within a Hubble time, and some might be progenitors of type Ia supernovae (SNe Ia). Systems with very short periods (0.2 d) are expected to merge within a Hubble time. To date, 38 of these systems have been found (Gianninas et al. 2015; Brown et al. 2013, and references therein). Their fate depends on the mass of both components. If the mass of the companion is high enough, then it might be a SN Ia progenitor. Once a He-WD has filled its Roche lobe, a typical scenario involves a C/O core WD companion that starts accreting He-rich material from the He-WD. Once a sufficiently large shell of He has been accreted, a detonation can be triggered in the He-shell resulting in underluminous type Ia SNe (Bildsten et al. 2007). However, some simulations have shown that the explosion of the He-shell almost inevitably triggers a subsequent detonation in the C/O core. Depending on the mass of the C/O core, this might produce SNe Ia explosions (Livne & Glasner 1990; Fink et al. 2007, 2010; Pakmor et al. 2013). This double-detonation scenario happens at a sub-Chandrasekhar mass and a C/O WD accretor of ~1.0 M is a good candidate for producing a normal SN Ia (Sim et al. 2010). If the mass of the companion is much lower than 1 M then extreme-helium stars, R Coronae Borealis (RCrB) stars, single hot subdwarfs, or single WDs can ultimately be formed, depending on the mass ratio of the components. This is why, in order to constrain the future product of the binary systems hosting these ELM WDs, the masses of both components need to be determined accurately.

The mass of the ELM WD is usually determined by comparing the atmospheric parameters (effective temperature, Teff, and surface gravity, log g) of the star with theoretical evolutionary tracks for different masses (Driebe et al. 1998; Panei et al. 2007; Kilic et al. 2010; Althaus et al. 2013). However, the precision of such mass determination is limited and also depends on the evolutionary track used. In addition, the mass of the companion depends on the inclination of the system, which remains unknown unless eclipses are seen (Steinfadt et al. 2010; Brown et al. 2011; Vennes et al. 2011; Parsons et al. 2011; Kilic et al. 2014; Hallakoun et al. 2015) or distortion effects (such as ellipsoidal variations) in the light curve allow the system’s parameters to be constrained (Kilic et al. 2011b; Hermes et al. 2014a). However, these are very rare cases, and for the vast majority of systems only a minimum companion mass can be determined or, using the mean inclination angle for a random sample, a most likely companion mass can be estimated.

Given all of these new ELM systems found (more than 70), HD 188112 might look like a common system among others, but this is not the case. First, when compared to the majority of ELM systems, its combination of hot temperature and low surface gravity makes it stand out on the log gTeff plane (see, e.g., Fig. 1 in Brown et al. 2013). This indicates that HD 188112 is at a rather early stage of its cooling process, which is why it is also classified as a low-mass sdB. Such “young” objects are sometimes referred to as pre-ELM WDs. In terms of fundamental parameters, only four systems so far are roughly similar to HD 188112: KIC 6614501 (Silvotti et al. 2012), SDSS J1625+3632 (Kilic et al. 2011a), SDSS J0815+2309 (Brown et al. 2013), and GALEX J0805-1058 (Kawka et al. 2015). Second, it is by far the brightest of the known ELM systems; it has a parallax measured by Hipparcos, and we can expect a Gaia measurement soon. The parallax allows a mass to be derived for the star that is independent of evolutionary models. Additionally, its brightness allows high-resolution, high signal-to-noise spectroscopy.

As an alternative to finding eclipses or modeling light-curves of distorted stars, if one assumes the system to be tidally locked, then its inclination can be measured via the rotational broadening (vrot sin i) of metallic lines (Geier et al. 2010). The surface rotation velocity of the star (vrot) is estimated from the orbital period (=rotation period) and then, by measuring the rotational broadening, the inclination can be derived. Given the fact that the optical spectrum of the star is devoid of metallic lines, except for a magnesium feature too weak for such a measurement, we turned to the UV range instead where at least some magnesium lines are expected to be observed. The UV spectra of HD 188112 turned out to be much richer in metallic lines than anticipated, allowing an in-depth chemical composition analysis of the star. In our preliminary analysis we presented abundances for Mg, Si, and Fe (Latour et al. 2015); in this paper, we extend the measured abundances to a dozen additional elements. This is of great interest since detailed abundance analyses of ELM WD are still scarce (Hermes et al. 2014b; Gianninas et al. 2014b; Kaplan et al. 2013). We present the UV observations in Sect. 2 and review some of the binary parameters of the system in Sect. 3. Section 4 is dedicated to the spectroscopic analysis, it includes a description of the models, the determination of the rotational broadening, and the abundances of all metallic species seen in the spectra. Finally, a discussion follows in Sect. 5.

thumbnail Fig. 1

Radial velocity curve, with measurements from different sources are indicated with different colors (see Sect. 3). The “optical old” refers to values published in Heber et al. (2003) and Edelmann et al. (2005).

Open with DEXTER

2. Observations

Ultraviolet spectra of HD 188112 were obtained with the Space Telescope Imaging Spectrograph (STIS) installed on the Hubble Space Telescope (HST) during cycle 20 (program ID 12865) on April 27, 2013.

First there is a set of 22 exposures of 120 s each taken with the E140H grating (R = 114 000) covering the wavelength range 1242–1440 Å. Because of the radial velocity variations of the star and the high resolution of the observations, a series of short exposures was necessary in order to avoid orbital smearing of the metallic lines1. This exposure time leads to ~3 km s-1 differences between subsequent spectra. The data retrieved from the HST Archive were reduced and processed with the CALSTIS pipeline. Radial velocities were measured for each exposure and the spectra were shifted to rest wavelength (information on the radial velocity measurements is presented in the following section). The 22 “shifted” spectra were then coadded into a single spectrum suitable for spectroscopic analysis.

Additionally, two near-UV observations were taken with the E230H grating (also R = 114 000) in Time-Tag mode. The wavelength range covered by these spectra is 2666–2932 Å. There is a short (1040.2 s, hereafter called NUV1) and a long (2860.2 s, hereafter called NUV2) exposure, taken at different orbital phases. Again, to avoid orbital smearing of the spectral lines, the retrieved data were processed following the STIS Time-Tag Analysis Guide2 in order to create subexposures that were then reduced using CALSTIS. This way, we created 8 spectra of 130 s for the NUV1 observation, and 21 spectra of 136.2 s for the NUV2. Radial velocities were measured for each spectrum that was then shifted to rest wavelength. The main spectral features of the near-UV spectra are a series of four Mg ii lines around 2797 Å, among which the two resonance lines also have interstellar (IS) components. The NUV1 observation was taken far away from the conjunction phase, so the IS components are distinct from the stellar lines. This is not the case though for the NUV2 observation; the first spectra of the series show blending of the stellar and IS resonance lines. For this reason, when we coadded the subexposures of NUV2, we included only the 11 spectra where the lines are not blended. As for NUV1, the 8 subexposures were combined together.

3. Binary system parameters

The HST spectra were obtained for spectroscopic analysis, but the radial velocities measured for the individual exposures are also useful to put new constraints on the orbital parameters determined in Heber et al. (2003) and Edelmann et al. (2005), or at least to confirm them. By combining our new values with the previous measurements taken between 2000 and 2004, we extended the temporal coverage to 13 years.

Radial velocities were measured in the UV spectra using strong unblended lines: the Si iv resonance lines (1393.7 and 1402.8 Å), Si iii 1296.7 Å, and Fe iii 1273.8 Å. For the NUV spectra, the three strongest Mg ii lines were used, except for six subexposures of NUV2 where the resonance lines were strongly blended with the IS lines. In this case, only the 2798.8 Å line was used. The mean values of individual measurements were used, which gave maximum standard deviations of 2.5 km s-1, but typically of ~1.5 km s-1. In addition to the previously published optical RVs, we also added five values kindly provided by J. Norris, from low-resolution spectra taken at the Anglo-Australian Telescope (AAT) in 2003.

A periodogram analysis was performed, including the 98 RV measurements and our resulting folded RV curve is shown in Fig. 1. The orbital parameters we obtained are in very good agreement with those stated in Heber et al. (2003); we have P = 0.60658584 (±5e8) d, K1 = 188.7 ± 0.2 km s-1, and γ0 = 26.6 ± 0.2 km s-1. This leads to a mass function (1)of f(m) = 0.4221 M. By allowing for an eccentric orbit during the fitting procedure, we derive an upper limit of e 0.00006. This low eccentricity value is consistent with a circular orbit.

The parallax of HD 188112 from the new reduction of the Hipparcos catalogue is 13.6 ± 1.7 mas (van Leeuwen 2007). We then used the V (10.22) and B (10.01) magnitudes from the TYCHO 2 catalogue (Høg et al. 2000) combined with the synthetic flux from a model atmosphere appropriate for HD 188112 in order to derive the stellar radius. We did not consider any reddening due to the proximity of the star. Both magnitudes could be satisfactorily reproduced without extinction. This allowed us to derive a spectroscopic mass of 0.245 M and a radius of 0.12 R. The main uncertainty results from the parallax measurement. The spectroscopic mass indicates a Mcomp ≥ 0.73 M and a most likely mass Mcomp(52°) = 1.21 M (Gray 2008).

The latest development in terms of evolutionary sequences for ELM WDs is due to Althaus et al. (2013). Their models predict that for masses in the range 0.180.4 M, the star is expected to undergo H-shell flashes (in CNO burning) during the early cooling phase. This can happen multiple times while the star is “trying” to contract and cool down, each time following tracks in a slightly different region of the log gTeff diagram (see Figs. 3 and 4 of Althaus et al. 2013). This makes the determination of the mass a bit more complicated since it is not known if the star is on its final or an intermediate cooling sequence. Moreover, the different cooling tracks of stars from different masses are found in similar places on the HR diagram. The approach of Althaus et al. (2013) is a statistical one, weighting the time spent by models of different masses in each region of the log gTeff plane. This way, they determine a mass of 0.211 ± 0.0175 M for HD 188112  as well as a cooling age of 329 ± 207 Myr. This mass is a bit lower than that determined by the parallax, but it nevertheless fits well within the spectroscopic mass uncertainty. Assuming 0.21 M for HD 188112 leads to Mcomp ≥ 0.70 M and a most probable mass Mcomp(52°) = 1.17 M. As mentioned before, the absence of NIR excess in the observed flux of the star excludes the possibility of a main-sequence companion, thus implying that it has to be a rather massive WD.

thumbnail Fig. 2

Fit of the Mg ii lines in the NUV1 spectrum for vrot sin i = 7.9 km s-1 and vmicro = 1.1 km s-1. The IS Mg lines are shown in grey and were not included in the fit.

Open with DEXTER

4. Spectroscopic analysis

4.1. Model atmospheres and synthetic spectra

The model atomspheres used for our analysis are computed following the hybrid non-local thermodynamical equilibrium (NLTE) approach discussed in Przybilla et al. (2006a) and Nieva & Przybilla (2007, 2008); it is a combination of three codes: Atlas, Detail, and Surface (hereafter ADS). The structure of the atmosphere, such as the temperature stratification and density, is based on line-blanketed, plane-parallel, homogeneous, and hydrostatic LTE-model atmospheres computed with the Atlas12 code (Kurucz 1996). The NLTE effects are accounted for in the computation of the atomic population numbers via updated versions of the Detail and Surface codes (Giddings 1981). Using the Atlas12 LTE atmospheric structure as input, Detail solves the coupled radiative transfer and statistical equilibrium equations to obtain population numbers in non-LTE. This is then used by Surface, which produces the final synthetic spectrum using more detailed line-broadening data. The ADS approach is successfully used to study some types of hot stars (e.g., OB-type stars and BA-type supergiants; Nieva & Przybilla 2012). It has also been tested in the sdB regime (Przybilla et al. 2006a). Compared to the fully NLTE approach, the hybrid approach is much less time-consuming in terms of computation. Given the atmospheric parameters of HD 188112, no strong NLTE effects are expected on the temperature–density stratification of the atmosphere and the ADS models are appropriate. It must be kept in mind that an NLTE computation of atomic populations requires specific model atoms. Such model atoms are available only for a limited number of atomic species (see Sect. 4.3). When no such models exist, the populations are computed assuming the LTE approximation.

Given the very low metallicity expected for the star, a metallicity of 3% (by numbers) solar was used in the computation of the Atlas12 LTE model atmospheres. Afterward, small grids of different abundances were computed for all the atomic species seen in our spectra. It should also be mentioned that the He abundance was kept fixed at the value of log N(He)/N(H) = 5.0 derived in Heber et al. (2003). In order to generate spectra including the lines of all visible elements, normalized synthetic spectra computed for the various elements were multiplied following the technique described in Irrgang et al. (2014). This approach is appropriate when blending of lines originating from different species is sparse, as is the case for HD 188112. The metallic lines are then fitted using a standard χ2 minimization technique.

4.2. Rotational velocity measurement

As a first step, we measured the rotational velocity (vrot sin i) using the Mg ii lines in the NUV spectra, as well as the Fe lines in the FUV spectra. At the same time, we also fitted the microturbulent velocity as it can contribute to line broadening, even though it is not expected to be very important for hot subdwarf stars because of their stable and radiative atomspheres. Microturbulence has to be taken into account during the model atmosphere computation, as opposed to the rotational broadening that is applied afterward to the synthetic spectra. So we computed models with microturbulence between 0 and 8 km s-1 only for Mg and Fe. Microturbulence affects the strong and weak lines in different ways, so the choice of Mg and Fe allowed lines of various strength to be fit in the NUV and FUV spectra. The rotational broadening, microturbulence, magnesium, and iron abundances were thus fitted simultaneously. The microturbulence could be constrained to vmicro 2 km s-1; its effect below this value is rather small and did not influence the resulting χ2 much. We measured a rotational broadening of vrot sin i = 7.9 ± 0.3 km s-1. The resulting fit for the Mg lines in the NUV1 spectrum and the residual of the fit are shown in Fig. 2. The IS lines are plotted in grey and were excluded from the fitting procedure.

If it is assumed that the system is in synchronous rotation, then the star’s rotation period is the same as the orbital period. Using the mass of the star to determine the radius, the equatorial rotation velocity can be derived: (2)Combined with the measured v sin i, it yields the inclination angle. For a mass of 0.21 M, as predicted by theoretical models, we find i = 57°± 4° and a companion mass of 1.0 M.

If we use the slightly higher sdB mass resulting from the Hipparcos parallax (0.245 M), we find i = 52°± 3° and a companion mass of 1.2 M3.

4.3. Metal abundance determinations

Our FUV spectrum reveals a number of lines from various elements, namely Si, Al, P, S, Ti, Cr, Mn, Fe, Ni, Zn, Ga, Sn, and Pb. The NUV spectrum is rather poor in metallic lines: besides the strong Mg lines, only a handful of weak Fe ii lines are visible. Our model atmospheres were computed using the parameters derived by Heber et al. (2003): Teff = 21 500 ± 500 K, log g = 5.66 ± 0.06, and log N(He)/N(H) = 5.0. Our grid also included the extreme parameters resulting from the uncertainties on Teff and log g, which were used to check the effects of changing those parameters on the metallic abundances. Microturbulence was kept to zero in our metal grids and v sin i = 7.9 km s-1 was used during the fitting procedures.

As mentioned earlier, in order to compute NLTE population numbers, detailed model atoms are needed. To date, the model atoms available for Detail cover C, N, O, Mg, Al, Si, S, and Fe, so these elements were treated in NLTE in Detail, while all others were computed using the LTE approximation. We note that the LTE approximation is expected to be adequate given the atmospheric parameters of the star (Przybilla et al. 2006b) and it is the usual approach for abundance analysis of white dwarfs and hot subdwarfs having similar or higher Teff. Nevertheless, we also computed LTE models for Mg, Al, Si, S, and Fe in order to compare abundances derived with both types of models. This is discussed in more detail in Sect. 5.1. Our resulting abundances are listed in Table 1 where the results obtained with NLTE and LTE models are indicated in two separate columns. The comparison between our final synthetic spectrum and the entire UV spectrum can be found in Fig. A.1.

Table 1

Abundances of metals detected in the spectra of HD 188112.

4.3.1. Upper limits for C, N, and O

No lines of the C, N, O elements are visible in the UV spectra. This is not very surprising in the cases of nitrogen and oxygen as no strong lines are predicted in our UV ranges.

For oxygen, the strongest predicted features are three O i lines at 1302, 1304, and 1305 Å. They were used to put an upper limit of log N(O)/N(H) = 6.3. For nitrogen nothing useful is predicted in the UV range, so we turned toward the optical instead, where we used the high-resolution FEROS spectra of HD 188112 obtained by Edelmann et al. (2005). The N ii lines at 3995 and 4630 Å were used to put an upper limit of log N(N)/N(H) = 5.0, which is, however, only slightly below the solar value.

The situation is different for carbon. Three strong lines are predicted in the UV range, including C ii resonance lines at 1334.5 and 1335.7 Å. The absence of these lines in our spectrum results in an extremely low upper limit of log N(C)/N(H) = 9.6. This is six orders of magnitude below the solar value. Carbon has essentially vanished from the surface of the star.

4.3.2. Mg, Al, Si, P, and S

thumbnail Fig. 3

Best fits for the P iii line in the FUV spectrum as well as for the Ca ii k line in an optical FEROS spectrum.

Open with DEXTER

Only four Mg lines are visible in the NUV spectra. Their abundances were fitted during the determination of the projected rotational velocity (see Fig. 2). We also fitted them in the NUV2 spectra, shown in the last panel of Fig. A.1. The fits with both spectra lead to very similar results, with an abundance of 1/100 solar. This is essentially the same abundance as was first estimated by Heber et al. (2003) based on the Mg ii λ4481 line visible in the optical spectrum.

The FUV spectrum shows the Al iii doublet at 1379, 1384 Å that matches an abundance of log N(Al)/N(H) = 7.37.

Only two very weak sulphur lines are visible at the blue end of the FUV spectra, but they were nevertheless used to constrain the abundance of this element to log N(S)/N(H) = −8.0.

One phosphorus line was used to fit: P iii λ1344.32. Because of the lack of a line list for this ion in Surface, we used the public codes TLUSTY and SYNSPEC4 to perform an independent fit of the line. TLUSTY allows the computation of NLTE stellar atmospheres, assuming plane-parallel geometry and hydrostatic and thermal equilibrium; SYNSPEC is then used to compute synthetic spectra, given an atmospheric structure (Lanz & Hubeny 2007). An NLTE model atmosphere was computed including hydrogen and helium only and as no model atom is available for P iii, we included it in an LTE way through Synspec. The resulting fit can be seen in Fig. 3.

Silicon is an interesting element since it shows strong lines originating from three different ionization stages. Lines originating from all three of them are well reproduced with an abundance of log N(Si)/N(H) = 7.3. This is a good sign that the ionization equilibrium for this element is correct given the hybrid NLTE approach and the effective temperature of the star. Interstellar components to Si iii λλ1301, 1303 are present, but are not blended with the photospheric lines due to the relatively high radial velocity of the star during the FUV exposures. The IS features seen in our spectra are broadened by the radial velocity correction applied to the individual exposures, which explains their unusual shape in Fig. A.1.

4.3.3. Ca

Calcium is an important element because it has been claimed that all ELM WDs with log g 5.9 show Ca in their optical spectra (Hermes et al. 2014b; Gianninas et al. 2014a). This was first illustrated by Kaplan et al. (2013) who combined results from different studies of ELM WDs in their Fig. 6. However, they erroneously indicated that HD 188112 showed Ca ii in its optical spectra, confusing it with the Mg ii feature mentioned in Heber et al. (2003). Nevertheless, since the tendency is true for all other ELM WDs analyzed so far, we went back to the optical spectra of HD 188112 to look for this expected Ca ii K line. After correcting a few spectra for their radial velocity, and scrupulously comparing them, we indeed found weak calcium features: one stable line at the predicted wavelength, thus of photospheric origin, and one moving around, being a very small IS component. As can be seen in Fig. 3, this line is very weak indeed, but present in the different spectra, thus real. It is understandable why it was first overlooked in Heber et al. (2003). It must also be kept in mind that most of the ELM WDs in which Ca has been measured have Teff < 16 000 K. At high temperatures, the strength of these lines rapidly decreases as Ca ii becomes ionized. Nevertheless, the tiny line we found in HD 188112 indicates an abundance of log N(Ca)/N(H) = 7.3. Calcium was also measured in another ELM WD (SDSS J0815+2309) having essentially the same temperature as HD 188112 but a slightly higher surface gravity (Gianninas et al. 2014a). Compared to HD 188112, the Ca abundance derived for J0815+2309 is a hundred times higher (log N(Ca)/N(H) = 5.27). As for other ELM WDs, the Ca abundance of HD 188112 is similar to the lowest abundances measured in the sample of Gianninas et al. (2014a).

4.3.4. The iron-group elements

Titanium, chromium, iron, and nickel show many strong lines from their main ionization stage (iii) that could be reproduced well by our synthetic spectra. Two well-defined lines of Mn iii are seen at 1283 and 1287 Å, and two lines of Zn iii are also visible at 1318 and 1328 Å. The resulting abundances are listed in Table 1.

Iron and nickel are the only iron-group elements that also show lines from their singly ionized stage. There is a small discrepancy in the abundances measured with the Fe iii lines in the FUV and the Fe ii lines mainly found in the NUV. Compared with the abundance of log N(Fe iii)/N(H)5 = 5.87 ± 0.07, the Fe ii lines need a slightly higher abundance of 5.6 ± 0.1 to be reproduced. Another investigation of UV iron lines in B-type stars using ADS models found a scatter of ~0.25 dex for FUV lines when abundances measured from different lines are compared (Schaffenroth 2015). The reason for this, possibly related to the ADS models, will be explored further in a future publication (Schaffenroth et al., in prep.). In any case, this scatter is similar to the discrepancy we found between the FUV and NUV lines of HD 188112. We thus adopt a mean abundance of −5.75 ± 0.15 for iron, as indicated in Table 1.

We also encountered a discrepancy between lines originating from both ionization stages of nickel, but this time with a significant difference. In the case of nickel, fitting the Ni ii lines results in an abundance of log N(Ni ii)/N(H) = 7.25, which is 0.55 dex below what is predicted by the main ionization stage (log N(Ni ii)/N(H) = 6.7 when using the ADS models). In other words, at the abundance derived from the Ni iii lines, the Ni ii lines are much too strong, as can be seen in Fig. A.1. Something to keep in mind is that, unlike iron, nickel is treated using the LTE approximation in ADS. As nickel model atoms are available for TLUSTY (but only for Ni iii to Ni vi), we built a small grid of NLTE models with various Ni abundances in order to fit the Ni iii lines. We obtained an abundance of log N(Ni iii)/N(H) = −6.6 ± 0.1 dex, essentially the same as derived with ADS. This is not surprising since Ni iii is the main ionization stage. It is also reassuring in the sense that we can expect the LTE abundances derived for the other iron-group elements via their main ionization stage to be reliable. However, owing to the lack of a Ni ii model atom for TLUSTY, this ion was treated in LTE and an abundance discrepancy similar to the one seen with the ADS models was observed.

thumbnail Fig. 4

Comparison between parts of the FUV spectrum featuring mainly the Si lines and an LTE model spectrum, in blue. The LTE model was computed using the NLTE abundances stated in Table 1. The Si ii, Fe ii, and S ii lines appear to be much stronger than their observed counterparts, while the Si iii lines are nicely reproduced.

Open with DEXTER

thumbnail Fig. 5

Summary of the determined chemical composition of HD 188112. The top panel shows the abundances relative to the solar ones (Asplund et al. 2009) and the bottom panel the absolute values. Down-pointing triangles indicate upper limits determined for C, N, and O. The stars in the bottom panel indicate abundances predicted by radiative levitation (P. Chayer, priv. comm.). No support is expected for oxygen and magnesium, thus the stars are at the bottom of the plot. The green crosses indicate the average abundances for sdB stars.

Open with DEXTER

4.3.5. Trans-iron elements: Ga, Sn, and Pb

We found lines originating from three trans-iron elements, gallium, tin, and lead, in the spectra of HD 188112. Such heavy elements are also seen in hot subdwarfs (O’Toole 2004; O’Toole & Heber 2006) as well as hot DA and DO WDs (Vennes et al. 2005; Werner et al. 2012; Rauch et al. 2013). However, they have not been detected yet in colder WDs, with the exception of the heavily polluted WD GD 362, in which a strontium line has been detected (Zuckerman et al. 2007). The metallic elements in this cool star are accreted from a debris disk surrounding the star. With its temperature of 21 500 K, HD 188112 is among the coolest WD/sdB stars for which trans-iron elements have been detected.

The Ga abundance of log N(Ga)/N(H) = 9.55 was determined using the λ1414.4 resonance line of Ga ii as well as two weak Ga iii lines at 1267 Å and 1293 Å. All of them are reproduced well with this abundance.

A sole line of lead is present in our spectra, the 1313 Å resonance line of Pb iv, which indicates an abundance of log N(Pb)/N(H) = 10.0, slightly above the solar value.

We see the resonance line of Sn iii at 1251.39 Å, as well as the resonance doublet of Sn iv λλ 1314, 1437. The fit of the Sn iii line gives an abundance of N(Sn)/N(H) = 10.6. However, at this abundance, the Sn iv lines are barely visible in our synthetic spectra. In order to reproduce them, we need an abundance of 8.4, which is more than 2 dex higher than what is indicated by the Sn iii line. This is reminiscent of the problem encountered with nickel, but much more dramatic. Given the fact that these three transitions are resonances, their oscillator strengths, taken from Morton (2000) are considered to be reliable. The Sn iii line is blended with a Cr iii line (1251.417 Å), which is however predicted to be rather weak at the Cr abundance determined. Indeed, a fit of the Sn iii line without the contribution of the Cr iii feature leads to an abundance only slightly higher, 10.5 dex, still far from the 8.4 dex needed to reproduce the Sn iv lines. As we did for nickel, we turned toward NLTE models in an attempt to solve the tin abundance discrepancy. An NLTE fit of the Sn lines was performed using this time the Tübingen NLTE model-atmosphere package (TMAP; Rauch & Deetjen 2003) because a simple NLTE Sn model atom for this code was built in the course of the spectral analysis of the hot DA white dwarf G191-B2B (Rauch et al. 2013; T. Rauch, private comm.). With NLTE models, the Sn iii line leads to an abundance of 10.0, and the Sn iv lines indicate 8.2 (1314.5 Å) and 8.5 (1437.5 Å). According to these results, NLTE effects apparently cannot solve the “tin problem”. However, because the model atom for tin is rather simple (3 and 6 NLTE levels for Sn iii and iv, respectively) and resonance lines can be very sensitive to NLTE effects, it would be premature to exclude these effects as the cause of the discrepancy and to blame it on the atomic data. Because Sn iii has the same ionization energy as Fe iii and Cr iii, it is expected to be the main ionization stage, as for the iron-group elements, so we adopt an abundance of log N(Sn)/N(H) = 10.6. This is based on our ADS results, but the upper limit includes the value indicated by TMAP.

5. Discussion

5.1. Non-LTE effects

The discrepancies found between the ionization stages of Ni and Sn show the same trend: the lowest ionization stage visible indicates a lower abundance. This is sometimes seen in much hotter stars and is a signature of NLTE effects when analyses are made using the LTE approximation. As stated in Dreizler & Werner (1993), in a consistent NLTE model, ionization equilibrium is shifted to higher stages. In our case, this could explain why the Ni ii lines are too strong in our models, while the Sn iv lines are too weak. However, these effects are not really expected to appear in an sdB at 21 500 K, although some line-by-line differences between LTE and NLTE treatment were found in HD 205805, an sdB star at 25 000 K and log g = 5.0 (Przybilla et al. 2006b). This is why we decided to compute additional LTE models for Mg, Al, Si, S, and Fe and to redo the fitting procedure in order to derive LTE abundances for these elements6. The resulting abundances are indicated in Table 1. For magnesium and aluminium, the LTE and NLTE abundances are similar, while there is a difference of 0.4 dex for the sulfur abundance. Most interesting is silicon; a simultaneous fit of all the Si lines leads to an abundance of log N(Si)/N(H) = 7.56 dex, but the quality of the fit is poor and it shows conspicuous hints of an incorrect ionization equilibrium. Instead, the lines originating from the different ions had to be fitted separately in order to be correctly reproduced, thus leading to discrepant abundances: log N(Si ii)/N(H) = 7.8 ± 0.1; log N(Si iii)/N(H) = 7.33 ± 0.05; and log N(Si iv)/N(H) = 7.1 ± 0.1.

The abundance of the main ionization stage (Si iii) remains the same as derived in NLTE, but the Si ii abundance is 0.5 dex lower. The Fe ii lines in the NUV spectrum also lead to an LTE abundance 0.5 dex lower than the value derived with NLTE models. To illustrate this effect on the spectral lines, we select four spectral ranges featuring Si lines and show the corresponding LTE spectrum in Fig. 4. The spectrum includes the abundances determined with the NLTE models. The Si ii lines (as well as few Fe ii lines and one S ii line) appear too strong, as can be seen with the Ni ii lines (Fig A.1).

It thus becomes clear that NLTE effects are present in HD 188112 and influence the line strengths, especially those of the non-dominant ionic species. For Si ii, S ii, and Fe ii the differences in strength of the modeled lines between LTE and NLTE lead to a ~0.5 dex shift in abundances. However, the LTE abundances derived for the main ionization stage (in this case doubly ionized) are consistent with the NLTE values. This result confirms that the abundances derived for the elements that could not be treated in NLTE can be regarded as reliable since the main ionization stage was used. The only exception is lead, for which only one Pb iv line was visible. In this case the derived abundance might be overestimated. The strongest gallium line is from Ga ii, but the LTE models simultaneously reproduced this line and the two weak Ga iii features. However, some mystery still remains about the Sn lines: it is not clear what is causing the 2 dex discrepancy between Sn iii and Sn iv. This “tin issue” is certainly worth further investigation in the future. Maybe more sophisticated NTLE model atoms (including more energy levels) would solve the problem, or perhaps the collisional radiative rates used for tin are not very reliable. There is also the possibility that the atomic data are inaccurate.

5.2. Metal abundances

Figure 5 summarizes the chemical composition of HD 188112 with upper limits for the C, N, O elements, and abundance values for 15 metallic elements. All of the elements, except lead, are found to have subsolar abundances, with the heavy elements (Ti and heavier) being the most abundant relative to solar (see top panel of Fig. 5). Lead is, however, slightly above the solar value. This abundance pattern suggests that diffusion mechanisms counteract gravitational settling in order to keep a relatively considerable amount of heavy elements in the photosphere.

5.2.1. Comparison with diffusion models

The bottom panel of Fig. 5 also features abundances predicted by radiative levitation (red stars) for the atomic species included in TOPbase7. These abundances were computed for a model with Teff = 21 000 K and log g = 5.6 according to the method described in Chayer et al. (1995) (P. Chayer, priv. comm.). In most cases (C, Al, Si, S, Fe) we measured abundances below the value predicted by radiative support, though the Fe abundance is quite close to the prediction. The two exceptions are Mg and Ca. Calcium is 1.5 dex higher than the predicted abundance while no support at all is expected for magnesium.

5.2.2. Comparison with ELM WDs

A UV spectrum has been analyzed for only one other ELM WD so far, namely that of GALEX J1717+6757 (Hermes et al. 2014b). This star has log g = 5.67, essentially identical to that of HD 188112, and Teff = 14 900 K, which gives a mass of ~0.19 M. Hermes et al. (2014b) detected lines from nine metals, and for all of them they derived an abundance higher (by rougly 0.5 to 2.0 dex) than that measured for HD 188112. In their discussion, Hermes et al. (2014b) point out that radiative levitation for N, O, Si, P, and Fe gives predictions that are in line with the abundance pattern seen in the star. One noteworthy exception is carbon for which strong radiative support was predicted, but the element is extremely underabundant (by a factor of 1000 with respect to the sun) in the atmosphere of J1717. Interestingly, this is reminiscent of the situation for HD 188112, where important radiative support is also predicted for carbon at 21 000 K, but the element is strongly depleted in the stellar atmosphere. Hermes et al. (2014b) argue that a possible explanation for carbon depletion might be a history of repeated CNO flashes, as predicted by evolutionary models with M 0.18 M. This might also apply to HD 188112, although we do not see the metal enrichment that might be associated with CNO-flash surface convection. Calcium is another element that defies radiative levitation predictions in J1717: more Ca is observed than predicted, just as it is for HD 188112. According to Fig. 5 in Hermes et al. (2014b) the abundance pattern of J1717 is not very different from those of WDs with ongoing accretion from debris discs, except for Mg and O. However, the authors conclude that debris-disk accretion is unlikely in the case of J1717 because of the discrepancy in oxygen abundance, and because accretion from a circumbinary disk is dynamically difficult to explain.

Two ELM WDs are known to be metal rich: PSR J1816+4510 (Teff = 16 000 K, log g = 4.9), the companion of a massive millisecond pulsar (Kaplan et al. 2013), and SDSS J0745+1949 (Teff = 8,380 K, log g = 6.2), a tidally distorted star (Gianninas et al. 2014b). The analyses of their optical spectra allowed for the abundance determination of a few elements. The abundances of Mg, Ca, Ti, Cr, and Fe turned out to be close to solar in the case of J0745 while He, Mg, Ca, Si, and Fe are about 10 times solar in J1816. The origin of the metals in these two stars is not fully understood. A recent shell flash is a possible explanation for the high metallicity of J1816, but the rather low mass inferred for J0745 (M = 0.16 M) does not suggest the star could have gone through such shell flash episodes. On the other hand, accretion of metals from a debris disk was discussed as a possibility for J0745, while the pulsar companion of J1816 prevents such an accretion.

Even though detailed abundances are known for only three other ELM-WDs, Gianninas et al. (2014a) measured abundances of He, Ca, and Mg in their sample of 61 stars (though such elements were visible only in a minority of stars). As mentioned previously, they noticed the presence of Ca lines in basically all of their stars having log g 6.0. HD 188112, which has low gravity, now also falls into this category thanks to the tiny line we identified in its optical spectra. The abundance we derived (log N(Ca)/N(H) = 7.3), though below the Gianninas et al. (2014a) detection threshold, is comparable to the lowest abundances measured in their cooler ELM WDs (we note here that the Ca ii lines get stronger at lower temperature). The same observation holds for Mg, where the abundance of HD 188112 is similar to the lowest values detected in the Gianninas et al. (2014a) sample.

5.2.3. Comparison with sdB stars

Given its fundamental parameters, HD 188112 is considered an sdB star, so it might be worth comparing its abundances with those of sdBs. To this end we also included in Fig. 5 the “upper averages” (including determinations and upper limits) of the Geier (2013) sample, consisting of a hundred sdBs. It thus appears that HD 188112 is a particularly metal-poor sdB, but also an especially cool one. PG1627+017 is a cool pulsating sdB (Teff~  2124 kK) that can be compared to HD 188112 in terms of effective temperature, but the lower log g of PG1627+017 places it onto the EHB, meaning that the star is likely to be more massive than HD 188112. Indeed, a mass of ~0.45 M was suggested by its seismic properties (Randall et al. 2006). The UV spectrum of PG1627+017 has allowed a detailed abundance determination, and the star is clearly more metal-rich than HD 188112 (see Fig. 6 in Blanchette et al. 2008). Carbon is the most extremely depleted element in HD 188112. Past studies of sdBs have shown a large star-to-star variation in atmospheric carbon. Very low carbon upper limits (log N(C)/N(C)) = 5.7 and 5.3) were determined for CD 24°731 and PG 1219+534, respectively (O’Toole & Heber 2006). The sdOB star Feige 110 is another example of a carbon poor star (Rauch et al. 2014; Heber et al. 1984). Lead is the most enriched element in HD 188112, and it has been identified in 22 out of the 27 hot subdwarf UV spectra examined by O’Toole (2004). Further abundance analysis of sdBs showed that lead enrichment varies from 10 to 1000 times solar (O’Toole & Heber 2006; Blanchette et al. 2008). These two studies also showed that trans-iron elements tend to be overabundant in the atmosphere of sdB stars. In this respect, the abundance pattern of HD 188112 is similar to other sdBs.

5.3. Binary system

In addition to the spectroscopic analysis, we attempted a better characterization of the HD 188112 binary system. By combining our new UV radial velocity measurement with the already published values, we refined the parameters of the binary system: period (0.60658584(5) d), amplitude (188.7 ± 0.2 km s-1), and systemic velocity (26.6 ± 0.2 km s-1). Access to a Hipparcos parallax for a star provides an independent mass measurement, i.e., independent of stellar evolution models. In our case, the spectroscopic mass of 0.245 M derived with the parallax distance agrees quite well with that indicated by evolutionary models of low-mass He-core WDs, which predict 0.211 ± 0.0175 M (Althaus et al. 2013). However, the uncertainty on the spectroscopic mass is rather large (due to the uncertainty on the parallax itself) and allows for a mass up to 0.32 M. Such a mass is at the lower limit of He-core burning for an sdB star (Han et al. 2002; Hu et al. 2007). Low-mass sdB stars are created when a massive progenitor (22.5 M) ignites helium under non-degenerate conditions. According to binary synthesis models, the second common envelope ejection channel can produce sdBs with masses down to 0.33 M in close binaries with a WD companion (Han et al. 2003). The position of HD 188112 in the log gTeff diagram would fit within the 0.33 M EHB band (see Fig. 7 from Silvotti et al. 2012). Although there is a possibility that HD 188112 is a He-core burning object, one must keep in mind that these low-mass sdBs occur much less frequently than the canonical (~0.47 M) sdBs and that the extreme upper limit of our spectroscopic mass coincides with the extreme lower limit of the He-core burning. Given the current parallax value, we thus consider it unlikely that HD 188112 is such an object. For the lower mass limit on HD 188112, both the parallax and evolutionary models give a value of ~0.19 M. In this case, the star is certainly a pre-ELM WD.

In our attempt to constrain the mass of the WD companion, we considered both the 0.21 and 0.245 M solutions indicated respectively by evolutionary models and the parallax. We derived a minimum mass of 0.7 M for the companion8. It is interesting to compare this minimum mass with the minimum mass distributions of sdB binaries with WD companions and of ELM WD binaries illustrated in Fig. 14 in Kupfer et al. (2015). When considering the distribution for sdB binaries, it appears that a 0.7 M companion is quite massive; only two known sdBs have a more massive WD companion9. After a survey was dedicated to finding massive companions to sdBs (Geier et al. 2011), it now seems that such systems are quite rare. On the other hand, when looking at the observed minimum mass distribution of ELM systems (data from Gianninas et al. 2014a), the companion masses are evenly distributed between 0.1 and 0.9 M with a few more massive outliers. In this context, the companion of HD 188112 would be on the massive side, but not outstandingly so.

We measured the rotational broadening of metallic lines in our UV spectra and derived vrot sin i = 7.9 ± 0.3 km s-1. In order to constrain the inclination of the system, we assumed it to be tidally locked, which led to inclination values between 49° and 61°. This range nicely overlaps with the most likely inclinations for a random sample. Our inclination range indicates companion masses between 0.92 and 1.33 M. The synchronization assumption thus seems reasonable, or at least does not lead to improbable results. Recent statistical investigations on the mass distribution of companions to ELM WDs found a mean mass of ~0.75 M, with σ~0.25 M (Andrews et al. 2014; Boffin 2015). Both studies found very similar results when assuming a Gaussian distribution to model the companion’s mass population and comparing it with the observed ELM WD systems of Brown et al. (2013) and Gianninas et al. (2014a). If the HD 188112 system is indeed synchronized, the companion would thus be on the massive side of these mass distributions. The possibility of a neutron star companion should not be discarded, as it would require an inclination angle only slightly below 49°, and would imply the star is rotating faster than synchronized. This last assumption is plausible considering the fact that evolutionary models predict that the rotation of the star should increase after the ejection of the common envelope, as the star is contracting toward the white dwarf regime. Given the orbital parameters of the HD 188112 system, the rotation period of the star can hardly be less than half the orbital period, in which case the inclination would be less than 25° and the companion mass above 6 M. A search for possible X-ray emission from the companion of HD 188112 might help to constrain its nature (Kilic et al. 2014).

It might be interesting to look at the possible future of the system. As the star contracts toward the cooling track it is likely to experience hydrogen shell flashes during which the star will expand and fill its Roche lobe (Istrate et al. 2014). Brief episodes of mass transfer might lead to mass loss of a few times 10-4 M during a shell flash (see Fig. 12 in Nelson et al. 2004). Because the hydrogen fuel is diminished both by burning and mass loss, the evolution of the star will speed up. In the long run, gravitational wave radiation will lead to a decay of the orbit resulting finally in a merger of the components after a few Hubble times. By the time of merging, the system has evolved into a degenerate He-core white dwarf and a massive companion, most likely a white dwarf of 0.7 to 1.3 M.

The outcome of the merger depends strongly on the mass of the white dwarf companion and its composition. If the mass is as low as 0.7 M, the white dwarf is of C/O composition and the merger most likely produces a RCrB star (Webbink 1984; Clayton et al. 2007). Mergers of He-core white dwarfs with more massive C/O white dwarfs may lead to SNe Ia (Pakmor et al. 2013; Dan et al. 2015). Some simulations showed that the He-shell ignition induced by the accretion of material onto a sub-Chandrasekhar C/O WD is very likely to trigger a detonation in the core, leading to a SN Ia explosion (double-detonation scenario, Fink et al. 2007, 2010; Moll & Woosley 2013). If the core fails to ignite10, then the helium explosion alone can produce a fainter type Ia SN (Bildsten et al. 2007). At the high end of the companion mass range (>1.1 M), the system would have a low mass ratio (q ≲ 0.20) and thus evolve via stable mass transfer into an AM Canum Venaticorum (AM CVn; Marsh et al. 2004; Dan et al. 2011), which is a potential SN Ia progenitor (Shen & Bildsten 2014). However, it is usually assumed that white dwarfs more massive than about 1.1 M have O/Ne/Mg cores and these objects, when accreting matter, are not likely to explode, but would rather collapse into a neutron star via electron capture (Nomoto & Kondo 1991; Shen & Bildsten 2014). On the other hand, recent investigations showed that near-Chandrasekhar mass (~1.11.3 M) white dwarfs may have a hybrid composition of C/O/Ne and could be prone to explode, leading to faint SNe of the observed type Iax (Chen et al. 2014; Kromer et al. 2015). In either case, the outcome of the system depends on the stability of the mass transfer, which in turn depends on the mass ratio of the components (Marsh et al. 2004). For mass ratio q in the range ~0.25–0.70 the mass transfer can be either stable or unstable (see Fig. 1 in Dan et al. 2011). For the HD 188112 system, q can be up to ~0.32 if the companion mass is close to the minimum value.

6. Conclusion

We began our analysis of the UV spectra of HD 188112 with the main goal of determining the projected rotational velocity, in order to put tighter constraints on the mass of the companion. However, the high-resolution spectra showed a large number of metallic lines that allowed us to perform a comprehensive abundance analysis of this pre-ELM WD. This is extremely interesting since the chemical compositions of only a few ELM WDs have been determined to date. Additionally, because it is a pre-ELM WD, HD 188112 represents an earlier evolutionary stage and comparing its metallic content with those of cooler, more compact, and thus older (in terms of cooling age) ELM WDs might help to determine the origin of their metals.

The spectral analysis was performed using hybrid NLTE model atmospheres, in which the NLTE effects are taken into account in the computation of the population numbers while the atmospheric structure is computed in LTE. The computation of NLTE populations was only possible for a limited number of species (C, N, O, Mg, Al, Si, S, Ca, and Fe); otherwise, the LTE approximation had to be used. By comparing both LTE and NLTE spectra for a few metallic species, we concluded that the spectral lines originating from the main ionization stages of an element were not affected, or were only marginally affected, by NLTE effects. Thus, they can be used reliably to derive abundances with LTE models. However, we clearly saw in our model spectra that NLTE effects affect lines of non-dominant ionic species. Fits of these lines led to discrepant abundances, by up to ~0.5 dex. As these effects are seen in HD 188112, they are likely to also be present in stars having a temperature higher than ~21 000 K. This is why we emphasize that it is more accurate to consider only lines originating from the dominant ionic species when using LTE models in order to derive appropriate abundances.

Based on our abundance determination and comparisons with other ELM WDs and hot subdwarf stars, HD 188112 appears to be an especially metal-poor object. Its abundance pattern suggests diffusion effects where gravitational settling is probably counteracted by radiative levitation. In constrast to other more metal-rich ELM WDs, HD 188112 does not require any “metal enrichment” mechanism to explain the presence of metals in its atmosphere. It is also unlikely that its metallicity will increase as it ultimately cools down and becomes more compact; radiative levitation should become weaker as Teff decreases and gravitational settling stronger as log g increases. However, given its mass, HD 188112 is likely to experience an H-shell flash that would temporarily increase the surface abundance of metals, especially He and N, although gravitational settling is expected to act quickly as the star contracts again. Considering the four detailed abundance analyses currently available for ELM/pre-ELM WDs, it seems that they show a rather wide variety of metallicities, from metal-rich to metal-poor, as is also observed for sdB stars.

Despite our efforts to constrain the mass of the HD 188112 companion, we are limited by the unknown inclination of the system and by the uncertainty on the mass of HD 188112 itself. By assuming tidally locked rotation, we could constrain the inclination to be between 49° and 61°. Unfortunately, we cannot currently test or verify whether the system is indeed tidally locked. Nevertheless, this assumption leads to inclinations close to the most likely ones for a random sample, so it does not appear to be inappropriate. This inclination range indicates a companion mass between 0.92 and 1.33 M, i.e., a rather massive WD, while the minimum mass is 0.7 M. For the mass of HD 188112, a value up to 0.32 M is possible according to the Hipparcos distance, in which case the star would be at the limit of the He-burning mass range. Gaia should provide us with a more precise distance.


Time-tag mode, as used for the NUV spectrum, would have been preferable, but the star is too bright in the far-UV.


Only the uncertainty on i is considered in these estimates of the companion mass.


We use the nomenclature Fe iii here and below to indicate the abundance derived by using lines from a specific ion.


We remind the reader that the hybrid approach of ADS only includes NLTE effects in the computation of the population numbers; in all cases (NLTE and LTE) the atmospheric structure is computed assuming the LTE approximation.


If we consider 0.19 M for HD 188112  the minimum mass of the companion decreases only slightly to 0.68 M.


HD 188112 was not included in the Kupfer et al. (2015) sample because of its low mass.


In the case of a very low-mass (and low-density) C/O WD companion or a high-mass WD companion with O/Ne composition, a secondary detonation may be less likely (Shen & Bildsten 2014).


This work was supported by the Deutsches Zentrum für Luft- und Raumfahrt (grant 50 OR 1315). Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with program no. 12865, cycle 20. We are most grateful to P. Chayer for sharing the results of his radiative levitation computation and T. Rauch for fitting the tin lines with TMAP model atmospheres. We also thank A. Istrate and H. Dreschel for interesting discussions. The work of F.K.R. is supported by the ARCHES prize of the German Ministry of Education and Research (BMBF) and by the DAAD/Go8 German-Australian exchange programme for travel support. WH and ST acknowledge support by project TRR 33 ’The Dark Universe’ of the German Research Foundation (DFG) This research has made use of ISIS functions provided by ECAP/Remeis observatory and MIT ( We thank John E. Davis for the development of the slxfig module used to prepare some of the figures in this paper.


  1. Althaus, L. G., Miller Bertolami, M. M., & Córsico, A. H. 2013, A&A, 557, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Andrews, J. J., Price-Whelan, A. M., & Agüeros, M. A. 2014, ApJ, 797, L32 [NASA ADS] [CrossRef] [Google Scholar]
  3. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bildsten, L., Shen, K. J., Weinberg, N. N., & Nelemans, G. 2007, ApJ, 662, L95 [NASA ADS] [CrossRef] [Google Scholar]
  5. Blanchette, J.-P., Chayer, P., Wesemael, F., et al. 2008, ApJ, 678, 1329 [NASA ADS] [CrossRef] [Google Scholar]
  6. Boffin, H. M. J. 2015, A&A, 575, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Brown, W. R., Kilic, M., Hermes, J. J., et al. 2011, ApJ, 737, L23 [NASA ADS] [CrossRef] [Google Scholar]
  8. Brown, W. R., Kilic, M., Allen de Prieto, C., Gianninas, A., & Kenyon, S. J. 2013, ApJ, 769, 66 [NASA ADS] [CrossRef] [Google Scholar]
  9. Chayer, P., Fontaine, G., & Wesemael, F. 1995, ApJS, 99, 189 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chen, M. C., Herwig, F., Denissenkov, P. A., & Paxton, B. 2014, MNRAS, 440, 1274 [NASA ADS] [CrossRef] [Google Scholar]
  11. Clayton, G. C., Geballe, T. R., Herwig, F., Fryer, C., & Asplund, M. 2007, ApJ, 662, 1220 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dan, M., Rosswog, S., Guillochon, J., & Ramirez-Ruiz, E. 2011, ApJ, 737, 89 [NASA ADS] [CrossRef] [Google Scholar]
  13. Dan, M., Guillochon, J., Brüggen, M., Ramirez-Ruiz, E., & Rosswog, S. 2015, MNRAS, 454, 4411 [NASA ADS] [CrossRef] [Google Scholar]
  14. Dreizler, S., & Werner, K. 1993, A&A, 278, 199 [Google Scholar]
  15. Driebe, T., Schoenberner, D., Bloecker, T., & Herwig, F. 1998, A&A, 339, 123 [Google Scholar]
  16. Edelmann, H., Heber, U., Altmann, M., Karl, C., & Lisker, T. 2005, A&A, 442, 1023 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Eisenstein, D. J., Liebert, J., Harris, H. C., et al. 2006, ApJS, 167, 40 [NASA ADS] [CrossRef] [Google Scholar]
  18. Fink, M., Hillebrandt, W., & Röpke, F. K. 2007, A&A, 476, 1133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Fink, M., Röpke, F. K., Hillebrandt, W., et al. 2010, A&A, 514, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Geier, S. 2013, A&A, 549, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Geier, S., Heber, U., Podsiadlowski, P., et al. 2010, A&A, 519, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Geier, S., Hirsch, H., Tillich, A., et al. 2011, A&A, 530, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Gianninas, A., Dufour, P., Kilic, M., et al. 2014a, ApJ, 794, 35 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gianninas, A., Hermes, J. J., Brown, W. R., et al. 2014b, ApJ, 781, 104 [NASA ADS] [CrossRef] [Google Scholar]
  25. Gianninas, A., Kilic, M., Brown, W. R., Canton, P., & Kenyon, S. J. 2015, ApJ, 812, 167 [NASA ADS] [CrossRef] [Google Scholar]
  26. Giddings, J. R. 1981, Ph.D. Thesis, University of London [Google Scholar]
  27. Gray, F. D. 2008, The Observation and Analysis of Stellar Photospheres (Cambridge University Press) [Google Scholar]
  28. Hallakoun, N., Maoz, D., Kilic, M., et al. 2015, ArXiv e-prints [arXiv:1507.06311] [Google Scholar]
  29. Han, Z., Podsiadlowski, P., Maxted, P. F. L., Marsh, T. R., & Ivanova, N. 2002, MNRAS, 336, 449 [NASA ADS] [CrossRef] [Google Scholar]
  30. Han, Z., Podsiadlowski, P., Maxted, P. F. L., & Marsh, T. R. 2003, MNRAS, 341, 669 [NASA ADS] [CrossRef] [Google Scholar]
  31. Heber, U. 2009, ARA&A, 47, 211 [NASA ADS] [CrossRef] [Google Scholar]
  32. Heber, U., Hamann, W.-R., Hunger, K., et al. 1984, A&A, 136, 331 [NASA ADS] [Google Scholar]
  33. Heber, U., Edelmann, H., Lisker, T., & Napiwotzki, R. 2003, A&A, 411, L477 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Hermes, J. J., Brown, W. R., Kilic, M., et al. 2014a, ApJ, 792, 39 [NASA ADS] [CrossRef] [Google Scholar]
  35. Hermes, J. J., Gänsicke, B. T., Koester, D., et al. 2014b, MNRAS, 444, 1674 [NASA ADS] [CrossRef] [Google Scholar]
  36. Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [NASA ADS] [Google Scholar]
  37. Hu, H., Nelemans, G., Østensen, R., et al. 2007, A&A, 473, 569 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Irrgang, A., Przybilla, N., Heber, U., et al. 2014, A&A, 565, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Istrate, A. G., Tauris, T. M., Langer, N., & Antoniadis, J. 2014, A&A, 571, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Kaplan, D. L., Bhalerao, V. B., van Kerkwijk, M. H., et al. 2013, ApJ, 765, 158 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kawka, A., Vennes, S., Oswalt, T. D., Smith, J. A., & Silvestri, N. M. 2006, ApJ, 643, L123 [NASA ADS] [CrossRef] [Google Scholar]
  42. Kawka, A., Vennes, S., O’Toole, S., et al. 2015, MNRAS, 450, 3514 [NASA ADS] [CrossRef] [Google Scholar]
  43. Kilic, M., Allen de Prieto, C., Brown, W. R., & Koester, D. 2007, ApJ, 660, 1451 [NASA ADS] [CrossRef] [Google Scholar]
  44. Kilic, M., Brown, W. R., Allen de Prieto, C., Kenyon, S. J., & Panei, J. A. 2010, ApJ, 716, 122 [NASA ADS] [CrossRef] [Google Scholar]
  45. Kilic, M., Brown, W. R., Allende Prieto, C., et al. 2011a, ApJ, 727, 3 [NASA ADS] [CrossRef] [Google Scholar]
  46. Kilic, M., Brown, W. R., Kenyon, S. J., et al. 2011b, MNRAS, 413, L101 [NASA ADS] [CrossRef] [Google Scholar]
  47. Kilic, M., Hermes, J. J., Gianninas, A., et al. 2014, MNRAS, 438, L26 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kromer, M., Ohlmann, S. T., Pakmor, R., et al. 2015, MNRAS, 450, 3045 [NASA ADS] [CrossRef] [Google Scholar]
  49. Kupfer, T., Geier, S., Heber, U., et al. 2015, A&A, 576, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Kurucz, R. L. 1996, in M.A.S.S., Model Atmospheres and Spectrum Synthesis, eds. S. J. Adelman, F. Kupka, & W. W. Weiss, ASP Conf. Ser., 108, 160 [Google Scholar]
  51. Lanz, T., & Hubeny, I. 2007, ApJS, 169, 83 [NASA ADS] [CrossRef] [Google Scholar]
  52. Latour, M., Irrgang, A., Heber, U., & Schaffenroth, V. 2015, in 19th European Workshop on White Dwarfs, eds. P. Dufour, P. Bergeron, & G. Fontaine, ASP Conf. Ser., 493, 9 [Google Scholar]
  53. Liebert, J., Bergeron, P., Eisenstein, D., et al. 2004, ApJ, 606, L147 [NASA ADS] [CrossRef] [Google Scholar]
  54. Livne, E., & Glasner, A. S. 1990, ApJ, 361, 244 [NASA ADS] [CrossRef] [Google Scholar]
  55. Marsh, T. R., Dhillon, V. S., & Duck, S. R. 1995, MNRAS, 275, 828 [NASA ADS] [CrossRef] [Google Scholar]
  56. Marsh, T. R., Nelemans, G., & Steeghs, D. 2004, MNRAS, 350, 113 [NASA ADS] [CrossRef] [Google Scholar]
  57. Moll, R., & Woosley, S. E. 2013, ApJ, 774, 137 [NASA ADS] [CrossRef] [Google Scholar]
  58. Morton, D. C. 2000, ApJS, 130, 403 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  59. Nelson, L. A., Dubeau, E., & MacCannell, K. A. 2004, ApJ, 616, 1124 [NASA ADS] [CrossRef] [Google Scholar]
  60. Nieva, M. F., & Przybilla, N. 2007, A&A, 467, 295 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Nieva, M. F., & Przybilla, N. 2008, A&A, 481, 199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Nieva, M.-F., & Przybilla, N. 2012, A&A, 539, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Nomoto, K., & Kondo, Y. 1991, ApJ, 367, L19 [NASA ADS] [CrossRef] [Google Scholar]
  64. O’Toole, S. J. 2004, A&A, 423, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. O’Toole, S. J., & Heber, U. 2006, A&A, 452, 579 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Pakmor, R., Kromer, M., Taubenberger, S., & Springel, V. 2013, ApJ, 770, L8 [NASA ADS] [CrossRef] [Google Scholar]
  67. Panei, J. A., Althaus, L. G., Chen, X., & Han, Z. 2007, MNRAS, 382, 779 [NASA ADS] [CrossRef] [Google Scholar]
  68. Parsons, S. G., Marsh, T. R., Gänsicke, B. T., Drake, A. J., & Koester, D. 2011, ApJ, 735, L30 [NASA ADS] [CrossRef] [Google Scholar]
  69. Przybilla, N., Butler, K., Becker, S. R., & Kudritzki, R. P. 2006a, A&A, 445, 1099 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Przybilla, N., Nieva, M. F., & Edelmann, H. 2006b, Balt. Astron., 15, 107 [NASA ADS] [EDP Sciences] [Google Scholar]
  71. Randall, S. K., Fontaine, G., Green, E. M., et al. 2006, ApJ, 643, 1198 [NASA ADS] [CrossRef] [Google Scholar]
  72. Rauch, T., & Deetjen, J. L. 2003, in Stellar Atmosphere Modeling, eds. I. Hubeny, D. Mihalas, & K. Werner, ASP Conf. Ser., 288, 103 [Google Scholar]
  73. Rauch, T., Werner, K., Bohlin, R., & Kruk, J. W. 2013, A&A, 560, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Rauch, T., Rudkowski, A., Kampka, D., et al. 2014, A&A, 566, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Schaffenroth, V. 2015, Ph.D. Thesis, Friedrich-Alexander-Universität Erlangen-Nürnberg [Google Scholar]
  76. Shen, K. J., & Bildsten, L. 2014, ApJ, 785, 61 [NASA ADS] [CrossRef] [Google Scholar]
  77. Silvotti, R., Østensen, R. H., Bloemen, S., et al. 2012, MNRAS, 424, 1752 [NASA ADS] [CrossRef] [Google Scholar]
  78. Sim, S. A., Röpke, F. K., Hillebrandt, W., et al. 2010, ApJ, 714, L52 [NASA ADS] [CrossRef] [Google Scholar]
  79. Steinfadt, J. D. R., Kaplan, D. L., Shporer, A., Bildsten, L., & Howell, S. B. 2010, ApJ, 716, L146 [NASA ADS] [CrossRef] [Google Scholar]
  80. van Kerkwijk, M. H., Bergeron, P., & Kulkarni, S. R. 1996, ApJ, 467, L89 [NASA ADS] [CrossRef] [Google Scholar]
  81. van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Vennes, S., Chayer, P., & Dupuis, J. 2005, ApJ, 622, L121 [NASA ADS] [CrossRef] [Google Scholar]
  83. Vennes, S., Thorstensen, J. R., Kawka, A., et al. 2011, ApJ, 737, L16 [NASA ADS] [CrossRef] [Google Scholar]
  84. Webbink, R. F. 1984, ApJ, 277, 355 [NASA ADS] [CrossRef] [Google Scholar]
  85. Werner, K., Rauch, T., Ringat, E., & Kruk, J. W. 2012, ApJ, 753, L7 [NASA ADS] [CrossRef] [Google Scholar]
  86. Zuckerman, B., Koester, D., Melis, C., Hansen, B. M., & Jura, M. 2007, ApJ, 671, 872 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Additional material

thumbnail Fig. A.1

Entirety of our UV spectra of HD 188112, compared with our best fitting model. The bottom of each panel features the difference between the observation and the model divided by the noise as χ. The portions of the spectrum plotted in grey are regions excluded from the fits. Most of them feature blueshifted IS lines. The IS lines are artificially broadened by the RV correction applied to the individual spectra and show two components due to a time gap in the series of FUV spectra. For this panel, the abundance of Sn was set to the value matching the Sn iii line.

Open with DEXTER

thumbnail Fig. A.1

continued. Here and in the following panels, the Sn abundance is set to the value matching the Sn iv lines. Nickel was included at the abundance determined by the Ni iii lines, thus the Ni ii lines appear too strong and were excluded from the fits (see Sect. 4.3.4). The broad feature at around 1345 Å is an artifact (blemish) from the STIS detector.

Open with DEXTER

thumbnail Fig. A.1

continued. Here starts the NUV2 spectrum. The scaling of the wavelength varies to include the spectral intervals devoid of fitted lines in unique panels. The iron abundance was set to the value derived with the Fe ii lines.

Open with DEXTER

All Tables

Table 1

Abundances of metals detected in the spectra of HD 188112.

All Figures

thumbnail Fig. 1

Radial velocity curve, with measurements from different sources are indicated with different colors (see Sect. 3). The “optical old” refers to values published in Heber et al. (2003) and Edelmann et al. (2005).

Open with DEXTER
In the text
thumbnail Fig. 2

Fit of the Mg ii lines in the NUV1 spectrum for vrot sin i = 7.9 km s-1 and vmicro = 1.1 km s-1. The IS Mg lines are shown in grey and were not included in the fit.

Open with DEXTER
In the text
thumbnail Fig. 3

Best fits for the P iii line in the FUV spectrum as well as for the Ca ii k line in an optical FEROS spectrum.

Open with DEXTER
In the text
thumbnail Fig. 4

Comparison between parts of the FUV spectrum featuring mainly the Si lines and an LTE model spectrum, in blue. The LTE model was computed using the NLTE abundances stated in Table 1. The Si ii, Fe ii, and S ii lines appear to be much stronger than their observed counterparts, while the Si iii lines are nicely reproduced.

Open with DEXTER
In the text
thumbnail Fig. 5

Summary of the determined chemical composition of HD 188112. The top panel shows the abundances relative to the solar ones (Asplund et al. 2009) and the bottom panel the absolute values. Down-pointing triangles indicate upper limits determined for C, N, and O. The stars in the bottom panel indicate abundances predicted by radiative levitation (P. Chayer, priv. comm.). No support is expected for oxygen and magnesium, thus the stars are at the bottom of the plot. The green crosses indicate the average abundances for sdB stars.

Open with DEXTER
In the text
thumbnail Fig. A.1

Entirety of our UV spectra of HD 188112, compared with our best fitting model. The bottom of each panel features the difference between the observation and the model divided by the noise as χ. The portions of the spectrum plotted in grey are regions excluded from the fits. Most of them feature blueshifted IS lines. The IS lines are artificially broadened by the RV correction applied to the individual spectra and show two components due to a time gap in the series of FUV spectra. For this panel, the abundance of Sn was set to the value matching the Sn iii line.

Open with DEXTER
In the text
thumbnail Fig. A.1

continued. Here and in the following panels, the Sn abundance is set to the value matching the Sn iv lines. Nickel was included at the abundance determined by the Ni iii lines, thus the Ni ii lines appear too strong and were excluded from the fits (see Sect. 4.3.4). The broad feature at around 1345 Å is an artifact (blemish) from the STIS detector.

Open with DEXTER
In the text
thumbnail Fig. A.1

continued. Here starts the NUV2 spectrum. The scaling of the wavelength varies to include the spectral intervals devoid of fitted lines in unique panels. The iron abundance was set to the value derived with the Fe ii lines.

Open with DEXTER
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.