author: M. Latour, U. Heber, A. Irrgang, V. Schaffenroth, S. Geier, W. Hillebrandt, F. K. Röpke, S. Taubenberger , M. Kromer, and M. Fink
date: Received 2015 September 25 / Accepted 2015 November 8
title: Quantitative spectral analysis of the sdB star HD 188112:
  A helium-core white dwarf progenitor

HD 188112 is a bright ($V$ = 10.2 mag) hot subdwarf B (sdB) star with a mass too low to ignite core helium burning and is therefore considered a pre-extremely low-mass (ELM) white dwarf (WD). ELM WDs (M $\protect\la$ 0.3 $M_{\rm \odot}$) are He-core objects produced by the evolution of compact binary systems. We present in this paper a detailed abundance analysis of  based on high-resolution Hubble Space Telescope (HST) near- and far-ultraviolet spectroscopy. We also constrain the mass of the star's companion. We use hybrid non-LTE model atmospheres to fit the observed spectral lines, and to derive the abundances of more than a dozen elements and the rotational broadening of metallic lines. We confirm the previous binary system parameters by combining radial velocities measured in our UV spectra with the previously published values. The system has a period of 0.60658584 days and a WD companion with M $\geq$ 0.70 $M_{\rm \odot}$. By assuming a tidally locked rotation combined with the projected rotational velocity ($v$ sin $i$ = 7.9 $\pm$ 0.3 km s$^{-1}$), we constrain the companion mass to be between 0.9 and 1.3 $M_{\rm \odot}$. We further discuss the future evolution of the system as a potential progenitor of an underluminous type Ia supernova. We measure abundances for Mg, Al, Si, P, S, Ca, Ti, Cr, Mn, Fe, Ni, and Zn, and for the trans-iron elements Ga, Sn, and Pb. In addition, we derive upper limits for the C, N, O elements and find  to be strongly depleted in carbon. We find evidence of non-LTE effects on the line strength of some ionic species such as Si ii and Ni ii. The metallic abundances indicate that the star is metal-poor, with an abundance pattern most likely produced by diffusion effects.

# 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 ($\sim$`<!-- -->`{=html}0.5 $M_{\rm \odot}$) 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 ($T_{\rm eff}$ $\ga$ 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 . The subsequent, more rapid evolution during helium shell burning brings the star above the EHB.

Given its fundamental parameters ($T_{\rm eff}$ = 21,500 K, log $g$ = 5.66),  was classified as an sdB; however, it lies at an odd position in the log $g$–$T_{\rm eff}$ diagram, i.e., significantly below the zero-age EHB . This is an indication that the star has too little mass to sustain He-core burning. The low-mass nature of  was confirmed in by comparing the trigonometric distance of the star (from the Hipparcos parallax) with the spectroscopic distance; a mass of $\sim$`<!-- -->`{=html}0.24 $M_{\rm \odot}$ is required for both values to match. Comparing the fundamental parameters of  with models for post red giant branch evolution led to a similar mass of $\sim$`<!-- -->`{=html}0.23 $M_{\rm \odot}$ (see Figure 4 in ). 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) . 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 . It is thus no surprise that  was found to be radial velocity (RV) variable with an orbital period P = 0.606585 d, a semi-amplitude K$_1$ = 188.3 km s$^{-1}$, a systemic velocity $\gamma_0$ = 26.6 km s$^{-1}$, and a miminum companion mass of 0.73 $M_{\rm \odot}$ . 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 , which does not show any deviation from that of a single star.

At the time of the discovery of , very few low-mass He-core objects were known. The first was discovered as the companion of a neutron star in a millisecond pulsar ; another serendipitous discovery followed a few years later, among the first data release of the Sloan Digital Sky Survey (SDSS) . No more than a handful of such objects had been spectroscopically confirmed when a dedicated search for what are now called extremely low-mass (ELM) WDs was undertaken in the SDSS data base (, 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 ($\la$ 0.2 d) are expected to merge within a Hubble time. To date, 38 of these systems have been found (, 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 . 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 . This double-detonation scenario happens at a sub-Chandrasekhar mass and a C/O WD accretor of $\sim$`<!-- -->`{=html}1.0 $M_{\rm \odot}$ is a good candidate for producing a normal SN Ia . If the mass of the companion is much lower than 1 $M_{\rm \odot}$ 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, $T_{\rm eff}$, and surface gravity, log $g$) of the star with theoretical evolutionary tracks for different masses . 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 or distortion effects (such as ellipsoidal variations) in the light curve allow the system's parameters to be constrained . 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),  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 $g$–$T_{\rm eff}$ plane (see, e.g., Fig. 1 in ). This indicates that  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 : KIC 6614501 , SDSS J1625+3632 , SDSS J0815+2309 , and GALEX J0805-1058 . 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 ($v_{\rm rot}$ sin $i$) of metallic lines . The surface rotation velocity of the star ($v_{\rm rot}$) 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  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 ; 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 . 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.

# Observations

Ultraviolet spectra of  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,$^{\rm }$ 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 lines[^1]. This exposure time leads to $\sim$`<!-- -->`{=html}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 Guide[^2] 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.

# 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 and , 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 $\sim$`<!-- -->`{=html}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. . The orbital parameters we obtained are in very good agreement with those stated in ; we have P = 0.60658584 ($\pm$ 5e$-$`<!-- -->`{=html}8) d, K$_1$ = 188.7 $\pm$ 0.2 km s$^{-1}$, and $\gamma_0$ = 26.6 $\pm$ 0.2 km s$^{-1}$. This leads to a mass function $$f(m)=\frac{M_{\rm{comp}}^{3}\rm{sin^3}\it i}{\left (M_{\rm{comp}}+M\right )^{2}} = \frac{PK^{3}}{2\pi G}$$ of $f(m)$ = 0.4221 $M_{\rm \odot}$. By allowing for an eccentric orbit during the fitting procedure, we derive an upper limit of $e \leq$ 0.00006. This low eccentricity value is consistent with a circular orbit.

The parallax of  from the new reduction of the Hipparcos catalogue is 13.6 $\pm$ 1.7 mas . We then used the V (10.22) and B (10.01) magnitudes from the TYCHO 2 catalogue combined with the synthetic flux from a model atmosphere appropriate for  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 $^{+0.075}_{-0.055}$ $M_{\rm \odot}$ and a radius of 0.12 $^{+0.02}_{-0.01}$ $R_{\rm \odot}$. The main uncertainty results from the parallax measurement. The spectroscopic mass indicates a $M_{\rm comp} \geq$ 0.73 $M_{\rm \odot}$ and a most likely mass $M_{\rm comp (52\degr)}$ = 1.21 $M_{\rm \odot}$ .

The latest development in terms of evolutionary sequences for ELM WDs is due to . Their models predict that for masses in the range 0.18$-$`<!-- -->`{=html}0.4 $M_{\rm \odot}$, 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 $g$–$T_{\rm eff}$ diagram (see figs. 3 and 4 of ). 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 is a statistical one, weighting the time spent by models of different masses in each region of the log $g$–$T_{\rm eff}$ plane. This way, they determine a mass of 0.211 $\pm$ 0.0175 $M_{\rm \odot}$ for   as well as a cooling age of 329 $\pm$ 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_{\rm \odot}$ for  leads to $M_{\rm comp} \geq$ 0.70 $M_{\rm \odot}$ and a most probable mass $M_{\rm comp (52\degr)}$ = 1.17 $M_{\rm \odot}$. 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.

# Spectroscopic analysis

## Model atmospheres and synthetic spectra

The model atomspheres used for our analysis are computed following the hybrid non-LTE (NLTE) approach discussed in and ; 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 . The NLTE effects are accounted for in the computation of the atomic population numbers via updated versions of the Detail and Surface codes . 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; ). It has also been tested in the sdB regime . Compared to the fully NLTE approach, the hybrid approach is much less time-consuming in terms of computation. Given the atmospheric parameters of , 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) = $-$`<!-- -->`{=html}5.0 derived in . 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 . This approach is appropriate when blending of lines originating from different species is sparse, as is the case for . The metallic lines are then fitted using a standard $\chi^2$ minimization technique.

## Rotational velocity measurement

As a first step, we measured the rotational velocity ($v_{\rm rot}$ 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 $v_{\rm micro}$ $\la$ 2 km s$^{-1}$; its effect below this value is rather small and did not influence the resulting $\chi^{2}$ much. We measured a rotational broadening of $v_{\rm rot}$ sin $i$ = 7.9 $\pm$ 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. . 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: $$v_{\rm rot}=\frac{2 \pi R}{P}.$$

Combined with the measured $v$ sin $i$, it yields the inclination angle. For a mass of 0.21 $M_{\rm \odot}$, as predicted by theoretical models, we find $i$ = 57$\degr$ $\pm$ 4$\degr$ and a companion mass of 1.0$^{+0.10}_{-0.08}$ $M_{\rm \odot}$.

If we use the slightly higher sdB mass resulting from the Hipparcos parallax (0.245 $M_{\rm \odot}$), we find $i$ = 52$\degr$ $\pm$ 3$\degr$ and a companion mass of 1.2$^{+0.13}_{-0.09}$ $M_{\rm \odot}$[^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 : $T_{\rm eff}$ = 21,500 $\pm$ 500 K, log $g$ = 5.66 $\pm$ 0.06, and log $N$(He)/$N$(H) = $-$`<!-- -->`{=html}5.0. Our grid also included the extreme parameters resulting from the uncertainties on $T_{\rm eff}$ 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 and it is the usual approach for abundance analysis of white dwarfs and hot subdwarfs having similar or higher $T_{\rm eff}$. 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 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. .

|  |  |  |
|:---|:--:|:--:|
| Element | Abundance |  |
| Z | log $N$(Z)/$N$(H) |  |
|  | NLTE | LTE |
| C ii | \<$-$`<!-- -->`{=html}9.6 | $-$ |
| N ii | \<$-$`<!-- -->`{=html}5.0 | $-$ |
| O i | \<$-$`<!-- -->`{=html}6.3 | $-$ |
| Mg ii | $-$`<!-- -->`{=html}6.40 $\pm$ 0.07 | $-$`<!-- -->`{=html}6.20 $\pm$ 0.07 |
| Al iii | $-$`<!-- -->`{=html}7.37 $\pm$ 0.05 | $-$`<!-- -->`{=html}7.39 $\pm$ 0.05 |
| Si ii-iii-iv | $-$`<!-- -->`{=html}7.3 $\pm$ 0.1 | $-$`<!-- -->`{=html}7.6 $\pm$ 0.1 |
| P iii | $-$ | $-$`<!-- -->`{=html}9.2 $\pm$ 0.1 |
| S ii | $-$`<!-- -->`{=html}8.0 $\pm$ 0.1 | $-$`<!-- -->`{=html}8.4 $\pm$ 0.1 |
| Ca ii | $-$`<!-- -->`{=html}7.3 $\pm$ 0.2 | $-$ |
| Ti iii | $-$ | $-$`<!-- -->`{=html}8.30 $\pm$ 0.05 |
| Cr iii | $-$ | $-$`<!-- -->`{=html}7.55 $\pm$ 0.05 |
| Mn iii | $-$ | $-$`<!-- -->`{=html}8.00 $\pm$ 0.05 |
| Fe ii-iii | $-$`<!-- -->`{=html}5.75 $\pm$ 0.15 | $-$`<!-- -->`{=html}6.1 $\pm$ 0.1 |
| Ni iii | $-$`<!-- -->`{=html}6.6 $\pm$ 0.1 | $-$`<!-- -->`{=html}6.7 $\pm$ 0.1 |
| Zn iii | $-$ | $-$`<!-- -->`{=html}7.72 $\pm$ 0.07 |
| Ga ii-iii | $-$ | $-$`<!-- -->`{=html}9.6 $\pm$ 0.1 |
| Sn iii | $-$ | $-$`<!-- -->`{=html}10.6 $^{+0.6}_{-0.1}$ |
| Pb iv | $-$ | $-$`<!-- -->`{=html}10.0 $\pm$ 0.2 |

Abundances of metals detected in the spectra of .

### 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 lines at 1302, 1304, and 1305 Å. They were used to put an upper limit of log $N$(O)/$N$(H) = $-$`<!-- -->`{=html}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  obtained by . The lines at 3995 and 4630 Å were used to put an upper limit of log $N$(N)/$N$(H) = $-$`<!-- -->`{=html}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 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) = $-$`<!-- -->`{=html}9.6. This is six orders of magnitude below the solar value. Carbon has essentially vanished from the surface of the star.

### Mg, Al, Si, P, and S

Only four Mg lines are visible in the NUV spectra. Their abundances were fitted during the determination of the projected rotational velocity (see Fig. ). We also fitted them in the NUV2 spectra, shown in the last panel of Fig. . 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 based on the $\lambda$`<!-- -->`{=html}4481 line visible in the optical spectrum.

The FUV spectrum shows the doublet at 1379, 1384 Å that matches an abundance of log $N$(Al)/$N$(H) = $-$`<!-- -->`{=html}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) = $-$`<!-- -->`{=html}8.0.

One phosphorus line was used to fit: $\lambda$`<!-- -->`{=html}1344.32. Because of the lack of a line list for this ion in Surface, we used the public codes TLUSTY and SYNSPEC[^4] 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 . An NLTE model atmosphere was computed including hydrogen and helium only and as no model atom is available for , we included it in an LTE way through Synspec. The resulting fit can be seen in Fig. .

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) = $-$`<!-- -->`{=html}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 $\lambda\lambda$`<!-- -->`{=html}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. .

### Ca

Calcium is an important element because it has been claimed that all ELM WDs with log $g$ $\la$ 5.9 show Ca in their optical spectra . This was first illustrated by who combined results from different studies of ELM WDs in their figure 6. However, they erroneously indicated that  showed in its optical spectra, confusing it with the feature mentioned in . Nevertheless, since the tendency is true for all other ELM WDs analyzed so far, we went back to the optical spectra of  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. , this line is very weak indeed, but present in the different spectra, thus real. It is understandable why it was first overlooked in . It must also be kept in mind that most of the ELM WDs in which Ca has been measured have $T_{\rm eff}$ \< 16,000 K. At high temperatures, the strength of these lines rapidly decreases as becomes ionized. Nevertheless, the tiny line we found in  indicates an abundance of log $N$(Ca)/$N$(H) = $-$`<!-- -->`{=html}7.3. Calcium was also measured in another ELM WD (SDSS J0815+2309) having essentially the same temperature as  but a slightly higher surface gravity . Compared to , the Ca abundance derived for J0815+2309 is a hundred times higher (log $N$(Ca)/$N$(H) = $-$`<!-- -->`{=html}5.27). As for other ELM WDs, the Ca abundance of  is similar to the lowest abundances measured in the sample of .

### 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 are seen at 1283 and 1287 Å, and two lines of are also visible at 1318 and 1328 Å. The resulting abundances are listed in Table .

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 lines in the FUV and the lines mainly found in the NUV. Compared with the abundance of log $N$()/$N$(H)[^5] = $-$`<!-- -->`{=html}5.87 $\pm$ 0.07, the lines need a slightly *higher* abundance of $-$`<!-- -->`{=html}5.6 $\pm$ 0.1 to be reproduced. Another investigation of UV iron lines in B-type stars using ADS models found a scatter of $\sim$`<!-- -->`{=html}0.25 dex for FUV lines when abundances measured from different lines are compared . The reason for this, possibly related to the ADS models, will be explored further in a future publication (Schaffenroth et al. 2016, in preparation). In any case, this scatter is similar to the discrepancy we found between the FUV and NUV lines of . We thus adopt a mean abundance of $-$`<!-- -->`{=html}5.75 $\pm$ 0.15 for iron, as indicated in Table .

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 lines results in an abundance of log $N$()/$N$(H) = $-$`<!-- -->`{=html}7.25, which is 0.55 dex *below* what is predicted by the main ionization stage (log $N$()/$N$(H) = $-$`<!-- -->`{=html}6.7 when using the ADS models). In other words, at the abundance derived from the lines, the lines are much too strong, as can be seen in Fig. . 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 to ), we built a small grid of NLTE models with various Ni abundances in order to fit the lines. We obtained an abundance of log $N$()/$N$(H) = $-$`<!-- -->`{=html}6.6 $\pm$ 0.1 dex, essentially the same as derived with ADS. This is not surprising since 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 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.

### Trans-iron elements: Ga, Sn, and Pb

We found lines originating from three trans-iron elements, gallium, tin, and lead, in the spectra of . Such heavy elements are also seen in hot subdwarfs as well as hot DA and DO WDs . 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 . The metallic elements in this cool star are accreted from a debris disk surrounding the star. With its temperature of 21,500 K,  is among the coolest WD/sdB stars for which trans-iron elements have been detected.

The Ga abundance of log $N$(Ga)/$N$(H) = $-$`<!-- -->`{=html}9.55 was determined using the $\lambda$`<!-- -->`{=html}1414.4 resonance line of as well as two weak 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 , which indicates an abundance of log $N$(Pb)/$N$(H) = $-$`<!-- -->`{=html}10.0, slightly above the solar value.

We see the resonance line of at 1251.39 Å, as well as the resonance doublet of $\lambda\lambda$ 1314, 1437. The fit of the line gives an abundance of $N$(Sn)/$N$(H) = $-$`<!-- -->`{=html}10.6. However, at this abundance, the lines are barely visible in our synthetic spectra. In order to reproduce them, we need an abundance of $-$`<!-- -->`{=html}8.4, which is more than 2 dex higher than what is indicated by the 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 are considered to be reliable. The line is blended with a line (1251.417 Å), which is however predicted to be rather weak at the Cr abundance determined. Indeed, a fit of the line without the contribution of the feature leads to an abundance only slightly higher, $-$`<!-- -->`{=html}10.5 dex, still far from the $-$`<!-- -->`{=html}8.4 dex needed to reproduce the 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; ) 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 (; T. Rauch, private communication). With NLTE models, the line leads to an abundance of $-$`<!-- -->`{=html}10.0, and the lines indicate $-$`<!-- -->`{=html}8.2 (1314.5 Å) and $-$`<!-- -->`{=html}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 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 has the same ionization energy as and , 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) = $-$`<!-- -->`{=html}10.6$^{+0.6}_{-0.1}$. This is based on our ADS results, but the upper limit includes the value indicated by TMAP.

# Discussion

## 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 , in a consistent NLTE model, ionization equilibrium is shifted to higher stages. In our case, this could explain why the lines are too strong in our models, while the 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 . 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 elements[^6]. The resulting abundances are indicated in Table . 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) = $-$`<!-- -->`{=html}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$()/$N$(H) = $-$`<!-- -->`{=html}7.8 $\pm$ 0.1,
log $N$()/$N$(H) = $-$`<!-- -->`{=html}7.33 $\pm$ 0.05,
and log $N$()/$N$(H) = $-$`<!-- -->`{=html}7.1 $\pm$ 0.1.
The abundance of the main ionization stage () remains the same as derived in NLTE, but the abundance is 0.5 dex lower. The 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. . The spectrum includes the abundances determined with the NLTE models. The lines (as well as few lines and one line) appear too strong, as can be seen with the lines (Fig ).

It thus becomes clear that NLTE effects are present in  and influence the line strengths, especially those of the non-dominant ionic species. For , , and the differences in strength of the modeled lines between LTE and NLTE lead to a $\sim$`<!-- -->`{=html}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 line was visible. In this case the derived abundance might be overestimated. The strongest gallium line is from but the LTE models simultaneously reproduced this line and the two weak features. However, some mystery still remains about the Sn lines: it is not clear what is causing the 2 dex discrepancy between and . 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.

## Metal abundances

Figure summarizes the chemical composition of  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. ). 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.

### Comparison with diffusion models

The bottom panel of Fig. also features abundances predicted by radiative levitation (red stars) for the atomic species included in TOPbase[^7]. These abundances were computed for a model with $T_{\rm eff}$ = 21,000 K and log $g$ = 5.6 according to the method described in (P. Chayer, private communication). 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.

### Comparison with ELM WDs

A UV spectrum has been analyzed for only one other ELM WD so far, namely that of GALEX J1717+6757 . This star has log $g$ = 5.67, essentially identical to that of , and $T_{\rm eff}$ = 14,900 K, which gives a mass of $\sim$`<!-- -->`{=html}0.19 $M_{\rm \odot}$. 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 . In their discussion, 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 , where important radiative support is also predicted for carbon at 21,000 K, but the element is strongly depleted in the stellar atmosphere. argue that a possible explanation for carbon depletion might be a history of repeated CNO flashes, as predicted by evolutionary models with $M \ga$ 0.18 $M_{\rm \odot}$. This might also apply to , 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 . According to Fig. 5 in 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 ($T_{\rm eff}$ = 16,000 K, log $g$ = 4.9), the companion of a massive millisecond pulsar , and SDSS J0745+1949 ($T_{\rm eff}$`<!-- -->`{=html}8̄,380 K, log $g$ = 6.2), a tidally distorted star . 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_{\rm \odot}$) 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, 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$ $\la$ 6.0. , 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) = $-$`<!-- -->`{=html}7.3), though below the detection threshold, is comparable to the lowest abundances measured in their cooler ELM WDs (we note here that the lines get stronger at lower temperature). The same observation holds for Mg, where the abundance of  is similar to the lowest values detected in the sample.

### Comparison with sdB stars

Given its fundamental parameters,  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. the "upper averages" (including determinations and upper limits) of the sample, consisting of a hundred sdBs. It thus appears that  is a particularly metal-poor sdB, but also an especially cool one. PG1627+017 is a cool pulsating sdB ($T_{\rm eff}$ $\sim$`<!-- -->`{=html}21$-$`<!-- -->`{=html}24 kK) that can be compared to  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 . Indeed, a mass of $\sim$`<!-- -->`{=html}0.45 $M_{\rm \odot}$ was suggested by its seismic properties . The UV spectrum of PG1627+017 has allowed a detailed abundance determination, and the star is clearly more metal-rich than  (see Fig. 6 in ). Carbon is the most extremely depleted element in . 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$_{\odot}$)) = $-$`<!-- -->`{=html}5.7 and $-$`<!-- -->`{=html}5.3) were determined for CD $-$`<!-- -->`{=html}24$\degr$`<!-- -->`{=html}731 and PG 1219+534, respectively . The sdOB star Feige 110 is another example of a carbon poor star . Lead is the most enriched element in HD188112, and it has been identified in 22 out of the 27 hot subdwarf UV spectra examined by . Further abundance analysis of sdBs showed that lead enrichment varies from 10 to 1000 times solar . 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  is similar to other sdBs.

## Binary system

In addition to the spectroscopic analysis, we attempted a better characterization of the  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 $\pm$ 0.2 km s$^{-1}$), and systemic velocity (26.6 $\pm$ 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 $^{+0.075}_{-0.055}$ $M_{\rm \odot}$ derived with the parallax distance agrees quite well with that indicated by evolutionary models of low-mass He-core WDs, which predict 0.211 $\pm$ 0.0175 $M_{\rm \odot}$ . 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_{\rm \odot}$. Such a mass is at the lower limit of He-core burning for an sdB star . Low-mass sdB stars are created when a massive progenitor (2$-$`<!-- -->`{=html}2.5 $M_{\rm \odot}$) 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_{\rm \odot}$ in close binaries with a WD companion . The position of  in the log $g$ - $T_{\rm eff}$ diagram would fit within the 0.33 $M_{\rm \odot}$ EHB band (see Fig. 7 from ). Although there is a possibility that  is a He-core burning object, one must keep in mind that these low-mass sdBs occur much less frequently than the canonical ($\sim$`<!-- -->`{=html}0.47 $M_{\rm \odot}$) 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  is such an object. For the lower mass limit on , both the parallax and evolutionary models give a value of $\sim$`<!-- -->`{=html}0.19 $M_{\rm \odot}$. 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_{\rm \odot}$ solutions indicated respectively by evolutionary models and the parallax. We derived a minimum mass of 0.7 $M_{\rm \odot}$ for the companion[^8]. 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 . When considering the distribution for sdB binaries, it appears that a 0.7 $M_{\rm \odot}$ companion is quite massive; only two known sdBs have a more massive WD companion[^9]. After a survey was dedicated to finding massive companions to sdBs , 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 ), the companion masses are evenly distributed between 0.1 and 0.9 $M_{\rm \odot}$ with a few more massive outliers. In this context, the companion of  would be on the massive side, but not outstandingly so.

We measured the rotational broadening of metallic lines in our UV spectra and derived $v_{\rm rot}$ sin $i$ = 7.9 $\pm$ 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$\degr$ and 61$\degr$. 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_{\rm \odot}$. 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 $\sim$`<!-- -->`{=html}0.75 $M_{\rm \odot}$, with $\sigma$ $\sim$`<!-- -->`{=html}0.25 $M_{\rm \odot}$ . 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 and . If the  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$\degr$, 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  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$\degr$ and the companion mass above 6 $M_{\rm \odot}$. A search for possible X-ray emission from the companion of  might help to constrain its nature .

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 . Brief episodes of mass transfer might lead to mass loss of a few times 10$^{-4}$ M$_\odot$ during a shell flash (see Fig. 12 in ). 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$_\odot$.

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$_\odot$, the white dwarf is of C/O composition and the merger most likely produces a RCrB star . Mergers of He-core white dwarfs with more massive C/O white dwarfs may lead to SNe Ia . 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, ). If the core fails to ignite[^10], then the helium explosion alone can produce a fainter type .Ia SN . At the high end of the companion mass range ($>$`<!-- -->`{=html}1.1 M$_\odot$), the system would have a low mass ratio ($q$ $\la$`<!-- -->`{=html}0.20) and thus evolve via stable mass transfer into an AM Canum Venaticorum (AM CVn) , which is a potential SN Ia progenitor . However, it is usually assumed that white dwarfs more massive than about 1.1 M$_\odot$ 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 . On the other hand, recent investigations showed that near-Chandrasekhar mass ($\sim$`<!-- -->`{=html}1.1$-$`<!-- -->`{=html}1.3 $M_{\rm \odot}$) 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 . 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 . For mass ratio $q$ in the range $\sim$`<!-- -->`{=html}0.25–0.70 the mass transfer can be either stable or unstable (see Fig. 1 in ). For the  system, $q$ can be up to $\sim$`<!-- -->`{=html}0.32 if the companion mass is close to the minimum value.

# Conclusion

We began our analysis of the UV spectra of  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,  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 $\sim$`<!-- -->`{=html}0.5 dex. As these effects are seen in , they are likely to also be present in stars having a temperature higher than $\sim$`<!-- -->`{=html}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,  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,  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 $T_{\rm eff}$ decreases and gravitational settling stronger as log $g$ increases. However, given its mass,  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  companion, we are limited by the unknown inclination of the system and by the uncertainty on the mass of  itself. By assuming tidally locked rotation, we could constrain the inclination to be between 49$\degr$ and 61$\degr$. 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_{\rm \odot}$, i.e., a rather massive WD, while the minimum mass is 0.7 $M_{\rm \odot}$. For the mass of , a value up to 0.32 $M_{\rm \odot}$ 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.

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

[^2]: http://www.stsci.edu/hst/stis/documents/isrs/200002.pdf

[^3]: Only the uncertainty on $i$ is considered in these estimates of the companion mass.

[^4]: http://nova.astro.umd.edu/

[^5]: We use the nomenclature here and below to indicate the abundance derived by using lines from a specific ion.

[^6]: 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.

[^7]: http://cdsweb.u-strasbg.fr/topbase/topbase.html

[^8]: If we consider 0.19 $M_{\rm \odot}$ for   the minimum mass of the companion decreases only slightly to 0.68 $M_{\rm \odot}$.

[^9]:  was not included in the sample because of its low mass.

[^10]: 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 .