Issue
A&A
Volume 633, January 2020
Article Number L5
Number of page(s) 9
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/201937343
Published online 10 January 2020

© ESO 2020

1. Introduction

Stellar black holes are hard to find if they are not in interacting binaries, that is, in X-ray binaries. Recently, Giesers et al. (2018) found a single-lined binary in the globular cluster NGC 3201 consisting of a turn-off star orbited by a detached compact object of minimum mass 4.36 ± 0.41 M in a highly eccentric orbit with a period of 166.8 d. Because this minimum mass exceeds the Oppenheimer-Volkoff limit, the discovery represents the first direct mass estimate of a detached black hole in a globular cluster.

Most recently, Liu et al. (2019) reported the discovery of the long-period (single-lined) spectroscopic binary system LB-1 (LS V+22 25). The visible component of this system is a bright B-type star with a prominent and broad Hα emission line very similar to a Be star. A fit of TLUSTY B-star model spectra lead to an effective temperature Teff = 18 104 ± 825 K, a surface gravity of log(g) = 3.43 ± 0.15, and a metallicity consistent with the solar value. Comparing these atmospheric parameters with stellar evolution models of B-type main sequence (MS) stars, the authors concluded that it is a young, massive star with M = 8 . 2 1.2 + 0.9 M $ M=8.2^{+0.9}_{-1.2}\,M_\odot $ and R = 9 ± 2 R located in the direction of the Galactic anti-center at a distance of d = 4.23 ± 0.24 kpc.

The absorption lines of the star show radial-velocity variations with a period P = 78.9 ± 0.3 d. Fitting a binary orbit to the radial-velocity curve yielded a semi-amplitude K1 = 52.8 ± 0.7 km s−1 and a very small eccentricity e = 0.03 ± 0.01. Adopting the above-mentioned mass for the visible B-star, the mass of the companion would have to be higher than 6 . 3 1.0 + 0.4 M $ 6.3^{+0.4}_{-1.0}\,M_\odot $. Since no spectral features of another star are detectable in the spectra, Liu et al. (2019) concluded that the companion cannot be a star, but instead has to be a black hole.

Based on the complex shape and substantial width (∼240 km s−1) of the Hα emission line, the authors concluded that it can only be a Keplerian disk. Since its radial velocity does not follow the motion of the star, but shows a very small variation (Kα = 6.4 ± 0.8 km s−1) in anti-phase, the authors propose that the disk is connected to the black hole companion. If so, the orbital inclination of the binary would be very small (io = 15–18°) and the black hole extremely massive ( M BH = 68 13 + 11 M $ M_{\mathrm{BH}}=68^{+11}_{-13}\,M_\odot $). This result strongly challenges stellar evolution models, which do not predict black holes of such high masses to be formed in metal-rich environments (Liu et al. 2019; Eldridge et al. 2019; Groh et al. 2019; Belczynski et al. 2019). The conclusion of Liu et al. (2019) rests on two main pillars: first, the B-star is a massive MS star; and second, the emission line is connected to the compact companion. However the latter was severely questioned by Abdul-Masih et al. (2019) and El-Badry & Quataert (2019) who argue that the radial velocity variation of the emission line is actually caused by the variability of the underlying absorption line of the visible star and is therefore not related to the compact companion. Because the Hα emission appears to be stationary, it might result from a circumbinary disk. Consequently, the measured mass function

f ( M ) M 2 sin 3 ( i o ) ( 1 + M 1 / M 2 ) 2 = ( 1 e 2 ) 3 / 2 K 1 3 P 2 π G = 1.20 ± 0.05 M , $$ \begin{aligned} f(M) \frac{M_2 \sin ^3(i_\mathrm{o} )}{(1+M_1/M_2)^2} = (1-e^2)^{3/2}\frac{K_1^3 P}{2 \pi G} = 1.20\pm 0.05\,M_\odot, \end{aligned} $$(1)

seems to remain as the only constraint on the mass M2 of the compact companion, which is then a function of two unknowns, namely the mass of the visible star M1 and the orbital inclination io. The spectral analysis by Abdul-Masih et al. (2019) revised the atmospheric parameters to Teff = 13 500 ± 700 K, log(g) = 3.3 ± 0.3, and a projected rotational velocity of v sin(i) = 7.5 ± 4.0 km s−1. Those parameter values might still be consistent with a massive B-type subgiant. However, the very low projected rotational velocity – although it can be explained by seeing the object relatively pole-on – is quite unusual for such an object and could also be indicative of a different nature. Indeed, some evolved low-mass stars are known that can mimic normal B-type stars spectroscopically (see, e.g., Hambly et al. 1997; Ramspeck et al. 2001); these can be distinguished from massive B stars by their slow rotation and chemical anomalies. Therefore, in order to unravel the nature of the visible star and hence its mass, we embarked on a quantitative spectroscopic analysis of LS V+22 25. This analysis shows that this object is very likely a rare stripped helium star in an evolutionary stage, at which it simply mimics a massive MS star. This has important consequences for the mass and potentially also for the nature of its compact companion.

2. Analysis

2.1. Quantitative spectroscopic analysis

Our spectral investigation is based on the publicly available high-resolution Keck/HIRES spectra described by Liu et al. (2019). The analysis strategy and the applied models are explained in detail in Irrgang et al. (2014). Briefly, model computations follow a so-called hybrid approach in which the structure of the atmosphere is computed in local thermodynamic equilibrium (LTE) with ATLAS12 (Kurucz 1996) while departures from LTE are accounted for by applying updated versions of DETAIL and SURFACE (Giddings 1981; Butler & Giddings 1985). The DETAIL code computes population numbers in non-LTE by numerically solving the coupled radiative transfer and statistical equilibrium equations. The SURFACE code uses the resulting departure coefficients and more detailed line-broadening data to compute the final synthetic spectrum. All three codes have been recently updated (Irrgang et al. 2018) by allowing for non-LTE effects on the atmospheric structure as well as by the implementation of the occupation probability formalism (Hubeny et al. 1994) for hydrogen and new Stark broadening tables for hydrogen (Tremblay & Bergeron 2009) and neutral helium (Beauchamp et al. 1997). Moreover, the concept of opacity sampling is consistently used throughout all computational steps, which is important when atmospheres with nonstandard chemical compositions are analyzed. The best-fitting model parameters are found by a simultaneous fit of all available spectra over their entire useful spectral range; see Fig. A.1 for a comparison.

Our atmospheric parameters (Teff = 12 720 ± 260 K, log(g) = 3.00 ± 0.08, see Table A.1) are, within error bars, consistent with those of Abdul-Masih et al. (2019) but quite different from those of Liu et al. (2019). Furthermore, the low values for the projected rotational velocity and for the microturbulence, which were adopted or measured by Liu et al. (2019) and Abdul-Masih et al. (2019), are more or less confirmed here (v sin(i) = 8.7 ± 0.2 km s−1, ξ ≤ 0.1 km s−1). While those values are already quite uncommon for a MS B-type star, the most striking result is revealed by the abundance analysis; see Fig. 1. The star shows an overwhelming signature of material that is processed by hydrogen burning via the CNO bi-cycle, namely extremely strong enrichment in helium and nitrogen at the expense of carbon and oxygen. Additionally, the heavier chemical elements also exhibit pronounced peculiarities. Magnesium and aluminum are considerably depleted while silicon, sulfur, argon, and iron are closer to “standard” values, but still systematically depleted in terms of mass fractions1. Atomic diffusion processes could possibly explain such an abundance pattern, in particular because the observed helium lines show an unusual shape: the observed line wings are too strong while the line cores are too weak (see Fig. A.2). This could be the signature of vertical stratification due to diffusion in the atmosphere, as observed in the spectra of some blue horizontal branch and hot subdwarf B stars (sdBs, see Schneider et al. 2018 and references therein). However, diffusion is unlikely to act in a star with such low gravity (log(g) = 3.00). Therefore, it is more plausible to attribute the mismatch of the helium lines to disk emission filling the line cores, which is the same effect that is observed for the hydrogen Balmer lines (Liu et al. 2019). However, much higher temperatures would be required to see this effect also in helium. An improved understanding of the disk structure is needed to solve this open question. For the time being, we adopt rather conservative systematic uncertainties on temperature and gravity to compensate for this source of uncertainty.

thumbnail Fig. 1.

Differential abundance pattern of LS V+22 25 with respect to the cosmic abundance standard by Nieva & Przybilla (2012) in terms of fractional particle numbers (upper panel) and mass fractions (lower panel). Solar abundances by Asplund et al. (2009) are used for Al, S, and Ar because cosmic abundances are not available for these elements. The error bars are the square roots of the quadratic sums of all given individual 1σ uncertainties. The gray-shaded area marks elements whose abundances show signatures of hydrogen burning via the CNO bi-cycle. To get rid of this effect, the invariant sum of the bi-cycle-catalysts C, N, and O (“CNO”) is also plotted. The abundance pattern of metals heavier than oxygen is nonstandard as well.

2.2. Analysis of the spectral energy distribution

To cross-check our spectroscopic results and to obtain a more comprehensive picture of LS V+22 25, we also fitted its spectral energy distribution using a two-component model consisting of a stellar spectrum plus a blackbody component. The latter was included to empirically account for the observed infrared excess, which might be associated with the disk structure. While the derived blackbody temperature ( T bb = 1140 100 + 120 $ T_{\mathrm{bb}}=1140^{+120}_{-100} $ K) indicates that this thermal source is much cooler than the stellar component, we do not want to overinterpret this result given that the disk is certainly much more complex than a simple blackbody. Nevertheless, by including this component, we are able to reproduce the entire spectral energy distribution from the near-ultraviolet to the far-infrared (see Fig. A.3) yielding results (see Table A.2) that are consistent with spectroscopy. Another notable outcome of this exercise is an estimate for the star’s angular diameter and for the interstellar reddening in its line of sight. Using the extinction law by Fitzpatrick (1999), we obtain an extinction coefficient R V = 3 . 64 0.20 + 0.22 $ R_V = 3.64^{+0.22}_{-0.20} $ and a color excess E(B − V) = 0.397 ± 0.024 mag. Those values become important below, when the reliability of the Gaia parallax is discussed.

3. The nature of LS V+22 25 and its unseen companion

3.1. B-type stars: not only a MS phenomenon

B-type spectra can be observed for stars of very different evolutionary stages. The most common ones are young, massive MS stars. However, the much less massive (∼0.47 M) and smaller (∼0.1–0.3 R) sdBs show relatively similar spectra. Most of them are likely core helium-burning stars located on or somewhat evolved away from the horizontal branch (Dorman et al. 1993; Han et al. 2002). However, in addition, even less massive (∼0.2–0.47 M) helium stars, which are direct progenitors of helium white dwarfs (He-WDs) can show B-type spectra for a period of time (∼106–107 yr; Driebe et al. 1998; Althaus et al. 2013; Istrate et al. 2014) during their evolution (e.g., Heber et al. 2003; Latour et al. 2016). Those objects are closely related to the cooler and less massive extremely low-mass white dwarfs, which have been discovered in substantial numbers in recent years (Brown et al. 2016; Pelisoli & Vos 2019).

Most low-mass B-type stars are much less luminous than their MS counterparts and their higher surface gravities lead to a stronger broadening of the spectral lines. Moreover, the majority of them also show chemical peculiarities in their atmospheres due to earlier mixing episodes, diffusion, isotopic anomalies, or stratification (e.g., O’Toole & Heber 2006; Geier 2013; Schneider et al. 2018). Finally, low-mass B-type stars are usually intrinsically very slow rotators (Geier & Heber 2012). Consequently, such objects can in general be easily distinguished from MS stars.

3.2. Helium stars

One common feature of all sdB types is that they are helium stars, which is difficult to explain using evolutionary scenarios involving only single stars. The helium core of a red giant has to be stripped either well before or at the moment at which helium-burning starts in the core. The stripping is done either by stable Roche lobe overflow (RLOF) to a companion star or via unstable common-envelope ejection (Han et al. 2002, 2003). Depending on the mass-transfer mechanism, very short-period post-common-envelope systems consisting of sdBs with small companions (P ∼ 0.05–30 d, e.g., Kupfer et al. 2015; Schaffenroth et al. 2019) or wide post-RLOF systems (P ∼ 300–1200 d, Vos et al. 2013) are formed.

Only a few candidates for core helium-burning sdB stars have been proposed to have neutron star or black hole companions (Geier et al. 2010), but an extended survey for such objects in close orbits could not confirm their existence (Geier et al. 2011, 2015). However, at least one pre-He-WD with a neutron star companion is known (Kaplan et al. 2013) and the more evolved He-WDs are regularly detected as visible companions of neutron stars.

Stripped helium stars with masses higher than ∼0.47 M are also predicted by binary evolution models (Götberg et al. 2018). However, these are much rarer and only a few of them have been observed so far, such as for example HD 49798 (Mereghetti et al. 2009). This enigmatic X-ray binary subdwarf O star (sdO) may be considered as the prototype for a stripped helium star (Brooks et al. 2017). Although HD 49798 has been assigned the spectral type sdO, it is unique among the hot subdwarfs because it is considerably more massive and luminous. It is a single-lined spectroscopic binary with an orbital period of 1.5 d, which was discovered as a soft X-ray source by the ROSAT satellite. The presence of X-rays indicates that a compact object accretes from the stellar wind. X-ray eclipses allowed dynamical measurement of the masses of the two binary components, revealing that HD 49798 is much more massive (1.5 M) than normal hot sudwarfs stars (∼0.47 M) and that its companion is either a white dwarf or a neutron star of 1.28 ± 0.05 M. HD 49798 has been revisited recently, using both spectroscopy and evolutionary modeling. Krtička et al. (2019) carried out a quantitative spectral analysis of optical and ultraviolet spectra to derive the abundance pattern. Similar to LS V+22 25, the atmospheric helium is strongly enriched and the CNO pattern is obvious. The abundances of heavier elements are also nonstandard: mostly overabundant with respect to the Sun. Brooks et al. (2017) computed the evolution of such a binary and were able to explain the observed abundances in helium, carbon, nitrogen, and oxygen as a result of stripping of the envelope of a 7 M star by its compact companion.

3.3. LS V+22 25: An intermediate-mass stripped helium star?

To constrain the mass, radius, and luminosity of LS V+22 25, we make use of its angular diameter, surface gravity, and parallax from Gaia (ϖ = 0.4403 ± 0.0856 mas), which places the star much closer than expected for a massive B-star. Liu et al. (2019) argued that the Gaia measurement could be erroneous because of the orbital motion of the star. However, the Gaia quality flags indicate that the astrometric solution is well-behaved: The uncertainty is smaller than 20%, 100 out of 103 astrometric observations are marked as good, the “renormalized unit weight error” (RUWE = 0.95, see Lindegren 2018) is perfectly acceptable, and possibly problematic correlations between position and parallax discussed by Liu et al. (2019) were shown to simply be the result of Gaia’s scanning pattern (Eldridge et al. 2019). Moreover, and contrary to Liu et al. (2019), the combination of our derived value for interstellar reddening towards the line-of-sight of LS V+22 25 with three-dimensional dust maps by Green et al. (2019) yields a distance estimate that is fully consistent with the parallactic distance; see Fig. A.4. Assuming the Gaia parallax to be correct, we can combine it with the photometrically derived angular diameter and the surface gravity from spectroscopy to determine the radius, mass, and luminosity of LS V+22 25 (see Table A.2), yielding R1 = 5.5 ± 1.1 R, M1 = 1.1 ± 0.5 M, and log(L1/L) = 2.95 ± 0.19. If LS V+2 25 was indeed a stripped helium star of that mass, it should be observed as a helium-rich sdO star (Götberg et al. 2018). However, and as already mentioned by Liu et al. (2019), the star could be caught in the short-period post-stripping transition phase during which it evolves on the Kelvin-Helmholtz timescale, which is of the order of 104 yr. The unseen companion could then have a mass as low as 2 . 5 0.5 + 0.4 M $ 2.5^{+0.4}_{-0.5}\,M_\odot $ and might therefore be either a black hole, a (massive) neutron star, or even a relatively unevolved low-mass MS star that is too faint to be detected (see Fig. A.5). For orbital inclinations smaller than 45°, the mass of the compact companion would exceed 5 . 0 0.7 + 0.6 M $ 5.0^{+0.6}_{-0.7}\,M_\odot $ and therefore reach values that are quite typical for known black holes in X-ray binaries.

If this were the correct scenario, we would see a stripped star in a very short-lived relaxation phase just after the stripping event. The presence of a disk around the companion consisting of the remnants of the envelope of the red-giant progenitor might then be the relic of mass transfer. Also, the almost perfectly circular orbit of the binary would fit to a previous episode of mass transfer which circularized the rather wide orbit.

3.4. LS V+22 25: A low-mass stripped pre-He-WD?

As outlined by Liu et al. (2019), the orbital motion in the binary system may have a non-negligible impact on the Gaia measurements. Based on the orbital parameters by Liu et al. (2019), the projected semimajor axis of the binary system is

a 1 sin ( i o ) = 1 2 π ( 1 e 2 ) 1 / 2 K 1 P = 0.383 ± 0.006 AU . $$ \begin{aligned} a_1 \sin (i_{\mathrm{o} }) = \frac{1}{2\pi } (1-e^2)^{1/2} K_1 P = 0.383\pm 0.006\,\mathrm{AU}. \end{aligned} $$(2)

Depending on the orientation of the orbital plane (e.g., the orbital inclination) and the scanning pattern of Gaia (e.g., the number of visits and their sampling of the binary orbit), this intrinsic binary motion may be comparable to Earth’s movement of 1 AU, and thus large enough to disturb the parallax measurement without producing a strong signal in the quality flags. If we thus allow LS V+22 25 to be even closer than suggested by the parallax and the reddening maps, another scenario may be considered, which involves a low-mass stripped pre-He-WD.

LS V+22 25 would not be the first pre-He-WD mimicking a MS B-type star. Already decades ago, a similar case was discussed (Luyten 1965). The star HZ 22 (also known as UX CVn) resembled an early-type MS star. However, variations in the radial velocity and in the light curve due to ellipsoidal deformation showed that it must be a single-lined spectroscopic binary system with a massive white dwarf companion. The resulting orbital period of about half a day (Young et al. 1972) was too short for the system to contain a normal MS B-star. Based on evolutionary models, HZ 22 was identified as a helium-core object of ∼0.39 M (Trimble 1973; Schönberner 1978). Shimanskii (2002) performed a detailed quantitative abundance analysis and found that helium, carbon, and iron are depleted, which is consistent with an old evolved star, that is, with the interpretation that this object is a proto-helium-core white dwarf.

Figure A.6 shows a (Teff, log(g)) diagram with MS evolutionary tracks (Ekström et al. 2012) as well as pre-He-WD tracks (Driebe et al. 1998). Therein, tracks for very low-mass pre-He-WDs can be clearly seen to overlap with the MS tracks of massive B-type stars. LS V+22 25 is located close to tracks with masses of ∼0.315 M and might be in a slightly earlier evolutionary phase than the prototype HZ 22, which is hotter and more compact. The evolutionary timescale of such a contracting object until it crosses the zero-age main sequence is about 105 yr (Driebe et al. 1998). The respective stellar radius would still be as large as ∼2.6 R. The minimum mass of the compact companion could be as low as 1.69 ± 0.04 M, which would be consistent with a neutron star or a black hole. A low-mass MS companion is less likely in this scenario because it would probably be luminous enough to leave imprints in the observed spectrum (see Fig. A.5), and these are not seen in the currently available data. For orbital inclinations smaller than 45°, the mass of the compact companion would exceed 3 . 95 0.09 + 0.10 M $ 3.95^{+0.10}_{-0.09}\,M_\odot $ which would rule out a neutron star companion. However, the spectroscopic distance of such a pre-He-WD would then be just 1.21 ± 0.12 kpc, that is, about half the Gaia and reddening distances.

4. Summary

We carried out a quantitative spectral analysis of LS V+22 25 and derive an effective temperature and a surface gravity that are both lower than previously reported. The inferred abundance pattern is characterized by the signature of hydrogen burning via the CNO bi-cycle. Moreover, heavier elements (magnesium, aluminum, silicon, sulfur, argon, and iron) deviate more or less from values expected for normal B-type MS stars. We also find a good fit of the spectral energy distribution based on a stellar model plus a blackbody component that was used as a simple proxy for the circumstellar disk. Making use of the Gaia parallax, the spectrophotometric angular diameter, and the spectroscopic surface gravity, we derive a stellar mass of 1.1 ± 0.5 M, which is much lower than that of a B subgiant. Together with the peculiar abundance pattern, this suggests that LS V+22 25 is not a normal MS B-type star but rather a stripped helium star.

Motivated by two prototypical low-mass hot stars, the sdO X-ray binary HD 49798 and the B + WD binary HZ 22, we considered two classes of stripped helium stars. The former scenario envisages envelope stripping of a massive star by a compact companion while the latter assumes that a low- to intermediate-mass red giant is stripped before the onset of the helium flash, leaving behind an inert helium star of low mass (<0.5 M). We consider the latter scenario unlikely because it is at variance with the distance determined from parallax and reddening. Moreover, the expected lifetimes are of the order of 105 to 106 years, which might be too long to sustain the observed disk if it is indeed circumbinary as suggested by El-Badry & Quataert (2019). The intermediate-mass scenario on the other hand is consistent with the distance estimates as well as with a young age as indicated by the star’s position in the Perseus spiral arm (see Fig. A.7), but would be in a rather short-lived evolutionary phase (104 yr), which might explain the presence of the circumstellar disk, that is, as a relic of the mass transfer. The abundances derived here will serve as crucial constraints for tailored evolutionary models.

Our mass estimate for LS V+22 25 is significantly smaller than that proposed by Liu et al. (2019). Consequently, for a wide range of orbital inclinations, the mass of the compact companion is considerably lower than the claimed value of M BH = 68 13 + 11 M $ M_{\mathrm{BH}}=68^{+11}_{-13}\,M_\odot $ and therefore not in conflict with the general picture of stellar evolution in our Galaxy. For inclinations close to 90°, the compact companion might not even be a black hole but a neutron star or a relatively unevolved low-mass MS star.

A possible progenitor system of LS V+22 25 has already been discovered recently: 2MASS J05215658+4359220 consists of a bright, rapidly rotating red giant star orbited by a noninteracting compact companion with a mass of 3 . 3 0.7 + 2.8 M $ 3.3_{-0.7}^{+2.8}\,M_\odot $, which is probably a neutron star or a black hole (Thompson et al. 2019). Just like LS V+22 25, the orbit is almost circular, but the most stunning similarity is the orbital period of ∼83 d. The object 2MASS J05215658+4359220 likely shows ellipsoidal variations and is close to filling its Roche lobe. This red giant will soon get stripped by the compact companion and might evolve to a configuration that resembles that of LS V+22 25.


1

We focus here on mass instead of number fractions because the former are invariant against stellar evolution, i.e., they do not change when hydrogen is converted to helium. Consequently, they are the natural measure for chemical peculiarities when stellar evolution is at play.

Acknowledgments

We thank the anonymous referee for a very quick and constructive report, Philipp Podsiadlowski for helpful discussions, Danny Lennon for emphasizing the importance of mass fractions, and John E. Davis for the development of the SLXFIG module used to prepare the figures in this paper. This research has made use of the Keck Observatory Archive (KOA), which is operated by the W. M. Keck Observatory and the NASA Exoplanet Science Institute (NExScI), under contract with the National Aeronautics and Space Administration. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.

References

  1. Abdul-Masih, M., Banyard, G., Bodensteiner, J., et al. 2019, ArXiv e-prints [arXiv:1912.04092] [Google Scholar]
  2. Althaus, L. G., Miller Bertolami, M. M., & Córsico, A. H. 2013, A&A, 557, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  4. Beauchamp, A., Wesemael, F., & Bergeron, P. 1997, ApJS, 108, 559 [NASA ADS] [CrossRef] [Google Scholar]
  5. Belczynski, K., Hirschi, R., Kaiser, E. A., et al. 2019, ApJ, submitted [arXiv:1911.12357] [Google Scholar]
  6. Bessell, M., & Murphy, S. 2012, PASP, 124, 140 [Google Scholar]
  7. Brooks, J., Kupfer, T., & Bildsten, L. 2017, ApJ, 847, 78 [NASA ADS] [CrossRef] [Google Scholar]
  8. Brown, W. R., Gianninas, A., Kilic, M., Kenyon, S. J., & Allende Prieto, C. 2016, ApJ, 818, 155 [NASA ADS] [CrossRef] [Google Scholar]
  9. Butler, K., & Giddings, J. R. 1985, Newsletter of Analysis of Astronomical Spectra, No. 9 (London: Univ. London) [Google Scholar]
  10. Cutri, R. M., Wright, E. L., Conrow, T., et al. 2014, VizieR Online Data Catalog: II/328 [Google Scholar]
  11. Dorman, B., Rood, R. T., & O’Connell, R. W. 1993, ApJ, 419, 596 [CrossRef] [Google Scholar]
  12. Driebe, T., Schoenberner, D., Bloecker, T., & Herwig, F. 1998, A&A, 339, 123 [Google Scholar]
  13. Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. El-Badry, K., & Quataert, E. 2019, MNRAS, submitted [arXiv:1912.04185] [Google Scholar]
  15. Eldridge, J. J., Stanway, E. R., Breivik, K., et al. 2019, MNRAS, submitted [arXiv:1912.03599] [Google Scholar]
  16. Evans, D. W., Riello, M., De Angeli, F., et al. 2018, A&A, 616, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
  18. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Geier, S. 2013, A&A, 549, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Geier, S., & Heber, U. 2012, A&A, 543, A149 [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. Geier, S., Kupfer, T., Heber, U., et al. 2015, A&A, 577, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Giddings, J. R. 1981, PhD Thesis, Univ. London [Google Scholar]
  25. Giesers, B., Dreizler, S., Husser, T.-O., et al. 2018, MNRAS, 475, L15 [NASA ADS] [CrossRef] [Google Scholar]
  26. Götberg, Y., de Mink, S. E., Groh, J. H., et al. 2018, A&A, 615, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Green, G. M., Schlafly, E., Zucker, C., Speagle, J. S., & Finkbeiner, D. 2019, ApJ, 887, 93 [NASA ADS] [CrossRef] [Google Scholar]
  28. Groh, J. H., Farrell, E., Meynet, G., et al. 2019, ApJ, submitted [arXiv:1912.00994] [Google Scholar]
  29. Hambly, N. C., Rolleston, W. R. J., Keenan, F. P., Dufton, P. L., & Saffer, R. A. 1997, ApJS, 111, 419 [NASA ADS] [CrossRef] [Google Scholar]
  30. Han, Z., Podsiadlowski, P., Maxted, P. F. L., Marsh, T. R., & Ivanova, N. 2002, MNRAS, 336, 449 [NASA ADS] [CrossRef] [Google Scholar]
  31. Han, Z., Podsiadlowski, P., Maxted, P. F. L., & Marsh, T. R. 2003, MNRAS, 341, 669 [NASA ADS] [CrossRef] [Google Scholar]
  32. Heber, U., Edelmann, H., Lisker, T., & Napiwotzki, R. 2003, A&A, 411, L477 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Henden, A. A., Levine, S., Terrell, D., & Welch, D. L. 2015, Amer. Astron. Soc. Meet. Abstr., 225, 336.16 [Google Scholar]
  34. Hou, L. G., & Han, J. L. 2014, A&A, 569, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Hubeny, I., Hummer, D. G., & Lanz, T. 1994, A&A, 282, 151 [NASA ADS] [Google Scholar]
  36. Irrgang, A., Wilcox, B., Tucker, E., & Schiefelbein, L. 2013, A&A, 549, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Irrgang, A., Przybilla, N., Heber, U., et al. 2014, A&A, 565, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Irrgang, A., Kreuzer, S., Heber, U., & Brown, W. 2018, A&A, 615, L5 [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. Krtička, J., Janík, J., Krtičková, I., et al. 2019, A&A, 631, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Kupfer, T., Geier, S., Heber, U., et al. 2015, A&A, 576, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Kurucz, R. L. 1996, in Model Atmospheres and Spectrum Synthesis, eds. S. J. Adelman, F. Kupka, & W. W. Weiss (San Francisco: ASP), 160 [Google Scholar]
  44. Latour, M., Heber, U., Irrgang, A., et al. 2016, A&A, 585, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Lindegren, L. 2018, Re-normalising the Astrometric Chi-square in Gaia DR2, GAIA-C3-TN-LU-LL-124, www.rssd.esa.int/doc_fetch.php?id=3757412 [Google Scholar]
  46. Liu, J., Zhang, H., Howard, A. W., et al. 2019, Nature, 575, 618 [NASA ADS] [CrossRef] [Google Scholar]
  47. Luyten, W. J. 1965, Proceedings of the 1st Conference on Faint Blue Stars (Minneapolis, MN: Univ. Minnesota) [Google Scholar]
  48. Maíz Apellániz, J., & Weiler, M. 2018, A&A, 619, A180 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Mereghetti, S., Tiengo, A., Esposito, P., et al. 2009, Science, 325, 1222 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  50. Mermilliod, J. C. 2006, VizieR Online Data Catalog: II/168 [Google Scholar]
  51. Nieva, M.-F., & Przybilla, N. 2012, A&A, 539, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. O’Toole, S. J., & Heber, U. 2006, A&A, 452, 579 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Pauli, E. M., Napiwotzki, R., Heber, U., Altmann, M., & Odenkirchen, M. 2006, A&A, 447, 173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Pelisoli, I., & Vos, J. 2019, MNRAS, 488, 2892 [NASA ADS] [CrossRef] [Google Scholar]
  55. Ramspeck, M., Heber, U., & Edelmann, H. 2001, A&A, 379, 235 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Schaffenroth, V., Barlow, B. N., Geier, S., et al. 2019, A&A, 630, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Schneider, D., Irrgang, A., Heber, U., Nieva, M. F., & Przybilla, N. 2018, A&A, 618, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Schönberner, D. 1978, A&A, 70, 451 [NASA ADS] [Google Scholar]
  59. Shimanskii, V. V. 2002, Astron. Rep., 46, 127 [NASA ADS] [CrossRef] [Google Scholar]
  60. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  61. Thompson, T. A., Kochanek, C. S., Stanek, K. Z., et al. 2019, Science, 366, 637 [NASA ADS] [CrossRef] [Google Scholar]
  62. Tremblay, P. E., & Bergeron, P. 2009, ApJ, 696, 1755 [NASA ADS] [CrossRef] [Google Scholar]
  63. Trimble, V. 1973, A&A, 23, 281 [NASA ADS] [Google Scholar]
  64. Vos, J., Østensen, R. H., Németh, P., et al. 2013, A&A, 559, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Young, A., Nelson, B., & Mielbrecht, R. 1972, ApJ, 174, 27 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Additional material

thumbnail Fig. A.1.

Comparison of best-fitting model spectrum (red line) with re-normalized observation (black line; HIRES spectrum taken on December 12, 2017). Light colors mark regions that have been excluded from fitting, e.g., due to data reduction artifacts or the presence of features that are not properly included in our models such as disk emission in the cores of the hydrogen and helium lines or possible stratification effects in the helium lines. Residuals χ are shown as well.

Table A.1.

Atmospheric parameters and abundances of LS V+22 25.

thumbnail Fig. A.2.

Comparison of best-fitting model spectrum (red line) with re-normalized observation (black line; HIRES spectrum taken on December 24, 2017) with special focus on the broad He I lines 4387.93 Å (left panel) and 4471.48 Å (right panel). Light colors mark regions that have been excluded from fitting. Residuals χ are shown as well. The resulting residual pattern, namely that the observed wings are too strong while the observed cores are too weak, could be indicative of vertical stratification; see, e.g., Schneider et al. (2018) for details, or disk emission filling the line cores.

Table A.2.

Stellar parameters derived from photometry and astrometry.

thumbnail Fig. A.3.

Comparison of synthetic and observed photometry: Top panel: spectral energy distribution. The colored data points are filter-averaged fluxes which were converted from observed magnitudes (the respective filter widths are indicated by the dashed horizontal lines), while the gray solid line represents the best-fitting model, i.e., it is based on the parameters from Table A.2, degraded to a spectral resolution of 6 Å. The individual contributions of the stellar (light blue) and blackbody (light red) component are shown as well. The panels at the bottom and on the side show the residuals χ, i.e., the difference between synthetic and observed data divided by the corresponding uncertainties, for magnitudes and colors, respectively. The photometric systems have the following color code: Johnson-Cousins (blue; Mermilliod 2006; Henden et al. 2015; calibration from Bessell & Murphy 2012); SDSS (golden; Henden et al. 2015); Gaia (cyan; Evans et al. 2018 with corrections and calibrations from Maíz Apellániz & Weiler 2018); 2MASS (red; Skrutskie et al. 2006); WISE (magenta; Cutri et al. 2014).

thumbnail Fig. A.4.

Interstellar reddening vs. distance for the line of sight towards LS V+22 25: The black solid line is the best-fitting sample of the “Bayestar19” dust map (Green et al. 2019; http://argonaut.skymaps.info/). The shaded areas illustrate 1σ uncertainty intervals resulting from the star’s Gaia parallax (ϖ = 0.4403 ± 0.0856 mas; vertical strip) and its measured color excess (the value of E(B − V) from Table A.2 has been converted to E(g − r); horizontal bar).

thumbnail Fig. A.5.

Mass of the secondary component as a function of the orbital inclination io: primary masses M1 of 1.1 ± 0.5 M (left panel; see Table A.2) and 0.315 ± 0.015 M (right panel; see Fig. A.6) are used to numerically solve Eq. (1) for M2. The width of the gray shaded region reflects all 1σ-uncertainties. The y-axis on the right-hand-side shows the luminosity of a zero-age main sequence (ZAMS) star with mass M2 based on evolutionary tracks for nonrotating stars of solar metallicity (Ekström et al. 2012). The light-blue shaded horizontal bar represents the inferred luminosity of the visible B-type primary, log(L1/L). For high orbital inclinations, a relatively unevolved low-mass MS companion could be outshined by its primary by more than one (right panel) to almost two (left panel) orders of magnitude, which could be faint enough to avoid imprints in the observed spectrum.

thumbnail Fig. A.6.

Position of LS V+22 25 (blue error bars) in the (Teff, log(g)) diagram: the black solid lines are pre-He-WD tracks from Driebe et al. (1998) labeled with their respective masses. The position of HZ 22 (red error bars; Shimanskii 2002), evolutionary tracks for nonrotating stars of solar metallicity (gray dashed lines; Ekström et al. 2012), and the corresponding ZAMS are shown for reference.

thumbnail Fig. A.7.

Three-dimensional orbit of LS V+22 25 in a Galactic Cartesian coordinate system in which the Galactic center (marked by a black +) lies at the origin, the Sun (marked by a black ⊙) is located on the negative x-axis, and the z-axis points to the Galactic north pole. The nine trajectories (red lines; arrows indicate the star’s current position) illustrate the effects of uncertainties in the parallax and proper motions from Gaia, and in the systemic radial velocity from Liu et al. (2019). These were computed back in time for 200 Myr using a standard three-component, axisymmetric model for the Galactic gravitational potential (see Irrgang et al. 2013 for details on the Milky Way mass model and on the orbit computations). The shape of the orbit (prograde, almost circular, small vertical oscillations) is typical for a thin-disk star (see, e.g., Pauli et al. 2006 for details on the properties of different kinematic groups). The thick blue solid lines schematically represent the loci of the spiral arms based on the polynomial logarithmic arm model of Hou & Han (2014). Interestingly, the current position of the star seems to be aligned with the Perseus spiral arm, which could be a possible sign of youth.

All Tables

Table A.1.

Atmospheric parameters and abundances of LS V+22 25.

Table A.2.

Stellar parameters derived from photometry and astrometry.

All Figures

thumbnail Fig. 1.

Differential abundance pattern of LS V+22 25 with respect to the cosmic abundance standard by Nieva & Przybilla (2012) in terms of fractional particle numbers (upper panel) and mass fractions (lower panel). Solar abundances by Asplund et al. (2009) are used for Al, S, and Ar because cosmic abundances are not available for these elements. The error bars are the square roots of the quadratic sums of all given individual 1σ uncertainties. The gray-shaded area marks elements whose abundances show signatures of hydrogen burning via the CNO bi-cycle. To get rid of this effect, the invariant sum of the bi-cycle-catalysts C, N, and O (“CNO”) is also plotted. The abundance pattern of metals heavier than oxygen is nonstandard as well.

In the text
thumbnail Fig. A.1.

Comparison of best-fitting model spectrum (red line) with re-normalized observation (black line; HIRES spectrum taken on December 12, 2017). Light colors mark regions that have been excluded from fitting, e.g., due to data reduction artifacts or the presence of features that are not properly included in our models such as disk emission in the cores of the hydrogen and helium lines or possible stratification effects in the helium lines. Residuals χ are shown as well.

In the text
thumbnail Fig. A.2.

Comparison of best-fitting model spectrum (red line) with re-normalized observation (black line; HIRES spectrum taken on December 24, 2017) with special focus on the broad He I lines 4387.93 Å (left panel) and 4471.48 Å (right panel). Light colors mark regions that have been excluded from fitting. Residuals χ are shown as well. The resulting residual pattern, namely that the observed wings are too strong while the observed cores are too weak, could be indicative of vertical stratification; see, e.g., Schneider et al. (2018) for details, or disk emission filling the line cores.

In the text
thumbnail Fig. A.3.

Comparison of synthetic and observed photometry: Top panel: spectral energy distribution. The colored data points are filter-averaged fluxes which were converted from observed magnitudes (the respective filter widths are indicated by the dashed horizontal lines), while the gray solid line represents the best-fitting model, i.e., it is based on the parameters from Table A.2, degraded to a spectral resolution of 6 Å. The individual contributions of the stellar (light blue) and blackbody (light red) component are shown as well. The panels at the bottom and on the side show the residuals χ, i.e., the difference between synthetic and observed data divided by the corresponding uncertainties, for magnitudes and colors, respectively. The photometric systems have the following color code: Johnson-Cousins (blue; Mermilliod 2006; Henden et al. 2015; calibration from Bessell & Murphy 2012); SDSS (golden; Henden et al. 2015); Gaia (cyan; Evans et al. 2018 with corrections and calibrations from Maíz Apellániz & Weiler 2018); 2MASS (red; Skrutskie et al. 2006); WISE (magenta; Cutri et al. 2014).

In the text
thumbnail Fig. A.4.

Interstellar reddening vs. distance for the line of sight towards LS V+22 25: The black solid line is the best-fitting sample of the “Bayestar19” dust map (Green et al. 2019; http://argonaut.skymaps.info/). The shaded areas illustrate 1σ uncertainty intervals resulting from the star’s Gaia parallax (ϖ = 0.4403 ± 0.0856 mas; vertical strip) and its measured color excess (the value of E(B − V) from Table A.2 has been converted to E(g − r); horizontal bar).

In the text
thumbnail Fig. A.5.

Mass of the secondary component as a function of the orbital inclination io: primary masses M1 of 1.1 ± 0.5 M (left panel; see Table A.2) and 0.315 ± 0.015 M (right panel; see Fig. A.6) are used to numerically solve Eq. (1) for M2. The width of the gray shaded region reflects all 1σ-uncertainties. The y-axis on the right-hand-side shows the luminosity of a zero-age main sequence (ZAMS) star with mass M2 based on evolutionary tracks for nonrotating stars of solar metallicity (Ekström et al. 2012). The light-blue shaded horizontal bar represents the inferred luminosity of the visible B-type primary, log(L1/L). For high orbital inclinations, a relatively unevolved low-mass MS companion could be outshined by its primary by more than one (right panel) to almost two (left panel) orders of magnitude, which could be faint enough to avoid imprints in the observed spectrum.

In the text
thumbnail Fig. A.6.

Position of LS V+22 25 (blue error bars) in the (Teff, log(g)) diagram: the black solid lines are pre-He-WD tracks from Driebe et al. (1998) labeled with their respective masses. The position of HZ 22 (red error bars; Shimanskii 2002), evolutionary tracks for nonrotating stars of solar metallicity (gray dashed lines; Ekström et al. 2012), and the corresponding ZAMS are shown for reference.

In the text
thumbnail Fig. A.7.

Three-dimensional orbit of LS V+22 25 in a Galactic Cartesian coordinate system in which the Galactic center (marked by a black +) lies at the origin, the Sun (marked by a black ⊙) is located on the negative x-axis, and the z-axis points to the Galactic north pole. The nine trajectories (red lines; arrows indicate the star’s current position) illustrate the effects of uncertainties in the parallax and proper motions from Gaia, and in the systemic radial velocity from Liu et al. (2019). These were computed back in time for 200 Myr using a standard three-component, axisymmetric model for the Galactic gravitational potential (see Irrgang et al. 2013 for details on the Milky Way mass model and on the orbit computations). The shape of the orbit (prograde, almost circular, small vertical oscillations) is typical for a thin-disk star (see, e.g., Pauli et al. 2006 for details on the properties of different kinematic groups). The thick blue solid lines schematically represent the loci of the spiral arms based on the polynomial logarithmic arm model of Hou & Han (2014). Interestingly, the current position of the star seems to be aligned with the Perseus spiral arm, which could be a possible sign of youth.

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.