Open Access
Issue
A&A
Volume 690, October 2024
Article Number A368
Number of page(s) 21
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202450778
Published online 22 October 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Hot subdwarf stars (sdO/Bs) have been discovered as a prominent population of faint blue stars at high Galactic latitudes, but they were subsequently found in all Galactic populations (e.g. Greenstein & Sargent 1974; Luo et al. 2020). They are much smaller and of lower mass than hot main sequence stars of similar spectral types and most of them are likely core-helium-burning extreme horizontal branch (EHB) stars (Heber 1986) with very thin hydrogen envelopes. State-of-the art in hot subdwarf research has been reviewed by Heber (2009, 2016).

The origin of sdO/B stars is still not fully understood, because they can only be formed if the progenitor loses its envelope almost entirely after passing the red-giant branch. While single-star scenarios are discussed, binary interactions are more likely (Pelisoli et al. 2020). Systematic surveys for radial-velocity (RV) variable stars revealed that about one third of the hot subdwarf stars are members of short-period, single-lined binaries (Maxted et al. 2001; see Geier et al. 2022 for a review) with white dwarf (WD), M-type main sequence or brown dwarf companions (Schaffenroth et al. 2022 and references therein). Those are formed after a common envelope phase, during which the companion becomes completely immersed in the red-giant envelope. Stable mass transfer to a main sequence companion in a wide binary and mergers involving He-WDs have been proposed as possible formation channels as well (Han et al. 2002, 2003; Chen et al. 2013; Zhang & Jeffery 2012). Double-lined hot subdwarf binaries with FGK-type main-sequence companions have indeed been discovered to be wide (e.g. Vos et al. 2018). Merger candidates among the single sdO/Bs have also been found (e.g. Dorsch et al. 2022; Werner et al. 2022).

Since most sdB stars are regarded as the progeny of low-mass stars, they can belong to all Galactic populations: the young thin disc, the older thick disc, or even the Galactic halo. Differences between these populations in metallicity and age are predicted to significantly affect hot subdwarf formation (e.g. Han 2008; Vos et al. 2020). Because the abundances in sdB atmospheres are altered by diffusion processes (e.g. Geier 2013), their population membership cannot be deduced from their abundance patterns, but only from their kinematic properties. While disc stars orbit the Galactic centre with moderate velocities and eccentricities, halo stars may have retrograde, very eccentric and highly inclined trajectories with respect to the Galactic plane. The kinematic method has been used to determine the population membership of rather bright samples of field sdBs (de Boer et al. 1997; Altmann et al. 2004; Kawka et al. 2015; Heber 2016; Martin et al. 2017; Bobylev & Bajkova 2019; Luo et al. 2020, 2021). However, the early studies were limited by the poor quality of the tangential velocities and the lack of information about the variability of their radial velocities. While Gaia measurements have remedied the first issue, the second one is still a drawback, even for the most recent studies. Previous studies concluded that most field sdO/Bs are members of the Galactic disc. The < 10% candidates for halo sdO/Bs can have very high Galactic restframe velocities close to the escape velocity of the Galaxy.

However, there are alternative ways to accelerate hot subdwarf stars to (very) high velocities, which were studied intensively after the discovery of the hypervelocity star (HVS) US 708; this is an He-rich sdO (He-sdO) with a velocity high enough to be unbound to the Galaxy, thereby we exclude an origin in the bound halo population (Hirsch et al. 2005). Perets (2009) proposed the disruption of a hierarchical triple system by the supermassive black hole in the centre of the Galaxy (a variant of the slingshot scenario proposed by Hills 1988) and the subsequent merger of an ejected close, helium white dwarf (He-WD) binary as a possible yet quite complicated scenario for US 708. Geier et al. (2015a), however, were able to exclude a Galactic centre origin of this star, calling for a different ejection scenario.

Possible alternatives are related to the supernova ejection scenario (Blaauw 1961), where a close binary is disrupted by the explosion of one component as supernova and the companion is ejected. In principle, core-collapse supernovae can eject He-stars (Renzo et al. 2019; Tauris 2015), but much higher velocities are possible if thermonuclear supernovae are involved. Close hot subdwarf binaries with massive WD companions turned out to be candidates for the progenitors of such thermonuclear type Ia supernovae (Maxted et al. 2000; Geier et al. 2007; Pelisoli et al. 2021) for the double-degenerate merger scenario (Iben & Tutukov 1984; Webbink 1984). Shen et al. (2018) proposed a scenario according to which the WD companion survives the SN Ia (see also Pakmor et al. 2021) and discovered several candidates with extremely high velocities (1000 − 3000 km s−1). These objects are more luminous than normal WDs, resemble hot subdwarfs of lower temperatures, and might remain in such a puffed-up non-equilibrium state for quite some time after the explosion (Bauer et al. 2019; Liu et al. 2021).

Also the exploding WDs themselves might survive the explosion and end up as fast moving objects with very peculiar abundance patterns (e.g. Jordan et al. 2012; Bravo et al. 2016). Several fast candidates with such very peculiar surface compositions have been discovered (Vennes et al. 2017; Raddi et al. 2018a,b, 2019; Gänsicke et al. 2020; El-Badry et al. 2023; Igoshev et al. 2023; Scholz 2024; Werner et al. 2024). Most of these fast objects also have luminosities similar to hot subdwarfs, but significantly lower temperatures than US 708. A possible evolutionary connection was proposed by Shen et al. (2018), according to which US 708 might represent the state after the puffed-up phase when the star is moving back to the WD cooling tracks (Bauer et al. 2019).

Hot subdwarfs are also good candidates for ejected companions in the single-degenerate SN Ia scenario, in which stable mass transfer from a donor star triggers the explosion of a WD (Whelan & Iben 1973; Nomoto 1982). The WD might either explode when reaching the Chandrasekhar mass or by the detonation of a helium layer on top of a C/O WD core (Woosley & Weaver 1994; Fink et al. 2010). The hot subdwarf companions ejected in such a scenario can reach velocities of several hundred km s−1 and exceed the escape velocity of the Galaxy (Justham et al. 2009; Wang & Han 2009; Geier et al. 2013; Neunteufel 2020; Neunteufel et al. 2019, 2021, 2022). Most of the predicted properties of those ejected He-star companions are consistent with the properties of US 708 (Geier et al. 2015a).

Several sdO/B binary progenitor candidates for this scenario with massive WD companions are known (Mereghetti et al. 2009, 2021; Vennes et al. 2012; Geier et al. 2013; Kupfer et al. 2022) as well as several other extremely close sdO/B+WD systems with lower mass companions (Kupfer et al. 2017a,b; Luo et al. 2024). Some of them are even found in a state of ongoing mass transfer (Kupfer et al. 2020a,b). Although those binaries do not fulfill all the criteria for SN Ia progenitor candidates, they provide evidence for a rich population of interacting sdB+WD systems. In contrast to that, no additional candidate for an ejected hot subdwarf companion similar to US 708 has been discovered yet.

The search for both the close binary progenitors and the surviving ejected companions of thermonuclear supernovae was one of the key science goals of a long-term survey project, which started in 2009 and is briefly summarised in the following sections. In this paper, we present the final results of the spectroscopic and kinematic survey of fast hot subdwarfs.

2. Previous results from the spectroscopic survey

A dedicated survey was performed to find hot subdwarf stars with massive compact companions, such as massive white dwarfs (> 1.0 M), neutron stars, or stellar-mass black holes: the “Massive Unseen Companions to Hot Faint Underluminous Stars from SDSS” project (MUCHFUSS). Hot subdwarf stars were selected from the Sloan Digital Sky Survey (SDSS, e.g. Ahn et al. 2012) by colour, followed by a visual inspection of the spectra. Hot subdwarf stars with high radial-velocity variations were selected as candidates for follow-up spectroscopy to derive the radial-velocity curves and the binary mass functions of the systems (see Geier et al. 2011a, 2015b).

Hot subdwarfs with extreme space velocities were studied in the Hyper-MUCHFUSS spin-off project. The initial survey was restricted to sdO/B stars with constant RVs exceeding 100 km s−1. Atmospheric parameters, spectroscopic distances and RVs were obtained from SDSS spectra. Proper motions were determined using position measurements from ground-based catalogues (Tillich et al. 2011). Twelve high-velocity sdO/Bs were discovered, either originating from the central bulge region of the Galaxy or the outer Galactic disc. While eleven of these stars were found to be bound to the Galaxy, one star was proposed to be unbound.

However, a more detailed analysis of this latter star PB 3877 revealed it to be a composite binary consisting of an sdB and a K-type main-sequence companion with an orbital period of the order of hundreds of days. Furthermore, a bound orbit could no longer be excluded. Since such a wide binary is quite fragile, any acceleration scenario requiring close binary interactions or close encounters with massive objects could be excluded, leaving a primordial or accreted halo origin as the only viable scenario for this star (Nemeth et al. 2016).

In contrast to that, a detailed follow-up study of US 708 (Geier et al. 2015a) based on improved spectroscopy and astrometry revealed that this star is moving even faster than estimated before (Hirsch et al. 2005). With a Galactic restframe velocity of ∼1000 km s−1, it turned out to be the fastest unbound star in the Galaxy known at the time (see also Neunteufel 2020). In addition, an origin close to the Galactic centre could be excluded, as well as the Hills ejection mechanism. Furthermore, the spectrum of US 708 shows significant rotational broadening of the spectral lines that is quite untypical for an He-sdO, but perfectly consistent with a prior phase of tidal spin-up in an ultracompact binary consisting of a low-mass hot subdwarf with ∼0.3 M and a WD of more than 1.0 M. These results are consistent with US 708 being the ejected companion of a thermonuclear SN.

Finally, detailed analyses of four more fast hot subdwarf candidates were performed: two from the sample of Tillich et al. (2011) and two more discovered in the MUCHFUSS survey (Geier et al. 2015b). All of these stars are bound to the Galaxy and three can be traced back to the Galactic centre or the inner disc. The He-sdO SDSS J205030.39−061957.8 turned out to be very similar to US 708 regarding its spectroscopic parameters. It originates from the outskirts of the Galactic disc, but it is also on a bound Galactic orbit (Ziegerer et al. 2017).

3. Revisiting and extending the sample

The results of the Gaia mission are a game changer, because they provide astrometry of unprecedented precision allowing for detailed kinematic studies. Large surveys have provided UV, optical and IR photometry, which allow us to construct detailed spectral energy distributions. Spectral databases (e.g. SDSS, LAMOST) provide multi-epoch spectroscopy. This wealth of new information motivated us to revisit and extend the sample of fast hot subdwarfs. We present the extension of the sample in Sect. 3. The quantitative spectral analysis is presented in Sect. 4. Multi-epoch spectra are used to check for radial velocity variability in Sect. 4.4. Spectral energy distributions are constructed and spectroscopic distances are derived in Sect. 5, and we then combine them with Gaia proper motions (PMs) to carry out a kinematic study (Sect. 6).

3.1. Target selection

To substantially increase the sample of high-velocity sdO/Bs and to find unbound objects similar to the prototype US 708, we performed an extended and more systematic survey, including archival spectroscopic and ground-based astrometric data. Our aim was to select candidates with high Galactic-restframe velocities, for which we wanted to obtain spectroscopic follow-up observations for detailed analyses. This part of the project was mostly conducted between 2015 and 2018 with the data and methods available to us at that time.

The starting point was again the regular data releases of the Sloan Digital Sky Survey (up to DR14), including newly classified hot subdwarfs (Geier et al. 2015b; Kepler et al. 2016, 2019). To further increase the input sample, we compiled a comprehensive catalogue of more than 5600 spectroscopically identified hot subdwarfs from the literature (Geier et al. 2017).

For more than 1900 sdO/Bs the spectroscopically determined fundamental parameters effective temperature Teff, surface gravity log g, helium abundance log n(He)/n(H) and radial velocity (RV) have been published in the literature. To improve our target selection, we performed quantitative spectral analyses of both the SDSS and the follow-up spectra of the most promising candidates for fast sdO/Bs. The method is described in Geier et al. (2011b). Spectroscopic distances have been determined as described in Ramspeck et al. (2001) assuming the canonical mass of 0.47 M for the subdwarfs.

Because our survey started well before the Gaia era, ground-based and hybrid ground- and space-based PMs have been taken from several catalogues (see Fig. 1). Since Galactic field sdO/B stars have typical distances of the order of several hundred pc to a few kpc, their PMs are quite small and often close to the detection limits of the ground-based surveys. Extreme care has to be taken using such PMs (see Ziegerer et al. 2015), which is why we applied a very conservative approach and only selected stars where all the independent PM measurements from the literature (each star has about 4 − 5) were significant and consistent with each other. Only in a few cases, where the majority of the other measurements were consistent and the data quality high, we allowed for one outlier at most.

thumbnail Fig. 1.

Proper motion measurements in right ascension (left panel) and declination (right panel). Gaia proper motions from Early Data Release 3 (EDR3, Gaia Collaboration 2021) are plotted against ground-based measurements from SDSS DR9 (Ahn et al. 2012, open circles), the PPMXL catalogue of positions and proper motions on the ICRS (Roeser et al. 2010, open squares), the Whole-Sky USNO-B1.0 Catalogue (Monet et al. 2003, open diamonds), the Absolute Proper motions Outside the Plane catalogue (APOP, Qi et al. 2015, open downward triangles), Gaia-PS1-SDSS (GPS1, Tian et al. 2017, open upward triangles), and Hot Stuff for One Year (HSOY, Altmann et al. 2017, crosses). PMs from the Fourth U.S. Naval Observatory CCD Astrograph Catalogue (UCAC4, Zacharias et al. 2013) were also used, but they are not plotted here, because only a few stars have measurements in this catalogue. The averages of the ground-based values from those surveys are plotted as filled diamonds. The ground-based measurements by Tillich et al. (2011) and Ziegerer et al. (2017) are plotted as open dark grey diamonds.

We calculated the Galactic-velocity components and the total Galactic-restframe velocities by combining the RV, PM and spectroscopic distance measurements. For the spectroscopic follow-up we finally selected the candidates with Galactic-restframe velocities vgrf > 300 km s−1.

3.2. Spectroscopic follow-up observations

A follow-up campaign was conducted to obtain higher-quality spectra and to improve the accuracy of the derived spectroscopic parameters. Second-epoch spectra are also needed to identify RV variable stars, since the close binary fraction of sdBs is very high. Reliable kinematic parameters can only be derived for stars with constant RV, and only single hot subdwarfs can be candidates for the ejected companions of thermonuclear supernovae.

The campaign was conducted from 2015 to 2018 in dedicated observing runs with 4m- to 10m telescopes equipped with medium-resolution spectrographs. In total more than 26 nights of observing time were granted for this project (see Table 1).

Table 1.

Summary of follow-up observations during the course of the MUCHFUSS project.

Additional spectra were retrieved from the SDSS DR17 (taken with the original SDSS spectrograph and the BOSS spectrograph, Abdurro’uf et al. 2022) and LAMOST DR8 databases. An overview of the spectroscopic data used for the spectroscopic analysis is given in Table A.1.

3.3. Cleaning the sample

We removed many false positives from the initial sample of fast sdO/B candidates as observational data were successively added. This happened for different reasons outlined in the following.

Stars with erroneous parameters. Over the last decade, several new ground-based proper-motion catalogues were published (see Sect. 3.1) and continually used to check the consistency of the individual measurements. Many stars turned out to have inconsistent PMs and were therefore removed from the sample. Finally, the highest quality, space-based PMs provided by Gaia allowed us to eliminate the remaining candidates with erroneous measurements.

Figure 1 shows a comparison of the ground-based proper motions of our final sample from the most relevant surveys with the proper motions from the Gaia mission. The results are, except for some outliers, generally consistent within the uncertainties confirming our ground-based target selection.

The spectroscopic follow-up allowed for an improved determination of the atmospheric parameters with respect to the ones determined from the SDSS spectra alone. In many cases this had a significant impact on the spectroscopic distances, which had previously often been overestimated. Many candidates turned out to be much closer than estimated and their Galactic velocity components are thus much slower. One of those stars, the He-sdO SDSS J141812.51−024426.8 was recently reanalysed by Werner et al. (2022).

Composites and blends. Higher resolution and quality spectra revealed spectroscopic signatures of cool MS companions in several of our programme stars, which initially ranked among our top candidates for fast sdO/Bs. The high space velocities of those stars are likely caused by inaccurate atmospheric parameters, which were determined by fitting single-star model spectra and led to an overestimation of the spectroscopic distances. However, we also found a composite system with a confirmed high space velocity in our survey (Nemeth et al. 2016).

To identify other composite binaries or blended sources we performed an analysis of the spectral energy distributions (SEDs) of all the stars from our sample, using a χ2-minimisation method to match filter-averaged magnitudes calculated from model spectra to observed magnitudes drawn automatically from photometric catalogues in the literature (see Heber et al. 2018 for a more detailed description and Table A.2 for a list of the photometric catalogues used). The method allows us to fit two stellar components and can be used to determine the effective temperatures of both stars, the surface ratio and the angular diameter of the primary in composite binaries (e.g. Dorsch et al. 2021). Since the SEDs include photometry up to the infrared regime, this method is very well suited to detecting cool companions, even if they are too faint to show significant spectroscopic features. In addition, comparing the surface ratio with theoretical predictions for the radii of hot subdwarfs and cool companions of the determined temperatures, significant mismatches allow us to identify physically unassociated blends.

In this paper, we only present the analysis of the sdO/Bs with SEDs indicative of single stars, which do not show any additional IR excess. The analysis of the composite binaries will be presented in another paper of this series (Heber et al., in prep.).

Proto-He-WDs. Another potential reason for misclassification as a very fast star might be the underlying assumption that the stars are in the (post-)EHB phase and that they all have the same mass of 0.47 M. It is known that red-giant-branch stars can be stripped via binary interactions even before they ignite core helium burning (e.g. Iben & Tutukov 1986; Benvenuto & De Vito 2004). Those objects subsequently cool down and evolve to become He-WDs. Stars in this proto-He-WD stage can have similar atmospheric parameters to sdB or blue horizontal branch (BHB) stars, but they have significantly lower masses (∼0.2 − 0.3 M), often smaller radii, and therefore lower luminosities (e.g. Heber et al. 2003; Kawka et al. 2015).

That such objects can be easily misclassified as much more luminous and therefore distant stars is impressively demonstrated by the putative top-hypervelocity BHB candidate SDSS J160429.12+100002.2. Detailed spectroscopic follow-up revealed that it is a proto-He-WD in a close binary, which is orbited by a very likely substellar companion (Irrgang et al. 2021), drastically reducing its estimated space velocity.

Although these proto-He-WDs appear to be much rarer than sdO/Bs in the (post-)EHB phases (especially at the hot end of the EHB, where their evolutionary times are one to two orders of magnitude shorter), our selection for apparently large space velocities may be biased towards such low-mass objects. We therefore conclude that the relative fraction of proto-He-WDs in our sample may be higher than in the field population.

Runaway main-sequence star. Finally, the B-type star PG 1610+062, an object originally suspected to be an RV-variable hot subdwarf rather than a fast moving one (Geier et al. 2015b), eventually turned out to be a runaway main sequence star with an exceptionally high ejection velocity from the disc, which challenges current ejection scenarios (Irrgang et al. 2019, 2021).

4. Homogeneous spectral analysis

The final selected objects showing diverse spectral types are classified, analysed and checked for RV variability in the following sections.

4.1. Spectral classification

A spectral classification scheme was established by Drilling et al. (2013) and was applied, for example by Lei et al. (2019) and Jeffery et al. (2021). This classification scheme adds a third dimension to the MK scheme by defining 40 helium classes to address the diverse helium line strengths in the optical spectra. We refrained from applying this scheme to our sample, but we condensed it into just nine spectral sub-classes. We can easily distinguish hydrogen-dominated from helium-dominated spectra. Helium lines may or may not be present for the H-dominated spectral types. In sdB stars He I lines may be present but no He II lines. Somewhat hotter H-strong-lined stars are classified as sdOB if they display He II 4686 Å in addition to He I lines. In sdO stars, the Balmer lines dominate the spectrum, but He II lines are stronger than the He I lines; the latter may even be absent from the hottest ones.

The helium-strong spectral types are more delicate to classify because both temperature and helium abundance have a strong impact on the spectral appearance. We may distinguish extremely He-strong-lined hot subdwarfs (He-sdB, He-sdOB and He-sdO) from intermediate ones by the presence or absence of any hydrogen Balmer lines. In the former case, they are classified as intermediate (iHe), in the latter as extreme helium (eHe) subdwarfs. The relative strength of He I and He II lines varies with temperature. While no He II lines are detectable in He-sdBs, He II lines are stronger than He I lines in He-sdOs. In He-sdOBs the strongest lines are from He I, but He II is detectable and may reach similar strength.

Hence, we distinguish three H-strong-lined classes (sdB, sdOB, and sdO) and six He-strong-lined classes (iHe-sdB, iHe-sdOB, iHe-sdO, eHe-sdB, eHe-sdOB, and eHe-sdO). In Table A.3 we give the spectral classifications of all programme stars1. In Fig. 2 we show examples for H-strong-lined types, while Fig. 3 shows examples of He-dominated spectra.

thumbnail Fig. 2.

Excerpts of BOSS spectral fits of He-poor hot subdwarfs. From top to bottom: SDSSJ124248.89+133632.6 (sdB), SDSSJ124310.58+343358.4 (sdOB), SDSSJ124819.08+035003.2 (sdO). In the upper panels the best-fit synthetic spectrum (red) is shown along with the observed one. The lower panels give residual χ, whereby the long-dashed horizontal lines mark the zero values and the short-dashed horizontal lines mark the deviations in terms of 1σ. Parts of the spectra, which have been excluded from the fits, are marked in grey.

thumbnail Fig. 3.

Same as Fig. 2, but for He-rich hot subdwarfs. From top to bottom: SDSSJ120352.24+235343.3 (iHe-sdOB), SDSSJ222515.34-011156.8 (eHe-sdOB), SDSSJ080833.76+180221.8 (eHe-sdO).

The final sample of fast hot subdwarfs comprises 53 stars. We classified 29 sdBs (55%), five sdOBs (9%), five sdOs (9%), seven iHe-sdOBs (13%), two iHe-sdOs (4%), two eHe-sdOBs (4%), two eHe-sdO (4%), and one BHB (2%). There are no He-sdB stars in the sample.

4.2. Model atmospheres

The chemical composition of sdB stars is known to be non-solar and varies considerably from star to star (Pereira 2011; Naslim et al. 2013; Geier 2013). Because individual abundances for the programme stars are not available, we used the average abundance pattern from Pereira (2011) as reported by Naslim et al. (2013) for the model atmosphere calculations.

Model atmospheres and synthetic spectra were calculated in a hybrid local thermodynamic equilibrium (LTE)/non-LTE (NLTE) approach with a suite of numerical codes. The ATLAS12 code (Kurucz 1996) was used to calculate the temperature and density stratification of the atmosphere in LTE, with the most recent Kurucz line lists2 used to incorporate the line-blanketing effect. The opacity sampling technique implemented in ATLAS12 allows non-solar chemical abundances to be used without pre-computing opacity distribution functions.

Once the model atmosphere calculations converged, detailed hydrogen and helium model atoms were set up to calculate occupation numbers for all relevant levels of hydrogen, He I, and He II. The DETAIL code (Giddings 1981; Butler & Giddings 1985) numerically solves the coupled equations for radiative transfer and statistical equilibrium. The emerging H/He spectrum was then computed with the SURFACE code (Butler & Giddings 1985) using detailed line-broadening tables for hydrogen and helium transitions (see Irrgang et al. 2018a for more details).

We calculated extensive grids covering the entire parameter range of hot subluminous stars. The backbone grid covers effective temperatures from 9000 K to 55 000 K, surface gravities log g = 4.6 to 6.6 and helium abundances from log n(He) = −5.05 to −0.041. Extensions to high Teff (75 000 K for log g = 5.2 to 6.6) and high helium abundances up to log n(He) = −0.001 (for Teff = 25 000 to 55 000 K) were computed to cover the parameters for all sub-types.

4.3. Quantitative spectral analyses

We replaced the previously used selective χ2 minimisation fitting procedure (Napiwotzki 1999; Hirsch 2009) by a global one developed by Irrgang et al. (2014). This means that we considered the whole useful spectral range (usually 3600 Å to 6700 Å) rather than pre-selecting spectral lines, removing artifacts, and simultaneously fitting the entire multi-parameter space (effective temperature, surface gravity, helium abundance, projected rotational velocity, and radial velocity). Available spectra of a star were simultaneously fit to determine the atmospheric parameters, projected rotation velocity and radial velocities. Example fits of BOSS spectra are shown (blue spectral parts only) in Figs. 2 and 3. Fits of WHT/ISIS spectra are shown in Fig. 4. Reduced χ is also plotted to demonstrate the quality of the fit. The resulting atmospheric parameters are listed in Table A.3.

thumbnail Fig. 4.

Same as Fig. 2, but showing fits of the WHT/ISIS spectra of hot subdwarfs. From top to bottom: SDSSJ130543.97+115840.8 (sdB), SDSSJ221728.35+121642.6 (sdOB), SDSSJ120352.24+235343.3 (iHe-sdOB).

The Teff − log g and Teff − log n(He)/n(H) diagrams (Figs. 5 and 6) generally resemble the corresponding diagrams of the field population (reviewed by Heber 2016). Furthermore, it can be seen that the overall distribution of parameters determined in our new homogeneous analysis matches the predicted location close to the EHB better than the preliminary analysis on which the target selection was based.

thumbnail Fig. 5.

Teff − log g diagram of sample of fast stars. The colour scales with the helium abundance from light orange to dark purple. The light grey symbols mark the results of the preliminary spectroscopic analysis on which the target selection was based. They are connected with the corresponding values from the new homogeneous analysis. The zero-age and terminal-age EHBs for a subsolar metallicity of −1.48 (solid grey lines) were interpolated from evolutionary tracks by Dorman et al. (1993). The helium main sequence (dashed line) is taken from Paczynski (1971).

thumbnail Fig. 6.

Teff − log n(He)/n(H) diagram of sample compared to the preliminary analysis (see Fig. 5). Solar helium abundance is marked by the solid horizontal line, while the dotted line marks the transition between intermediate and extreme helium abundance.

4.4. Multi-epoch radial velocities

Radial velocities of all the stars have been measured along with the atmospheric parameters from the available medium-resolution spectra. The low-resolution spectra obtained with GTC/OSIRIS are not suited to measuring sufficiently accurate RVs. We also encountered problems with spectra taken during our observing run with WHT/ISIS in July, 2017 (see Table 1), where the seeing reached down to the astonishingly small value of 0.3 arcsec. Since the smallest slit provided by ISIS is of 0.6 arcsec and the telescope was not defocused to match the slit width, significant movement of the stars within the slit lead to erroneous RVs, which had to be discarded.

Because a large fraction of the sdB stars is known to be radial-velocity variable, multi-epoch measurements are needed to identify the variables in the sample. Therefore, we made use of the individual observations and the combined SDSS spectra from the final data release of the fourth phase of the Sloan Digital Sky Survey (SDSS-IV, DR17, Abdurro’uf et al. 2022), as well as those from our follow-up spectroscopy.

To estimate the probability that a star is variable in radial velocity, we calculated the statistical p value as described by Maxted et al. (2001). Assuming χ2 statistics, we determined the goodness of fit of a constant radial velocity to the observed radial velocities. Thus the p value is the probability of obtaining the observed χ2 value or higher from random fluctuations of a constant value. To allow for the systematic errors, we added an external error in quadrature to all radial velocity uncertainties prior to calculating statistics. For spectra taken with the SDSS and BOSS spectrographs Yanny et al. (2009) and Bolton et al. (2012) estimated systematic radial-velocity uncertainties of 4 km s−1 and 4.5 km s−1, respectively. Here we adopted 5 km s−1 for all spectra, irrespective of their spectral resolution. If log p < −4, we consider the variations not to be caused by random fluctuation, but the star to be a radial-velocity-variable close binary (Maxted et al. 2001).

For all stars, multi-epoch observations are available, and the number of individual measurements and the resulting p values are given in Table A.3. The timespans between the individual RVs are quite different, and we therefore, divided the sample into four categories. Objects in the best category, which is number 1, have several spectra from SDSS, which were taken at distinct epochs separated by one day or more. In addition, category 1 stars have follow-up observations taken with different instruments usually several years after the SDSS epochs. Stars of category number 2 have SDSS spectra only, but these were again taken at several distinct epochs separated by at least one day. Stars of category number 3 only have individual SDSS spectra, which were taken back-to-back usually within less than one hour, but their longer term RV variability can be checked with follow-up spectra taken with different instruments. Finally. the least constraining category, number 4, only has individual SDSS spectra taken back-to-back within less than one hour of each other. Since the typical orbital periods of sdBs range from hours to days, category 4 stars can only probe the closest binary systems with the shortest periods. Non-detections of RV variations in category 4 stars are therefore not suited to excluding close binarity in general.

The p values listed in Table A.3 are high for most stars with no indication of RV variability (see also Fig. 7). Only two stars in the sample show evidence of significant RV variations. These are the sdOB J180313.45+234000.1 and the sdB J102439.43+383917.9, both considered members of category 1. The former has log p = −4.43, which is marginally above our detection threshold. The four SDSS spectra, however, show a higher log p indicative of constant RV and the deviation stems only from one follow-up spectrum taken with FORS1. We therefore treated this star as a candidate for RV variability only. The latter, on the other hand, shows both a lower log p = −5.87 and a more consistent variation among the nine individual spectra taken with SDSS and ISIS. We therefore classify J102439.43+383917.9 as a close binary.

thumbnail Fig. 7.

log p distribution of sample. The significance level log p < −4.0 is marked by the vertical line.

5. Spectral energy distribution and spectroscopic distances

Photometric measurements are obtained by querying several public databases. The optical spectral range is well covered by observations in the Sloan filters and broad band Gaia measurements. Most stars have UV observations from the GALEX mission, while infrared magnitudes are only available for a few stars (e.g. Table A.2 for a list of the photometric surveys used).

We fitted the observed magnitudes with synthetic SEDs computed with ATLAS12. The stellar parameters were set to the values determined from spectroscopy. The angular diameter Θ, and the interstellar colour excess E(44 − 55) remained as free parameters to match the observed SED (see Heber et al. 2018 for more details). To model interstellar extinction, the function given by Fitzpatrick et al. (2019) was used, assuming a standard extinction coefficient of R(55) = 3.02. For illustration we show the fits for nine stars of different spectral types in Fig. 8.

thumbnail Fig. 8.

Fits of SEDs of selected programme stars. Each plot consists of two panels; the upper one compares the observed fluxes to the synthetic SED. To ease the slope of the distribution the flux is multiplied by the wavelength to the power of three. Photometric fluxes are displayed as coloured data points (for details see Table A.2) with their respective uncertainties and filter widths (dashed lines). The best-fit models are drawn as full grey lines. The lower panels give reduced χ to demonstrate the quality of the fit. Example fits are selected to demonstrate relevant issues. In each row, the stars are arranged with increasing temperature. Top row: He-weak sdB stars: cool sdB J123953.52+062853.0, which is the most distant programme star; J135651.26+155810.4, the sdB with the largest (positive) heliocentric RV; and J123137.56+074621.7. Middle row: sdOB stars: J164419.44+452326.7; the low gravity J130543.97+115840.8; and the highly reddened J180313.45+234000.1. Bottom row: He-sdO stars: iHe-sdOB J204358.55–065025.8; the best observed iHe-sdOB star J022422.21+000313.5; and the eHe-sdO J205030.39–061957.8, which was considered a spectroscopic twin to US 708 by Ziegerer et al. (2017). The strong colour excess of J180313.45+234000.1 and J205030.39–061957.8 is most obvious from the strength of the 2200 Å UV feature effecting the GALEX NUV band.

Adopting the theoretically predicted degenerate core mass (0.47 M) for the helium flash in low-mass stars, with an uncertainty of ±0.05 M to allow for different metallicities, the radii R of our sample stars are calculated using the surface gravity and R = G M / g $ R=\sqrt{GM/g} $. The spectroscopic distances are determined using the angular diameter d = 2R/Θ (see Table A.5).

6. Astrometry and kinematic analysis

Due to the large distances of the stars in our sample, most Gaia parallaxes are of insufficient accuracy to derive meaningful distances from astrometry only (see Table A.4). Instead we calculated the spectroscopic distances as described in Sect. 5 using the improved spectroscopic parameters derived from the follow-up analysis. Proper motions and astrometric correlation coefficients were taken from Gaia Early Data Release 3 (EDR3; GaiaGaia Collaboration 2021; see Table A.4).

The Galactic velocity components (U points away from the Galactic centre, V points in the direction of Galactic rotation, and W points toward Galactic north) and the total Galactic restframe velocities have been updated accordingly.

Four stars of our sample have Gaia parallaxes with uncertainties smaller than 20%, which is accurate enough to compare them directly to the spectroscopic distances (see Table 2). The parallax distances are consistent with the spectroscopic distances within the uncertainties for all those stars confirming the validity of our approach (see Table 2).

Table 2.

Comparison of spectroscopic versus parallax ϖ based distances. The inversion accounts for asymmetry and zero-point offset (Lindegren et al. 2020).

The Galactic orbits were calculated as described in Irrgang et al. (2018b) using Model I of Irrgang et al. (2013) for the Galactic potential. The orbits were traced back and the eccentricities of the Galactic orbits e as well as the components of the angular momentum in Z-direction LZ determined. The uncertainties in the input parameters were propagated with Monte Carlo simulations assuming Gaussian distributions while also accounting for asymmetric errors and the correlations between the proper motion components.

The kinematic parameters of our sample are provided in Table A.6. Besides the Galactic velocity components U, V and W as well as the total Galactic rest frame velocity vgrf corrected for the motion of the local standard of rest, we provide e and LZ and the time of flight td since the last crossing of the Galactic disc within 100 pc of Z = 0.

Figure 9 shows the distribution of distances from the Galactic plane while the Toomre and the LZ − e diagrams are shown in Fig. 10. Most stars are situated well above the Galactic disc in the halo, but all of them are bound to the Galaxy. The prototype US 708 remains the only unbound HVS sdO/B. The sample represents a kinematically extreme population in the Galactic halo and contains only very few stars, which could also belong to the disc population.

thumbnail Fig. 9.

Distribution of distances of our sample from the Galactic plane. The dashed vertical lines mark the scale height of the Galactic thin disc; the solid vertical lines mark the scale height of the thick disc.

thumbnail Fig. 10.

Left panel: Toomre diagram of sample. The 2σ contours of the thin disc (grey), the thick disc (green), and the halo (blue) are constructed from the U, V, and W distributions of Anguiano et al. (2020). Objects with halo kinematics are marked in blue; the few objects also consistent with thick-disc kinematics are marked in green. Right panel: e − LZ diagram of the sample using the same colour coding for the stars. The green and grey 2σ contours in the upper left denote the locations of thin- and thick-disc stars, respectively, based on the GALAH + Gaia EDR3 sample of Buder et al. (2021). The dashed lines mark prograde (positive LZ) and retrograde (negative LZ) directions of Galactic rotation, respectively.

7. The sample of fast hot subdwarf stars: Ejected or born in the halo

After a thorough analysis of our sample of fast hot subdwarfs, we can now evaluate whether the properties of the sample are consistent with a genuine halo population or whether there is evidence of a (sub-)population of stars ejected from the disc.

Although bound to the Galaxy, the extreme kinematics of our sample might still hint at the SN Ia scenario or other ejection scenarios from the disc. Since the core-collapse and the thermonuclear supernova scenarios require either the explosion of a massive star or the prior evolution of an intermediate-mass close binary (e.g. Geier et al. 2013), runaways from these scenarios must originate from young stellar populations, likely in the Galactic disc. Since the lifetime of sdO/B stars is of the order of 100 Myr at most, a maximum possible time of flight since the ejection is defined. We calculated the flight times from the last disc crossing for our sample, the location of the crossings in Galactic coordinates, and the velocity at the moment of the crossing (or “ejection velocity”, vej). Out of 53 stars, 36 (68%) could have been ejected from the disc within the uncertainties (see Table A.7), assuming a disc radius of 15 kpc. This means that the SN ejection scenario is likely excluded for the remaining 32%.

As shown in Fig. 11, most Galactic plane crossing locations are at Rdisc < 15 kpc; only two stars cross further out. The disc crossings are quite concentrated towards the central region of the Galaxy and fall off sharply at larger galactocentric distances. There is also no difference in the crossing locations between stars, which fulfill the criteria for an ejection and the rest of the sample. This is expected in a halo population, because the density of stars in the halo is known to be highest in the central bulge region and to drop off in the radial direction (e.g. Bland-Hawthorn & Gerhard 2016 and references therein). In addition, we know that the 32% stars of our sample are very likely halo stars.

thumbnail Fig. 11.

Upper panel: Points of last Galactic disc crossings of stars in our sample in the XY plane. The two stars with the largest distances have been omitted here for better visualisation. Lower panel: Distribution of the galactocentric radii at disc crossing for all stars of the sample. Stars that could have been ejected from the disc are marked in dark grey, and the others are in light grey.

The distribution of disc-crossing velocities is shown in Fig. 12. Unlike the (hyper-)runaway A/B-type stars studied by Irrgang et al. (2021), the present sample does not contain stars with high ejection velocities at large Rdisc. On the contrary, the ejection velocities are in fact highest towards the Galactic centre and decrease with galactocentric distance. This is in part a bias introduced by our selection of stars by their high current space velocity. Since the gravitational attraction of the Galaxy decreases with galactocentric distance, the stars coming from the central part of the Galaxy require higher disc-crossing velocities than those from further out in the disc to remain fast at their present positions.

thumbnail Fig. 12.

Ejection or disc-crossing velocities plotted against galactocentric distance.

The halo scenario is further corroborated by the diversity of spectral types in our sample. Although we did not derive spectroscopic distances of the composite binaries, we identified their locations at high Galactic latitudes in our survey and their faintness point to a halo membership (except the well-studied system PB 3877, Nemeth et al. 2016). Our extended sample, therefore, consists of almost all known types of hot subdwarfs including composite binaries, which is another indication that we have a halo population rather than a sample of ejected stars, which would have to be less diverse and should consist exclusively of single stars. In addition, other members of this halo population with similarly extreme kinematics to our fast sdO/Bs, but in earlier evolutionary stages on the main sequence and the red giant branch have been discovered (Scholz et al. 2015; Hattori et al. 2018; Li et al. 2021).

On the other hand, the distribution of the stars in the Teff − log g and the Teff − log n(He)/n(H) diagrams (see Figs. 5 and 6) resembles the typical field population of sdO/Bs (e.g. Heber 2016; Luo et al. 2021; Geier et al. 2022) rather than the markedly different distribution seen in the old halo globular cluster ω Cen (Latour et al. 2018, 2023). It is, however, not clear whether ω Cen is a good proxy for the halo population of sdO/Bs as such, because other old GCs (e.g. NGC 6752) contain a different mixture of hot subdwarf types (Latour et al. 2023).

The extremely small fraction of RV-variable close binaries in our sample is especially for the hydrogen-rich sub-types completely different than what is observed in the disc population (see Geier et al. 2022 and references therein). We only found one star in our sample to be significantly variable (J102439.43+383917.9), next to a single candidate for RV variability (J180313.45+234000.1). Both stars are kinematically more likely to be members of the thick disc rather than the halo population – we do not find any RV variable close binaries in our sample of hot subdwarf stars with halo kinematics. This is consistent with the lack of such binaries in ω Cen (Latour et al. 2018). A possible explanation might be the age dependency of the contribution of the different binary formation channels. For older populations, the relative contribution of mergers is predicted to rise, leading to a higher fraction of single stars (Han 2008; Lisker & Han 2008). This strongly indicates that the lack of hot subdwarfs in close binary systems in GCs is an age effect rather than an environmental effect.

We therefore conclude that we neither found another unbound HVS sdO/B similar to the prototype US 708 nor strong candidates for bound ejected companions of SN Ia in our spectroscopic survey. Of the 53 stars in our sample, 48 are very likely halo sdO/Bs with extreme kinematics and five are more likely members of the thick-disc population.

8. Discussion and outlook

The lack of other good sdO/B-type candidates for the single-degenerate SN Ia ejection scenario and the unclear relation to other fast (pre-)WDs (e.g. El-Badry et al. 2023) makes the prototype US 708 even more special.

However, it might still be premature to draw strong conclusions on the birth rate of such objects and the number of corresponding SN Ia or other sub-types of stellar explosions. The observed lack of other candidates might be related to our search strategy. Neunteufel et al. (2021, 2022) simulated a model population of He-stars ejected by SNe Ia and determined the observable properties of this sample, such as their positions and velocity components. The authors concluded that most of the ejected stars observed should be located in and around the Galactic bulge, where the stellar density is the highest. In addition, their simulation predicted that many objects with high absolute space velocities would have rather small proper motions, and they would thus only be detectable through their large RVs; this requires spectral observations. Furthermore, most of those objects should orbit on bound trajectories with velocities around 500 km s−1 and should have masses around 0.8 M – more massive than typical sdO/Bs. In these simulations, US 708 represents an extreme case of a low-mass He-star ejected from the tightest possible binary progenitor, which was already suggested by Geier et al. (2015a).

Comparing these predictions with our selection criteria reveals some potentially quite significant biases. By selecting objects from the northern SDSS survey, we missed the Galactic-bulge region, which is best accessible from the south. Focusing on objects with high tangential velocities and therefore relatively high proper motions instead of RVs, we did not impose a particularly strong selection criterion to single out the ejected sdO/B population. Probably most importantly, we selected sdO/Bs of all types, whereas the predicted higher mass He-stars should be He-sdOs with quite high luminosities. Since those stars have much shorter lifetimes on the EHB (∼107 yr) than the bulk of the low-mass sdO/Bs (∼108 yr, e.g. Han et al. 2002), they are rarer. If ejected from the disc at moderate velocities, they will not fly very deep into the halo and would be unable to reach our search area at high Galactic latitudes before turning into WDs. We therefore suspect that our selection criteria were likely not specific enough to find the best candidates for ejected companions in the thermonuclear SN scenario.

We expect the situation to change for the better in the years to come. Using Gaia astrometry and photometry, we are now able to compile extensive, all-sky catalogues of hot subluminous stars (e.g. Culpan et al. 2022), which will be used as input target lists for the next generation of spectroscopic surveys, in particular those covering the yet neglected southern sky (4MOST, de Jong et al. 2019; SDSS V, Kollmeier et al. 2017). In addition, we will use simulations such as those of Neunteufel et al. (2021, 2022) to define better selection criteria for targeted spectroscopic follow-up campaigns. Based on those spectra, which will reveal both the nature and the RV of tens of thousands of hot subluminous stars, we will be able to continue our search for the ejected companions of SN Ia. We are thus confident that we will find the hidden siblings of US 708.


1

J174211.75+643009.9 turns out to be a BHB star.

Acknowledgments

S.G. was supported by the Deutsche Forschungsgemeinschaft (DFG) through a Heisenberg fellowship (GE2506/8-1). V.S. was supported by the Deutsche Forschungsgemeinschaft (DFG) through grant GE2506/9-1. I.P. acknowledges support by the Deutsche Forschungsgemeinschaft (DFG) through grant GE2506/12-1 and from a Royal Society University Research Fellowship (URF\R1\231496). M.D. is supported by the German Aerospace Center (DLR) through grant 50OR2304. P.N. acknowledges support from the Grant Agency of the Czech Republic (GAČR 18-20083S). N.R. is supported by the Deutsche Forschungsgemeinschaft (DFG) through grant RE3915/2-1. R.R. acknowledges support from Grant RYC2021-030837-I funded by MCIN/AEI/10.13039/501100011033 and by “European Union NextGeneration EU/PRTR”. This work was partially supported by the AGAUR/Generalitat de Catalunya grant SGR-386/2021 and by the Spanish MINECO grant PID2020-117252GB-I00. We thank M. Pritzkuleit for his comments and suggestions. Based on observations at the La Silla-Paranal Observatory of the European Southern Observatory for programmes number 095.D-0055(A), 095.D-0588(A), 297.D-5004, 099.D-0208(A), 0102.D-0092(A). Based on observations with the William Herschel Telescope operated by the Isaac Newton Group at the Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias on the island of La Palma, Spain. Based on observations collected at the Centro Astronómico Hispano Alemán (CAHA) at Calar Alto, operated jointly by the Max-Planck Institut für Astronomie and the Instituto de Astrofísica de Andalucía (CSIC). Based on observations made with the Gran Telescopio Canarias (GTC), installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias, in the island of La Palma. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Website is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington. Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS-IV acknowledges support and resources from the Center for High Performance Computing at the University of Utah. The SDSS website is www.sdss4.org. SDSS-IV is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, Center for Astrophysics | Harvard & Smithsonian, the Chilean Participation Group, the French Participation Group, Instituto de Astrofísica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe (IPMU)/University of Tokyo, the Korean Participation Group, Lawrence Berkeley National Laboratory, Leibniz Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astrophysik (MPA Garching), Max-Planck-Institut für Extraterrestrische Physik (MPE), National Astronomical Observatories of China, New Mexico State University, New York University, University of Notre Dame, Observatário Nacional/MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Autónoma de México, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University.

References

  1. Abbott, T. M. C., Adamov, M., Aguena, M., et al. 2021, ApJS, 255, 20 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35 [CrossRef] [Google Scholar]
  3. Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2012, ApJS, 203, 21 [Google Scholar]
  4. Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [Google Scholar]
  5. Altmann, M., Edelmann, H., & de Boer, K. S. 2004, A&A, 414, 181 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Altmann, M., Roeser, S., Demleitner, M., Bastian, U., & Schilbach, E. 2017, A&A, 600, A4 [Google Scholar]
  7. Anguiano, B., Majewski, S. R., Hayes, C. R., et al. 2020, AJ, 160, 43 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bauer, E. B., White, C. J., & Bildsten, L. 2019, ApJ, 887, 68 [Google Scholar]
  9. Benvenuto, O. G., & De Vito, M. A. 2004, MNRAS, 352, 249 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bianchi, L., Shiao, B., & Thilker, D. 2017, VizieR Online Data Catalog: II/335 [Google Scholar]
  11. Blaauw, A. 1961, Bull. Astron. Inst. Neth., 15, 265 [NASA ADS] [Google Scholar]
  12. Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529 [Google Scholar]
  13. Bobylev, V. V., & Bajkova, A. T. 2019, Astron. Lett., 45, 580 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bolton, A. S., Schlegel, D. J., Aubourg, E., et al. 2012, AJ, 144, 144 [NASA ADS] [CrossRef] [Google Scholar]
  15. Bravo, E., Gil-Pons, P., Gutierrez, J. L., & Doherty, C. L. 2016, A&A, 589, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Buder, S., Sharma, S., Kos, J., et al. 2021, MNRAS, 506, 150 [NASA ADS] [CrossRef] [Google Scholar]
  17. Butler, K., & Giddings, J. R. 1985, Newsletter of Analysis of Astronomical Spectra (London: Univ. London) [Google Scholar]
  18. Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2017, VizieR Online Data Catalog: II/349 [Google Scholar]
  19. Chen, X., Han, Z., Deca, J., & Podsiadlowski, Ph. 2013, MNRAS, 434, 18 [Google Scholar]
  20. Culpan, R., Geier, S., Reindl, N., et al. 2022, A&A, 662, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. de Boer, K. S., Aguilar Sanchez, Y., Altmann, M., et al. 1997, A&A, 327, 577 [NASA ADS] [Google Scholar]
  22. de Jong, R. S., Agertz, O., Agudo Berbel, A., et al. 2019, The Messenger, 175, 3 [NASA ADS] [Google Scholar]
  23. Dorman, B., Rood, R. T., & O’Connell, R. W. 1993, ApJ, 419, 596 [NASA ADS] [CrossRef] [Google Scholar]
  24. Dorsch, M., Jeffery, C. S., Irrgang, A., Woolf, V., & Heber, U. 2021, A&A, 653, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Dorsch, M., Reindl, N., Pelisoli, I., et al. 2022, A&A, 658, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Drilling, J. S., Jeffery, C. S., Heber, U., Moehler, S., & Napiwotzki, R. 2013, A&A, 551, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Drlica-Wagner, A., Carlin, J. L., Nidever, D. L., et al. 2021, ApJS, 256, 2 [NASA ADS] [CrossRef] [Google Scholar]
  28. El-Badry, K., Shen, K. J., Chandra, V., et al. 2023, IEEE Open J. Antennas Propag., 6, 28 [Google Scholar]
  29. Fink, M., Röpke, F. K., Hillebrandt, W., et al. 2010, A&A, 514, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Fitzpatrick, E. L., Massa, D., Gordon, K. D., Bohlin, R., & Clayton, G. C. 2019, ApJ, 886, 10 [NASA ADS] [CrossRef] [Google Scholar]
  31. Gänsicke, B. T., Koester, D., Raddi, R., Toloza, O., & Kepler, S. O. 2020, MNRAS, 496, 4079 [CrossRef] [Google Scholar]
  32. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Geier, S. 2013, A&A, 549, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Geier, S., Nesslinger, S., Heber, U., et al. 2007, A&A, 464, 299 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Geier, S., Hirsch, H., Tillich, A., et al. 2011a, A&A, 530, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Geier, S., Maxted, P. F. L., Napiwotzki, R., et al. 2011b, A&A, 526, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Geier, S., Marsh, T. R., Wang, B., et al. 2013, A&A, 554, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Geier, S., Fürst, F., Ziegerer, E., et al. 2015a, Science, 347, 1126 [Google Scholar]
  39. Geier, S., Kupfer, T., Heber, U., et al. 2015b, A&A, 577, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Geier, S., Østensen, R. H., Nemeth, P., et al. 2017, A&A, 600, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Geier, S., Dorsch, M., Pelisoli, I., et al. 2022, A&A, 661, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Giddings, J. R. 1981, Ph.D. Thesis, University of London, UK [Google Scholar]
  43. Greenstein, J. L., & Sargent, A. I. 1974, ApJS, 28, 157 [Google Scholar]
  44. Han, Z. 2008, A&A, 484, L31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Han, Z., Podsiadlowski, P., Maxted, P. F. L., Marsh, T. R., & Ivanova, N. 2002, MNRAS, 336, 449 [Google Scholar]
  46. Han, Z., Podsiadlowski, P., Maxted, P. F. L., & Marsh, T. R. 2003, MNRAS, 341, 669 [NASA ADS] [CrossRef] [Google Scholar]
  47. Hattori, K., Valluri, M., Bell, E. F., & Roederer, I. U. 2018, ApJ, 866, 121 [CrossRef] [Google Scholar]
  48. Heber, U. 1986, A&A, 155, 33 [NASA ADS] [Google Scholar]
  49. Heber, U. 2009, ARA&A, 47, 211 [Google Scholar]
  50. Heber, U. 2016, PASP, 128, 082001 [Google Scholar]
  51. Heber, U., Edelmann, H., Lisker, T., & Napiwotzk, R. 2003, A&A, 411, L477 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Heber, U., Irrgang, A., & Schaffenroth, J. 2018, Open Astron., 27, 35 [Google Scholar]
  53. Henden, A. A., Levine, S., Terrell, D., & Welch, D. L. 2015, AAS, 225, 336 [Google Scholar]
  54. Hills, J. G. 1988, Nature, 331, 687 [Google Scholar]
  55. Hirsch, H. 2009, Ph.D. Thesis, Friedrich-Alexander University Erlangen-Nürnberg [Google Scholar]
  56. Hirsch, H. A., Heber, U., O’Toole, S. J., & Bresolin, F. 2005, A&A, 444, L61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Iben, I., & Tutukov, A. V. 1984, ApJS, 54, 335 [NASA ADS] [CrossRef] [Google Scholar]
  58. Iben, I., & Tutukov, A. V. 1986, ApJ, 311, 742 [NASA ADS] [CrossRef] [Google Scholar]
  59. Igoshev, A. P., Perets, H., & Hallakoun, N. 2023, MNRAS, 518, 6223 [Google Scholar]
  60. Irrgang, A., Wilcox, B., Tucker, E., & Schiefelbein, L. 2013, A&A, 549, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Irrgang, A., Przybilla, N., Heber, U., et al. 2014, A&A, 565, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Irrgang, A., Kreuzer, S., Heber, U., & Brown, W. 2018a, A&A, 615, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Irrgang, A., Kreuzer, S., & Heber, U. 2018b, A&A, 620, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Irrgang, A., Geier, S., Heber, U., Kupfer, T., & Fürst, F. 2019, A&A, 628, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Irrgang, A., Geier, S., Heber, U., et al. 2021, A&A, 650, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Jeffery, C. S., Miszalski, B., & Snowdon, E. 2021, MNRAS, 501, 623 [Google Scholar]
  67. Jordan, G. C., Perets, H. B., Fisher, R. T., & van Rossum, D. R. 2012, ApJ, 761, L23 [NASA ADS] [CrossRef] [Google Scholar]
  68. Justham, S., Wolf, C., Podsiadlowski, Ph., & Han, Z. 2009, A&A, 493, 1081 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Kawka, A., Vennes, S., O’Toole, S. J., et al. 2015, MNRAS, 450, 3514 [NASA ADS] [CrossRef] [Google Scholar]
  70. Kepler, S. O., Pelisoli, I., Koester, D., et al. 2016, MNRAS, 455, 3413 [NASA ADS] [CrossRef] [Google Scholar]
  71. Kepler, S. O., Pelisoli, I., Koester, D., et al. 2019, MNRAS, 486, 2169 [NASA ADS] [Google Scholar]
  72. Kollmeier, J. A., Zasowski, G., Rix, H. W., et al. 2017, ArXiv e-prints [arXiv:1711.03234] [Google Scholar]
  73. Kuijken, K., Heymans, C., Dvornik, A., et al. 2019, A&A, 625, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Kupfer, T., Ramsay, G., van Roestel, J., et al. 2017a, ApJ, 851, 28 [NASA ADS] [CrossRef] [Google Scholar]
  75. Kupfer, T., van Roestel, J., Brooks, J., et al. 2017b, ApJ, 835, 131 [NASA ADS] [CrossRef] [Google Scholar]
  76. Kupfer, T., Bauer, E. B., Burdge, K. B., et al. 2020a, ApJ, 898, L25 [Google Scholar]
  77. Kupfer, T., Bauer, E. B., Marsh, T. R., et al. 2020b, ApJ, 891, 45 [Google Scholar]
  78. Kupfer, T., Bauer, E., van Roestel, J., et al. 2022, ApJ, 925, L12 [NASA ADS] [CrossRef] [Google Scholar]
  79. Kurucz, R. L. 1996, Model Atmospheres and Spectrum Synthesis (San Francisco: ASP) [Google Scholar]
  80. Latour, M., Randall, S. K., Calamida, A., Geier, S., & Moehler, S. 2018, A&A, 618, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Latour, M., Hämmerich, S., Dorsch, M., et al. 2023, A&A, 677, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Lawrence, A., Warren, S. J., Almaini, O., et al. 2012, VizieR Online Data Catalog: II/314 [Google Scholar]
  83. Lei, Z., Bu, Y., Zhao, J., Németh, P., & Zhao, G. 2019, PASJ, 71, 41 [NASA ADS] [CrossRef] [Google Scholar]
  84. Li, Y.-B., Luo, A.-L., Lu, Y.-L., et al. 2021, ApJS, 252, 3 [NASA ADS] [CrossRef] [Google Scholar]
  85. Lindegren, L., Bastian, U., Biermann, M., et al. 2020, A&A, 649, A4 [Google Scholar]
  86. Lisker, T., & Han, Z. 2008, ApJ, 680, 1042 [NASA ADS] [CrossRef] [Google Scholar]
  87. Liu, Z.-W., Röpke, F. K., Zeng, Y., & Heger, A. 2021, A&A, 654, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. López-Sanjuan, C., Vázquez Ramió, H., & Xiao, K. 2024, A&A, 683, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Luo, Y.-P., Nemeth, P., & Li, Q. 2020, ApJ, 898, 64 [NASA ADS] [CrossRef] [Google Scholar]
  90. Luo, Y.-P., Nemeth, P., Wang, K., Wang, X., & Han, Z. 2021, ApJS, 256, 28 [NASA ADS] [CrossRef] [Google Scholar]
  91. Luo, C., Li, J., Zheng, C., et al. 2024, ArXiv e-prints [arXiv:2404.04835] [Google Scholar]
  92. Marocco, F., Eisenhardt, P. R. M., & Fowler, J. W. 2021, ApJS, 253, 8 [NASA ADS] [CrossRef] [Google Scholar]
  93. Martin, P., Jeffery, C. S., Naslim, N., & Woolf, V. M. 2017, MNRAS, 467, 68 [NASA ADS] [Google Scholar]
  94. Maxted, P. F. L., Marsh, T. R., & North, R. C. 2000, MNRAS, 317, L4 [Google Scholar]
  95. Maxted, P. F. L., Heber, U., Marsh, T. R., & North, R. C. 2001, MNRAS, 326, 1391 [CrossRef] [Google Scholar]
  96. McMahon, R. G., Banerji, M., & Gonzalez, E. 2021, VizieR Online Data Catalog: II/367 [Google Scholar]
  97. Mereghetti, S., Tiengo, A., Esposito, P., et al. 2009, Science, 325, 1222 [NASA ADS] [CrossRef] [Google Scholar]
  98. Mereghetti, S., Pintore, F., Rauch, T., et al. 2021, MNRAS, 504, 920 [NASA ADS] [CrossRef] [Google Scholar]
  99. Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 [Google Scholar]
  100. Napiwotzki, R. 1999, A&A, 350, 101 [Google Scholar]
  101. Naslim, N., Jeffery, C. S., Hibbert, A., & Behara, N. T. 2013, MNRAS, 434, 1920 [NASA ADS] [CrossRef] [Google Scholar]
  102. Nemeth, P., Ziegerer, E., Irrgang, A., et al. 2016, ApJ, 821, L13 [NASA ADS] [CrossRef] [Google Scholar]
  103. Neunteufel, P. 2020, A&A, 641, A52 [EDP Sciences] [Google Scholar]
  104. Neunteufel, P., Yoon, S.-C., & Langer, N. 2019, A&A, 627, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Neunteufel, P., Kruckow, M., Geier, S., & Hamers, A. S. 2021, A&A, 646, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  106. Neunteufel, P., Preece, H., Kruckow, M., et al. 2022, A&A, 663, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  107. Nomoto, K. 1982, ApJ, 257, 780 [Google Scholar]
  108. Onken, C. A., Wolf, C., Bessell, M. S., et al. 2019, PASA, 36, 33 [NASA ADS] [CrossRef] [Google Scholar]
  109. Paczynski, B. 1971, Acta Astron., 21, 1 [NASA ADS] [Google Scholar]
  110. Pakmor, E., Zenati, Y., Perets, H. B., & Toonen, S. 2021, MNRAS, 503, 4734 [Google Scholar]
  111. Pelisoli, I., Vos, J., Geier, S., Schaffenroth, V., & Baran, A. S. 2020, A&A, 642, A180 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Pelisoli, I., Neunteufel, P., Geier, S., et al. 2021, Nat. Astron., 5, 1052 [NASA ADS] [CrossRef] [Google Scholar]
  113. Pereira, C. 2011, Ph.D. Thesis, Queen’s University Belfast, UK [Google Scholar]
  114. Perets, H. B. 2009, ApJ, 698, 1330 [NASA ADS] [CrossRef] [Google Scholar]
  115. Qi, Z., Yu, L., Bucciarelli, B., et al. 2015, AJ, 150, 137 [Google Scholar]
  116. Raddi, R., Hollands, M. A., Gänsicke, B. T., et al. 2018a, MNRAS, 479, L96 [NASA ADS] [CrossRef] [Google Scholar]
  117. Raddi, R., Hollands, M. A., Koester, D., et al. 2018b, ApJ, 858, 3 [NASA ADS] [CrossRef] [Google Scholar]
  118. Raddi, R., Hollands, M. A., Koester, D., et al. 2019, MNRAS, 489, 14 [Google Scholar]
  119. Ramspeck, M., Heber, U., & Edelmann, H. 2001, A&A, 379, 235 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  120. Renzo, M., Zapartas, E., de Mink, S. E., et al. 2019, A&A, 624, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  121. Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  122. Roeser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440 [Google Scholar]
  123. Schaffenroth, V., Pelisoli, I., Barlow, B. N., et al. 2022, A&A, 666, A182 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Schlafly, E. F., Meisner, A. M., & Green, G. M. 2019, ApJS, 240, 30 [Google Scholar]
  125. Scholz, R. D. 2024, A&A, 685, A162 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  126. Scholz, R. D., Heber, U., Heuser, C., et al. 2015, A&A, 574, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  127. Shen, K. J., Boubert, D., Gänsicke, B. T., et al. 2018, ApJ, 865, 15 [Google Scholar]
  128. Tauris, T. M. 2015, MNRAS, 448, L6 [Google Scholar]
  129. Tian, H.-J., Gupta, P., Sesar, B., et al. 2017, ApJS, 232, 4 [NASA ADS] [CrossRef] [Google Scholar]
  130. Tillich, A., Heber, U., Geier, S., et al. 2011, A&A, 527, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  131. Vennes, S., Kawka, A., O’Toole, S. J., Nemeth, P., & Burton, D. 2012, ApJ, 759, L25 [NASA ADS] [CrossRef] [Google Scholar]
  132. Vennes, S., Nemeth, P., Kawka, A., et al. 2017, Science, 357, 680 [Google Scholar]
  133. Vos, J., Nemeth, P., Vuckovic, M., Østensen, R. H., & Parsons, S. 2018, MNRAS, 473, 693 [NASA ADS] [CrossRef] [Google Scholar]
  134. Vos, J., Bobrick, A., & Vuckovic, M. 2020, A&A, 641, A163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  135. Wang, B., & Han, Z. 2009, A&A, 508, L27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  136. Webbink, R. F. 1984, ApJ, 277, 355 [NASA ADS] [CrossRef] [Google Scholar]
  137. Whelan, J., & Iben, I. 1973, ApJ, 186, 1007 [CrossRef] [Google Scholar]
  138. Werner, K., Reindl, N., Dorsch, M., et al. 2022, A&A, 658, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  139. Werner, K., Reindl, N., Rauch, T., et al. 2024, A&A, 682, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  140. Woosley, S. E., & Weaver, T. A. 1994, ApJ, 423, 371 [Google Scholar]
  141. Xu, Z., & Zhaoji, J. 2005, VizieR Online Data Catalog: II/262 [Google Scholar]
  142. Yanny, B., Rockosi, C., Newberg, H. J., et al. 2009, AJ, 137, 4377 [Google Scholar]
  143. Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44 [Google Scholar]
  144. Zhang, X., & Jeffery, C. S. 2012, MNRAS, 419, 452 [NASA ADS] [CrossRef] [Google Scholar]
  145. Ziegerer, E., Volkert, M., Heber, U., et al. 2015, A&A, 576, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  146. Ziegerer, E., Heber, U., Geier, S., et al. 2017, A&A, 601, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Additional material

The Tables A.1 and A.2 give an overview of the available spectra and list the photometric surveys. Table A.3 lists the classifications and spectroscopic parameters, Table A.4 lists the astrometric data from Gaia EDR3, and Table A.5 the results of the SED analysis. Tables A.6 and A.7 provide the kinematic parameters and the parameters at the last Galactic disc crossing, respectively.

Table A.1.

Spectroscopic data overview

Table A.2.

Photometric surveys used for the SED analysis

Table A.3.

Classifications and spectroscopic parameters

Table A.4.

Astrometric parameters

Table A.5.

Results of the SED analysis

Table A.6.

Kinematic parameters of the fast star candidates

Table A.7.

Parameters of the fast star candidates at the last Galactic disc passage

All Tables

Table 1.

Summary of follow-up observations during the course of the MUCHFUSS project.

Table 2.

Comparison of spectroscopic versus parallax ϖ based distances. The inversion accounts for asymmetry and zero-point offset (Lindegren et al. 2020).

Table A.1.

Spectroscopic data overview

Table A.2.

Photometric surveys used for the SED analysis

Table A.3.

Classifications and spectroscopic parameters

Table A.4.

Astrometric parameters

Table A.5.

Results of the SED analysis

Table A.6.

Kinematic parameters of the fast star candidates

Table A.7.

Parameters of the fast star candidates at the last Galactic disc passage

All Figures

thumbnail Fig. 1.

Proper motion measurements in right ascension (left panel) and declination (right panel). Gaia proper motions from Early Data Release 3 (EDR3, Gaia Collaboration 2021) are plotted against ground-based measurements from SDSS DR9 (Ahn et al. 2012, open circles), the PPMXL catalogue of positions and proper motions on the ICRS (Roeser et al. 2010, open squares), the Whole-Sky USNO-B1.0 Catalogue (Monet et al. 2003, open diamonds), the Absolute Proper motions Outside the Plane catalogue (APOP, Qi et al. 2015, open downward triangles), Gaia-PS1-SDSS (GPS1, Tian et al. 2017, open upward triangles), and Hot Stuff for One Year (HSOY, Altmann et al. 2017, crosses). PMs from the Fourth U.S. Naval Observatory CCD Astrograph Catalogue (UCAC4, Zacharias et al. 2013) were also used, but they are not plotted here, because only a few stars have measurements in this catalogue. The averages of the ground-based values from those surveys are plotted as filled diamonds. The ground-based measurements by Tillich et al. (2011) and Ziegerer et al. (2017) are plotted as open dark grey diamonds.

In the text
thumbnail Fig. 2.

Excerpts of BOSS spectral fits of He-poor hot subdwarfs. From top to bottom: SDSSJ124248.89+133632.6 (sdB), SDSSJ124310.58+343358.4 (sdOB), SDSSJ124819.08+035003.2 (sdO). In the upper panels the best-fit synthetic spectrum (red) is shown along with the observed one. The lower panels give residual χ, whereby the long-dashed horizontal lines mark the zero values and the short-dashed horizontal lines mark the deviations in terms of 1σ. Parts of the spectra, which have been excluded from the fits, are marked in grey.

In the text
thumbnail Fig. 3.

Same as Fig. 2, but for He-rich hot subdwarfs. From top to bottom: SDSSJ120352.24+235343.3 (iHe-sdOB), SDSSJ222515.34-011156.8 (eHe-sdOB), SDSSJ080833.76+180221.8 (eHe-sdO).

In the text
thumbnail Fig. 4.

Same as Fig. 2, but showing fits of the WHT/ISIS spectra of hot subdwarfs. From top to bottom: SDSSJ130543.97+115840.8 (sdB), SDSSJ221728.35+121642.6 (sdOB), SDSSJ120352.24+235343.3 (iHe-sdOB).

In the text
thumbnail Fig. 5.

Teff − log g diagram of sample of fast stars. The colour scales with the helium abundance from light orange to dark purple. The light grey symbols mark the results of the preliminary spectroscopic analysis on which the target selection was based. They are connected with the corresponding values from the new homogeneous analysis. The zero-age and terminal-age EHBs for a subsolar metallicity of −1.48 (solid grey lines) were interpolated from evolutionary tracks by Dorman et al. (1993). The helium main sequence (dashed line) is taken from Paczynski (1971).

In the text
thumbnail Fig. 6.

Teff − log n(He)/n(H) diagram of sample compared to the preliminary analysis (see Fig. 5). Solar helium abundance is marked by the solid horizontal line, while the dotted line marks the transition between intermediate and extreme helium abundance.

In the text
thumbnail Fig. 7.

log p distribution of sample. The significance level log p < −4.0 is marked by the vertical line.

In the text
thumbnail Fig. 8.

Fits of SEDs of selected programme stars. Each plot consists of two panels; the upper one compares the observed fluxes to the synthetic SED. To ease the slope of the distribution the flux is multiplied by the wavelength to the power of three. Photometric fluxes are displayed as coloured data points (for details see Table A.2) with their respective uncertainties and filter widths (dashed lines). The best-fit models are drawn as full grey lines. The lower panels give reduced χ to demonstrate the quality of the fit. Example fits are selected to demonstrate relevant issues. In each row, the stars are arranged with increasing temperature. Top row: He-weak sdB stars: cool sdB J123953.52+062853.0, which is the most distant programme star; J135651.26+155810.4, the sdB with the largest (positive) heliocentric RV; and J123137.56+074621.7. Middle row: sdOB stars: J164419.44+452326.7; the low gravity J130543.97+115840.8; and the highly reddened J180313.45+234000.1. Bottom row: He-sdO stars: iHe-sdOB J204358.55–065025.8; the best observed iHe-sdOB star J022422.21+000313.5; and the eHe-sdO J205030.39–061957.8, which was considered a spectroscopic twin to US 708 by Ziegerer et al. (2017). The strong colour excess of J180313.45+234000.1 and J205030.39–061957.8 is most obvious from the strength of the 2200 Å UV feature effecting the GALEX NUV band.

In the text
thumbnail Fig. 9.

Distribution of distances of our sample from the Galactic plane. The dashed vertical lines mark the scale height of the Galactic thin disc; the solid vertical lines mark the scale height of the thick disc.

In the text
thumbnail Fig. 10.

Left panel: Toomre diagram of sample. The 2σ contours of the thin disc (grey), the thick disc (green), and the halo (blue) are constructed from the U, V, and W distributions of Anguiano et al. (2020). Objects with halo kinematics are marked in blue; the few objects also consistent with thick-disc kinematics are marked in green. Right panel: e − LZ diagram of the sample using the same colour coding for the stars. The green and grey 2σ contours in the upper left denote the locations of thin- and thick-disc stars, respectively, based on the GALAH + Gaia EDR3 sample of Buder et al. (2021). The dashed lines mark prograde (positive LZ) and retrograde (negative LZ) directions of Galactic rotation, respectively.

In the text
thumbnail Fig. 11.

Upper panel: Points of last Galactic disc crossings of stars in our sample in the XY plane. The two stars with the largest distances have been omitted here for better visualisation. Lower panel: Distribution of the galactocentric radii at disc crossing for all stars of the sample. Stars that could have been ejected from the disc are marked in dark grey, and the others are in light grey.

In the text
thumbnail Fig. 12.

Ejection or disc-crossing velocities plotted against galactocentric distance.

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.