# Candidate bow shock stars: Spectral analysis and environment

D. Russeil, A. Zavagno, J. C. Bouret, C. Adami

► To cite this version:

D. Russeil, A. Zavagno, J. C. Bouret, C. Adami. Candidate bow shock stars: Spectral analysis and environment. Astronomy & Astrophysics - A&A, 2025, 698, pp.A64. 10.1051/0004-6361/202452262. hal-05096953

HAL Id: hal-05096953

https://hal.science/hal-05096953v1

Submitted on 4 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Candidate bow shock stars: Spectral analysis and environment

D. Russeil<fnref n="1" />,*, A. Zavagno<fnref n="1" />,<fnref n="2" />, J. C. Bouret<fnref n="1" />, and C. Adami<fnref n="1" />

Received 16 September 2024 / Accepted 4 April 2025

## ABSTRACT

**Aims.** Infrared bow shocks are arc-shaped structures located ahead of a star and generally observed at mid- to far-IR wavelengths. They are thought to result from the interaction of the stellar wind with the ambient interstellar medium and are typically (but not always) related to runaway stars. However, the formation of bow shocks seems to be dominated by local environmental factors rather than stellar motion. In this context, we aim to probe the links between bow-shock driving stars and their environment.

**Methods.** We observed 47 bow shock driving star candidates with the Multi-purpose InSTRument for Astronomy at Low-resolution (MISTRAL) spectro-imager at Haute-Provence Observatory (OHP) in the 420–800 nm range to perform spectral classification of the candidate stars. In parallel, we evaluated the transverse motion of stars from GAIA DR3 in order to determine whether they are runaways. We then characterised the bow shock environmental conditions.

**Results.** We find that among the 47 candidates we have 3 unclassifiable stars (suspected to be G- or K-type stars), 3 M- or K-type stars, 2 A-type stars, 10 O stars, and 29 B (mainly giant and supergiant) stars. We find that 17 stars (among the 37 with determined transverse velocity) are runaways, among which only 7 have their transverse velocity aligned to the bow-shock axis. This suggests that runaway is not the only origin for bow shock formation. We note the diversity of environments where bow shocks are observed: stellar associations, a cluster, and H II regions. For most stars, the origin of the bow shock is not clear; however, the 11 bow shocks observed in the Cygnus OB stellar association suggest that the ISM conditions in such regions favour bow shock observability. We also identify that the bow shock ahead of the star ionising the H II region Sh2-135 could be produced by a photoevaporated flow of about 16 km/s coming from the H II region molecular cloud’s interface. Finally, for six stars we were able to identify the cluster from which they were ejected and determine the ejection process (dynamical ejection from star cluster or binary supernovae scenarios).

**Conclusions.** The formation of bow shocks seems to be dominated by local environmental factors rather than stellar motion.

**Key words.** circumstellar matter – stars: general – stars: winds, outflows – HII regions

## 1. Introduction

Infrared stellar bow-shock nebulae (SBNs) are arc-shaped structures located ahead of a star and generally observed at 24 µm (and sometimes at 8 µm) by the Spitzer Space Telescope or at 22 µm by the Wide-field Infrared Survey Explorer (WISE) Telescope (e.g. Churchwell et al. 2009; Carey et al. 2009; Wright et al. 2010; Kobulnicky et al. 2016; Jayasinghe et al. 2019). They are expected to arise from the interaction between the ambient interstellar medium (ISM) and the star wind. Bow shocks caused by a runaway star and in situ bow shocks are the two primary categories.

In the in situ scenario, a bow shock is produced at the interaction layer with the stellar wind as hot gas from an H II star-forming region (Henney & Arthur 2019a) rushes towards the star. For instance, this process is expected to explain the 24 µm bow shock observed close to the exciting star but inside the circular photodissociation region (PDR), traced at 8 µm, of the H II region RCW120 (Jayasinghe et al. 2019) and the bow shocks facing the very massive Westerlund 2 cluster in RCW49 (Povich et al. 2008; Jayasinghe et al. 2019). More precisely, in RCW120, the bubble rim is opened on its upper northeastern side while, due to its location at the bottom part of the bubble, the ionising star strongly irradiates the nearby molecular cloud’s interface, entraining the ablated gas and dust in a photoevaporative flow

directed towards the open side of the bubble. Passing by the ionising star, it interacts with the stellar wind to form the observed 24 µm bow shock structure.

When the star is a runaway, the bow shock is caused by the star’s supersonic motion relative to the ambient medium. Large volumes of gas and dust are swept up by the runaway OB star’s powerful wind during this process, and the material piles up ahead to produce the bow shock. The dust and gas swept along by the bow shock are then heated by the stellar and shock-excited radiation, and the dust re-radiates the energy as a mid- to far-infrared excess emission. In this way, the bow shock develops as an arc-shaped structure, with bows pointing in the same direction as the stellar motion in the local rest frame. The wind is confined by the ram pressure of the ISM at distances from the star determined by momentum balance (e.g. Martinez et al. 2023). Nevertheless, only a small fraction of runaway stars produce detectable bow shocks (e.g. Gvaramadze et al. 2011), indicating that the ISM properties have a significant impact on the appearance of bow shocks. In particular, Kobulnicky & Chick (2022) conclude that the formation of SBNs is dominated by local environmental factors rather than stellar motion.

In this paper, we present a spectroscopic follow-up in the 420–800 nm range of a sample of 47 bow-shock driving star candidates in order to characterise their O or B spectral type and then verify their runaway status using Gaia DR3 astrometric data. This is done in the context of probing the environmental relations of bow-shock driving stars. Section 2 of this paper

<footnote n="1">¹ Aix-Marseille Univ., CNRS, CNES, LAM, 13388 Marseille, France</footnote>

<footnote n="2">² Institut Universitaire de France, 1 rue Descartes, 75005 Paris, France</footnote>

<footnote>* Corresponding author.</footnote>
presents the observations and data reduction. Section 3 presents the ancillary datasets used for this study and Sect. 4 the spectral classification. Section 5 presents the results obtained on extinction and infrared excess and Sect. 6 the ones on the physical parameters. The results are discussed in Sect. 7. Conclusions and perspectives are given in Sect. 8.

## 2. Observation and data reduction

Kobulnicky et al. (2016) have created the largest Spitzer/WISE bow shock catalogue, which has subsequently been completed by Jayasinghe et al. (2019). They identify the putative central star at the origin of the identified bow shocks from position and Two Micron All-Sky Survey (2MASS) colours. These catalogues also provide the measured standoff distance (distance between the central star and the arc-shaped nebula apsis), the position angle of the bow shocks, and an assigned environment code based on the local surroundings (FB=facing bright-rimmed cloud, FH=facing H II region, H=inside H II region, and I=isolated).

For this paper, we selected the stellar bow shock candidates that are observable from Haute-Provence Observatory (OHP) with the Multi-purpose InSTRument for Astronomy at Low-resolution (MISTRAL) instrument (star brighter than V = 16 mag and declination larger than -5°), with right ascension between 18h30 and 23h, and close to the Galactic plane (b < 2.5°). Combining the Kobulnicky et al. (2016) and Jayasinghe et al. (2019) catalogues, these constraints led to a sample of 70 stars of which 47 were observed (among which 32 are isolated, 6 are inside the H II region, 5 face an H II region, and 4 face a bright-rimmed cloud). They belong mainly to the Jayasinghe et al. (2019) catalogue.

MISTRAL is the new Faint Object Spectroscopic Camera mounted at the folded Cassegrain focus of the 1.93m telescope whose full technical description can be found in Schmitt et al. (2024). Our primary goal being to establish and specify the spectral type of the stars driving the SBN, we then used the blue disperser, covering the 420 nm to 800 nm range. With this disperser, the spectral resolution is R ~ 700 at 600 nm, giving us low-resolution spectra (instrumental FWHM ~ 8.5 Å). The typical exposure time was 30 min to reach a S/N of 10, but because the spectral response of the receiver falls drastically below ~490 nm, for most stars the most useful part of the spectrum is above 490 nm. The spectral calibration was carried out from an XeHgAr spectral lamp observed just before and just after each science exposure. The data reduction was performed using the dedicated pipeline<fnref n="1" /> which is based on the Automated SpectroPhotometric Image REDuction (Lam et al. 2023) software. The spectra were rebinned to a 1 Å step and then normalised. An example of a spectrum is shown in Schmitt et al. (2024).

Raw data are available online<fnref n="2" /> in a database hosted by the Centre de donnéeS Astrophysiques de Marseille (CeSAM) service<fnref n="3" /> at LAM.

## 3. Complementary datasets

We complemented the spectral type characterisation with photometric information. For that, we crossmatched our sample (the nearest object within a 3'' cone search radius) with the Gaia Early Data Release 3 (EDR3; Gaia Collaboration 2021),

IGAPS (Monguió et al. 2020), Pan-STARRS-DR1 (Chambers et al. 2016), 2MASS (Cutri et al. 2003), Spitzer/GLIMPSE (Churchwell et al. 2009), and WISE (Cutri et al. 2013) point source catalogues. For U, B, and V magnitudes, we performed crossmatching with the Simbad<fnref n="4" /> database. The IGAPS point source catalogue gives u, g, r, i, and Hα magnitudes up to limiting magnitudes of 21.5, 22.4, 21.5, 20.4, and 20.5, respectively (Monguió et al. 2020; Drew et al. 2005; Groot et al. 2009). For u, g, and r, the saturation magnitudes are 14.5, 14, and 13, respectively, and the seeing is between 1.5" and 1.1". Pan-STARRS-DR1 provides g, r, i, z, and y broad band magnitudes up to limiting magnitudes of 23.3, 23.2, 23.1, 22.3, and 21.4, respectively. The J, H, and Ks magnitudes are from the 2MASS full-sky survey (seeing ~2.5"), the 3.6, 4.5, and 5.8, 8 µm magnitudes are from the Spitzer/GLIMPSE survey, and the WISE survey gives point source magnitudes at 3.4, 4.6, 12, and 22 µm (W1 to W4 bands). We also looked at the MIPS 24 µm (Gutermuth & Heyer 2015) point source catalogue data, but only two objects have a counterpart. From Gaia EDR3, we retrieved the proper motions and the parallaxes (with errors) in addition to the G magnitude, the colour excess (E(BP-RP)), and the modelled stellar temperature.

## 4. Spectral classification

The spectral classification is commonly performed from lines in the 400–500 nm domain but, due to the low S/N we have in this spectral range, we performed the spectral classification from lines in the red part of the spectrum. Unfortunately, for OB stars, only a few papers (Conti 1974; Didelon 1982; Lennon et al. 1993; Jaschek et al. 1994; Leone & Lanzafame 1998; Martins et al. 2005b; Kobulnicky et al. 2012; Martins & Palacios 2021) give predicted and/or observed equivalent widths (EWs) for H, HeI, and HeII lines in the red domain. In our spectra, the Hα, Hβ, HeI/HeII667.8, HeI706.5, HeI587.6, HeI541.1, HeI492.1, and HeII468.6 lines are the most prominent. We then followed a classification scheme based mainly on these lines. The first step in our classification scheme was to distinguish O from B stars. For that, we assumed that O stars (down to B0.5) are characterised by prominent HeII lines (e.g. Evans et al. 2004). Once an O star was identified, we used the HeI706.5 over HeII541.1 EW ratio from Martins & Palacios (2021, their Fig. 11) to specify the spectral type and the HeII 468.6 EW (figures 3 and 5 in Martins 2018) to define the luminosity class. For B stars, we based the classification on the HeI706.5 EW (Fig. 2 from Jaschek et al. 1994). The typical uncertainty on the HeI706.5 EW suggests a typical uncertainty of one luminosity class and up to two spectral types. In addition, for a given HeI706.5 EW, up to three solutions are possible, so to find the spectral type and luminosity class we first used the HeI668.0 EW (Fig. 1 from Jaschek et al. 1994) and the HeI587.6 and/or HeI501.6 EWs (Leone & Lanzafame 1998), and then complemented these with the Hα and Hβ EWs from Didelon (1982). We note that, frequently, the Balmer EWs do not converge to the same spectral type as suggested by HeI lines. To specify the classification, when the S/N allows it (for 12 stars out of 47), for B stars we looked at the OII464.0–465.0, NII 463.1, and SiIII 455.2–456.8–457.5 features that are typically observed in B supergiants (e.g. Fig. 3.5 from Gray & Corbally 2009) and, with a CIII/OII blend, are present from late O-type to early B-type stars (with a maximum at B1Ia).

To test our classification scheme, the spectral classification was carried out on stars of known type. We complemented the

<footnote n="1">¹ http://wwwobs-hp.fr/guide/mistral/MISTRAL_spectrographCamera.shtml#H5</footnote>

<footnote n="2">² https://data.lam.fr/mistral/home</footnote>

<footnote n="3">³ https://www.lam.fr/service/cesam/presentation/</footnote>

<footnote n="4">⁴ https://simbad.unistra.fr/simbad/</footnote>
comparison stars we observed by retrieving OB stellar spectra (and the given spectral classification) from the MILES, STELIB, and X-shooter spectral libraries<fnref n="5" />. Tables A.1 and A.2 summarise our EW measurements and classification. We note that, for O stars, the spectral type we find is in good agreement with the one given in the library. For B stars, the agreement is not as good, with a typical difference of two spectral types. In addition, because of its early B-type, HD 34816 has been classified following the O-type and the B-type schemes, illustrating the large uncertainty for B0/0.5 stars.

We then performed the spectral classification for our sample. The measurements and results are shown in Table A.3. We evaluated the EWs, repeating the measures for slightly different wavelength ranges, as the main source of EW uncertainty is the choice of the continuum level. Several stars have a spectral type in the literature (from Simbad or Chick et al. 2020) that is globally in agreement with our classification. We find that, among the 47 candidates, we have: 3 cool stars, 1 A-type star, 10 O stars and 31 B (mainly giant and supergiant) stars. Three stars are not classifiable because no line other than Hα is identifiable (noisy spectrum). We performed a crossmatching with the recent Gaia DR3 astrophysical parameters table<fnref n="6" />, which gives a spectral type for 46 of our 47 stars. This allows us to confirm the later spectral type of cold stars, also suggesting that the three unclassified stars (#3, #4, and #7) are probably cold.

Then, as already observed (e.g. Kobulnicky et al. 2019; Chick et al. 2020; Kobulnicky & Chick 2022) and as expected by prevailing models (e.g. Meyer et al. 2014; Henney & Arthur 2019a; Raga et al. 2022), we also find that bow shocks are mainly driven by hot stars.

In parallel, we note that only two stars (#26 and #28) show a clear Hα line in emission, which, given our low spectral resolution, suggests that most of the stars in our sample must have low mass-loss rates, as Marcolino et al. (2009) show that Hα is in emission for $log \dot{M} \ge -6.30$ only. For lower rates, the Hα profiles are in absorption and, therefore, only a small portion of the Hα profile is refilled by the wind, making it undetectable within the broad photospheric profile. At lower wind densities (roughly, below a few times $10^{-8} M_{\odot} \text{yr}^{-1}$ in the case of O/early B stars), the Hα line is even unsuitable for determining a mass-loss rate (e.g. Peppel 1984; Mokiem et al. 2007; Puls et al. 2008).

## 5. Extinction and infrared excess

To subsequently calculate the stellar luminosity needed in Sect. 6.2, we determined the $A_V$ extinction. To do this, we used the stellar spectral energy distribution (SED). Because of the effective temperature of the hot stars (between 10 000 and 50 000 K), the theoretical photospheric emission should follow $F_V \propto v^{\beta}$ with $\beta \sim -1.8$ to $-2$. To determine the expected value for $\beta$, we fitted synthetic spectra, extracted from the Pollux<fnref n="7" /> database, for stars with a temperature between 12 000 and 35 000 K and log(g) between 2.3 and 4.2. We find that the expected $\beta$ is between $-1.90$ and $-1.54$ and, therefore, is not too sensitive to stellar parameters. We then performed a power-law fit to the near-IR part (1 to 8 µm) of the SED and adjusted the extinction so that the data points (especially the one in the optical part with $\lambda < 1$ µm) of the extinction-corrected SED followed the power law (Fig. B.1). This method being appropriate for hot

<figure><img src="image_2.png" /><figcaption>Fig. 1. 2MASS colour-colour diagram. The red, green, and blue lines represent the main sequence, supergiant and giant locations, respectively. These were built from the PARSEC evolutionary tracks<fnref n="8" /> with Solar metallicity. The dashed grey lines are 10 mag length reddening vectors. The upward and downward triangles are the classified cool stars (#1, #9, and #30) and the suspected cool ones (#3, #4, and #7), respectively.</figcaption></figure>

stars only, the cool stars #1, #9, and #30 were excluded from this process, while the $A_V$ values of stars #3, #4, and #7 (suspected to be cold stars, as well) must be taken with caution. We find $\beta$ between $-2.09$ and $-1.58$, in agreement with the theoretical expected range for hot stars. The $A_V$ values are listed in Table A.4.<fnref n="8" />

The infrared excess from OB stars is commonly considered to be the contribution from the ionised stellar wind or circumstellar dust (Hovhannessian & Hovhannessian 2001; Siebenmorgen et al. 2018; Deng et al. 2022). The mechanism explaining the IR excess is either the free-free emission (which gives a power-law continuum emission from IR to radio bands) or dust thermal emission. Deng et al. (2022), from a large OB star sample, show that about 8% are identified to have circumstellar IR excess emission. They show that the IR excess can be better explained by free-free emission in the stellar wind for ~33% of them, by synchrotron emission for ~32% of them, and by both models for ~30% of them, while the wind-dust model works better for only ~4% of them. They also show that for the wind-dust objects, a parsec-scale dusty envelope with a low dust temperature is identified, which implies a large-scale circumstellar dust halo. Because the bow shock around a star is expected to be due to wind interaction with the interstellar medium, we propose to check if our star sample shows such IR excess. From the 2MASS colour-colour plot (Fig. 1) we can note that, within the error bars, our sample exhibits no near-infrared excess and shows typical

<footnote n="5">⁵ http://svo2.cab.inta-csic.es/theory/libraries/</footnote>

<footnote n="6">⁶ See https://gea.esac.esa.int/archive/documentation/GDR3/</footnote>

<footnote n="7">⁷ https://pollux.oreme.org/</footnote>

<footnote n="8">⁸ http://stev.oapd.inaf.it/cgi-bin/cmd</footnote>
reddening less than 10 mag in agreement with the $A_V$ measure-
ments. In parallel, we see that cool stars are clearly located in the
upper part of the diagram.

In addition, looking at the SED (Fig. B.1), we note that most
of the stars exhibit an IR excess above 10 µm. However, this
excess can be attributed to the bow shock itself. We checked the
WISE 22 µm and/or the MIPS 24 µm images to see if the stars
effectively emit at these wavelengths. Only 11 stars (#6, #7, #11,
#20, #21, #22, #24, #28, #33, #34, and #45) have a point source
counterpart on the image, while eight (#10, #12, #13, #23, #29,
#32, #41, and #44) are clearly dominated or contaminated by the
bow-shock emission. Lastly, stars #20, #24, and #28 are the best
mid-IR excess stellar candidates. Among them, star #28 is the
only one to present an $H\alpha$ line in emission corroborating the
presence of a wind.

**6. Physical parameters**

Because our goal is to determine as independently as possible the
different quantities involved in the momentum balance equation
which characterises SBNs, we derived the physical parameters
for the stars, bow shocks, and their local environment. The SBNs
and the stellar transverse velocity vectors are shown in Fig. C.1,
while the physical quantities discussed below are summarised in
Table A.4.

6.1. *Stellar distance and transverse velocity*

With regard to the distance of the star, we adopted the distance
estimated from parallaxes by Bailer-Jones et al. (2021) or, if not
available (as was the case for star #25), the usual spectropho-
tometric distance, following Russeil (2003). We then have star,
and hence bow-shock, distances ranging from 555 pc to 9328 pc,
with 76.6% of them between 1300 and 5000 pc. We used the
proper motions and the parallaxes retrieved from Gaia EDR3 to
calculate the Galactic transverse velocity components of the star
with respect to the local interstellar medium, following Russeil
et al. (2020). To ensure good quality of the transverse velocity
vector (Vt), it was calculated only for stars with a parallax $\pi > 0$,
a relative uncertainty $\sigma_\pi / \pi \le 0.2$, and RUWE$^9 \le 1.4$ (Fabricius
et al. 2021).

With these constraints, we could calculate Vt for 37 of the
47 stars and find Vt between 2.9 and 120.6 km/s (median value
23 km/s). Because the radial velocity is not known for most of
the stars, the transverse velocity, Vt, is then a lower limit for the
true stellar velocity. However we adopt it as the space velocity
of the star with respect to the ambient medium. Mackey et al.
(2021), from simulations, show that the arcuate shape of the SBN
is prominent only if it is seen edge-on (inclination angle greater
than 45°). Since they were identified on the basis of their arcuate
morphology, we can assume that the peculiar velocity is domi-
nated by the transverse component rather than the radial one. To
check this hypothesis, we considered the 20 stars for which the
radial heliocentric velocity is given in the SIMBAD$^{10}$ database,
Gaia DR3 database, or Chick et al. (2020). After heliocentric
to Local Standard of Rest (LSR) conversion, the peculiar radial
velocity was calculated and quadratically added to Vt. We find
that the ratio between the full space velocity and Vt is on aver-
age 1.54, with only three stars showing a ratio larger than 2 (stars

#1, #33, and #36). We can then conclude that Vt is a representa-
tive measure of the stellar velocity (with respect to the ambient
medium).

The assumption that bow shocks require supersonic star-ISM
relative velocity suggests that runaway stars are good candidates.
The velocity threshold adopted to define a runaway star depends
on the authors, but also on whether the considered velocity is
the full, the radial, or the transverse one. The velocity threshold
is usually between 25 km/s and 40 km/s (Blaauw 1961; Tetzlaff
et al. 2011; Kobulnicky & Chick 2022). Following Kobulnicky
& Chick (2022), we adopted a threshold of 25 km/s and find
that 17 stars (among the 37 with a determined Vt) are runaways,
among which 15 are isolated. In addition, one can note that only
five stars have a Vt vector pointing at more than 45° away from
the Galactic plane (among which only two are runaways). This
has already been observed for the morphology and kinematics
position angles of previous samples of SBN central stars
(Kobulnicky et al. 2016; Kobulnicky & Chick 2022).
Kobulnicky & Chick (2022) explain that by a selection
effect, arguing that OB stars with a velocity vector perpendic-
ular to the Galactic plane quickly reach the less dense region
of the disk, which minimises the probability of producing a
detectable bow shock.

6.2. *Stellar mass-loss rate*

From the spectral type, we adopted the stellar temperature
($T_{\text{eff}}$) and mass from Schmidt-Kaler (1982) and de Jager &
Nieuwenhuijzen (1987) for B stars and from Martins et al.
(2005a) for O stars. With the $A_V$ determined in Sect. 5, we deter-
mined the Ks and G absolute magnitudes (the $A_V$ to $A_{Ks}$ and
$A_G$) conversion factors are from Wang & Chen 2019). From the
MIST$^{11}$ database, we evaluated (retrieving the values for solar
metallicity and log g between 2 and 4.5) the bolometric correc-
tions (BC) to apply to the Ks and G absolute magnitudes. The BC
appears to mainly depend on $T_{\text{eff}}$, while the log g dependence
can be taken into account as a typical uncertainty of 0.1 mag.
The stellar luminosity was calculated adopting the Solar bolo-
metric magnitude from Willmer (2018) and the values obtained
from Ks and G magnitudes were averaged (if there was only one
value, that was adopted). We then used relations from Vink et al.
(2001) to evaluate the mass-loss rate ($\dot{M}$) from the adopted lumi-
nosity, mass, and $T_{\text{eff}}$ (and assuming Solar metallicity). For our
sample, $\dot{M}$ is between $2 \ 10^{-14} M_{\odot} \text{yr}^{-1}$ and $1.8 \ 10^{-5} M_{\odot} \text{yr}^{-1}$.
This is in agreement with the $\dot{M}$ found by Peri et al. (2012),
Kobulnicky et al. (2018) and Kobulnicky et al. (2019) for other
samples of SBN star candidates. We have to keep in mind that,
while there is a general agreement between Vink et al. (2001)
and $\dot{M}$ observed by various techniques (in particular the one
developed by Kobulnicky et al. 2018), $\dot{M}$ evaluations suffer from
uncertainties that impact the quantities derived in Sect. 7. For
example, with an $\dot{M}$ on average a factor of 2.7 below the Vink
et al. (2001) prescription, as Kobulnicky et al. (2019) find, we
can estimate that Vt would be 60% smaller when it is calculated
from Eq. (1).

6.3. *Density of the local interstellar medium*

Kobulnicky et al. 2018 developed a method to measure the local
interstellar medium based on the infrared surface brightness of
the nebula. Here we propose to estimate the ambient density in a
way independent of the shock itself. For this we used the HI4PI

<footnote n="9">⁹ RUWE is the Renormalised Unit Weight Error. It is a measure of the reliability of the astrometric solution.</footnote>

<footnote n="10">¹⁰ See http://simbad.cds.unistra.fr/simbad/</footnote>

<footnote n="11">¹¹ See https://waps.cfa.harvard.edu/MIST/</footnote>
datacubes (HI4PI Collaboration 2016). These have an angular resolution of 16.2' and a spectral resolution of 1.28 km/s. From the adopted distance of the star, we first evaluated the $V_{LSR}$ (the radial velocity in the local standard of rest) of the local interstellar medium, assuming that it follows the rotation curve from Russeil et al. (2017). We then measured the column density (assuming optical thin HI emission), integrating the HI profile on a 10 km/s velocity range (characteristic velocity dispersion for the warm neutral gas in the inner Galaxy; Marasco et al. 2017) around the $V_{LSR}$ and in a circular aperture of 720'' radius (this size was chosen to be representative of an area of constant column density and comes from a compromise between the average size of the sources, which is 9'', the spatial resolution of the survey and the size of the background structures observed there), centred on the star position. Assuming a spherical geometry, we then deduced the HI number density ($cm^{-3}$), which was expected to be an upper limit for the ambient density. We find a mean density of $\langle n_{up} \rangle = 32.2 cm^{-3}$, which is in agreement with the average values found by Peri et al. (2012), Peri et al. (2015), and Kobulnicky et al. (2018) (a slightly larger mean value of 102 $cm^{-3}$ is reported in Kobulnicky et al. 2019).

We also estimated the ambient density using the HI4PI total column density map (HI4PI Collaboration 2016), from which, following Ochsendorf et al. (2014), we evaluated the average HI density along the line of sight by dividing the total column density at the position of the star by its distance. This gives us a lower estimate for the local interstellar medium ($\langle n_{low} \rangle = 1.9 cm^{-3}$).

Figure 2 shows that, if the two methods give different values they are, nevertheless, quite well correlated. Both methods assume a homogeneous medium but correspond to different scaling; the correlation can then suggest that most stars are in a similar environment, although density inhomogeneities are difficult to evaluate. Five stars, with $n_{low}$ larger than 3 $cm^{-3}$, strongly depart from the tendency. These stars are among those at the smallest Galactic longitudes (between 27° and 35°) and are less than 2 kpc away. The two stars with the highest density are also the closest ones (less than 610 pc). This shows that our approaches have significant limitations for stars that are closer than 2 kpc. Indeed, $n_{up}$ is badly estimated because the relationship between kinematic distance and velocity does not work for close objects, and $n_{low}$ is poorly calculated for nearby stars since the background emission is likely to dominate the HI integrated emission. In Fig. 2 we also compare the different estimates to the expected relations between density and the Galactic plane distance, z (e.g. Dickey & Lockman 1990; Loup et al. 1993). Globally, $n_{low}$ and $n_{up}$ follow the general trend expected for an exponential disk, while discrepancies can be attributed to local density increases due to star-forming regions, as illustrated by the stars in the well-known Cygnus X region and for which Emig et al. (2022) measure electron densities between $10 cm^{-3}$ and $400 cm^{-3}$.

The densities that we evaluate here using the different arc regimes presented in Fig. 2 of Henney & Arthur (2019a) (which combine the density, the ambiant velocity, and the stellar mass), suggest to us that only three stars (#5, #12, and #40) would possibly develop a radiation bow wave, instead of a wind bow shock, if $n_{up}$ is adopted.

## 6.4. Geometrical quantities

Similarly to Kobulnicky et al. (2018), we re-evaluated the stand-off distance, $R_0$, defined as the distance between the stellar candidate and the bow shock apex. It is measured directly from

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Local ISM density evaluation in the direction of the BS. The upper plot shows $n_{up}$ versus $n_{low}$. The fitted dashed line follows $n_{up} = 16.8 n_{low} + 2.3$. The middle and lower plots show $n_{low}$ and $n_{up}$ versus the distance above the Galactic plane, z, respectively. The dashed line draws the Dickey & Lockman 1990 and Loup et al. 1993 tendencies. In all the plots, the red points are the stars in the Cygnus X region ($1 \sim 79°$).</figcaption></figure>

mid-infrared images (MIPS 24 µm, WISE 22 µm, or Herschel 70 µm). Since estimating the distance to the apex is made difficult when the star is off-centred relative to the bow geometry and/or when the bow morphology is irregular, we estimate an uncertainty of 20% on the measurement of $R_0$. In our sample of 47 sources, we find $R_0$ in the range 0.01–1.5 pc. We also estimate the distance from the stellar position to the bow shock at 90° ($R_{90}$) of the apex direction. We find a mean $R_0 / R_{90}$ ratio of 0.57, in good agreement with values found by Meyer et al. (2014) for modelled bow shocks around massive runaway stars.

In addition, we measured the position angle (PA) difference between the SBN morphological axis and the Vt direction ($\Delta_{PA}$). We find 9 stars (of which 7 are runaways) with Vt well aligned with the SBN ($\Delta_{PA} < 45°$), 6 stars (of which 3 are runaways) with Vt pointing in the SBN opposite direction ($\Delta_{PA} = 180° \pm 45°$), and 21 stars (of which 10 are runaways) in the intermediate $\Delta_{PA}$. So we have a majority of non-aligned cases. This is in agreement with Kobulnicky & Chick (2022) who, analysing $\Delta_{PA}$ for a larger sample (267 SBNs), identify two populations: a small one (31 % of their sample) of highly aligned (Gaussian width of 25°) and the other one with a majority of random non-aligned cases. They
<figure><img src="image_4.png" /><figcaption>Fig. 3. Log-log plots of the dependencies of R₀ relative to Vt (upper), mass-loss rate (middle), and n<sub>low</sub> density (bottom). The dashed lines are the fitted tendencies.</figcaption></figure>

conclude that the local environment factors rather than the stellar motion dominate the formation of BSs.

On the other hand, from the theoretical expression of the standoff distance (recalled in Eq. (1)) it is expected that R₀ ∝ Vt⁻¹, R₀ ∝ M¹/², and R₀ ∝ n⁻¹/². Plotting these quantities (Fig. 3), we find a slope (in log scale) of 0.18 (correlation coefficient r = 0.78) and -0.69 (correlation coefficient r = -0.57) for R₀ versus M and n<sub>low</sub>, respectively. Obviously, the R₀ versus n<sub>low</sub> correlation is poor and so not significant. The discrepancy with the theoretical expression is most important for R₀ versus Vt, as we do not find the expected decreasing behaviour. This seems to suggest that, in most cases, the relative velocity of the ambient ISM may come more from the motion of the ISM itself than from the velocity of the star in the ISM. However, the expected velocity trend may also simply be blurred because we compared objects that have different mass-loss values and that are located in different density environments. So M seems to be the main parameter that correlates with R₀, as already supported by Kobulnicky et al. (2018, 2019).

## 7. Discussion

The relative velocity between the central star and the interstellar medium is often invoked to explain the creation of an observable bow shock. As underlined by, for example, Kobulnicky & Chick (2022), for this, it is expected that this relative velocity is provided either by the supersonic motion of the star with respect to the ISM (called a runaway star) and/or a bulk ISM flow due to photoevaporative or champagne flows in H II regions. In this context, the theoretical momentum-driven SBN relation leads to the

following expression for the standoff distance (e.g. Kobulnicky et al. 2010):

$$R_0 = \sqrt{\frac{V_w \dot{M}}{4\pi\rho_a V_a^2}}, \quad (1)$$

where V<sub>w</sub> and M are the stellar wind velocity and mass-loss rate, respectively, and V<sub>a</sub> and ρ<sub>a</sub> are the ambient relative velocity flow and local ISM density, respectively (ρ<sub>a</sub> = μn<sub>a</sub>, with the mean ISM gas mass per H atom μ = 2.36 × 10⁻²⁴ g, Povich et al. 2008). The adopted terminal wind velocities, V<sub>w</sub>, were from Martins et al. (2005a), Prinja & Massa (1998), and Lamers et al. (1995) for O-, B-, and A-type stars, respectively.

Adopting the established values for the mass-loss rate, wind terminal velocity, and standoff distance, this relation can be used in different ways. For stars for which we are not able to determine Vt, the relation is used to determine it from n<sub>low</sub> and n<sub>up</sub>. For stars with Vt aligned with the SBN, then the transverse component of the ambient ISM velocity flow, V<sub>a</sub>, corresponds to Vt, and the relation can be used to re-evaluate the ambient density. For stars with Vt not aligned with the SBN apex direction, Vt is no longer representative of the relative ambient medium velocity. We then used the relation to evaluate V<sub>a</sub> and subsequently the expected ISM flow velocity, V<sub>flow</sub> (with V<sub>flow</sub> = V<sub>a</sub> - Vt × cos (Δ<sub>PA</sub>)). All these quantities are summarised in Table A.5.

However, other processes can also lead to the formation of bows. For example, Henney & Arthur (2019a) recall that, whereas bow shocks are usually attributed to the star's wind - ISM interaction (wind-supported bow shocks, WBSs) there are regimes where the bows can be supported by the radiation (radiation-supported bow waves, RBWs, and radiation-supported bow shocks, RBSs). In particular, Henney & Arthur (2019b) underline that WBSs dominate for fast-moving stars and for low-density environments, RBWs are important for B stars and weak-wind O stars in moderate density environments (n > 100 cm⁻³), and RBSs apply to slow-moving O stars in dense environments. In bow waves, the distance (R<sub>BW</sub>) where dust grains are stopped due to radiation pressure should exceed the bow-shock standoff distance, R₀. However, Ochsendorf et al. (2014) showed that for O stars earlier than O6 and early type-B supergiants, it is expected that R<sub>BW</sub> < R₀, while for dwarf stars with spectral type later than O6, R<sub>BW</sub> > R₀. In addition, they showed that creating dust waves is most efficient around relatively slowly moving stars. For example, Ochsendorf et al. (2014) argued that the arc structure engulfing σ Ori AB (with a space velocity of only 15.1 km/s) is a dust wave in which the radiation pressure has stalled the dust carried along by the IC 434 photoevaporative flow of ionised gas coming from the dark cloud L1630. In parallel, Ngoumou et al. (2013) observed an arc-shaped feature around a high-velocity (100 km/s) B star near the stellar cluster Trumpler 14. They noted that the star is not centred with respect to the bow, but showed that it is at the edge of a molecular clump. They then argued that the arc is a bow shock produced by the stellar wind while the star is grazing along the surface of the clump, showing, from simulation, that the non-axisymmetry of the feature is consistent with the presence of a density gradient perpendicular to the bow shock. More generally, Wilkin (2000) theoretically showed that flows at oblique angles relative to a moving star can produce asymmetric bow-shock nebulae.

In addition, as underlined by Kobulnicky & Chick (2022), since Vt is calculated assuming a circular orbit model, it does not take into account circular velocity departures due to the cloud-to-cloud velocity dispersion of the ISM and streaming motion along
the arms. Indeed, arm perturbation can introduce radial and azimuthal gas streaming motions (e.g. Roberts 1969; Englmaier 2000; Ramón-Fox & Bonnell 2018) which translate into average peculiar motions towards the Galaxy centre and opposite to Galactic rotation. Such velocity departures are well known for the Perseus arm for Galactic longitudes between 90° and 160° (e.g. Roberts 1972; Brand & Blitz 1993). Typical velocity departures are about 10 km/s (e.g. Anderson et al. 2012; Wienen et al. 2015), but can reach up to 40 km/s (Brand & Blitz 1993). They can explain some of the cases of misalignment between the SBNs and Vt.

In our sample, for most of the stars, their transverse velocity is not aligned with their SBN axis. Therefore, to better understand the formation of SBNs, in the following, we study the properties of SBNs with respect to their environmental conditions.

## 7.1. Bow shocks in the Cygnus region

Eleven OB stars and their associated bow shocks in Cygnus X were studied by Kobulnicky et al. (2010). In this context, twelve stars (#31 to #42) of our sample can be associated with the Cygnus star-forming region, of which three (#31, #39, and #40) have already been studied by Kobulnicky et al. (2010). They classified the SBNs around these stars as ambiguous, probable, and doubtful for #31, #39, and #40, respectively, and they showed that some SBNs can be an interaction or illumination of a nearby cloud.

In our sample, the stars have distances between 1.4 and 1.9 kpc (except star #31, which has a slightly larger distance of 2.15 kpc), in agreement with the distances of the Cygnus OB associations, especially Cyg OB2 (Berlanas et al. 2019 suggest that the main Cyg OB2 association is located at a distance of 1.76 kpc, with a foreground population located around 1.35 kpc), which is the most distinct and most highly concentrated of them (Quintana & Wright 2021).

Quintana & Wright (2021) show that the Cygnus OB associations do not have the kinematic coherence expected for OB associations, leading them to identify six new OB associations (Quintana & Wright 2022), among which the C, E (Cyg OB2), and D associations (see Fig. 1 in Quintana & Wright 2022) show clearly different motions. Indeed, associations D and E move in the direction of increasing Galactic coordinates, while C moves in the opposite direction. In this context, we note that three stars (#37, #40, and #42) are moving in the direction of increasing Galactic coordinates, while six stars (#33, #34, #35, #36, #39, and #41) are moving in the opposite direction, consistent with these stellar association motions.

As already noted by Kobulnicky et al. (2010), the orientation of the SBNs and of Vt appear random and have no preferential alignment with respect to Cyg OB2, as expected if the stars are runaways originating from this association or if Cyg OB2 is at the origin of the local ISM pressure. Except for star #31, which is a runaway, the others have Vt between 3.4 and 13.4 km/s, more representative of the motion of the stars within their parent association (Ward & Kruijssen 2018 measure stellar velocity dispersion in associations between 3 and 13 km/s).

Runaway star #31 has its Vt vector pointing back to Cyg OB2. Taking into account its velocity and its angular distance to Cygnus OB2 (~3.2°), we estimate a travel time to its current position of 2.4 Myr, in agreement with the typical age of ~2.5 Myr for Cygnus OB2 (Negueruela et al. 2008). It is, therefore, very likely that this star originates from Cygnus OB2.

Because none of the other stars are runaways (although stars #33 and #36 seem to be runaways when considering their 3D velocity), the bows should be formed due to a bulk ISM flow, as is also supported by the non-isolated status assigned for nine of them by Jayasinghe et al. (2019). The ISM in the Cygnus region is very filamentary and Emig et al. (2022) measure electron densities of the filamentary structures between $10 \text{ cm}^{-3}$ and $400 \text{ cm}^{-3}$ (with a median value of $35 \text{ cm}^{-3}$). For three of the four aligned objects (#31, #39, and #40), the density evaluated from Eq. 1 is in this range, while for #41 (a B6.5I) a much lower density ($4 \times 10^{-3} \text{ cm}^{-3}$) is required to account for the measured $R_0$. However, #41 is quite unusual, as its SBN is better seen at Herschel-70 µm and WISE-12 µm than at WISE-22 µm which, in addition to the extremely low density estimated, means its nature is still to be clarified. For the two objects (#32 and #38) with no Vt, the estimated velocity is between 2 and 40 km/s. However, for #38 the SBN seems doubtful as the star is not centred on the arc and the bow follows a larger-scale feature.

Lastly, concerning the six non-aligned objects (stars #33 to #37 and #42), for four of them (#34 to #37), depending on the density, an estimated $V_{flow}$ in agreement with the typical ISM velocity (which is between ~6 km/s for typical photoevaporative outflow velocity, e.g. McLeod et al. 2015, and 30 km/s for gas velocity in the Champagne phase H II region, e.g. Tenorio-Tagle 1979) is found. Such photoevaporative flows could be invoked to explain the SBNs for stars #34 to #36, as they are less than 4' (~1.6 pc) away from a bright rim and in the direction in which their bow points. For star #33, the estimated $V_{flow}$ is very large (45 or 153 km/s for $n_{up}$ and $n_{low}$, respectively) and no origin for this flow is identifiable. One could assume that, as this star is a runaway, thanks to its radial velocity component, the direction of the bow could be due to density gradient.

## 7.2. Bow shocks in H II regions

In RCW 49 and M17, Povich et al. (2008) detected SBNs pointing towards the exciting clusters around OB stars located in the PDR and at the periphery of the cluster, providing very obvious examples of SBNs in H II regions. In addition, 24 µm emission is frequently observed within H II regions and near the exciting stars (e.g. Deharveng et al. 2010), occasionally exhibiting an arcuate shape, as seen in RCW 82 (Pomarès et al. 2009) and RCW 120 (Zavagno et al. 2007).

Five stars in our sample are associated with H II regions. Stars #6 and #44 are found inside the optical H II regions Sh2-66 and Sh2-135, respectively. While stars #12 and #29 are situated on the edge of radio H II regions catalogued as GB6J1910+0814 and GB6J1952+2610 (also identified as the IR bubble N131 by Churchwell et al. 2006), respectively, by Gregory et al. (1996), star #23 appears to be located at the edge of the optical H II region Sh2-86.

Star #44, at a projected distance of 2.5 pc from the main pillar's tip (Fig. 4 – upper panel) is the main ionising source of the H II region Sh2-135 (Lahulla 1985) and shows its SBN facing a very bright rim with pillars. The density calculated from Eq. 1 (Table A.5) is around $283 \text{ cm}^{-3}$ larger than $n_{up}$. However, from the radio flux (given by Anderson et al. 2016), we estimate (using Eq. 4) in Das et al. 2018 an electron density of ~39 $\text{cm}^{-3}$. Adopting this last value, a $V_{flow}$ of 16 km/s is then required to find the measured $R_0$. Since Pismis et al. (1986) reported that Sh2-135 is a blister-like H II region with the ionising gas flowing away from the CO cloud with a line of sight velocity of 4 km/s, the SBN could be produced by a photoevaporative flow coming from the rim.
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** IPHAS-Hα (left) and MIPS 24 µm (right) images of Sh2-135 (upper panel) and Sh2-66 (lower panel). In the upper panel, the arrow is the transverse velocity vector (arbitrary length) for star #44. In the lower panel, star #6 is indicated by the cross. The images are oriented according to Galactic coordinates.</figcaption></figure>

Star #6 is located in Sh2-66 and has its SBN facing a dark cloud (Fig. 4 – lower panel). Because of a bad RUWE, Vt is not estimated and the parallactic distance (8.6 kpc) is very uncertain. This distance is very different from the spectrophotometric distance (3.6 kpc, Forbes 1989) and favours the far kinematic distance (near and far distances are 2.9 or 11.4 kpc, respectively, calculated from the Fich et al. 1990 Hα LSR velocity of 46.8 km/s for the H II region). The star is at the centre of the region and is probably the main exciting source, while its projected distance to the dark cloud is ~4.2 pc (for a distance of 8.6 kpc). From the radio flux (given by Anderson et al. 2015), we estimate an electron density of ~25.5 cm⁻³, in agreement with nₐp. Due to its similar velocity (V<sub>LSR</sub> = 44.9 km/s, Eden et al. 2012), the dark cloud (GRS G030.44 +00.44) can be linked to the H II region. The PDR and the dark cloud are barely distinguishable from the background in the near-IR and mid-IR wavelengths. This suggests that the H II region evolutionary state is rather advanced. The calculated Vₐ ~ 48 km/s is quite large to be due to a photoevaporative flow alone and so better suggests the star to be a runaway.

The SBN of star #12 looks more like an enhancement in the radio H II region GB6 B1907+0809, classified as a diffuse H II region by Lockman et al. (1996). With an LSR velocity of 62 km/s (Rigby et al. 2019), the region's near and far kinematic distances are 4.1 or 8.2 kpc, respectively. The distance of #12 is ~555 pc, making it a foreground star. In addition, it is an A5Ib star that cannot ionise the H II region.

Star #29 appears located on the edge of the IR bubble N131 (Churchwell et al. 2006). The bubble morphology suggests that

this could be the PDR of an H II region that is not optically visible (Fig. C.2) and whose ionising source is unknown. According to Hou & Gao 2014 and Zhang et al. 2016, N131 is approximately 8.5 kpc away, while #29 is only 4 kpc away. The large Vt (55.3 km/s) suggests that the star is a runaway and, because the SBN and the Vt are well aligned and point outwards from the centre of the bubble, one can suspect at first glance that it has travelled from the centre of the bubble (at 5.7'). But, in this context, its travel time would be too short (0.25 Myr) with respect to the typical age for such a star (between 3.5 and 13 Myr, McEvoy et al. 2015). One can conclude that the star and its SBN are probably seen by chance in the direction of the edge of N131.

Star #23 is on the border of the optical H II region Sh2-86 (Fig. C.3). The SBN points towards the exciting cluster NGC 6823 and the distance of #23 is in agreement with that of NGC 6823 (2 kpc, e.g. Riaz et al. 2012). Located in an optically dark filament identified as a molecular filament by Kohno et al. (2022), we favour nₐp and then estimate a low velocity (Vt ~ 7.7 km/s) for the star. However, the nature of the SBN must be made clear, as it is straight rather than arched, and seems to be connected to a larger scale extension.

### 7.3. Bow shock in the cluster Berkeley 97

Star #46 (HD 240016) is observed in the direction of the open cluster Berkeley 97. The distance for Berkeley 97 is 1.8 kpc (Tadross 2008), 2.5 kpc (Dias et al. 2021), or 3.2 kpc (Poggio et al. 2021). Covered only by WISE (low spatial resolution), about 16 stars belonging to Berkeley 97 (Hunt & Reffert 2023) are observed in the SBN footprint. Among them is the star HD 240015, separated by 50.6" from #46. HD 240015 has a similar magnitude to #46 and is classified as a B0.5II star by Dinçel et al. (2024). Both stars have similar distances (Bailer-Jones et al. 2021); so, as in Dinçel et al. (2024), we can assume that #46 also belongs to the cluster and then that they are the two most massive stars of the cluster. This cluster, with log(age[yr]) ~ 7.1 (Dinçel et al. 2024), does not show an H II region and seems to be in a cleared environment, suggesting that n<sub>low</sub> is to be favoured. This is in agreement with recent results from Dinçel et al. (2024), who show that the cluster is surrounded by an H I void.

For #46 (poor RUWE value), Vt cannot be evaluated; nonetheless, its radial velocity (~32 km/s) indicates that it is a runaway, even though its motion may be severely influenced by binarity and internal motion in the cluster. The nature of the SBN is not clear either, because the cluster seems to move in the opposite direction and no interstellar medium is noted close to it.

In this context, and because #46 is not centred relative to the SBN, the SBN could be due either to a joint contribution from #46 and HD 240015 or due to the combined and unresolved contribution of the stars observed in the SBN footprint. Another alternative could be that the mid-IR emission comes from the wind interaction with any supernovae remnant as, given the age of the cluster, the most massive stars should have become supernovae. This scenario would be in agreement with Dinçel et al. (2024), who speculate that the H I void could underline a dissipated supernova remnant (SNR) and that some part of the dust dragged inward by the reverse shock interacts with massive stars to form the observed SBN. Higher spatial resolution analysis of the SBN is then needed to clarify its nature.

### 7.4. Bow shocks around cool stars

Scherer et al. (2020) point out that cool stars drive super-sonic winds and may show bow shocks if the star has velocity
<table><caption>**Table 1.** Possible parent cluster.</caption><thead><tr><th>Star</th><th>β (°)</th><th>T<sub>t</sub> (Myr)</th><th>Cluster name</th><th>Distance (kpc)</th><th>log(age(yr))</th><th>References</th></tr></thead><tbody><tr><td>#5</td><td>0.24</td><td></td><td>Theia 12</td><td>0.62±1.2 × 10<sup>-3</sup></td><td>6.48</td><td>1</td></tr><tr><td>#11</td><td>1.57</td><td>2.96</td><td>HSC 358</td><td>3.00±0.033</td><td>8.36</td><td>1</td></tr><tr><td>#13</td><td>2.26</td><td>1.62</td><td>MWSC 3093-3096</td><td>5.87±0.07</td><td>6.1</td><td>2</td></tr><tr><td>#14</td><td>2.36</td><td>2.11</td><td>Berkeley 45</td><td>3.28±0.53</td><td>8.60</td><td>1, 2, 3</td></tr><tr><td>#16</td><td>3.27</td><td>2.73</td><td>Kronberger 79</td><td>4.08±1.02</td><td>7.71</td><td>1, 2</td></tr><tr><td>#43</td><td>3.90</td><td>4.54</td><td>UBC 1169</td><td>4.55±0.084</td><td>8.46</td><td>1</td></tr><tr><td>#47</td><td>1.21</td><td>1.79</td><td>Teutsch 54</td><td>3.81±0.056</td><td>8.02</td><td>1</td></tr></tbody></table>

**Notes.** References are: (1) Hunt & Reffert (2023), (2) Just et al. (2023), and (3) Cantat-Gaudin et al. (2020). When there are several references, the mean value for the distance and log(age) are given.

with respect to the ambient ISM. In parallel, Cox et al. (2012) pointed out that evolved stars lose material, via a dusty wind, during their asymptotic giant branch (AGB) and supergiant phase, developing extended envelopes which can then form bow shocks, as around Betelgeuse which is the best example of such phenomenon (Mohamed et al. 2012).

Here, the spectral classification suggests that stars #1, #9, and #30 are cool stars, while all of them display a Vt that is not aligned with the SBN. For #1 and #9, the SBN is thick, with an irregular shape and a thin surrounding secondary arc (better seen at IRAC-8 µm), suggesting that it could better underline a photodissociation region. In addition, their SBNs present radio counterparts (HRDSG027.334+00.176 and HRDSG034.031-00.59 for #1 and #9, respectively) classified as H II regions by Anderson et al. (2011). From the main radio recombination line velocity measured by Anderson et al. (2011), HRDSG027.334+00.176 and HRDSG034.031-00.59 would be at the near kinematic distances of 4.6 and 3.4 kpc, respectively. With a distance of 1.06 and 1.65 kpc, respectively, #1 and #9 then seem to be foreground stars, not related to the background H II regions. For star #30, because it is not centred relative to the SBN and its Vt is not aligned with the SBN, one suspects that it is probably not the SBN driver. Based solely on an adequate direction and value of their tangential velocities, the 2MASS sources 19533053+2808202 or 19533020+2808216 would be better candidates.

### 7.5. Bow shocks of velocity aligned stars

In addition to the stars in the Cygnus region and H II regions, we have four star candidates with $\Delta_{PA} < 45^\circ$ and four others with $45^\circ < \Delta_{PA} < 52^\circ$ (see Table A.4). One can consider here that they have their Vt and their SBN aligned. These are #5, #10, #11, #13, #14, #16, #43, and #47. Among them, #13, #14, #16, #43, and #47 are runaways (Vt > 25 km/s). One notes also that, except for star #10 (n ~ 34 cm$^{-3}$), their measured parameters suggest a low density (between 0.003 and 3.4 cm$^{-3}$), in agreement with the weakly structured local environment aspect seen on the *Spitzer* 8 µm and mid-IR images, WISE 22 µm, and/or *Spitzer* 24 µm images (Fig. C.1).

There are two main scenarios to explain the presence of runaway or walkaway stars (stars with velocities of 5–30 km/s, de Mink et al. 2014): the binary supernova scenario (Blaauw 1961; Hoogerwerf et al. 2001) and the dynamical ejection from star clusters scenario (e.g. Hoogerwerf et al. 2001; Fujii & Portegies Zwart 2011). We then searched the possible origin for these stars by simply looking for clusters (with similar distances to the stars) present on the straight line drawn backwards. For that we used

the stellar cluster catalogues from Hunt & Reffert (2023), Just et al. (2023), and Cantat-Gaudin et al. (2020). Table 1 summarises our findings, where $\beta$ is the angular separation distance between the SBN and the candidate parental cluster from which the travel time (noted $T_t$) is calculated.

Star #5 has a low velocity (Vt = 8.22 km/s) and is identified as a member of the cluster Theia 12 (Kounkel & Covey 2019) by Hunt & Reffert (2023). It is located at the periphery of the cluster, which has a radius of 2.5° (Hunt & Reffert 2023). Due to the large extension of the cluster, it is difficult to know if the star is being ejected from the cluster or if its velocity is only representative of its motion in the cluster.

For most of the other stars we can find a possible parental cluster. Assuming that the star has been ejected early, the flight time, $T_t$, is expected to represent the star's age. We find that $T_t$ is between 1.62 and 4.54 Myr, in agreement with the typical age of O- and B-type stars (e.g. Weidner & Vink 2010; Daszyńska-Daszkiewicz & Miszuda 2019). For star #13, two clusters, MWSC 3093 and MWSC 3096 (Kharchenko et al. 2016), in the extreme star-forming complex W51, could be the ones it originates from. For this star, the cluster age and $T_t$ are in agreement with the dynamical ejection scenario. For stars #11, #14, #16, #43, and #47 the cluster age is roughly two orders of magnitude larger than $T_t$. These stars could have been ejected later after their formation via the binary supernova process, since the age of the clusters suggests that the most massive stars should have already exploded into supernovae (e.g. Crowther 2012, Eldridge & Stanway 2022). Finally, for star #10, no parental cluster candidate is found. This star is moving in the direction of the Galactic plane and is located in an elongated structure (possibly a cloud) visible from 8 µm to 250 µm. This could explain the higher density of its surrounding ISM, which would favour SBN formation.

### 7.6. Other bow shocks

There remain 18 stars, of which five have no calculated Vt, while the others have $\Delta_{PA} > 52^\circ$.

Among them, stars #3, #4, and #7 have no spectral type, and then no characteristic can be established.

For the five stars with no Vt (#18, #25, #26, #28, and #45) we estimate the expected velocity needed to produce the SBN. Depending on the density, these stars can be runaways or not, except for #45 which, whatever the density, has a low velocity and #26, which shows an unrealistic velocity (larger than 200 km/s). For #18, because of its large Galactic plane distance of 168.9 pc, we would favour $n_{low}$, hence a runaway status with a velocity of 86.5 km/s.
For the other thirteen stars with misaligned Vt (#2, #3, #4, #7, #8, #15, #17, #19, #20, #21, #22, #24, and #27), eight of them (#2, #4, #15, #17, #19, #20, #21, and #22) are runaways (with Vt between 31.3 and 80.7 km/s), for which it is therefore difficult to understand the misalignment of their bow shock relative to their transverse velocity direction. For #20, the suggested $V_{\text{flow}}$ is unrealistic (>100 km/s), but its SBN clearly points towards the star-forming region GAL54.4+01 (Day et al. 1972) (which is at a similar distance of 7.3±1 kpc, Anderson et al. 2012), while its Vt direction suggests it could have been ejected from it.

Five stars (#24, #25, #26, #27, and #28) are observed in the Vulpecula OB1 association region, which is a region of active star formation (Billot et al. 2010). Three of them (#26, #27, and #28) have their SBNs facing Sh2-87 (distance ~ 2.3 kpc; Xue & Wu 2008), while the others (#24 and #25) have their SBNs pointing towards pillars which themselves point in the opposite direction towards stars belonging to the Vulpecula OB1 association (distance ~ 1.6 kpc; Mel'nik & Dambis 2017). But #25, #26, and #28 have a distance placing them farther away, and so not related to these structures. Conversely, stars #24 and #27, due to their distance (~2 kpc) and their Vt < 20 km/s, may belong to the Vulpecula OB1 association.

Six stars exhibit a SBN with a remarkable morphology (see Fig. C.1): #2 shows two cirrus-like arcs on either side of the star, #7 has a quasi-circular aspect (it is also classified as a bubble MWP1G031117+007214 by Simpson et al. 2012), #17 is surrounded by irregular structures of the ISM, #22 has an isolated long, thin and distorted SBN (it shows the largest $R_0$ ~ 1.5 pc), #26 has four small radial finger-shaped extensions which extend from the SBN, one of which passes through the star, and #28 has an SBN drawing three parallel arcs. Several explanations can be invoked to explain these shapes: the SBN may be formed by a density gradient perpendicular to the Vt (e.g. Ngoumou et al. 2013), or it is a simple enhancement of the local interstellar medium, or it is due to previous stellar mass ejections (#26 and #28 are long-period variable candidates from Lebzelter et al. 2023), or it is caused by dust grains interacting with the interstellar magnetic field while disturbed by the stellar wind collision (e.g. Katushkina et al. 2018).

Three stars are found in the direction of SNRs: #17, #21, and #22 fall in the direction of the SNRs G53.4+0.03 (Anderson et al. 2017), G55.0+0.3 (Green et al. 2022), and G56.56-0.75 (Anderson et al. 2017), respectively. The SNR-SBN link is difficult to establish because the SNR's distance is not well constrained. In addition, their SBNs do not point towards the SNR's centre, making their link uncertain. Furthermore, their stellar Vt does not point away from the SNR's centre, indicating that these stars were not expelled during the supernova explosion. In this context, at least for #17, the SBN could simply be an ISM substructure shaped by the SNR.

## 8. Conclusions and perspectives

We led a 400–800 nm spectral survey of a sample of 47 candidate bow-shock stars. We complemented these observations with Gaia ERD3 data to find their transverse velocity relative to the interstellar medium. We find that, of the 47 candidates, there are 3 unclassifiable stars, 3 cool stars, 2 A-type stars, 10 O stars, and 29 B (mainly giant and supergiant) stars. We find that 17 stars (among the 37 with determined Vt) are runaways, of which 15 are isolated, but only 7 have their transverse velocity aligned with the SBN axis. This confirms that runaway is not the only origin for SBN formation.

We estimated, as independently as possible, the different quantities involved in the momentum balance equation which characterises the SBN. The uncertainties are large and only the relationship between $R_0$ and $\dot{M}$ seems robust, confirming that the stellar wind is the most important parameter. The misalignment between the SBN axis and the transverse velocity direction of the star is observed even for runaway stars, suggesting that the environment (ambient density inhomogeneity, ISM flows, etc.) has a strong influence on the appearance and asymmetry of the bow shocks.

In this context, we note the diversity of environments where SBNs are observed: stellar associations, clusters, H II regions, and possibly SNRs. For stars with transverse velocity aligned with the BS, whether or not they are runaways, we were able, for most of them, to identify a possible parental cluster and we confirm that these stars are mainly located in a low-density medium. As already observed, numerous SBNs are present in the Cygnus OB stellar association, suggesting that the ISM conditions in such regions favour SBN observability. Few SBNs are observed in H II regions. For one of them (Sh2-135), the SBN could be explained by a ~16 km/s photoevaporated flow coming from the H II region molecular cloud's interface. It appears then that it is necessary to specify the conditions for a flow to form an SBN: flow velocity and density, star distance from the interface, and so on. To validate the photoevaporative flow mechanism in SBN creation, statistical analysis of SBNs in H II regions and direct flow velocity determination are necessary. A more general approach would be to quantify and identify what causes the emergence of SBNs in star-forming regions: turbulence, magnetic field, density fluctuations, medium nature, and so on.

Finally, if few SBNs are observed in the direction of cool stars they seem to not be related.

A direct improvement of our study would be to have higher mid-resolution images for objects with only WISE-22 µm observations and to update and specify the transverse velocity value and direction through future GAIA releases.

**Acknowledgements.** The authors thank the anonymous refree for its his/her insightful remarks and comments. The authors acknowledge Interstellar Institute's program "I16" and the Paris-Saclay University's Institut Pascal for hosting discussions that nourished the development of the ideas behind this work. In particular, the authors thank P. Lesaffre for discussions during I16 and M. Dennefeld for his contribution to the observations. Based on observations made at Observatoire de Haute Provence (CNRS), France, with MISTRAL on the T193 telescope. This research also has made use of the MISTRAL database, operated at CeSAM (LAM), Marseille, France. AZ thanks the support of the Institut Universitaire de France. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France and "Aladin sky atlas" developed at CDS, Strasbourg Observatory, France.

## References

Anderson, L. D., Bania, T. M., Balser, D. S., & Rood, R. T. 2011, ApJS, 194, 32

Anderson, L. D., Bania, T. M., Balser, D. S., & Rood, R. T. 2012, ApJ, 754, 62

Anderson, L. D., Armentrout, W. P., Johnstone, B. M., et al. 2015, ApJS, 221, 26

Anderson, L. D., Armentrout, W. P., Johnstone, B. M., et al. 2016, VizieR Online Data Catalog: Radio observations of Galactic WISE H II regions (Anderson+, 2015), VizieR On-line Data Catalog: J/ApJS/221/26. Originally published in: 2015ApJS...221...26A

Anderson, L. D., Wang, Y., Bihr, S., et al. 2017, A&A, 605, A58

Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147

Berlanas, S. R., Herrero, A., Comerón, F., et al. 2018, A&A, 612, A50

Berlanas, S. R., Wright, N. J., Herrero, A., Drew, J. E., & Lennon, D. J. 2019, MNRAS, 484, 1838

Berlanas, S. R., Herrero, A., Comerón, F., et al. 2020, A&A, 642, A168

Billot, N., Noriega-Crespo, A., Carey, S., et al. 2010, ApJ, 712, 797

Blaauw, A. 1961, Bull. Astron. Inst. Netherlands, 15, 265
Brand, J., & Blitz, L. 1993, A&A, 275, 67

Cantat-Gaudin, T., Anders, F., Castro-Ginard, A., et al. 2020, A&A, 640, A1

Lame, S. J., Noriega-Crespo, A., Mizuno, D. R., et al. 2009, PASP, 121, 76

Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, arXiv e-prints [arXiv:1612.05560]

Chick, W. T., Kobulnicky, H. A., Schurhammer, D. P., et al. 2020, ApJS, 251, 29

Churchwell, E., Povich, M. S., Allen, D., et al. 2006, ApJ, 649, 759

Churchwell, E., Babler, B. L., Meade, M. R., et al. 2009, PASP, 121, 213

Comerón, F., & Pasquali, A. 2012, A&A, 543, A101

Conti, P. S. 1974, ApJ, 187, 539

Cox, N. L. J., Kerschbaum, F., van Marle, A. J., et al. 2012, A&A, 537, A35

Crowther, P. 2012, Astron. Geophys., 53, 430

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246

Cutri, R. M., Wright, E. L., Conrow, T., et al. 2013, Explanatory Supplement to the AllWISE Data Release Products

Das, S. R., Tej, A., Vig, S., et al. 2018, A&A, 612, A36

Daszyńska-Daszkiewicz, J., & Miszuda, A. 2019, ApJ, 886, 35

Day, G. A., Caswell, J. L., & Cooke, D. J. 1972, Aust. J. Phys. Astrophys. Suppl., 25, 1

Deharveng, L., Schuller, F., Anderson, L. D., et al. 2010, A&A, 523, A6

de Jager, C., & Nieuwenhuijzen, H. 1987, A&A, 177, 217

de Mink, S. E., Sana, H., Langer, N., Izzard, R. G., & Schneider, F. R. N. 2014, ApJ, 782, 7

Deng, D., Sun, Y., Wang, T., Wang, Y., & Jiang, B. 2022, ApJ, 935, 175

Dias, W. S., Monteiro, H., Moitinho, A., et al. 2021, MNRAS, 504, 356

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215

Didelon, P. 1982, A&AS, 50, 199

Dinçel, B., Sheth, S., Specht, L., et al. 2024, A&A, 691, A63

Drew, J. E., Greimel, R., Irwin, M. J., et al. 2005, MNRAS, 362, 753

Eden, D. J., Moore, T. J. T., Plume, R., & Morgan, L. K. 2012, MNRAS, 422, 3178

Eldridge, J. J., & Stanway, E. R. 2022, ARA&A, 60, 455

Emig, K. L., White, G. J., Salas, P., et al. 2022, A&A, 664, A88

Englmaier, P. 2000, Rev. Mod. Astron., 13, 97

Evans, C. J., Howarth, I. D., Irwin, M. J., Burnley, A. W., & Harries, T. J. 2004, MNRAS, 353, 601

Fabricius, C., Luri, X., Arenou, F., et al. 2021, A&A, 649, A5

Fich, M., Treffers, R. R., & Dahl, G. P. 1990, AJ, 99, 622

Forbes, D. 1989, A&AS, 77, 439

Fujii, M. S., & Portegies Zwart, S. 2011, Science, 334, 1380

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Georgelin, Y. M., Georgelin, Y. P., & Roux, S. 1973, A&A, 25, 337

Gray, R. O., & Corbally, C. J. 2009, Stellar Spectral Classification (Princeton University Press)

Green, S., Mackey, J., Kavanagh, P., et al. 2022, A&A, 665, A35

Gregory, P. C., Scott, W. K., Douglas, K., & Condon, J. J. 1996, ApJS, 103, 427

Groot, P. J., Verbeek, K., Greimel, R., et al. 2009, MNRAS, 399, 323

Gutermuth, R. A., & Heyer, M. 2015, AJ, 149, 64

Gvaramadze, V. V., Kniazev, A. Y., Kroupa, P., & Oh, S. 2011, A&A, 535, A29

Hardorp, J., Theile, I., & Voigt, H. H. 1964, Hamburger Sternw. Warner & Swasey Obs., C03, 0

Henney, W. J., & Arthur, S. J. 2019a, MNRAS, 486, 3423

Henney, W. J., & Arthur, S. J. 2019b, MNRAS, 486, 4423

HI4PI Collaboration (Ben Bekhti, N., et al.) 2016, A&A, 594, A116

Hoogerwerf, R., de Bruijne, J. H. J., & de Zeeuw, P. T. 2001, A&A, 365, 49

Hou, L. G., & Gao, X. Y. 2014, MNRAS, 438, 426

Hovhannessian, R. K., & Hovhannessian, E. R. 2001, Astrophysics, 44, 454

Hunt, E. L., & Reffert, S. 2023, A&A, 673, A114

Jaschek, M., Andrillat, Y., Houziaux, L., & Jaschek, C. 1994, A&A, 282, 911

Jayasinghe, T., Dixon, D., Povich, M. S., et al. 2019, MNRAS, 488, 1141

Just, A., Piskunov, A. E., Klos, J. H., Kovaleva, D. A., & Polyachenko, E. V. 2023, A&A, 672, A187

Katushkina, O. A., Alexashov, D. B., Gvaramadze, V. V., & Izmodenov, V. V. 2018, MNRAS, 473, 1576

Kharchenko, N. V., Piskunov, A. E., Schilbach, E., Röser, S., & Scholz, R. D. 2016, A&A, 585, A101

Kobulnicky, H. A., & Chick, W. T. 2022, AJ, 164, 86

Kobulnicky, H. A., Gilbert, I. J., & Kiminki, D. C. 2010, ApJ, 710, 549

Kobulnicky, H. A., Lundquist, M. J., Bhattacharjee, A., & Kerton, C. R. 2012, AJ, 143, 71

Kobulnicky, H. A., Chick, W. T., Schurhammer, D. P., et al. 2016, ApJS, 227, 18

Kobulnicky, H. A., Chick, W. T., & Povich, M. S. 2018, ApJ, 856, 74

Kobulnicky, H. A., Chick, W. T., & Povich, M. S. 2019, AJ, 158, 73

Kohno, M., Nishimura, A., Fujita, S., et al. 2022, PASJ, 74, 24

Kounkel, M., & Covey, K. 2019, AJ, 158, 122

Lahulla, J. F. 1985, A&AS, 61, 537

Lame, M. C., Smith, R. J., Arcavi, I., et al. 2023, AJ, 166, 13

Lamers, H. J. G. L. M., Snow, T. P., & Lindholm, D. M. 1995, ApJ, 455, 269

Lebzelter, T., Mowlavi, N., Lecoeur-Taibi, I., et al. 2023, A&A, 674, A15

Lennon, D. J., Dufton, P. L., & Fitzsimmons, A. 1993, A&AS, 97, 559

Leone, F., & Lanzafame, A. C. 1998, A&A, 330, 306

Lockman, F. J., Pisano, D. J., & Howard, G. J. 1996, ApJ, 472, 173

Loup, C., Forveille, T., Omont, A., & Paul, J. F. 1993, A&AS, 99, 291

Mackey, J., Green, S., Moutzouri, M., et al. 2021, MNRAS, 504, 983

Maŕz Apellániz, J., Sota, A., Arias, J. I., et al. 2016, ApJS, 224, 4

Marasco, A., Fraternali, F., van der Hulst, J. M., & Oosterloo, T. 2017, A&A, 607, A106

Marcolino, W. L. F., Bouret, J. C., Martins, F., et al. 2009, A&A, 498, 837

Martinez, J. R., del Palacio, S., & Bosch-Ramon, V. 2023, A&A, 680, A99

Martins, F. 2018, A&A, 616, A135

Martins, F., & Palacios, A. 2021, A&A, 645, A67

Martins, F., Schaerer, D., & Hillier, D. J. 2005a, A&A, 436, 1049

Martins, L. P., González Delgado, R. M., Leitherer, C., Cerviño, M., & Hauschildt, P. 2005b, MNRAS, 358, 49

McEvoy, C. M., Dufton, P. L., Evans, C. J., et al. 2015, A&A, 575, A70

McLeod, A. F., Dale, J. E., Ginsburg, A., et al. 2015, MNRAS, 450, 1057

Mel'nik, A. M., & Dambis, A. K. 2017, MNRAS, 472, 3887

Meyer, D. M. A., Mackey, J., Langer, N., et al. 2014, MNRAS, 444, 2754

Mohamed, S., Mackey, J., & Langer, N. 2012, A&A, 541, A1

Mokiem, M. R., de Koter, A., Evans, C. J., et al. 2007, A&A, 465, 1003

Monguió, M., Greimel, R., Drew, J. E., et al. 2020, A&A, 638, A18

Neguerauela, I., Marco, A., Herrero, A., & Clark, J. S. 2008, A&A, 487, 575

Ngoumou, J., Preibisch, T., Ratzka, T., & Burkert, A. 2013, ApJ, 769, 139

Ochsendorf, B. B., Cox, N. L. J., Krijt, S., et al. 2014, A&A, 563, A65

Peppel, U. 1984, A&AS, 57, 107

Peri, C. S., Benaglia, P., Brookes, D. P., Stevens, I. R., & Isequilla, N. L. 2012, A&A, 538, A108

Peri, C. S., Benaglia, P., & Isequilla, N. L. 2015, A&A, 578, A45

Pismis, P., Hasse, I., & Moreno, M. A. 1986, Rev. Mexicana Astron. Astrofis., 13, 131

Poggio, E., Drimmel, R., Cantat-Gaudin, T., et al. 2021, A&A, 651, A104

Pomarès, M., Zavagno, A., Deharveng, L., et al. 2009, A&A, 494, 987

Povich, M. S., Benjamin, R. A., Whitney, B. A., et al. 2008, ApJ, 689, 242

Prinja, R. K., & Massa, D. L. 1998, in *Astronomical Society of the Pacific Conference Series*, 131, Properties of Hot Luminous Stars, ed. I. Howarth, 218

Puls, J., Vink, J. S., & Najarro, F. 2008, A&A Rev., 16, 209

Quintana, A. L., & Wright, N. J. 2021, MNRAS, 508, 2370

Quintana, A. L., & Wright, N. J. 2022, MNRAS, 515, 687

Raga, A. C., Cantó, J., & Noriega-Crespo, A. 2022, Rev. Mexicana Astron. Astrofis., 58, 395

Ramón-Fox, F. G., & Bonnell, I. A. 2018, MNRAS, 474, 2028

Riaz, B., Martín, E. L., Tata, R., et al. 2012, MNRAS, 419, 1887

Rigby, A. J., Moore, T. J. T., Eden, D. J., et al. 2019, A&A, 632, A58

Roberts, W. W. 1969, ApJ, 158, 123

Roberts, Jr., W. W. 1972, ApJ, 173, 259

Russeil, D. 2003, A&A, 397, 133

Russeil, D., Zavagno, A., Mège, P., et al. 2017, A&A, 601, L5

Russeil, D., Zavagno, A., Nguyen, A., et al. 2020, A&A, 642, A21

Scherer, K., Baalmann, L. R., Fichtner, H., et al. 2020, MNRAS, 493, 4172

Schmidt-Kaler, T. 1982, Landolt-Börnstein, New Series, Gr. 6, Vol. 2, Subvol. B

Schmitt, J., Adami, C., Dennefeld, M., et al. 2024, A&A, 687, A198

Siebenmorgen, R., Scicluna, P., & Krelowski, J. 2018, A&A, 620, A32

Simpson, R. J., Povich, M. S., Kendrew, S., et al. 2012, MNRAS, 424, 2442

Tadross, A. L. 2008, MNRAS, 389, 285

Tenorio-Tagle, G. 1979, A&A, 71, 59

Tetzlaff, N., Neuhäuser, R., & Hohle, M. M. 2011, MNRAS, 410, 190

Uranova, T. A. 1977, Astron. Tsirk., 946, 4

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574

Wang, S., & Chen, X. 2019, ApJ, 877, 116

Ward, J. L., & Kruijssen, J. M. D. 2018, MNRAS, 475, 5659

Weidner, C., & Vink, J. S. 2010, A&A, 524, A98

Wienen, M., Wyrowski, F., Menten, K. M., et al. 2015, A&A, 579, A91

Wilkin, F. P. 2000, ApJ, 532, 400

Willmer, C. N. A. 2018, ApJS, 236, 47

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Xue, R., & Wu, Y. 2008, ApJ, 680, 446

Zavagno, A., Pomarès, M., Deharveng, L., et al. 2007, A&A, 472, 835

Zhang, C.-P., Li, G.-X., Wyrowski, F., et al. 2016, A&A, 585, A117
## Appendix A: Tables

<table><caption>**Table A.1.** B-type star classification test.</caption><thead><tr><th>Star name / Library ident.</th><th>Spec. type.<br>(from the lib.)</th><th>EW706.5<br>Å</th><th>EW667.8<br>Å</th><th>EW587.6<br>Å</th><th>EWHα<br>Å</th><th>Spec. type<br>(from our scheme)</th></tr></thead><tbody><tr><td>HD184927 / STELIB</td><td>B2V</td><td>0.695</td><td>1.322</td><td>1.479</td><td>3.765</td><td>B0.5III</td></tr><tr><td>HD135485 / MILES</td><td>B3V</td><td>0.405</td><td>0.516</td><td>0.729</td><td>5.356</td><td>B3V</td></tr><tr><td>HD187879 / MILES</td><td>B1III</td><td>0.391</td><td>0.907</td><td>0.930</td><td>3.295</td><td>B0.5III</td></tr><tr><td>HD206165 / MILES</td><td>B2Ib</td><td>0.699</td><td>1.092</td><td>1.104</td><td>0.961</td><td>B3Ia</td></tr><tr><td>HD013267 / MILES</td><td>B5Ia</td><td>0.487</td><td>0.705</td><td>0.787</td><td>0.603</td><td>B6Ia</td></tr><tr><td>HD176301 / MILES</td><td>B7III-IV</td><td></td><td>0.296</td><td>0.357</td><td>5.431</td><td>B7Ib</td></tr><tr><td>HD174959 / MILES</td><td>B6IV</td><td>0.187</td><td>0.248</td><td>0.329</td><td>5.634</td><td>B7V</td></tr><tr><td>HD138764 / MILES</td><td>B6IV</td><td>0.090</td><td>0.200</td><td>0.341</td><td>5.945</td><td>B9V</td></tr><tr><td>HD164432 / MILES</td><td>B2IV</td><td>0.523</td><td>0.681</td><td>0.682</td><td>3.804</td><td>B3III</td></tr><tr><td>HD160762 / MILES</td><td>B3IV</td><td>0.479</td><td>0.525</td><td>0.683</td><td>4.501</td><td>B2V</td></tr><tr><td>HD180163 / MILES</td><td>B2.5IV</td><td>0.316</td><td>0.503</td><td>0.636</td><td>3.733</td><td>B4V</td></tr><tr><td>HD147550 / xshooter</td><td>B9V</td><td></td><td></td><td>0.052</td><td>8.456</td><td>B8V</td></tr><tr><td>HD34816 / xshooter</td><td>B0.5V</td><td>0.442</td><td>0.661</td><td>0.695</td><td>2.750</td><td>B2V</td></tr></tbody></table>

<table><caption>**Table A.2.** O-type star classification test.</caption><thead><tr><th>Star name / Library ident.</th><th>Spec. type.<br>(from the lib.)</th><th>EW706.5<br>Å</th><th>EW541.2<br>Å</th><th>EW468.6<br>Å</th><th>Spec. type<br>(from our scheme)</th></tr></thead><tbody><tr><td>HD057060 / MILES</td><td>O7Ia</td><td>0.379</td><td>0.619</td><td>-1.243</td><td>O7Ia</td></tr><tr><td>HD057061 / MILES</td><td>O9Ib</td><td>0.599</td><td>0.519</td><td>0.317</td><td>O9II</td></tr><tr><td>HD209975 / MILES</td><td>O8Ib</td><td>0.733</td><td>0.457</td><td>0.285</td><td>O9Ib</td></tr><tr><td>HD46223 / STELIB</td><td>O5e</td><td>0.194</td><td>0.925</td><td>0.608</td><td>O5III</td></tr><tr><td>HD41161 / MISTRAL</td><td>O8Vn*</td><td>0.406</td><td>0.470</td><td>0.719</td><td>O7.5V</td></tr><tr><td>HD216898 / MISTRAL</td><td>O9V*</td><td>0.587</td><td>0.779</td><td>0.702</td><td>O8V</td></tr><tr><td>HD34816 / xshooter</td><td>B0.5V</td><td>0.442</td><td>0.095</td><td>0.214</td><td>B0Ib</td></tr></tbody></table>

Note: *: Spectral type from Simbad⁴.
<table><caption>Table A.3. Equivalent Widths and Spectral classification of the stellar sample.</caption><thead><tr><th rowspan="2">nb.</th><th rowspan="2">MWP ident.</th><th rowspan="2">Hα</th><th rowspan="2">HeH668.0</th><th rowspan="2">He706.6</th><th rowspan="2">He587.6</th><th rowspan="2">He541.1</th><th colspan="2">Equivalent Widths (Å)</th><th rowspan="2">Hβ</th><th rowspan="2">He468.6</th><th rowspan="2">He454.2</th><th rowspan="2">He447.1</th><th rowspan="2">He438.8</th><th rowspan="2">Classif.</th><th rowspan="2">Notes</th></tr><tr><th>He501.6</th><th>He492.1</th></tr></thead><tbody><tr><td>1</td><td>2G0723335+007832</td><td>1.626±2.9 10<sup>-2</sup></td><td>0.8794±0.108</td><td>0.76±0.110</td><td>0.960±0.131</td><td></td><td></td><td></td><td>2.005±0.130</td><td></td><td></td><td></td><td></td><td>M-type</td><td>M</td></tr><tr><td>2</td><td>2G07287881+0062158</td><td>2.265±0.124</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B</td><td>Noisy spectrum; G</td></tr><tr><td>3</td><td>2G07287072+0010280</td><td>2.784±0.209</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B</td><td>Noisy spectrum; K</td></tr><tr><td>4</td><td>2G0300472+0017596</td><td>4.893±0.100</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B</td><td>CW/OU/NII; O9.5III; O</td></tr><tr><td>5</td><td>2G0302964+0057613</td><td>2.593±6.0 10<sup>-2</sup></td><td>0.639±6.7 10<sup>-2</sup></td><td>0.496±8.5 10<sup>-2</sup></td><td>0.665±6.4 10<sup>-2</sup></td><td>0.653±0.119</td><td>0.133±0.165</td><td>0.539±0.108</td><td>6.999±0.576</td><td>3.25±0.186</td><td>1.094±0.116</td><td></td><td></td><td>A5b</td><td>Noisy spectrum; K</td></tr><tr><td>6</td><td>2G0304619+0039234</td><td>1.993±7.6 10<sup>-2</sup></td><td>0.729±0.109</td><td>0.360±0.167</td><td>0.383±0.180</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1II</td><td>CW/OU/NII; O9.5III; O</td></tr><tr><td>7</td><td>2G0311175+0072601</td><td>0.619±0.148</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1II</td><td>CW/OU/NII; O9.5III; O</td></tr><tr><td>8</td><td>2G0311438+0087763</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>M-type</td><td>M</td></tr><tr><td>9</td><td>2G0340267+0045756</td><td>1.969±0.190</td><td>0.793±0.145</td><td>0.801±0.136</td><td>0.901±0.119</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>B; O</td></tr><tr><td>10</td><td>2G0372791+0077658</td><td>2.264±0.169</td><td>0.961±4.6 10<sup>-2</sup></td><td>0.726±0.206</td><td>0.823±8.9 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>B; O</td></tr><tr><td>11</td><td>2G0390920+0069154</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>No Hel line; F</td></tr><tr><td>12</td><td>2G0423947+0038527</td><td>5.315±0.208</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>No Hel line; F</td></tr><tr><td>13</td><td>2G0469130+0061119</td><td>no Halpha</td><td>1.127±0.181</td><td>0.581±0.146</td><td>0.798±0.105</td><td>0.593±0.159</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1.5I</td><td>B; O</td></tr><tr><td>14</td><td>2G0480954+0019923</td><td>2.065±0.114</td><td>0.829±9.9 10<sup>-2</sup></td><td>0.711±9.7 10<sup>-2</sup></td><td>0.557±0.195</td><td>0.752±6.9 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>OB; O</td></tr><tr><td>15</td><td>2G0508365+0040783</td><td>2.702±0.151</td><td>0.478±0.143</td><td>1.249±0.186</td><td>0.948±0.145</td><td>0.929±8.9 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>OB; O</td></tr><tr><td>16</td><td>2G0534177+0009895</td><td>1.228±0.225</td><td>1.142±8.4 10<sup>-2</sup></td><td>0.877±0.125</td><td>0.864±0.119</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>O</td></tr><tr><td>17</td><td>2G0541619+0094689</td><td>2.944±0.107</td><td>0.835±0.108</td><td>0.955±0.170</td><td>0.593±0.281</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>B; O</td></tr><tr><td>18</td><td>2G0541763+0068755</td><td>1.857±0.140</td><td>0.634±0.108</td><td>0.959±0.170</td><td>0.700±0.150</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>B0.5III</td></tr><tr><td>19</td><td>2G0543235+0094471</td><td>2.007±0.102</td><td>0.597±5.2 10<sup>-2</sup></td><td>0.891±9.1 10<sup>-2</sup></td><td>0.706±0.150</td><td></td><td>1.032±0.217</td><td></td><td>2.136±0.183</td><td>2.850±8.7 10<sup>-2</sup></td><td>-1.570±2.8 10<sup>-2</sup></td><td>-1.366±9.0 10<sup>-2</sup></td><td>0.837±0.194</td><td>B1</td></tr><tr><td>20</td><td>2G0543235+0094471</td><td>2.596±0.161</td><td>0.697±9.1 10<sup>-2</sup></td><td>0.461±0.135</td><td>0.544±8.9 10<sup>-2</sup></td><td></td><td>1.368±8.8 10<sup>-2</sup></td><td>0.232±8.2 10<sup>-2</sup></td><td>2.736±7.9 10<sup>-2</sup></td><td>1.293±0.277</td><td>0.901±0.520</td><td>0.590±0.116</td><td>0.837±0.194</td><td>B1</td></tr><tr><td>21</td><td>2G0553236+0036329</td><td>1.062±0.285</td><td>0.768*±8.1 10<sup>-2</sup></td><td>0.410±0.102</td><td>0.716±0.188</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>CW/OU/NII; B; B</td></tr><tr><td>22</td><td>2G0565716+0082862</td><td>2.01±0.167</td><td>0.951±6.3 10<sup>-2</sup></td><td>0.837±0.230</td><td>0.602±0.183</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>O; O</td></tr><tr><td>23</td><td>2G0592997+0012519</td><td>2.792±0.113</td><td>0.741±8.4 10<sup>-2</sup></td><td>0.598±0.234</td><td>0.803±7.4 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B2I</td><td>CW/OU/NII; B3; B</td></tr><tr><td>24</td><td>2G0600615+0018340</td><td>2.793±0.113</td><td>0.741±8.4 10<sup>-2</sup></td><td>0.598±0.234</td><td>0.803±7.4 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>CW/OU/NII; O8(n); O (f)</td></tr><tr><td>25</td><td>2G0601249+0011967</td><td>2.85±0.114</td><td>0.607±5.9 10<sup>-2</sup></td><td>0.435±0.189</td><td>0.575±6.6 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>CW/OU/NII; O9.5III; O</td></tr><tr><td>26</td><td>2G0605078+0028576</td><td>-0.099±0.108</td><td>1.10±8.1 10<sup>-2</sup></td><td>0.910±5.6 10<sup>-2</sup></td><td>0.704±0.197</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>Lum; Class from Hβ; Oe; B</td></tr><tr><td>27</td><td>2G0609757+0034412</td><td>2.542±0.220</td><td>0.595±7.5 10<sup>-2</sup></td><td>0.410±0.100</td><td>0.518±7.4 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>OB; B</td></tr><tr><td>28</td><td>2G0619180+0013914</td><td>-5.17±0.369</td><td>0.742±9.8 10<sup>-2</sup></td><td>0.451±0.366</td><td>0.814±0.179</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>O</td></tr><tr><td>29</td><td>2G0629904+0043101</td><td>2.729±0.266</td><td>0.949±0.147</td><td>0.5749±0.264</td><td>0.281±0.265</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1II</td><td>K</td></tr><tr><td>30</td><td>2G0647586+0028772</td><td>1.162±8.2 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1III; O</td><td>K</td></tr><tr><td>31</td><td>2G0771850+0060886</td><td>3.507±0.128</td><td>0.511±7.1 10<sup>-2</sup></td><td>0.629±0.123</td><td>1.152±0.149</td><td></td><td>1.430±0.269</td><td></td><td>0.884±0.198</td><td></td><td></td><td></td><td></td><td>B1V; B2V; O; O</td><td>B1V; B2V; O; O</td></tr><tr><td>32</td><td>2G0770694+00195568</td><td>2.790±0.122</td><td>0.720±6.3 10<sup>-2</sup></td><td>0.496±0.246</td><td>0.635±7.8 10<sup>-2</sup></td><td></td><td>0.252±0.111</td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>CW/OU/NII; SHII; B</td></tr><tr><td>33</td><td>2G0773612+0011690</td><td>1.310±8.0 10<sup>-2</sup></td><td>1.015±4.5 10<sup>-2</sup></td><td>0.698±0.172</td><td>0.865±6.4 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1.5V; B</td><td>B1V; B</td></tr><tr><td>34</td><td>2G0777958+0097869</td><td>1.578±4.1 10<sup>-2</sup></td><td>0.915±3.8 10<sup>-2</sup></td><td>0.686±0.217</td><td>1.064±6.1 10<sup>-2</sup></td><td></td><td>0.115±0.107</td><td></td><td>0.519±0.141</td><td></td><td></td><td></td><td></td><td>B1</td><td>CW/OU/NII; O7; O</td></tr><tr><td>35</td><td>2G0782866+0077837</td><td>2.756±5.4 10<sup>-2</sup></td><td>0.733±7.2 10<sup>-2</sup></td><td>0.582±8.7 10<sup>-2</sup></td><td>0.727±6.7 10<sup>-2</sup></td><td></td><td>0.102±0.253</td><td>0.643±0.267</td><td>0.704±0.187</td><td>0.417±0.348</td><td></td><td></td><td></td><td>B1II; O9.5Vn; O</td></tr><tr><td>36</td><td>2G0784216+0081171</td><td>2.748±0.135</td><td>0.823±5.5 10<sup>-2</sup></td><td>0.636±0.107</td><td>0.761±7.0 10<sup>-2</sup></td><td></td><td>0.238±0.233</td><td></td><td>0.238±0.233</td><td></td><td></td><td></td><td></td><td>B1II; O9.5Vn; O</td></tr><tr><td>37</td><td>2G079207+00120627</td><td>1.805±4.2 10<sup>-2</sup></td><td>0.418 ± 2.110<sup>-2</sup></td><td>0.625±6.9 10<sup>-2</sup></td><td>0.790±6.4 10<sup>-2</sup></td><td></td><td>0.178±0.175</td><td>0.811±0.173</td><td>0.577±0.150</td><td></td><td></td><td></td><td></td><td>B2I</td><td>B1V; B2V; O; O</td></tr><tr><td>38</td><td>2G0795135+0008924</td><td>2.119±8.1 10<sup>-2</sup></td><td>0.177* ± 2.410<sup>-2</sup></td><td>0.673±0.104</td><td>0.757±0.117</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B3I</td><td>B1V; B2II</td></tr><tr><td>39</td><td>2G0798222+0009549</td><td>2.777±0.143</td><td>0.966±9.3 10<sup>-2</sup></td><td>0.625±0.155</td><td>0.586±0.187</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B3III</td><td>B3III</td></tr><tr><td>40</td><td>2G0802624+00191414</td><td>1.785±3.3 10<sup>-2</sup></td><td>0.512±4.9 10<sup>-2</sup></td><td>0.503±5.4 10<sup>-2</sup></td><td>0.757±2.2 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B3III</td><td>B3III; O</td></tr><tr><td>41</td><td>2G0809410+0001786</td><td>2.726±7.6 10<sup>-2</sup></td><td>0.680±0.102</td><td>0.339±9.9 10<sup>-2</sup></td><td>0.372±2.48</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B6.5I</td><td>O</td></tr><tr><td>42</td><td>2G0814164+0037181</td><td>non gauss. em.</td><td>0.691±0.106</td><td>0.650±0.147</td><td>0.483±7.6 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B4I</td><td>O; O</td></tr><tr><td>43</td><td>2G1043461+0022039</td><td>2.898±5.2 10<sup>-2</sup></td><td>0.482±9.6 10<sup>-2</sup></td><td>1.090±0.166</td><td>0.244±0.109</td><td></td><td>1.078±0.138</td><td></td><td>0.449±0.142</td><td></td><td></td><td></td><td></td><td>B1</td><td>CW/OU/NII; BOb; O</td></tr><tr><td>44</td><td>2G1045666+0028085</td><td>2.889±5.5 10<sup>-2</sup></td><td>0.647±7.6 10<sup>-2</sup></td><td>0.511±8.0 10<sup>-2</sup></td><td>0.501±6.6 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1</td><td>CW/OU/NII; O9.5V; O</td></tr><tr><td>45</td><td>2G1046762+0036328</td><td>4.025±7.4 10<sup>-2</sup></td><td>0.350±0.130</td><td>0.259±0.114</td><td>0.176±4.9 10<sup>-2</sup></td><td></td><td>0.328±6.6 10<sup>-2</sup></td><td>0.432±0.171</td><td>0.192±8.5 10<sup>-2</sup></td><td>0.469±0.164</td><td></td><td></td><td></td><td>B1V; B2V; O; O</td></tr><tr><td>46</td><td>2G1066374+0037825</td><td>1.396±6.2 10<sup>-2</sup></td><td>0.939±7.2 10<sup>-2</sup></td><td>0.542±6.9 10<sup>-2</sup></td><td>0.626±4.6 10<sup>-2</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>B1B; B</td><td>CW/OU/NII; B1.5II; B</td></tr><tr><td>47</td><td>2G1091158+0068191</td><td>non gauss. em.</td><td>1.130±9.9 10<sup>-2</sup></td><td>0.968±0.108</td><td>0.727±0.127</td><td></td><td>0.466±0.159</td><td></td><td>0.961±0.186</td><td>0.837±0.178</td><td></td><td></td><td></td><td>B2I</td><td>CW/OU/NII; B1.5II; B</td></tr></tbody></table>

Notes: Hell6680 line are indicated by a "*". In the last column: "SHII" means that the SHII 4552-5568-4575 lines are detected. "CIII/OU/NII" means that the lines CIII 4650 / OII 4640-4650 blend and NII 4631 are well detected (but not resolved). In the last column the spectral types are from Simbad (the references are: Hardorp et al. 1964; Georgelin et al. 1973; Uranova 1977, citate/Nesterov 95; Kobelnicky et al. 2010; Comerón & Pasquali 2012). Maiz Apellániz et al. 2016; Berlanas et al. 2018, 2020; Chick et al. 2020; Dinçel et al. 2024) or (in bold) from GAIA DR3 astrophysical parameters Table (see https://gea.esac.esa.int/archive/documentation/GDR3/).
<table><caption>Table A.4. Adopted and measured physical quantities</caption><thead><tr><th>nb.</th><th>RAJ2000</th><th>DECJ2000</th><th>A<sub>v</sub> (mag)</th><th>Vwind (km/s)</th><th>Teff (K)</th><th>Mass (M<sub>☉</sub>)</th><th>Distance (kpc)</th><th>log(L) (L<sub>☉</sub>)</th><th>log(M) (M<sub>☉</sub>yr<sup>-1</sup>)</th><th>R<sub>0</sub> (pc)</th><th>Vt (km/s)</th><th>PA (deg.)</th><th>ΔP<sub>4</sub> (deg.)</th><th>n<sub>low</sub> cm<sup>-3</sup></th><th>n<sub>up</sub> cm<sup>-3</sup></th></tr></thead><tbody><tr><td>1</td><td>280.141</td><td>-4.899</td><td>5.5</td><td>690.0</td><td>17579.23</td><td>23.0</td><td>1.06<sup>+0.04</sup><sub>-0.03</sub></td><td>4.466</td><td>-7.451</td><td>0.06±0.01</td><td>10.8±5.7</td><td>181.6</td><td>78.0</td><td>5.5</td><td>76.7</td></tr><tr><td>2</td><td>281.114</td><td>-4.763</td><td>7.0</td><td></td><td></td><td></td><td>4.64<sup>+0.03</sup><sub>-0.03</sub></td><td></td><td></td><td>0.31±0.04</td><td>67.6±8.4</td><td>263.0</td><td>126.4</td><td>0.9</td><td>10.5</td></tr><tr><td>3</td><td>280.839</td><td>-3.712</td><td></td><td></td><td></td><td></td><td>2.45<sup>+0.17</sup><sub>-0.18</sub></td><td></td><td></td><td>0.20±0.02</td><td>23.3±2.2</td><td>11.8</td><td>142.3</td><td>2.6</td><td>34.5</td></tr><tr><td>4</td><td>281.387</td><td>-2.487</td><td>7.2</td><td></td><td></td><td></td><td>1.72<sup>+0.084</sup><sub>-0.087</sub></td><td></td><td></td><td>0.20±0.02</td><td>48.8±3.3</td><td>28.4</td><td>279.9</td><td>4.0</td><td>54.0</td></tr><tr><td>5</td><td>282.170</td><td>-2.608</td><td>6.0</td><td>180.0</td><td>8609.93</td><td>13.0</td><td>0.61<sup>+0.005</sup><sub>-0.006</sub></td><td>1.990</td><td>-12.946</td><td>0.02±0.01</td><td>8.2±0.1</td><td>207.3</td><td>45.7</td><td>8.9</td><td>98.1</td></tr><tr><td>6</td><td>281.383</td><td>-2.019</td><td>4.1</td><td>2063.67</td><td>33570.0</td><td>29.6</td><td>8.60<sup>+0.006</sup><sub>-0.006</sub></td><td>6.017</td><td>-4.731</td><td>1.23±0.05</td><td></td><td>320.6</td><td>0.6</td><td>0.8</td><td>24.8</td></tr><tr><td>7</td><td>281.386</td><td>-1.283</td><td>7.5</td><td></td><td></td><td></td><td>2.67<sup>+0.18</sup><sub>-0.20</sub></td><td></td><td></td><td>0.25±0.03</td><td>10.1±1.2</td><td>264.8</td><td>2.3</td><td>33.2</td><td></td></tr><tr><td>8</td><td>282.826</td><td>-1.991</td><td>5.5</td><td>380.0</td><td>16595.86</td><td>9.6</td><td>2.92<sup>+0.16</sup><sub>-0.18</sub></td><td>3.822</td><td>-8.414</td><td>0.17±0.03</td><td>18.8±1.2</td><td>234.6</td><td>1.5</td><td>13.7</td><td></td></tr><tr><td>9</td><td>283.399</td><td>+0.954</td><td></td><td></td><td></td><td></td><td>1.66<sup>+0.18</sup><sub>-0.20</sub></td><td></td><td></td><td>0.68±0.02</td><td>51.8±8.7</td><td>191.6</td><td>78.3</td><td>3.4</td><td>43.4</td></tr><tr><td>10</td><td>285.537</td><td>+3.513</td><td>6.3</td><td>1100.0</td><td>21379.62</td><td>24.0</td><td>4.27<sup>+0.18</sup><sub>-0.18</sub></td><td>4.841</td><td>-6.601</td><td>0.19±0.04</td><td>22.2±2.4</td><td>307.9</td><td>46.5</td><td>1.0</td><td>12.0</td></tr><tr><td>11</td><td>286.294</td><td>+5.164</td><td>5.4</td><td>490.0</td><td>15995.58</td><td>22.0</td><td>2.53<sup>+0.07</sup><sub>-0.07</sub></td><td>4.080</td><td>-8.343</td><td>0.23±0.02</td><td>22.8±1.4</td><td>185.4</td><td>51.9</td><td>2.1</td><td>28.5</td></tr><tr><td>12</td><td>287.551</td><td>+8.236</td><td>4.0</td><td>180.0</td><td>8609.93</td><td>13.0</td><td>0.55<sup>+0.065</sup><sub>-0.065</sub></td><td>1.653</td><td>-13.691</td><td>0.04±0.01</td><td>3.0±0.1</td><td>210.4</td><td>111.2</td><td>8.5</td><td>80.1</td></tr><tr><td>13</td><td>288.783</td><td>+12.702</td><td>7.2</td><td>1853.0</td><td>29991.62</td><td>17.5</td><td>5.06<sup>+0.006</sup><sub>-0.006</sub></td><td>5.237</td><td>-6.259</td><td>0.26±0.05</td><td>120.6±19.6</td><td>286.2</td><td>8.2</td><td>0.9</td><td>14.9</td></tr><tr><td>14</td><td>290.085</td><td>+13.369</td><td>7.7</td><td>865.0</td><td>19364.21</td><td>23.5</td><td>3.86<sup>+0.034</sup><sub>-0.034</sub></td><td>4.384</td><td>-7.628</td><td>0.11±0.04</td><td>73.6±10.0</td><td>299.6</td><td>37.9</td><td>1.1</td><td>28.6</td></tr><tr><td>15</td><td>290.959</td><td>+15.477</td><td>7.4</td><td>490.0</td><td>15995.58</td><td>22.0</td><td>3.33<sup>+0.034</sup><sub>-0.034</sub></td><td>4.091</td><td>-8.319</td><td>0.62±0.03</td><td>34.6±3.4</td><td>92.7</td><td>165.8</td><td>1.3</td><td>28.6</td></tr><tr><td>16</td><td>291.629</td><td>+15.710</td><td>6.9</td><td>490.0</td><td>15995.58</td><td>22.0</td><td>4.00<sup>+0.031</sup><sub>-0.031</sub></td><td>4.151</td><td>-8.187</td><td>0.08±0.04</td><td>81.6±8.2</td><td>283.4</td><td>16.4</td><td>1.0</td><td>26.4</td></tr><tr><td>17</td><td>292.436</td><td>+18.196</td><td>5.8</td><td>1100.0</td><td>21379.62</td><td>24.0</td><td>4.69<sup>+0.040</sup><sub>-0.040</sub></td><td>4.722</td><td>-6.863</td><td>0.25±0.05</td><td>31.3±3.6</td><td>126.5</td><td>126.2</td><td>0.9</td><td>18.9</td></tr><tr><td>18</td><td>292.026</td><td>+19.255</td><td>4.0</td><td>1642.0</td><td>26835.4</td><td>20.0</td><td>9.32<sup>+2.021</sup><sub>-2.021</sub></td><td>5.016</td><td>-6.976</td><td>0.38±0.09</td><td></td><td>95.9</td><td>0.3</td><td>7.4</td><td></td></tr><tr><td>19</td><td>292.274</td><td>+19.143</td><td>5.5</td><td>490.0</td><td>15995.58</td><td>22.0</td><td>6.43<sup>+0.027</sup><sub>-0.027</sub></td><td>3.991</td><td>-8.539</td><td>0.15±0.06</td><td>102.8±19.5</td><td>222.1</td><td>68.9</td><td>0.5</td><td>10.0</td></tr><tr><td>20</td><td>292.112</td><td>+19.399</td><td>5.5</td><td>2161.32</td><td>31846.0</td><td>22.04</td><td>6.49<sup>+0.027</sup><sub>-0.027</sub></td><td>5.420</td><td>-5.938</td><td>0.73±0.06</td><td>79.1±9.2</td><td>61.8</td><td>211.0</td><td>0.5</td><td>6.9</td></tr><tr><td>21</td><td>293.185</td><td>+19.978</td><td>2.9</td><td>370.0</td><td>13551.89</td><td>20.0</td><td>2.73<sup>+0.070</sup><sub>-0.070</sub></td><td>3.662</td><td>-9.290</td><td>0.14±0.03</td><td>34.8±1.4</td><td>21.6</td><td>252.5</td><td>1.3</td><td>37.4</td></tr><tr><td>22</td><td>294.919</td><td>+20.502</td><td>3.5</td><td>2442.53</td><td>39249.0</td><td>37.35</td><td>7.28<sup>+0.092</sup><sub>-0.092</sub></td><td>5.490</td><td>-5.895</td><td>1.50±0.07</td><td>80.7±8.9</td><td>131.2</td><td>92.8</td><td>0.4</td><td>9.5</td></tr><tr><td>23</td><td>295.430</td><td>+23.267</td><td>8.0</td><td>1100.0</td><td>21379.62</td><td>24.0</td><td>1.65<sup>+0.03</sup><sub>-0.03</sub></td><td>4.416</td><td>-7.540</td><td>0.16±0.02</td><td></td><td>149.9</td><td>24.2</td><td>36.4</td><td>44.2</td></tr><tr><td>24</td><td>296.184</td><td>+23.852</td><td>3.4</td><td>690.0</td><td>17579.23</td><td>23.0</td><td>1.98<sup>+0.089</sup><sub>-0.089</sub></td><td>4.066</td><td>-8.334</td><td>0.20±0.02</td><td>15.2±0.5</td><td>199.1</td><td>52.8</td><td>1.9</td><td>36.4</td></tr><tr><td>25</td><td>295.931</td><td>+24.058</td><td>3.7</td><td>490.0</td><td>15995.58</td><td>22.0</td><td>5.69<sup>+1.042</sup><sub>-0.042</sub></td><td>4.817</td><td>-6.714</td><td>0.26±0.05</td><td></td><td>220.5</td><td>0.6</td><td>12.4</td><td></td></tr><tr><td>26</td><td>296.526</td><td>+24.187</td><td>7.0</td><td>1100.0</td><td>21379.62</td><td>24.0</td><td>8.99<sup>+2.032</sup><sub>-0.032</sub></td><td>5.659</td><td>-4.793</td><td>0.34±0.09</td><td></td><td>58.6</td><td>0.4</td><td>7.9</td><td></td></tr><tr><td>27</td><td>296.839</td><td>+24.562</td><td>3.0</td><td>370.0</td><td>12589.25</td><td>19.0</td><td>2.07<sup>+0.066</sup><sub>-0.066</sub></td><td>3.699</td><td>-9.213</td><td>0.11±0.02</td><td>11.5±0.4</td><td>329.1</td><td>76.2</td><td>36.1</td><td></td></tr><tr><td>28</td><td>297.169</td><td>+25.479</td><td>6.7</td><td>1100.0</td><td>21379.62</td><td>24.0</td><td>5.62<sup>+0.051</sup><sub>-0.051</sub></td><td>4.942</td><td>-6.378</td><td>0.78±0.05</td><td></td><td>240.0</td><td>0.6</td><td>14.8</td><td></td></tr><tr><td>29</td><td>298.055</td><td>+26.253</td><td>6.3</td><td>1000.0</td><td>24490.63</td><td>17.4</td><td>4.06<sup>+0.038</sup><sub>-0.038</sub></td><td>4.287</td><td>-7.614</td><td>0.18±0.04</td><td>55.3±5.0</td><td>102.0±2.8</td><td>7.7</td><td>19.2</td><td></td></tr><tr><td>30</td><td>298.381</td><td>+28.140</td><td></td><td></td><td></td><td></td><td>1.93<sup>+0.026</sup><sub>-0.026</sub></td><td></td><td></td><td>0.18±0.02</td><td>102.8±2.8</td><td>328.7</td><td>78.4</td><td>1.7</td><td>38.6</td></tr><tr><td>31</td><td>307.342</td><td>+37.929</td><td>7.5</td><td>2396.4</td><td>38192.0</td><td>34.53</td><td>2.14<sup>+0.084</sup><sub>-0.084</sub></td><td>4.775</td><td>-7.432</td><td>0.11±0.02</td><td>48.6±2.4</td><td>274.7</td><td>18.6</td><td>1.5</td><td>27.7</td></tr><tr><td>32</td><td>304.688</td><td>+39.406</td><td>2.2</td><td>2103.99</td><td>30789.0</td><td>20.55</td><td>1.39<sup>+0.07</sup><sub>-0.07</sub></td><td>4.562</td><td>-7.828</td><td>0.13±0.01</td><td></td><td>107.3</td><td>2.5</td><td>33.8</td><td></td></tr><tr><td>33</td><td>305.725</td><td>+39.208</td><td>4.0</td><td>1100.0</td><td>21379.62</td><td>24.0</td><td>1.65<sup>+0.04</sup><sub>-0.04</sub></td><td>5.078</td><td>-6.078</td><td>0.20±0.02</td><td>11.7±0.4</td><td>91.8</td><td>2.4</td><td>39.2</td><td></td></tr><tr><td>34</td><td>306.242</td><td>+39.458</td><td>5.5</td><td>1905.62</td><td>30626.0</td><td>24.19</td><td>1.63<sup>+0.030</sup><sub>-0.030</sub></td><td>5.344</td><td>-6.213</td><td>0.68±0.02</td><td>7.3±0.2</td><td>44.4</td><td>220.1</td><td>2.6</td><td>39.6</td></tr><tr><td>35</td><td>306.823</td><td>+39.742</td><td>4.3</td><td>2269.31</td><td>35019.0</td><td>27.45</td><td>1.63<sup>+0.034</sup><sub>-0.034</sub></td><td>5.034</td><td>-6.800</td><td>0.36±0.02</td><td>13.4±0.4</td><td>56.7</td><td>229.2</td><td>2.4</td><td>51.1</td></tr><tr><td>36</td><td>306.604</td><td>+40.028</td><td>5.0</td><td>1810.82</td><td>32274.0</td><td>31.54</td><td>1.66<sup>+0.032</sup><sub>-0.032</sub></td><td>4.635</td><td>-7.846</td><td>0.22±0.02</td><td>8.8±2.0</td><td>210.3</td><td>109.0</td><td>2.3</td><td>45.4</td></tr><tr><td>37</td><td>307.064</td><td>+40.735</td><td>4.4</td><td>690.0</td><td>17579.23</td><td>23.0</td><td>1.90<sup>+0.032</sup><sub>-0.032</sub></td><td>4.381</td><td>-7.640</td><td>0.14±0.02</td><td>3.4±0.1</td><td>297.0</td><td>190.7</td><td>2.2</td><td>49.3</td></tr><tr><td>38</td><td>308.493</td><td>+40.328</td><td>7.5</td><td>490.0</td><td>15995.58</td><td>22.0</td><td>1.72<sup>+0.020</sup><sub>-0.020</sub></td><td>4.571</td><td>-7.258</td><td>0.50±0.02</td><td></td><td>287.7</td><td>2.0</td><td>50.1</td><td></td></tr><tr><td>39</td><td>308.730</td><td>+40.579</td><td>8.8</td><td>800.0</td><td>20183.66</td><td>14.8</td><td>1.58<sup>+0.019</sup><sub>-0.019</sub></td><td>4.268</td><td>-7.607</td><td>0.17±0.02</td><td>4.4±0.3</td><td>345.6</td><td>2.0</td><td>53.1</td><td></td></tr><tr><td>40</td><td>307.126</td><td>+42.010</td><td>1.5</td><td>500.0</td><td>18281.0</td><td>12.2</td><td>1.43<sup>+0.004</sup><sub>-0.004</sub></td><td>4.179</td><td>-7.756</td><td>0.11±0.01</td><td>11.2±0.5</td><td>145.5</td><td>31.8</td><td>2.7</td><td>70.4</td></tr><tr><td>41</td><td>309.707</td><td>+41.424</td><td>6.7</td><td>370.0</td><td>12182.65</td><td>18.5</td><td>1.51<sup>+0.003</sup><sub>-0.003</sub></td><td>2.970</td><td>-10.823</td><td>0.17±0.01</td><td>10.2±0.4</td><td>258.0</td><td>34.0</td><td>2.0</td><td>53.3</td></tr><tr><td>42</td><td>309.747</td><td>+42.044</td><td>5.5</td><td>380.0</td><td>14689.26</td><td>21.0</td><td>1.44<sup>+0.002</sup><sub>-0.002</sub></td><td>3.856</td><td>-8.844</td><td>0.16±0.01</td><td>4.8±0.1</td><td>178.1</td><td>83.3</td><td>2.1</td><td>55.1</td></tr><tr><td>43</td><td>334.137</td><td>+59.392</td><td>5.0</td><td>2320.85</td><td>36077.0</td><td>29.59</td><td>4.10<sup>+0.015</sup><sub>-0.015</sub></td><td>5.378</td><td>-6.062</td><td>0.60±0.04</td><td></td><td>296.3</td><td>45.7</td><td>0.7</td><td>12.6</td></tr><tr><td>44</td><td>335.511</td><td>+58.719</td><td>2.8</td><td>1960.37</td><td>31607.0</td><td>25.83</td><td>2.13<sup>+0.038</sup><sub>-0.038</sub></td><td>4.958</td><td>-7.053</td><td>0.17±0.02</td><td>7.0±1.2</td><td>315.0</td><td>1.3</td><td>13.2</td><td></td></tr><tr><td>45</td><td>337.578</td><td>+60.642</td><td>3.3</td><td>200.0</td><td>10495.42</td><td>16.0</td><td>2.27<sup>+0.016</sup><sub>-0.016</sub></td><td>4.280</td><td>-7.914</td><td>0.33±0.02</td><td></td><td>89.4</td><td>1.1</td><td>21.6</td><td></td></tr><tr><td>46</td><td>339.844</td><td>+59.007</td><td>2.9</td><td>490.0</td><td>18281.0</td><td>12.2</td><td>2.63<sup>+0.017</sup><sub>-0.017</sub></td><td>4.431</td><td>-7.198</td><td>0.28±0.03</td><td></td><td>67.8</td><td>0.9</td><td>8.9</td><td></td></tr><tr><td>47</td><td>343.939</td><td>+60.411</td><td>5.3</td><td>800.0</td><td>18071.74</td><td>18.9</td><td>3.74<sup>+0.015</sup><sub>-0.015</sub></td><td>4.389</td><td>-7.495</td><td>0.26±0.04</td><td></td><td>118.3</td><td>0.7</td><td>6.6</td><td></td></tr></tbody></table>
<table><caption>**Table A.5.** Calculated quantities.</caption><thead><tr><th>nb.</th><th>V<sub>flow</sub><br>(for n<sub>low</sub>)<br>(km/s)</th><th>V<sub>flow</sub><br>(for n<sub>up</sub>)<br>(km/s)</th><th>V<sub>a</sub><br>(for n<sub>low</sub>)<br>(km/s)</th><th>V<sub>a</sub><br>(for n<sub>up</sub>)<br>(km/s)</th><th>density<br>cm<sup>-3</sup></th></tr></thead><tbody><tr><td>2</td><td>63.8</td><td>46.9</td><td>24.1</td><td>7.2</td><td></td></tr><tr><td>5</td><td></td><td></td><td></td><td></td><td>0.003</td></tr><tr><td>6</td><td></td><td></td><td>269.1</td><td>47.8</td><td></td></tr><tr><td>8</td><td>19.4</td><td>13.9</td><td>8.5</td><td>2.8</td><td></td></tr><tr><td>10</td><td></td><td></td><td></td><td></td><td>34.32</td></tr><tr><td>11</td><td></td><td></td><td></td><td></td><td>0.18</td></tr><tr><td>12</td><td>1.1</td><td>1.1</td><td>0.01</td><td>0.01</td><td></td></tr><tr><td>13</td><td></td><td></td><td></td><td></td><td>2.40</td></tr><tr><td>14</td><td></td><td></td><td></td><td></td><td>0.72</td></tr><tr><td>15</td><td>36.8</td><td>34.1</td><td>3.3</td><td>0.7</td><td></td></tr><tr><td>16</td><td></td><td></td><td></td><td></td><td>0.18</td></tr><tr><td>17</td><td>94.7</td><td>34.9</td><td>76.2</td><td>16.6</td><td></td></tr><tr><td>18</td><td></td><td></td><td>86.6</td><td>18.8</td><td></td></tr><tr><td>19</td><td>-21.4</td><td>-33.6</td><td>16.0</td><td>3.7</td><td>0.01</td></tr><tr><td>20</td><td>209.4</td><td>106.9</td><td>141.2</td><td>38.6</td><td></td></tr><tr><td>21</td><td>14.9</td><td>11.7</td><td>4.0</td><td>0.8</td><td></td></tr><tr><td>22</td><td>88.9</td><td>21.6</td><td>85.4</td><td>18.0</td><td></td></tr><tr><td>23</td><td></td><td></td><td>32.7</td><td>7.7</td><td></td></tr><tr><td>24</td><td>0.6</td><td>-6.9</td><td>9.8</td><td>2.2</td><td></td></tr><tr><td>25</td><td></td><td></td><td>68.2</td><td>15.6</td><td></td></tr><tr><td>26</td><td></td><td></td><td>933.6</td><td>210.4</td><td></td></tr><tr><td>27</td><td>2.4</td><td>-1.7</td><td>5.1</td><td>1.1</td><td></td></tr><tr><td>28</td><td></td><td></td><td>52.4</td><td>10.7</td><td></td></tr><tr><td>29</td><td></td><td></td><td></td><td></td><td>0.55</td></tr><tr><td>31</td><td></td><td></td><td></td><td></td><td>7.57</td></tr><tr><td>32</td><td></td><td></td><td>40.2</td><td>11.0</td><td></td></tr><tr><td>33</td><td>152.8</td><td>45.4</td><td>143.2</td><td>35.8</td><td></td></tr><tr><td>34</td><td>52.4</td><td>17.6</td><td>46.8</td><td>12.0</td><td></td></tr><tr><td>35</td><td>59.3</td><td>19.9</td><td>50.5</td><td>11.0</td><td></td></tr><tr><td>36</td><td>25.1</td><td>7.8</td><td>22.2</td><td>5.0</td><td></td></tr><tr><td>37</td><td>31.5</td><td>9.3</td><td>28.1</td><td>5.9</td><td></td></tr><tr><td>38</td><td></td><td></td><td>11.1</td><td>2.2</td><td></td></tr><tr><td>39</td><td></td><td></td><td></td><td></td><td>77.48</td></tr><tr><td>40</td><td></td><td></td><td></td><td></td><td>13.02</td></tr><tr><td>41</td><td></td><td></td><td></td><td></td><td>0.004</td></tr><tr><td>42</td><td>4.0</td><td>0.3</td><td>4.6</td><td>0.9</td><td></td></tr><tr><td>43</td><td></td><td></td><td></td><td></td><td>3.42</td></tr><tr><td>44</td><td></td><td></td><td></td><td></td><td>283.55</td></tr><tr><td>45</td><td></td><td></td><td>6.7</td><td>1.5</td><td></td></tr><tr><td>46</td><td></td><td></td><td>32.0</td><td>10.1</td><td></td></tr><tr><td>47</td><td></td><td></td><td></td><td></td><td>0.44</td></tr></tbody></table>
Appendix B: Figures of the stellar SED

<figure><img src="image_6.png" /><figcaption>**Fig. B.1.** Spectral energy distribution (SED) of the SBN star candidates classified as OB stars. The log of the power law $F_v \propto v^\beta$ fit is shown and the slope which is the $\beta$ value (and its uncertainty) is given at the top of the panels.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Fig. B.1. Continued.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>Fig. B.1. Continued.</figcaption></figure>
**Appendix C: SBN images with transverse velocity vectors and additional figures.**

<figure><img src="image_9.png" /><figcaption>**Fig. C.1.** WISE 22µm or *Spitzer*-MIPS 24 µm images of the studied SBN. The candidate star position and the transverse velocity vector (arbitrary length) are shown as circles and arrows, respectively. The coordinates around the frames are the Galactic coordinates (in degrees).</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>Fig. C.1. Continued.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>Fig. C.1. Continued.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Fig. C.1. Continued.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. C.2.** IRAC-8µm (left) and MIPS 24µm (right) images of N131. The displayed arrow is the transverse velocity vector (arbitrary length) for star #29. The images are oriented according to Galactic coordinates.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Fig. C.3.** DSS-Red (left) and MIPS 24µm (right) images of Sh2-86. The star #6 is indicated by the cross. The images are oriented according to Galactic coordinates.</figcaption></figure>