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

HAL
open science

# Study of a sample of faint Be stars in the exofield of CoRoT. III. Global spectroscopic characterization and astrophysical parameters of the central stars

J. Zorec, A. M. Hubert, C. Martayan, Y. Frémat

► To cite this version:

J. Zorec, A. M. Hubert, C. Martayan, Y. Frémat. Study of a sample of faint Be stars in the exofield of CoRoT. III. Global spectroscopic characterization and astrophysical parameters of the central stars. Astronomy & Astrophysics - A&A, 2023, 676, 10.1051/0004-6361/202346018. insu-04848725

HAL Id: insu-04848725

https://insu.hal.science/insu-04848725v1

Submitted on 19 Dec 2024

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Study of a sample of faint Be stars in the exofield of CoRoT

## III. Global spectroscopic characterization and astrophysical parameters of the central stars*

J. Zorec<fnref n="1" />, A. M. Hubert<fnref n="2" />, C. Martayan<fnref n="3" />, and Y. Frémat<fnref n="4" />

Received 30 January 2023 / Accepted 12 June 2023

## ABSTRACT

**Context.** The search and interpretation of non-radial pulsations from Be star light curves observed with the CoRoT satellite requires high-quality stellar astrophysical parameters.

**Aims.** The present work is devoted to the spectroscopic study of a sample of faint Be stars observed by CoRoT in the fourth long run (LRA02).

**Methods.** The astrophysical parameters were determined from the spectra in the λλ4000–4500 Å wavelength domain observed with the VLT/FLAMES instruments at ESO. Spectra were fitted with models of stellar atmospheres using our GIRFIT package. Spectra obtained in the λλ6400–7200 Å wavelength domain enabled the confirmation or, otherwise, a first identification of Be star candidates.

**Results.** The apparent parameters ($T_{eff}$, log g, Vsin i) for a set of 19 B and Be stars were corrected for the effects induced by the rapid rotation. These allowed us to determine: (1) stellar masses that are in agreement with those measured for detached binary systems; (2) spectroscopic distances that agree with the Gaia parallaxes; and (3) centrifugal/gravity equatorial force ratios of ~0.6–0.7, which indicate that our Be stars are subcritical rotators. A study of the Balmer Hα, Hγ and Hδ emission lines produced: (1) extents of the circumstellar disk (CD) emitting regions that agree with the interferometric inferences in other Be stars; (2) R-dependent exponents $n(R) = \ln[\rho(R)/\rho_o]/\ln(R_o/R)$ of the CD radial density distributions; and (3) CD base densities $\rho_o$ similar to those inferred in other recent works.

**Conclusions.** The Hγ and Hδ emission lines are formed in CD layers close to the central star. These lines produced a different value of the exponent $n(R)$ than assumed for Hα. Further detailed studies of Hγ and Hδ emission lines could reveal the physical properties of regions where the viscous transport of angular momentum to the remaining CD regions is likely to originate from. The subcritical rotation of Be stars suggests that their huge discrete mass-ejections and concomitant non-radial pulsations might have a common origin in stellar envelope regions that become unstable to convection due to rotation. If it is proven that the studied Be stars are products of binary mass transfer phases, the errors induced on the estimated $T_{eff}$ by the presence of stripped sub-dwarf O/B companions are not likely to exceed their present uncertainties.

**Key words.** stars: early-type – stars: emission-line, Be – stars: fundamental parameters – stars: oscillations

## 1. Introduction

The photometric data obtained by the CoRoT satellite (Weiss et al. 2004; Baglin et al. 2006) concern a wide range of variable objects, including Be stars, which are rapidly, albeit subcritical, rotating objects (Cranmer 2005; Zorec et al. 2016). During their main sequence lifespan, these stars build up a circumstellar disk (hereafter CD), (Jaschek et al. 1981; Collins 1987), whose presence is revealed by line emissions and continuum flux excess, mainly in the visible and infrared spectral domains (Rivinius et al. 2013).

The CoRoT photometric data of Be stars are particularly relevant because they carry signatures of their non-radial pulsation modes. As pulsating objects, Be stars lie in the instability region of the Hertzsprung-Russell (HR) diagram between (or even partially shearing) the instability regions of β Cephei stars and the slowly pulsating-B stars (SPBs, Miglio et al. 2007; Walczak et al. 2015; Saio et al. 2017; Szewczuk & Daszyńska-Daszkiewicz 2017; Burssens et al. 2020). Unlike Be stars, SPBs have often been considered slow rotators. However, long ground-based photometric surveys (e.g., Salmon et al. 2014) and recent four-year data provided by the Kepler mission (Pedersen et al. 2021, and references therein) have enabled the identification of rapid rotators among B-star pulsators (previously considered SPBs), where properties of their internal structure and pulsations might perhaps become an approximation for Be stars (Moravveji et al. 2016).

<footnote n="1">¹ Sorbonne Université, CNRS, UPMC, UMR7095 Institut d'Astrophysique de Paris, 98bis Bd. Arago, 75014 Paris, France
e-mail: zorec@iap.fr</footnote>

<footnote n="2">² GEPI, Observatoire de Paris, PSL Research University, CNRS UMR 8111, 5 place Jules Janssen, 92190 Meudon, France
e-mail: Anne-Marie.Hubert@obspm.fr</footnote>

<footnote n="3">³ European Organization for Astronomical Research in the Southern Hemisphere, Alonso de Cordova 3107, Vitacura, Santiago de Chile, Chile</footnote>

<footnote n="4">⁴ Royal Observatory of Belgium, 3 Av. Circulaire, 1180 Bruxelles, Belgium</footnote>

<footnote>* The complete version of Tables D.1 and D.2 are only available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/676/A81</footnote>
The analysis of the CoRoT photometric data of Be stars carried out by Huat et al. (2009) and Semaan et al. (2018) revealed complex frequency spectra exhibiting multiple-periodicity and groups of closely separated frequencies, as well as isolated frequencies. This same behavior was also observed in data of Be stars obtained in other space missions: MOST (Walker et al. 2003, 2005a,b), SMEI (Jackson et al. 2004; Goss et al. 2011), KEPLER (Borucki et al. 2010; Kurtz et al. 2015; Pápics et al. 2017), BRITE-Constellation (Weiss et al. 2014; Baade et al. 2016), and TESS (Ricker et al. 2016; Labadie-Bartz et al. 2022).

The groups of closely separated frequencies can be considered as combinations of g-mode pulsations, which cannot be produced by spots in the stellar surface (Kurtz et al. 2015). One of the most important findings in the CoRoT data of HD 49330 (Huat et al. 2009), later also observed in another 7 Be stars analyzed by (Semaan et al. 2018), is the presence of light outbursts accompanied by changes mainly of the photometric power spectrum. They are supposed to be manifestations of discontinuous mass-ejections, similar to those reported by Hubert & Floquet (1998) using HIPPARCOS photometry, as well as in stars of the LMC and SMC observed by the MACHO Cook et al. (1995); Keller et al. (2002) and the OGLE Mennickent et al. (2002), de Wit et al. (2006) surveys. Although these ejections are supposed to provide the required mass to build up a CD as described by the viscous decretion model (VDD), (Lee et al. 1991; Ghoreyshi et al. 2021; Haubois et al. 2012), the mechanisms causing the ejections proper are still heavily debated. On the one hand, from theoretical grounds the combination of prograde g-modes of stellar pulsation could lead to sporadic mass ejections at circumstances when stellar rotation is near critical (Osaki 1986; Kee et al. 2016; Saio et al. 2017). On the other hand, based on spectroscopic studies, Baade et al. (2016) suggested a possible relation between a stellar oscillation called "difference frequency", which somehow should result from the combination of two other observed close pulsation frequencies, with sporadic mass ejection events. A series of correlations (e.g., Semaan et al. 2018), theoretical predictions of pulsating modes (e.g., Saio 2013), and recent results (e.g., Pápics et al. 2017) confirm that Be stars should be considered as genuine non-radial pulsators.

Two preceding papers by Semaan et al. (2013, 2018) were dedicated to spectroscopic studies and analysis of non-radial pulsation modes in faint Be stars observed in the initial run (IRA01) and the first two long runs (LRC01 and LRA01) of the CoRoT mission. The present work deals with a sample of faint Be stars belonging to the fourth long run (LRA02), whose main objective is to achieve a global spectroscopic characterization for them. The specific aim of the present paper is to determine reliable stellar astrophysical parameters of a new set of faint Be stars. In this work we consider Be stars to behave as single objects.

Although the astrophysical parameters are here determined using spectra obtained in the 4000–4500 Å wavelength domain, spectra in a wavelength region near the $H\alpha$ line were also obtained in order to confirm the Be star nature of the studied stars. To take full advantage of the obtained spectroscopic data, we also inquire the physical structure of the CD in the studied Be stars. Our approach is based on first principles of radiation transfer that control the formation of Balmer lines. It enables to probe the extents, and the density distribution in the CDs of the studies Be stars.

The present work is organized as follows. Section 2 briefly identifies the sky region of the studied Be stars and specifies the source where the CoRoT data were extracted. Section 3 presents the spectroscopic observations of the program objects. Section 4 gives a short description of the spectral characteristics of the

selected Be stars. Section 5 describes the methods used to determine the primary and secondary apparent stellar astrophysical parameters and their treatment for rapid rotation effects. Some correlations between the Balmer line emission intensities and the near-infrared photometric color indices are shown in Sect. 6. A method to describe the physical characteristic of CD of the program stars is detailed in Sect. 7. A short summary of the results obtained is given in Sect. 8.

## 2. CoRoT observations

The telescope of the CoRoT satellite was pointed to the “anticentre” of the Milky Way to observe the LRA02 field whose coordinates are $\alpha = 103°52, δ = -04°38$ and Roll = 06°00. Half of the field was dedicated to observation of bright stars devoted to programs of seismology and the other half to observation of faint stars with magnitudes $11 < V < 17$ mag to additional programs, such as the search for exoplanets. The upper part of the LRA02 field was very close to the lower part of the IRA01 field, as shown in Deleuil et al. (2016, Fig. III.1.1).

The LRA02 exofield contains stars within coordinates $06^h47^m30^s < \alpha < 06^h53^m50^s$ and $-05°51'30'' < δ < -03°06'30''$. In this field, 11 448 stars were observed over 112 days from 16 November 2008 to 8 March 2009. CoRoT data, namely, the light curves with several levels of corrections and headers, of all faint stars are available at CDS<fnref n="1" />.

## 3. GIRAFFE observations

### 3.1. *Observations and reduction of spectra*

A spectroscopic campaign with the VLT/FLAMES instrumentation at ESO<fnref n="2" /> was undertaken to obtain spectra of the faint Be stars identified in the LRA02 exofield of CoRoT. Observations were conducted in 2010 (November and December) and in 2011 (January). At least three spectra were obtained for each target; one spectrum with 3600 s of exposure time in the MEDUSA MODE at medium resolution, $R = 6400$, in the LR2 setup (396.4–456.7 nm), hereafter referred to as the “blue spectrum”, and two or three successive spectra with $R = 8600$ in the LR6 setup (643.8–718.4 nm), hereafter referred to as the “red spectra”, having 780 s of exposure time each. A total of 2209 stars have been observed in 18 fields, of which 2189 ones in two blue and red setups, and 19 stars in only one. We thus dispose spectra for about 19% of the stars that were photometrically observed by CoRoT in the LRA02 exofield. GIRAFFE observations have been conducted at a nearly full moon. Consequently, the stellar spectra are most often polluted by telluric absorption lines and the solar spectrum due to the scattered moonlight. We made the reduction of spectra with the use of the ESO GIRAFFE pipeline in the same way as in Semaan et al. (2013, Paper I) for Be stars studied in the first exofields of CoRoT. Then, we used the IRAF package to separate all the spectra contained in each exposure made in 18 GIRAFFE fields (labeled from *a* to *r*). The reduced spectra were calibrated in wavelength, corrected from heliocentric velocity and cleaned for cosmics. For each star, the red spectra were summed up to increase their S/N.

<footnote n="1">¹ http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=B/corot</footnote>

<footnote n="2">² Program ID: 086.D-0212(A), (PI: C. Neiner).</footnote>
<table><caption>**Table 1.** Be star candidates in the LRA02 exofield of CoRoT.</caption><thead><tr><th>CoRoT ID</th><th>N°</th><th>V (mag)</th><th>R (mag)</th><th>Gaia DR3 ID</th><th>G (mag)</th><th>SpT</th></tr></thead><tbody><tr><td>CoRoT 103000272</td><td>1</td><td>13.244</td><td>13.124</td><td>3105803586041696768</td><td>13.138</td><td>B8 V</td></tr><tr><td>CoRoT 103032255</td><td>2</td><td>13.551</td><td>13.352</td><td>3105810496648942720</td><td>13.400</td><td>B7 II</td></tr><tr><td>CoRoT 110655185</td><td>3</td><td>15.670</td><td>15.300</td><td>3102172212671872256</td><td>15.401</td><td></td></tr><tr><td>CoRoT 110655384</td><td>4</td><td>16.073</td><td>15.750</td><td>3101951348274296704</td><td>15.862</td><td></td></tr><tr><td>CoRoT 110655437</td><td>5</td><td>15.590</td><td>15.640</td><td>3102171250599205248</td><td>15.813</td><td></td></tr><tr><td>CoRoT 110662847</td><td>6</td><td>14.114</td><td>13.780</td><td>3105813279788040320</td><td>13.905</td><td>B3 V</td></tr><tr><td>CoRoT 110663174</td><td>7</td><td>16.464</td><td>15.900</td><td>3102281133042324096</td><td>15.983</td><td></td></tr><tr><td>CoRoT 110663880</td><td>8</td><td>14.513</td><td>13.985</td><td>3105801253879634944</td><td>14.023</td><td>B8 I</td></tr><tr><td>CoRoT 110672515</td><td>9</td><td>16.062</td><td>15.563</td><td>3101912109453608576</td><td>15.567</td><td></td></tr><tr><td>CoRoT 110681176</td><td>10</td><td>12.058</td><td>11.936</td><td>3102777420808111104</td><td>11.985</td><td>B8 V</td></tr><tr><td>CoRoT 110688151</td><td>11</td><td>11.311</td><td>11.150</td><td>3102105206881278848</td><td>11.242</td><td>B6 III</td></tr><tr><td>CoRoT 110747131</td><td>12</td><td>14.180</td><td>14.470</td><td>3101888744831685888</td><td>14.586</td><td></td></tr><tr><td>CoRoT 110751872</td><td>13</td><td>11.402</td><td>11.330</td><td>3102285324930474624</td><td>11.374</td><td>B5 III</td></tr><tr><td>CoRoT 110751876</td><td>14</td><td>11.443</td><td>11.310</td><td>3105805651925926784</td><td>11.324</td><td></td></tr><tr><td>CoRoT 110752156</td><td>15</td><td>15.875</td><td>15.341</td><td>3102336314781346816</td><td>15.246</td><td></td></tr><tr><td>CoRoT 110827583</td><td>16</td><td>15.331</td><td>14.960</td><td>3101956326137090560</td><td>15.087</td><td></td></tr><tr><td>CoRoT 300002611</td><td>17</td><td>16.395</td><td>15.916</td><td>3105402234235943424</td><td>15.643</td><td></td></tr><tr><td>CoRoT 300002834</td><td>18</td><td>16.004</td><td>15.445</td><td>3105481124194678656</td><td>15.441</td><td></td></tr><tr><td>CoRoT 300003290</td><td>19</td><td>15.471</td><td>15.058</td><td>3105385054366798336</td><td>15.175</td><td></td></tr></tbody></table>

**Notes.** Column 2: star number; Col. 3: V magnitude from the AAVSO photometric all-sky survey (APAAS) reported in the UCAC4 Catalog; Col. 4: magnitude in the R band from the CoRoT faint stars database (Vizier) On-line Data Catalog: B/CoRoT; Col. 5: Gaia DR3 ID; Col. 6: magnitude G from the Gaia EDR3 Catalog; Col. 7: Spectral types determined by spectroscopy in Sebastian et al. (2012).

## 3.2. Setting up the faint Be star sample in the LRA02 exofield of CoRoT

In a first step, we made the selection of emission-line stars by inspecting the $H\alpha$ line in the red spectrum. We identified about 40 emission-line stars and confirmed two of them: CoRoT 110751872 (TYC 4808-1005-1) and CoRoT 110751876 (EM*AS 136). Moreover, we observed that the $H\alpha$ line profiles of many fainter B/A stars are often polluted by night sky emission lines (e.g., Hanuschik 2003) which could bring on confusion with regard to the true nature of these objects. We finally selected only 20 Be stars. Other stars are cool or much fainter stars that have noisy spectra polluted by OH lines (Osterbrock et al. 1996) in their $H\alpha$ line profile.

The selected CoRoT Be star candidates are listed in Table 1. They have been identified in UCAC4 catalog (Fourth U.S. Naval Observatory CCD Astrograph Catalog) and Gaia DR3 catalog (Gaia Collaboration 2023). Information on GIRAFFE spectra and conditions of their exposure times are given in Table A.1. The mean signal-to-noise ratios (S/Ns) range from 40 to 400. The fraction of moonlight that contaminates the spectra is also mentioned.

In a second step, we investigated the blue GIRAFFE spectra of all Be star candidates to identify the photospheric lines in order to determine their spectral type and the corresponding photospheric parameters. Most spectra were obtained for a full moon (see Table A.1). The exposure times in the blue wavelength range are longer than those in the red one. Therefore, the spectra of the faintest stars ($15 < V < 17$ mag) were more heavily polluted by sky background in the LR2 setup than in the LR6 one. The contamination due to the moonlight was estimated using the processed data of our objects and of their sky environments, which are available in the ESO Science Spectrum archive. We

subtracted the moonlight contribution to the flux of each faint star. Thanks to these corrections, we could confirm the B spectral type of all stars in Table 1 (19 stars) and the Be character for 18 of them.

## 4. Short description of the spectral characteristics of the program Be stars

We present a short description of the spectral characteristics of each Be star candidate in the LRA02 exofield of CoRoT. Several quantities related with Balmer emission lines and shell features are reported. We note that diffuse interstellar bands (DIBs) are often detected in spectra. Most of them are identified at $\lambda\lambda 4232, 4430, 4501, 6614, 6660, 6993, 7060$, and $7120 \AA$ (Jenniskens & Desert 1994; Galazutdinov et al. 2000). The spectral region around the $H\alpha$ line from $\lambda 6500$ to $6700 \AA$ of the program stars is shown in Fig. 1.

**CoRoT 103000272.** $H\alpha$ shows a faint, double-peaked emission, which is nearly symmetrical ($V \ge R$, $\Delta_{\text{peaks}} = 50 \text{ km s}^{-1}$) and centrally superposed on the broad photospheric line profile. No emission is detected in the blue spectrum. We note the presence of very weak metallic lines (Fe II) in the absorption as well as several DIBs.

**CoRoT 103032255.** $H\alpha$ shows a strong, single and symmetrical emission line profile. The photospheric red lines He I $\lambda\lambda 6678$ and $7065 \AA$ are rather weak. In the blue domain, a very weak emission is present at each side of the core of the photospheric $H\gamma$ line profile.

**CoRoT 110655185.** $H\alpha$ shows a symmetrical and double-peaked emission line profile ($\Delta_{\text{peaks}} = 110 \text{ km s}^{-1}$) of moderate intensity superimposed on the photospheric line profile; other wavelength ranges in the red spectrum are polluted by telluric lines. A weak emission is present at each side of the core of the photospheric $H\gamma$ and $H\delta$ line profiles.

**CoRoT 110655384.** $H\alpha$ depicts a very weak, double-peaked emission with $V \ge R$ and $\Delta_{\text{peaks}} = 290 \text{ km s}^{-1}$ superimposed on the broad photospheric line profile. Weak nebular lines (sky lines) are visible on each $H\alpha$ side. No emission is detected in blue Balmer lines.

**CoRoT 110655437.** from its CoRoT light curve, this faint object is known as an eclipsing binary whose period is 3.96685 d (Klagyivik et al. 2017). In the red spectrum, $H\alpha$ displays a complex line profile with an asymmetric double-peaked emission ($V \gg R$, $\Delta_{\text{peaks}} = 367 \text{ km s}^{-1}$) disturbed by two narrow absorptions or shell features, the first one at $RV = +155 \text{ km s}^{-1}$, the second at $RV = -540 \text{ km s}^{-1}$. He I $\lambda\lambda 6678$ and $7065$ show a double-like structure with a sharp absorption component located at the center of a broad one. Weak nebular lines are present on each side of $H\alpha$. It is difficult to search signatures of a binary component in the blue spectrum, which is noisy and has a lower spectral resolution than the red one. Other comments on this object are presented in Sect. 8.2.

**CoRoT 110662847.** $H\alpha$ displays a symmetrical, double-peaked emission ($\Delta_{\text{peaks}} = 210 \text{ km s}^{-1}$) superimposed on the photospheric line profile. Red He I lines $\lambda\lambda 6678$ and $7065$ are conspicuous. No emission is detected in the blue lines.

**CoRoT 110663174.** $H\alpha$ shows a nearly symmetrical, double-peaked emission line profile ($V \ge R$, $\Delta_{\text{peaks}} = 120 \text{ km s}^{-1}$) of moderate intensity superposed on the photospheric line
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Spectra in the Hα (λ6562 Å) line region of the studied Be stars. This figure also indicates the positions of the He I 6678 line and the diffuse interstellar absorption band at λ6614 Å. The Hα emission component of circumstellar origin is marked in red. The small infilling emission in CoRoT 300002834 is probably not of circumstellar origin.</figcaption></figure>
profile. The blue spectrum of this faint star is very noisy and dis-
turbed by the presence of the full moon in the sky at the epoch
of the observation. After correction of this spectrum from sky
lines contribution, the core of the photospheric Hγ line profile is
disturbed by emission.

CoRoT 110663880. Hα displays a strong emission line profile with a complex, possibly triple structure in its core. The blue spectrum of this faint star is very noisy and strongly affected by the presence of the full moon in the sky at the epoch of the observation. After correction of this spectrum from sky lines contribution, the core of the photospheric Hγ and Hδ line profiles appeared to be disturbed by a double weak emission, itself disturbed by a central shell feature.

CoRoT 110672515. Hα displays a symmetrical, double-peaked emission (Δ<sub>peaks</sub> = 225 km s<sup>-1</sup>) with a deep H shell feature (RV = +34 km s<sup>-1</sup>). Red He I lines are conspicuous; we also note the presence of a metallic shell, mainly of Fe II lines, which depict a weak, double emission on each side of an absorption feature. Hγ and Hδ line profiles are affected by emission on each side of their photospheric core. A shell is also present (H and probably metals) in the blue domain.

CoRoT 110681176. Hα shows a very weak emission, which is nearly symmetrical (V ≤ R, Δ<sub>peaks</sub> = 295 km s<sup>-1</sup>) in the core of the photometric line profile. Blue and red He I lines are weak. No emission is observed in the blue Balmer lines.

CoRoT 110688151. Among the Be stars of the sample, this one displays a significant contamination of its spectra due to lines of circumstellar origin. Hα depicts a strong, closely double-peaked emission which is nearly symmetrical (V ≤ R, Δ<sub>peaks</sub> = 120 km s<sup>-1</sup>). Red He I lines are conspicuous. In the blue domain, the photospheric Hγ and Hδ line profiles are disturbed in their core by a double-peaked emission. Moreover, the strongest Fe II lines of multiplets 27, 28, 37, 38 in the blue spectrum, and of multiplets 40 and 74 in the red one are present in emission with a double structure. Metallic shell lines highly pollute the blue He I and Mg II lines commonly used to determine stellar astrophysical parameters.

CoRoT 110747131. Hα depicts a weak, closely doubled-
peaked emission, which is nearly symmetrical (V ≥ R; Δ<sub>peaks</sub> =
70 km s<sup>-1</sup>), in the core of the photospheric line profile. No
emission is detected in blue lines.

CoRoT 110751872. Hα shows a symmetrical, closely
double-peaked emission of moderate intensity (Δ<sub>peaks</sub> =
100 km s<sup>-1</sup>) superimposed on the photospheric profile. Red He I
lines are conspicuous. No emission is detected in blue lines.

CoRoT 110751876. This star displays lines formed in the circumstellar medium that contaminate its spectra. Hα shows a rather strong, closely double-peaked emission line (Δ<sub>peaks</sub> = 140 km s<sup>-1</sup>), which is symmetrical. In the blue domain, Balmer lines are affected by emission on each side of their photospheric core. Blue He I absorption lines are deep and strong. As in CoRoT 110688151, the stronger lines of the main Fe II multiplets are present with a double emission profile in spectra, which makes particularly difficult the fitting of blue He I and Mg II lines.

CoRoT 110752156. Hα shows a very weak photospheric line profile disturbed on each side by a double-peaked emission (V ≤ R; Δ<sub>peaks</sub> = 525 km s<sup>-1</sup>). Moreover, this complex line

presents a deep and sharp core in absorption (H shell, RV = +35 km s⁻¹). Red He I lines are conspicuous, He I 6678 depicts a double structure that could be explained as due to binarity, DIP or other facts. The blue spectrum of this faint star is strongly affected by the presence of the full moon in the sky at the epoch of the observation. After correction from sky lines contribution, this spectrum remains noisy, however He I lines are conspicuous as well as a shell H feature in Hy.

CoRoT 110827583. Hα displays a weak, symmetrical and double-peaked emission (Δ<sub>peaks</sub> = 170 km s<sup>-1</sup>) in the core of the broad photospheric line profile. A lot of deep and narrow absorption lines, probably of telluric origin, pollute the red spectrum of this faint object. Moreover, weak nebular lines (OH) are also present on each side of the core of Hα. No emission is detected in the blue domain. After the subtraction of moonlight, the blue spectrum corresponds to a B type star.

CoRoT 300002611. Hα is a weak, symmetrical and double-peaked emission (Δ<sub>peaks</sub> = 150 km s<sup>-1</sup>) superimposed to the broad and weak photospheric line. Nebular lines are conspicuous in the vicinity of Hα. Spectra are noisy, the S/N is poor in the bluer part of the blue domain, however weak blue He I lines are detected in this star. Moreover, a weak emission seems to be present at each side of the core of the photospheric Hγ line profile, as well as in Hδ, but with lower intensity.

CoRoT 300002834. Very weak and narrow emission lines (OH) are observed over the wings of the Hα photospheric line; nebular lines [S II] 6717 and 6731 seem also to be present. The emission component in the Hα core could not have a circumstellar origin. Thus, the object is not considered a Be star. We note that fainter early-type stars in the same GIRAFFE field as well as in some others have similar characteristics in their red spectrum.

CoRoT 300003290. Hα shows a strong, single-peaked emission. Very weak sky lines are observed on each side of the line. The photospheric Hγ line is disturbed by a weak, narrow and single-peaked emission in its core. The same is observed in Hδ, but to a lesser degree. As the star is faint, and consequently the blue spectrum is rather noisy, it is difficult to identify whether circumstellar Fe II lines are present or not.

**5. Astrophysical parameters**

The astrophysical parameters derived here are meant to help
to estimate the masses, ages, and radii of the program stars.
They will enable calculations of rotational frequencies to test
the pulsation-instability properties of stars and compare them to
theoretical predictions. The determination of the astrophysical
parameters implies several steps.

First, a correction from the “veiling” effect of spectra dis-
playing emission lines is needed, which is induced by the
disk. Then, the veiling-corrected spectral energy distribution
(SED) allows us to access to the “apparent” stellar astrophysical
parameters. The apparent parameters are those of classical plane-
parallel models of stellar atmospheres that correctly represent
the spectral characteristics of the studied spectral wavelength
region. In our case, this corresponds to the λλ4000–4500 Å
interval, emitted by the observed stellar hemisphere, which is
deformed by the rapid stellar rotation and has non-uniform sur-
face distributions of gravity and effective temperature. Finally,
using models of stellar atmospheres that account for the effects
due to rapid rotation, mostly stellar geometrical deformation and
gravity darkening (GD), the apparent astrophysical parameters
are translated into parent non-rotating counterpart astrophysical parameters (pnrc), which are ultimately used to determine the stellar mass and its age. In this step, we also derive the Vsin i parameter corrected for the GD effect, known as Stoeckley's underestimation. AA avoids bulleted lists in the man text.

## 5.1. The veiling effect

The observed continuum energy distribution of Be stars corresponds to the stellar energy emitted by the central stars, modified by the presence of the CD. According to these cases, the positive or negative flux excess produced by the circumstellar environment can represent a significant fraction of the total flux in the visible spectral domain. The apparent residual intensities of spectral lines are meant to be genuine photospheric components that are usually analyzed to infer the astrophysical stellar parameters; thus, they may not only be affected by emission and absorption components produced in the CD, but also can be over- or under-normalized by the continuum energy distribution modified by the flux excess of circumstellar origin. These perturbations have been known from the Fifties, and are called the “veiling effect” (de Jager & Neven 1959; Hawkins 1961). Since then, this effect has been widely studied in the framework of T Tauri stars (Basri & Batalha 1990), but also for Be stars (see, e.g., Hubert-Delplace et al. 1982; Ballereau et al. 1995; Semaan et al. 2013). In order to avoid complex calculations of radiation transfer effects in circumstellar media, which in Be stars are still poorly mastered, we propose taking into account the veiling effect using heuristic expressions that describe the relation between line spectrum $F^{\text{obs}}(\lambda)$ and the continuum spectrum $F_c^{\text{obs}}(\lambda)$ emitted by the star+disk system, with $F^*(\lambda)$ and $F_c^*(\lambda)$ attributed to the stellar photosphere. The expression used in this work to carry out the correction of spectra by the veiling effect was discussed in Semaan et al. (2013). Some details regarding this correction and the values of the corresponding veiling factors are given in Appendix B.

## 5.2. Apparent astrophysical parameters

To derive the apparent astrophysical parameters, plane-parallel model atmospheres were computed for effective temperatures ranging from $T_{\text{eff}} = 8000-55000$ K and gravities $log g = 2.5$ to 4.5 dex. We computed the temperature structure of the atmospheres with the ATLAS9 FORTRAN program (Kurucz 1993; Castelli et al. 1997). For effective temperatures of $15000$ K $\le T_{\text{eff}} \le 27000$ K, the non-LTE level populations were calculated for the considered atoms using TLUSTY (Hubeny & Lanz 1995) by keeping the temperature and density distributions fixed. We adopted the solar chemical abundances given by Grevesse & Sauval (1998). For model atmospheres with effective temperatures lower than $T_{\text{eff}} = 15000$ K, full LTE was assumed, while those having $T_{\text{eff}} > 27500$ K are models from the OSTAR2002 non-LTE grid (Lanz & Hubeny 2003). The grid of fluxes we use during the fitting procedure was built with SYNSPEC. Other computing details of models are given in Frémat et al. (2005, 2006).

The apparent effective temperature and surface gravity are based on the fit of spectral regions not affected by the CE emission/absorption, which encompass wings of hydrogen lines as well as those of helium, carbon, oxygen, silicon, and Mg II lines in the spectral domain from 4000 to 4500 Å. The normalized spectra of the studied stars, $\Phi^*(\lambda)$, namely, those corrected for the veiling effect with Eq. (B.1), are Doppler-shifted according to the radial velocity. The fit of the veiling-corrected spectra is carried out by comparing them with a grid of synthetic spectra

using the robust least squares method employed in the MINUIT minimization package developed at CERN. This package was transcribed into a FORTRAN language and is named GIRFIT (Frémat et al. 2006). The synthetic spectra were convoluted with a Gaussian function for the required instrumental resolution.

In the most general form, the fitting procedure of spectra with MINUIT considers as an unknown the parametric set ($T_{\text{eff}}$, $log g$, $V\sin i$, $V_z$), where $V_z$ is the radial velocity not corrected for the Earth movement around the Sun. The fitting procedure is controlled with the $\chi^2$ test computed over zones in the 4000–4500 Å wavelength range where the line circumstellar emissions and shell absorptions and the diffuse interstellar absorption bands are carefully avoided. For three studied stars (stars Nos. 4, 13, and 14) we determined the $V\sin i$ using the Fourier transform (FT) method, which was then considered as a fixed quantity. The FT method was applied to the lines He I 4026 and 4472 and Mg II 4481 Å. These objects are the only ones in our sample for which the FT $V\sin i$ parameters issued from the three spectral lines are: on the one hand, mutually consistent within the average GIRFIT uncertainty and, on the other hand, consistent with the value obtained with the GIRFIT fitting method.

We note that in spite of the robust character of the MINUIT algorithm, fits of equal quality according to the $\chi^2$ test with slightly different values of the set ($T_{\text{eff}}$, $log g$, $V\sin i$) can be obtained; that is to say: the noise introduces degeneracies in the combination of different parameters, in particular, between $log g$ and $V\sin i$. This phenomenon is worth to be considered in some detail since $M_V$ and $log L/L_\odot$ rely sensitively on the value of $log g$, mainly in late B spectral types. We then introduced a refining protocol to determine the astrophysical parameters that is detailed in the following section.

## 5.3. Refined astrophysical parameters

In the refining operation of astrophysical parameters, we used the following BCD (Barbier-Chalonge-Divan) calibrations: $M_V = M_V(\lambda_1, D)$, (Zorec & Briot 1991), $T_{\text{eff}} = T_{\text{eff}}(\lambda_1, D)$, (Zorec et al. 2009), and $log g = log g(\lambda_1, D)$, (Zorec 1986). To this panoply of relations, we added the expression for the color excess $E(B-V) = [(B-V)_{\text{obs}} - (B-V)_*]$, where the intrinsic color index $(B-V)_*$ is a function of ($T_{\text{eff}}$, $log g$). These come from the synthetic UVB colors updated by Castelli F. in 2011 of the model atmospheres published in Castelli & Kurucz (2003). We employed the photometric quantities $V_{\text{obs}}$ and $(B-V)_{\text{obs}}$ derived from the (G, BP, RP) magnitudes of the Gaia photometric system (Riello et al. 2021), which represent a uniform photometric data base for our program stars. This uniformity is less obvious among the UVB photometric data proper.

We assumed that the corrections $\Delta T_{\text{eff}}$ and $\Delta log g$ meant to be applied to the first determination of ($T_{\text{eff}}$, $log g$) can be obtained from the deviations $\Delta M_V$ and $\Delta E(B-V)$ that best approach the condition $\lim(\Delta d/d) \to 0$ derived from Pogson's formula (Pogson 1856):

$$ \Delta d/d = 0.46[\Delta M_V + R_V \Delta E(B-V)]. \quad (1) $$

In Eq. (1), d is the stellar distance and $R_V = A_V/E(B-V) = 3.1$ because no other values of $R_V$ toward the direction of our program stars are known. In this treatment, the color excess $E(B-V)_T = (B-V)_{\text{obs}} - (B-V)(T_{\text{eff}}, \log g)$ has two components:

$$ E(B-V)_T = E(B-V)_{\text{ISM}} + E(B-V)_{\text{CD}}, \quad (2) $$
<table><caption>**Table 2.** Primary and secondary apparent astrophysical parameters.</caption><thead><tr><th rowspan="2">Star</th><th><i>T</i><sub>eff</sub></th><th>log <i>g</i></th><th>V sin <i>i</i></th><th>V<sub>Z</sub></th><th>log L/L<sub>☉</sub></th><th>M/M<sub>☉</sub></th><th>t/t<sub>MS</sub></th><th>R/R<sub>☉</sub></th><th>V<sub>c</sub></th></tr><tr><th>±σ<i>T</i><sub>eff</sub><br>(K)</th><th>±σlog <i>g</i><br>(dex)</th><th>±σ<i>V</i>sin <i>i</i><br>(km s<sup>-1</sup>)</th><th>±σ<i>V</i><sub>Z</sub><br>(km s<sup>-1</sup>)</th><th>±σlog L/L<sub>☉</sub></th><th>±σ<i>M</i>/M<sub>☉</sub></th><th>±σ<i>t</i>/t<sub>MS</sub></th><th>±σ<i>R</i>/R<sub>☉</sub></th><th>±σ<i>V</i><sub>c</sub><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>1</td><td>11450 ± 570</td><td>3.84 ± 0.12</td><td>120 ± 13</td><td>38.2 ± 2.9</td><td>2.318 ± 0.101</td><td>3.44 ± 0.17</td><td>0.88 ± 0.06</td><td>3.66 ± 0.35</td><td>348 ± 11</td></tr><tr><td>2</td><td>13670 ± 690</td><td>3.56 ± 0.12</td><td>230 ± 25</td><td>75.8 ± 3.2</td><td>3.076 ± 0.105</td><td>5.07 ± 0.30</td><td>1.03 ± 0.01</td><td>6.15 ± 0.62</td><td>324 ± 09</td></tr><tr><td>3</td><td>11180 ± 600</td><td>3.58 ± 0.13</td><td>220 ± 25</td><td>106.0 ± 3.2</td><td>2.584 ± 0.110</td><td>3.77 ± 0.22</td><td>1.02 ± 0.02</td><td>5.22 ± 0.55</td><td>304 ± 10</td></tr><tr><td>4</td><td>13610 ± 659</td><td>4.22 ± 0.11</td><td>370 ± 40</td><td>69.1 ± 11.9</td><td>2.264 ± 0.094</td><td>3.60 ± 0.17</td><td>0.30 ± 0.16</td><td>2.43 ± 0.21</td><td>442 ± 13</td></tr><tr><td>5</td><td>16030 ± 740</td><td>3.70 ± 0.11</td><td>230 ± 22</td><td>117.0 ± 8.5</td><td>3.288 ± 0.092</td><td>6.01 ± 0.29</td><td>0.97 ± 0.03</td><td>5.71 ± 0.49</td><td>368 ± 11</td></tr><tr><td>6</td><td>15590 ± 779</td><td>3.47 ± 0.11</td><td>250 ± 25</td><td>47.1 ± 6.6</td><td>3.513 ± 0.102</td><td>6.62 ± 0.42</td><td>1.02 ± 0.00</td><td>7.83 ± 0.75</td><td>330 ± 08</td></tr><tr><td>7</td><td>14030 ± 730</td><td>3.88 ± 0.13</td><td>270 ± 30</td><td>56.6 ± 2.7</td><td>2.752 ± 0.105</td><td>4.50 ± 0.24</td><td>0.85 ± 0.07</td><td>4.02 ± 0.41</td><td>381 ± 13</td></tr><tr><td>8</td><td>22830 ± 1187</td><td>3.76 ± 0.15</td><td>200 ± 39</td><td>45.9 ± 3.5</td><td>4.095 ± 0.127</td><td>10.70 ± 0.84</td><td>0.89 ± 0.06</td><td>7.13 ± 0.89</td><td>442 ± 16</td></tr><tr><td>9</td><td>14210 ± 720</td><td>3.37 ± 0.10</td><td>250 ± 27</td><td>57.2 ± 0.7</td><td>3.423 ± 0.097</td><td>6.21 ± 0.38</td><td>1.03 ± 0.00</td><td>8.50 ± 0.75</td><td>307 ± 07</td></tr><tr><td>10</td><td>10430 ± 550</td><td>3.59 ± 0.12</td><td>220 ± 24</td><td>76.0 ± 2.1</td><td>2.399 ± 0.103</td><td>3.38 ± 0.18</td><td>1.02 ± 0.02</td><td>4.85 ± 0.47</td><td>299 ± 09</td></tr><tr><td>11</td><td>14300 ± 700</td><td>4.00 ± 0.13</td><td>300 ± 33</td><td>61.5 ± 2.0</td><td>2.652 ± 0.105</td><td>4.33 ± 0.23</td><td>0.72 ± 0.11</td><td>3.45 ± 0.35</td><td>405 ± 14</td></tr><tr><td>12</td><td>11650 ± 580</td><td>3.58 ± 0.14</td><td>200 ± 20</td><td>46.8 ± 4.6</td><td>2.673 ± 0.113</td><td>3.98 ± 0.24</td><td>1.02 ± 0.02</td><td>5.33 ± 0.61</td><td>309 ± 11</td></tr><tr><td>13</td><td>13750 ± 690</td><td>3.82 ± 0.13</td><td>260 ± 29</td><td>78.8 ± 0.3</td><td>2.777 ± 0.107</td><td>4.50 ± 0.24</td><td>0.90 ± 0.06</td><td>4.31 ± 0.44</td><td>368 ± 12</td></tr><tr><td>14</td><td>22160 ± 1140</td><td>4.00 ± 0.15</td><td>325 ± 36</td><td>60.5 ± 4.3</td><td>3.714 ± 0.121</td><td>8.75 ± 0.58</td><td>0.65 ± 0.13</td><td>4.88 ± 0.59</td><td>484 ± 19</td></tr><tr><td>15</td><td>21830 ± 1500</td><td>3.58 ± 0.17</td><td>285 ± 31</td><td>57.2 ± 3.8</td><td>4.206 ± 0.156</td><td>11.11 ± 1.17</td><td>0.99 ± 0.04</td><td>8.87 ± 1.34</td><td>401 ± 17</td></tr><tr><td>16</td><td>10790 ± 550</td><td>3.97 ± 0.15</td><td>220 ± 24</td><td>68.9 ± 3.4</td><td>2.021 ± 0.116</td><td>2.95 ± 0.17</td><td>0.76 ± 0.11</td><td>2.93 ± 0.34</td><td>360 ± 14</td></tr><tr><td>17</td><td>12000 ± 600</td><td>4.04 ± 0.14</td><td>350 ± 38</td><td>57.4 ± 6.4</td><td>2.192 ± 0.112</td><td>3.30 ± 0.19</td><td>0.66 ± 0.13</td><td>2.88 ± 0.33</td><td>386 ± 14</td></tr><tr><td>18</td><td>13490 ± 670</td><td>3.57 ± 0.11</td><td>225 ± 24</td><td>86.8 ± 6.4</td><td>3.034 ± 0.097</td><td>4.94 ± 0.27</td><td>1.03 ± 0.01</td><td>6.03 ± 0.54</td><td>325 ± 09</td></tr><tr><td>19</td><td>17200 ± 1270</td><td>3.40 ± 0.16</td><td>100 ± 12</td><td>72.1 ± 6.9</td><td>3.851 ± 0.135</td><td>8.27 ± 0.76</td><td>1.02 ± 0.00</td><td>9.48 ± 1.38</td><td>336 ± 11</td></tr></tbody></table>

where $E(B-V)_{ISM}$ is the interstellar extinction proper and $E(B-V)_{CD}$ corresponds to the reddening produced on the Paschen continuum by the CD. We estimated $E(B-V)_{CD}$ with the empirical relation given as a function of the equivalent width $W_\alpha$ of the H$\alpha$ line emission component used in Gkouvelis et al. (2016), which is based on previous determinations by Raddi et al. (2013) and Dachs et al. (1988). Because the magnitude excess $\Delta V_{CD}$ in the V magnitude due to the CD, the reddening $E(B-V)_{CD}$ does not conform to the law $R_V = A_V/E(B-V) = 3.1$ that is otherwise valid for the interstellar medium, we must dissociate the circumstellar reddening from the interstellar one and apply it directly to $V_{obs}$. While $E(B-V)_{ISM}$ was estimated using Eq. (2) and the noted relation between $E(B-V)_{CD}$ and $W_\alpha$, the magnitude excess $\Delta V_{CD}$ was obtained from Eq. (8) in Zorec & Briot (1991).

When dealing with Eq. (1), knowledge of *d* is never required. We simply parameterized the ratio $\Delta d/d$ from -0.5 to +0.5, and for each $\Delta d/d$ we tested a series of values $\Delta M_V$ and $\Delta E(B-V)$, which (via the above BCD calibrations) enable us to calculate the corrections $\Delta T_{eff}$ and $\Delta \log g$. The chosen fit with GIRFIT then corresponds to the smallest possible value found of $|\Delta d/d|$ that leads to the best $\chi^2$ test. From this fit, we get not only the apparent astrophysical parameters ($T_{eff}$, $\log g$, $V \sin i$), but also $M_V$ and $E(B-V)_{ISM}$. In most cases, it is $|\langle\Delta d/d\rangle| \le 0.05$.

Specifically, BCD calibrations of astrophysical parameters are displayed in tabular forms and can be reverted to: $M_V = M_V(T_{eff}, \log g)$ and $E = E(B-V)(T_{eff}, \log g)$, from which the required relations readily follow:

$$ \begin{aligned} \Delta M_V &= (\partial M_V / \partial T_{eff}) \Delta T_{eff} + (\partial M_V / \partial \log g) \Delta \log g, \\ \Delta E &= (\partial E / \partial T_{eff}) \Delta T_{eff} + (\partial E / \partial \log g) \Delta \log g, \end{aligned} \quad (3) $$

needed to derive $\Delta T_{eff}$ and $\Delta \log g$ from $\Delta M_V$ and $\Delta E(B-V)$. The best fits of the observed spectra previously corrected for the veiling effect obtained with GIRFIT, subject to the noted iteration procedure, are shown in Appendices as Figs. C.1–C.4. The apparent astrophysical parameters derived are given in

Table 2. We refer to them as the “primary apparent astrophysical parameters”.

Due to the noted subtleties that have to be taken into account to calculate the spectroscopic distances of Be stars, it is worth specifying Pogson’s formula, as it is used in the present work:

$$ \left. \begin{aligned} \log d_{\text{spect.}} &= 0.2 (V_* - M_V) + 1.0 \\ V_* &= (V_{\text{obs}} - \Delta V_{\text{CD}}) - 3.1 E(B-V)_{\text{ISM}} \end{aligned} \right\}. \quad (4) $$

Once the final couple ($T_{eff}$, $\log g$) was adopted, we calculated the uncertainties associated to the estimated distances marked in tables of results, which depend only on the measurement errors of $V_{obs}$ and $(B-V)_{obs}$, $E(B-B)_{ISM}$ and $I_\alpha$. Each of these parameters is represented as $X + \epsilon_X$, where the error $\epsilon_X$ is simulated through a Monte Carlo procedure. The errors $\epsilon_X$ are assumed having a normal distribution characterized by the respective standard dispersion $\sigma_X$ reported in our tables. The reported values of $d_{spect.}$ in Table 3 are the averages $\langle d_{spect.} \rangle$ of 10⁴ estimates, and $\sigma_d$ represents their standard deviation. The same procedure was adopted also for $\langle d_{Gaia} \rangle = \langle 1/(\pi + \epsilon_\pi) \rangle$, where $\epsilon_\pi$ is considered to have a normal distribution with the standard deviation published for each Gaia DR3 parallax $\pi$ (Gaia Collaboration 2016, 2021; Lindegren et al. 2021). The errors associated to $M_V$ in Table 3 are also due to $V_{obs}$, $(B-V)_{obs}$, $E(B-B)_{ISM}$ and $I_\alpha$.

## 5.4. Secondary apparent astrophysical parameters

We use the term “secondary apparent astrophysical parameters” to apply to all those that tightly depend on the “primary apparent astrophysical parameters” and were not obtained from the fitting operation of spectra. Some of them also rely on models of stellar evolution. These are: the BCD parameters ($\lambda_1$, $D$), the bolometric luminosity, $\log L/L_\odot$, the stellar mass, $M/M_\odot$, equatorial radius, $R/R_\odot$, the equatorial linear critical velocity, $V_c$, and the fractional stellar age, $t/t_{MS}$ ($t$ is the age and $t_{MS}$ is the time that a non rotating star spends in the main sequence evolutionary phase).
<table><caption>**Table 3.** Indirect apparent quantities, interstellar reddening, and distances.</caption><thead><tr><th>Star</th><th>D<br>(dex)</th><th>λ<sub>1</sub><br>(Å)</th><th>MK</th><th>V<sub>obs</sub><br>(mag)</th><th>(B − V)<sub>obs</sub><br>(mag)</th><th>M<sub>V</sub> ± σ<sub>V</sub><br>(mag)</th><th>E(B − V)<sub>ISM</sub> ± σ<sub>EISM</sub><br>(mag)</th><th>d<sub>spect.</sub> ± σ<sub>d</sub><br>(pc)</th><th>d<sub>GAIA</sub> ± σ<sub>d</sub><br>(pc)</th></tr></thead><tbody><tr><td>1</td><td>0.415</td><td>45.8</td><td>B9 IV</td><td>13.259</td><td>0.442</td><td>−0.40 ± 0.12</td><td>0.525 ± 0.049</td><td>2556 ± 179</td><td>2540 ± 86</td></tr><tr><td>2</td><td>0.308</td><td>34.6</td><td>B7 III</td><td>13.526</td><td>0.452</td><td>−1.90 ± 0.12</td><td>0.554 ± 0.047</td><td>5425 ± 365</td><td>5175 ± 391</td></tr><tr><td>3</td><td>0.425</td><td>36.1</td><td>B9 III</td><td>15.698</td><td>0.782</td><td>−1.13 ± 0.13</td><td>0.841 ± 0.052</td><td>6858 ± 469</td><td>6572 ± 1936</td></tr><tr><td>4</td><td>0.310</td><td>61.6</td><td>B6-7 III</td><td>16.048</td><td>0.582</td><td>+0.17 ± 0.12</td><td>0.692 ± 0.049</td><td>5584 ± 377</td><td>5433 ± 1241</td></tr><tr><td>5</td><td>0.252</td><td>39.9</td><td>B5 III</td><td>16.233</td><td>0.966</td><td>−1.97 ± 0.14</td><td>1.119 ± 0.051</td><td>8768 ± 576</td><td>8774 ± 2741</td></tr><tr><td>6</td><td>0.242</td><td>34.9</td><td>B5 III</td><td>14.105</td><td>0.612</td><td>−2.65 ± 0.13</td><td>0.760 ± 0.049</td><td>7509 ± 498</td><td>7393 ± 1114</td></tr><tr><td>7</td><td>0.313</td><td>45.0</td><td>B6 IV</td><td>16.391</td><td>0.950</td><td>−0.97 ± 0.14</td><td>1.046 ± 0.052</td><td>6574 ± 439</td><td>6464 ± 1757</td></tr><tr><td>8</td><td>0.144</td><td>50.0</td><td>B2 III-IV</td><td>14.479</td><td>1.015</td><td>−3.09 ± 0.18</td><td>1.167 ± 0.052</td><td>6092 ± 403</td><td>5972 ± 698</td></tr><tr><td>9</td><td>0.268</td><td>32.5</td><td>B6 III</td><td>15.872</td><td>0.791</td><td>−2.61 ± 0.12</td><td>0.886 ± 0.051</td><td>13862 ± 925</td><td>13603 ± 4457</td></tr><tr><td>10</td><td>0.472</td><td>37.7</td><td>B9.5 III</td><td>12.059</td><td>0.301</td><td>−0.80 ± 0.12</td><td>0.363 ± 0.050</td><td>2227 ± 162</td><td>2184 ± 84</td></tr><tr><td>11</td><td>0.302</td><td>50.5</td><td>B6 III-IV</td><td>11.284</td><td>0.167</td><td>−0.66 ± 0.12</td><td>0.221 ± 0.046</td><td>1758 ± 124</td><td>1706 ± 60</td></tr><tr><td>12</td><td>0.399</td><td>35.7</td><td>B9 III</td><td>14.834</td><td>0.698</td><td>−1.26 ± 0.12</td><td>0.791 ± 0.051</td><td>5344 ± 366</td><td>5326 ± 751</td></tr><tr><td>13</td><td>0.322</td><td>43.0</td><td>B6-7 IV</td><td>11.418</td><td>0.174</td><td>−1.11 ± 0.14</td><td>0.274 ± 0.045</td><td>2135 ± 148</td><td>2073 ± 80</td></tr><tr><td>14</td><td>0.164</td><td>57.4</td><td>B2.5 V</td><td>11.411</td><td>0.344</td><td>−2.17 ± 0.14</td><td>0.510 ± 0.047</td><td>2489 ± 169</td><td>2510 ± 129</td></tr><tr><td>15</td><td>0.143</td><td>44.6</td><td>B2 III</td><td>15.855</td><td>1.211</td><td>−3.59 ± 0.19</td><td>1.434 ± 0.053</td><td>10000 ± 648</td><td>9476 ± 2950</td></tr><tr><td>16</td><td>0.445</td><td>55.7</td><td>B9 IV</td><td>15.304</td><td>0.641</td><td>+0.21 ± 0.13</td><td>0.698 ± 0.052</td><td>3857 ± 276</td><td>3829 ± 400</td></tr><tr><td>17</td><td>0.379</td><td>54.5</td><td>B7-8 V</td><td>16.454</td><td>1.444</td><td>+0.02 ± 0.12</td><td>1.521 ± 0.055</td><td>2185 ± 142</td><td>2260 ± 186</td></tr><tr><td>18</td><td>0.316</td><td>34.8</td><td>B7 III</td><td>15.929</td><td>1.057</td><td>−1.81 ± 0.12</td><td>1.189 ± 0.052</td><td>6482 ± 424</td><td>6255 ± 1463</td></tr><tr><td>19</td><td>0.202</td><td>36.7</td><td>B4 III</td><td>15.557</td><td>0.910</td><td>−3.21 ± 0.14</td><td>1.054 ± 0.051</td><td>12481 ± 819</td><td>13504 ± 4243</td></tr></tbody></table>

The BCD quantities ($\lambda_1$, $D$) are determined from the adopted ($T_{eff}$, log $g$) parameters, which in turn lead to the BCD spectral classification of the studied Be stars (Zorec & Briot 1991; Zorec et al. 2009). This set of secondary apparent parameters is given in Table 3.

Adopting the bolometric correction BC($T_{eff}$) by Flower (1996), to avoid possible conflicts with the zero point of bolometric corrections, we determined the bolometric “apparent” absolute magnitude $M_{bol}$ following the prescription from Torres (2010):

$$M_{bol} = M_{bol}^{\odot} + [M_{V} - M_{V}^{\odot}] + [BC(T_{eff}) - BC(\odot)]. \quad (5)$$

From this, we can derive the bolometric luminosity, log $L/L_\odot$. Once log $L/L_\odot$ and $T_{eff}/T_{eff}^{\odot}$ are determined, we estimate the stellar radius $R/R_\odot$. Finally, by interpolating in the grids of stellar evolutionary models without rotation given by Georgy et al. (2013) for solar metallicity $Z = 0.014$, we obtain the stellar masses $M/M_\odot$ and its fractional ages $t/t_{MS}$. All these secondary apparent astrophysical parameters are reported in Table 2 with their uncertainties.

To complete the presentation of the program objects, we show their HR diagram in Fig. 2, where we have included the instability zones calculated by Walczak et al. (2015) for non rotating objects using the OPLIB new Los Alamos opacities. Similar calculations can be found in Miglio et al. (2007), Burssens et al. (2020) for other opacities, and rotating stars in Saio et al. (2017) and Szewczuk & Daszyńska-Daszkiewicz (2017), where the widening of the instability regions due to stellar rotation can be seen.

## 5.5. Astrophysical parameters corrected for effects carried by the rapid rotation

We consider that the stars are rigid rotators characterized by a uniform surface rate of angular velocity ratio $\Omega/\Omega_c$, where $\Omega_c$ is the surface angular velocity at critical rotation. The apparent astrophysical parameters obtained in Sect. 5.2 can then be given

according to the following formal expression:

$$P_{\text{app}} = P_{\text{pncr}} F_{P}(M, t, \Omega/\Omega_c, i), \quad (6)$$

where $P$ stands for parameters as effective temperature, $T_{eff}$, surface gravity, $g$, bolometric luminosity, $L$, and projected rotational velocity, $V\sin i$. On the left-hand side of Eq. (6), these quantities are considered “apparent”, that is, they are issued from models of stellar atmospheres and evolution without rotation that fit the observed spectral domain. On the right side of Eq. (6), the parameters $P$ correspond to the parent-non-rotating-counterparts (pncr; Frémat et al. 2005), which are multiplied by functions that represent the effects carried by a rapid rotation in an object of actual mass, $M$, and age, $t$, whose surface angular velocity is $\Omega$ and its rotational axis is seen according to an aspect angle, $i$.

The functions $F_P(M, t, \Omega, i)$ in Eq. (6) were calculated with FASTROT (Frémat et al. 2005) as follows. For a series of stellar masses, ages, and rotational rates, we calculated 2D models of stellar structure assuming rigid rotation over the whole star (Zorec et al. 2011; Zorec & Royer 2012). The geometrical deformations of stars induced by the rapid rotation produce a surface gravity that depends on the latitude and a non uniform distribution of the effective temperature, which is calculated by the gravity darkening formula given by Espinosa Lara & Rieutord (2011). Using local plane-parallel model atmospheres, we recomposed the apparent stellar spectra in the $\lambda\lambda4000-4500$ Å wavelength region for objects of several masses, ages, angular velocity rates $\Omega/\Omega_c$, and inclination angles, $i$. The system of equations in Eq. (6) is given in a tabular form whose solutions for the unknown quantities ($M, t, \Omega/\Omega_c, i$) are obtained by iteration. Since only $T_{eff}$, log $g$ and $V\sin i$ as entry apparent parameters can be considered independent, we solved these equations for a series of imposed values of $\Omega/\Omega_c$. Due to gravitational darkening, the observed (apparent) $V\sin i$ are underestimated (Stoeckley 1968; Townsend et al. 2004; Frémat et al. 2005). The amount of this underestimation is automatically taken into account in the solution of Eq. (6). The evolutionary tracks of rotating stars used were calculated by
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** HR diagram of the studied stars. Evolutionary tracks are from Georgy et al. (2013) for $\Omega/\Omega_c = 0.0$ (blue) and for rotation rates at the ZAMS $\Omega/\Omega_c = 0.95$ (red). Blue points: stars with their apparent parameters; red points: Stars with their pncr parameters obtained for $\Omega/\Omega_c = 0.95$. Instability strips of $\beta$ Ceph- (---) and SPB-type (---) pulsators adapted from Walczak et al. (2015) for OPLIB opacities and metallicity of $Z = 0.015$.</figcaption></figure>

Georgy et al. (2013) and are given in terms of effective temperatures and bolometric luminosity averaged over the surface of the geometrically deformed rotating objects. This imposes that at each iteration step, the corresponding averaged quantities over the stellar surface are calculated from the running pncr parameters. Moreover, at each iteration step we need to determine the rotational velocity of stars in the ZAMS to identify the right evolutionary track needed to interpolate the stellar mass and age at the running stellar age, which is consistent with the imposed rotational rates $\Omega/\Omega_c$. Once a cycle of iterations is over, we calculate the rotational frequency, $v_r$, given by:

$$v_r = 0.02 \left[ \frac{V_{eq}(\Omega/\Omega_c, M, t)}{R_{eq}(\Omega/\Omega_c, M, t)} \right] \text{cycles/day}, \quad (7)$$

where $V_{eq}$ is the equatorial linear velocity given in km s$^{-1}$ and $R_{eq}$ is the equatorial radius in solar units of the star distorted by its rotation.

Each entry parameter in Eq. (6) is considered with its uncertainty assumed having a normal (Gaussian) distribution. We produced $10^4$ Monte Carlo trials of all entry parameters and each time we sought a new solution for the system using Eq. (6). Sometimes, the set of entry parameters ($T_{eff} \pm \epsilon_{T_{eff}}$, $\log g \pm \epsilon_{\log g}$, $\log L/L_{\odot} \pm \epsilon_{\log L}$, $V\sin i \pm \epsilon_{V\sin i}$) may not correspond to realistic objects. In such cases, the functions $F_P(M, t, \Omega, i)$ cannot produce a reliable solution and the trial is considered as having no solution. Although the uncertainties of the entry parameters are assumed to have normal distributions, the distributions of the solution quantities do not always have normal distributions. So, the averages of parameters with non symmetrical distribution do not coincide with the mode of the distribution (mode = parameter corresponding to the maximum of the distribution). In order

to illustrate the obtained differences in the estimated parameters, we produced two tables of results: one of them gives the average values of solutions with their classical standard $1\sigma$ deviations, while the other is for modes, where the “sigmas” correspond to the uncertainties related to the identification of the modes. Similar solutions to the system Eq. (6) have been previously done in Floquet et al. (2000), Neiner et al. (2003), Frémat et al. (2005), Zorec et al. (2005), Zorec et al. (2016), Vinicius et al. (2006), Martayan et al. (2006, 2007), Semaan et al. (2013). In Appendix D, Tables D.1 and D.2 give the averages and modes of distributions of astrophysical parameters and rotational frequencies corrected for rotational effects assuming $\Omega/\Omega_c = 0.95$, respectively. These tables give: the pncr effective temperature, $T_{eff}$, $\log g$, true $V\sin i$, pncr bolometric luminosity in solar units, $\log L/L_{\odot}$, mass $M/M_{\odot}$, critical linear equatorial velocity, $V_c$, estimated inclination angle, $i$, true fractional age, $t/t_{MS}$, the true stellar age, $t$, and the rotational frequency, $v_r$, in cycles/day. All these quantities are given with their respective $1\sigma$ dispersion calculated from the obtained distribution of solutions. Similar solutions obtained for imposed values $\Omega/\Omega_c = 0.8, 0.9, 0.999$, and 1.0 are given in Tables D.1 and D.2, respectively, available at the CDS.

# 6. Some empirical relations concerning the CD in the studied Be stars

The close relation between the spectroscopic distances and those measured by Gaia, provide us with some confidence on the reliability of the inferred astrophysical parameters. We can then try to benefit from some information on the properties of CDs of the studied stars, drawn from the emission in spectral lines. Thanks to the derived astrophysical parameters, good estimates of the neat emission line profiles and of the amount of emission filling up the photospheric component of the here observed Balmer lines $H\alpha$, $H\gamma$, and $H\delta$ can be obtained, at least when it is possible according to the S/N of the spectra.

## 6.1. Total emission in the observed Balmer lines

In cases when the emission in a line clearly exceeds the continuum over its entire wavelength range, we can simply add to this latter the amount of emission that fills up the photospheric absorption component derived from the stellar astrophysical parameters. If this filling up is partial, here the apparent photospheric absorption profile underlying the emission is called the pseudo-photospheric absorption. We estimated the pseudo-photospheric absorption and add to this last the difference between the pseudo-photospheric and the photospheric proper absorption line profile. The pseudo-photospheric (pph) absorption line profile is easily fitted using the relation $\psi(\lambda) = \exp\{-1/[a + b \times (\lambda - \lambda_o)^c]\}$ introduced by Ballereau et al. (1995), Chauville et al. (2001), where the constants $a, b$, and $c$ are determined using only three points in the line profile. In Fig. 3, the red lines show the total $H\alpha$ line-emission components reduced to the continuum intensity level $I_{\lambda}/I_c = 1$ (violet dotted line). Thus, the total emission profile (red line) makes up the observed profile proper, to which the residual intensity between the photospheric profile and the pseudo-photospheric absorption profile has been added. The blue dashed lines correspond to the photospheric absorption line profile broadened according to the measured $V\sin i$ parameter. The green dashed lines are the pseudo-photospheric absorption line profiles. The equivalent width of the total $H\alpha$ emission corresponds then to the area below the red profile and the $I_{\lambda}/I_c = 1$ line, that is,
<figure><img src="image_5.png" /><figcaption>Fig. 3. Smoothed Hα line-emission components (red profiles). The emission residual intensities are relative to the photospheric line calculated according to the astrophysical parameters ($T_{eff}$, log g, Vsin i) of each object. The figure shows the photospheric component broadened by rotation (blue profiles). It includes the empirical pseudo-photospheric component (green profiles, see Sect. 6.1) in cases of partial emission feeling up of the absorption photospheric line. For the stars 103000272, 110681176, 110747131, and 110827583, the fit of the observed pseudo-photospheric absorption component closely reproduces the actual photospheric line.</figcaption></figure>

$W_{H\alpha} = W_{H\alpha}^{em} + W_{H\alpha}^{ph} - W_{H\alpha}^{pph}$. In Table 4, we give the equivalent widths $W_{H\alpha}$, $W_{H\alpha}^{em}$, $W_{H\alpha}^{ph}$, and $W_{H\alpha}^{pph}$.

The line emission profiles in the Hγ and Hδ lines are determined through the residual intensities between the emission proper and the photospheric line profile. We also assume that the emission in these lines is weak enough to take the pseudo-photospheric profile as the genuine photospheric absorption. In Fig. 4, we reproduce the emission line profiles in Hγ and Hδ lines. Some emission in these lines can also be present in other stars, but they are too noisy to carried out a clear extraction. We attempted to determine the equivalent width for all detectable emissions; they are given in Tables 4 and B.1 and were previously used in Sect. 5.1 to calculate the veiling factor.

## 6.2. Other measurements on the emission line profiles

An estimate of the total extent of the CD line emission formation region is frequently done using a relation which assumes that the disk is optically thin and that it has a rotation law that is dependent of the distance $R/R_o$ ($R_o$ is radius of the base of the disk) from the central star (Huang 1972). Following the reasoning in Arcos et al. (2017) we can write Huang’s radius as

$$ (R/R_o)_H = (2.7V\sin i/\Delta_p)^{1/j}, \qquad (8) $$

where $\Delta_p$ is the distance in km s⁻¹ between the emission peaks, and j = 1/2 if the rotation is Keplerian or j = 1 if the disk conserves the angular momentum. The factor of 2.7 in Eq. (8)
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Smoothed Hγ (red) and Hδ (blues) line-emission components. The emission residual intensities are relative to each photospheric line. Emissions in these Balmer lines are seen also in some other stars, but they are too noisy and not exploitable.</figcaption></figure>

<table><caption>**Table 4.** Equivalent widths and normalized fluxes of the emission components in the Hα and Hδ lines.</caption><thead><tr><th>Star</th><th>W<sup>em</sup><sub>Hα</sub><br>(Å)</th><th>W<sup>ph</sup><sub>Hα</sub><br>(Å)</th><th>W<sup>pph</sup><sub>Hα</sub><br>(Å)</th><th>I<sub>Hα</sub><br>(Å)</th><th>W<sup>em</sup><sub>Hδ</sub><br>(Å)</th><th>I<sub>Hδ</sub><br>(Å)</th></tr></thead><tbody><tr><td>1</td><td>1.8</td><td>9.0</td><td>8.9</td><td>0.6</td><td>0.0</td><td>0.0</td></tr><tr><td>2</td><td>10.8</td><td>5.7</td><td>0.0</td><td>7.7</td><td>0.0</td><td>0.0</td></tr><tr><td>3</td><td>7.3</td><td>8.0</td><td>3.7</td><td>3.8</td><td>1.2</td><td>0.4</td></tr><tr><td>4</td><td>4.3</td><td>8.2</td><td>4.6</td><td>3.6</td><td>0.9</td><td>0.3</td></tr><tr><td>5</td><td>2.0</td><td>5.1</td><td>3.7</td><td>2.0</td><td>0.7</td><td>0.3</td></tr><tr><td>6</td><td>3.5</td><td>4.8</td><td>0.0</td><td>4.7</td><td>0.0</td><td>0.0</td></tr><tr><td>7</td><td>14.1</td><td>6.6</td><td>0.0</td><td>10.0</td><td>0.5</td><td>0.2</td></tr><tr><td>8</td><td>34.8</td><td>3.4</td><td>0.0</td><td>39.0</td><td>0.7</td><td>0.8</td></tr><tr><td>9</td><td>20.8</td><td>5.2</td><td>0.0</td><td>12.8</td><td>1.2</td><td>0.5</td></tr><tr><td>10</td><td>0.8</td><td>9.4</td><td>0.0</td><td>3.0</td><td>0.0</td><td>0.0</td></tr><tr><td>11</td><td>35.8</td><td>6.9</td><td>0.0</td><td>21.2</td><td>0.7</td><td>0.3</td></tr><tr><td>12</td><td>1.4</td><td>7.5</td><td>0.0</td><td>3.2</td><td>0.0</td><td>0.0</td></tr><tr><td>13</td><td>14.3</td><td>6.6</td><td>2.9</td><td>8.5</td><td>0.1</td><td>0.0</td></tr><tr><td>14</td><td>24.5</td><td>4.0</td><td>0.0</td><td>27.8</td><td>0.7</td><td>0.70</td></tr><tr><td>15</td><td>0.0:</td><td>3.2</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>16</td><td>1.9</td><td>10.7</td><td>0.0</td><td>3.9</td><td>0.0</td><td>0.0</td></tr><tr><td>17</td><td>5.1</td><td>9.0</td><td>5.5</td><td>3.2</td><td>1.3</td><td>0.4</td></tr><tr><td>18</td><td>0.0:</td><td>5.9</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>19</td><td>15.3</td><td>4.1</td><td>0.0</td><td>12.6</td><td>1.3</td><td>1.7</td></tr></tbody></table>

differs from that used in Arcos et al. (2017), because stellar deformations were calculated following Frémat et al. (2005) for $\Omega/\Omega_c = 0.88$. In Table 5, we give the measured emission peak separations and the disks extents $R/R_o$ determined assuming Keplerian rotation ($j = 0.5$). To avoid confusion, we note that Huang’s radius ($R/R_o$)$_H$ in Table 5 is calculated with $\Delta_p$, while ($R/R_o$)$_{hm}$ represents a kind of Huang’s CD extent determined with the full-widths-at-half-maximum (FWHM), $\Delta_{hm}$, of the respective Balmer emission lines. The radius ($R/R_o$)$_H$ is meant to represent the extent of the entire region in the CD where the emission is raised, if the CD is optically thin, while ($R/R_o$)$_{hm}$

represents just a part of this line-forming region. The width of the base of emission line should in principle correspond to ($R/R_o$)$_{base} = 1$. However, the wings of emission lines are affected by other broadening mechanisms apart from rotation, so that the relation in Eq. (8) cannot identify well defined distances in the CD when using the FWHM or the width at the base of emission lines. Relations with the $V\sin i$ of the peak separations $\Delta_p$ and line widths at half intensity $\Delta_{hm}$ of the emission components in the Hα, Hγ and Hδ lines are shown in Fig. 5.

### 6.3. Near-IR circumstellar color excess and the Balmer line emissions

To complete the presentation of observational characteristics of the program Be stars, we include in this section a short discussion of their near-IR photometric data. These allow us, on the one hand, to situate them in the color-color diagram and show they do not belong to the B[e] class of stars known as Herbig HAc/Be, and on the other hand, to obtain some relations commonly used to present the Be phenomenon (Ashok et al. 1984; Dachs et al. 1988; Hanuschik 1989; Kastner & Mazzali 1989). However, it is not in the scope of the present work to extract information on the CDs from the photometric data.

Be stars present a thermal IR flux excess mainly due to bound-free, free-free transitions, and electron scattered radiation in their CD. To characterize the energy distribution in the visible domain, we could use the $(B-V)$ color index from the UBV Johnson-Cousins photometric system derived from the $(G, G_{BP}, G_{RP})$ Gaia photometry using the published transformation relations into the UBV photometric system (Riello et al. 2021) in the way that has widely been studied by Moujtahid et al. (1998, 1999). However, this color index is heavily marred by the interstellar reddening so that the genuine color excess due to the CD is marred by the uncertainties of the $E(B-V)_{ISM}$ color excess determinations. We then used the near-IR photometry in the J, H, and K bands to try to characterize the color excess produced by the Be star CD. To this end, the 2MASS photometry (Cutri et al. 2003) was used.
<table><caption>**Table 5.** Separation of emission peaks $\Delta_p$, line emission FWHM, $\Delta_{hm}$, and Huang's radius of the CD according to $\Delta_p$ and $\Delta_{hm}$.</caption><thead><tr><th rowspan="2">Star</th><th colspan="4">Hα</th><th colspan="4">Hγ</th><th colspan="4">Hδ</th></tr><tr><th>Δ<sub>p</sub><br>(km s<sup>-1</sup>)</th><th>Δ<sub>hm</sub><br>(R/R<sub>o</sub>)<sub>H</sub><br>(R/R<sub>o</sub>)<sub>hm</sub></th><th></th><th></th><th>Δ<sub>p</sub><br>(km s<sup>-1</sup>)</th><th>Δ<sub>hm</sub><br>(R/R<sub>o</sub>)<sub>H</sub><br>(R/R<sub>o</sub>)<sub>hm</sub></th><th></th><th></th><th>Δ<sub>p</sub><br>(km s<sup>-1</sup>)</th><th>Δ<sub>hm</sub><br>(R/R<sub>o</sub>)<sub>H</sub><br>(R/R<sub>o</sub>)<sub>hm</sub></th><th></th><th></th></tr></thead><tbody><tr><td>1</td><td>-</td><td>180</td><td>-</td><td>4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>2</td><td>-</td><td>230</td><td>-</td><td>9</td><td>103</td><td>255</td><td>45</td><td>7</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>3</td><td>111</td><td>312</td><td>33</td><td>4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>4</td><td>296</td><td>523</td><td>13</td><td>4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>5</td><td>355</td><td>571</td><td>4</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>6</td><td>209</td><td>412</td><td>13</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>7</td><td>126</td><td>332</td><td>40</td><td>6</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>8</td><td>158</td><td>391</td><td>13</td><td>2</td><td>201</td><td>339</td><td>8</td><td>3</td><td>205</td><td>305</td><td>8</td><td>4</td><td>4</td></tr><tr><td>9</td><td>225</td><td>483</td><td>11</td><td>2</td><td>225</td><td>473</td><td>11</td><td>3</td><td>241</td><td>512</td><td>10</td><td>2</td><td>2</td></tr><tr><td>10</td><td>253</td><td>392</td><td>6</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>11</td><td>122</td><td>371</td><td>52</td><td>6</td><td>182</td><td>408</td><td>24</td><td>5</td><td>225</td><td>491</td><td>15</td><td>3</td><td>3</td></tr><tr><td>12</td><td>71</td><td>214</td><td>68</td><td>8</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>13</td><td>95</td><td>294</td><td>67</td><td>7</td><td>182</td><td>335</td><td>18</td><td>5</td><td>213</td><td>351</td><td>13</td><td>5</td><td>5</td></tr><tr><td>14</td><td>142</td><td>403</td><td>42</td><td>5</td><td>221</td><td>449</td><td>17</td><td>4</td><td>233</td><td>504</td><td>16</td><td>3</td><td>3</td></tr><tr><td>16</td><td>170</td><td>334</td><td>14</td><td>4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>17</td><td>154</td><td>369</td><td>40</td><td>7</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>19</td><td>-</td><td>199</td><td>-</td><td>11</td><td>40</td><td>201</td><td>60</td><td>2</td><td>71</td><td>20</td><td>18</td><td>2</td><td>2</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Relations with the $V\sin i$ of the peak separations $\Delta_p$ and line widths at half intensity $\Delta_{hm}$ of the emission components in the Hα (red), Hγ (blue) and Hδ (green) lines. A black circle surrounds the binary star N° 5.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>Fig. 6. Dereddened JHK color-color diagram of the program Be stars. The JHK magnitudes are from 2MASS. In this diagram are also shown the model sequences of the intrinsic colors corresponding to log g = 1.0, 3.0, and 5.0 dex, and the corner occupied by the bluest Herbig HAE/Be stars. A black circle surrounds the binary star N° 5.</figcaption></figure>

The observed (J - H)obs and (H - J)obs colors were corrected from interstellar extinction using the relations between the color excesses E(V - J) = 2.23 E(B - V), E(V - H) = 2.55 E(B - V), and E(V - K) = 2.76 E(B - V) and the reddening relation, A_V = 3.1E(B - V), established by J. S. Mathys (1999, priv. comm.) and published in Cox (2000), Kinman & Castelli (2002). The determination of the color excess E(B - V) was detailed in Sect. 5.4. The JHK dereddened color-color diagram of the studied stars is shown in Fig. 6 (red stars), where we added the model sequences of intrinsic colors for log g = 1.0 to 5.0 dex by Castelli & Kurucz (2003, updated in 2011), and the region occupied by the bluest Herbig HAE/Be stars (Hernández et al. 2005) to show that our program stars were not mistakenly chosen among Herbig HAE/Be objects. The stars numbered 17 and 19 have much larger JHK color uncertainties than the remaining ones, which might explain their extreme blue (H - K)_o colors.

The color excesses proper in the IR that characterize the flux excess produced by the CD are defined as

$$ \begin{cases} \Delta(J - H)_{\text{CD}} = (J - H)_o - (J - H)_* \\ \Delta(H - K)_{\text{CD}} = (H - K)_o - (H - K)_* \end{cases} \quad (9) $$

where (J - H)* and (H - K)* are the respective intrinsic colors of stars that underlay the CD. The intrinsic colors are inferred from the synthetic colors indices calculated with LTE models of stellar atmospheres by Castelli & Kurucz (2003), using the apparent (T<sub>eff</sub>, log g) parameters. Also, (J - H)_o and (H - K)_o are the observed color indices corrected for the ISM extinction.

Having at our disposal observations in the Hα line – as well as, for some stars, also in the Hγ and Hδ reliable line emissions – we can look for a correlation between the flux IR color excess and the intensity of the emission component in the Balmer line emission components. Because the IR flux excess and the Balmer line emissions are produced in rather different regions of the same CD, such correlations can be of interest in attempting to control the consistency of physical inputs meant to describe the structure of Be star disks and their formation mechanisms. Marginally, we note that the sought correlations are marred by uncertainties due mainly to the non-simultaneity of observations carried of the IR colors and of the Balmer lines, as well as of the

physical origins. In fact, since the perturbations produced by the CD on the continuum spectrum and on the line emissions are not produced in the same layers, there may be line emission intensities produced rather far from the star (roughly 2 R<sub>o</sub> ≲ R ≲ 10 R<sub>o</sub>) that correspond to positive or negative flux excesses coming from disk layers near the star, generally R ≲ 2 R<sub>o</sub>. In spite of these shortcoming, we can attempt to establish some relations with the data at disposal in the way were previously attempted in Ballereau et al. (1995) between the flux excess ΔV in the magnitude V and the Hγ line emission, which revealed a quite reliable and useful correlation, which is described in Sect 5.1 as a way of correcting the observed spectra from the veiling effect in the λλ 4000–4500 Å spectral range.

The emission intensities at Hα, Hγ, and Hδ lines are given by:

$$ I_{H\alpha,H\gamma,H\delta} = W_{H\alpha,H\gamma,H\delta}^{\text{em}} \times \left[ \frac{F_{H\alpha,H\gamma,H\delta}^{\text{c}}}{F_{H\alpha,H\gamma,H\delta}^{\text{c}}(22\,500, 4.0)} \right], \quad (10) $$

where W<sub>Hα,Hγ,Hδ</sub><sup>em</sup> are the equivalent widths of the emission components in the Hα, Hγ, and Hδ lines, F<sub>Hα,Hγ,Hδ</sub><sup>c</sup> are the fluxes of the continuum spectrum at the center of the lines, and F<sub>Hα,Hγ,Hδ</sub><sup>c</sup>(22 500, 4.0) are the normalizing fluxes given by a model with T<sub>eff</sub> = 22 500 K and log g = 4.0 (roughly a B2V type star). In Table 4, we give the equivalent widths W<sub>Hα,Hγ</sub><sup>em</sup> of the emission component in the Hα and Hγ lines, and the respective intensities I<sub>Hα</sub> and I<sub>Hγ</sub>. In Fig. 7, we show the trends determined by I<sub>Hα</sub>, I<sub>Hγ</sub>, and I<sub>Hδ</sub> against the color index excess Δ(J - H)<sub>CD</sub>, which are rather well defined, while the correlations of the line emission intensities with Δ(H - K)<sub>CD</sub>, not shown here, look much more scattered. In Table 6, we give the observed color indices (J - H)<sub>obs</sub> and (H - K)<sub>obs</sub>, the adopted stellar intrinsic colors (J - H)* and (H - K)*, and the color index excess Δ(J - H)<sub>CD</sub>.

# 7. Physical characteristics of CDs derived from the interpretation of the Balmer line emission components

In this section we explore the properties of the CDs derived from the spectroscopic characteristics of the Hα, Hγ, and Hδ emission lines. We produce a simple representation of the CD that may contrast in simplicity with more detailed calculations as those by Carciofi & Bjorkman (2006, 2008), Sigut & Jones (2007), Sigut et al. (2009), McGill et al. (2013), Catanzaro (2013), Kurfürst et al. (2018). The simplifications we use are based on first physical principles and are introduced in as physically consistent a way as possible.

## 7.1. Assumptions about the geometry and temperature of the CD

To describe the geometrical characteristics of the CD, we use two reference systems centered on the star. In the star system (X, Y, Z), the Z-axis contains the rotation axis, while (X, Y) are in the stellar equatorial plane with X on the back sky plane. The system (x, y, z) has the z-axis directed towards the observer, and x coincides with X. The angle between Z and z is the inclination angle, i, of the star-disk system. The radius vector, R, is measured on the (X, Y) plane. Between Y and R, there is the azimuth angle θ.

We take the CD driven by a quasi-Keplerian speed of rotation V<sub>K</sub>(R) = V<sup>0</sup><sub>K</sub>(R<sub>o</sub>/R)<sup>γ/2</sup>. We include an expansion (or contraction) velocity parameter V<sub>exp</sub> to account for the emission line
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** Relations between the normalized line emission intensities $I_{H\alpha}$ (red), $I_{H\gamma}$ (blue), and $I_{H\delta}$ (green) and the color dereddened color excess $\Delta(J-K)_{CD}$ mag due to the CD. A black circle surrounds the binary star $N^5$.</figcaption></figure>

<table><caption>**Table 6.** JHK color indices of the studied stars.</caption><thead><tr><th>Star</th><th>(J-H)obs</th><th>(H-K)obs</th><th>(J-H)<sub>*</sub></th><th>(H-K)<sub>*</sub></th><th>Δ(J-H)<sub>CD</sub></th><th>Δ(H-K)<sub>CD</sub></th></tr></thead><tbody><tr><td>1</td><td>0.143</td><td>0.102</td><td>-0.026</td><td>-0.028</td><td>0.001</td><td>0.019</td></tr><tr><td>2</td><td>0.154</td><td>0.119</td><td>-0.044</td><td>-0.045</td><td>0.021</td><td>0.047</td></tr><tr><td>3</td><td>0.264</td><td>0.137</td><td>-0.022</td><td>-0.027</td><td>0.018</td><td>-0.012</td></tr><tr><td>4</td><td>0.198</td><td>0.228</td><td>-0.045</td><td>-0.043</td><td>0.022</td><td>0.126</td></tr><tr><td>5</td><td>0.482</td><td>0.223</td><td>-0.060</td><td>-0.057</td><td>0.184</td><td>0.045</td></tr><tr><td>6</td><td>0.192</td><td>0.150</td><td>-0.056</td><td>-0.055</td><td>0.005</td><td>0.046</td></tr><tr><td>7</td><td>0.600</td><td>0.156</td><td>-0.048</td><td>-0.046</td><td>0.313</td><td>-0.018</td></tr><tr><td>8</td><td>0.549</td><td>0.247</td><td>-0.093</td><td>-0.086</td><td>0.270</td><td>0.088</td></tr><tr><td>9</td><td>0.512</td><td>0.128</td><td>-0.047</td><td>-0.048</td><td>0.275</td><td>-0.010</td></tr><tr><td>10</td><td>0.090</td><td>0.077</td><td>-0.012</td><td>-0.019</td><td>-0.014</td><td>0.020</td></tr><tr><td>11</td><td>0.103</td><td>0.120</td><td>-0.050</td><td>-0.047</td><td>0.082</td><td>0.121</td></tr><tr><td>12</td><td>0.454</td><td>0.112</td><td>-0.027</td><td>-0.031</td><td>0.228</td><td>-0.023</td></tr><tr><td>13</td><td>0.082</td><td>0.037</td><td>-0.046</td><td>-0.045</td><td>0.040</td><td>0.024</td></tr><tr><td>14</td><td>0.153</td><td>0.185</td><td>-0.090</td><td>-0.083</td><td>0.080</td><td>0.161</td></tr><tr><td>15</td><td>0.420</td><td>0.296</td><td>-0.089</td><td>-0.082</td><td>0.050</td><td>0.078</td></tr><tr><td>16</td><td>0.138</td><td>0.164</td><td>-0.018</td><td>-0.021</td><td>-0.067</td><td>0.038</td></tr><tr><td>17</td><td>0.491</td><td>0.245</td><td>-0.032</td><td>-0.032</td><td>0.037</td><td>-0.042</td></tr><tr><td>18</td><td>0.345</td><td>0.164</td><td>-0.043</td><td>-0.044</td><td>0.008</td><td>-0.042</td></tr><tr><td>19</td><td>0.409</td><td>0.226</td><td>-0.066</td><td>-0.063</td><td>0.138</td><td>0.067</td></tr></tbody></table>

**Notes.** Stars numbered 7, 9, and 12 have color indices obtained from the (G, BP, RP) magnitudes of the Gaia photometry.

profile asymmetry. The integration of the hydrostatic equilibrium equation for the circumstellar environment leads to the axisymmetric density profile of the disk given in Kurfürst et al. (2018), of which we adopt the simplified expression valid for $Z/R \ll 1$:

$$ \rho(R, z) = \rho_o D(R) \exp \left\{ - \left[ Z/\sqrt{2h(R)} \right]^2 \right\}, \quad (11) $$

which still describes the main geometrical characteristics of disks and enables easier mathematical handling. In Eq. (11), $\rho_o$ is the base density of the disk at the stellar equator, $D(R)$ is a function meant to describe the CD density distribution in the equatorial plane, and $h(R) = [V_s/V_K]R$ is the Z-scale height of the disk with $V_s$ as the sound speed and $V_K$ as its quasi-Keplerian rotation. Explicitly

$$ \begin{aligned} h(R) &= k \left[ (R_o/R_\odot)^\gamma (T_{\text{eff}}/T_\odot)/(M/M_\odot) \right]^{1/2} (R/R_o)^\beta \\ k &= \left[ R/(\sqrt{2}\mu G) (R_\odot T_{\text{eff}}^\odot/M_\odot) \right]^{1/2} \end{aligned} \quad (12) $$

with $\mathcal{R}$ the constant of gases; $\mu$ the average molecular weight; gravitational constant $G$; $R_\odot$, $T_{\text{eff}}^\odot$, and $M_\odot$ are the solar radius, effective temperature, and mass, respectively. In the following, we use $\gamma = 1$ to otherwise describe deviations from a strict Keplerian rotation. Whether the temperature of CD follows the geometrical dilution of the stellar bolometric radiation or whether it is isothermal, we take either $\beta = (2 + \gamma - 1/2)/2 = 1.25$ or $\beta = (2 + \gamma)/2 = 1.5$, respectively.

The function $D(R)$ can be specified only if a theory is used to describe the formation and further evolution of the CD, as done in VDD models (Narita et al. 1994; Okazaki 2001; Haubois et al. 2012; Ghoreyshi et al. 2021; Marr et al. 2021). In this work, we assume that the R-dependent density profile of the disk is parameterized as $D(R) = (R_o/R)^n$, where the exponent $n$ is a function of R. The density exponent, $n$, that is reported currently in the literature, actually represents the global slope of the disk density variation with R. However, from Haubois et al. (2012), it is clear that depending on the evolution phase of the circumstellar disk, $n(R)$ is a strong function of R, mainly in proximity to the star, where emission component in lines such as Hγ and Hδ are formed. Anticipating the possibility that the regions where the emission of Hα and those of Hγ and Hδ do not respond to the same value of the density exponent $n$, as in Zorec et al. (2007), we adopted the following expression for $n(R)$:

$$ n(R) = \left[ \frac{(n_2 - n_1)}{\pi} \right] \arctan \left\{ Q \times \left[ \frac{R}{R_o} - \frac{R_D}{R_o} \right] \right\} + \frac{1}{2}(n_1 + n_2), \quad (13) $$

whose value changes from $n(R) = n_1$ within $R_o \lesssim R \lesssim R_D$ to $n(R) = n_2 > n_1$ for $R \gtrsim R_D$. In Eq. (13), $Q \lesssim 10$ is a constant chosen freely. However, no significant changes are carried on the determination of parameters if we adopt a simple step relation $n(R) = n_1$ at $R_o \le R \le R_D$ and $n(R) = n_2$ for $R \gtrsim R_D$.

From the theory of VDD (Lee et al. 1991; Okazaki 2001, 2007), we know that in isothermal disks at the steady state, it is $n = 3.5$ over the entire disk radial extent. The results obtained by Zorec et al. (2007) based on the study of Fe II emission lines in the blue spectral range of Be stars, suggest that $n_1 \lesssim 1.0$ for $R \lesssim R_D \sim 3R_o$, while studies on the SED in the near- and far-IR of Be stars (Waters 1986; Vieira et al. 2015) suggest $n_2 \approx 2.5-3.9$ for $R \gtrsim R_D$. The modeling of the SED in the 3500–10 500 Å wavelength interval of a set of Be stars by Moujtahid (1998),
Moujtahid et al. (2000a,b) and Chauville et al. (2001) suggested an average exponent $<n> = 1.6 \pm 0.2$. Granada et al. (2010) used $n_1 = 0.5$ in a study of emissions in the hydrogen Humphry's series, which form in CD regions near the central star. Studies of the temporal evolution of VDD density distributions under different formation scenarios (Haubois et al. 2012) end up concluding that $n$ varies from 0.5 to 4.5, but sometimes $n < 0$ in the base of the disk. The function $n(R)$ on the distance from the star, it is sensitive to the kinematic viscous coefficient and depends on the disk decretion history (Haubois et al. 2012; Vieira et al. 2017; Ghoreyshi et al. 2018, 2021; Marr et al. 2021).

In a large series of papers, disk models for Be stars were assumed to be isothermal (Marlborough 1969; Poeckert & Marlborough 1978a,b; Hummel 1994; Catanzaro 2013). More refined discussions on the distribution of the temperature in the CD of Be stars obtained high non-uniform temperature distributions, whose global behavior depends on the physical inputs and assumptions made in the models. Millar & Marlborough (1998, 1999), Carciofi & Bjorkman (2006, 2008), Sigut & Jones (2007), McGill et al. (2013) noted that not far from the star, the disk temperature decreases in the equatorial region and increases somewhat at larger distances and higher Z-coordinates. On the contrary, by including viscous heating, Kurfürst et al. (2018) have found that the temperature increases in the equatorial regions and decreases towards larger distances in both R and Z.

However, it has long been known that the source function of Balmer lines is strongly dominated by photoionization processes (Thomas 1957, 1965; Jefferies 1968; Mihalas 1978; Hubeny & Mihalas 2014). This means that the production and destruction of line photons is dominated by the stellar radiation field. The source function is thus dissociated from local temperature and density in the line formation region. To formalize a marginal dependencies of the source function with the CD temperature, we assume that the local temperature is determined by the geometrical dilution of the stellar bolometric flux $[T(R)/T_{\text{eff}}]^4 = (1/2)[1 - (1 - (R_o/R)^2)^{1/2}]$, which is close to the relation found by Moujtahid (1998) $T/T_{\text{eff}} = 0.19 + 0.61/R$ in a study of the observed visible energy distribution in 21 Be stars. Moreover, as the line emission efficiency depends on the electron density squared, the emission rapidly decreases with Z, reducing thus the extent of effective zone perpendicular to the equatorial plane for the emission of photons. We assume then a uniform temperature in the Z-direction.

Furthermore, we reduce the CD line formation region into an equivalent cylinder having an average half-height $\bar{H}$ and a radial extent $\Delta R$ that extends from a given radius, $R$. Here, $\Delta R$ is the distance over which it is recovered 99% of the disk radial opacity integrated from $R$ to $R + \Delta R$. Using Eq. (11), it follows that the opacity of the disk in the Z-direction averaged over the $\Delta R$ distance, is the same as for a cylinder with height $\bar{H} = 3\sqrt{2}h(R)$, where $h(R)$ is the weighted average of the disk density scale-height, $h(R)$, over the distance, $\Delta R$

$$ \overline{h(R)/R_o} \approx [1 + \beta(\Delta R/R)] h(R), \quad \text{for } \Delta R/R \lesssim 1, \quad (14) $$

where $h(R)$ is given by Eq. (12).

Although in our discussion, $\bar{H}$ is considered a free fitting parameter, the value obtained from the analytical expression is adequate enough in almost all cases to obtain the best fit of a given observed line emission. In Tables 7–9, we refer the model $\bar{H}$ to the value given by $3\sqrt{2}h(R)$.

<table><caption>**Table 7.** Parameters obtained from the fit of the observed Hα emission lines for their CD formation region.</caption><thead><tr><th rowspan="2">Star</th><th colspan="2">S°λ/F<sub>c</sub></th><th colspan="3">H&#773;/R<sub>o</sub></th><th rowspan="2">V<sub>rot</sub><br>km s<sup>-1</sup></th><th rowspan="2">V<sub>exp</sub><br>km s<sup>-1</sup></th></tr><tr><th>model</th><th>fit</th><th>τ<sub>λo</sub></th><th>R/R<sub>o</sub></th><th>ΔR/R<sub>o</sub></th><th>model</th><th>fit</th></tr></thead><tbody><tr><td>1</td><td>0.053</td><td></td><td>0.30</td><td>5.5</td><td>2.3</td><td>1.06</td><td>270</td><td>-35</td></tr><tr><td>2</td><td>0.064</td><td>0.100</td><td>0.80</td><td>4.5</td><td>2.0</td><td>0.98</td><td>1.00</td><td>0</td></tr><tr><td>3</td><td>0.052</td><td></td><td>1.10</td><td>7.5</td><td>3.6</td><td>1.78</td><td></td><td>0</td></tr><tr><td>4</td><td>0.063</td><td>0.100</td><td>1.05</td><td>4.0</td><td>1.6</td><td>0.59</td><td></td><td>0</td></tr><tr><td>5</td><td>0.072</td><td>0.135</td><td>1.40</td><td>2.2</td><td>1.7</td><td>0.42</td><td></td><td>-25</td></tr><tr><td>6</td><td>0.072</td><td></td><td>1.05</td><td>4.3</td><td>2.2</td><td>0.98</td><td></td><td>0</td></tr><tr><td>7</td><td>0.065</td><td>0.160</td><td>0.90</td><td>4.0</td><td>1.7</td><td>0.81</td><td>0.90</td><td>-5</td></tr><tr><td>8</td><td>0.094</td><td></td><td>1.70</td><td>4.0</td><td>6.1</td><td>1.23</td><td></td><td>-55</td></tr><tr><td>9</td><td>0.066</td><td>0.200</td><td>0.76</td><td>4.9</td><td>2.4</td><td>1.18</td><td></td><td>5</td></tr><tr><td>10</td><td>0.047</td><td>0.100</td><td>0.40</td><td>2.0</td><td>0.8</td><td>0.34</td><td></td><td>-10</td></tr><tr><td>11</td><td>0.067</td><td></td><td>1.60</td><td>8.1</td><td>5.5</td><td>1.82</td><td></td><td>30</td></tr><tr><td>12</td><td>0.054</td><td></td><td>2.00</td><td>3.3</td><td>1.0</td><td>0.59</td><td></td><td>-18</td></tr><tr><td>13</td><td>0.064</td><td></td><td>1.40</td><td>6.9</td><td>3.9</td><td>1.53</td><td></td><td>0</td></tr><tr><td>14</td><td>0.091</td><td></td><td>3.00</td><td>3.5</td><td>3.6</td><td>0.76</td><td></td><td>0</td></tr><tr><td>16</td><td>0.050</td><td>0.070</td><td>0.35</td><td>5.0</td><td>1.8</td><td>0.85</td><td></td><td>-5</td></tr><tr><td>17</td><td>0.056</td><td>0.120</td><td>0.43</td><td>5.0</td><td>2.1</td><td>0.85</td><td></td><td>-2</td></tr><tr><td>19</td><td>0.076</td><td></td><td>2.10</td><td>5.0</td><td>3.6</td><td>1.44</td><td></td><td>0</td></tr></tbody></table>

<table><caption>**Table 8.** Parameters obtained from the fit of the observed Hγ emission lines for their CD formation region.</caption><thead><tr><th rowspan="2">Star</th><th colspan="2">S°λ/F<sub>c</sub></th><th colspan="3">H&#773;/R<sub>o</sub></th><th rowspan="2">V<sub>rot</sub><br>km s<sup>-1</sup></th><th rowspan="2">V<sub>exp</sub><br>km s<sup>-1</sup></th></tr><tr><th>model</th><th>fit</th><th>τ<sub>λo</sub></th><th>R/R<sub>o</sub></th><th>ΔR/R<sub>o</sub></th><th>model</th><th>fit</th></tr></thead><tbody><tr><td>2</td><td>0.182</td><td></td><td>0.06</td><td>2.1</td><td>1.4</td><td>0.42</td><td>230</td><td>-10</td></tr><tr><td>8</td><td>0.227</td><td>0.500</td><td>0.09</td><td>2.2</td><td>2.0</td><td>0.47</td><td>240</td><td>-10</td></tr><tr><td>9</td><td>0.182</td><td>0.500</td><td>0.01</td><td>1.6</td><td>1.8</td><td>0.38</td><td>150</td><td>-17</td></tr><tr><td>11</td><td>0.186</td><td>0.400</td><td>0.50</td><td>2.1</td><td>1.6</td><td>0.34</td><td>230</td><td>10</td></tr><tr><td>13</td><td>0.188</td><td>0.210</td><td>0.07</td><td>3.1</td><td>1.0</td><td>0.51</td><td>200</td><td>-10</td></tr><tr><td>14</td><td>0.233</td><td></td><td>0.04</td><td>3.5</td><td>3.6</td><td>0.96</td><td>320</td><td>0</td></tr><tr><td>19</td><td>0.211</td><td></td><td>0.14</td><td>3.7</td><td>2.1</td><td>0.89</td><td>300</td><td>-20</td></tr></tbody></table>

<table><caption>**Table 9.** Parameters obtained from the fit of the observed Hδ emission lines for their CD formation region.</caption><thead><tr><th rowspan="2">Star</th><th colspan="2">S°λ/F<sub>c</sub></th><th colspan="3">H&#773;/R<sub>o</sub></th><th rowspan="2">V<sub>rot</sub><br>km s<sup>-1</sup></th><th rowspan="2">V<sub>exp</sub><br>km s<sup>-1</sup></th></tr><tr><th>model</th><th>fit</th><th>τ<sub>λo</sub></th><th>R/R<sub>o</sub></th><th>ΔR/R<sub>o</sub></th><th>model</th><th>fit</th></tr></thead><tbody><tr><td>8</td><td>0.263</td><td></td><td>0.04</td><td>2.0</td><td>2.1</td><td>0.42</td><td>250</td><td>-130</td></tr><tr><td>9</td><td>0.230</td><td>0.190</td><td>1.20</td><td>2.4</td><td>1.3</td><td>0.51</td><td>210</td><td>-40</td></tr><tr><td>11</td><td>0.224</td><td>0.500</td><td>0.20</td><td>1.9</td><td>1.7</td><td>0.30</td><td>230</td><td>8</td></tr><tr><td>13</td><td>0.230</td><td></td><td>0.03</td><td>2.8</td><td>1.0</td><td>0.47</td><td>200</td><td>-30</td></tr><tr><td>14</td><td>0.275</td><td></td><td>0.04</td><td>3.2</td><td>2.5</td><td>0.64</td><td>320</td><td>-20</td></tr><tr><td>19</td><td>0.253</td><td></td><td>0.07</td><td>3.2</td><td>1.5</td><td>0.72</td><td>320</td><td>-20</td></tr></tbody></table>

## 7.2. The source function

In the frame of an equivalent two-level atom with continuum, the main characteristic of the source functions for the first Balmer lines is that they are strongly dominated by photionizations and radiative recombinations, where collisional ionization and excitation rates can be neglected (Thomas 1957, 1965; Jefferies 1968; Mihalas 1978; Hubeny & Mihalas 2014). The non-LTE source function of these lines are then determined only by the radiation field of the central star. The source function of the first Balmer lines with continuum in an equivalent emitting ring with uniform
<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Regions in the disk that contribute in different ways to the observed line emission, according to our simplified representation of the star-disk system.</figcaption></figure>

physical properties can then be approached with the following expressions (Mihalas 1978; Cidale & Ringuelet 1989; Hubeny & Mihalas 2014)

$$S_{\lambda}(\tau_o) = \begin{cases} [\eta/(1+\eta)]^{1/2} B_{\lambda}^* & \text{for } \tau_o < 1 \\ [\eta/(1+\eta)]^{1/2} B_{\lambda}^* \tau_o^{1/2} & \text{for } \tau_o \ge 1, \end{cases} \quad (15)$$

where we do not consider the thermalization of the source function that should be produced for $\tau_o \ge 10-10^2$. In Eq. (15), $\tau_o$ is the optical depth in the central wavelength of the Balmer $\lambda$ line; $\eta$ is the “sink” term in the line source function, $S_\lambda$; $B^*$ is the “source” factor of $S_\lambda$. This approximation was already used in Vinicius et al. (2006) and Arias et al. (2007). The expressions for the radiation sink and the source factors $\eta$ and $B_*$ are given in (Thomas 1965; Hubeny & Mihalas 2014).

In Tables 7–9, the source function ratios $S^\circ_\lambda/F_*= [\eta/(1+\eta)]^{1/2}B^*/F_*$ are given for the H$\alpha$, H$\gamma$, and H$\delta$ lines of the studied stars. Within the same framework of radiation-dominated transitions the non-LTE departure coefficients were calculated.

### 7.3. Formal solution of the radiation transfer for the CD

From the simplifications adopted on the geometrical shape of the CD, as a radiation emitting region characterize by a uniform opacity $\tau_Z$ in the direction perpendicular to the equator, can be considered to behave as an equivalent thin ring, having a total height of $2\bar{H}$ and radius $R$.

Although the simulation of the formation of emission lines from a flat disk of uniform height does not offer particular difficulties (see Horne & Marsh 1986), calculations as a function of the inclination angle, $i$, in the frame of an equivalent ring are a little more delicate. The shape of the emission ring-region projected towards the observer is sketched in Fig. 8. The radiation field coming from the ring-star system can be represented using the formal solution of the equation of radiation transfer at each point $(x, y)$ projected on the background plane of the sky. For optical depths $0 \le \tau_o < 1$, according to Eq. (15), the source function $S_\lambda = S^\circ_\lambda$ of the studied Balmer lines is constant. The formal solution of the radiation transfer equation applied to the star-disk system applied to the geometrical configuration given in Fig. 8

is as follows:

$$ \left. \begin{aligned} I_{\lambda}^{u}(x, y, V - V_{r}) &= I_{\lambda}^{*}(x, y) \exp[-\tau_{\lambda}^{f}(x, y, V - V_{r})/\mu(x, z)] \\ & \quad + S^{\circ}_{\lambda}\{1 - \exp[-\tau_{\lambda}^{f}(x, y, V - V_{r})/\mu(x, y)]\} \\ & \quad (\text{stellar emission absorbed by the front shell}) + (\text{emission from the front shell}) \\ I_{\lambda}^{b}(x, y, V - V_{r}) &= S^{\circ}_{\lambda}\{1 - \exp[-\tau_{\lambda}^{r}(x, y, V - V_{r})/\mu(x, y)]\} \\ & \quad \times \exp[-\tau_{\lambda}^{f}(x, y)/\mu(x, y)] \\ & \quad + S^{\circ}_{\lambda}\{1 - \exp[-\tau_{\lambda}^{f}(x, y, V - V_{r})/\mu(x, y)]\} \\ & \quad = (\text{emission from rear shell absorbed by front shell}) + (\text{emission by front shell}) \end{aligned} \right\} \quad (16) $$

where $I_\lambda^*$ is the contribution to the field of radiation by the central star. In these expressions, $\mu(x, y)$ is the cosine between the normal to the ring at a test point $(x, y)$ and the sight direction, $V$ is the wavelength displacement in the spectral line in velocity units, and $V_r$ is the radial velocity at the point $(x, y)$, defined below in Sect. 7.4 (Eq. (21)).

For optical depths $\tau_o \ge 1$, the formal solution for the radiation field introduces such integrals as:

$$ \begin{aligned} I_{\lambda} &= S^{\circ}_{\lambda} \int_{1}^{\tau_o} \tau_o^{1/2} \exp(-\tau_{\lambda}) d\tau_{\lambda}, \\ \tau_{\lambda} &= \tau_o \Phi_{\lambda}, \\ I_{\lambda} &= \Phi_{\lambda} [\tau_o^{1/2} f(\Phi_{\lambda} \tau_o) - f(\tau_o)], \end{aligned} \quad (17) $$

where $\Phi_\lambda$ is the intrinsic line absorption profile and $f(t)$ is determined by the incomplete gamma function, $\gamma(t, 3/2)$.

These relations can be used to compose the field of radiation due to overlapping absorbing and emitting layers of the front and rear parts of the equivalent ring. In the optical depths, a distinction is made between $\tau_\lambda^r(x, y, V)$ and $\tau_\lambda^f(x, y, V)$ of the “rear” and “front” layers, respectively, to caution that the radiation field is affected by the Doppler effect. This is due to the possibility of several apparent radial velocity compositions produced by the rotation and expansion velocities of the projected shells relative to the observer.

Finally, the radiation flux in the lines emitted by the star-disk system is given by:

$$F(x, y) = \int_{S} I_{\lambda}(x, y)dxdy, \quad (18)$$

where $S$ is the surface of the star-disk system projected onto the sky background, which in addition to the ring surface shown in Fig. 8, encompasses the projected emitting surface of effective cylinder from $R$ to $R + \Delta R$.

The expressions in Eq. (16) can be adapted to several configurations that depend on the value of the height $\bar{H}: \bar{H} \le R_o$ or $\bar{H} > R_o$ and on the value of the inclination angle $i_{RH} = \arctan(\Delta R/2\bar{H})$. We previously used this procedure in Floquet et al. (2000), Arias (2004), Vinicius et al. (2006) and Arias et al. (2007).

### 7.4. Line opacity

We considered the opacity of the studied Balmer lines parameterized according to their value in the $Z$-direction, which at the respective central wavelength, $\lambda_o$, and in the non-LTE approximation is:

$$ \begin{aligned} \tau_{\lambda_o} &= \alpha_{\lambda_o} \overline{N_2(R)_Z} \{1 - (b_n/b_2) \exp[-(hc/\lambda_o kT)]\}, \\ \alpha_{\lambda_o} &= (\sqrt{\pi} e^2)/(m_e c)(\lambda_o f)/V_D, \end{aligned} \quad (19) $$
where $\overline{N_2(\bar{R})}$ represents the population of the atomic level $n = 2$ in the Z-direction, averaged over the radial distance from R to R + $\Delta R$, f is the oscillator strength of the given transition, $V_D$ is the total Doppler velocity of its intrinsic absorption line profile that encompasses not only thermal but other possible macroscopic velocity components of statistical nature, and $b_n$ represents the non-LTE departure coefficients for atomic levels n. The term $(b_n/b_2) \exp[-(hc/\lambda_o kT)] \ll 1$ can be neglected in most cases.

Writing the line opacity as:

$$ \tau_Z = \tau_{\lambda_o} \Phi(\Delta\lambda/\Delta\lambda_D), \quad (20) $$

the intrinsic Balmer line profile $\Phi(\lambda)$ formally corresponds to the composite expression for a Stark profile given in Gray (2008) or Hubeny & Mihalas (2014). We calculated it using the tables of Stark profiles given by Stehlé & Hutcheon (1999). However, in most cases, the gas density in the disk is low, so that the Stark line broadening is small.

The opacity in the Z-direction, $\tau_Z$, is adopted in this work as a free parameter. For aspect angles, $i \le i_{RH} = \arctan(\Delta R/2\bar{H})$, the optical depth of the shell as a function of the inclination is then parameterized as $\tau(i) = \tau_Z/\cos i$ and then as $\tau(i) = \tau_Z \tan i_{RH}/\sin i$ when $i \ge i_{RH}$. Apart from these i-angle dependencies of the calculated line emission characteristics, stronger inclination angle dependencies are carried by the apparent shape of the entire emission ring-region projected towards the observer, as sketched out in Fig. 8.

The emission intensity at an observed wavelength $\lambda$ corresponds to the sum of the emissions at wavelength offsets $\delta\lambda$ with respect to the central absorption of the atom, which, in turn, is displaced due to the macroscopic movement of the medium with respect to the observer. The absorption (or emission) intensity with which the atom contributes at $\lambda$ is finally determined by the profile $\Phi(\delta\lambda) = \Phi(\lambda - \lambda_o - \Delta\lambda)$, where $\Delta\lambda = \lambda_o \frac{V_r}{c}$ and $V_r$ is the total macroscopic radial velocity of the gas at a given point $(R, \theta)$ in the ring:

$$ V_r = [V_{\text{rot}} \sin \theta + V_{\text{exp}} \cos \theta] \sin i, \quad (21) $$

where $V_{\text{rot}}$ is the rotation velocity of the ring and $V_{\text{exp}} = V_{\text{exp}}(R_o/R)^\alpha$ is its “expansion” velocity that is used only to account for the line profile asymmetries. The quantity $V_{\text{rot}}$ is also only to fit at best the width of emission lines, which encompass the actual line broadening by the disk rotation integrated over its formation region and a host of other line broadening mechanisms. Hence, we do not expect to have well defined correlations between $V_{\text{rot}}$ and $V_{\text{rot}}^*/\sqrt{R}$.

For the line profile calculations we use velocities instead of wavelengths, so that the argument in the profile $\Phi(\lambda - \lambda_o - \Delta\lambda)$ is readily transformed into velocity displacements $[(\lambda - \lambda_o - \Delta\lambda)]/\Delta\lambda_D = (V - V_r)/V_D$, where $V - V_r$ measures in velocity units the offset $\lambda - \lambda_o$ in the observed profile.

The use of equivalent emitting (absorbing) shell with radius R, also requires a consideration of the Doppler displacements due to the changes of the macroscopic velocity fields in the Z-direction. These contributions are called “shear” components. Following Horne & Marsh (1986), it can be shown that the integrated effect on the line profile of the shear velocity components carries an enlarged line Doppler width $\Delta\lambda_D \rightarrow \Delta_D$,

$$ \Delta_D = \Delta\lambda_D \left[ 1 + \left( \frac{\lambda_o}{\Delta\lambda_D} \frac{V_{sh}}{c} \right)^2 \right]^{\frac{1}{2}}, \quad (22) $$

where the factor in braces acts as a supplementary broadening agent of the line profile. In Eq. (22), $V_{\text{sh}}$ is the contribution of the shear velocity component given by

$$ V_{\text{sh}} = \left[ \frac{1}{2} V_{\text{rot}} \sin \theta + \alpha V_{\text{exp}} \cos \theta \right] \frac{\bar{H}}{R} \cos \theta \tan i \sin i. \quad (23) $$

Shear effects are anisotropic with strong dependence on the inclination angle, i, and azimuth, $\theta$. These effects has almost no incidence on the optically thin lines, but it saturates somewhat the emission peaks as soon as $\tau_{\lambda_o} \ge 1$. In what follows, we assume $\alpha = 0$. Disks where $\alpha \ne 0$, namely, those with gradients of expansion velocities, can be studied using separate equivalent rings which may then reproduce multi-peak emission line profiles (Arias et al. 2007).

### 7.5. Fit of the line emission components and determination of disk parameters

In our approach the main free parameters are: (1) $\tau_Z$, the opacity of the ring in the Z-direction; (2) the radius $R/R_o$ of the ring; and (3) the rotation and expansion velocity parameters $V_{\text{rot}}$ and $V_{\text{exp}}$, respectively. At each iteration step, the extent $\Delta R$ of the emitting region in the disk is calculated as a function of the chosen value of $R/R_o$. We initiate the iterative fitting procedure of the line emission profiles by adopting the calculated value of $[\eta/(1+\eta)]^{1/2}B^*/F_c$ and of $\bar{H}/R_o$ that corresponds to the chosen value of R. In the search for finer adjustments, these quantities can be considered as free parameters. This enables us to test the relevance of the assumptions made on the models of CDs relative to $\eta^{1/2}B^*/F_c$ and $\bar{H}/R_o$. The inclination angle of the star-disk system is taken from Table D.2 already determined with the remaining astrophysical parameters of stars.

All fitting parameters were obtained through a trial-and-error procedure, which is not arbitrary because the global intensity of the emission, the central absorption component in emission lines, and the width of the line wings are sensitive to the ratios between the searched quantities. Nevertheless, we cannot discuss their uniqueness – although it would be possible perhaps if a robust method of adjustment such as the “annealing method” were used (Metropolis et al. 1953). The parameters displayed in Tables 7–9 must then be considered first-order estimates. We note that there are additional uncertainties affecting them which are related to the shape of the emission line profiles proper, mainly when it comes to the extracted ones in Hγ and Hδ lines.

### 7.6. Extent of the line emission regions

The Gaussian fit of the $H\alpha$ emission intensity maps obtained with long baseline interferometry for about 20 Be stars (Quirrenbach et al. 1997; Tycner et al. 2004, 2005, 2006, 2008; Koubský et al. 2010; Delaa et al. 2011; Touhami et al. 2013; de Almeida et al. 2020), produce an average observational estimate of the extent of the emission region $\langle E/R_o \rangle = 7.4 \pm 3.5$. This value must be compared with our $R + \Delta R$ determinations, which for 17 program Be stars lead to the average $\langle (R + \Delta R)/R_o \rangle = 7.4 \pm 2.8$ that can be considered in good agreement with interferometrical estimates.

### 7.7. Density distribution in the CD

Only a theory for the CD formation and evolution can determine the function $D(R)$ in Eq. (11) as in VDD model that depends on the chosen viscosity. In this work, we simply assume an analytic
<table><caption>**Table 10.** Mean non-LTE departure coefficients in the line-forming regions, the coefficient $n_1$ for the inner disk density distribution (Eq. (13)), the base density, $\rho_o$, of the CD and the CD mass in the Balmer line formation region, $M_{CD}$.</caption><thead><tr><th>Star</th><th>b<sub>2</sub><br>(Hα)</th><th>b<sub>2</sub><br>(Hγ)</th><th>b<sub>2</sub><br>(Hδ)</th><th>n<sub>1</sub></th><th>ρ<sub>o</sub> ± σ<sub>ρ<sub>o</sub></sub><br>(g cm<sup>-3</sup>)</th><th>M<sub>CD</sub><br>(M<sub>☉</sub>)</th></tr></thead><tbody><tr><td>1</td><td>0.04</td><td>—</td><td>—</td><td>—</td><td>1.5 ± 0.2 × 10<sup>-11</sup></td><td>2.0 × 10<sup>-08</sup></td></tr><tr><td>2</td><td>0.25</td><td>1.25</td><td>—</td><td>4.8</td><td>6.4 ± 1.3 × 10<sup>-11</sup></td><td>2.1 × 10<sup>-08</sup></td></tr><tr><td>3</td><td>0.01</td><td>—</td><td>—</td><td>—</td><td>1.0 ± 0.2 × 10<sup>-10</sup></td><td>4.0 × 10<sup>-07</sup></td></tr><tr><td>4</td><td>0.36</td><td>—</td><td>—</td><td>—</td><td>6.0 ± 0.3 × 10<sup>-11</sup></td><td>2.4 × 10<sup>-08</sup></td></tr><tr><td>5</td><td>1.78</td><td>—</td><td>—</td><td>—</td><td>2.0 ± 0.1 × 10<sup>-13</sup></td><td>1.0 × 10<sup>-09</sup></td></tr><tr><td>6</td><td>0.55</td><td>—</td><td>—</td><td>—</td><td>5.7 ± 0.7 × 10<sup>-12</sup></td><td>7.7 × 10<sup>-08</sup></td></tr><tr><td>7</td><td>0.42</td><td>—</td><td>—</td><td>—</td><td>2.7 ± 1.3 × 10<sup>-12</sup></td><td>5.2 × 10<sup>-09</sup></td></tr><tr><td>8</td><td>2.05</td><td>4.26</td><td>4.31</td><td>4.8</td><td>1.5 ± 0.4 × 10<sup>-11</sup></td><td>6.5 × 10<sup>-08</sup></td></tr><tr><td>9</td><td>0.23</td><td>1.68</td><td>1.18</td><td>1.1</td><td>3.7 ± 3.7 × 10<sup>-12</sup></td><td>5.0 × 10<sup>-08</sup></td></tr><tr><td>10</td><td>0.62</td><td>—</td><td>—</td><td>—</td><td>6.0 ± 0.3 × 10<sup>-12</sup></td><td>2.0 × 10<sup>-10</sup></td></tr><tr><td>11</td><td>0.03</td><td>1.38</td><td>1.54</td><td>7.1</td><td>1.7 ± 0.5 × 10<sup>-10</sup></td><td>9.5 × 10<sup>-8</sup></td></tr><tr><td>12</td><td>0.29</td><td>—</td><td>—</td><td>—</td><td>4.0 ± 0.3 × 10<sup>-13</sup></td><td>1.6 × 10<sup>-09</sup></td></tr><tr><td>13</td><td>0.05</td><td>0.73</td><td>0.92</td><td>6.1</td><td>8.7 ± 1.0 × 10<sup>-11</sup></td><td>9.1 × 10<sup>-08</sup></td></tr><tr><td>14</td><td>2.53</td><td>2.53</td><td>3.11</td><td>0.5</td><td>1.2 ± 0.5 × 10<sup>-12</sup></td><td>4.3 × 10<sup>-09</sup></td></tr><tr><td>16</td><td>0.04</td><td>—</td><td>—</td><td>—</td><td>1.1 ± 0.2 × 10<sup>-11</sup></td><td>1.5 × 10<sup>-08</sup></td></tr><tr><td>17</td><td>0.07</td><td>—</td><td>—</td><td>—</td><td>1.4 ± 0.2 × 10<sup>-11</sup></td><td>1.0 × 10<sup>-08</sup></td></tr><tr><td>19</td><td>0.53</td><td>1.12</td><td>1.55</td><td>5.9</td><td>1.9 ± 0.4 × 10<sup>-11</sup></td><td>1.5 × 10<sup>-07</sup></td></tr></tbody></table>

<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Correlation of Huang's radii derived with $V_{rot}$ against line formation extends $R + \Delta R$. In the ordinates the values are from Table 5.</figcaption></figure>

expression for $D(R)$ to be valid for the studied line formation region.

In recent interpretations of the photometric, spectroscopic, and polarimetric variations of $\omega$ CMa and 66 Oph (Ghoreyshi et al. 2021; Marr et al. 2021), $D(R)$ is characterized in Eq. (11) by a single global exponents, $n$, that describe the several disk formation phases. These values range from some $n \lesssim 4.0$ at CD formation phases to $n \gtrsim 2.0$ at disk dissipation epochs. Nevertheless, much higher values of $n$ may describe the disk formation controlled by periodic mass injection rates, where the injection and dissipation timescales are similar, as shown by Haubois et al. (2012).

For the $H\alpha$ line formation region, we made the conservative assumption $n_2 = 3.0$ for all program stars, because in our method we cannot determine it otherwise. Since the formation region of the $H\alpha$ emission component is situated at radii of $R \gtrsim 4R_o$, the value of $n_1$ has no incidence on the determination of the hydrogen disk base density $\rho_o$ g cm$^{-3}$.

In stars where we could extract reliable emission profiles for the $H\gamma$ and $H\delta$ lines, we estimated the exponent $n_1$ that characterizes the inner layers of the CD in a relative way by demanding that the values of $\rho_o$ derived from the optical depths of all three lines be the same and equal to that found with $H\alpha$. In Table 10, we give the values of $n_1$ that were thus derived. These values are, on average, $\langle n_1 \rangle = 4.3 \pm 2.5$ – similar to the estimates by Vieira et al. (2017), where $n$ is nevertheless meant to represent the entire radial density distribution at CD forming phases. We can also note that the average $[n_1 + n_2(= 3)]/2 = 3.7 \pm 1.3$, which is close to the values of $n$ estimated by several authors which used a single $n$ to characterize the entire CD, namely, Silaj et al. (2010), Arcos et al. (2017).

The large variety of $n_1$ values estimated in this work may then correspond to different evolution phases of disks, from their formation to their likely steady state, or dissipation. In cases where spectroscopic observations of Be stars exist only for

specific dates, a method such as the one developed in the present work, which is of relative simple use, may still bring useful insights on the nature of CDs. In the present case, the higher values of $n_1$ suggest that in many Be stars the disk formation could respond to periodic mass injection rates, as discussed in Haubois et al. (2012). It is worth noting that further detailed studies of the $H\gamma$ and $H\delta$ emission lines at different mass ejection episodes could reveal the physical properties of regions where the viscous transfer of angular momentum is probably organized to the rest of the CD regions and responsible for its Keplerian rotation.

From the radial optical depth $\tau_Z$, the non-LTE departure coefficients averaged over the respective line formation regions, and assuming that the He abundance in number is roughly $N_{He} \approx N_{H}/4$, we readily determine the CD base density, $\rho_o$, which is given in Table 10. These values are of the same order of magnitude than adopted in other more detailed studies (Silaj et al. 2010; Catanzaro 2013; Arcos et al. 2017; Vieira et al. 2017).

The integration of Eq. (11) together with the use of the obtained values of $\rho_o$, lead us to an estimate of the total mass, $M_{CD}$, gathered in the formation regions of the studied Balmer line emissions that are listed in Table 10.

## 7.8. Additional comments on the CDs

Our CD emission model is challenged with the star $n = 5$. This object is a binary system, where it is likely that contributions by two distinct regions in the disk would be needed to account for the entire $H\alpha$ line emission profile in a more reliable way.

Arcos et al. (2017) obtained a correlation at 61% between Huang's radii and the radii $R_{90}$, which represents the distance in the CD up to which 90% of the line emission is formed. We note that the $V_{sin i}$ parameters used to this end were not measured, but obtained as one among several parameters used to fit the $H\alpha$ emission line profiles. We obtain a similar correlation at 68% between Huang's radii calculated with the fitting quantity $V_{rot}$ and then distance $R + \Delta R$ (see Fig. 9). However, because in both approaches Huang's radii are calculated with fitted rotational velocities, doubts can be cast on the physical reliability of these correlations.

Then, $V_{exp}$ is a fitting parameter needed to account for the line emission asymmetry. It resumes perhaps the effects carried by the known quasi-cyclic V/R variations attributed to the global one-armed disk density oscillations (Kato 1983; Okazaki 1991, 1997).

Finally, in this work we did not attempt to interpret the correlations implying emission lines and photometric measurements, simply because it would require to include a study of the formation of the continuum energy distribution with approaches like
those carried by Moujtahid (1998), Moujtahid et al. (2000a,b),
which are not in the scope of the present paper.

8. Discussions

8.1. Comments on the pncr parameters

A rotating star of mass, *M*, behaves as an object with lower effective mass *M* − Δ*M*, where Δ*M* is proportional to the pressure-weighted average of the centrifugal to the gravitational force over the entire object (Sackmann 1970; Clement 1979). Consequently, the rotating star produces a lower bolometric luminosity during roughly the first half of the main sequence (MS) evolutionary phase and has longer evolving time scales than its non-rotating counterpart. Rotational mixing carries fresh H-fuel into the convective core, which still prolongs its life in the MS phase (Maeder & Meynet 2000). This mixing also transports helium (and other H-burning products) into the radiative envelope, which lowers its opacity and contributes to the enhancement of the emitted stellar luminosity over that produced by the non-rotating object with same mass, taking place near to the second half of the MS phase (Maeder 2009).

The apparent luminosity depends mainly on the geometri-
cal deformation of the star, the gravitational darkening effect
(GD) and the inclination angle at which is seen the object.
Generally, the GD effect produces an average lower apparent
effective temperature and the geometrical deformation of the
star leads to an apparent average lower effective gravity. This
explains that the pncr effective temperatures and gravities are
larger than the apparent ones. Concerning the bolometric lumi-
nosity, the picture is less clear and strongly depends on the
combination of the astrophysical parameters, mainly on the rota-
tion rate and on the stellar aspect angle. Roughly, the ratio of the
apparent to the pncr bolometric luminosity is L<sup>app</sup>(i, Ω)/L<sup>pncr</sup> ≈
[S(i, Ω)/S<sub>pncr</sub>][⟨T<sup>app</sup><sub>eff</sub>(i, Ω)⟩/T<sup>pncr</sup>]<sup>4</sup>, where “S” is the area of the
apparent stellar hemisphere and that of the non-rotating counter-
part, respectively. In this relation we take the effective tempera-
ture of the rotating object averaged over the deformed observed
stellar hemisphere. See details on the behavior of L<sup>app</sup>(i, Ω) and
T<sup>app</sup><sub>eff</sub>(i, Ω) in Zorec et al. (2016). The conditions L<sup>pncr</sup> ≥ L<sup>app</sup>(i, Ω)
or L<sup>pncr</sup> ≤ L<sup>app</sup>(i, Ω) are controlled by the competition between
S(i, Ω) and T<sup>app</sup><sub>eff</sub>(i, Ω). The apparent Vsin i parameter is under-
estimated because the equatorial regions, which have the larger
surface linear rotational velocities, reduce their contribution to
the rotational line broadening due to the GD effect. A compari-
son between the pncr parametric set (Vsin i, T<sub>eff</sub>, log g, log L/L<sub>☉</sub>)
obtained for Ω/Ω<sub>c</sub> = 0.95, with the apparent one is shown in
Fig. 10.

A test on the reliability of the pncr parameters can be done
by comparing the pncr stellar masses with the masses obtained
using detached binary systems. The pncr stellar masses are
meant to represent the actual ones. This comparison is shown
in Fig. 11. The blue points are for masses determined using
the apparent astrophysical parameters and the red ones are for
masses determined with the pncr parameters. In this figure,
the average (log M/M<sub>☉</sub>, log T<sub>eff</sub>) relations are from Harmanec
(1988, green curves) and from Xiong et al. (2023, brown
curves), with the latter obtained in the frame of the LAMOST
medium-resolution spectroscopic survey. We see that the pncr
masses closely lie, within the uncertainty strips, around the
mean relation of masses measured in detached binary systems.

In principle, an O-, B- or A-type star may acquire the
required properties to display the Be phenomenon in several

<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Comparison between the pncr parametric set (Vsin i, T<sub>eff</sub>, log g, log L/L<sub>☉</sub>) obtained for Ω/Ω<sub>c</sub> = 0.95 with the apparent one. The average uncertainties of the plotted quantities are indicated. The diagonal in red represents the identity relation.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 11.** Comparison of Be stellar masses obtained in this paper with those of non-Be stars against the effective temperature. Blue dots represent the masses obtained from the apparent astrophysical parameters. Red dots represent the pncr stellar masses determined for Ω/Ω_c = 0.95. The green curve and the dotted green lines represent the average empirical relation curve with the error margins obtained by Harmanec (1988, __________). The red curve and the dotted brown lines represent the average empirical relation with the error margins obtained from LAMOST data Xiong et al. (2023, __________). The typical error bar that characterizes the apparent and pncr mass determinations is also shown. The points with a red circle identify the star N° 5.</figcaption></figure>
ways. It may be due to: formation and pre-main sequence evolution circumstances of single stars; internal angular momentum evolution properties during the MS phase (see references in Zorec 2023); a rejuvenating process as a consequence that one of the components has gained mass and angular momentum during mass transfer phases in binary systems (Packet 1981; Pols et al. 1991; de Mink et al. 2013). Depending on the several circumstances that may characterize the mass transfer phenomenon, the resulting Be stars should be binaries or single stars (Shao & Li 2014; Hastings et al. 2021). If they are binaries, their companions can be stripped sub-dwarf O/B (sdO/B) He stars (Be<sub>He</sub>), white dwarfs (Be<sub>WD</sub>, neutron stars (Be<sub>NS</sub>), or black holes (Be<sub>BH</sub>). Single Be stars can be merger objects (Be<sub>merger</sub>), disrupted Be<sub>NS</sub> or Be<sub>BH</sub> binaries. A significant number of recent observational data show that Be stars are binaries indeed, where the secondary is probably a stripped sdO/B object (Bodensteiner et al. 2020; El-Badry et al. 2022, and references therein). The paradigm of mass transfer may perhaps explain the nature of the star N° 5 studied in this paper and that we briefly discuss in Sect. 8.2.

From the above-cited set of possible companions of Be stars in binary system, those whose presence was neglected in the present work but may in principle have some effect on the $T_{\text{eff}}$ determination include: white dwarfs and the stripped sdO/B companions. However, as white dwarfs have radii of $R_{\text{WD}}/R_{\odot} = 0.01(M_{\odot}/M_{\text{WD}})^{1/3}$, masses of $0.15 \le M_{\text{WD}}/M_{\text{WD}} \le 1.4$, and effective temperatures of $8000 \le T_{\text{eff}}^{\text{WD}} \le 40000$ K, they can hardly carry any effect on the estimated $T_{\text{eff}}$ of our program Be stars. This may not be the case with hot sub-dwarf companions. According to Heber (2009), hot sub-dwarfs have most frequently effective temperatures in the range of $20000 \le T_{\text{eff}} \le 40000$ K, and radii in the interval $0.15 \le R_{\text{sd}}/R_{\odot} \le 0.25$. If we neglect their presence, the extreme errors that we will risk to commit on the estimate of the apparent effective temperatures of Be stars, can be estimated using the following extreme approximate bolo-metric luminosity relations, namely, both stars are seen or the sdO/B eclipses the Be star: $(s_B + s_{\text{sd}})T_{\text{Be}}^4 = s_B T_B^4 + s_{\text{sd}} T_{\text{sd}}^4$ and $s_B T_{\text{Be}}^4 = (s_B - s_{\text{sd}})T_B^4 + s_{\text{sd}} T_{\text{sd}}^4$, where $T_{\text{Be}}$, $T_B$ and $T_{\text{sd}}$ are: the effective temperature obtained with GIRFIT of the studied Be star, the actual effective temperature expected for the Be star, and the effective temperature of the sub-dwarf, respectively; the $s_B$ and $s_{\text{sd}}$ represent the areas of the Be and sub-dwarf respective apparent stellar surfaces. Because $r = s_{\text{sd}}/s_B \ll 1$, we have for both cases $T_B \approx T_{\text{Be}} \times \{[1 - r(T_{\text{sd}}/T_{\text{Be}})^4]\}$, which applied individually to our program objects enables us to estimate the effective temperature and the radius of the sub-dwarf star that carry on the $T_{\text{eff}}$ estimate an error on the order of the uncertainty characterizing our effective temperature determinations. The average effective temperature and radius of sdO/B stars thus calculated are: $\langle T_{\text{sd}} \rangle = 43600 \pm 13000$ K and $\langle R_{\text{sd}}/R_{\odot} \rangle = 0.27 \pm 0.06$, respectively. According to these results, the sdO/B stars in our sample should belong to the largest and hottest population of sdO/B stars. Such objects would certainly leave spectral signatures in the wavelength region studied in this work by He II λλ 4100, 4200, 4338, C IV λ 4442 and N III λ 4379 lines (Drilling et al. 2013), which are not obvious to identify in our spectra with certainty. We may then conclude that if these sdO/B are proven to exist, their effect on the $T_{\text{eff}}$ determinations will be at most of the order of uncertainties inherent to the spectral resolutions and S/Ns used in this work. However, star N° 5 decidedly requires a detailed study of the actual nature of the secondary, which (according to the characteristics of eclipses) could be larger than those considered in the above estimates.

It is worth recalling that the Be phenomenon proper remains linked to the physical characteristics of the O-, B- or A-type star, whatever the origin of Be stars as either single or merger byproducts. However, when dealing with binary systems, precise determinations of the astrophysical parameters of Be stars still requires that from now on it would be wise to carry out a thorough inquiry on the presence of companions and on their characteristics, as well as on their perturbation power on the energy distribution of the primary in the studied spectral interval.

## 8.2. The binary star N° 5

The astrophysical parameters of the star N° 5 were determined assuming that the object was a single star, although the light curve clearly indicates that it is a close eclipsing binary system (Klagyivik et al. 2017). Neglecting the presence of the companion and the perturbations that it may introduce to the SED of the primary would lead us to infer the improbable inclination angle for the primary $i \sim 50^\circ$ if $\Omega/\Omega_c = 0.95$. In fact, we could expect that the most probable inclination should be close to $i \sim 90^\circ$, because we are dealing with an eclipsing binary. However, assuming $i \sim 90^\circ$ and consequently fixing the value of the rotation parameter of the primary at $V \sin i = 337$ km s$^{-1}$, the GIRFIT fitting procedure produces bad adjustments of the synthetic spectra with the observed ones, and ends up with the following estimates of the apparent astrophysical parameters: $T_{\text{eff}} = 14000 \pm 700$ K and $\log g = 3.20 \pm 0.16$.

Moreover, the Hα emission line has likely a shell aspect. Other shell signatures are also found in the unique spectrum at our disposal (He I 6678 and 7065 Å). Unfortunately, the scarce data we have on this star do not help in determining whether the shell features are produced by a circumstellar disk or whether they are signatures attributed to the secondary component.

Due to the combined effect of gravitational and centrifugal forces, the circumstellar disk of the Be component will tend to be on or around its equatorial plane and probably near or over the binary orbital plan. However, we cannot give in to the temptation of speculating on a possible complex circumbinary envelope deploying strongly outside the binary plane, produced by a mass transfer episode between the components, such as those recently predicted by Smallwood et al. (2023). Such complex circumbinary gaseous structures can lead to misinterpretations of the position (inclination) of the circumstellar disk of the Be star.

## 8.3. Considering the rotational velocity rates

Regarding the corrections of effects carried by the rapid rotation on the apparent astrophysical parameters, we note that if we had enough independent apparent parameters as in the case of the calculations presented in Frémat et al. (2005) or Zorec et al. (2016), we could produce an estimate of the actual rotation frequency for each star individually. Since in this paper we do not have the required number of independent entry parameters, we ought to parameterize the ratios $\Omega/\Omega_c$.

Thus, we considered only the two extreme cases $\Omega/\Omega_c = 0.0$ and 1.0, for which the average velocity ratios corrected for GD effects are $\langle V/V_c \rangle = 0.85 \pm 0.20$ and $\langle V \sin i/V_c \rangle = 0.91 \pm 0.20$, respectively. In turn, these ratios correspond to $\langle \Omega/\Omega_c \rangle = 0.96 \pm 0.10$ and $\langle \Omega/\Omega_c \rangle = 0.99 \pm 0.05$, or to equatorial force ratios $\langle \eta = F_{\text{centrifugal}}/F_{\text{gravity}} \rangle = 0.64 \pm 0.33$ and $\langle \eta \rangle = 0.77 \pm 0.29$, respectively. Thus, the force ratios clearly indicate that on average our stars have quite under-critical rotations, even though the average $\Omega/\Omega_c$ ratios seem to suggest otherwise (see Zorec et al. 2016; Zorec 2023).
The obtained rotation rates can be translated into ratios of rotational frequencies that are in the interval $\langle v_r(\Omega/\Omega_c = 1.0)/v_r(\Omega/\Omega_c = 0.95) \rangle = 1.10 \pm 0.02$. This range can be taken as an indicator of the uncertainties affecting the rotational frequencies due to possible undue choices of the individual $\Omega/\Omega_c$ values.

## 9. Concluding remarks

A set of faint 19 B-type stars of which 18 are Be stars, observed by the satellite CoRoT in the fourth long run (LRA02), was studied by spectroscopy. Using spectra in the $\lambda\lambda 4000-4500$ Å obtained during a spectroscopic campaign with the VLT/FLAMES instruments at ESO. When necessary, the SEDs of Be stars were corrected for the veiling effect due to the circumstellar disk emission. All studied stars were assumed to be single stars or their SEDs had not been perturbed in the analyzed spectral range by the presence of possible stripped sdO/B components produced during the mass transfer phase in binary systems.

The fit of spectra with models of stellar atmospheres using our GIRFIT package based on the MINUIT minimization algorithm developed by CERN, produced apparent astrophysical parameters of stars, which were then corrected for the effects induced by the rapid rotation according to parameterized ratios $\Omega/\Omega_c$. The corrected astrophysical parameters enabled us to estimate true rotation frequencies that will be used in a following paper dealing with the analysis of non-radial pulsations derived from the light curves observed by CoRoT. We have concluded that it is difficult to decide whether the estimated effective temperatures in this work are systematically affected by the presence of sub-dwarf companions, because the expected mis-estimate of the effective temperatures depend on the actual characteristics of the companion sdO/B stars, and that they could be on the order of uncertainties related to the resolution and S/Ns of the spectra analyzed.

The masses of stars derived with the pncr astrophysical parameters and models of stellar evolution are in agreement with the masses measured in detached binary systems. Moreover, the obtained astrophysical parameters faithfully reproduce the stellar distances measured by Gaia.

We note however, that the inclination angle and other astrophysical parameters of the star N° 5 remain uncertain because in this case, the perturbation due to the secondary could be not negligible. Unfortunately, we do not have enough data to study this issue and decide whether some spectral features are genuine spectroscopic shell characteristics or they are signatures due to the presence a binary component.

Using the pncr astrophysical parameters lead us to conclude that the studied Be stars are rapid rotators, but the relation of the equatorial centrifugal force to the gravity indicate that they are far from being critical rotators. Critical or near-critical rotation is a necessary ingredient to favor massive mass ejections related with the non-radial pulsations (Kee et al. 2016). Nevertheless, the inferred $\Omega/\Omega_c$ ratios are large enough to favor the appearance of more or less extended regions in the stellar envelope and in the sub-photospheric regions that become unstable to convection (Clement 1979; Maeder et al. 2008; Cantiello et al. 2009). Since convection favors the setting out of differential rotation with a concomitant production of internal magnetic fields, magnetic dynamos, and cyclic activities (Zorec et al. 2011; Zorec 2023; Keszthelyi 2023), we may ask whether instabilities related with these magnetic fields could be the cause of significant upheavals

in the outer stellar structure, where the correlated non-radial pulsations and huge sporadic mass ejections could both be phenomena triggered by these disturbances. In the framework of this phenomenology, the study of non-radial pulsations in Be stars may acquire a new major area of interest centered on the study of the internal structure of rapidly rotating stars.

Spectra were obtained also in the region of the $H\alpha$ line to confirm the Be nature of objects in the chosen set of stars. For some objects, emission components could also be extracted in the $H\gamma$ and $H\delta$ lines. Using a simplified representation of the line emitting regions in a flat CD, we inferred indications on the physical properties of disks surrounding the observed Be stars. A non-LTE approach was employed to calculate the source function of lines and their opacity based on first principles that control the level populations of atoms. Entrusting the fit obtained of the observed $H\alpha$, $H\gamma$, and $H\delta$ emission lines profiles with model CD disks oriented according to the inclination angles derived from the astrophysical parameters, we have discussed the density distribution in the CDs. We thus found that the exponent of the density distribution in the radial direction of the circumstellar disk is a function of the distance, $R$.

Although the VDD models also predict exponents of the radial distribution of the density in the CDs that depend on the distance, our results point out on the necessity to carry out controls of the predicted exponents using detailed analysis of spectral lines that preferentially form in regions close to the central object. These regions are fed with material by sporadic huge mass ejections and variable stellar winds that provide the material to build up the CD. The spectral lines formed in these regions may reveal the actual physical properties of layers which make the transition from the star to the CD and where is organized the transport of the angular momentum that in fine determines the Keplerian rotation of disks in Be stars.

The obtained extents of the $H\alpha$ line emission formation regions are in agreement with similar quantities estimated elsewhere by interferometry in other Be stars. From the derived line opacity, we inferred the CD base densities that closely approach those estimated with much detailed models developed recently by other authors. Assuming that the snapshot image brought by our spectra of the line emission profiles can be attributed to a kind of stationary structure of the CD, the disk base density enabled us to obtain a rough estimate of the mass contained in the region of the line formation regions. These estimates are also in good agreement with other published estimates in the literature, where they were derived considering that the CD are subject to a dynamical evolution as well as that a quasi stationary picture actually results from several mass injection phenomena and to their dissipation.

**Acknowledgements.** This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This research has made use of the service of the ESO Science Archive Facility. We are thankful to the referee for his/her attentive reading of the manuscript and for his/her many suggestions and corrections that improved the presentation of our results.

## References

Arcos, C., Jones, C. E., Sigut, T. A. A., Kanaan, S., & Curé, M. 2017, ApJ, 842, 48

Arias, M. L. 2004, PhD thesis, National University of La Plata, Argentina

Arias, M. L., Zorec, J., & Frémat, Y. 2007, ASP Conf. Ser., 361, 419

Ashok, N. M., Bhatt, H. C., Kulkarni, P. V., & Joshi, S. C. 1984, MNRAS, 211, 471
Baade, D., Rivinius, T., Pigulski, A., et al. 2016, A&A, 588, A56

Baglin, A., Auvergne, M., Barge, P., et al. 2006, ESA SP, 1306, 33

Ballereau, D., Chauville, J., & Zorec, J. 1995, A&AS, 111, 423

Basri, G., & Batalha, C. 1990, ApJ, 363, 654

Bodensteiner, J., Shenar, T., & Sana, H. 2020, A&A, 641, A42

Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977

Burssens, S., Simón-Díaz, S., Bowman, D. M., et al. 2020, A&A, 639, A81

Cantiello, M., Langer, N., Brott, I., et al. 2009, A&A, 499, 279

Carciofi, A. C., & Bjorkman, J. E. 2006, ApJ, 639, 1081

Carciofi, A. C., & Bjorkman, J. E. 2008, ApJ, 684, 1374

Castelli, F., & Kurucz, R. L. 2003, IAU Symp., 210, A20

Castelli, F., Gratton, R. G., & Kurucz, R. L. 1997, A&A, 318, 841

Catanzaro, G. 2013, A&A, 550, A79

Chauville, J., Zorec, J., Ballereau, D., et al. 2001, A&A, 378, 861

Cidale, L. S., & Ringuelet, A. E. 1989, PASP, 101, 417

Clement, M. J. 1979, ApJ, 230, 230

Collins, George W., I. 1987, IAU Colloq., 92, 3

Cook, K. H., Alcock, C., Allsman, H. A., et al. 1995, ASP Conf. Ser., 83, 221

Cox, A. N. 2000, Allen's Astrophysical Quantities (UK: The University of London, The Athlone Press)

Cranmer, S. R. 2005, ApJ, 634, 585

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources (NASA/IPAC)

Dachs, J., Engels, D., & Kiehling, R. 1988, A&A, 194, 167

de Almeida, E. S. G., Meilland, A., Domiciano de Souza, A., et al. 2020, A&A, 636, A110

de Jager, C., & Neven, L. 1959, The Observatory, 79, 102

de Mink, S. E., Langer, N., Izzard, R. G., Sana, H., & de Koter, A. 2013, ApJ, 764, 166

de Wit, W. J., Lamers, H. J. G. L. M., Marquette, J. B., & Beaulieu, J. P. 2006, A&A, 456, 1027

Delaa, O., Stee, P., Meilland, A., et al. 2011, A&A, 529, A87

Deleuil, M., Moutou, C., Cabrera, J., et al. 2016, in The CoRoT Legacy Book: the Adventure of the Ultra High Precision Photometry from Space (France: EDP Sciences), 117

Drilling, J. S., Jeffery, C. S., Heber, U., Moehler, S., & Napiwotzki, R. 2013, A&A, 551, A31

El-Badry, K., Conroy, C., Quataert, E., et al. 2022, MNRAS, 516, 3602

Espinosa Lara, F., & Rieutord, M. 2011, A&A, 533, A43

Floquet, M., Hubert, A. M., Hirata, R., et al. 2000, A&A, 362, 1020

Flower, P. J. 1996, ApJ, 469, 355

Frémat, Y., Zorec, J., Hubert, A. M., & Floquet, M. 2005, A&A, 440, 305

Frémat, Y., Neiner, C., Hubert, A. M., et al. 2006, A&A, 451, 1053

Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1

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

Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1

Galazutdinov, G. A., Musaev, F. A., Krelowski, J., & Walker, G. A. H. 2000, PASP, 112, 648

George, C., Ekström, S., Granada, A., et al. 2013, A&A, 553, A24

Ghoreyshi, M. R., Carciofi, A. C., Rímulo, L. R., et al. 2018, MNRAS, 479, 2214

Ghoreyshi, M. R., Carciofi, A. C., Jones, C. E., et al. 2021, ApJ, 909, 149

Gkouvelis, L., Fabregat, J., Zorec, J., et al. 2016, A&A, 591, A140

Goss, K. J. F., Karoff, C., Chaplin, W. J., Elsworth, Y., & Stevens, I. R. 2011, MNRAS, 411, 162

Granada, A., Arias, M. L., & Cidale, L. S. 2010, AJ, 139, 1983

Gray, D. F. 2008, The Observation and Analysis of Stellar Photospheres (Cambridge: Cambridge University Press)

Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161

Hanuschik, R. W. 1989, Ap&SS, 161, 61

Hanuschik, R. W. 2003, A&A, 407, 1157

Harmanec, P. 1988, Bull. Astron. Inst. Czechoslovakia, 39, 329

Hastings, B., Langer, N., Wang, C., Schootemeijer, A., & Milone, A. P. 2021, A&A, 653, A144

Haubois, X., Carciofi, A. C., Rivinius, T., Okazaki, A. T., & Bjorkman, J. E. 2012, ApJ, 756, 156

Hawkins, F. M. 1961, MNRAS, 122, 285

Heber, U. 2009, ARA&A, 47, 211

Hernández, J., Calvet, N., Hartmann, L., et al. 2005, AJ, 129, 856

Horne, K., & Marsh, T. R. 1986, MNRAS, 218, 761

Huang, S.-S. 1972, ApJ, 171, 549

Huat, A.-L., Hubert, A.-M., Baudin, F., et al. 2009, A&A, 506, 95

Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875

Hubeny, I., & Mihalas, D. 2014, Theory of Stellar Atmospheres (Princeton: Princeton Series in Astrophysics)

Hubert, A. M., & Floquet, M. 1998, A&A, 335, 565

Hubert-Delplace, A. M., Hubert, H., Ballereau, D., & Chambon, M. T. 1982, IAU Symp., 98, 195

Hummel, W. 1994, A&A, 289, 458

Jackson, B. V., Buffington, A., Hick, P. P., et al. 2004, Sol. Phys., 225, 177

Jaschek, M., Slettebak, A., & Jaschek, C. 1981, Be Star Newslett., 4, 9

Jefferies, J. T. 1968, Spectral line formation (Waltham: Blaisdell Publishing Company)

Jenniskens, P., & Desert, F. X. 1994, A&AS, 106, 39

Kastner, J. H., & Mazzali, P. A. 1989, A&A, 210, 295

Kato, S. 1983, PASJ, 35, 249

Kee, N. D., Owocki, S., Townsend, R., & Müller, H. R. 2016, ASP Conf. Ser., 506, 47

Keller, S. C., Bessell, M. S., Cook, K. H., Geha, M., & Syphers, D. 2002, AJ, 124, 2039

Keszthelyi, Z. 2023, Galaxies, 11, 40

Kinman, T., & Castelli, F. 2002, A&A, 391, 1039

Klagyivik, P., Deeg, H. J., Cabrera, J., Csizmadia, S., & Almenara, J. M. 2017, A&A, 602, A117

Koubský, P., Hummel, C. A., Harmanec, P., et al. 2010, A&A, 517, A24

Kurfürst, P., Feldmeier, A., & Krtička, J. 2018, A&A, 613, A75

Kurtz, D. W., Shibahashi, H., Murphy, S. J., Bedding, T. R., & Bowman, D. M. 2015, MNRAS, 450, 3015

Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. Kurucz CD-ROM Cambridge (Cambridge: Smithsonian Astrophysical Observatory), 13

Labadie-Bartz, J., Carciofi, A. C., Henrique de Amorim, T., et al. 2022, AJ, 163, 226

Lanz, T., & Hubeny, I. 2003, ASP Conf. Ser., 288, 157

Lee, U., Osaki, Y., & Saio, H. 1991, MNRAS, 250, 432

Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2

Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Berlin: Springer Berlin Heidelberg)

Maeder, A., & Meynet, G. 2000, ARA&A, 38, 143

Maeder, A., Georgy, C., & Meynet, G. 2008, A&A, 479, L37

Marlborough, J. M. 1969, ApJ, 156, 135

Marr, K. C., Jones, C. E., Carciofi, A. C., et al. 2021, ApJ, 912, 76

Martayan, C., Frémat, Y., Hubert, A. M., et al. 2006, A&A, 452, 273

Martayan, C., Frémat, Y., Hubert, A. M., et al. 2007, A&A, 462, 683

McGill, M. A., Sigut, T. A. A., & Jones, C. E. 2013, ApJS, 204, 2

Mennickent, R. E., Pietrzyński, G., Gieren, W., & Szewczyk, O. 2002, A&A, 393, 887

Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., & Teller, E. 1953, J. Chem. Phys., 21, 1087

Miglio, A., Montalbán, J., & Dupret, M.-A. 2007, Commun. Asteroseismol., 151, 48

Mihalas, D. 1978, Stellar Atmospheres (San Francisco: W.H. Freeman and Company)

Millar, C. E., & Marlborough, J. M. 1998, ApJ, 494, 715

Millar, C. E., & Marlborough, J. M. 1999, ApJ, 516, 276

Moravveji, E., Townsend, R. H. D., Aerts, C., & Mathis, S. 2016, ApJ, 823, 130

Moujtahid, A. 1998, PhD thesis, Université Paris VI, France

Moujtahid, A., Zorec, J., Hubert, A. M., Garcia, A., & Burki, G. 1998, A&AS, 129, 289

Moujtahid, A., Zorec, J., & Hubert, A. M. 1999, A&A, 349, 151

Moujtahid, A., Zorec, J., & Hubert, A. M. 2000a, ASP Conf. Ser., 214, 506

Moujtahid, A., Zorec, J., & Hubert, A. M. 2000b, ASP Conf. Ser., 214, 510

Narita, S., Kiguchi, M., & Hayashi, C. 1994, PASJ, 46, 575

Neiner, C., Hubert, A. M., Frémat, Y., et al. 2003, A&A, 409, 275

Okazaki, A. T. 1991, PASJ, 43, 75

Okazaki, A. T. 1997, A&A, 318, 548

Okazaki, A. T. 2001, PASJ, 53, 119

Okazaki, A. T. 2007, ASP Conf. Ser., 361, 230

Osaki, Y. 1986, PASP, 98, 30

Osterbrock, D. E., Fulbright, J. P., Martel, A. R., et al. 1996, PASP, 108, 277

Packet, W. 1981, A&A, 102, 17

Pápicos, P. I., Tkachenko, A., Van Reeth, T., et al. 2017, A&A, 598, A74

Pedersen, M. G., Aerts, C., Pápicos, P. I., et al. 2021, Nat. Astron., 5, 715

Poeckert, R., & Marlborough, J. M. 1978a, ApJ, 220, 940

Poeckert, R., & Marlborough, J. M. 1978b, ApJS, 38, 229

Pogson, N. 1856, MNRAS, 17, 12

Pols, O. R., Cote, J., Waters, L. B. F. M., & Heise, J. 1991, A&A, 241, 419

Quirrenbach, A., Bjorkman, K. S., Bjorkman, J. E., et al. 1997, ApJ, 479, 477

Raddi, R., Drew, J. E., Fabregat, J., et al. 2013, MNRAS, 430, 2169

Ricker, G. R., Vanderspek, R., Winn, J., et al. 2016, SPIE Conf. Ser., 9904, 9904AB

Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3

Rivinius, T., Carciofi, A. C., & Martayan, C. 2013, A&ARv, 21, 69

Sackmann, I. J. 1970, A&A, 8, 76

Saio, H. 2013, in Lecture Notes in Physics, eds. M. Goupil, K. Belkacem, C. Neiner, F. Lignières, & J. J. Green (Berlin: Springer Verlag), 865, 159
Saio, H., Ekström, S., Mowlavi, N., et al. 2017, MNRAS, 467, 3864

Salmon, S. J. A. J., Montalbán, J., Reese, D. R., Dupret, M. A., & Eggenberger, P. 2014, A&A, 569, A18

Sebastian, D., Guenther, E. W., Schaffenroth, V., et al. 2012, A&A, 541, A34

Semaan, T., Hubert, A. M., Zorec, J., et al. 2013, A&A, 551, A130

Semaan, T., Hubert, A. M., Zorec, J., et al. 2018, A&A, 613, A70

Shao, Y., & Li, X.-D. 2014, ApJ, 796, 37

Sigut, T. A. A., & Jones, C. E. 2007, ApJ, 668, 481

Sigut, T. A. A., McGill, M. A., & Jones, C. E. 2009, ApJ, 699, 1973

Silaj, J., Jones, C. E., Tycner, C., Sigut, T. A. A., & Smith, A. D. 2010, ApJS, 187, 228

Smallwood, J. L., Martin, R. G., & Lubow, S. H. 2023, MNRAS, 520, 2952

Stehlé, C., & Hutcheon, R. 1999, A&AS, 140, 93

Stoeckley, T. R. 1968, MNRAS, 140, 141

Szewczuk, W., & Daszyńska-Daszkiewicz, J. 2017, MNRAS, 469, 13

Thomas, R. N. 1957, ApJ, 125, 260

Thomas, R. N. 1965, Some Aspects of Non-equilibrium Thermodynamics in the Presence of a Radiation field (Colorado: University of Colorado Press)

Torres, G. 2010, AJ, 140, 1158

Touhami, Y., Gies, D. R., Schaefer, G. H., et al. 2013, ApJ, 768, 128

Townsend, R. H. D., Owocki, S. P., & Howarth, I. D. 2004, MNRAS, 350, 189

Tycner, C., Hajian, A. R., Armstrong, J. T., et al. 2004, AJ, 127, 1194

Tycner, C., Lester, J. B., Hajian, A. R., et al. 2005, ApJ, 624, 359

Tycner, C., Gilbreath, G. C., Zavala, R. T., et al. 2006, AJ, 131, 2710

Tycner, C., Jones, C. E., Sigut, T. A. A., et al. 2008, ApJ, 689, 461

Vieira, R. G., Carciofi, A. C., & Bjorkman, J. E. 2015, MNRAS, 454, 2107

Vieira, R. G., Carciofi, A. C., Bjorkman, J. E., et al. 2017, MNRAS, 464, 3071

Vinicius, M. M. F., Zorec, J., Leister, N. V., & Levenhagen, R. S. 2006, A&A, 446, 643

Walczak, P., Fontes, C. J., Colgan, J., Kilcrease, D. P., & Guzik, J. A. 2015, A&A, 580, L9

Walker, G., Matthews, J., Kuschnig, R., et al. 2003, PASP, 115, 1023

Walker, G. A. H., Kuschnig, R., Matthews, J. M., et al. 2005a, ApJ, 635, L77

Walker, G. A. H., Kuschnig, R., Matthews, J. M., et al. 2005b, ApJ, 623, L145

Waters, L. B. F. M. 1986, A&A, 162, 121

Weiss, W. W., Aerts, C., Aigrain, S., et al. 2004, ESA SP, 538, 435

Weiss, W. W., Rucinski, S. M., Moffat, A. F. J., et al. 2014, PASP, 126, 573

Xiong, J., Liu, C., Li, J., et al. 2023, AJ, 165, 30

Zorec, J. 1986, Thèse d'État, Structure et rotation différentielle dans les étoiles B avec et sans émission (Université Paris VII, France)

Zorec, J. 2023, Galaxies, 11, 54

Zorec, J., & Briot, D. 1991, A&A, 245, 150

Zorec, J., & Royer, F. 2012, A&A, 537, A120

Zorec, J., Frémat, Y., & Cidale, L. 2005, A&A, 441, 235

Zorec, J., Arias, M. L., Cidale, L., & Ringuelet, A. E. 2007, A&A, 470, 239

Zorec, J., Cidale, L., Arias, M. L., et al. 2009, A&A, 501, 297

Zorec, J., Frémat, Y., Domiciano de Souza, A., et al. 2011, A&A, 526, A87

Zorec, J., Frémat, Y., Domiciano de Souza, A., et al. 2016, A&A, 595, A132
## Appendix A: Table containing the information on FLAMES/GIRAFFE spectra and conditions of their exposure time.

<table><caption>**Table A.1.** Log of the stars observed with FLAMES/GIRAFFE</caption><thead><tr><th>CoRoT-ID</th><th>N°</th><th>RA</th><th>DEC</th><th>Run</th><th>Grat</th><th>F</th><th>EXP</th><th>DATE-OBS</th><th>AM</th><th>S/N</th><th>Moon</th></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>sec.</td><td></td><td></td><td></td><td></td></tr></thead><tbody><tr><td rowspan="2">103000272</td><td rowspan="2">1</td><td rowspan="2">06 50 39.09</td><td rowspan="2">-03 18 14.37</td><td rowspan="2">LRA2-b</td><td>LR2</td><td>1</td><td>3685</td><td>2010-11-23T06:11:00</td><td>1.13</td><td>170</td><td>98%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-11-23T07:25:16</td><td>1.07</td><td>140</td><td>98%</td></tr><tr><td rowspan="2">103032255</td><td rowspan="2">2</td><td rowspan="2">06 51 23.79</td><td rowspan="2">-03 14 46.51</td><td rowspan="2">LRA2-c</td><td>LR2</td><td>1</td><td>3685</td><td>2010-12-24T04:18:40</td><td>1.12</td><td>170</td><td>90%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-12-24T05:32:01</td><td>1.07</td><td>160</td><td>90%</td></tr><tr><td rowspan="2">110655185</td><td rowspan="2">3</td><td rowspan="2">06 49 06.10</td><td rowspan="2">-05 00 21.86</td><td rowspan="2">LRA2-m</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-03T02:18:40</td><td>1.29</td><td>65</td><td>B/H</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-12-26T04:20:50</td><td>1.09</td><td>-</td><td>72%</td></tr><tr><td rowspan="2">110655384</td><td rowspan="2">4</td><td rowspan="2">06 49 06.80</td><td rowspan="2">-05 20 41.77</td><td rowspan="2">LRA2-m</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-03T02:18:40</td><td>1.29</td><td>90</td><td>B/H</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-12-26T04:20:50</td><td>1.09</td><td>80</td><td>72%</td></tr><tr><td rowspan="2">110655437</td><td rowspan="2">5</td><td rowspan="2">06 49 06.98</td><td rowspan="2">-05 02 27.30</td><td rowspan="2">LRA2-m</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-03T02:18:40</td><td>1.29</td><td>65</td><td>B/H</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-12-26T04:20:50</td><td>1.09</td><td>75</td><td>72%</td></tr><tr><td rowspan="2">110662847</td><td rowspan="2">6</td><td rowspan="2">06 50 10.73</td><td rowspan="2">-03 20 55.79</td><td rowspan="2">LRA2-b</td><td>LR2</td><td>1</td><td>3685</td><td>2010-11-23T06:11:00</td><td>1.13</td><td>85</td><td>98%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-11-23T07:25:16</td><td>1.07</td><td>90</td><td>98%</td></tr><tr><td rowspan="2">110663174</td><td rowspan="2">7</td><td rowspan="2">06 50 11.94</td><td rowspan="2">-04 40 53.22</td><td rowspan="2">LRA2-k</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-20T04:25:44</td><td>1.09</td><td>62</td><td>100%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2011-01-22T02:50:53</td><td>1.08</td><td>450</td><td>92%</td></tr><tr><td rowspan="2">110663880</td><td rowspan="2">8</td><td rowspan="2">06 50 14.36</td><td rowspan="2">-03 23 59.50</td><td rowspan="2">LRA2-b</td><td>LR2</td><td>1</td><td>3685</td><td>2010-11-23T06:11:00</td><td>1.13</td><td>170</td><td>98%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-11-23T07:25:16</td><td>1.07</td><td>140</td><td>98%</td></tr><tr><td rowspan="2">110672515</td><td rowspan="2">9</td><td rowspan="2">06 51 14.24</td><td rowspan="2">-05 23 51.02</td><td rowspan="2">LRA2-q</td><td>LR2</td><td>1</td><td>3685</td><td>2010-11-27T06:15:00</td><td>1.09</td><td>59</td><td>69%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-11-27T07:28:53</td><td>1.06</td><td>60</td><td>69%</td></tr><tr><td rowspan="2">110681176</td><td rowspan="2">10</td><td rowspan="2">06 52 08.32</td><td rowspan="2">-03 37 52.33</td><td rowspan="2">LRA2-f</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-17T01:36:53</td><td>1.27</td><td>300</td><td>91%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2011-01-17T02:48:54</td><td>1.11</td><td>250</td><td>91%</td></tr><tr><td rowspan="2">110688151</td><td rowspan="2">11</td><td rowspan="2">06 53 00.94</td><td rowspan="2">-04 40 18.17</td><td rowspan="2">LRA2-l</td><td>LR2</td><td>1</td><td>3685</td><td>2010-12-24T06:16:48</td><td>1.09</td><td>350</td><td>90%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-12-24T07:28:55</td><td>1.24</td><td>250</td><td>90%</td></tr><tr><td rowspan="2">110747131</td><td rowspan="2">12</td><td rowspan="2">06 50 27.04</td><td rowspan="2">-05 26 36.86</td><td rowspan="2">LRA2-q</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-27T06:15:00</td><td>1.09</td><td>91</td><td>69%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2011-01-27T07:28:53</td><td>1.06</td><td>60</td><td>69%</td></tr><tr><td rowspan="2">110751872</td><td rowspan="2">13</td><td rowspan="2">06 51 17.51</td><td rowspan="2">-04 36 59.97</td><td rowspan="2">LRA2-k</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-20T04:25:44</td><td>1.09</td><td>450</td><td>100%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2011-01-22T02:50:53</td><td>1.08</td><td>450</td><td>92%</td></tr><tr><td rowspan="2">110751876</td><td rowspan="2">14</td><td rowspan="2">06 51 17.52</td><td rowspan="2">-03 19 29.71</td><td rowspan="2">LRA2-c</td><td>LR2</td><td>1</td><td>3685</td><td>2010-12-23T04:18:40</td><td>1.12</td><td>450</td><td>90%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-12-24T05:32:01</td><td>1.07</td><td>300</td><td>90%</td></tr><tr><td rowspan="2">110752156</td><td rowspan="2">15</td><td rowspan="2">06 51 18.19</td><td rowspan="2">-04 07 40.62</td><td rowspan="2">LRA2-h</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-21T02:36:14</td><td>1.09</td><td>60</td><td>98%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2011-01-21T03:58:07</td><td>1.07</td><td>50</td><td>98%</td></tr><tr><td rowspan="2">110827583</td><td rowspan="2">16</td><td rowspan="2">06 48 58.53</td><td rowspan="2">-00 13 06.25</td><td rowspan="2">LRA2-m</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-03T02:18:40</td><td>1.29</td><td>130</td><td>B/H</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-12-26T04:20:50</td><td>1.09</td><td>100</td><td>72%</td></tr></tbody></table>
<table><caption>**Table A.1. Continued.**</caption><thead><tr><th>CoRoT-ID</th><th>N°</th><th>RA</th><th>DEC</th><th>Run</th><th>Grat</th><th>F</th><th>EXP</th><th>DATE-OBS</th><th>AM</th><th>S/N</th><th>Moon</th></tr></thead><tbody><tr><td rowspan="2">300002611</td><td rowspan="2">17</td><td rowspan="2">06 48 19.69</td><td rowspan="2">-03 43 27.41</td><td rowspan="2">LRA2-d</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-16T06:08:08</td><td>1.29</td><td>65</td><td>B/H</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-11-23T08:03:29</td><td>1.09</td><td>50</td><td>98%</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>2</td><td>780</td><td>2010-11-23T08:17:23</td><td>1.10</td><td>-</td><td>98%</td></tr><tr><td rowspan="2">300002834</td><td rowspan="2">18</td><td rowspan="2">06 48 25.86</td><td rowspan="2">-03 26 21.14</td><td rowspan="2">LRA2-a</td><td>LR2</td><td>1</td><td>3685</td><td>2010-12-26T04:57:58</td><td>1.07</td><td>70</td><td>72%</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2011-01-10T03:07:30</td><td>1.12</td><td>60</td><td>B/H</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>2</td><td>780</td><td>2011-01-10T03:21:24</td><td>1.10</td><td>-</td><td>B/H</td></tr><tr><td rowspan="2">300003290</td><td rowspan="2">19</td><td rowspan="2">06 48 36.90</td><td rowspan="2">-03 53 45.79</td><td rowspan="2">LRA2-d</td><td>LR2</td><td>1</td><td>3685</td><td>2011-01-16T06:08:08</td><td>1.29</td><td>60</td><td>B/H</td></tr><tr><td>LR6</td><td>1</td><td>780</td><td>2010-11-23T08:03:29</td><td>1.09</td><td>50</td><td>98%</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td>2</td><td>780</td><td>2010-11-23T08:17:23</td><td>1.10</td><td>-</td><td>98%</td></tr></tbody></table>

Notes: F = ; Grant = ; Run = ; RA = right ascension' DEC = declination; EXP = exposure time in seconds; AM = air mass; S/N = signal-to-noise ratio ; Moon = Moon's illumination ; B/H = moon below Horizon
## Appendix B: Comparison between an observed and the corresponding veiling corrected spectrum

According to Semaan et al. (2013), the expression to correct the spectra for the veiling effect is

$$F^*(\lambda)/F_c^*(\lambda) = \left\{ \frac{[F_{\text{obs}}(\lambda)/F_c^{\text{obs}}(\lambda)][1+r(\lambda)] - r(\lambda)}{r(\lambda) \cdot \bar{E}_c(\lambda)/[F_c^*(\lambda)A_c(\lambda)]} \right\}, \quad (\text{B.1})$$

where $F^{\text{obs}}(\lambda)$ is the observed line spectrum; $F_c^{\text{obs}}(\lambda)$ is the continuum spectrum emitted by the star+disk system; $F^*(\lambda)$ and $F_c^*(\lambda)$ are the line and continuum spectra due to the stellar photosphere, respectively. The term “1 + r” is called the “veiling factor.” $E_c(\lambda)$ and $A_c(\lambda)$ respectively represent the amount of continuum emission and absorption due to the CD that can be estimated using the empirical correlations in (Ballereau et al. 1995). The veiling factor in Eq. (B.1) is a function of $\lambda$ that varies by $\Delta r/r \gtrsim 0.04$ from $\lambda = 4000$ Å to $\lambda = 4500$ Å as it can be inferred from relations given in Moujtahid et al. (1999, 2000a) and using the required parameters calculated in Zorec et al. (2016). The term $r(\lambda)$ roughly behaves as $r(\lambda) = a \times (\lambda/\lambda_{\text{Hy}})^3 + b$. Frémat et al. (2006) showed that by neglecting the veiling effect in Be stars displaying strong line emissions, the determination of astrophysical parameters ($T_{\text{eff}}$, log g) may suffer from substantial uncertainties. The veiling factors r at $\lambda = 4500$ Å determined using the empirical correlations in Ballereau et al. (1995) and the emission intensities $I_{\text{Hy}}$ in the Hγ line are given in Table B.1. The emission intensity in the Hγ line is defined as $I_{\text{Hy}} = W_{\text{Hy}}^{\text{em}} \times \{[F_{\text{Hy}}^{\text{c}}(T_{\text{eff}}, \log g)]/[F_{\text{Hy}}^{\text{c}}(22500, 4.0)]\}$, where $W_{\text{Hy}}^{\text{em}}$ represents the equivalent width of the emission component of the Hγ line and $F_{\text{Hy}}^{\text{c}}$ is the flux of the continuum spectrum at Hγ. In Table B.1, stars N° 1, 6, 12, and 16 have very low emission and accordingly it is $r \approx 0$. An example of the effect introduced on spectra by the correction by the veiling effect is shown in Fig. B.1.

As an example of spectra before and after correction made for the veiling effect are shown in Fig. B.1. They correspond to the Be star ID 300002611 (N° 17). The correction is made over the entire spectral range of $\lambda\lambda 4000 - 4500$ Å. The spectra in Fig. B.1 are shown to be arbitrarily shifted in the ordinates. In this figure is also shown the difference between the observed and veiling corrected spectra to appreciate the change in the equivalent widths of spectral lines and the almost negligible variation introduced on the energy distribution in the treated wavelength interval.

<figure><img src="image_14.png" /><figcaption>**Fig. B.1.** Comparison between normalized observed and veiling corrected spectra of the program star N° 17 (ID 300002611) in the wavelength interval $\lambda\lambda 4000 - 4500$ Å.</figcaption></figure>

<table><caption>**Table B.1.** Equivalent width and normalized flux of the emission component in the Hγ line, and the veiling factor $1 + r$.</caption><thead><tr><th>Star</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>W<sup>em</sup><sub>Hy</sub></td><td>0.00</td><td>0.31</td><td>0.74</td><td>0.38</td><td>0.51</td><td>0.00</td><td>0.77</td><td>1.49</td><td>1.06</td><td>0.06</td></tr><tr><td>I<sub>Hy</sub></td><td>0.00</td><td>0.12</td><td>0.22</td><td>0.14</td><td>0.23</td><td>0.00</td><td>0.32</td><td>1.83</td><td>0.47</td><td>0.01</td></tr><tr><td>1 + r</td><td>1.00</td><td>1.04</td><td>1.14</td><td>1.05</td><td>1.08</td><td>1.00</td><td>1.14</td><td>1.36</td><td>1.23</td><td>1.00</td></tr></tbody></table><table><thead><tr><th>Star</th><th>11</th><th>12</th><th>13</th><th>14</th><th>15</th><th>16</th><th>17</th><th>18</th><th>19</th></tr></thead><tbody><tr><td>W<sup>em</sup><sub>Hy</sub></td><td>1.24</td><td>0.00</td><td>0.37</td><td>0.96</td><td>0.20</td><td>0.00</td><td>1.89</td><td>0.31</td><td>1.36</td></tr><tr><td>I<sub>Hy</sub></td><td>0.53</td><td>0.00</td><td>0.15</td><td>0.96</td><td>0.18</td><td>0.00</td><td>0.59</td><td>0.12</td><td>1.73</td></tr><tr><td>1 + r</td><td>1.28</td><td>1.00</td><td>1.05</td><td>1.20</td><td>1.02</td><td>1.00</td><td>1.50</td><td>1.04</td><td>1.32</td></tr></tbody></table><p>Equivalent widths W<sup>em</sup><sub>Hy</sub> and intensities I<sub>Hy</sub> are given in Å</p>

## Appendix C: Spectra in the blue region and the obtained best fits

## Appendix D: Astrophysical parameters and rotational frequencies corrected for rotational effects of the observed Be stars assuming the angular velocity ratio $\Omega/\Omega_c = 0.95$
<figure><img src="image_15.png" /><figcaption>Fig. C.1. Observed spectra corrected for veiling effects of the program stars from N° 1 to 6 in the wavelength interval λλ4000 – 4500 Å.</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>Fig. C.2. Observed spectra corrected for veiling effects of the program stars from N° 7 to 12 in the wavelength interval λλ 4000 – 4500 Å.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>Fig. C.3. Observed spectra corrected for veiling effects of the program stars from N° 13 to 18 in the wavelength interval λλ 4000 – 4500 Å.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>Fig. C.4. Observed spectra corrected for veiling effects of the program star N° 19 in the wavelength interval λλ 4000–4500 Å.</figcaption></figure>

<table><caption>**Table D.1.** Averaged astrophysical parameters and rotational frequencies corrected for rotational effects of the observed Be stars assuming the angular velocity ratio $\Omega/\Omega_c = 0.95$</caption><thead><tr><th>Star</th><th>⟨T<sub>eff</sub>⟩ ± σ<br>K</th><th>⟨log g⟩ ± σ<br>dex</th><th>⟨Vsin i⟩ ± σ<br>km s<sup>-1</sup></th><th>⟨log L/L<sub>☉</sub>⟩ ± σ<br>dex</th><th>⟨M/M<sub>☉</sub>⟩ ± σ</th><th>⟨V<sub>c</sub>⟩ ± σ<br>km s<sup>-1</sup></th><th>⟨i⟩ ± σ<br>deg.</th><th>⟨t/t<sub>MS</sub>⟩ ± σ</th><th>⟨age⟩ ± σ<br>yr</th><th>⟨v<sub>r</sub>⟩ ± σ<br>c/d</th></tr></thead><tbody><tr><td>1</td><td>11980 ± 590</td><td>3.92 ± 0.15</td><td>134 ± 13</td><td>2.209 ± 0.085</td><td>2.90 ± 0.35</td><td>368 ± 25</td><td>26 ± 9</td><td>0.366 ± 0.256</td><td>2.11E+08 ± 8.01E+07</td><td>1.61 ± 0.25</td></tr><tr><td>2</td><td>15500 ± 880</td><td>3.71 ± 0.14</td><td>254 ± 25</td><td>3.167 ± 0.108</td><td>4.79 ± 0.47</td><td>366 ± 21</td><td>58 ± 18</td><td>0.703 ± 0.211</td><td>9.36E+07 ± 2.74E+07</td><td>0.89 ± 0.15</td></tr><tr><td>3</td><td>12620 ± 720</td><td>3.72 ± 0.15</td><td>234 ± 22</td><td>2.674 ± 0.110</td><td>3.26 ± 0.38</td><td>341 ± 21</td><td>57 ± 16</td><td>0.479 ± 0.264</td><td>1.84E+08 ± 7.22E+07</td><td>0.97 ± 0.17</td></tr><tr><td>4</td><td>15850 ± 840</td><td>4.40 ± 0.12</td><td>392 ± 36</td><td>2.428 ± 0.103</td><td>3.13 ± 0.68</td><td>513 ± 28</td><td>68 ± 14</td><td>0.199 ± 0.275</td><td>3.66E+09 ± 1.15E+09</td><td>3.12 ± 0.49</td></tr><tr><td>5</td><td>17740 ± 810</td><td>3.83 ± 0.12</td><td>258 ± 23</td><td>3.321 ± 0.093</td><td>5.19 ± 0.20</td><td>408 ± 21</td><td>50 ± 14</td><td>0.292 ± 0.225</td><td>4.68E+07 ± 2.35E+07</td><td>1.10 ± 0.17</td></tr><tr><td>6</td><td>18030 ± 880</td><td>3.65 ± 0.12</td><td>282 ± 24</td><td>3.652 ± 0.106</td><td>5.92 ± 0.81</td><td>378 ± 19</td><td>67 ± 14</td><td>0.533 ± 0.313</td><td>4.43E+07 ± 1.82E+07</td><td>0.71 ± 0.12</td></tr><tr><td>7</td><td>15880 ± 890</td><td>4.02 ± 0.15</td><td>296 ± 29</td><td>2.840 ± 0.112</td><td>4.13 ± 0.29</td><td>428 ± 27</td><td>58 ± 10</td><td>0.214 ± 0.227</td><td>5.85E+07 ± 4.01E+07</td><td>1.62 ± 0.30</td></tr><tr><td>8</td><td>24140 ± 1220</td><td>3.86 ± 0.16</td><td>214 ± 39</td><td>4.005 ± 0.115</td><td>9.05 ± 0.89</td><td>473 ± 32</td><td>33 ± 14</td><td>0.394 ± 0.297</td><td>1.60E+07 ± 7.73E+06</td><td>1.07 ± 0.20</td></tr><tr><td>9</td><td>16700 ± 940</td><td>3.58 ± 0.11</td><td>276 ± 23</td><td>3.597 ± 0.099</td><td>5.56 ± 0.41</td><td>355 ± 18</td><td>73 ± 14</td><td>0.657 ± 0.213</td><td>6.85E+07 ± 1.83E+07</td><td>0.61 ± 0.10</td></tr><tr><td>10</td><td>11790 ± 640</td><td>3.73 ± 0.14</td><td>230 ± 22</td><td>2.493 ± 0.105</td><td>3.08 ± 0.13</td><td>335 ± 20</td><td>57 ± 16</td><td>0.606 ± 0.196</td><td>2.69E+08 ± 6.57E+07</td><td>1.02 ± 0.17</td></tr><tr><td>11</td><td>16280 ± 880</td><td>4.15 ± 0.14</td><td>327 ± 32</td><td>2.755 ± 0.115</td><td>4.10 ± 0.30</td><td>460 ± 29</td><td>60 ± 14</td><td>0.097 ± 0.150</td><td>4.21E+07 ± 3.36E+07</td><td>2.06 ± 0.38</td></tr><tr><td>12</td><td>12930 ± 680</td><td>3.70 ± 0.15</td><td>217 ± 19</td><td>2.714 ± 0.106</td><td>3.41 ± 0.47</td><td>341 ± 21</td><td>51 ± 12</td><td>0.504 ± 0.296</td><td>1.62E+08 ± 6.26E+07</td><td>0.96 ± 0.16</td></tr><tr><td>13</td><td>15580 ± 850</td><td>3.96 ± 0.15</td><td>286 ± 28</td><td>2.867 ± 0.113</td><td>4.11 ± 0.25</td><td>412 ± 26</td><td>58 ± 18</td><td>0.304 ± 0.243</td><td>7.02E+07 ± 4.21E+07</td><td>1.43 ± 0.27</td></tr><tr><td>14</td><td>24400 ± 1240</td><td>4.12 ± 0.15</td><td>340 ± 38</td><td>3.731 ± 0.121</td><td>8.14 ± 1.34</td><td>541 ± 36</td><td>50 ± 18</td><td>0.057 ± 0.147</td><td>8.96E+06 ± 6.64E+06</td><td>1.78 ± 0.32</td></tr><tr><td>15</td><td>24450 ± 1610</td><td>3.72 ± 0.14</td><td>305 ± 34</td><td>4.256 ± 0.161</td><td>9.74 ± 1.16</td><td>449 ± 25</td><td>57 ± 15</td><td>0.483 ± 0.309</td><td>1.75E+07 ± 7.51E+06</td><td>0.78 ± 0.18</td></tr><tr><td>16</td><td>11800 ± 620</td><td>4.08 ± 0.15</td><td>235 ± 23</td><td>2.030 ± 0.108</td><td>2.61 ± 0.22</td><td>394 ± 26</td><td>46 ± 12</td><td>0.158 ± 0.197</td><td>1.89E+08 ± 1.11E+08</td><td>2.08 ± 0.36</td></tr><tr><td>17</td><td>14150 ± 760</td><td>4.23 ± 0.14</td><td>362 ± 32</td><td>2.386 ± 0.111</td><td>3.11 ± 0.17</td><td>454 ± 28</td><td>74 ± 14</td><td>0.015 ± 0.049</td><td>1.09E+09 ± 0.92E+09</td><td>2.37 ± 0.40</td></tr><tr><td>18</td><td>15200 ± 820</td><td>3.72 ± 0.12</td><td>249 ± 24</td><td>3.115 ± 0.100</td><td>4.58 ± 0.49</td><td>365 ± 19</td><td>56 ± 16</td><td>0.651 ± 0.228</td><td>9.77E+07 ± 3.20E+07</td><td>0.90 ± 0.14</td></tr><tr><td>19</td><td>17860 ± 870</td><td>3.49 ± 0.16</td><td>113 ± 13</td><td>3.708 ± 0.116</td><td>6.07 ± 0.81</td><td>351 ± 23</td><td>23 ± 9</td><td>0.625 ± 0.297</td><td>4.78E+07 ± 1.70E+07</td><td>0.61 ± 0.10</td></tr></tbody></table>

<table><caption>**Table D.2.** Modes of astrophysical parameters and rotational frequencies corrected for rotational effects of the observed Be stars assuming the angular velocity ratio $\Omega/\Omega_c = 0.95$</caption><thead><tr><th>Star</th><th>⟨T<sub>eff</sub>⟩ ± σ<br>K</th><th>⟨log g⟩ ± σ<br>dex</th><th>⟨Vsin i⟩ ± σ<br>km s<sup>-1</sup></th><th>⟨log L/L<sub>☉</sub>⟩ ± σ<br>dex</th><th>⟨M/M<sub>☉</sub>⟩ ± σ</th><th>⟨V<sub>c</sub>⟩ ± σ<br>km s<sup>-1</sup></th><th>⟨i⟩ ± σ<br>deg.</th><th>⟨t/t<sub>MS</sub>⟩ ± σ</th><th>⟨age⟩ ± σ<br>yr</th><th>⟨v<sub>r</sub>⟩ ± σ<br>c/d</th></tr></thead><tbody><tr><td>1</td><td>11980 ± 70</td><td>3.98 ± 0.02</td><td>134 ± 1</td><td>2.208 ± 0.010</td><td>3.12 ± 0.24</td><td>378 ± 3</td><td>24 ± 0</td><td>0.581 ± 0.230</td><td>2.51E+08 ± 8.17E+06</td><td>1.65 ± 0.03</td></tr><tr><td>2</td><td>15500 ± 100</td><td>3.76 ± 0.02</td><td>254 ± 2</td><td>3.167 ± 0.013</td><td>5.14 ± 0.06</td><td>366 ± 2</td><td>55 ± 1</td><td>0.869 ± 0.018</td><td>9.81E+07 ± 2.43E+06</td><td>0.88 ± 0.02</td></tr><tr><td>3</td><td>12620 ± 80</td><td>3.78 ± 0.02</td><td>234 ± 2</td><td>2.717 ± 0.013</td><td>3.10 ± 0.55</td><td>349 ± 2</td><td>53 ± 1</td><td>0.431 ± 0.109</td><td>1.72E+08 ± 1.40E+07</td><td>0.95 ± 0.02</td></tr><tr><td>4</td><td>15850 ± 100</td><td>4.45 ± 0.01</td><td>392 ± 4</td><td>2.428 ± 0.012</td><td>3.12 ± 0.20</td><td>525 ± 3</td><td>62 ± 1</td><td>0.034 ± 0.132</td><td>6.28E+09 ± 2.58E+07</td><td>3.11 ± 0.05</td></tr><tr><td>5</td><td>17740 ± 90</td><td>3.88 ± 0.01</td><td>258 ± 2</td><td>3.321 ± 0.011</td><td>5.29 ± 0.21</td><td>416 ± 2</td><td>47 ± 0</td><td>0.282 ± 0.093</td><td>4.31E+07 ± 4.56E+06</td><td>1.12 ± 0.02</td></tr><tr><td>6</td><td>18020 ± 100</td><td>3.70 ± 0.01</td><td>282 ± 2</td><td>3.651 ± 0.012</td><td>5.46 ± 0.17</td><td>386 ± 2</td><td>61 ± 1</td><td>0.288 ± 0.085</td><td>3.36E+07 ± 5.72E+06</td><td>0.73 ± 0.01</td></tr><tr><td>7</td><td>15890 ± 110</td><td>4.08 ± 0.02</td><td>296 ± 3</td><td>2.840 ± 0.013</td><td>4.05 ± 0.13</td><td>439 ± 3</td><td>54 ± 1</td><td>0.109 ± 0.084</td><td>2.66E+07 ± 8.33E+06</td><td>1.66 ± 0.03</td></tr><tr><td>8</td><td>24140 ± 140</td><td>3.92 ± 0.02</td><td>214 ± 4</td><td>4.050 ± 0.013</td><td>9.15 ± 0.42</td><td>486 ± 3</td><td>33 ± 0</td><td>0.436 ± 0.049</td><td>1.94E+07 ± 1.10E+06</td><td>1.05 ± 0.02</td></tr><tr><td>9</td><td>17090 ± 100</td><td>3.62 ± 0.01</td><td>286 ± 2</td><td>3.597 ± 0.011</td><td>5.62 ± 0.35</td><td>355 ± 2</td><td>89 ± 0</td><td>0.604 ± 0.100</td><td>6.48E+07 ± 5.28E+06</td><td>0.64 ± 0.01</td></tr><tr><td>10</td><td>11790 ± 80</td><td>3.79 ± 0.02</td><td>230 ± 2</td><td>2.492 ± 0.012</td><td>3.10 ± 0.08</td><td>335 ± 2</td><td>54 ± 1</td><td>0.624 ± 0.044</td><td>2.84E+08 ± 7.76E+06</td><td>1.02 ± 0.02</td></tr><tr><td>11</td><td>16280 ± 100</td><td>4.21 ± 0.02</td><td>327 ± 3</td><td>2.754 ± 0.013</td><td>4.30 ± 1.00</td><td>460 ± 3</td><td>58 ± 1</td><td>0.007 ± 0.016</td><td>5.09E+06 ± 3.68E+06</td><td>2.02 ± 0.04</td></tr><tr><td>12</td><td>12930 ± 80</td><td>3.76 ± 0.02</td><td>217 ± 2</td><td>2.714 ± 0.013</td><td>3.10 ± 0.28</td><td>350 ± 2</td><td>51 ± 0</td><td>0.254 ± 0.133</td><td>1.25E+08 ± 1.52E+07</td><td>1.00 ± 0.02</td></tr><tr><td>13</td><td>15580 ± 100</td><td>4.02 ± 0.02</td><td>286 ± 3</td><td>2.866 ± 0.013</td><td>4.17 ± 0.14</td><td>412 ± 3</td><td>54 ± 1</td><td>0.252 ± 0.078</td><td>5.51E+07 ± 7.99E+06</td><td>1.45 ± 0.03</td></tr><tr><td>14</td><td>24400 ± 150</td><td>4.18 ± 0.02</td><td>340 ± 4</td><td>3.730 ± 0.014</td><td>7.52 ± 0.05</td><td>555 ± 4</td><td>46 ± 1</td><td>0.000 ± 0.000</td><td>8.95E+06 ± 6.60E+06</td><td>1.79 ± 0.04</td></tr><tr><td>15</td><td>24450 ± 180</td><td>3.78 ± 0.02</td><td>305 ± 3</td><td>4.254 ± 0.019</td><td>9.43 ± 0.66</td><td>449 ± 3</td><td>50 ± 1</td><td>0.316 ± 0.127</td><td>1.66E+07 ± 1.62E+06</td><td>0.79 ± 0.02</td></tr><tr><td>16</td><td>11800 ± 70</td><td>4.13 ± 0.02</td><td>235 ± 2</td><td>2.028 ± 0.013</td><td>2.57 ± 0.05</td><td>404 ± 3</td><td>44 ± 0</td><td>0.017 ± 0.045</td><td>3.05E+07 ± 2.31E+07</td><td>2.16 ± 0.04</td></tr><tr><td>17</td><td>14150 ± 90</td><td>4.29 ± 0.02</td><td>362 ± 3</td><td>2.430 ± 0.013</td><td>3.10 ± 0.06</td><td>466 ± 3</td><td>89 ± 0</td><td>0.000 ± 0.000</td><td>1.09E+09 ± 0.92E+09</td><td>2.21 ± 0.04</td></tr><tr><td>18</td><td>14880 ± 100</td><td>3.76 ± 0.01</td><td>249 ± 2</td><td>3.115 ± 0.012</td><td>4.20 ± 0.15</td><td>365 ± 2</td><td>53 ± 1</td><td>0.556 ± 0.080</td><td>1.10E+08 ± 4.77E+06</td><td>0.87 ± 0.02</td></tr><tr><td>19</td><td>17860 ± 100</td><td>3.55 ± 0.02</td><td>113 ± 1</td><td>3.755 ± 0.014</td><td>5.65 ± 0.28</td><td>360 ± 2</td><td>23 ± 0</td><td>0.405 ± 0.143</td><td>4.18E+07 ± 4.75E+06</td><td>0.59 ± 0.01</td></tr></tbody></table>