# Benchmarking the fundamental parameters of Ap stars with optical long-baseline interferometric measurements

K. Perraut, M. Cunha, A. Romanovskaya, D. Shulyak, T. Ryabchikova, V.

Hocdé, N. Nardetto, D. Mourard, A. Meilland, F. Morand, et al.

► To cite this version:

K. Perraut, M. Cunha, A. Romanovskaya, D. Shulyak, T. Ryabchikova, et al., Benchmarking the fundamental parameters of Ap stars with optical long-baseline interferometric measurements. Astronomy and Astrophysics - A&A, 2020, 642, pp.A101. 10.1051/0004-6361/202038753. hal-02963288

HAL Id: hal-02963288

https://hal.science/hal-02963288

Submitted on 10 Oct 2020

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

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

<figure><img src="image_1.png" /></figure>
# Benchmarking the fundamental parameters of Ap stars with optical long-baseline interferometric measurements*

K. Perraut<fnref n="1" />, M. Cunha<fnref n="2" />,<fnref n="3" />, A. Romanovskaya<fnref n="4" />, D. Shulyak<fnref n="5" />, T. Ryabchikova<fnref n="4" />, V. Hocdé<fnref n="6" />, N. Nardetto<fnref n="6" />, D. Mourard<fnref n="6" />, A. Meilland<fnref n="6" />, F. Morand<fnref n="6" />, I. Tallon-Bosc<fnref n="7" />, C. Farrington<fnref n="8" />, and C. Lanthermann<fnref n="9" />

Received 25 June 2020 / Accepted 31 August 2020

## ABSTRACT

**Context.** The variety of physical processes at play in chemically peculiar stars makes it difficult to determine their fundamental parameters. In particular, for the magnetic ones, called Ap stars, the strong magnetic fields and the induced spotted stellar surfaces may lead to biased effective temperatures when these values are derived through spectro-photometry.

**Aims.** We propose to benefit from the exquisite angular resolution provided by long-baseline interferometry in the visible to determine the accurate angular diameters of a number of Ap stars, and thus estimate their radii by a method that is as independent as possible of atmospheric models.

**Methods.** We used the visible spectrograph VEGA at the CHARA interferometric array to complete the sample of Ap stars currently observable with this technique. We estimated the angular diameter and radius of six new targets. We estimated their bolometric flux based solely on observational spectroscopic and photometric data to derive nearly model-independent luminosities and effective temperatures.

**Results.** We extend to 14 the number of Ap stars for which interferometric angular diameters have been measured. The fundamental parameters we derived for the complete Ap sample are compared with those obtained through a self-consistent spectroscopic analysis. Based on a model fitting approach of high-resolution spectra and spectro-photometric observations over a wide wavelength range, this method takes into account the anomalous chemical composition of the atmospheres and the inhomogeneous vertical distribution for different chemical elements. Regarding both the radii and the effective temperatures, the derived values from our interferometric observations and from self-consistent modelling are consistent within better than $2\sigma$ for nine targets out of ten. We thus benchmark nine Ap stars for effective temperatures ranging from 7200 and 9100 K, and luminosities ranging between 7 $L_\odot$ and 86 $L_\odot$.

**Conclusions.** These results will be key for the future derivation of accurate radii and other fundamental parameters of fainter peculiar stars for which both the sensitivity and the angular resolution of the current interferometers are not sufficient. Within the context of the observations of Ap stars with the Transiting Exoplanet Survey Satellite (TESS), these interferometric measurements are crucial for testing the mechanism of pulsation excitation at work in these peculiar stars. In particular, our interferometric measurements provide accurate locations in the Hertzsprung-Russell diagram for hot Ap stars among which pulsations may be searched for with TESS, putting to test the blue edge of the theoretical instability strip. These accurate locations could be used to derive masses and ages of these stars through a specific grid of models, and to test correlations between the properties of these peculiar stars and their evolutionary state.

**Key words.** methods: observational – techniques: high angular resolution – techniques: interferometric – stars: fundamental parameters

## 1. Introduction

Magnetic fields are detected in a large variety of stars throughout the spectral type range and are known to play a crucial role during star formation and evolution. They have an impact on numerous physical phenomena at play within the stars and their environments such as accretion, mass loss, turbulence, and pulsations, and on some fundamental stellar properties such

as rotational speed and stellar chemical abundances (Donati & Landstreet 2009). One of the best examples of how stellar magnetic fields impact stellar surface properties is found in the chemically peculiar Ap stars. These are stars of spectral types A and B, characterized by strong and large-scale organized magnetic fields typically of several kG, and abundance inhomogeneities of a few orders of magnitude with respect to solar values, leading to spotted surfaces (Preston 1974). They also exhibit low rotational speeds (≤100 km s⁻<fnref n="1" />) supposed to be due to magnetic braking (Mathys 2017). Owing to their abnormal

<footnote n="1">¹ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France
e-mail: karine.perraut@univ-grenoble-alpes.fr</footnote>

<footnote n="2">² Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal</footnote>

<footnote n="3">³ School of Physics and Astronomy, University of Birmingham, Birmingham, B15 2TT, UK</footnote>

<footnote n="4">⁴ Institute of Astronomy, Russian Academy of Sciences, Pyatnitskaya 48, 119017 Moscow, Russia</footnote>

<footnote n="5">⁵ Max-Planck Institute fur Sonnensystemforschung, Justus-von-Liebig-Weg 3, 37077 Göttingen, Germany</footnote>

<footnote n="6">⁶ Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Laboratoire Lagrange, France</footnote>

<footnote n="7">⁷ Univ. Lyon, Univ. Lyon1, ENS de Lyon, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, 69230 Saint-Genis-Laval, France</footnote>

<footnote n="8">⁸ CHARA Array, Mount Wilson Observatory, 91023 Mount Wilson CA, USA</footnote>

<footnote n="9">⁹ Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, 3001 Leuven, Belgium</footnote>

<footnote>* CHARA/VEGA observations.</footnote>
<table><caption>**Table 1.** Stellar properties for the complete Ap sample, with $T_{\text{eff}}$ and $L$ derived from Kochukhov & Bagnulo (2006).</caption><thead><tr><th>#</th><th>HD</th><th><i>B</i></th><th><i>V</i></th><th>Sp. type</th><th><i>T</i><sub>eff</sub> [K]</th><th><i>L</i> [<i>L</i><sub>☉</sub>]</th></tr></thead><tbody><tr><td>1</td><td>4778</td><td>6.16</td><td>6.15</td><td>A3Vp</td><td>9980</td><td>32.4</td></tr><tr><td>2</td><td>24712<sup>(a)</sup></td><td>6.00</td><td>6.32</td><td>A9p</td><td>7190</td><td>7.4</td></tr><tr><td>3</td><td>108662</td><td>5.18</td><td>5.2</td><td>A0p</td><td>10 500</td><td>64.6</td></tr><tr><td>4</td><td>108945</td><td>5.47</td><td>5.4</td><td>A2pv</td><td>8950</td><td>51.3</td></tr><tr><td>5</td><td>118022<sup>(b)</sup></td><td>4.94</td><td>4.96</td><td>Ap</td><td>9060</td><td>28.8</td></tr><tr><td>6</td><td>120198</td><td>5.66</td><td>5.7</td><td>B9pe</td><td>10 540</td><td>44.7</td></tr><tr><td>7</td><td>128898<sup>(c)</sup></td><td>3.19</td><td>3.43</td><td>A7p</td><td>7670</td><td>11.0</td></tr><tr><td>8</td><td>137909<sup>(d)</sup></td><td>3.68</td><td>3.97</td><td>F2p</td><td>7430</td><td>27.5</td></tr><tr><td>9</td><td>153882</td><td>6.26</td><td>6.3</td><td>B9p</td><td>9730</td><td>100</td></tr><tr><td>10</td><td>176232<sup>(e)</sup></td><td>5.89</td><td>6.14</td><td>A7p</td><td>7930</td><td>20.9</td></tr><tr><td>11</td><td>188041<sup>(f)</sup></td><td>5.81</td><td>5.6</td><td>A5p</td><td>8430</td><td>35.5</td></tr><tr><td>12</td><td>201601<sup>(g)</sup></td><td>4.68</td><td>4.94</td><td>A9p</td><td>7620</td><td>12.6</td></tr><tr><td>13</td><td>204411<sup>(f)</sup></td><td>5.37</td><td>5.3</td><td>A4pv</td><td>8750</td><td>89.1</td></tr><tr><td>14</td><td>220825</td><td>4.94</td><td>4.9</td><td>A2pv</td><td>9080</td><td>22.4</td></tr></tbody></table>

**Notes.** The objects previously observed are in italics, and the corresponding references are given in the Notes.

**References.** (a)Perraut et al. (2016); (b)Perraut et al. (2015); (c)Bruntt et al. (2008); (d)Bruntt et al. (2010); (e)Perraut et al. (2013); (f)Romanovskaya et al. (2019b); (g)Perraut et al. (2011).

stellar surfaces, the determination of their fundamental parameters by (spectro)-photometry can be biased. Using long-baseline interferometry to provide accurate angular diameters of these particular stars appears to be a promising approach to overcoming the difficulties in deriving accurate global parameters. However, this approach is also very challenging because of the stars' small angular size (less than 1 mas).

The first interferometric observations of the largest Ap star $\alpha$ Cir by Bruntt et al. (2008) were performed with the VLTI (Glindemann et al. 2003). Then several campaigns on the brightest and the largest Ap stars were conducted with the CHARA array (Ten Brummelaar et al. 2010) since this interferometric array with its baselines $B$ as long as 330 m and its visible instrumentation allows targets as small as $\lambda/2B = 0.2$ mas to be partially resolved. To date, accurate angular diameters have been determined by optical long-baseline interferometry for eight Ap stars (given in italics in Table 1). In this paper we report on the angular diameter measurements of six additional Ap stars in the northern hemisphere, which extends the Ap sample for which we derive stellar radii by interferometry to 14 (hereafter referred to as the complete Ap sample). Our objective is to cover a wide range of luminosities and effective temperatures, and thus to probe the theoretical instability strip where the sub-category of Ap stars exhibiting pulsations, the rapidly oscillating Ap (roAp) stars, are expected (Balmforth et al. 2001; Cunha 2002; Cunha et al. 2019). The effective temperatures of our complete Ap sample range from about 7200 to 10 550 K and the luminosities from about 7 to 100 $L_{\odot}$ (Table 1). Ultra-precise photometric light curves will be available for most known Ap stars once the NASA Transiting Exoplanet Survey Satellite (TESS) mission (Ricker et al. 2014) completes its first two years of science observations. Having accurate stellar parameters and the necessary data to search for the presence of pulsations in the sample of stars considered in this work is an important step towards testing theoretical predictions for the driving of pulsations in roAp stars (Cunha et al. 2020).

Our aim is to derive stellar radii and effective temperatures that are as independent of stellar models as possible. Once

validated these models could be used to derive radii of fainter and/or farther stars for which the sensitivity and/or the angular resolution of interferometry are currently not sufficient. When looking at the Renson catalogue of Ap, HgMn, and Am stars (Renson & Manfroid 2009), most of the targets are located in the southern hemisphere and the histograms of the stars for both hemispheres is centred on Vmag = 9 (Fig. 1, left), which is out of reach of the sensitivity of the current visible interferometric instrumentations whose best performance corresponds to a limiting magnitude of about 8–8.5 in R band (Ireland et al. 2008; Ligi et al. 2013). In addition, when looking at the expected angular sizes of the targets that are brighter than this limit, we find that only a few of them are above 0.2 mas (Fig. 1, right), a value that is already very challenging to measure.

We present the interferometric observations and the angular diameters and radii derived from our last CHARA/VEGA observations in Sect. 2. We then use (spectro-)photometric data to compute the bolometric fluxes based solely on observational data in Sect. 3. In Sect. 4 we compute the stellar effective temperatures and luminosities and discuss our results in Sect. 5.

## 2. Last interferometric observations of Ap stars with VEGA/CHARA

### 2.1. Observations

During the 2016–2018 period, we observed six Ap stars with different telescope triplets of the CHARA array (Ten Brummelaar et al. 2010) and the VEGA visible instrument (Mourard et al. 2009). We operated VEGA in parallel with the K-band CLIMB beam combiner acting as a fringe tracker (Sturmann et al. 2010). We used the medium spectral resolution of VEGA ($R \sim 6000$) and recorded fringes on spectral bands centred around $\sim 700$ or $\sim 730$ nm, depending on the targets. We interleaved our target observations with calibrator observations, following the sequence calibrator-target-calibrator, with 40 or 60 blocks of 2500 short exposures (of 10 ms) per star. We swapped from one star to another every 30 min to correctly calibrate the instrumental transfer function. We used the SearchCal software (Bonneau et al. 2006) proposed by the JMMC<fnref n="1" /> to find suitable calibrators. The observation log is given in Appendix A along with the spatial frequency coverage for all targets.

### 2.2. Data processing

We used the standard VEGA reduction pipeline (Mourard et al. 2009) to compute the squared visibilities over spectral bands of 15 or 20 nm according to the width of water absorption lines in the spectrum. Such large bandwidths lead to a mean spectral resolution of about 40, and our visibility measurements are thus mainly sensitive to the stellar photosphere. We followed the same processing as described in our previous papers (Perraut et al. 2011, 2013, 2015, 2016) to compute raw squared visibility for each block of 2500 individual frames. Even though we recorded data with telescope triplets, we cannot retrieve closure phases on our targets since the signal-to-noise ratio on this observable is strongly limited by the saturation limit of the VEGA photon-counting detectors (see Fig. 4 in Mourard et al. 2012).

### 2.3. Angular diameters and radii

We followed the same approach as in our previous papers to calibrate the target squared visibilities $V^2$ by using the angular

<footnote n="1">¹ Available at www.jmmc.fr/searchcal</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Left: Vmag histograms of the Renson targets (Renson & Manfroid 2009) for the southern and northern hemispheres. Right: expected angular diameter as a function of the brightest Ap stars of the northern (orange circles) and southern (blue triangles) hemispheres. The yellow area corresponds to the sensitivity and the angular resolution of the visible instrumentations of the CHARA array.</figcaption></figure>

<table><caption>**Table 2.** Uniform-disc angular diameters in the R band of our calibrators from the JSDC catalogues.</caption><thead><tr><th>Calibrator</th><th>V-K</th><th>θ<sub>UD</sub> [mas]</th></tr></thead><tbody><tr><td>HD 4335</td><td>-0.08</td><td>0.184 ± 0.013</td></tr><tr><td>HD 108765</td><td>0.25</td><td>0.257 ± 0.019</td></tr><tr><td>HD 118214</td><td>-0.04</td><td>0.230 ± 0.016</td></tr><tr><td>HD 120874</td><td>0.26</td><td>0.194 ± 0.014</td></tr><tr><td>HD 151862</td><td>0.05</td><td>0.19 ± 0.015</td></tr><tr><td>HD 156653</td><td>0.076</td><td>0.18 ± 0.016</td></tr><tr><td>HD 217186</td><td>0.21</td><td>0.188 ± 0.013</td></tr><tr><td>HD 223855</td><td>0.05</td><td>0.170 ± 0.012</td></tr><tr><td>HD 224103</td><td>-0.12</td><td>0.163 ± 0.011</td></tr><tr><td>HD 224926</td><td>-0.31</td><td>0.222 ± 0.016</td></tr></tbody></table>

diameters of the calibrators provided by the JSDC catalogues (Bonneau et al. 2006; Bourges et al. 2017) to compute the transfer function. We adopted a conservative uncertainty of 7% on the uniform angular diameter (Table 2). Our calibrators are stars of A-B spectral types. For these early-type stars, since the JSDC diameter estimation exhibits discrepancies with the surface brightness–colour relations (SBCRs) for $V-K < 0$, we adopted the uniform diameters provided by the JSDC2 catalogue which shows a better consistency with SBCRs (see Fig. 3d in Challouf et al. 2014). The calibrated squared visibility curves versus spatial frequencies are given in Fig. 2.

We determined the uniform-disc angular diameter $\theta_{UD}$ of our targets in the R band using the model fitting tool LITpro<fnref n="2" /> (Tallon-Bosc et al. 2008), and the limb-darkened angular diameters $\theta_{LD}$ considering the linear limb-darkening coefficients in the R band from the Claret tables (Claret & Bloemen 2011) when using an effective temperature range in agreement with the values in Table 1 and a surface gravity log g between 4 and 4.5.

Stellar radii R (in solar radius, $R_{\odot}$) can be derived from our angular diameter measurements $\theta_{LD}$ (in mas), and the parallaxes $\pi_{P}$ (in seconds of arc) through the formula

$$R = \frac{\theta_{LD}}{9.305 \times \pi_{P}}. \qquad (1)$$

We used the parallaxes $\pi_{P}$ from the second Gaia DR2 release (Gaia Collaboration 2018) and computed the error bars on $\pi_{P}$ using the formula provided by Lindegren et al. (2018) at the IAU colloquium:

$$\sigma_{\text{ext}} = \sqrt{k^2 \sigma_i^2 + \sigma_s^2}. \qquad (2)$$

Here $\sigma_i$ is the internal error provided in the Gaia DR2 catalogue, $k = 1.08$, and $\sigma_s = 0.021$ mas is the systematic error for $G < 13$.

All the angular diameters, the parallaxes we consider, and the radii are given in Table 4.

## 3. Bolometric fluxes

We follow the same approach as in our previous papers to derive the bolometric fluxes from observational spectral energy distributions (SEDs; e.g. Perraut et al. 2011, 2013, 2015, 2016). So, in addition to the six Ap targets presented in the last section, we also computed the SED of HD 188041 and HD 204411, which was not done in Romanovskaya et al. (2019b).

### 3.1. Building the spectral energy distributions

To compute the bolometric fluxes from the SED, we gathered the following spectro-photometric data:

- For wavelengths in the range $1265 \AA < \lambda < 3180 \AA$, we used rebinned spectra from the Sky Survey Telescope obtained at the IUE Newly Extracted Spectra (INES) data archive<fnref n="3" />. We removed all bad pixels and measurements with negative flux using the quality flag (Garhart et al. 1997);

- We collected the photometric data from TD1 (Thompson et al. 1978) at 1565, 1965, 2365, and 2740 Å;

- For wavelengths in the visible range we used flux-calibrated ELODIE (Prugniel & Soubiran 2001; Prugniel et al. 2007) and/or MILES (Sánchez-Blázquez et al. 2006) spectra, along with spectroscopic observations with the Boller & Chivens long-slit spectrograph mounted at the 2.1 m telescope at the Observatorio Astronómico Nacional (OAN) at San Pedro Mártir in Mexico (see Romanovskaya et al. 2019b for the data description and processing).

<footnote n="2">² www.jmmc.fr/litpro_page.htm</footnote>

<footnote n="3">³ http://sdc.cab.inta-csic.es/cgi-ines/IUEdbsMY</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Squared VEGA visibilities (circles) and uniform-disc models for the best fits (orange solid lines) and for the angular diameters at ±1σ (grey dotted lines).</figcaption></figure>

– For wavelengths in the range 3300 Å < λ < 7850 Å we also collected the photometric data from the Adelman catalogue (Adelman et al. 1989) and converted the flux as described in our Appendix B. We used these data to calibrate the visible spectrum of HD 108662 in flux;

– At longer wavelengths (in the *J*, *H*, and *K* bands) we collected the photometric data from the 2MASS All Sky Catalog of point sources<fnref n="4" /> (Cutri et al. 2003) and corrected for the filter responses using the zeropoints from Cohen et al. (2003).

At the lower wavelength end of the spectral distribution, we
performed a linear interpolation on logarithmic scale between

912 and 1842 Å, assuming zero flux at 912 Å. At the higher wave-
length end of the spectral distribution, we performed a linear
interpolation on logarithmic scale using the 2MASS magnitudes
and assuming zero flux at 1.6 × 10⁶ Å.

All the datasets we considered are found in Table 3, and the
corresponding spectral distributions are given in Fig. 3. The vis-
ible spectra appear to have the same shape for all our targets, but
scaling in amplitude with the *V* magnitude; they all are in agree-
ment with the Adelman measurements when available. For the
lower wavelengths (between 1265 and 3180 Å) the few spectra we
collected exhibit complex features; when both datasets are avail-
able, the INES ones and the TD1 ones are in good agreement.
Finally, all the Adelman datasets display a flux that does not vary

<footnote n="4">⁴ http://irsa.ipac.caltech.edu/applications/Gator/</footnote>
<table><caption>**Table 3.** Spectroscopic and photometric datasets used for the bolometric flux computation with the corresponding spectral ranges.</caption><thead><tr><th>Target</th><th>INES</th><th>TD1</th><th>Adelman <sup>(+)</sup></th><th>Visible spectra</th><th>2MASS</th><th>f<sub>bol</sub> [erg cm<sup>-2</sup> s<sup>-1</sup>] This work</th><th>f<sub>bol</sub> [erg cm<sup>-2</sup> s<sup>-1</sup>] from modelling</th></tr></thead><tbody><tr><td>HD 4778</td><td>[1265 ... 3180]</td><td>Y</td><td>–</td><td>[3900 ... 6800] <sup>(a)</sup></td><td>J, H, K</td><td>0.97 ± 0.10 × 10<sup>-7</sup></td><td>–</td></tr><tr><td>HD 108662</td><td>[1265 ... 3180]</td><td>–</td><td>[3300 ... 7100]</td><td>[4000 ... 6800] <sup>(a)</sup></td><td>J, H, K</td><td>2.24 ± 0.21 × 10<sup>-7</sup> <sup>(•)</sup></td><td>2.53 ± 0.15 × 10<sup>-7</sup></td></tr><tr><td>HD 108945</td><td>[1265 ... 3180]</td><td>Y</td><td>[3300 ... 7100]</td><td>[3900 ... 6800] <sup>(a)</sup></td><td>J, H, K</td><td>1.68 ± 0.15 × 10<sup>-7</sup></td><td>–</td></tr><tr><td>HD 120198</td><td>[1265 ... 3180]</td><td>Y</td><td>[3300 ... 7100]</td><td>[3725 ... 7210] <sup>(b)</sup></td><td>J, H, K</td><td>1.62 ± 0.13 × 10<sup>-7</sup></td><td>–</td></tr><tr><td>HD 153882</td><td>–</td><td>Y</td><td>[3390 ... 7850]</td><td>[3535 ... 7405] <sup>(c)</sup></td><td>J, H, K</td><td>0.81 ± 0.06 × 10<sup>-7</sup></td><td>0.79 ± 0.04 × 10<sup>-7</sup></td></tr><tr><td>HD 188041</td><td>[1977 ... 3180]</td><td>Y</td><td>[3390 ... 7850]</td><td>[3535 ... 7405] <sup>(c)</sup></td><td>J, H, K</td><td>1.34 ± 0.16 × 10<sup>-7</sup></td><td>1.34 ± 0.25 × 10<sup>-7</sup></td></tr><tr><td>HD 204411</td><td>–</td><td>Y</td><td>[3300 ... 6800]</td><td>[3720 ... 7215] <sup>(b)</sup></td><td>J, H, K</td><td>1.90 ± 0.15 × 10<sup>-7</sup></td><td>1.84 ± 0.15 × 10<sup>-7</sup></td></tr><tr><td>HD 220825</td><td>–</td><td>–</td><td>[3300 ... 7850]</td><td>[3725 ... 7215] <sup>(b)</sup></td><td>J, H, K</td><td>2.30 ± 0.17 × 10<sup>-7</sup></td><td>–</td></tr></tbody></table>

**Notes.** The bolometric fluxes obtained from our work and from self-consistent atmosphere modelling are given in the last two columns. See text for details. <sup>(+)</sup>Adelman data are photometric values collected over the optical range. Visible spectra: <sup>(a)</sup>ELODIE; <sup>(b)</sup>OAN; <sup>(c)</sup>MILES. <sup>(•)</sup>We accounted for the 1.32 scale factor of the IUE flux, as explained in Sect. 3.1.

by more than 4% over the 3400–3700 Å wavelength range (hereafter called constant flux). This corresponds to the Balmer jump, which is very pronounced for the A-type stars. We noticed that the IUE flux for HD 108662 is too low compared to this Adelman constant flux (Fig. 3, right, first row), which might come from flux losses of the IUE satellite for bright sources. To match IUE flux with Adelman UV flux, the IUE flux needs to be scaled up by a factor of 1.32 (Romanovskaya et al. 2020b).

## 3.2. Estimating the bolometric fluxes

Three targets have a complete SED (HD 108662, HD 108945, and HD 120198) and their bolometric flux values can be determined by computing the area under the curve of this distribution. When spectra and photometric data are both available in a spectral range, we always used the spectra to compute the area under the SED.

For the objects for which the (spectro-)photometric data do not cover the full spectral range, we used the spectro-photometry dataset coming from the complete target sample to devise a way to estimate, for each star, the flux in wavelength ranges where data are missing. This was done by establishing the correlations between the different flux data, i.e. the TD1 flux integrated over different wavelength intervals, the Adelman constant flux, and the B and V magnitudes, as described in Appendix B. In practice we proceeded as follows:

- for HD 4778, for which no Adelman data are available, we used the linear relation we found between the constant flux and the B magnitude (Fig. B.1) to estimate a constant flux of $0.99 \times 10^{-11} \text{ erg cm}^{-2} \text{ s}^{-1} \text{ Å}^{-1}$ (white stars in Fig. 3, top left). Using the linear relation derived from the V magnitude led to the same constant flux within an accuracy better than 4%. We then computed the bolometric flux from flux integration.

- for HD 153882 and HD 204411 the INES spectra are missing, but TD1 data are available (blue triangles in Fig. 3, first column, third and fourth rows, respectively). For all the targets for which we have TD1 and INES spectra, we computed the integrated flux from 1975 to 3180 Å and showed that it strongly linearly correlates with the TD1 flux at 2365 Å (regression coefficient of 0.97; blue symbols and line on Fig. B.2). We used the derived linear law to determine an integrated flux of $7.9 \times 10^{-9} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 153882 and of $1.46 \times 10^{-8} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 204411 over the [1975 Å; 3180 Å] range. We also estimated this integrated flux by considering the correlation with the Adelman constant flux for the targets for which

we have both INES and Adelman datasets (regression coefficient of 0.7; orange symbols and line in Fig. B.2). With this approach, we estimated an integrated flux over the [1975 Å; 3180 Å] range of $9.7 \times 10^{-9} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 153882 and of $1.58 \times 10^{-8} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 204411. The agreement between the two approaches is better than 8% for HD 204411 and of the order of 20% for HD 153882. In the following we used the flux values obtained with the first method since the linear regression coefficient in that case is very close to 1.

- for HD 153882, HD 188041, and HD 204411, the INES spectra over the [1265 Å; 1975 Å] range are missing, but TD1 data are available. We proceeded as before. The flux integrated between 1265 and 1975 Å was correlated with the TD1 flux at 1965 Å (regression coefficient of 0.90; Fig. B.3). We then used this linear law to estimate an integrated flux over the [1265 Å; 1975 Å] range of $2.5 \times 10^{-9} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 153882, $5.1 \times 10^{-9} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 188041, and $7.1 \times 10^{-9} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 204411.

- for HD 220825 the situation is different since we do not have TD1 data and we could not apply the same methods as before. We estimated the integrated flux for the [1975 Å; 3180 Å] range by using its linear correlation with the Adelman constant flux (Fig. B.2). We obtained an integrated flux of $3.1 \times 10^{-8} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 220825. We studied for all our targets the correlation between the integrated fluxes over the [1265 Å; 1975 Å] range and over the [1975 Å; 3180 Å] range and derived a linear law (regression coefficient of 0.82; Fig. B.4). We used this linear law to estimate an integrated flux over the [1265 Å; 1975 Å] range of $1.3 \times 10^{-8} \text{ erg cm}^{-2} \text{ s}^{-1}$ for HD 220825.

As in our previous papers, we estimated the uncertainty associated with the bolometric flux by considering the following conservative uncertainties for each contribution: 3% uncertainty on the flux computed from the ELODIE, MILES, or OAN spectrum; 10% on the flux computed from the combined INES spectra; and 15% on the flux derived from interpolations and/or estimations from correlations. These largest uncertainties are on the fluxes over the lowest wavelength range, which generally correspond to about 12–15% of the total flux, and always to less than 25%.

## 3.3. Comparison with atmospheric models

The bolometric fluxes we obtain for our new targets are listed in Table 3. We derived uncertainties on the bolometric fluxes
<figure><img src="image_4.png" /></figure>

**Fig. 3.** Spectral energy distributions: UV-EUV INES (grey) and visible (orange) spectra, TD1 (blue triangles), Adelman (purple stars) and 2MASS (red circles) photometries. The dashed lines correspond to the interpolations. The insets are zoomed-in images of the UV-EUV and visible domains.

between 7 and 12%. We compared our values with the bolometric fluxes derived from a self-consistent spectroscopic analysis, when available. The latter approach is based on model fitting of spectro-photometric observations over a wide wavelength range, including the same INES, TD1, Adelman, and 2MASS datasets as ours. This method includes the modelling of the anomalous

chemical composition of the atmospheres of these peculiar stars and takes into account the inhomogeneous vertical distribution for different chemical elements. Model atmosphere calculations are conducted with the LLMODELS code (Shulyak et al. 2004). Through an iterative process, a theoretical SED is produced and compared to the observed one; the radius is a free parameter
<table><caption>**Table 4.** Limb-darkened angular diameters in the R band derived from interferometric measurements, Gaia DR2 parallaxes, radii, effective temperatures, and luminosities derived from interferometric measurements for the complete Ap sample.</caption><thead><tr><th rowspan="2">#</th><th rowspan="2">HD</th><th rowspan="2">θ<sub>LD</sub><br>[mas]</th><th rowspan="2">π<sub>P</sub><br>[mas]</th><th colspan="3">Interferometric observations</th><th colspan="2">Modelling</th><th rowspan="2">⟨B<sub>3</sub>⟩<br>[kG]</th></tr><tr><th>R<br>[R<sub>☉</sub>]</th><th>L<br>[L<sub>☉</sub>]</th><th>T<sub>eff</sub><br>[K]</th><th>R<br>[R<sub>☉</sub>]</th><th>T<sub>eff</sub><br>[K]</th></tr></thead><tbody><tr><td>1</td><td>4778</td><td>0.204 ± 0.008</td><td>9.32 ± 0.10</td><td>2.36 ± 0.12</td><td>34.9 ± 4.3</td><td>9135 ± 400</td><td></td><td></td><td>2.6 <sup>(a)</sup></td></tr><tr><td>2</td><td>24712</td><td>0.335 ± 0.009</td><td>20.62 ± 0.09</td><td>1.75 ± 0.05</td><td>7.6 ± 1.2</td><td>7235 ± 280</td><td>1.75 ± 0.04</td><td>7250 ± 100 <sup>(b)</sup></td><td>3.1 <sup>(c)</sup></td></tr><tr><td>3</td><td>108662</td><td>0.328 ± 0.009</td><td>13.54 ± 0.24</td><td>2.59 ± 0.12</td><td>38.1 ± 4.9</td><td>8880 ± 330</td><td>2.09 ± 0.10</td><td>10 225 ± 100 <sup>(d)</sup></td><td>3.3 <sup>(e)</sup></td></tr><tr><td>4</td><td>108945</td><td>0.315 ± 0.006</td><td>12.00 ± 0.15</td><td>2.82 ± 0.09</td><td>36.5 ± 4.2</td><td>8430 ± 270</td><td></td><td></td><td>0.1 <sup>(a)</sup></td></tr><tr><td>5</td><td>118022</td><td>0.346 ± 0.006</td><td>17.16 ± 0.25</td><td>2.17 ± 0.06</td><td>28.9 ± 3.0</td><td>9100 ± 190</td><td>2.16 ± 0.02</td><td>9140 ± 30 <sup>(f)</sup></td><td>3.0 <sup>(f)</sup></td></tr><tr><td>6</td><td>120198</td><td>0.226 ± 0.008</td><td>10.62 ± 0.078</td><td>2.28 ± 0.10</td><td>44.9 ± 4.3</td><td>9865 ± 370</td><td></td><td></td><td>1.6 <sup>(e)</sup></td></tr><tr><td>7</td><td>128898</td><td>1.105 ± 0.037</td><td>60.35 ± 0.14 <sup>(*)</sup></td><td>1.967 ± 0.066</td><td>10.51 ± 0.60</td><td>7420 ± 170</td><td>1.94 ± 0.01</td><td>7500 ± 130 <sup>(g)</sup></td><td>2.0 <sup>(h)</sup></td></tr><tr><td>8</td><td>137909</td><td>0.699 ± 0.017</td><td>28.60 ± 0.69 <sup>(*)</sup></td><td>2.63 ± 0.09</td><td>25.3 ± 2.9</td><td>7980 ± 180</td><td>2.47 ± 0.07</td><td>8100 ± 150 <sup>(i)</sup></td><td>5.4 <sup>(j)</sup></td></tr><tr><td>9</td><td>153882</td><td>0.193 ± 0.019</td><td>5.98 ± 0.053</td><td>3.46 ± 0.37</td><td>70.8 ± 6.5</td><td>8980 ± 600</td><td>3.35 ± 0.03</td><td>9150 ± 100 <sup>(k)</sup></td><td>3.8 <sup>(k)</sup></td></tr><tr><td>10</td><td>176232</td><td>0.275 ± 0.009</td><td>13.37 ± 0.08</td><td>2.21 ± 0.08</td><td>16.9 ± 1.4 <sup>(+)</sup></td><td>7900 ± 190 <sup>(+)</sup></td><td>2.35 ± 0.06</td><td>7550 ± 50 <sup>(l)</sup></td><td>1.5 <sup>(l)</sup></td></tr><tr><td>11</td><td>188041</td><td>0.247 ± 0.005</td><td>11.72 ± 0.11</td><td>2.26 ± 0.05</td><td>30.5 ± 4.2</td><td>9000 ± 360</td><td>2.39 ± 0.05</td><td>8770 ± 150 <sup>(m)</sup></td><td>3.6 <sup>(m)</sup></td></tr><tr><td>12</td><td>201601</td><td>0.564 ± 0.017</td><td>28.77 ± 0.25</td><td>2.11 ± 0.07</td><td>11.0 ± 0.93</td><td>7253 ± 235</td><td>1.98 ± 0.05</td><td>7550 ± 150 <sup>(j)</sup></td><td>4.0 <sup>(j)</sup></td></tr><tr><td>13</td><td>204411</td><td>0.328 ± 0.004</td><td>8.33 ± 0.12</td><td>4.23 ± 0.11</td><td>85.6 ± 9.2</td><td>8520 ± 220</td><td>4.42 ± 0.07</td><td>8300 ± 150 <sup>(m)</sup></td><td>&lt;0.8 <sup>(n)</sup></td></tr><tr><td>14</td><td>220825</td><td>0.339 ± 0.005</td><td>20.44 ± 0.23</td><td>1.78 ± 0.03</td><td>17.2 ± 1.7</td><td>8790 ± 230</td><td></td><td></td><td>2.0 <sup>(e)</sup></td></tr></tbody></table>

**Notes.** The last three columns give the radii and the effective temperatures as derived from self-consistent modelling and when considering the Gaia DR2 parallaxes, and the magnetic field strengths. Four targets have not been modelled yet. <sup>(*)</sup>For these two targets we keep the parallaxes from HIPPARCOS because of the binarity and/or the excessive brightness. <sup>(+)</sup>We report the average values of the effective temperatures and luminosities provided in Perraut et al. (2013).

**References.** <sup>(a)</sup>Glagolevskij (2019); <sup>(b)</sup>Shulyak et al. (2009); <sup>(c)</sup>Ryabchikova et al. (2007); <sup>(d)</sup>Romanovskaya et al. (2020b); <sup>(e)</sup>Kochukhov (2017); <sup>(f)</sup>Romanovskaya et al. (2019a); <sup>(g)</sup>Bruntt et al. (2008); <sup>(h)</sup>Kochukhov et al. (2009); <sup>(i)</sup>Bruntt et al. (2010); <sup>(j)</sup>Shulyak et al. (2013); <sup>(k)</sup>Romanovskaya et al. (2020a); <sup>(l)</sup>Nesvacil et al. (2013); <sup>(m)</sup>Romanovskaya et al. (2019b); <sup>(n)</sup>Ryabchikova et al. (2005).

together with the surface gravity and the effective temperature; the stellar parallax is used to scale the model flux to the distance of the star to be compared with the observed flux at the Earth's orbit; and the iterations continue until the fundamental parameters used in the SED fitting converge. Thus, the abundance pattern is consistent with the physical parameters of the final stellar atmosphere. This approach is described in detail in several papers (e.g. Shulyak et al. 2009; Romanovskaya et al. 2019b). As shown in Table 3, the two methods provide bolometric fluxes that are in excellent agreement. We use our values, derived solely from data, combined with our angular diameters to compute the fundamental parameters of our Ap targets.

# 4. Results on the complete sample of Ap stars

## 4.1. Fundamental parameters

We used our bolometric fluxes $f_{bol}$ (Table 3) and the parallaxes $\pi_P$ from Gaia DR2 (Gaia Collaboration 2018) to determine the luminosities of our targets as

$$L = 4\pi f_{bol} \frac{C^2}{\pi P^2}, \qquad (3)$$

where C stands for the conversion factor from parsecs to metres.

We used the measured angular diameters and bolometric fluxes to estimate the effective temperatures of our new targets with the equation

$$\sigma T_{\text{eff}}^4 = 4 f_{\text{bol}} / \theta_{\text{LD}}^2, \qquad (4)$$

where $\sigma$ stands for the Stefan-Boltzmann constant ($5.67 \times 10^{-5}$ erg cm$^{-2}$ s$^{-1}$ K$^{-4}$).

All our results are given in Table 4. We also report the fundamental parameters of the Ap stars that were observed previously via interferometry. We recomputed the values when considering the Gaia DR2 parallaxes, except for the two brightest targets HD 128898 and HD 137909, the latter being a binary star. For both we kept the HIPPARCOS parallaxes (van Leeuwen 2007) that are more accurate than the Gaia DR2 values.

## 4.2. Positions in the Hertzsprung-Russell (HR) diagram

We used our fundamental parameters to put the stars of the complete Ap sample in the HR diagram along with evolutionary tracks from Girardi et al. (2000), when considering a solar metallicity for all the stars (i.e. an initial metallicity $Z = 0.019$ and an initial helium abundance $Y = 0.273$), and models with overshooting (Fig. 4). We also overplotted the blue and red edges of the theoretical instability strip (blue dash-dotted lines), as computed by Cunha et al. (2019). Our sample nicely covers this instability strip. It can be clearly seen that none of the roAp stars in our sample (denoted by stars in Fig. 4) is located close to the blue edge of this theoretical instability strip since all of them have effective temperatures $\log T_{\text{eff}} < 3.90$. This is in line with the position of the overall sample of roAp stars known to date, as discussed by Cunha et al. (2019), and potentially indicates that the excitation mechanism at work in roAp stars is not fully understood.

# 5. Discussion and concluding remarks

Thanks to our interferometric measurements of angular diameters of Ap stars, we can compare the derived interferometric radii with those computed from self-consistent modelling, which are available for ten stars from our complete sample
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Complete sample of Ap stars in the Hertzsprung-Russell diagram. Stars denote the known rapidly oscillating Ap stars, while circles denote the Ap stars for which pulsations have not been detected yet. The dashed blue lines denote the blue and red edges of the theoretical instability strip. The thin dash-dotted lines correspond to evolutionary tracks for masses ranging from 1.4 $M_⊙$ to 2.5 $M_⊙$ and the thick black line denotes the ZAMS. The colour-coding gives the stellar magnetic field strength $B_s$. The numbers close to the symbols refer to Table 4.</figcaption></figure>

(Table 4). We recomputed these modelled radii using the Gaia DR2 parallaxes for all stars except HD 128898 and HD 137909, as explained in the notes of Table 4. Radii and effective temperatures derived from interferometric observations and self-consistent modelling are consistent within better than $2\sigma$ for nine of the ten targets (Fig. 5). The two radii determinations are in tension at a $4\sigma$ level for HD 108662, which is resolved by our observations and for which the scatter of the visibility points is the largest (Fig. 2, top right), i.e. the difference between the measured visibilities and the uniform-disc modelled visibilities can reach up to $\pm 2.6$ times the error on our visibility measurements. For this object such a scatter in visibilities remains when considering the 2017 dataset alone or the 2018 dataset alone. When looking at the weather conditions, for all observations, the Fried parameter is quoted to be larger than 9 cm, which corresponds to good conditions for such a source as bright as 5.2 in V and K bands. HD 108662 has the same calibrator as HD 108945, whose visibility curve does not exhibit such a large scatter. All these findings suggest that the observed scatter in the HD 108662 data might have a stellar origin. This star has a known companion at 140" as revealed in the Washington Double Star catalogue (WDS; Mason et al. 2001). This companion (HD 108651) is a spectroscopic binary. Even if hierarchical systems are commonly observed for these spectral types (HD 108662 is an A0 star), HD 108662A is not known as a spectroscopic binary (no note in WDS; no entry in Pourbaix et al. 2004). To go a step further, we used the CANDID tool (Gallenne et al. 2015) to estimate the $3\sigma$ detection limit of a companion within 50 mas around HD 108662A. We defined two sets of observations, one for the 2017 observations (recorded in April and May) and one for the observations recorded in April 2018. Despite the lack of closure phases, we can rule out a companion two magnitudes fainter than the star (i.e. $V \sim 7.2$). Based on previous studies, even if HD 108662 is quoted in Simbad as a variable star of $\alpha^2$CVn type and, as such, might exhibit stellar spots, it seems difficult

to explain such a scatter in visibility as being due to a stellar spot; Shulyak et al. (2014) present several simulations of spotted surfaces of chemically peculiar stars and mention that the impact of a spot becomes noticeable once the squared visibility is below 0.3, provided that the visibility accuracy is better than 1%. In our case, even if the squared visibilities are as small as 0.25 along the longest baselines mainly oriented along a position angle of 100–130° (see the spatial frequency coverage in Fig. A.1), the visibility differences we observed are much larger than 1%. The too-large radius of HD 108662 we derived from the interferometric angular diameter leads to a too-small effective temperature (see object 3 in Fig. 5) which does not allow the visible spectrum, and especially the H lines, to be correctly fitted (Romanovskaya et al. 2020b).

Benchmarking the Ap stars that have been already self-consistently modelled is of strong interest for deriving accurate radii, and then other fundamental parameters, of fainter Ap stars for which both the sensitivity and the angular resolution of the current interferometers are not sufficient. As clearly illustrated in this work with HD 153882, one of the faintest and smallest targets in our sample, deriving an accurate angular diameter for targets as small as 0.2 mas is very challenging. We have to cope with calibration issues, and thus with large error bars on the radius and the effective temperatures. Thus, for most of the Ap stars, indirect spectroscopic analysis would be the relevant approach to derive accurate fundamental parameters.

In order to further test the excitation mechanism at work in roAp stars, which is not yet fully understood, it is important to look for pulsations in hotter Ap stars whose location in the HR diagram is known accurately. That is the case of the hotter stars in our sample (indicated by circles in Fig. 4) whose effective temperatures range between 8000 and 10000 K (i.e. $3.90 < \log T_{\text{eff}} < 4.0$). Seven of these stars will be observed by the TESS satellite during the currently ongoing second year of science operations, and it will be possible to establish if any of these
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Comparison between the radii (**left**) and the effective temperatures (**right**) derived from interferometric measurements and those derived from the self-consistent spectroscopic modelling. The dash lines give the 1:1 relations. The numbers close to the symbols refer to Table 4.</figcaption></figure>

hotter stars exhibits pulsations with amplitudes within reach of the TESS observations. A detailed modelling of the stars, aimed at deriving their masses and ages based on a comprehensive grid of models, will be presented in a future work.

With a precise location of the Ap stars in the HR diagram, one can also consider searching for evidence of a correlation between the magnetic field strength and stellar age. This question is still a matter of debate. Looking at the distribution of the Ap magnetic stars in the HR diagram, Hubrig et al. (2000) claimed, for stars with masses smaller than $3 M_\odot$, that magnetic fields only appear in stars that have already completed at least approximately 30% of their main sequence lifetime. However, this is contradicted by Bagnulo et al. (2003) who have observed the Ap star HD 66318. They report the largest magnetic field detected in an Ap star, but this target has only completed about 10–20% of its main sequence life. This is also in contradiction with the findings of the study of Kochukhov & Bagnulo (2006) on the evolutionary state of these magnetic stars, and with the distribution of roAp stars of a TESS sample presented in Fig. 1 of Cunha et al. (2019), where about 20 targets are very close to the zero age main sequence (ZAMS). Our accurate determination of the fundamental parameters of HD 220825 (object 14) is in direct line with such a distribution and clearly sets this $2 M_\odot$ magnetic Ap star very close to the ZAMS. More recently, based on a study of stars in open clusters, Landstreet (2011) argued that the typical magnetic field declines rapidly as the Ap stars evolve on the main sequence. With this in mind, we collected the surface magnetic field strength of our targets in the literature and reported them in Table 4. Looking at the position on the main sequence of each target with respect to the ZAMS represented by the thick black line in Fig. 4, no apparent correlation emerges in our results; for instance, the youngest star of our sample, HD 220825 (object 14), exhibits a stellar magnetic field smaller than 3 kg, while HD 137909 (object 8) located at the middle of the main sequence has a magnetic field more than two times higher. More generally, a decrease in the magnetic field strength with age would translate into green symbols close to and along the ZAMS, and black symbols farther away at the end of the main sequence. Perhaps we could not see a clear correlation because the number of stars in our study is limited and

the star-to-star variations do not allow us to identify a possible statistical correlation between the two properties. It could also be that the correlation found by Landstreet (2011) is not present for stars with masses lower than $2.5 M_\odot$. The stars in their study were typically of higher mass than the stars in our study, with their lowest mass bin including stars in the range of $2-3 M_\odot$. In the future, the use of self-consistent atmospheric modelling will allow us to increase the number of stars in the mass range studied here, and thus to verify whether the evolutionary dependence found by the authors is also present for the lower mass Ap stars.

**Acknowledgements.** The authors thank the anonymous referee for her/his fruitful comments. This work is based upon observations obtained with the Georgia State University Center for High Angular Resolution Astronomy Array at Mount Wilson Observatory. The CHARA Array is supported by the National Science Foundation under Grant No. AST-1636624 and AST-1715788. Institutional support has been provided from the GSU College of Arts and Sciences and the GSU Office of the Vice President for Research and Economic Development. 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 SIMBAD and VIZIER3 databases at CDS, Strasbourg (France) and of electronic bibliography maintained by the NASA/ADS system. This work has made use of the Jean-Marie Mariotti Center Aspro, SearchCal, and LITPRO services (available at http://www.jmmc.fr/). His work was supported by the "Programme National de Physique Stellaire" (PNPS) of CNRS/INSU co-funded by CEA and CNES. This work was supported by the "Action Spécifique pour la Haute Résolution Angulaire" (ASHRA) of CNRS/INSU co-funded by CNES. K.P. acknowledges funding from LABEX OSUG@2020. M.S.C. is supported by national funds through FCT – Fundação para a Ciência e a Tecnologia – in the form of a work contract and through the research grants UIDB/04434/2020, UIDP/04434/2020 and PTDC/FIS-AST/30389/2017, and by FEDER – Fundo Europeu de Desenvolvimento Regional through COMPETE2020 – Programa Operacional Competitividade e Internacionalização (grant: POCI-01-0145-FEDER-030389). A.R. acknowledges funding from RFBR, project number 19-32-90147.

## References

* Adelman, S. J. 1973, ApJ, 183, 95

* Adelman, S. J., Pyper, D. M., Shore, S. N., White, R. E., & Warren, W. H., J. 1989, A&AS, 81, 221

* Bagnulo, S., Landstreet, J. D., Lo Curto, G., Szeifert, T., & Wade, G. A. 2003, A&A, 403, 645
Balmforth, N. J., Cunha, M. S., Dolez, N., Gough, D. O., & Vauclair, S. 2001, MNRAS, 323, 362

Bonneau, D., Clausse, J.-M., Delfosse, X., et al. 2006, A&A, 456, 789

Bourges, L., Mella, G., Lafrasse, S., et al. 2017, VizieR Online Data Catalog: II/346

Bruntt, H., North, J. R., Cunha, M., et al. 2008, MNRAS, 386, 2039

Bruntt, H., Kervella, P., Mérand, A., et al. 2010, A&A, 512, A55

Challouf, M., Nardetto, N., Mourard, D., et al. 2014, A&A, 570, A104

Claret, A., & Bloemen, S. 2011, VizieR Online Data Catalog: J/A+A/352/075

Cohen, M., Wheaton, W. A., & Megeath, S. T. 2003, AJ, 126, 1090

Cunha, M. S. 2002, MNRAS, 333, 47

Cunha, M. S., Antoci, V., Holdsworth, D. L., et al. 2019, MNRAS, 487, 3523

Cunha, M., Alentiev, D., Brandao, I., & Perraut, K. 2020, MNRAS, in press

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

Donati, J. F., & Landstreet, J. D. 2009, ARA&A, 47, 333

Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1

Gallenne, A., Mérand, A., Kervella, P., et al. 2015, A&A, 579, A68

Garhart, M. P., Smith, M. A., Turnrose, B. E., Levay, K. L., & Thompson, R. W. 1997, IUE NASA Newsletter, 57, 1

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371

Glagolevskij, Y. V. 2019, Astrophys. Bull., 74, 66

Glindemann, A., Algomedo, J., Amestica, R., et al. 2003, SPIE Conf. Ser., 4838, 89–100

Hubrig, S., North, P., & Mathys, G. 2000, ApJ, 539, 352

Ireland, M. J., Mérand, A., ten Brummelaar, T. A., et al. 2008, SPIE Conf. Ser., 7013

Kochukhov, O. 2017, A&A, 597, A58

Kochukhov, O., & Bagnulo, S. 2006, A&A, 450, 763

Kochukhov, O., Shulyak, D., & Ryabchikova, T. 2009, A&A, 499, 851

Landstreet, J. D. 2011, ASP Conf. Ser., 449, 249

Ligi, R., Mourard, D., Nardetto, N., & Clausse, J. M. 2013, J. Astron. Instrum., 2, 1340003

Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, https://www.cosmos.esa.int/web/Gaia/dr2-known-issues#AstrometryConsiderations

Mason, B. D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. G., & Worley, C. E. 2001, AJ, 122, 3466

Mathys, G. 2017, A&A, 601, A14

Mourard, D., Clausse, J. M., Marcotto, A., et al. 2009, A&A, 508, 1073

Mourard, D., Harmanec, P., Stencel, R., et al. 2012, A&A, 544, A91

Nesvacil, N., Shulyak, D., Ryabchikova, T. A., et al. 2013, A&A, 552, A28

Perraut, K., Brandão, I., Mourard, D., et al. 2011, A&A, 526, A89

Perraut, K., Borgniet, S., Cunha, M., et al. 2013, A&A, 559, A21

Perraut, K., Cunha, M., Brandão, I., et al. 2015, A&A, 579, A85

Perraut, K., Brandão, I., Cunha, M., et al. 2016, A&A, 590, A117

Pourbaix, D., Tokovinin, A. A., Batten, A. H., et al. 2004, A&A, 424, 727

Preston, G. W. 1974, ARA&A, 12, 257

Prugniel, P., & Soubiran, C. 2001, A&A, 369, 1048

Prugniel, P., Soubiran, C., Koleva, M., & Le Borgne, D. 2007, VizieR Online Data Catalog: III/251

Renson, P., & Manfroid, J. 2009, A&A, 498, 961

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, Proc. SPIE, 9143, 914320

Romanovskaya, A., Ryabchikova, T., & Shulyak, D. 2019a, ASP Conf. Ser., 518, 173

Romanovskaya, A., Ryabchikova, T., Shulyak, D., et al. 2019b, MNRAS, 488, 2343

Romanovskaya, A. M., Ryabchikova, T. A., & Shulyak, D. V. 2020a, Astron. Lett., 46, 331

Romanovskaya, A. M., Ryabchikova, T. A., & Shulyak, D. V. 2020b, INASAN Sci. Rep., 5, 219

Ryabchikova, T., Leone, F., & Kochukhov, O. 2005, A&A, 438, 973

Ryabchikova, T., Sachkov, M., Weiss, W. W., et al. 2007, A&A, 462, 1103

Sánchez-Blázquez, P., Peletier, R. F., Jiménez-Vicente, J., et al. 2006, MNRAS, 371, 703

Shulyak, D., Tsymbal, V., Ryabchikova, T., Stütz, C., & Weiss, W. W. 2004, A&A, 428, 993

Shulyak, D., Ryabchikova, T., Mashonkina, L., & Kochukhov, O. 2009, A&A, 499, 879

Shulyak, D., Ryabchikova, T., & Kochukhov, O. 2013, A&A, 551, A14

Shulyak, D., Paladini, C., Li Causi, G., Perraut, K., & Kochukhov, O. 2014, MNRAS, 443, 1629

Sturmann, J., Ten Brummelaar, T., Sturmann, L., & McAlister, H. A. 2010, SPIE Conf. Ser., 7734, 773403

Tallon-Bosc, I., Tallon, M., Thiébaut, E., et al. 2008, SPIE Conf. Ser., 7013, 70134H

Ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2010, SPIE Conf. Ser., 7734

Thompson, G. I., Nandy, K., Jamar, C., et al. 1978, VizieR On-line Data Catalog: II/59B

van Leeuwen, F. 2007, A&A, 474, 653
## Appendix A: Log of VEGA observations

The log of the VEGA observations is given in Table A.1 and the corresponding spatial frequency coverage (UV-planes) in Fig. A.1.

<table><caption>**Table A.1.** Observing log.</caption><thead><tr><th>Target</th><th>Date</th><th>UT time</th><th>Telescopes</th><th>Calibrators</th></tr></thead><tbody><tr><td rowspan="8">HD 4778</td><td>2017-07-26</td><td>12:07</td><td>E1W1W2</td><td>HD 4335</td></tr><tr><td>2017-09-16</td><td>10:34</td><td>S1E1</td><td>HD 4335</td></tr><tr><td>2017-09-16</td><td>11:09</td><td>S1E1</td><td>HD 4335</td></tr><tr><td>2017-09-16</td><td>11:55</td><td>S1E1</td><td>HD 4335</td></tr><tr><td>2017-09-16</td><td>12:34</td><td>S1E1</td><td>HD 4335</td></tr><tr><td>2017-10-16</td><td>6:07</td><td>E1E2W1</td><td>HD 4335</td></tr><tr><td>2017-10-16</td><td>7:38</td><td>E1E2W1</td><td>HD 4335</td></tr><tr><td rowspan="6">HD 108662</td><td>2017-04-12</td><td>8:24</td><td>E1E2W1</td><td>HD 108765</td></tr><tr><td>2017-05-12</td><td>4:30</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2017-05-12</td><td>5:27</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-25</td><td>4:13</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-25</td><td>5:51</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-25</td><td>6:19</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td rowspan="6">HD 108945</td><td>2018-04-24</td><td>4:28</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-24</td><td>5:05</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-24</td><td>5:59</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-27</td><td>3:49</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-27</td><td>4:14</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td>2018-04-27</td><td>5:06</td><td>E1W1W2</td><td>HD 108765</td></tr><tr><td rowspan="6">HD 120198</td><td>2015-06-02</td><td>5:53</td><td>E1E2W1</td><td>HD 120874, HD 118214</td></tr><tr><td>2015-06-02</td><td>6:39</td><td>E1E2W1</td><td>HD 120874, HD 118214</td></tr><tr><td>2016-06-13</td><td>6:03</td><td>E1W1W2</td><td>HD 118214</td></tr><tr><td>2016-07-26</td><td>4:25</td><td>E1E2W1</td><td>HD 118214</td></tr><tr><td>2016-07-27</td><td>4:24</td><td>E1E2W1</td><td>HD 118214</td></tr><tr><td>2016-07-27</td><td>5:01</td><td>E1E2W1</td><td>HD 118214</td></tr><tr><td rowspan="4">HD 153882</td><td>2015-09-01</td><td>4:31</td><td>E2W1</td><td>HD 151862, HD 156653</td></tr><tr><td>2015-09-01</td><td>5:23</td><td>E2W1</td><td>HD 151862, HD 156653</td></tr><tr><td>2016-07-28</td><td>5:33</td><td>E1E2W1</td><td>HD 151862, HD 156653</td></tr><tr><td>2017-07-26</td><td>4:52</td><td>E1W1W2</td><td>HD 156653</td></tr><tr><td rowspan="10">HD 220825</td><td>2016-08-20</td><td>10:50</td><td>E1E2W1</td><td>HD 224103</td></tr><tr><td>2016-08-25</td><td>10:50</td><td>E1E2W1</td><td>HD 223855, HD 224103</td></tr><tr><td>2016-09-30</td><td>4:03</td><td>E1W1W2</td><td>HD 217186, HD 224926</td></tr><tr><td>2016-09-30</td><td>5:05</td><td>E1W1W2</td><td>HD 217186, HD 224926</td></tr><tr><td>2016-09-30</td><td>5:29</td><td>E1W1W2</td><td>HD 217186, HD 224926</td></tr><tr><td>2016-09-30</td><td>6:00</td><td>E1W1W2</td><td>HD 217186, HD 224926</td></tr><tr><td>2017-07-26</td><td>9:21</td><td>E1W1W2</td><td>HD 217186, HD 224926</td></tr><tr><td>2017-07-26</td><td>10:13</td><td>E1W1W2</td><td>HD 217186, HD 224926</td></tr><tr><td>2017-07-30</td><td>10:02</td><td>E1E2W1</td><td>HD 224926</td></tr><tr><td>2017-07-30</td><td>10:52</td><td>E1E2W1</td><td>HD 224926</td></tr></tbody></table>

**Notes.** The data are available via the Jean-Marie Mariotti Center OiDB service (available at http://oidb.jmmc.fr).
<figure><img src="image_7.png" /><figcaption>Fig. A.1. Maps of the baselines in the UV-plane of the CHARA/VEGA observations.</figcaption></figure>

## Appendix B: Flux estimations

Adelman constant flux. For the Adelman data (Adelman 1973), we convert into erg cm⁻² s⁻¹ A⁻¹ the values f_λ provided in the catalogue with respect to the flux at 5556 Å (f_5556) and per unit of frequency as follows:

$$ \text{Flux}(\lambda) = 3.46 \times 10^{-9} \times \left(\frac{5556}{\lambda}\right)^2 \times 10^{-0.4(V_{\text{mag}}-0.026+f_{\lambda}-f_{5556})}. \quad (B.1) $$

When looking at the Adelman datasets displayed in Fig. 3 (purple stars), we observed that for all the stars, the flux over the

3400–3700 Å wavelength range is stable within 4% (referred to as constant flux in this paper). This constant flux corresponds to the Balmer jump and hence depends on the temperature and the density in the region responsible for the absorption. It is correlated with the B and V magnitudes, while it can change by a factor of more than two between stars of similar effective temperatures, which can explain the scatter in Fig. B.1.

Integrated flux in the bluest part of the spectra. We gathered the (spectro-)photometric data for the complete Ap sample and looked for correlations between integrated flux and photometric data (Figs. B.2–B.4).
<figure><img src="image_8.png" /><figcaption>**Fig. B.1.** Adelman constant flux as a function of Bmag (blue circles) and of Vmag (orange diamonds). The dash-dotted lines display the best linear laws for each dataset with the corresponding equations and regression coefficients at the top.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. B.3.** Integrated flux over [1265 Å; 1975 Å] based on INES spectra as a function of TD1 flux at 1965 Å (blue pentagons). The dash-dotted line displays the best linear law for this dataset with the corresponding equation and regression coefficient at the top.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. B.2.** Integrated flux over [1975 Å; 3180 Å] based on INES spectra as a function of TD1 flux at 2365 Å (blue stars) and of Adelman constant flux (orange triangles). The dash-dotted lines display the best linear laws for each dataset with the corresponding equations and regression coefficients at the top.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. B.4.** Integrated flux over [1265 Å; 1975 Å] based on INES spectra as a function of the integrated flux over [1975 Å; 3180 Å] (blue squares). The dash-dotted line displays the best linear law for this dataset with the corresponding equation and regression coefficient at the top.</figcaption></figure>