Astronomy & Astrophysics manuscript no. Vioque_et_al_2018 c ESO 2018 August 3, 2018

Gaia DR2 study of Herbig Ae/Be stars M. Vioque1, 2?, R. D. Oudmaijer1, D. Baines3, I. Mendigutía4, and R. Pérez-Martínez2

1 School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK. 2 Ingeniería de Sistemas para la Defensa de España (Isdefe), XMM/Newton Science Operations Centre, ESA-ESAC Campus, PO Box 78, 28691 Villanueva de la Cañada, Madrid, Spain. 3 Quasar Science Resources for ESA-ESAC, ESAC Science Data Center, PO Box 78, 28691 Villanueva de la Cañada, Madrid, Spain. 4 Centro de Astrobiología (CSIC-INTA), Departamento de Astrofísica, ESA-ESAC Campus, PO Box 78, 28691 Villanueva de la Cañada, Madrid, Spain.

Accepted for publication in Astronomy & Astrophysics.

ABSTRACT

Aims. We use Gaia Data Release 2 (DR2) to place 252 Herbig Ae/Be stars in the HR diagram and investigate their characteristics and properties. Methods. For all known Herbig Ae/Be stars with parallaxes in Gaia DR2, we collected their atmospheric parameters and photometric and extinction values from the literature. To these data we added near- and mid-infrared photometry, collected Hα emission line properties such as equivalent widths and line profiles, and their binarity status. In addition, we developed a photometric variability indicator from Gaia’s DR2 information. Results. We provide masses, ages, luminosities, distances, photometric variabilities and infrared excesses homogeneously derived for the most complete sample of Herbig Ae/Be stars to date. We find that high mass stars have a much smaller infrared excess and have much lower optical variabilities compared to lower mass stars, with the break at around 7M .Hα emission is generally correlated with infrared excess, with the correlation being stronger for infrared emission at wavelengths tracing the hot dust closest to the star. The variability indicator as developed by us shows that ∼25% of all Herbig Ae/Be stars are strongly variable. We observe that the strongly variable objects display doubly peaked Hα line profiles, indicating an edge-on disk. Conclusions. The fraction of strongly variable Herbig Ae stars is close to that found for A-type UX Ori stars. It had been suggested that this variability is in most cases due to asymmetric dusty disk structures seen edge-on. The observation here is in strong support of this hypothesis. Finally, the difference in dust properties occurs at 7M , while various properties traced at UV/optical wavelengths differ at a lower mass, 3M . The latter has been linked to different accretion mechanisms at work whereas the differing infrared properties and photometric variabilities are related to different or differently acting (dust-)disk dispersal mechanisms. Key words. stars: Herbig Ae/Be – stars: Hertzsprung-Russell diagram – stars: formation – stars: pre-main sequence – stars: emission- line – infrared: stars

1. Introduction nier et al. 2005) that may suggest a different accretion scenario with the disk reaching directly into the star (Kraus 2015). Other Herbig Ae/Be stars (HAeBes) are Pre-Main Sequence stars spectro-photometric (Mendigutía et al. 2011b; Cauley & Johns- (PMS) of intermediate mass, spanning the range between low Krull 2015 and Patel et al. 2017) and spectro-polarimetric stud- mass T-Tauri stars and the embedded Massive Young Stellar Ob- ies (Vink et al. 2002) also point in the direction that the accretion jects (MYSOs). They are optically bright so they are much eas- physics change within the Herbig Ae/Be stars mass range. In ad- ier to observe and to study than MYSOs and it is expected that dition, Herbig Be stars are more likely to be found in binaries within the mass range of HAeBes a change in accretion mech- than Herbig Ae stars (Baines et al. 2006). anism from the magnetically controlled accretion acting for T- Tauri stars (see Bouvier et al. 2007) to an, as yet, unknown An important indicator of their PMS nature, together with mechanism for high mass stars occurs. Indeed, there is evidence emission lines, is the infrared (IR) excess that also traces the that the magnetically driven accretion model is valid for Her- Herbig Ae/Be forming environment. The IR excess profile have been classified into two groups differentiated by a flat or rising arXiv:1808.00476v1 [astro-ph.SR] 1 Aug 2018 big Ae stars but not for several Herbig Be stars (Fairlamb et al. 2015; Ababakr et al. 2017; Oudmaijer 2017; Grady et al. 2010; shape of the continuum (Meeus et al. 2001). This difference has Schöller et al. 2016). Moreover, there are multiple evidences that a geometric origin depending on the presence of flaring outer Herbig Ae and T-Tauri stars behave more similarly than Herbig disks and puffed-up inner disks (Dullemond & Dominik 2004a, Be stars, and Herbig Ae and Herbig Be stars have different obser- 2004b, 2005), and the presence of gaps in the disk (Maaskant vational properties. Examples of this are the different outer gas et al. 2013; Honda et al. 2015). The IR excess of HAeBes is dispersal rates (higher for Herbig Be stars, Fuente et al. 1998), expected to be characteristic and different from the IR excess of the higher incidence of clustering scenarios for Herbig Be stars other similar objects like for example, ordinary Be stars (Finken- (Testi et al. 1999), and the evidences of Herbig Be stars hosting zeller & Mundt 1984). denser and larger inner gaseous disks (Ilee et al. 2014; Mon- Herbig Ae/Be stars are known to present irregular photomet- ric variations, with a typical timescale from days to weeks (Eiroa ? e-mail: [email protected] et al. 2002; Oudmaijer et al. 2001) and of the order of one mag-

Article number, page 1 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 nitude in the optical. This variability is typically understood as largely understudied. We therefore keep objects with F-type clas- due to variable extinction, due to for example rotating circum- sification in Chen et al. (2016) in the sample. Similarly, no upper stellar disks, or as an effect of rotation on cold photospheric spots limit in mass was imposed, leaving the separation between Mas- and also pulsation due to the source crossing the instability strip sive Young Stellar Objects (MYSOs) and HAeBes to the optical in the HR diagram (Marconi & Palla 1998). An extreme case brightness of the sources1. of large non-periodic photometric and polarimetric variations is Then, we crossmatched the sources with Gaia DR2. Detec- observed in UX Ori type stars (UXORs) with amplitudes up to tions were considered to be matched with the catalogue when 2−3 mag. Many of them are catalogued as HAeBes and their ex- their coordinates agreed to within 0.5 arcsecond. If more than treme variability is explained by eclipsing dust clouds in nearly one match was found we took the closest one. If no match was edge-on sources and the scattering radiation in the circumstellar found within 0.5 arcsecond, successive crossmatches with larger environment (see Grinin 2000 and references therein; Natta et al. apertures were performed up to 2 arcsecond. In these latter cases 1997 and Natta & Whitney 2000). an individual inspection of the crossmatch was applied. Finally, Infrared photometric variability, related to disk structure a comparison between the Johnson V band magnitudes and the variations, is not always correlated with the optical variability Gaia filters was done for each source in order to discard possi- (Eiroa et al. 2002) which implies that different mechanisms re- ble incorrect matches. This provides us with parallaxes for 254 garding both the disk structure and accretion underlie the final HAeBes. observed variability. Spectroscopic variability is also present in As Lindegren et al. (2018) point out, not all Gaia DR2 par- Herbig Ae/Be stars (Mendigutía et al. 2011a). allaxes are of the same quality, and some values - despite their With the advent of the second data release of Gaia (DR2, sometimes very small error bars - appear erroneous (e.g. Linde- Gaia Collaboration et al. 2016, Gaia Collaboration et al. 2018b), gren et al. 2018). We included the following constraint in astro- providing parallaxes to over 1.3 billion objects (Lindegren et al. metric quality following the indications in Appendix C of Linde- 2018), including the majority of known Herbig Ae/Be stars, the gren et al. (2018) and what was applied in Gaia Collaboration et time is right for a new study on the properties of the class. Gaia al. (2018). This constraint will remove from the sample objects DR2 contains a five dimensional astrometric solution (α, δ, µα, with spurious parallaxes: µδ and parallax ($)) up to G . 21 (white G band, described in Evans et al. 2018). Almost all of the known Herbig Ae/Be u < 1.2 × max(1, e−0.2(G−19.5)) (1) stars have parallaxes in Gaia DR2, which allowed luminosities to be derived and 252 HAeBes to be placed in the HR diagram, where G is the Gaia G band and u is the unit weight error, a tenfold increase on earlier studies using Hipparcos data alone. defined as the square root of the ratio of the astrometric_chi2_al The paper is organized as follows: In Sect. 2, we describe the and (astrometric_n_good_obs_al − 5) columns (Lindegren et al. data acquisition of not only the parallaxes, but also optical and 2018, their Equation C.2). 228 of our sources satisfy this condi- infrared photometry, effective temperatures, extinction values, tion. Hα emission line information and binarity. In Sect. 3 we derive Some objects are found to be very close to this condition, the stellar luminosities and place the objects in an Hertzsprung- PDS 144S, PV Cep and V892 Tau, and as we will show later, Russell (HR) diagram, while we also present a method to de- they would appear significantly below the Main Sequence in the rive a statistical assessment of the objects’ variability in Gaia’s HR-diagram. Given that the Lindegren condition is presented as database. In addition, we homogeneously derive masses and ages a guideline rather than a rule by the Gaia astrometry team, we for all the sources, together with near- and mid-infrared ex- decided to treat these three objects as if they satisfy Eq. 1 as cesses. In Sect. 4 we carry out an analysis of the data and present well. various correlations and interdependencies, which we discuss in We will refer to the set of astrometrically well behaved the context of intermediate mass star formation in Sect. 5. We sources as the high quality sample and to those that do not sat- conclude in section Sect. 6. isfy Eq. 1 as the low quality sample. We will not be able to place 2 sources in the HR diagram because lacking of appropriate pa- rameters (Sect. 2.2). In addition, we will move 5 more sources to 2. Data acquisition the low quality sample in Sect. 3.1 because of different reasons. Summarizing, there are 218 objects (228 − 3 − 2 − 5) in the final 2.1. Construction of the sample high quality sample and 34 in the low quality one. Information We have gathered the majority of Herbig Ae/Be stars known and about the objects in different samples is presented in separated proposed to date from different works (272, see Chen et al. 2016 tables at the end of the paper. The high quality sample will be the for a compilation of most of them). Chen et al. (2016) based one taken into account in further considerations unless otherwise their sample mostly on the work of Zhang et al. (2006) which in specified. turn is based on the work of The et al. (1994) and Vieira et al. Distances are not obtained by straightforwardly inverting the (2003). In addition, we included a few HAeBes from Alecian et parallax. The conversion of one parameter to the other one is not al. (2013); Baines et al. (2006); Carmona et al. 2010; Fairlamb et strictly trivial because of the non-linearity of the inverse func- al. 2015; Hernández et al. (2005); Manoj et al. (2006) and Sartori tion (see for example Bailer-Jones 2015). In the case of Gaia et al. (2010) that are not present in the aforementioned papers. DR2, Bailer-Jones et al. (2018) proposed distance values using Although Herbig Ae/Be stars have long been considered, a weak distance prior that follows a galactic model. Their dis- by definition, of type A or B, there should be some flexibility tances begin to differ from the distances obtained through sim- in this constraint as the physical boundary between Herbig Ae ple inversion for sources with large errors, σ$/$ & 0.5. Hence, stars and intermediate mass T-Tauris is fairly unstudied. This is 1 The MYSOs are typically infrared-bright and optically faint (Lums- because spectral types of T-Tauri stars are typically K-M with den et al. 2013). However, a number of optically visible objects are some G-type objects while Herbig Ae/Be stars are, quite un- known to have passed all selection criteria such as the early type ob- surprisingly, limited to A and B spectral type. Hence, pre-Main jects PDS 27 and PDS 37 that are also classified as Herbig Be stars Sequence stars of intermediate spectral types have often been (Ababakr et al. 2015).

Article number, page 2 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars in our initial (high and low quality) sample only a small subset number of Herbig Ae/Be stars that can be placed in the HR di- of 12 Herbig Ae/Be stars will suffer substantially from this ef- agram and for which we can derive stellar luminosities, masses, fect. Following the indications in Luri et al. (2018) on how to ages, IR excesses and variabilities in Sect. 3 is therefore reduced treat the Gaia parallaxes we decided to apply a simpler exponen- to 252 objects. tially decreasing prior to estimate distances. For completeness, we should note that the parallaxes provided by Gaia DR2 have a regional, not gaussian systematic error as large as 0.1mas and a 2.3. Infrared photometry global zero point error of about −0.029mas that are not included All the sources were crossmatched with the Two-Micron All in the gaussian random errors provided in the Gaia archive (see Sky Survey (2MASS, see Skrutskie et al. 2006) and with the Arenou et al. 2018 and Lindegren et al. 2018). Hence, the un- Wide-Field Infrared Survey Explorer ‘AllWISE’ all-sky cata- certainty in the parallaxes is slightly underestimated. The final logue (hereafter WISE, see Cutri et al. 2013). Both these sur- errors in the high quality sample range from 0.016 to 0.37mas. veys contain hundreds of millions of stars, guaranteeing a large Herbig Ae/Be stars have been historically confused with overlap with Gaia. We used a 3 arcsecond aperture for the cross- classical Be stars, with which they share many characteristics match. The few sources that did not lie within that 3 arcsecond (Rivinius et al. 2013; Klement et al. 2017; Grundstrom & Gies threshold were studied individually and, if present, their IR pho- 2006). Indeed, the nature of some of the objects in our sample tometry was included. This provides values and uncertainties for is still under debate. An interesting example in this respect is the J, H and Ks bands (1.24, 1.66 and 2.16µm respectively) and HD 76534, a B2Ve object that appears in listings of Be stars (eg. for the W1, W2, W3 and W4 bands (3.4, 4.6, 12, and 22µm re- Oudmaijer & Drew 1997) and Herbig Be stars alike (Fairlamb et spectively) for most of the HAeBes. Note that for some sources al. 2015). The latest dedicated study puts the object in the Herbig some of the bands may be missing or just be upper limits. We Be category (Patel et al. 2017). To assess the effect of ambiguous double-checked all infrared matches with the dereddened opti- classifications in our study we will also, next to the full sample, cal photometry and found no inconsistencies. consider the subset of Herbig Ae/Be stars in Table 1 of The et al. (1994). This catalogue contains all historically known, and best studied, Herbig Ae/Be stars. 98/254 of our initial sources with 2.4. Hα equivalent width and emission line profile parallaxes are present in this table (their best candidates). We collected all the Hα equivalent widths (EW) we could find in the literature for the Herbig Ae/Be stars. Not only the inten- 2.2. Atmospheric parameters, photometry and extinction sity of the line but also the shape contains very useful informa- values tion. Therefore, when possible, information about the shape of We obtained atmospheric parameters and photometric and ex- the Hα line was included. We have classified the Hα line profile tinction values for all the sources from the literature. These were as single-peaked (s), double-peaked (d) and showing a P-Cygni mainly Alecian et al. 2013; Carmona et al. 2010; Chen et al. profile (P), both regular or inverse. EW and line shape informa- 2016; Fairlamb et al. 2015; Hernández et al. 2004; Hernández et tion are presented in Table 3 and Table 6 for the high quality and al. 2005; Manoj et al. 2006; Montesinos et al. 2009; Mendigutía low quality samples respectively. Many Herbig Ae/Be stars are et al. 2012; Sartori et al. 2010; Vieira et al. 2003 and the APASS quite variable in their Hα emission and its EW may significantly change at short timescales (e.g. Costigan et al. 2014). This is also Data Release 9. Whenever the effective temperature (Teff) was not available it was derived from the spectral type with the effec- the case for the line shape, although spot checks on objects that tive temperature calibration tables of Gray & Corbally (2009). A have more than one Hα observation listed in the literature ap- 1 subspectral type uncertainty was assigned in all cases. When pear to indicate that there are not many changes in line profile classification (see also for example Aarnio et al. 2017), although not listed in the literature, AV values were derived from the ob- served photometry and using the intrinsic colours of Pecaut & changes between single-peaked and double-peaked profiles in a given star are also observed (Mendigutía et al. 2011a). We do Mamajek (2013). An RV = 3.1 was used in all cases in which AV was derived although other studies like Hernández et al. (2004) note the additional complication that emission line shapes are or Manoj et al. (2006) have suggested that a larger value of for often difficult to unambiguously classify. Regarding the Hα EWs compiled, we note that our main ref- example RV = 5 could be more appropriate for HAeBes where local extinction dominates the total extinction. This is a topic for erences (Fairlamb et al. 2015 and Mendigutía et al. 2011a) pro- future investigations using Diffuse Interstellar Bands (as done vide the non-photospheric contribution of the EW, while most by for example Oudmaijer et al. 1997). The relevant data of each other authors state the observed EW, which includes the pho- source is presented in Table 2 and Table 5 for the high quality tospheric contribution. This photospheric absorption peaks for and low quality samples respectively. A0-A1 type objects, with EW values of ∼ +10Å (See e.g. Fig. HAeBes usually show photometric variability. Thus, for ob- 7 of Joner & Hintz 2015) but is only ∼ +2Å for B0 objects. We jects with multi-epoch photometry available, we selected the used the Joner & Hintz (2015) results to correct those EWs that brightest set to determine the extinction towards the objects and were not corrected for absorption. thus their intrinsic brightnesses. As we will also show later, the We have Hα EWs for 218/252 of the HAeBes and line variability is often caused by irregular extinction, using those profiles for 197 of these: 31% are single-peaked, 52% double- data with minimum extinction introduces the smallest errors in peaked and 17% P-Cygni (of which the vast majority are of the determination of the stellar parameters. For this reason, we regular P-Cygni type). This is in agreement with Finkenzeller only used simultaneous photometry when deriving AV values. & Mundt (1984) who found that out of 57 HAeBes, 25% were All the photometric values were corrected for extinction using single-peaked, 50% showed double-peaked Hα profiles and 20% the reddening law of Cardelli et al. (1989). presented a P-Cygni profile (both regular and inverse). The main Two sources, V833 Ori and GSC 1829-0331, do not have references for the EW values are Baines et al. (2006); Fairlamb enough simultaneous photometry available to derive extinctions et al. (2017); Hernández et al. (2004); Mendigutía et al. (2011a) for them and hence they were excluded for the sample. The total and Wheelwright et al. (2010). Main references for the line pro-

Article number, page 3 of 39 A&A proofs: manuscript no. Vioque_et_al_2018

6.0 5.5 5.0 4.5 20M 4.0 14M ] 3.5

L 9M [ 3.0

)

L 2.5

( 5M

g 2.0 3.4M o l 1.5 2M 1.0 0.5 1.5M 0.0 Binary 1M -0.5 4.6 4.5 4.4 4.3 4.2 4.1 4.0 3.9 3.8 3.7 log(Teff) [K] Fig. 1. 223 Herbig Ae/Be stars in the HR diagram satisfying Eq. 1 constraint. In most cases vertical error bars are dominated by parallax uncer- tainties. Sources with a white dot have been classified as binaries. The mass of each Pre-Main Sequence track (Bressan et al. 2012) is indicated on the righthand side. An isochrone (Bressan et al. 2012; Marigo et al. 2017) of 2.5 Myr is also shown for reference as a dashed line.

6.0 6.0 5.5 5.5 5.0 5.0 4.5 20M 4.5 20M 4.0 14M 4.0 14M ] ] 3.5 3.5

L 9M L 9M [ 3.0 [ 3.0

) )

L 2.5 L 2.5

( 5M ( 5M

g 2.0 3.4M g 2.0 3.4M o o l 1.5 l 1.5 2M 2M 1.0 1.0 0.5 29 low quality objects 1.5M 0.5 1.5M 0.0 223 high quality objects 1M 0.0 In The et al. 1994 sample 1M -0.5 -0.5 4.6 4.5 4.4 4.3 4.2 4.1 4.0 3.9 3.8 3.7 4.6 4.5 4.4 4.3 4.2 4.1 4.0 3.9 3.8 3.7 log(Teff) [K] log(Teff) [K] Fig. 2. Left: 223 high quality and 29 low quality Herbig Ae/Be stars in the HR diagram after the cut in astrometric quality described in Eq. 1. Right: 218 Herbig Ae/Be stars in the final high quality sample after removing the 5 problematic objects described in Sect. 3.1. In red those objects present in Table 1 of The et al. (1994). The mass of each Pre-Main Sequence track (Bressan et al. 2012) is indicated on the righthand side. An isochrone (Bressan et al. 2012; Marigo et al. 2017) of 2.5 Myr is also shown for reference as a dashed line.

files are van den Ancker et al. (1998); Baines et al. (2006); and Baines et al. (2006) who performed spectro-astromety of 45 Mendigutía et al. (2011a); Vieira et al. (2003) and Wheelwright HAeBes and 31 HAeBes respectively, probing companions in et al. (2010). The rest of the references can be found in Table 3 the ≈ 0.1 - 2 arcsec range, and Leinert et al. (1997) who per- and Table 6. formed speckle interferometry of 31 objects, sampling separa- tions of order 0.1 arcsec. 81/252 HAeBes (∼ 32%) of our set are catalogued as binary systems, a fraction in agreement with the 2.5. Binarity Duchêne (2015) findings if we take into account the large num- ber of faint Herbig Ae/Be stars which have never been studied More than half of the Herbig Ae/Be stars are known to be in bi- for binarity. The binary status of each HAeBe is presented in Ta- nary systems (Duchêne 2015). The true number is likely much ble 2 and Table 5 for the high quality and low quality samples larger, as there have been a small number of targeted surveys for respectively; main references were Baines et al. (2006); Leinert binarity of HAeBe stars, the largest are Wheelwright et al. (2010)

Article number, page 4 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars et al. (1997) and Wheelwright et al. (2010), we refer to Table 2 MWC 314, MWC 623 and MWC 930 on the other hand ap- and Table 5 for a complete list. pear quite luminous and very much to the right of the Main Se- Baines et al. (2006) found a typical wide (few hundred au) quence, something very unusual for high mass PMS objects. An separation in the binary systems. Wheelwright et al. (2010) de- individual inspection reveals that these objects are more likely to tected no binaries closer than 30 au and established a range of be evolved giants and they appear in the literature as such (e.g. ≈ 40 - 4000 au in their data. for MWC 314: Carmona et al. 2010, for MWC 623: Lee et al. 2016, for MWC 930: Miroshnichenko et al. 2014). Deciding on the nature of the various Herbig Ae/Be candidates in our master 3. Derived quantities sample is not our intention and it is beyond the scope of this pa- per that is essentially a statistical study. However, these objects 3.1. Luminosity and Hertzsprung-Russell diagram occupy a special place in the HR diagram that is consistent with Using the parallaxes, atmospheric parameters and extinction val- both a pre- as well as a post-Main Sequence nature while there is ues we derived the luminosity for the 252 HAeBes with par- much information regarding these objects supporting their post- allaxes employing a similar method to Fairlamb et al. (2015), Main Sequence nature. We therefore decided to err on the cau- which is similar to Montesinos et al. (2009) and van den An- tious side and remove these as well from further analysis. cker et al. (1997). In short, it first consists of using values of The final HR diagram without these 2+3 problematic objects Teff and surface gravity (log(g)) to select an atmosphere model is presented in Fig. 2 on the right. In addition, in this graph, we from Castelli & Kurucz (2004) (referred to as CK-models here- highlight the sample of The et al. 1994 bonafide HAeBes in red. after) for each star to be used for its intrinsic spectral energy This final high quality sample of 218 objects will be the one distribution (SED). Solar metallicity CK-models were used in we will use in the following plots and studies. The information all cases but for BF Ori, RR Tau, SV Cep, XY Per and WW Vul concerning the 34 discarded objects in the low quality sample for which the metallicities are known not to be solar from the can be found in the tables at the end of the paper. spectroscopic work of Montesinos et al. (2009). When possible, the log(g) values were estimated from the luminosity class; oth- In this last high quality HR diagram we see that there are erwise they were taken as 4.00 (typical values range from 3.5 to many more low mass HAeBes than high mass HAeBes (69% 4.5). Uncertainties in log(g) and metallicity can be neglected in of the sources are below 4M ). This is most likely because of our study as their effect on the model SED and derived quantities the Initial Mass Function (IMF). This trend of more objects for ∼ is negligible. lower masses discontinues below 2M . This is roughly the We then scaled the model to the dereddened photometric mass corresponding to the boundary between Main Sequence A- Johnson V band. The energy distribution is then integrated over and F-type stars, and thus the traditional lower mass boundary at frequency to get the total flux. The final luminosities, presented which the Herbig Ae/Be stars were originally selected. in Table 2 and Table 5 for the high quality and low quality sam- For lower masses the sources are more spread in tempera- ples respectively, are then obtained by means of the total flux and ture, occupying larger parts of the PMS tracks, while, instead, the parallax. All sources of uncertainty were taken into account the high mass objects tend to be predominately located close to at this step including using different CK models for the different the Zero Age Main Sequence (ZAMS). This is likely because the temperatures within the T eff uncertainty range. higher the mass the faster the PMS evolution is. This fast evolu- The 223 Herbig Ae/Be stars satisfying Eq. 1 constraint are tion could explain why high mass objects at low temperatures plotted in the resulting HR diagram in Fig. 1. This number is an (and thus low surface gravities) are hardly present in Fig. 1 or increase of more than a factor of ten compared to the previous, the sample. Hipparcos-based, study by van den Ancker et al. (1998). PMS We will encounter more examples below where high mass evolutionary tracks from Bressan et al. (2012) are also plotted and low mass objects have different properties. in Fig. 1 plus a 2.5Myr isochrone (Bressan et al. 2012, Marigo et al. 2017), all of them with solar metallicities (Z = 0.01 and Y = 0.267). 3.2. Mass and age Before we analyze this sample, we also plot the HR diagram for all 252 objects with parallaxes (the high and low quality sam- Using the isochrones, the mass and age of the Herbig Ae/Be stars ples together, hence including those which failed the Lindegren were estimated. We used a hundred PARSEC isochrones with so- quality selection criteria) in Fig. 2 on the left. Many of these lar metallicity (Bressan et al. 2012; Marigo et al. 2017) from 0.01 badly astrometrically behaved sources are located in unphysical to 20Myr, from which we only use the PMS tracks. To each Her- positions, significantly below the Main Sequence, validating our big Ae/Be star we assigned the closest two isochrone points in approach of removing those from our analyses. the HR diagram; the solar metallicity isochrones did not match Returning to the HR diagram in Fig. 1, there are still several seven sources from the high quality sample in the HR diagram outliers that do not seem to be PMS objects. GSC 5360-1033 and isochrones with lower metallicities were used in those cases. and UY Ori appear way below the Main Sequence, just like the As each point is associated with a mass (M) and an age, we com- lower quality objects that were removed earlier. However, the puted for each HAeBe an average of those values weighted by Gaia DR2 data of these two objects appear of good quality. Re- the distance to the points. The result is an estimate of age and garding GSC 5360-1033 the situation is not clear, an ambiguous mass for 236/252 HAeBes. These values are presented in Ta- spectral type or photometry for this object, or an incorrect esti- ble 3 and Table 6 for the high quality and low quality samples mation or the extinction may be the reason for the unexpected respectively. Uncertainties were derived from the error bars in location of the object. For UY Ori, Fairlamb et al. 2015 assigned the HR diagram (Fig. 1 and Fig. 2) keeping a minimum error a spectral type of B9 to this object, but the photometry listed in of 5%. We compared many of our masses and ages with those simbad indicates a large variability. Pending more certainty, we of Alecian et al. (2013) and Reiter et al. (2018). We found that decided to remove these two objects from the high quality sam- our determinations of these parameters are consistent with the ple and place them in the low quality sample. results of the previous authors.

Article number, page 5 of 39 A&A proofs: manuscript no. Vioque_et_al_2018

3.3. Infrared excesses In the process of deriving the luminosity it is also possible to de- 0.4 rive the infrared excess. We have logarithmically interpolated the Herbig Ae/Be stars UBVRI standard stars different dereddened observed fluxes from the J band (1.24µm) Faint UBVRI standard stars to the W4 band (22µm) and defined the infrared excess (E) as: 0.3

(F − F∗) 0.2 E = e [λ1,λ2] (2) F∗

Fe is the total flux underneath the observed dereddened pho- 0.1 tometry (the infrared photometry has also been dereddened) and Frequency density F∗ is the total photospheric flux below the CK model. λ1 and λ2 define the range of wavelengths of interest and the total fluxes 0.0 in the numerator just refer to that range. This measure expresses 0 5 10 15 20 25 30 35 40 45 50 55 60 V the excess in terms of the total luminosity of the object. For ex- i ample, all things being equal, if we have two stars with the same Fig. 3. Distribution of the variability indicator for Herbig Ae/Be stars amount of dust surrounding them, with one of them brighter, the and two catalogues of photometric standards; one of bright sources infrared re-radiated emission will be larger, but the IR excess (Landolt 2009) and one of faint sources (Clem & Landolt 2016). As as defined here, would be the same as it is a relative measure. a class, the Herbig Ae/Be stars are more variable than the photometric The same or a very similar indicator was used by Cote & Waters standards. (1987, their Eq. 8), Waters et al. (1987, Eq. 3), and more recently by Banzatti et al. (2018) in their Sect. 2.3. of all Gaia objects in a brightness interval of ±0.1 magnitude Uncertainties in the infrared excesses were derived using the around the G band value of the object (i.e. to Aa,G ∈(a ,a ), with uncertainties in the observed fluxes and the uncertainties in the a 1 2 a indexing the Gaia catalogue and being a1 = Gi − 0.1mag and temperature (which affect the CK models) of each object. a2 = Gi + 0.1mag). The equation is as follows: We have split the total infrared excess in two, a Near Infrared Excess (1.24 − 3.4µm, roughly the 2MASS region) and a Mid − − Infrared Excess (3.4 22µm, the WISE region). The values for Ai Aa,Ga∈(a1,a2) these excesses are presented in Table 3 and Table 6 for the high Vi = (4) σ[A ] ∈ quality and low quality samples respectively. The total infrared a Ga (a1,a2) excess (1.24 − 22µm) is the sum of the two. G is the Gaia white G band magnitude and σ is the standard In addition, we also computed the infrared excess at each deviation. In essence, we subtract the error to flux ratio of each individual band (J, H, Ks, W1, W2, W3 and W4) as the flux ratio HAeBe, weighted by the number of observations, to the mean of between the dereddened observed monochromatic flux and the the same expression (Aa, Eq. 3) for the sources in the Gaia cata- expected flux according to the CK model. The values for these logue within ±0.1mag of the Herbig star in the G band. Then we excesses are presented in Table 4 and Table 7 for the high quality divide by the standard deviation of Aa s of that Gaia subset. This and low quality samples respectively results in a variability indicator which measures the variability (in standard deviations, σ) for each Herbig Ae/Be star compared 3.4. Variability information to the mean of field objects of the same brightness. For completeness, we note that it is necessary to impose Gaia DR2 does not provide a general variability indicator for more constraints to exclude the cases in which a larger error is all sources. Here, we use Gaia’s repeated observations to extract not due to intrinsic variability. Following Deason et al. (2017), photometric variability information. Gaia DR2 used a total of 22 Appendix C of Lindegren et al. (2018) and what was done in months of observations and each source was observed repeatedly Gaia Collaboration et al. (2018) we require Nobs > 70 and more in a non periodic fashion. Data Release 2 provides the average than 8 visibility periods (i.e., groups of observations separated photometry, the uncertainty on this value and the number of ob- by at least four days), plus the Eq. 1 constraint that limits the as- servations. All things being equal, the photometric “error” will trometric quality (and hence the variability indicator just can be be larger for a photometrically variable object than for a stable derived for sources in the high quality sample). In order to also object. Here we aim to quantify the variability of the objects. We limit the photometric quality we included the following criterion start with the “Variability Amplitude” (Ai) for a certain source i presented in Gaia Collaboration et al. (2018): as presented in Deason et al. (2017): − 2 − 2 p 1.0 + 0.015(GBP GRP) < EF < 1.3 + 0.06(GBP GRP) (5) Ai = Nobs,i e(Fi)/Fi (3) where EF is the flux excess factor and GBP and GRP the Gaia where Nobs is the number of CCD crossings, F and e(F) are blue and red passbands respectively. Note that these constraints the flux and flux error respectively. This quantity is powerful may inevitably exclude many of the very variable HAeBes as in identifying objects that show larger flux variations than ex- they also trace larger errors and hence variability. These con- pected for a stable star. However, in order to statistically assess straints will also be biased to discarding binaries and faint the level of variability, we introduce a variability indicator Vi, sources in crowded areas (Lindegren et al. 2018; Gaia Collab- which quantifies how much more variable an object is compared oration et al. 2018). to stable objects of the same brightness. In short it compares the The variability indicator values for the 193 sources satisfying Variability Amplitude from Eq. 3 of a given object (i) to that the previous conditions are presented in Table 3.

Article number, page 6 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars

3.75 3.75 Age 0.25Myr Teff 8150K (up to A5) 3.50 3.50 0.25Myr < Age 1.25Myr 8150K < Teff 9500K (A5-A1) m m 3.25 1.25Myr < Age 4Myr 3.25 9500K < Teff 10700K (A1-B9)

2 3.00 4Myr < Age 7.5Myr 2 3.00 10700K < Teff 16500K (B9-B3) 2 2 16500K < Teff (hotter than B3) o 2.75 7.5Myr < Age o 2.75 t t

4 2.50 4 2.50 2 2 . .

1 2.25 1 2.25

m 2.00 m 2.00 o o r r f f

1.75 1.75 s s

s 1.50 s 1.50 e e c c

x 1.25 x 1.25 e e 1.00 1.00 d d e e

r 0.75 r 0.75 a a r r f 0.50 f 0.50 n n I I 0.25 0.25

0.00 0.00 101 10 2 10 1 100 101 Mass [M ] Age [Myr]

Fig. 4. Left: IR excess in the range 1.24 − 22µm vs. estimated mass of the objects. The most massive objects (more massive than ∼ 7M ) barely show an infrared excess. Right: IR excess in the range 1.24 − 22µm vs. estimated age. Ages and effective temperatures are respectively colour coded in the legend. The symbols stand for the Hα line profiles: circles (double-peaked), triangles (single-peaked), stars (P-Cygni profile) and diamonds (no information). Note that although it is not necessarily a one-to-one correlation, lower ages correspond to higher masses.

In Fig. 3, we show the Vi distribution of Herbig Ae/Be stars variability indicator values larger than 2. We would therefore ex- and compare it to the Vi distribution of bright photometric stan- pect only 3 of our 193 Herbig Ae/Be stars for which we could dards from Landolt (2009) and faint photometric standards taken determine Vi to be strongly variable, at amplitudes of 0.5 magni- from Clem & Landolt (2016). If Eq. 4 would have not been used tudes in the V band or higher. However, we find 45 more, indi- these two latter samples would have had a different mean in cating that a large fraction of Herbig Ae/Be stars exhibit strong the distribution of Ai. The Herbig Ae/Be stars appear to show, variations. on average, a larger variability indicator value than the standard In addition, it is interesting to compare our variability in- stars, being the break at ∼ Vi = 2. We performed a two-sample dicator with the variability catalogues published alongside the Kolmogorov-Smirnov (KS) statistical test to study whether Her- Gaia DR2 general catalogue (Holl et al. 2018). Just 1 every 3000 big Ae/Be stars can be drawn from those two samples of stan- objects passed the Gaia DR2 stringent selection criteria for vari- dard stars. The result shows that we can reject that hypothesis to ability. 10/252 objects in our list fall in this category and appear within a 0.001 significance and hence this variability indicator as variable in those catalogues. Of the 5 of those that have de- differentiates them as a group. rived Vi values they are larger than Vi = 5. In order to assess the relation between our variability indica- tor (G band variability) and variability in the V band we com- pared the magnitude variations in the V band as presented in the International Variable Star Index VSX (Watson et al. 2006) 4. Data analysis with our variability indicator values. We found that we are trac- 4.1. Infrared excesses ing variabilities as small as ∼ 0.5mag with the Vi = 2 cut-off. In Eiroa et al. (2002) 7/23 (30%) PMS objects homogeneously ob- In Fig. 4 the total infrared excess (1.24−22µm) vs. the estimated served for variability have variabilities above 0.5mag. In our case mass and age of the sources is plotted. There appears a difference 48/193 sources have values above Vi = 2 (25%) and hence can be in IR properties between high and low mass stars. Whereas low considered as strongly variable. Of those 48, 17 are catalogued mass stars show a range of infrared excess, the higher mass stars as UXOR type (Mendigutía 2011c; Oudmaijer et al. 2001; Poxon in general only present very low levels of excess. A similar be- 2015). There are 5 other UXORs in our sample with Vi values, haviour is seen when the excess is plotted as function of age; 4 of them have reported optical variabilities smaller than 0.5mag the excess for the youngest objects is smallest. This is probably in the V band. The other one is BO Cep. This object has been re- readily explained by the fact that the more massive PMS objects ported to have a periodic variability with a single prominent peak in the HR diagram have the lowest ages by virtue of their rapidly with a period of ∼ 10 days (Gürtler et al. 1999). The regular non evolving isochrones, so trends in mass will automatically also be periodic variability of the object is smaller than 0.5mag which present in those with age. To study trends as a function of age, it explains why this UXOR has not been detected by our variability would be necessary to consider subsamples with a narrow range indicator. Supporting this, it appears as UXOR in Poxon (2015) in mass. We therefore consider that the main result of this exer- but not in Oudmaijer et al. (2001) or Mendigutía (2011c). cise is that high mass objects have a very low infrared excess, To put the variability indicator into perspective, we find that and that there appears to be a break at ∼ 7M from where almost 6 out of 411 photometric standards from Landolt (2009) have no sources with significant excess appear.

Article number, page 7 of 39 A&A proofs: manuscript no. Vioque_et_al_2018

each band (the values for these excesses are presented in Table 4 and Table 7 for the high quality and low quality samples respec-

m tively). In all cases, there is a general and consistent increase of 0 2 10 the Hα EW from sources with very little IR excess to those with 2

o higher IR excess. t

4 1 . 10 In Table 1 we show that the Hα emission line equivalent 3 width is more correlated with the infrared excess at shorter wave- m

o lengths than at larger wavelengths, with the correlation peaking r 2 f 10 at 2.16µm (K band). s s s

e An obvious question might be whether there is a causal cor- c 3 x 10 Single-peaked H profile relation between the Hα emission and presence of emission due e Double-peaked H profile d to dust around these objects. The various excesses at various e r P-Cygni H profile wavebands are correlated with each other (Fig. 5), and as a con- a 4

r 10 f No line shape information sequence the IR excess at many wavelengths also correlate with n I 10 3 10 2 10 1 100 the EW. However, the correlation with Hα is strongest at the Ks Infrared excess from 1.24 to 3.4 m band which traces the hot dust in the inner disk, suggesting the accretion mechanism or wind activity as traced by Hα is related Fig. 5. IR excess in the range 3.4 − 22µm (Mid IR excess) vs. IR excess to the inner parts of the dusty disk (see also Manoj et al. 2006). in the range 1.24 − 3.4µm (Near IR excess). The symbols stand for As presented in Table 1, the correlation rises from a minimum at the Hα line profiles: circles (double-peaked), triangles (single-peaked), 1.24µm (effectively tracing the stellar photosphere) up to 3.4µm stars (P-Cygni profile) and diamonds (no information). A linear fit in the and then goes down again to the same minimum at 22µm (W4 log space is shown in blue (log(Mid IRexcess) = 1.16 log(Near IRexcess) + 0.23, r = 0.88). band), where dust in the outer disk is found. In fact, Mendigutía et al. (2012) discovered the same correlation between IR excess and accretion rate and they found that it is no longer present be- yond 20µm. Fig. 5 splits the total infrared excess into two, a near-infrared For comparison purposes, in Fig. 7 the Ks band is plotted in and a mid-infrared part. It demonstrates that the excess at both the upper panel and the W4 band in the lower. It is noteworthy wavelength ranges are highly correlated with each other (the lin- that for the Ks band, where we have the strongest correlation, ear fit in logarithmic space that can be seen in the plot has a cor- small EWs are almost only present in sources with little IR ex- relation coefficient of r = 0.88). Therefore, it is not unexpected cess and, in consonance with Sect. 4.1, for a given Hα EW value that the ∼ 7M break is also present at near- and mid-infrared. low mass stars (M < 7M ) tend to have higher IR excesses. However, these trends are weaker or non existent in the case of the W4 band, where we have the weaker correlation. This 4.2. Hα equivalent width reinforces the idea that the Hα emission is correlated with the inner parts of the disk. Note that in both panels the average ex- Fig. 6 shows the EW as a function of mass and age respectively. cess is still lower for the higher mass objects. The emission line As the definition of a Herbig Ae/Be star includes the presence strengths will also be subject of a follow-on study using accre- of emission, which is mostly from the Hα line, it may not come tion rates (Wichittanakom et al. in preparation). as a surprise that essentially all measurements are negative (i.e. tracing emission). The EWs show a large range of values, which appears to in- 4.4. Variability crease with increasing mass and decrease with increasing age We conclude this section by studying the variability of the ob- (studied by Manoj et al. 2006). The older objects typically have jects and its correlation with the various properties discussed so lower EW’s than younger objects. It is tempting to read an evo- far, including the Hα line profiles we took from the literature. lutionary effect in this finding - after all it would be expected that The left panel of Fig. 8 presents the variability indicator as the accretion (and therefore emission) would decrease when the function of the total (near plus mid) IR excess. As described in PMS objects are closer to the MS. However, we should recall Sect. 3.4 the variability indicator states the number of standard that there is a strong correlation between the age and the mass deviations a certain source is separated from the mean of the of the stars, so we may well be looking at a mass effect instead. Gaia objects of the same brightness. No, or hardly any variabil- As the EW is a relative measurement with respect to the stellar ity is present at the lowest IR excesses but sources can be both continuum, a larger EW for otherwise similar objects indicates a variable and non-variable at the higher IR excesses, consistent stronger emission line. The observed trend towards higher tem- with van den Ancker et al. (1998) based on a smaller sample. peratures/masses and thus higher luminosities implies that the The figure next to it shows the variability as a function of lines become even stronger than the EW alone would seem to mass. As high mass stars in this sample generally do not have a imply. strong IR excess, we find that mostly the lower mass and cooler objects display high variabilities, with the break also around 7 4.3. Hα EW and infrared excess M , corresponding to a Main Sequence spectral type of around B3. Although cooler objects tend to have larger variabilities (also The correlation between Hα emission, measured by its equiva- observed by van den Ancker et al. 1998), we can observe how the lent width, and near- and mid-infrared excess is studied in Fig. 7 range in temperatures for variable sources is wide in the right and Table 1 for each one of the infrared bands (J, H, Ks, W1, panel of Fig. 8, and that there are in fact many Herbig Be stars W2, W3 and W4). In this case, we computed the infrared excess with very strong variabilities. Hence, this is more likely a trend as the flux ratio between the dereddened observed monochro- with mass and not with temperature. We do note that although we matic flux and the expected flux according to the CK model at detect photometric variability from the Vi = 2 value, the Vi = 5

Article number, page 8 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars

-600 -600 Teff 8150K (up to A5) Mass 1.8M

8150K < Teff 9500K (A5-A1) 1.8M < Mass 2.2M

9500K < Teff 10700K (A1-B9) 2.2M < Mass 3.4M

-500 10700K < Teff 16500K (B9-B3) -500 3.4M < Mass 6M

16500K < Teff (hotter than B3) 6M < Mass

-400 -400 ] ] Å Å [ [

W -300 W -300 E E

H H

-200 -200

-100 -100

0 0 101 10 2 10 1 100 101 Mass [M ] Age [Myr]

Fig. 6. Left: Hα EW vs. estimated mass. Right: Hα EW vs. estimated age. Effective temperatures and masses are respectively colour coded in the legend. The symbols stand for the Hα line profiles: circles (double-peaked), triangles (single-peaked), stars (P-Cygni profile) and diamonds (no information). value is a better separation boundary for the observed trends in sible due to the more than a factor of ten increase in available both panels of Fig. 8. distances to these objects compared to Hipparcos. With the Gaia The challenge is to identify which property lies at the cause DR2 data, the majority of known Herbig Ae/Be stars could be of the variability, is it the mass of the objects, their age or in- placed in the HR diagram. We found the following: frared excess emission or something else? An important clue is that many objects with strong variability (above Vi = 2) and – There are more low mass objects than high mass objects, line shape information have doubly peaked Hα profiles ( 31 out with the high mass objects mostly located close to the Main of 43; 72±7%, 68% confidence interval). In general, double- Sequence. peaked emission line profiles are due to rotating disks, so the data – High mass objects have in general very small infrared ex- are suggestive of an edge-on disk-type orientation and structure cesses and low variability, the properties appear to differ (from the remaining 12 objects they all have a P-Cygni profile around 7M . and none have a single-peaked profile). The number of variable – Hα emission is generally correlated with infrared excess, objects with doubly peaked line profiles is significantly differ- with the correlation stronger for IR emission at wavelengths ent from the full sample, in which only half of the targets with tracing the hot dust closest to the star. known line classifications have a double-peaked profile (of the – More massive and younger objects have higher Hα EWs. sources with derived variability indicator and known line pro- – When split at 7M into “low” and “high” mass samples, the file 79 out of 155; 51±4% are double-peaked and 48 our of 155; Hα - IR excess correlations hold for both mass ranges, with 31±4% are single-peaked). These fractions are significantly dif- the average excess lower for the higher mass objects. ferent, and we therefore suspect that the variable sources are – Photometric variability can be traced back to those objects mostly oriented edge-on, and that the line-of-sight inclination with double-peaked Hα emission and large infrared ex- to the objects could be a decisive factor in the cause of the vari- cesses. ability. This is in agreement with the trend observed in the left – All catalogued UXORs in the sample with detected variabil- panel of Fig. 8. Sources with large amounts of circumstellar ma- ities above 0.5mag in the V band appear as strongly variable terial show large infrared excesses and high or low levels of vari- (above Vi = 2) with the exception of BO Cep (discussed in ability depending on the inclination of their disk whilst sources Sect. 3.4). with little material around have low infrared excesses and low variabilities in all cases (also discussed in van den Ancker et al. Here, we will discuss these findings and their implications 1998). for the formation of intermediate mass stars.

5.2. Selection effects 5. Discussion Let us first investigate the various selection effects and biases 5.1. General findings that could potentially affect the results. In the above we have determined fundamental parameters such Quality parallaxes: It could be argued that the quality of the as- as temperature, mass, age, IR excess, variability and luminosity trometric data has an effect on the findings. The parallax errors for a large sample of Herbig Ae/Be stars which was made pos- occupy a comparatively small range, from ∼0.016-0.37 mas, but

Article number, page 9 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 because of the large spread in distances, the relative uncertainties can be very large. To investigate whether this has a detrimental effect on the results, we repeated the analysis with only the ob- 102 jects with the very best parallaxes ($/σ$ > 10). This, of course, M 7M (2.16 m) limits the sample and 182/218 objects remain in the high qual- M < 7M (2.16 m) ity sample. These 182 objects are less luminous, which may be expected as in general they have larger parallaxes and are thus 101

closer. As a result they will be less massive and have larger ages ] K than the objects in the entire sample. This, as a consequence of C F / the trends described in previous sections, implies that these ob- d 0

e 10 v jects also show larger infrared excesses and variabilities as well r e s

as smaller Hα EWs (see Fig. 4, Fig. 6 and Fig. 8). However, we b o

find that essentially all correlations also hold for the higher qual- F [ 5 0 101 2 10 M 7M 1(022 m) 10 ity parallax sample, and if anything, they appear stronger. For s

s M < 7M (22 m) example almost all of the high mass sources that have large in- e 104 c

frared excesses and variabilities in figures Fig. 4 and Fig. 8 have x e $/σ$ < 10. The inclusion of lower quality parallaxes induces 103 R an extra scatter in the results, but the larger sample and wider I coverage in luminosity aids in reinforcing them. 102

Quality identification as Herbig Ae/Be star: Another potential 101 source of error is source misclassification. We have used the 100 largest sample of Herbig Ae/Be stars published to date (Chen et 0 1 2 al. 2016 with some added from Alecian et al. (2013); Baines et 10 10 10 al. (2006); Carmona et al. 2010; Fairlamb et al. 2015; Hernández H EW [Å] (absolute value) et al. (2005); Manoj et al. (2006) and Sartori et al. (2010)). The defining characteristics of HAeBes are not unique to the class, Fig. 7. 2.16µm (blue markers) and 22µm (red markers) infrared excesses and can often also be found in other types of stars such as clas- defined as Fobserved/FCK vs. Hα equivalent width (absolute value). Note sical Be stars, which display Hα emission and a near-infrared that this IR excess indicator is a flux ratio and not the one described in excess (e.g. Rivinius et al. 2013) or evolved stars which can Eq. 2 where we integrated under the SED. Dots are Herbig Ae/Be stars have spectral types A and B, display hydrogen recombination with M < 7M and triangles are Herbig Ae/Be stars with M > 7M . emission and be surrounded by dusty shells and disks such as Lines are linear fits to the data, dashed for HAeBes with M > 7M the Luminous Blue Variables and B[e] stars (Davies et al. 2007; and in solid colours for HAeBes with M < 7M ; black solid lines are the linear fits for all the sources (equations and correlation coefficients Oudmaijer, et al. 1998 on HD 87643). It is therefore inevitable for these fits to all the sources for all the infrared bands can be seen in that some sources will have been misclassified. It would be fair Table 1). Note the difference in the scale of the vertical axis between the to say that the more “classical” Herbig Ae/Be stars going back to two panels. the Herbig (1960) and The et al. (1994) papers have been stud- ied in more detail and are better established as young pre-Main Table 1. Correlation between IR excess and Hα EW at different wave- Sequence stars. lengths. We therefore studied the The et al. sample of objects (their Table 1, 85 sources out of our 218) separately and find that all Band Correlation A B correlations do hold for this “gold standard” sample as well. We coefficient (r) do find that on average these objects have a larger Hα EW and J (1.24µm) 0.41 0.15 ± 0.03 0.025 ± 0.034 have larger IR excesses than the full sample. These properties are H (1.66µm) 0.56 0.32 ± 0.03 0.0024 ± 0.0478 the defining characteristics of a Herbig Ae/Be star, and it may Ks (2.16µm) 0.60 0.48 ± 0.05 0.046 ± 0.066 not be surprising that the first objects to be proposed as Herbig W1 (3.4µm) 0.57 0.64 ± 0.07 0.15 ± 0.10 Ae/Be stars are on average more extreme in these properties. Yet, W2 (4.6µm) 0.57 0.78 ± 0.09 0.24 ± 0.12 again, as with the higher quality parallax sample, the trends are W3 (12µm) 0.52 0.93 ± 0.12 0.79 ± 0.16 still present in this sub-sample. W4 (22µm) 0.41 0.71 ± 0.12 2.05 ± 0.17 Mass distribution of the sample: The known Herbig Ae/Be stars have mostly been found serendipitously, and a large-scale sys- Notes. Correlation between IR excess (defined as a flux ratio, Fobserved/FCK) and Hα EW at different wavelengths for all the sources. tematic search for them has yet to be carried out. Yet, an inter- The coefficients are defined by: log(F /F ) = A log(|EW|) + B. esting question is how representative the present sample is for observed CK The Ks band, with the higher correlation, and the W4 band are in bold; the class. To this end, we consider the mass distribution of the both are shown in Fig. 7. objects. There are more or less the same number of low mass, A-type objects than higher mass B-type objects. There are more Herbig Be stars than might be expected from the Initial Mass goals is to draw an increased and well-selected sampled of Her- Function, however, the B-type objects are brighter and are sam- big Ae/Be stars from the Gaia catalogues. pled from a larger volume, as also attested by their smaller par- allaxes. We will therefore expect a larger fraction of Herbig Be Binarity: One may think that binarity may affect the observed stars in the sample. When limiting our sample in distance, we photometry and for example produce fake levels of variability in obtain a Herbig Ae/Herbig Be ratio that is close to the IMF. As our variability indicator. This is because binary sources tend to far as the mass distribution is concerned, we may say that the be more astrometrically and photometrically irregular. We stud- current sample is representative of the class. One of our future ied the group of binaries against the group of isolated sources

Article number, page 10 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars and overall we find that the known binaries are slightly brighter these massive objects. Thus, the observation here is consistent than the objects that have not been reported to be a binary. This with the classical scenario that the Kelvin-Helmholtz contrac- is probably a selection effect in that brighter objects were more tion timescale is much smaller for massive objects compared to likely to be included in the binary surveys. We compared the the free-fall timescale of the collapsing parental cloud. In this brightnesses of binaries and non-binaries in the Baines et al. scenario, the massive young stars only become visible once they (2006), Wheelwright et al. (2010) and Leinert et al. (1997) stud- are on, or close to, the Main Sequence - the so-called birthline. ies separately and find that within the surveys there are indeed We will discuss more on this later, but note that with this inter- no brightness differences between binaries and non-binaries. pretation one still would expect a range of infrared excesses in Returning to the Gaia sample; all other properties but in- any sample. This is consistent with what we find for massive ob- frared excesses, including variability, are similar. We do find that jects (larger than 7 M ); a large number of low excess stars, but binaries have in general slightly larger IR excesses. With the ben- still a few with noticeable excess (see Fig. 4). efit of hindsight, this is perhaps something that could have been Moving to the lower mass objects, which do display a large expected. Most of the binaries are distant binaries with separa- range of infrared excess emission, an immediate question to ask tions larger than 0.1 arcsec (Gaia’s resolution). Indeed, no bi- is whether we can detect any evolutionary effect in the sense that naries are found closer than 30au so it is expected that binarity objects that are further evolved have smaller infrared excesses, does not play a significant role in the optical photometry. At the as one expected from the progressive dust dispersal, and as sug- same time, companions could potentially contribute to the in- gested by Fuente et al. (1998). For example, if the inside out frared emission whose fluxes have been measured with apertures clearing model of disk evolution is correct, we should see a trend larger than the typical separations. Given that we do detect slight at each PMS track from high excess to little excess. differences in IR-excess between binaries and non-binaries, a However, it appears Herbig Ae/Be stars do not show any con- preliminary inference would be that the companions may con- sistent evolution of the infrared excess from high to low excess at tribute to the infrared flux in some cases. any mass range. There are many objects appearing younger than 2.5 Myr or even 1 Myr at all mass ranges with little IR excess. Arguably the lack of an evolutionary effect can be explained by the size of the error bars on for example the luminosity. The 5.3. Infrared excess as function of mass evolutionary timescales vary strongly with mass (and thus lu- Fig. 4 shows the infrared excess as function of mass and age re- minosity), masking any trend of infrared excess emission with spectively. There is a marked difference in the infrared excess age. Here, we would highlight that many young Herbig Ae stars observed towards high and low mass objects. Herbig Be stars show little excess. By looking at these objects in the right panel more massive than ∼7 M in general appear to have little to no of Fig. 2 it is not difficult to find sources with error bars small excess, while the lower mass objects show a wide range of ex- enough to discard the contribution of uncertainty to the prob- cesses. There is also a trend with age with the youngest objects lem. Finally, the contamination by binaries as discussed in Sect. having the smallest infrared excess. Although it would be tempt- 5.2 can play a role in here as many HAeBes can still remain as ing to assume a causal relation between age and presence of dust, undetected binaries. and try to explain why the youngest objects have the smallest We should also note that the underlying assumption of the amount of dust around them, we suspect the stellar mass is the evolutionary calculations is that the conditions under which the dominant factor. The durations of the PMS evolutionary tracks stars form are uniform, the accretion rates a smooth function of are progressively shorter for higher masses, and an underlying time resulting in an overall similar evolution for all stars. How- relation between mass and infrared excess would therefore also ever, the final configuration is undoubtedly affected by inhomo- appear as a correlation between age and infrared excess. geneities, varying accretion rates and even the masses of the ini- Either way, the lack of dusty emission from high mass ob- tial clouds. Nevertheless, looking for real evolutionary effects in jects is puzzling, as we might expect the more massive objects the spectral energy distributions requires selecting subsamples to be formed in more massive clouds and therefore be more em- of objects that are located at or close to the same mass tracks. bedded. A natural conclusion would be that at any time of their This may require even more precise parallaxes than can be pro- PMS evolution, these young objects would be surrounded by vided by Gaia at the moment in many cases. It also requires pre- more dust than their lower mass counterparts and therefore, at cise determinations of the atmospheric parameters and extinction any stage, they would have a stronger infrared emission. A coun- values. A proper statistical study with high quality parameters of terargument is that the Herbig Be stars are predominately found the evolutionary properties of the HAeBes as they move towards closer to the ZAMS and are therefore more evolved, having dis- the Main Sequence is hence still pending and planned for the persed their circumstellar material. Supporting this idea, Alonso- future. Albi et al. (2009) found, from their compilation of millimetre ob- servations of 44 objects, that Herbig Be stars have much weaker 5.4. Variability in terms of the UXOR phenomenon millimetre emission than their later type counterparts. In addi- tion, they found that the masses of the disks around Herbig Be The variability indicator that was developed specifically for the stars traced at mm wavelengths are usually 5-10 times lower than Gaia data demonstrates that the class of Herbig Ae/Be stars is those around lower mass stars, with the boundary also around 7 more variable than the general population of stars. Fig. 8 shows M . These authors suggest that the disk dispersal is more effi- that the lower mass objects are much more photometrically vari- cient and faster in high mass objects above 7 M . Indeed, the able than the higher mass objects, for which the variability ap- disk dispersal times are a steep, declining function with stellar pears to cease beyond ∼7 M . The photometrically variable ob- mass, from millions of years for the lower mass stars to tens of jects contain most of the so-called UXOR variables reported thousands of years for the highest mass young stars of 10 M in the literature. Using the compilation of UXOR variables by and higher (Gorti et al. 2009). Mendigutía (2011c); Oudmaijer et al. (2001) and Poxon (2015), The latter timescales are comparable to the evolutionary we find that 17 out of the 48 strongly variable objects - those timescales as for example computed by Bressan et al. (2012) for with variability indicator values larger than 2, representing vari-

Article number, page 11 of 39 A&A proofs: manuscript no. Vioque_et_al_2018

60 60

55 55

50 50

45 45

40 40 T 8150K (up to A5) 35 Single-peaked H profile 35 eff 8150K < Teff 9500K (A5-A1)

i Double-peaked H profile i 30 30 9500K < Teff 10700K (A1-B9) V P-Cygni H profile V 10700K < Teff 16500K (B9-B3) 25 No line shape information 25 16500K < Teff (hotter than B3) 20 20

15 15

10 10 V = 5 V = 5 5 i 5 i Vi = 2 Vi = 2 0 0 10 3 10 2 10 1 100 101 Infrared excess from 1.24 to 22 m Mass [M ]

Fig. 8. Left: Variability indicator vs. IR excess in the range 1.24 − 22µm. It can be seen how objects with the lower IR excess do not show high variability. Right: Variability indicator vs. estimated mass. It can be seen how the most massive objects (more massive than ∼ 7M ) do barely show variability. Line profiles and temperatures are colour coded in the legend in the left and right panel respectively. The symbols stand for the Hα line profiles: circles (double-peaked), triangles (single-peaked), stars (P-Cygni profile) and diamonds (no information). The Vi = 2 and Vi = 5 values are stressed for clarity. ations of 0.5 magnitudes (in the V band) or higher - are classified traced by our variability indicator. Given that doubly peaked line as UXORs. The remaining 5 UXORs with variability indicator profiles are most easily explained by at least part of the emission values present in the sample have documented variabilities be- originating in a rotating disk leads us to conclude that the photo- low 0.5 magnitudes with the exception of BO Cep (discussed in metrically variable objects are seen edge-on and surrounded by Sect. 3.4). a disk-like structure. It is true that outflows or winds not limited The defining characteristic of the UXOR phenomenon is not to the disk can produce double-peaked Hα profiles (Kurosawa only the photometric variability but also the reddening and blue- et al. 2006; Tambovtseva et al. 2014). Supporting the hypothesis ing associated during the variations. The explanation put forward of edge-on disks being the main cause of photometric variabil- for this behaviour is the obscuration of the star by a rotating, in- ity, we find in variability the same separation at ∼7 M between homogeneous, dusty edge-on disk. The objects first become red- low and high mass objects we found when studying IR excesses, der when dust obscures the object, and can even become blue which suggests that photometric variability and IR excess have at their faintest phases, when the direct light from the stars is the same cause. In addition, sources with high IR excesses have blocked and, predominately blue light is scattered into the line both high and low variability levels, which can be understood as of sight. As the polarization - resulting from scattered light - depending on the disk inclination, while sources with lower IR also peaks during the faintest phases (e.g. Grinin 2000), the ob- excesses show little variability in all cases (left panel of Fig. 8, scuring disk hypothesis is favoured. Interestingly, observational discussed in Sect. 4.4). This would also explain the few high evidence other than the polarization supporting this conclusion mass strongly variable objects that can be seen in the right panel has been relatively sparse. of Fig. 8; they are mostly the ones with high IR excess in the With the large sample of Herbig Ae/Be stars, and the large left panel of Fig. 4 (discussed before in Sect. 5.3). Given that an number of UXORs among them, we can repeat a similar exper- edge-on orientation is the major and main ingredient of the dust iment using the Hα line as a proxy for the inclination of the obscuration hypothesis, these results lend very strong support to circumstellar disks. We will consider the line profiles of the Hα it using a large sample of Herbig Ae/Be stars. emission in tandem with the variability indicator. Fig. 8 shows The large fraction of objects with double-peaked line pro- that all but twelve of the strongly variable objects with docu- files or variability is in agreement with the model predictions mented line profiles (above Vi = 2, those with ∆V > 0.5mag) by Natta & Whitney (2000) who worked out how many Herbig have double-peaked Hα emission. In fact, the five objects for Ae/Be stars would undergo the UXOR phenomenon considering which no line profile is listed have, to our knowledge, no re- the scale heights of dusty disks and under which inclinations the ported profiles. The occurrence of double-peaked profiles in the photometric variability would still be visible. They conclude that highly variable sample is significantly higher than for the other around half of the Herbig Ae stars could be UXORs. In our high objects (see Sec. 4.4). It is significant that the other twelve ob- quality sample we have 85 A type stars with variability indica- jects have P-Cygni profiles and none of them show a single- tor values and just 16 of them were previously listed as UXORs, peaked profile. The P-Cygni profile is often related to episodic but again most of them have been largely unstudied. However, of energetic phenomena and it is not unexpected that it is also the 25 A type stars with variabilities above Vi = 2, 13 are known

Article number, page 12 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars

UXORs. This means that for the Herbig Ae stars for which we ary PMS tracks. However, these object are present in Gaia DR2 detect variability at the Vi = 2 level, ∼ 52% are known UXORs although they are yet uncatalogued as HAeBes. In a sense this is (and just two have P-Cygni profiles). Moreover, this implies that a situation similar to that outlined for the low infrared excesses we are retrieving ∼ 81% of known A-type UXORs with our vari- observed toward the Herbig Be stars that are mostly located close ability indicator and hence, assuming that all the 25 A type stars to the Main Sequence. This could be explained by the fact that with variabilities above Vi = 2 are of UXOR type, there should the objects would be embedded and thus optically invisible or be about 31/85 UXORs in the sample. In turn this would imply faint in earlier phases of their evolution. that ∼ 37% of all Herbig Ae stars belong to the UXOR class. Further observations of optically fainter objects will be nec- If we also take into account that we have potentially removed essary to settle this issue. Additional progress can be made by some UXORs, possibly the most variable ones, from considera- connecting the PMS evolutionary tracks with radiative transfer tion when applying the constraints described in Sec. 3.4 we get codes to provide synthetic observations (as e.g. Davies et al. to values close to the 50% predicted by Natta & Whitney (2000). 2011, or Zhang et al. 2014 for Massive Young Stellar Objects) Finally, Davies et al. (2018) recently studied in detail the extended to optical wavelengths in the Herbig Be mass range. UXOR object CO Ori, which has single-peaked Hα emission. Related to the “missing” high mass stars in the HR diagram, it Consequently, they found that the inclination of its disk is of will be important to fill the historic, and entirely man-made, gap ∼ 30◦ (i.e. it is nearly face-on). In this particular case, whether if between the Herbig Ae stars and the T-Tauri stars. The latter the disk is still causing the UXOR phenomenon or if it is caused are confined to have spectral types G-K-M, and typically Herbig through fluctuations in the circumstellar material outside the disk Ae/Be stars, in this case by definition, have spectral types A and is still uncertain. We could not derive a variability indicator value B. We are missing out the F-type stars, resulting in an incomplete for this object to assess its variability. Inspired by this example, coverage of the HR diagram for pre-Main Sequence stars. we took a look to the other UXORs in our sample with single- peaked profiles, they all have variabilities below Vi = 2 in our variability indicator (HD 100546, HD 142527, HD 98922 and IL 5.6. On the difference between Herbig Ae and Herbig Be Cep), this suggesting a category of low variability UXORs with stars nearly face-on disks. Nonetheless, the results presented in this From the above it appears that the dusty disks surrounding Her- section strongly support the idea that most UXORs are caused by big Ae and Herbig Be stars are different, with the break in IR ex- edge-on disks, which are responsible of large photometric vari- cess occurring at 7M (around B3 spectral type), a value which abilities. was also found by Alonso-Albi et al. (2009) from their compi- lation of millimetre emission tracing the outer parts of the dusty 5.5. Missing objects in the HR diagram disks. As discussed, given the much stronger radiation field from B-stars, both in intensity and photon-energies, the most straight- When inspecting the right panel of Fig. 2, it appears that most forward explanation for the much less massive disks of higher Herbig Be stars are located relatively close to the Main Se- mass objects is a more efficient disk dispersal mechanism (see quence, whereas the lower mass Herbig Ae stars occupy a larger e.g. Gorti et al. 2009). This also explains why the same 7M part of their evolutionary tracks, contracting to higher tempera- break is seen in variability (Fig. 8). As described in Sect. 5.4, the tures at constant luminosity. In other words, the late type Herbig high levels of variability in some sources are caused by edge-on Be and Herbig Ae stars at high luminosities (and low surface dusty disks. A more efficient disk dispersal mechanism beyond gravities) that would occupy the tracks towards the locations of 7M would result in these sources showing no strong variability B-type stars on the Main Sequence are missing. It is only due in our indicator. It also explains why the objects with the lower to the use of Gaia parallaxes, expanding the number of Herbig IR excesses are not strongly variable while the rest can have both Ae/Be stars with well established luminosities that we can make high and low variability values. this observation. Other studies of large samples of Herbig Ae/Be stars indi- In our discussion earlier, we mentioned the fact that these cate a break in properties at a much lower mass of 3M , around objects could still be heavily embedded in their parental clouds, the B7 spectral boundary. Fairlamb et al. (2015) studied the ac- preventing them from being optically visible when evolving on cretion rates, which are proportional to the mass of the objects, their way to the Main Sequence. There is evidence for optically and found a different slope for lower mass than for higher mass invisible, but infrared bright objects at locations in these regions objects. Ababakr et al. (2017), extending the work of Mottram of the HR-diagram. For example, Pomohaci et al. (2017) were et al. (2007), found a distinct difference in spectro-polarimetric the first to spectrally type an infrared bright Massive Young properties across the Hα line between the Herbig Ae and late Be Stellar Object based on the, rare, absorption spectrum at near- type stars on the one hand and earlier Herbig Be type objects infrared wavelengths (higher order Brackett lines are in absorp- on the other hand. These authors also point out the similarity in tion for this object, while Brγ is in emission). They found that the the Hα spectro-polarimetry of the Herbig Ae stars and T-Tauri object could be fitted with that of an A-type giant star. Had this stars. Finally, Mendigutía et al. (2011a) noted the difference in object been optically visible, it would have occupied the empty Hα variability; Herbig Ae and late Be stars are largely variable, region in the right panel of Fig.2. To this, we add the early B- whereas Herbig Be stars are not. Later, Fang et al. (2013) showed type Herbig Be stars/infrared bright MYSOs PDS 27 and PDS that T-Tauri stars display even more variable Hα emission - again 37 (Ababakr et al. 2015). They are found in the upper regions of hinting at a similar accretion mechanism for the T-Tauri stars and the HR diagram, slightly off the Main Sequence. They are opti- Herbig Ae stars. cally visible, but not overly bright at V ∼ 13mag, and have not How can we reconcile the fact that some studies show a dif- been included in many (optical) magnitude limited catalogues. ferent break in properties than others? It is worth pointing out Thus, there are several examples that might lead us to conclude that the latter investigations consider regions much closer to the that the - implicit - optical brightness limit of any catalogue of star than the dusty emission. Fairlamb et al. (2015) derives accre- Herbig Ae/Be stars would prevent the inclusion of massive pre- tion rates from the UV-excess, which trace the shocked material Main Sequence stars on the horizontal portions of the evolution- on the stellar surface, Ababakr et al. (2017)’s spectropolarimetry

Article number, page 13 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 traces the free electrons in ionized material at distances of order another mechanism, possibly boundary layer accretion, for stellar radii from the stars. The spectro-polarimetric properties of the B-type objects. The differing IR and variability properties the B-type stars can be explained by stable circumstellar disks, are related to different or differently acting (dust-)disk dis- while the line properties for T-Tauri and Herbig Ae objects are persal mechanisms, which occurs at much larger size scales consistent with magnetically controlled accretion. Likewise, the than the accretion traced by hydrogen recombination line Hα emission traces the ionized zones close to the star, such as emission. the accretion columns and circumstellar disks and the variabil- ity is explained due to the accretion columns orbiting the central Finally, the findings presented in this paper signal just the star (e.g Kurosawa et al. 2008). beginning in unveiling the formation of intermediate mass stars Earlier, we showed that the IR fluxes and Hα properties are using Gaia. Gaia presents us with an excellent opportunity to largely correlated, but that the IR fluxes are smaller for the earlier search and identify new Herbig Ae/Be stars, resulting in a well- type objects. We therefore conclude this section with the obser- selected and properly characterized sample. The results pre- vation that the infrared and millimetre emission trace the circum- sented here will assist greatly in identifying new Herbig Ae/Be stellar disks and originates much further from the stars than the objects from the more than a billion stars with astrometric pa- UV, hydrogen recombination emission and free electrons which rameters in Gaia. This is the subject of our follow-on study, the trace the accretion onto the stars. The break in accretion mech- STARRY project. anism appears to occur around 3M , whereas the disk dispersal Acknowledgements. The STARRY project has received funding from the Eu- becomes significant at higher masses, 7M . ropean Union’s Horizon 2020 research and innovation programme under MSCA ITN-EID grant agreement No 676036. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www. 6. Conclusions cosmos.esa.int/gaia), processed by the Gaia Data Processing and Anal- ysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/ In this paper we have collated the largest astrometric dataset of consortium). Funding for the DPAC has been provided by national institutions, Herbig Ae/Be stars. We present parallaxes for the vast majority in particular the institutions participating in the Gaia Multilateral Agreement. This research also made use of Astropy, a community-developed core Python of known Herbig Ae/Be stars and have gathered atmospheric pa- package for Astronomy (Astropy Collaboration, 2013), the TOPCAT tool (Taylor rameters, optical and infrared photometry, extinction values, Hα 2005) and the VizieR catalogue access tool and the SIMBAD database, operated emission line information, binary statistics, and devised an ob- at CDS, Strasbourg, France. This research was made possible through the use of jective measure for the photometric variability. From these we the AAVSO Photometric All-Sky Survey (APASS), funded by the Robert Mar- tin Ayers Sciences Fund. We thank the referee for his/her insightful comments derive luminosities which allow us to place the objects in an HR which have improved the paper. diagram, containing over ten times more objects than previously possible.

Thus, we homogeneously derived luminosities, distances, References masses, ages, variabilities and infrared excesses for the most Aarnio, A. N., Monnier, J. D., Harries, T. J., et al. 2017, ApJ, 848, 18 complete sample of Herbig Ae/Be stars to date. We investigated Ababakr, K. M., Fairlamb, J. R., Oudmaijer, R. D., & van den Ancker, M. E. the various properties and reach the following conclusions: 2015, MNRAS, 452, 2566 Ababakr, K. M., Oudmaijer, R. D., & Vink, J. S. 2016, MNRAS, 461, 3089 Ababakr, K. M., Oudmaijer, R. D., & Vink, J. S. 2017, MNRAS, 472, 854 1. The Gaia photometric variability indicator as developed here Acke, B., van den Ancker, M. E. & Dullemond, C. P. 2005, A&A, 436, 209 indicates that 48/193 or ∼25% of all Herbig Ae/Be stars are Alecian, E., Wade, G. A., Catala, C., et al. 2013, MNRAS, 429, 1001 strongly variable. We find that the presence of variability cor- Alonso-Albi, T., Fuente, A., Bachiller, R., et al. 2009, A&A, 497, 117 relates very well with the Hα line profile. The variable ob- van den Ancker, M. E., The, P. S., Feinstein, A., et al. 1997, A&AS, 123, 63 van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 1998, A&A, 330, jects display doubly peaked profiles, indicating an edge-on 145 disk. It had been suggested that this variability is in most Arenou, F., Luri, X., Babusiaux, C., et al. 2018, submitted to A&A. cases due to asymmetric dusty disk structures seen edge-on. Aspin, C. 1998, A&A, 335, 1040 The observation here is the most compelling confirmation Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., et al. 2018, submitted to ApJ. Bailer-Jones, C. A. L. 2015, PASP, 127, 994 of this hypothesis. Most sources catalogued as UXORs in Baines, D., Oudmaijer, R. D., Porter, J. M., & Pozzo, M. 2006, MNRAS, 367, the sample appear as strongly variable with double-peaked 737 profiles. The fraction of strongly variable A-type objects is Banzatti, A., Garufi, A., Kama, M., et al. 2018, A&A, 609, L2 close to that found for the A-type objects with the UXOR Biller, B., Lacour, S., Juhász, A., et al. 2012, ApJ, 753, L38 phenomenon. Boehm, T., & Catala, C. 1995, A&A, 301, 155 Boersma, C., Peeters, E., Martín-Hernández, N. L., et al. 2009, A&A, 502, 175 2. High mass stars do not display an infrared excess and show Borges Fernandes, M., Kraus, M., Lorenz Martins, S., & de Araújo, F. X. 2007, no strong photometric variability. Several suggestions have MNRAS, 377, 1343 been put forward to explain this. These include fast evo- Bouvier, J., Alencar, S. H. P., Harries, T. J., Johns-Krull, C. M., & Romanova, lutionary timescales and fast dust dispersion timescales for M. M. 2007, Protostars and Planets V, 479 Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 high mass objects. We do note that the break is around 7 M , Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 which is intriguingly similar to other statistical studies re- Carmona, A., van den Ancker, M. E., Audard, M., et al. 2010, A&A, 517, A67 lated to dusty disks around Herbig Ae/Be stars which sign- Castelli, F., & Kurucz, R. L. 2004, arXiv:astro-ph/0405087 post a different or more efficient disk dispersal mechanism Cauley, P. W., & Johns-Krull, C. M. 2015, ApJ, 810, 5 Chelli, A., Cruz-Gonzalez, I., & Reipurth, B. 1995, A&AS, 114, 135 for high mass objects. Chen, P. S., Shan, H. G., & Zhang, P. 2016, New A, 44, 1 3. Whereas the break in IR properties and photometric vari- Clem, J. L., & Landolt, A. U. 2016, AJ, 152, 91 abilities occurs at 7 M , various Hα line properties includ- Connelley, M. S., Reipurth, B. & Tokunaga, A. T. 2008, AJ, 135, 2496. ing mass accretion rates, spectropolarimetric properties and Corporon, P., & Lagrange, A.-M. 1999, A&AS, 136, 429 Costigan, G., Vink, J. S., Scholz, A., Ray, T., & Testi, L. 2014, MNRAS, 440, emission line variability seem to differ at a lower mass, 3 M . 3444 The latter has been linked to different accretion mechanisms Coulson, I. M., & Walther, D. M. 1995, MNRAS, 274, 977 at work; magnetospheric accretion for the A-type objects and Coté, J., & Waters, L. B. F. M. 1987, A&A, 176, 93

Article number, page 14 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars

Cutri, R. M., Wright, E. L., Conrow, T., et al. 2013, Explanatory Supplement to Marigo, P., Girardi, L., Bressan, A., et al. 2017, ApJ, 835, 77 the AllWISE Data Release Products, by R. M. Cutri et al. , Marconi, M., & Palla, F. 1998, ApJ, 507, L141 Davies, B., Oudmaijer, R. D., & Sahu, K. C. 2007, ApJ, 671, 2059 Marston, A. P. & McCollum, B. 2008, A&A, 477, 193 Davies, B., Hoare, M. G., Lumsden, S. L., et al. 2011, MNRAS, 416, 972 Mayer, A., Deschamps, R., & Jorissen, A. 2016, A&A, 587, A30 Davies, C. L., Kreplin, A., Kluska, J., et al. 2018, MNRAS, 474, 5406 Meeus, G., Waters, L. B. F. M., Bouwman, J., et al. 2001, A&A, 365, 476 Deason, A. J., Belokurov, V., Erkal, D., Koposov, S. E., & Mackey, D. 2017, Mendigutía, I., Eiroa, C., Montesinos, B., et al. 2011a, A&A, 529, A34 MNRAS, 467, 2636 Mendigutía, I., Calvet, N., Montesinos, B., et al. 2011b, A&A, 535, A99 Doering, R. L., & Meixner, M. 2009, AJ, 138, 780 Mendigutía, I. 2011c, Ph.D. thesis, Universidad Autónoma de Madrid Duchêne, G. 2015, Ap&SS, 355, 291 Mendigutía, I., Mora, A., Montesinos, B., et al. 2012, A&A, 543, A59 Dullemond, C. P., & Dominik, C. 2004a, A&A, 417, 159 Millour, F., Chesneau, O., Borges Fernandes, M., et al. 2009, A&A, 507, 317 Dullemond, C. P., & Dominik, C. 2004b, A&A, 421, 1075 Miroshnichenko, A. S., Fremat, Y., Houziaux, L., et al. 1998, Astronomy and Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971 Astrophysics Supplement Series, 131, 469 Dunhill, A. C., Cuadra, J., & Dougados, C. 2015, MNRAS, 448, 3545 Miroshnichenko, A. S., Gray, R. O., Vieira, S. L. A., Kuratov, K. S., & Bergner, Dunkin, S. K., Barlow, M. J., & Ryan, S. G. 1997, MNRAS, 290, 165 Y. K. 1999, A&A, 347, 137 Eiroa, C., Oudmaijer, R. D., Davies, J. K., et al. 2002, A&A, 384, 1038 Miroshnichenko, A. S., Chentsov, E. L., Klochkova, V. G., et al. 2000, Astron- Evans, D. W., Riello, M., De Angeli, F., et al. 2018, ArXiv e-prints , omy and Astrophysics Supplement Series, 147, 5 arXiv:1804.09368 Miroshnichenko, A. S., Bjorkman, K. S., Chentsov, E. L., et al. 2002, A&A, 388, Fairlamb, J. R., Oudmaijer, R. D., Mendigutía, I., Ilee, J. D., & van den Ancker, 563 M. E. 2015, MNRAS, 453, 976 Miroshnichenko, A. S., Levato, H., Bjorkman, K. S., et al. 2004, A&A, 417, 731 Fairlamb, J. R., Oudmaijer, R. D., Mendigutia, I., Ilee, J. D., & van den Ancker, Miroshnichenko, A. S., Manset, N., Zharikov, S. V., et al. 2014, Advances in M. E. 2017, MNRAS, 464, 4721 Astronomy, 2014, E7 Fang, M., Kim, J. S., van Boekel, R., et al. 2013, ApJS, 207, 5 Monnier, J. D., Millan-Gabet, R., Billmeier, R., et al. 2005, ApJ, 624, 832 Arellano Ferro, A., & Giridhar, S. 2003, A&A, 408, L29 Montesinos, B., Eiroa, C., Mora, A., & Merín, B. 2009, A&A, 495, 901 Finkenzeller, U., & Mundt, R. 1984, A&AS, 55, 109 Morrell, N., & Levato, H. 1991, ApJS, 75, 965 Folsom, C. P., Wade, G. A., Kochukhov, O., et al. 2008, Contributions of the Mottram, J. C., Vink, J. S., Oudmaijer, R. D., & Patel, M. 2007, MNRAS, 377, Astronomical Observatory Skalnate Pleso, 38, 245 Frasca, A., Miroshnichenko, A. S., Rossi, C., et al. 2016, A&A, 585, A60 1363 Friedemann, C., Guertler, J., & Loewe, M. 1996, A&AS, 117, 205 Nakano, M., Sugitani, K., Watanabe, M., et al. 2012, AJ, 143, 61 Fuente, A., Martin-Pintado, J., Bachiller, R., Neri, R., & Palla, F. 1998, A&A, Natta, A., Grinin, V. P., Mannings, V., & Ungerechts, H. 1997, ApJ, 491, 885 334, 253 Natta, A., & Whitney, B. A. 2000, A&A, 364, 633 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1. Oudmaijer, R. D., & Drew, J. E. 1997, A&A, 318, 198 Gaia Collaboration, Babusiaux, C., van Leeuwen, F., et al. 2018, accepted for Oudmaijer, R. D., Busfield, G., & Drew, J. E. 1997, MNRAS, 291, 797 publication by A&A. Oudmaijer, R. D., Proga, D., Drew, J. E., et al. 1998, MNRAS, 300, 170 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, ArXiv e-prints , Oudmaijer, R. D., & Drew, J. E. 1999, MNRAS, 305, 166 arXiv:1804.09365. Oudmaijer, R. D., Palacios, J., Eiroa, C., et al. 2001, A&A, 379, 564 Gorti, U., Dullemond, C. P., & Hollenbach, D. 2009, ApJ, 705, 1237 Oudmaijer, R. D. 2017, The B[e] Phenomenon: Forty Years of Studies, 508, 175 Grady, C. A., Hamaguchi, K., Schneider, G., et al. 2010, ApJ, 719, 1565 Patel, P., Sigut, T. A. A., & Landstreet, J. D. 2017, ApJ, 836, 214 Gray, R. O., & Corbally, C., J. 2009, Stellar Spectral Classification by Richard Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9 O. Gray and Christopher J. Corbally. Princeton University Press, 2009. ISBN: Pogodin, M. A., Hubrig, S., Yudin, R. V., et al. 2012, Astronomische 978-0-691-12511-4 Nachrichten, 333, 594 Grinin, V. P. & Rostopchina, A. N. 1996, Astronomy Reports, 40, 171 Polster, J., Korcáková,ˇ D., Votruba, V., et al. 2012, A&A, 542, A57 Grinin, V. P. 2000, Disks, Planetesimals, and Planets, 219, 216 Pomohaci, R., Oudmaijer, R. D., Lumsden, S. L., Hoare, M. G., & Mendigutía, Grundstrom, E. D., & Gies, D. R. 2006, ApJ, 651, L53 I. 2017, MNRAS, 472, 3624 Gürtler, J., Friedemann, C., Reimann, H.-G., et al. 1999, Astronomy and Astro- Poxon, M. 2015, Journal of the American Association of Variable Star Observers physics Supplement Series, 140, 293 (JAAVSO), 43, 35 Hamaguchi, K., Choi, M., Corcoran, M. F., et al. 2008, ApJ, 687, 425-432 Reiter, M., Calvet, N., Thanathibodee, T., et al. 2018, ApJ, 852, 5 Herbig, G. H. 1960, ApJS, 4, 337 Rivinius, T., Carciofi, A. C., & Martayan, C. 2013, A&A Rev., 21, 69 Herbig, G. H. & Bell, K. R. 1988, Third catalog of emission-line stars of the Sartori, M. J., Gregorio-Hetem, J., Rodrigues, C. V., Hetem, A., Jr., & Batalha, Orion population. C. 2010, AJ, 139, 27-38 Hernández, J., Calvet, N., Briceño, C., Hartmann, L., & Berlind, P. 2004, AJ, Schöller, M., Pogodin, M. A., Cahuasquí, J. A., et al. 2016, A&A, 592, A50 127, 1682 Schütz, O., Meeus, G., Carmona, A., Juhász, A., & Sterzik, M. F. 2011, A&A, Hernández, J., Calvet, N., Hartmann, L., et al. 2005, AJ, 129, 856 533, A54 Holl, B., Audard, M., Nienartowicz, K., et al. 2018, accepted for publication by Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 A&A. Spezzi, L., Alcalá, J. M., Covino, E., et al. 2008, ApJ, 680, 1295 Honda, M., Maaskant, K., Okamoto, Y. K., et al. 2015, ApJ, 804, 143 Tambovtseva, L. V., Grinin, V. P. & Weigelt, G. 2014, A&A, 562, A104 Hou, W., Luo, A.-L., Hu, J.-Y., et al. 2016, Research in Astronomy and Astro- Taylor, M. B. 2005, Astronomical Data Analysis Software and Systems XIV, physics, 16, 138 347, 29 Ilee, J. D., Fairlamb, J., Oudmaijer, R. D., et al. 2014, MNRAS, 445, 3723 Testi, L., Palla, F., & Natta, A. 1999, A&A, 342, 515 Joner, M. D., & Hintz, E. G. 2015, AJ, 150, 204 Thé, P. S., de Winter, D., & Perez, M. R. 1994, A&AS, 104, 315 Klement, R., Carciofi, A. C., Rivinius, T., et al. 2017, A&A, 601, A74 Torres, C. A. O., Quast, G. R., da Silva, L., & de La Reza, R. 2000, IAU Sympo- Kraus, S., Calvet, N., Hartmann, L., et al. 2012, ApJ, 746, L2 sium, 200, 118 Kraus, S. 2015, Ap&SS, 357, 97 Vieira, S. L. A., Corradi, W. J. B., Alencar, S. H. P., et al. 2003, AJ, 126, 2971 Kubát, J., Saad, S. M., Kawka, A., et al. 2010, A&A, 520, A103 Vieira, R. G., Gregorio-Hetem, J., Hetem, A., et al. 2011, A&A, 526, A24 Kucerová,ˇ B., Korcáková,ˇ D., Polster, J., et al. 2013, A&A, 554, A143 Vink, J. S., Drew, J. E., Harries, T. J., & Oudmaijer, R. D. 2002, MNRAS, 337, Kurosawa, R., Harries, T. J. & Symington, N. H. 2006, MNRAS, 370, 580 356 Kurosawa, R., Romanova, M. M., & Harries, T. J. 2008, MNRAS, 385, 1931 Waters, L. B. F. M., Cote, J., & Lamers, H. J. G. L. M. 1987, A&A, 185, 206 Lazareff, B., Berger, J.-P., Kluska, J., et al. 2017, A&A, 599, A85 Watson, C. L., Henden, A. A., & Price, A. 2006, Society for Astronomical Sci- Landolt, A. U. 2009, AJ, 137, 4186 ences Annual Symposium, 25, 47 Lee, C.-D., Chen, W.-P. & Liu, S.-Y. 2016, A&A, 592, A130 Wheelwright, H. E., Oudmaijer, R. D., & Goodwin, S. P. 2010, MNRAS, 401, Leinert, C., Richichi, A., & Haas, M. 1997, A&A, 318, 472 Lindegren, L., Hernandez, J., Bombrun, A., et al. 2018, accepted for publication 1199 by A&A. Wheelwright, H. E., Vink, J. S., Oudmaijer, R. D., & Drew, J. E. 2011, A&A, Liu, Q. Z., van Paradijs, J., & van den Heuvel, E. P. J. 2000, A&AS, 147, 25 532, A28 Lumsden, S. L., Hoare, M. G., Urquhart, J. S., et al. 2013, ApJS, 208, 11 Zhang, P., Yang, H. T., & Liu, J. 2006, Ap&SS, 305, 11 Luri, X., Brown, A. G. A., Sarro, L. M., et al. 2018, accepted for publication by Zhang, Y., Tan, J. C., & Hosokawa, T. 2014, ApJ, 788, 166 A&A. Zhang, Q., Claus, B., Watson, L., et al. 2017, ApJ, 837, 53 Maaskant, K. M., Honda, M., Waters, L. B. F. M., et al. 2013, A&A, 555, A64 Zuckerman, B., Melis, C., Song, I., et al. 2008, ApJ, 683, 1085 Maheswar, G., Manoj, P., & Bhatt, H. C. 2002, A&A, 387, 1003 Malkov, O. Y., Oblak, E., Snegireva, E. A., & Torra, J. 2006, A&A, 446, 785 Manoj, P., Bhatt, H. C., Maheswar, G., & Muneer, S. 2006, ApJ, 653, 657

Article number, page 15 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 1 6 4 4 5 7 5 8 V Binary 7.32 Yes 8.90 Yes 9.90 - 9.82 - 8.98 Yes 7.78 Yes 6.69 - 9.24 Yes 6.52 Yes 7.23 - 8.91 - 12.49 - 12.44 - 10.05 Yes 10.51 - 11.23 - 11.52 - 11.13 Yes 14.94 - 14.75 - 14.00 - 11.59 - 14.13 - 12.35 - 11.48 - 13.08 - 10.73 - 12.00 - 12.01 - 14.32 - 10.27 - 13.65 - 14.37 - 10.05 - V 28 35 26 29 62 42 46 12 24 12 50 50 24 35 20 12 20 19 29 23 42 28 13 54 29 23 35 17 34 12 29 20 25 29 29 23 20 06 20 11 23 15 19 07 29 23 29 23 26 29 ...... 00 7.65 - 050 000 030 020 040 000 050 000 050 000 030 000 050 000 050 000 050 000 . 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...... 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 0 + − + − + − + − + − + − + − + − + − 43 13 27 19 73 06 83 12 14 83 41 94 64 78 92 89 58 51 06 11 37 07 53 48 61 ...... 0 4 2 1 1 1 0 0 2 3 0 3 4 0 1 0 2 1 0 0 2 0 4 5 2 000 330 940 000 000 210 000 000 000 ...... 0 0 0 0 0 0 0 0 0 ) (mag) (mag) 19 21 39 44 47 27 36 12 25 14 31 31 10 14 10 09 51 39 18 12 33 45 49 38 23 13 26 11 13 09 18 18 46 35 12 06 12 08 35 19 42 32 42 33 12 09 34 40

...... 028 005 060 051 088 056 055 038 093 046 024 003 069 046 049 036 035 011 073 050 015 013 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 ...... 0 0 0 0 . . . . L 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 61 17 01 22 54 96 76 50 85 87 58 85 53 38 47 17 70 79 82 18 87 71 24 31 Log(L) A ...... 1 4 3 2 1 2 0 1 2 0 1 3 1 1 0 2 3 0 0 2 3 3 3 4 623 287 467 243 983 790 371 581 329 606 119 ...... 0 1 0 1 0 0 1 1 1 1 1 ff e 750 790 250 250 200 200 210 300 100 100 150 80 370 200 600 380 900 300 250 100 250 250 250 250 150 80 150 80 80 160 250 250 500 500 250 250 250 250 500 500 250 250 250 250 4500 5000 1600 800 1700 500 1500 800 1000 1000 5000 3000 3000 3000 5000 3000 5000 3000 3000 800 5000 3000 5000 3000 4500 5000 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 9500 6250 8970 6630 6650 6250 6750 8900 9800 6750 9750 7000 6250 6250 6170 7250 9750 9750 8000 6750 8200 16000 10250 24500 10000 11500 17000 19500 11000 19500 19500 16500 19500 19500 24500 6 4 1 0 4 3 7 5 7 5 3 1 0 0 1 6 1 1 9 1 ...... 30 27 13 11 64 53 14 12 12 11 52 43 46 39 82 68 12 11 35 30 31 27 10 9 70 69 57 44 15 13 2 2 2 2 3 3 4 4 3 3 5 5 . . 1 1 8 7 1 1 2 3 350 240 620 440 470 290 580 390 170 130 370 250 180 140 250 190 210 160 180 130 + − + − + − + − + − + − + − + − + − + − + − + − + − + − 0 0 + − + − + − + − + − + − + − + − + − − + + − + − + − + − + − + − + − + − + − + − + − 9 6 1 5 1 3 0 3 4 1 ...... 20 . 057 162 062 140 046 2110 027 4040 029 773 083 295 035 1003 053 389 018 335 020 374 042 404 066 2130 082 163 042 833 046 3000 029 942 023 1409 021 386 023 558 030 1780 037 2390 040 696 035 706 032 1770 039 386 019 2760 027 2070 039 104 051 110 028 411 058 108 062 705 036 1650 057 382 090 142 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 140 113 390 137 290 384 987 570 984 670 474 293 131 191 154 055 701 588 792 538 346 432 411 535 592 338 456 597 089 431 226 405 574 618 036 ...... 30:33:04 6 57:21:31 0 47:13:44 1 31:24:57 3 42:18:01 0 69:44:36 2 70:03:45 2 11:25:39 2 14:30:28 0 24:44:54 6 29:57:55 0 57:26:48 1 09:33:37 1 09:44:46 1 55:23:00 0 61:38:23 0 43:49:32 0 + + + + + + + + + + + + + + + + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( Be star belonging to the high quality sample of 218 sources. / 30 549 03:29:19.8 41 3731 20:24:15.7 + + NameAB Aur 04:55:45.9 RA DEC Parallax Distance T AK Sco 16:54:44.8 -36:53:19 7 AS 310 18:33:21.2 -04:58:06 0 AS 470 21:36:14.2 AS 477 21:52:34.1 BD BD BF Ori 05:37:13.3 -06:35:01 2 BH Cep 22:01:42.9 BO Cep 22:16:54.1 CO Ori 05:27:38.3 CPM 25 06:23:56.3 CQ Tau 05:35:58.5 DG Cir 15:03:23.8 -63:22:59 1 GSC 1876-0892 06:07:15.4 GSC 3975-0579 21:38:08.5 GSC 6546-3156 07:24:17.5 -26:16:05 0 GSC 8143-1225 07:59:11.6 -50:22:47 2 GSC 8581-2002 08:44:23.6 -59:56:58 1 GSC 8645-1401 12:17:47.5 -59:43:59 0 GSC 8994-3902 13:19:04.0 -62:34:10 0 HBC 217 06:40:42.2 HBC 222 06:40:51.2 HBC 334 02:16:30.7 HBC 442 05:34:14.2 -05:36:54 2 HBC 7 00:43:25.3 HBC 705 20:51:02.7 HD 100453 11:33:05.5 -54:19:29 9 HD 100546 11:33:25.3 -70:11:41 9 HD 101412 11:39:44.4 -60:10:28 2 HD 104237 12:00:04.9 -78:11:35 9 HD 114981 13:14:40.7 -38:39:06 1 HD 130437 14:50:50.2 -60:17:10 0 HD 132947 15:04:56.0 -63:07:53 2 HD 135344 15:15:48.9 -37:08:56 7 Main parameters of each Herbig Ae Table 2.

Article number, page 16 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 9 4 4 4 2 4 2 2 6 6 10 11 12 13 14 V Binary 8.63 Yes 8.40 - 7.10 Yes 8.64 Yes 8.27 Yes 8.67 - 8.17 Yes 6.75 - 8.80 Yes 8.64 Yes 4.81 - 6.85 - 8.16 - 4.24 - 8.70 - 7.56 Yes 7.38 Yes 5.96 Yes 7.33 Yes 9.82 - 9.38 - 9.91 Yes 9.54 - 8.72 Yes 9.68 Yes 9.32 - 9.88 - 9.27 - 10.10 - 11.04 - 10.10 - 10.69 - 10.40 - 10.02 Yes 10.86 - 10.01 - V 12 25 14 14 12 34 34 00 12 34 14 00 29 23 12 34 12 26 23 15 34 09 24 12 21 30 34 12 27 06 12 14 23 15 ...... 050 060 050 000 010 000 030 040 050 000 080 090 060 000 020 040 020 050 000 040 050 000 050 000 050 000 050 000 050 000 050 000 000 020 050 000 059 000 ...... 0 0 . . 0 0 . . 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 ...... 0 0 . . 0 0 . . 0 0 0 0 . . . . 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 31 95 47 00 02 00 50 38 53 05 09 20 11 57 06 16 80 ...... 0 0 0 0 1 0 2 0 0 1 0 0 1 0 0 0 1 230 000 010 400 000 500 000 550 000 000 100 140 000 000 000 000 000 000 300 ...... 0 0 0 2 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 ) (mag) (mag) 15 14 09 12 14 12 20 28 26 07 12 22 19 17 43 34 12 22 09 14 29 14 23 08 13 10 29 18 17 12 27 40 24 10 13 04 12 09 11 08 15 14 42 23

...... 033 035 032 010 032 027 026 005 043 045 038 011 039 044 055 032 036 012 058 063 057 066 053 031 093 091 085 060 ...... 0 0 . . . . 0 0 . . 0 0 . . 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 . . . . 0 0 0 0 L 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 74 46 09 74 22 31 58 23 58 05 07 51 29 56 94 97 66 84 27 81 55 47 Log(L) A ...... 4 0 2 2 2 1 2 4 1 2 3 1 1 1 1 2 1 0 1 0 4 3 786 773 216 963 939 971 367 199 391 153 461 116 178 522 ...... 0 0 1 0 0 0 1 1 1 1 1 1 1 1 ff 38 80 e 120 120 250 250 250 250 250 250 250 250 250 250 850 600 250 250 800 900 900 300 250 250 800 900 800 900 200 200 220 130 900 300 250 250 250 250 500 500 500 500 120 120 900 300 150 150 500 500 500 500 380 380 250 250 + − 1500 1500 1000 1500 5000 3000 3000 800 1500 800 2100 2900 1000 1000 3000 800 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 5430 6380 7750 9750 6500 7500 7500 9000 9800 9250 9500 7380 9800 8500 9000 8380 9800 6400 9750 9500 7880 10250 10700 10700 10700 10000 11000 10500 28000 13000 19500 16500 11500 14000 34000 16500 4 3 6 6 0 9 0 9 1 2 4 2 5 2 7 5 7 2 9 0 3 4 3 9 5 6 2 6 3 4 ...... 91 88 10 9 15 14 19 17 19 17 37 29 66 49 94 64 44 37 12 11 15 13 10 9 31 25 24 20 18 15 21 18 2 2 1 1 . . 2 1 2 1 7 6 2 2 3 3 2 2 6 8 1 2 1 1 6 7 2 2 5 5 1 1 150 110 100 70 130 90 630 440 280 170 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 0 0 + − + − + − + − + − + − + − + − − + − + − + − + − + − + − + + − + − + − + − + − + − 8 7 3 3 1 4 8 8 8 5 0 7 2 0 3 ...... 63 . 15 166 30 122 11 101 37 106 10 89 11 429 ...... 077 135 053 134 044 110 044 1340 047 157 053 148 058 155 080 158 070 150 058 950 064 114 059 1100 061 159 047 266 045 361 032 524 059 432 062 421 078 669 098 697 047 721 053 359 056 396 048 354 082 455 068 438 063 398 074 399 027 4040 077 1400 ...... 0 . . . . 0 0 . 0 . . . 0 . . . . 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 02 15 85 32 38 20 ...... 365 424 039 717 356 744 437 296 632 024 775 879 281 759 771 907 313 374 476 401 379 784 524 821 191 277 510 502 134 616 ...... 08:42:10 0 15:47:15 3 50:59:48 1 11:20:20 2 11:17:41 2 10:01:51 2 16:39:57 1 16:30:57 1 10:19:20 1 02:27:47 2 02:28:10 2 00:25:08 2 68:09:48 2 + + + + + + + + + + + + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( NameHD 135344B 15:15:48.4 -37:09:16 7 RA DEC Parallax Distance T HD 139614 15:40:46.4 -42:29:54 7 HD 141569 15:49:57.7 -03:55:17 9 HD 141926 15:54:21.8 -55:19:44 0 HD 142527 15:56:41.9 -42:19:24 6 HD 142666 15:56:40.0 -22:01:40 6 HD 143006 15:58:36.9 -22:57:16 6 HD 144432 16:06:57.9 -27:43:10 6 HD 149914 16:38:28.6 -18:13:14 6 HD 150193 16:40:17.9 -23:53:45 6 HD 155448 17:12:58.8 -32:14:34 1 HD 158643 17:31:25.0 -23:57:46 8 HD 163296 17:56:21.3 -21:57:22 9 HD 169142 18:24:29.8 -29:46:50 8 HD 17081 02:44:07.3 -13:51:32 9 HD 174571 18:50:47.2 HD 176386 19:01:38.9 -36:53:27 6 HD 179218 19:11:11.3 HD 199603 20:58:41.8 -14:29:00 11 HD 200775 21:01:36.9 HD 235495 21:21:27.5 HD 244314 05:30:19.0 HD 244604 05:31:57.3 HD 245185 05:35:09.6 HD 249879 05:58:55.8 HD 250550 06:02:00.0 HD 259431 06:33:05.2 HD 287823 05:24:08.0 HD 288012 05:33:04.8 HD 290380 05:23:31.0 -01:04:24 2 HD 290409 05:27:05.5 HD 290500 05:29:48.1 -00:23:43 2 HD 290764 05:38:05.3 -01:15:22 2 HD 290770 05:37:02.4 -01:37:21 2 HD 305298 10:33:05.0 -60:19:51 0 HD 313571 18:01:07.2 -22:15:04 0 continued. Table 2.

Article number, page 17 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 2 4 2 5 5 2 5 1 5 1 4 1 15 16 17 18 19 30 V Binary 7.62 - 9.89 Yes 9.09 Yes 8.17 Yes 8.12 Yes 8.27 Yes 6.09 Yes 8.03 - 8.45 - 9.67 Yes 8.84 Yes 7.90 - 4.59 Yes 7.94 Yes 8.38 Yes 9.07 - 4.36 - 4.93 Yes 8.50 Yes 8.64 - 6.91 - 6.67 Yes 8.42 - 6.85 Yes 8.31 - 9.82 - 8.07 Yes 8.51 Yes 9.26 - 8.51 Yes 9.12 - 8.19 - 8.47 - 10.99 - 11.11 - 10.65 - V 28 06 15 17 17 16 25 29 21 00 50 50 25 15 12 24 24 34 12 16 12 34 23 15 19 13 24 18 50 21 25 00 35 17 23 15 25 22 23 15 12 03 25 08 19 12 50 50 ...... 020 000 050 000 050 040 050 000 050 000 050 000 020 010 030 020 050 000 050 000 050 000 030 010 0 0 0 0 0 0 . . 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...... 0 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 06 38 32 39 00 81 15 38 52 82 83 17 13 46 21 00 51 57 79 68 03 08 37 90 ...... 0 2 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 010 000 060 000 000 000 620 890 000 000 000 780 ...... 0 0 0 0 0 0 0 0 0 0 0 0 ) (mag) (mag) 14 05 24 21 20 21 11 13 10 02 30 24 12 08 18 22 24 29 14 15 12 10 35 20 15 24 46 25 19 14 30 22 31 20 12 02 24 14 32 17 33 38 34 18 13 09 11 04 11 09 11 08 10 08 16 11 64 53 19 15 10 07 20 17

...... 049 039 046 023 060 036 087 074 0 0 0 0 0 0 . . 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 L 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 27 19 51 72 36 45 04 13 43 83 24 04 45 69 17 19 72 42 83 88 89 04 46 97 44 51 19 84 60 99 85 81 Log(L) A ...... 1 3 2 0 1 1 1 3 2 2 1 2 2 3 2 2 1 1 2 2 3 4 2 0 2 2 3 2 4 2 2 4 980 789 193 549 . . . . 0 1 1 3 ff e 200 200 750 750 250 250 120 250 110 100 900 300 250 250 220 300 830 830 500 500 250 250 800 900 800 900 160 220 600 380 450 450 250 250 200 200 500 500 250 250 500 500 500 500 250 250 250 250 800 1000 1100 1300 1000 1000 3000 800 1000 700 2900 800 3000 800 3400 4200 3000 800 5000 3000 1000 1000 1500 1500 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 8250 9500 7000 5900 9800 7000 7600 9750 9500 7830 8900 9450 7000 9250 15750 11930 10700 10700 10500 12500 19000 13000 10000 10000 15000 11200 20000 16500 10500 13600 16500 21000 16500 19500 11500 25500 4 2 9 2 6 2 6 2 9 8 9 0 1 2 8 4 ...... 11 10 13 12 35 27 29 24 19 16 24 19 18 16 80 43 19 16 22 18 24 20 16 13 41 33 16 14 39 33 10 9 55 46 53 45 82 66 3 3 7 7 5 5 4 4 2 2 0 1 2 2 6 7 200 130 120 90 280 120 100 70 100 80 280 190 570 350 210 150 510 350 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − − + − + − + + − + − + − + − + − + − + − + − + − 8 5 5 0 3 1 3 5 ...... 10 163 30 650 36 312 10 338 30 44 17 222 12 2010 ...... 076 161 063 1300 055 1070 047 311 036 300 048 356 052 387 067 160 098 435 072 474 064 408 099 363 062 411 058 428 081 405 058 620 058 930 046 1090 064 380 049 165 043 318 049 668 037 396 037 911 028 1023 043 1910 034 1830 033 1168 034 3100 ...... 0 ...... 0 . 0 . 0 . 0 0 ...... 0 . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 13 27 16 95 68 01 50 ...... 182 716 908 210 328 805 580 240 288 103 447 751 428 335 463 603 049 897 630 051 139 489 526 089 971 465 511 844 217 ...... 29:50:37 6 21:31:51 3 05:38:43 2 24:57:37 6 25:19:57 6 17:03:29 2 04:07:17 3 05:05:04 0 + + + + + + + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( NameHD 31648 04:58:46.3 RA DEC Parallax Distance T HD 319896 17:31:05.9 -35:08:29 0 HD 323771 17:34:04.6 -39:23:41 0 HD 34282 05:16:00.5 -09:48:35 3 HD 344261 19:21:53.5 HD 34700 05:19:41.4 HD 35187 05:24:01.2 HD 35929 05:27:42.8 -08:19:39 2 HD 36112 05:30:27.5 HD 36408 05:32:14.1 HD 36917 05:34:47.0 -05:34:15 2 HD 36982 05:35:09.8 -05:27:53 2 HD 37258 05:36:59.3 -06:09:16 2 HD 37357 05:37:47.1 -06:42:30 1 HD 37371 05:38:09.9 -00:11:01 2 HD 37490 05:39:11.1 HD 37806 05:41:02.3 -02:43:01 2 HD 38087 05:43:00.6 -02:18:45 2 HD 38120 05:43:11.9 -04:59:50 2 HD 39014 05:44:46.3 -65:44:08 22 HD 41511 06:04:59.1 -16:29:04 4 HD 45677 06:28:17.4 -13:03:11 1 HD 46060 06:30:49.8 -09:39:15 1 HD 50083 06:51:45.8 HD 50138 06:51:33.4 -06:57:59 2 HD 56895B 07:18:31.8 -11:11:34 6 HD 58647 07:25:56.1 -14:10:44 3 HD 59319 07:28:36.8 -21:57:49 1 HD 68695 08:11:44.6 -44:05:09 2 HD 76534 08:55:08.7 -43:28:00 1 HD 85567 09:50:28.5 -60:58:03 0 HD 87403 10:02:51.4 -59:16:55 0 HD 87643 10:04:30.3 -58:39:52 0 HD 94509 10:53:27.2 -58:25:24 0 HD 95881 11:01:57.6 -71:30:48 0 HD 96042 11:03:40.5 -59:25:59 0 continued. Table 2.

Article number, page 18 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 1 3 2 3 3 2 1 2 32 V Binary 5.61 - 8.44 - 6.77 Yes 6.78 Yes 9.21 Yes 9.39 - 9.78 - 11.87 - 11.55 - 13.00 - 11.02 - 13.76 - 10.35 - 10.61 - 13.40 - 10.42 - 12.36 Yes 14.16 - 14.63 - 11.44 Yes 10.60 Yes 13.29 - 15.20 - 12.63 - 13.47 - 12.31 - 12.85 - 11.61 Yes 12.10 - 12.03 Yes 10.56 - 10.03 - 12.60 Yes 10.62 - 10.84 - 14.43 - 3 3 V . . 18 00 18 26 18 26 23 15 14 18 23 15 22 37 10 15 16 18 12 34 34 36 16 18 17 16 14 24 24 12 34 12 60 34 18 26 18 26 23 33 15 17 26 29 60 34 26 29 29 23 18 26 29 13 . . 2 2 ...... 050 020 010 000 050 040 050 020 020 010 010 030 010 010 040 030 0 0 ...... 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 1 + − + − + − + − + − + − + − + − . 4 00 48 84 18 24 16 52 70 20 60 36 02 65 05 64 64 89 03 63 25 31 94 03 06 62 88 13 ...... 0 3 3 1 5 3 0 2 3 2 0 2 2 4 3 0 8 5 4 4 2 0 5 3 3 5 0 900 090 330 230 800 570 540 470 ...... 0 0 0 0 0 0 3 8 ) (mag) (mag) 3 3 . . 12 10 10 08 44 43 34 37 39 21 32 22 25 29 29 15 13 18 11 11 24 23 17 26 35 28 18 19 13 15 09 13 23 13 32 16 77 87 37 37 37 38 12 12 35 34 25 21 36 40 72 83 38 44 48 37 28 33 15 08

. . . . 1 1 ...... 093 020 072 060 059 052 048 041 058 055 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − L 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 9 + − + − + − + − + − . 3 47 06 12 49 93 07 03 86 89 71 09 13 04 57 91 16 08 10 49 52 94 59 89 58 12 62 32 21 62 19 Log(L) A ...... 1 2 4 3 2 2 5 3 0 0 2 1 1 2 1 1 2 1 5 4 4 4 4 3 4 4 4 4 3 1 174 544 032 717 922 . . . . . 1 1 3 1 1 ff e 200 200 500 500 250 250 250 250 500 500 250 250 500 500 340 380 500 500 800 900 500 580 220 130 250 250 900 300 750 750 250 200 5100 4800 3900 4500 3900 4500 3000 800 2000 3900 3000 800 1000 1000 1000 1000 800 1000 1500 800 9200 9200 3900 4500 3900 4500 1500 1500 5900 4000 4500 5000 9200 9200 4500 5000 5000 3000 3900 4500 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 8900 8500 7990 8500 7830 7380 9800 6750 10500 10500 10500 13000 14000 10700 10500 15750 26800 29000 29000 16500 32900 16500 14000 14000 15000 11500 19800 29000 29000 24500 26000 24500 19800 24500 19500 29000 2 1 1 9 . . . . 55 53 28 25 83 66 31 27 13 11 12 10 63 50 41 34 18 16 34 29 16 14 33 29 22 18 15 13 . . 2 2 3 2 100 80 460 280 390 210 340 210 180 120 290 140 640 450 160 110 170 90 110 80 320 200 260 150 600 400 140 110 220 140 250 180 400 290 300 200 310 200 120 90 + − + − + − + − + − + − + − + − + − + − + − + − + − 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 8 1 . . 07 . 10 57 17 1150 23 780 14 221 20 630 . . . . . 039 184 034 689 070 161 036 1121 039 1170 098 1750 069 1660 068 1170 032 3850 027 805 078 308 056 1230 044 885 045 713 016 794 033 755 025 605 026 857 089 800 066 1640 092 1200 055 2910 085 375 024 1810 064 1340 029 2170 023 3160 053 1770 058 1720 062 1030 057 387 0 ...... 0 . . 0 . . . 0 . . 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 52 55 52 49 03 . . . . . 411 448 207 879 834 310 487 805 118 238 250 771 119 395 258 319 651 163 217 507 729 110 663 536 682 424 265 486 493 942 585 ...... 62:08:45 1 32:31:54 3 18:39:26 1 10:09:34 1 47:12:10 1 66:26:13 1 50:15:10 1 58:13:11 1 41:40:44 0 60:50:43 0 55:50:05 0 60:24:01 0 15:16:52 0 39:29:50 0 + + + + + + + + + + + + + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( NameHD 9672 01:34:37.9 RA -15:40:35 17 DEC Parallax Distance T HD 97048 11:08:03.2 -77:39:17 5 HD 98922 11:22:31.7 -53:22:11 1 HR 5999 16:08:34.3 -39:06:19 6 HT CMa 07:02:42.5 -11:26:12 0 HU CMa 07:04:06.7 -11:26:08 0 Hen 3-1121 15:58:09.6 -53:51:18 0 Hen 3-1121S 15:58:09.7 -53:51:35 0 Hen 3-1191 16:27:15.1 -48:39:27 0 Hen 3-823 12:48:42.4 -59:54:35 0 Hen 3-847 13:01:17.8 -48:53:19 1 Hen 3-938 13:52:42.8 -63:32:49 0 IL Cep 22:53:15.6 IP Per 03:40:47.0 KK Oph 17:10:08.1 -27:15:19 4 LKHa 260 18:19:09.4 -13:50:41 0 LKHa 338 06:10:47.1 -06:12:51 1 LkHa 208 06:07:49.5 LkHa 215 06:32:41.8 LkHa 257 21:54:18.8 LkHa 259 23:58:41.6 LkHa 324 21:03:54.2 LkHa 339 06:10:57.8 -06:14:40 1 MQ Cas 00:09:37.6 MWC 1021 20:29:26.9 MWC 1080 23:17:25.6 MWC 137 06:18:45.5 MWC 297 18:27:39.5 -03:49:52 2 MWC 342 20:23:03.6 MWC 593 17:49:10.2 -24:14:21 0 MWC 655 22:38:31.8 MWC 657 22:42:41.8 MWC 878 17:24:44.7 -38:43:51 0 MWC 953 18:43:28.4 -03:46:17 0 NSV 2968 06:26:53.9 -10:15:35 0 NV Ori 05:35:31.4 -05:33:09 2 continued. Table 2.

Article number, page 19 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 27 23 V Binary 9.80 Yes 10.81 - 10.6310.38 - - 10.23 - 14.03 - 13.11 - 12.44 - 11.78 - 12.15 - 13.40 - 13.13 - 12.20 - 12.84 - 13.31 - 13.79 - 13.26 - 12.06 - 13.00 - 10.02 - 12.15 - 14.55 - 12.02 - 14.45 - 12.34 - 14.04 - 13.15 - 12.85 Yes 13.54 - 14.18 - 10.88 - 13.42 - 12.80 - 14.38 - 15.17 - 12.40 - V 26 00 34 12 34 12 19 16 60 34 25 25 20 27 09 10 18 26 12 34 13 13 15 00 18 24 12 19 26 29 29 23 50 50 13 13 17 28 35 25 34 12 26 29 29 20 18 26 ...... 060 040 030 000 010 000 030 020 070 070 040 030 040 010 040 000 010 020 030 030 050 040 070 070 0 0 0 0 . . 0 0 0 0 0 0 . . 0 0 0 0 . . 0 0 . . 0 0 . . 0 0 . . . . 0 0 0 0 . . 0 0 0 0 0 0 . . 0 0 . . 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 00 91 20 44 64 90 67 43 29 99 03 00 31 21 57 96 20 81 84 48 73 28 07 12 ...... 0 0 0 1 2 0 0 1 5 2 5 0 3 1 3 2 1 5 4 1 1 4 4 7 170 230 070 220 510 110 600 270 520 860 760 600 ...... 1 1 2 1 3 1 2 6 0 0 1 1 ) (mag) (mag) 12 01 24 10 34 17 18 16 10 09 69 82 10 07 15 14 15 16 23 24 17 13 31 36 16 18 16 25 27 19 37 39 20 18 18 24 14 16 35 42 39 31 09 10 37 35 37 39 13 15 16 12 12 11 35 17 46 50 23 16 32 36 20 21

...... 064 057 088 082 074 016 075 054 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 . . 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 L 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 79 28 82 58 60 80 38 37 96 95 72 14 96 79 05 15 21 60 12 43 87 24 87 00 95 44 88 60 61 59 21 70 Log(L) A ...... 0 1 1 2 0 2 1 1 0 1 2 5 1 1 4 4 5 1 2 3 2 2 2 4 1 0 1 1 3 0 5 2 734 274 948 193 . . . . 1 1 0 1 ff e 250 100 900 300 900 300 180 160 250 250 160 220 100 150 250 250 500 500 800 900 500 500 250 100 250 250 250 250 500 500 380 180 250 250 500 500 900 300 250 100 1500 1500 9200 9200 2000 2000 3900 4500 2000 2000 1500 1500 3500 3500 3000 3000 1000 1500 4500 5000 5000 3000 1000 1000 3500 3500 4500 5000 3900 4500 2000 2000 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 6750 9800 9800 8150 7830 6550 6750 9750 8520 6250 9800 6750 10250 10500 14000 10700 10500 10750 15250 10500 15000 19800 14000 29000 17000 26000 17500 30000 13000 24500 19500 18500 17500 24500 29000 15000 8 3 8 3 2 7 ...... 14 12 53 44 19 18 64 51 51 42 60 46 62 53 92 76 10 9 65 54 42 37 36 31 43 37 38 32 56 37 33 28 7 7 6 7 4 4 140 90 140 110 370 260 750 540 600 400 460 310 210 150 140 110 330 240 140 110 170 140 430 300 360 230 160 120 190 120 640 410 240 160 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − − + − + + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 6 8 2 . . . 10 850 13 144 10 2470 19 380 . . . . 027 414 050 399 037 887 019 763 035 1490 048 853 040 847 069 684 020 1316 024 1475 033 2550 029 4630 035 398 031 1074 019 1130 054 2890 044 2550 036 345 034 1820 027 869 031 1590 023 2880 027 951 017 2130 016 2440 030 2950 056 1930 033 804 027 1810 083 1080 063 1410 044 643 . . . 0 ...... 0 . . 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 11 93 19 60 . . . . 0 411 505 119 308 649 160 173 447 754 668 328 043 514 922 881 118 275 891 515 146 610 307 046 457 395 266 416 238 531 864 639 551 ...... − 29:41:46 2 20:14:48 0 10:01:13 1 29:25:17 0 00:42:33 2 02:54:17 0 + + + + + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( NamePDS 002 01:17:43.5 RA -52:33:31 2 DEC Parallax Distance T PDS 004 03:39:00.6 PDS 021PDS 022 06:02:14.9 06:03:37.1 -10:00:59 1 -14:53:03 1 PDS 025 06:54:27.9 -25:02:16 1 PDS 123 05:50:54.8 PDS 124 06:06:58.5 -05:55:07 1 PDS 126 06:13:37.3 -06:25:02 1 PDS 129 06:31:03.6 PDS 130 06:49:58.6 -07:38:52 0 PDS 133 07:25:05.0 -25:45:50 0 PDS 134 07:32:26.6 -21:55:36 0 PDS 138 11:53:13.2 -62:05:21 0 PDS 174 05:06:55.5 -03:21:13 2 PDS 211 06:10:17.3 PDS 24 06:48:41.7 -16:48:06 0 PDS 241 07:08:38.8 -04:19:05 0 PDS 27 07:19:35.9 -17:39:18 0 PDS 277 08:23:11.8 -39:07:01 2 PDS 286 09:06:00.0 -47:18:58 0 PDS 290 09:26:11.1 -52:42:27 1 PDS 297 09:42:40.3 -56:15:34 0 PDS 324 10:57:24.2 -62:53:13 0 PDS 33 08:48:45.7 -40:48:21 1 PDS 34 08:49:58.5 -45:53:06 0 PDS 344 11:40:32.8 -64:32:06 0 PDS 361S 13:03:21.5 -62:13:26 0 PDS 37 10:10:00.3 -57:02:07 0 PDS 389 15:14:47.0 -62:17:00 1 PDS 415N 16:18:37.2 -24:05:18 6 PDS 431 16:54:59.2 -43:21:50 0 PDS 469 17:50:58.1 -14:16:12 0 PDS 477 18:00:30.3 -16:47:26 PDS 520 18:30:06.2 PDS 543 18:48:00.7 PDS 69 13:57:43.9 -39:58:47 1 continued. Table 2.

Article number, page 20 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 2 2 2 2 4 2 3 2 1 2 21 33 21 21 24 25 26 V Binary 9.53 - 8.87 - 9.26 Yes 7.73 - 8.90 - 9.10 Yes 11.51 - 10.92 Yes 11.36 - 10.97 Yes 10.43 Yes 10.34 Yes 12.04 - 13.15 Yes 15.34 - 10.63 Yes 12.93 - 13.84 - 11.07 - 10.90 Yes 11.42 - 10.21 - 10.82 Yes 10.03 - 12.73 Yes 10.21 Yes 10.53 Yes 14.69 - 11.75 Yes 10.61 - 12.99 - 13.76 - 12.41 Yes 11.41 Yes 11.82 - 11.39 Yes V 29 20 28 12 20 19 23 15 07 10 12 12 12 34 12 26 60 34 35 17 39 20 19 24 19 16 18 24 12 12 12 34 25 03 15 17 18 24 17 16 12 35 33 16 24 12 24 12 17 16 18 26 12 12 28 12 ...... 080 050 070 030 020 040 050 000 020 030 050 070 060 070 060 050 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 . . 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 21 55 96 01 89 96 98 40 91 75 75 85 08 72 26 50 03 90 07 00 00 78 90 97 05 88 91 81 ...... 1 1 0 2 1 0 1 0 8 5 3 1 4 3 2 2 0 1 3 3 4 3 1 5 3 4 2 0 500 480 320 000 690 210 650 740 ...... 1 0 1 0 0 2 4 0 ) (mag) (mag) 18 13 19 12 10 09 40 22 16 18 10 07 14 23 16 20 81 89 34 20 42 29 27 22 11 09 24 29 11 10 15 24 10 09 12 03 15 14 18 24 15 17 21 18 24 32 46 38 24 13 39 23 21 21 33 34 10 08 16 10

...... 080 080 064 044 048 050 045 022 057 057 046 039 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 . . 0 0 0 0 0 0 L 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 36 01 86 80 58 56 41 90 17 15 80 53 15 96 48 33 04 40 06 29 36 00 45 86 13 51 60 76 95 21 Log(L) A ...... 1 2 0 3 3 1 1 2 5 2 3 1 1 2 2 2 1 1 3 2 3 2 2 3 2 3 2 4 1 1 073 962 773 918 444 899 ...... 1 0 0 0 1 0 ff e 310 150 200 200 150 80 200 200 500 500 800 900 250 250 250 250 200 200 600 380 160 110 180 160 200 200 800 900 250 250 250 250 160 220 750 750 750 750 300 600 250 250 250 250 200 200 200 200 3000 800 1500 1500 9200 9200 5000 3000 1000 1500 1000 1500 800 1000 4500 4500 1500 800 1500 800 800 1000 3900 4500 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 6750 6250 9000 8500 8500 9500 8900 6010 8150 7750 9000 7830 9750 9500 8000 8000 10000 10250 10700 11000 10700 15750 13800 10000 16500 16000 19800 19500 13000 13000 15000 22500 11500 11500 15000 29000 0 3 0 8 7 2 3 4 3 1 0 2 ...... 50 39 49 41 47 40 13 11 13 12 46 39 14 13 82 59 39 34 36 31 10 9 35 27 35 31 33 29 40 35 76 49 16 14 12 11 14 13 18 16 9 8 4 3 4 4 9 8 8 9 3 3 250 160 100 70 320 190 210 130 280 160 540 340 250 160 280 170 120 90 240 160 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − − + − + + − + − + − + − + − + − + − + − + − + − 5 3 5 9 8 5 ...... 15 136 10 1300 12 1080 12 390 16 482 . . . . . 068 627 046 773 039 368 061 1480 031 932 020 344 044 408 051 324 051 386 067 870 088 1100 031 910 053 391 088 695 030 860 029 843 043 366 045 341 025 893 022 922 030 872 061 2470 068 1400 028 569 091 1280 042 410 066 980 081 152 054 1550 046 420 037 537 ...... 0 . . . 0 . . 0 . . . . . 0 . . . . 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 33 62 78 55 04 . . . . . 581 284 713 604 066 904 452 077 586 122 828 092 555 412 158 182 728 925 116 081 141 196 634 755 670 436 990 556 587 381 861 ...... 23:53:51 1 26:22:27 1 73:40:27 2 05:44:17 1 40:39:37 0 44:37:30 0 41:21:51 1 41:21:28 0 43:47:25 1 43:47:26 1 01:43:48 2 00:08:40 2 66:06:35 1 66:06:54 1 62:15:36 1 61:13:32 0 61:54:40 1 43:18:10 0 00:08:40 2 61:58:51 1 + + + + + + + + + + + + + + + + + + + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( NamePX Vul 19:26:40.3 RA DEC Parallax Distance T RR Tau 05:39:30.5 RY Ori 05:32:09.9 -02:49:47 2 SAO 185668 17:43:55.6 -22:05:45 0 SAO 220669 08:55:45.9 -44:25:14 1 SV Cep 22:21:33.2 T Ori 05:35:50.5 -05:28:35 2 TY CrA 19:01:40.8 -36:52:34 7 UX Ori 05:04:30.0 -03:47:14 3 V1012 Ori 05:11:36.5 -02:22:48 2 V1295 Aql 20:03:02.5 V1478 Cyg 20:32:45.5 V1493 Cyg 20:52:04.6 V1685 Cyg 20:20:28.2 V1686 Cyg 20:20:29.3 V1787 Ori 05:38:09.3 -06:49:17 2 V1818 Ori 05:53:42.6 -10:24:01 1 V1977 Cyg 20:47:37.5 V2019 Cyg 20:48:04.8 V346 Ori 05:24:42.8 V350 Ori 05:40:11.8 -09:42:11 2 V351 Ori 05:44:18.8 V361 Cep 21:42:50.2 V373 Cep 21:43:06.8 V374 Cep 23:05:07.5 V380 Ori 05:36:25.4 -06:42:58 2 V388 Vel 08:42:17.3 -40:44:10 0 V431 Sct 18:29:25.7 -06:04:37 0 V594 Cas 00:43:18.3 V594 Cyg 21:20:23.4 V599 Ori 05:38:58.6 -07:16:46 2 V669 Cep 22:26:38.7 V718 Sco 16:13:11.6 -22:29:07 6 V921 Sco 16:59:06.8 -42:42:08 0 VV Ser 18:28:47.9 VX Cas 00:31:30.7 continued. Table 2.

Article number, page 21 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 1 V Binary 9.67 Yes 12.86 - 10.64 - V 26 29 22 12 12 12 ...... 0 0 0 0 0 0 + − + − + − 56 95 47 . . . 3 0 1 ) (mag) (mag) 39 44 13 09 13 12

...... 0 0 0 0 0 0 L + − + − + − 68 42 00 Log(L) A . . . 3 1 2 ff e 200 200 200 200 4500 5000 + − + − + − 8970 9750 24500 16 15 35 27 300 200 + − + − + − , Mass and Age. The references for binarity are: (1) Baines et al. (2006); (2) Wheelwright et al. (2010); L 048 1860 037 504 089 459 . . . 0 0 0 ± ± ± 467 984 169 . . . 21:12:31 1 38:58:55 2 + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( and V taken from the following sources in order of choice: Fairlamb et al. 2015; Montesinos et al. 2009; Hernández et al. 2004; Mendigutía et al. et al. (2017); (19) Mayer et al. (2016); (20) Folsom et al. (2008); (21) Corporon & Lagrange (1999); (22) Doering & Meixner (2009); (23) Chelli et V ff ,A ff e NameWRAY 15-1435 16:13:06.7 -50:23:20 0 RA DEC Parallax Distance T WW Vul 19:25:58.8 XY Per A 03:49:36.3 Atmospheric parameters T continued. 2012; Carmona et al. 2010;If Chen not et al. available 2016; they Alecian were et derived al. as 2013; described Sartori et in al. Sect. 2010; 2.2. Manoj See et Sect. al. 2006; 3 Hernández for et derivation al. of 2005; Vieira et al. 2003; APASS Data Release 9 and the SIMBAD database. Table 2. References. (3) Leinert et al. (1997);(10) (4) Liu Maheswar et et al. al.(17) (2000); (2002); Morrell (11) (5) & Biller Wheelwright Levato et etal. (1991); al. al. (1995); (18) (2011); (2012); (24) Lazare (6) (12) Miroshnichenko Alecian et(2009); Schütz et (31) al. et al. Frasca (2002); al. (2013); et (25) (7) (2011); al. Friedemann Hamaguchi (13) (2016); et et (32) Boersma al. al. Marston et (1996); (2008); & al. (26) (8) McCollum (2009); Kraus (2008); Dunhill (14) (33) et et Zhang Malkov al. al. et et (2015); (2012); al. (9) (27) al. (2017). Coulson Torres (2006); et & (15) al. Walther Arellano (1995); (2000); Ferro (28) & Aspin Giridhar (1998); (29) (2003); Connelley (16) et Kubát al. et (2008); al. (30) (2010); Millour et al.

Article number, page 22 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 7 6 3 0 7 8 9 4 8 7 2 5 3 1 4 5 7 4 0 3 7 4 0 1 9 4 8 5 ...... 16 3 13 06 14 2 12 1 33 23 44 08 11 06 21 08 45 49 22 32 27 40 32 20 29 67 54 32 46 04 64 73 15 2 20 07 62 36 ...... 1 0 . . . . 0 1 4 1 . . . . 3 0 . . . . 5 0 3 3 0 2 . . . 2 2 . 0 1 2 3 . . . . 0 1 0 1 1 1 1 0 . . . . 053 068 077 026 + − + − + − + − 0 0 0 0 0 0 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − − + − + − + − + − + − + − + − + − + − + − + − + − + 0 0 0 0 5 4 6 8 + − + − + − + − + − + − + − + + − − + − + − + − + − + − + − 4 7 3 2 5 5 1 0 4 5 7 6 1 9 + − + − ...... 4 7 7 5 4 2 4 8 0 0 8 12 09 53 37 48 05 75 38 25 09 76 10 56 12 06 ...... 10 17 17 0 0 6 4 5 4 5 6 1 0 1 0 0 0 0 277 046 . . 0 0 ) (Myr) 9 0 6 8 1 1 6 8 8 4 2 7 1 4 2 2 2 4 8 8 ...... 37 12 29 14 66 56 16 08 14 08 10 12 11 11 59 34 25 09 49 19 11 11 36 21 37 19 15 10 20 24 19 11 60 65

3 1 4 2 ...... 4 3 . . 1 0 1 0 . . 2 1 5 2 ...... 3 2 ...... 4 3 2 1 064 071 094 094 063 063 092 092 087 087 070 070 090 090 061 061 . . 0 0 0 0 . . 0 0 0 0 0 0 . . . . 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . 0 . 0 . 0 0 0 0 0 0 0 0 . . + − + − + − + − + − + − + − + − + + − − 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 M + − + − + − + − + − + + − + − + − − + − + − + − + − + + − − + − + − 4 5 5 4 9 0 2 4 0 3 + − + − + − + − + − + − + − + − Mass Age ...... 05 10 09 80 22 28 37 06 60 47 30 04 04 51 53 71 15 ...... 2 1 1 251 849 736 807 215 275 880 401 ...... 1 1 i - 5 - 2 ∗ ∗ ∗ ∗ ∗ -* 2 V 39 33 06 25 09 . . . . . 0.50 1 0.35 6 2.93 1 0.01 5 0.98 9 0.24 3 9.76 4 0.21 2 -0.17 1 -0.04 2 -0.21 11 0 0 9 18.42 2 12.47 1 11 27 − - -0.19 2 - - 1 - 0.04 2 - -0.12 1 - 0.63 1 - 1.45 1 - 1.15 8 - -0.06 9 - 0.18 3 - 1.22 7 6 4 5 3 3 8 3 3 3 6 α 31 21 21 28 26 21 21 22 21 21 21 s s d d d d P P P s s d d d d d d P P P - - 0.43 13 - - 0.22 2 - s - - 49.58 1 1 5 5 5 1 1 1 5 1 1 5 5 1 6 6 1 5 1 3 3 3 3 1 27 11 15 15 27 23 22 3 3 6 7 1 7 2 5 6 4 5 1 ...... 0 5 72 65 19 22 14 63 52 01 38 . . EW H 1 1 1 1 2 . 1 . 1 . 1 . . . . 1 1 . 4 1 . 10 25 34 45 45 060 065 . . 4 1 . . 0 0 0 0 0 0 0 . 0 . 0 α 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± 0 0 ± ± ± ± ± ± ± ± ± ± ± 1 0 5 5 4 6 7 1 6 8 7 9 ± ± ± ± ...... ± ± 2 8 . . 5 1 2 5 24 03 03 89 20 50 80 10 55 200 . . 27 39 27 . 41 15 . . 55 . 14 . . . 61 42 21 09 78 − − − − . . . . 49 38 − 5 8 6 7 4 − − − − − − − − 451 773 5 − − 13 10 21 14 . . − − − − − 16 18 14 − 5 4 − − − − − − − − − 3 4 2 2 3 3 1 3 3 2 . . 08 10 09 13 14 14 05 11 25 12 − − − − − − − − ...... 0 1 045 037 032 032 015 019 033 019 013 012 077 042 012 015 043 027 049 041 044 063 092 056 038 043 071 047 016 035 019 026 027 036 016 024 m) (Å) line shape ( . . 0 0 . . 0 0 ...... 0 . 0 ...... 0 0 ...... 0 . . 0 . . − + 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ − + − + − + + − + − 10 10 10 10 10 10 10 10 7 − + − + − − + − + − + + − − + + − − + + − + − + − + − + − − + + + − . · · · · · · · · 0 22 37 27 29 27 15 ) ) ) ) ) ) ) ) . . . . . 0 0 0 0 0 − 157 197 162 120 309 184 151 078 168 196 139 170 072 049 173 173 162 38 62 69 30 30 50 04 12 67 52 24 61 41 72 73 82 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . − + − + − + − + − + − + − + − + 64 65 22 51 09 63 95 24 ...... (0 (3 (0 (0 (1 (6 (0 2 2 2 2 2 11 11 08 12 13 25 14 08 14 13 12 10 92 28 07 11 − − − − − ...... 055 013 043 021 033 079 053 056 028 016 068 036 055 071 031 042 011 015 009 012 022 032 009 010 040 062 035 053 011 012 054 031 013 015 m) (3 ...... 0 0 . . 0 0 . . . 0 0 . . . . 0 0 . 0 0 . . 0 0 0 0 . . . . 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − + − + − + − − + − + − + + − 10 10 10 10 10 + − + − + − + − + − − − + + − − + + − + − + + − + − + − + − + − + − 4 · · · · · . 17 20 31 19 56 28 20 32 ) ) ) ) ) 3 ...... 177 072 103 154 113 121 022 059 190 050 146 165 020 100 212 031 030 57 27 46 27 04 10 31 20 78 28 − ...... 0 0 0 0 0 0 0 0 0 0 + − + − − + + − + − 24 . 36 29 20 53 45 . . . . . Be star belonging to the high quality sample of 218 sources. (1 / 41 3731 - (0 30 549 0 + + GSC 3975-0579 0 GSC 6546-3156 0 GSC 8143-1225 0 GSC 8581-2002 (4 GSC 8645-1401 0 GSC 8994-3902 - - - - 0.04 8 HBC 217 0 HBC 222 0 HBC 705HD 100453 (0 0 HD 100546HD 101412HD 104237HD 114981 0 0 (0 0 HBC 334HBC 7 0 (0 HBC 442 0 HD 130437 (0 HD 132947 - - HD 135344 - - NameAB Aur Near IR excess Mid IR excess 0 H AK Sco 0 AS 310 - - BF Ori 0 BH Cep 0 AS 470BD BD 0 CQ Tau 0 BO CepCO OriGSC 1876-0892 0 0 0 AS 477 0 CPM 25DG Cir 0 0 Other parameters of each Herbig Ae Table 3.

Article number, page 23 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 8 5 5 4 6 4 8 3 3 7 4 6 1 2 2 7 8 3 5 ...... 16 1 54 11 4 26 20 20 15 20 37 15 2 54 24 12 3 43 52 67 53 28 37 37 49 54 31 01 14 11 13 2 45 91 11 1 25 55 64 54 44 27 40 41 45 20 29 57 77 47 − ...... 0 1 9 3 ...... 1 1 1 0 . . . . 0 1 ...... 0 1 0 1 . . 1 3 0 . . 1 057 027 + − + − + − + − + − + − 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . 0 0 0 0 0 . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + + − − + − + − + − + − + − 0 0 4 9 5 8 7 9 + − + − + − + + − + − + − − + − + − + − + − + − + + − − + − + − + − + − + − 10 7 8 9 2 4 5 6 7 6 + − ...... · 4 3 6 7 7 3 6 59 04 17 66 74 48 89 98 44 22 05 41 33 43 89 56 42 43 93 ) ...... 10 14 4 0 0 1 1 5 0 4 0 1 4 0 9 7 4 2 0 7 8 041 76 53 . . . 0 0 0 + − 37 . (2 ) (Myr) 8 3 1 0 9 2 4 2 4 9 0 2 3 5 ...... 30 10 30 10 49 14 10 18 16 10 11 11 39 13 16 13 08 08 18 09 18 30 39 29 81 47 14 30 79 22 13 13 10 10 12 10 14 08

2 2 2 1 . . . . 6 3 . . . . 0 1 ...... 1 0 ...... 1 1 ...... 2 2 . . . . 093 085 085 085 082 069 072 072 074 074 093 095 095 093 069 069 092 092 075 075 0 0 0 0 0 0 0 0 . . 0 0 0 0 . . 0 0 . . 0 0 0 0 0 0 0 0 . 0 . 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 ...... 0 0 . . + − − + − + − + − + + − + − 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 M + − + − + − + − + − + − + − + + + − − + − − + + − − + − + + − − + − + − + − 3 7 5 9 8 3 2 + − + − + − + − + − + − + − + − + − + − Mass Age ...... 22 69 98 94 30 97 35 00 03 56 05 60 98 25 92 22 53 90 61 ...... 1 891 833 691 704 383 432 481 860 386 493 ...... 1 i - 1 ∗ ∗ V 30 01 . . 0.23 7 0.11 1 0.32 1 0.20 1 0.00 2 0.52 12 0.08 1 0.07 5 0.16 1 0.63 2 0.07 1 0.00 2 0.38 1 0.35 5 0.20 1 0.24 1 2.05 1 0.65 19 0.38 1 -0.17 1 -0.09 17 -0.04 2 -0.19 1 -0.18 1 -0.13 2 -0.16 2 -0.13 3 0 3 2 3 3 3 3 3 6 7 7 6 3 8 α 31 31 21 21 21 12 21 13 12 31 21 22 21 26 31 14 21 21 s s s s d d d d d d d P s s s s s s s s s s s d d d d P P P - - 0.27 3 - - -0.21 2 - - 0.53 2 - - 0.31 2 1 1 1 1 7 1 1 2 1 1 1 1 1 1 3 1 1 1 3 3 1 1 6 1 15 14 12 32 13 12 7 4 4 3 7 9 7 1 6 5 1 6 6 5 7 2 4 ...... 0 52 41 10 68 94 55 35 . EW H . 1 1 1 1 . 1 . . . 5 1 1 1 1 2 1 1 . 1 1 . 1 1 15 20 70 15 . . . 2 . 0 0 0 0 0 0 0 050 α 0 0 . 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 0 ± ± ± ± ± ± ± 8 1 7 8 6 9 5 0 5 4 7 6 6 8 4 1 6 ± ± ± ± ...... 7 . ± 5 6 24 26 60 30 80 80 88 . 63 29 58 12 28 76 . 38 37 15 14 14 . 12 10 21 . . . . 13 59 04 79 21 − − . 47 . . . 3 3 − − − − − − − − − − − − − − − − 13 22 10 46 11 235 − − 22 . 11 10 − − − − − 8 − − − − - 1 - 3 2 3 3 2 3 3 06 13 19 06 10 10 − − − − − − − ...... 016 011 033 040 067 032 030 040 012 047 020 007 007 051 027 015 028 022 050 089 016 011 055 061 016 010 060 077 039 049 007 013 020 054 012 016 013 017 030 024 005 015 046 038 m) (Å) line shape ( . . . . 0 0 ...... 0 0 ...... 0 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − + − + − 10 10 10 10 10 10 10 − + − + − + − + − + + − − + − − + − + − − + + − + − + − + + − + − + − + − + − + − + − + · · · · · · · 22 32 11 53 ) ) ) ) ) ) ) . . . 0 0 0 − 208 136 114 121 221 171 211 161 029 186 023 199 039 221 260 256 175 357 145 193 333 158 35 15 47 12 21 22 80 80 13 04 67 74 10 07 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . − − + + − + − + + − + − + − 04 62 31 75 21 41 26 ...... (4 (0 (0 (1 (8 (0 3 2 2 3 2 2 14 06 03 11 10 04 10 07 15 10 − − − − − − ...... 014 030 009 010 032 034 022 023 012 019 048 042 011 029 020 023 024 018 017 008 013 020 011 010 018 065 058 006 044 050 024 015 024 018 009 010 010 016 022 029 026 044 049 043 m) (3 ...... 0 0 0 0 ...... 0 . 0 . . 0 0 ...... 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ − + − + − + − + + − 10 10 10 10 10 10 − − + + − − − + − + − + + − + − + − + − − + + − + − + + + − + − − + − + − + + − + − + − 4 · · · · · · . 24 33 07 24 59 ) ) ) ) ) ) 3 . . . . . 152 038 013 162 016 175 028 032 041 173 286 169 065 133 128 115 245 072 248 106 225 196 15 10 61 66 34 23 35 49 43 20 16 00 − ...... 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − − + − + 24 . 87 24 20 00 92 97 ...... (1 HD 200775 0 HD 235495 0 HD 244314 0 HD 250550 0 HD 244604HD 249879 0 0 HD 245185HD 259431 0 0 HD 287823 0 HD 288012 - - - - -0.15 2 HD 305298HD 313571 (0 (1 HD 290770 0 HD 290380HD 290500HD 290764 0 0 0 HD 290409 0 NameHD 135344B Near IR excess 0 Mid IR excess H HD 139614 0 HD 141569 (1 HD 179218 0 HD 141926 (1 HD 149914HD 150193HD 155448 0 0 HD 17081HD 174571 (0 - (0 HD 199603 - - (0 - - -0.09 4 HD 158643HD 163296HD 169142 0 HD 176386 0 0 0 HD 142527HD 144432 0 0 HD 143006 0 HD 142666 0 continued. Table 3.

Article number, page 24 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 8 3 4 4 1 2 2 3 4 6 8 0 1 4 4 2 6 2 2 ...... 14 1 14 10 12 05 19 18 12 18 24 20 11 08 17 12 23 44 15 2 07 17 90 50 55 26 17 09 18 13 47 17 12 2 87 93 65 34 11 07 64 60 − 2 0 3 2 . . . . 0 1 ...... 9 0 0 1 . . 3 0 . . . . 0 1 . . . . 1 1 . . . . 0 1 . . 037 022 023 028 045 051 052 010 050 071 + − + − + − 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 . . . . 0 0 . . 0 0 . . + − + − + − + − + − + − + − + − + − 0 0 0 0 0 0 0 0 0 0 5 8 3 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 10 3 6 2 5 1 3 8 4 6 + − + − + − + − + − ...... · 8 7 6 6 8 1 4 0 40 46 99 22 08 73 69 30 86 10 56 22 09 63 84 96 28 28 ) ...... 12 1 1 0 0 1 0 1 0 0 0 1 0 0 0 0 0 0 0 047 171 217 020 280 79 54 ...... 0 0 0 0 0 0 0 + − 90 . (1 ) (Myr) 1 8 9 9 2 0 2 0 4 8 9 6 3 8 2 4 4 4 7 8 ...... 46 32 11 08 33 19 31 26 51 37 53 39 65 53 11 7 50 27 13 09 32 13 25 25 15 15 11 08 94 75 36 38 42 29 14 11 63 67 55 33 79 38 43 24 23 10

...... 1 1 . . 1 0 . . . . 3 1 ...... 1 0 1 0 3 2 3 1 ...... 1 0 . . 3 2 . . 1 0 . . 092 092 072 072 067 067 + − 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + + − − + − 0 0 0 0 0 0 M + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 7 7 9 2 0 6 8 7 6 4 + − + − + − Mass Age ...... 10 66 92 56 71 20 88 84 85 11 21 37 00 17 53 87 81 46 32 51 50 78 ...... 1 343 833 450 . . . i - 2 ∗ V 22 . 0.71 1 0.30 3 0.56 3 1.55 8 0.49 3 0.02 2 0.33 5 0.46 4 0.05 11 0.62 4 0.24 1 0.65 7 0.17 6 5.31 18 0.10 5 0.05 20 0.01 1 0.00 5 0.91 3 -0.17 3 -0.19 5 -0.10 5 -0.09 2 -0.26-0.04 2 1 -0.23 6 2 - -0.19 3 - -0.12 3 3 7 7 3 8 8 8 6 7 6 3 8 6 3 α 22 22 21 22 21 21 22 21 21 31 31 34 31 21 s s s s s d d d d d d d d s s s s s s d d d d d d P P - - - 3 - - 0.78 2 -P - - -0.21 9 7 6 3 1 1 1 1 1 6 1 1 1 3 3 1 7 1 1 1 9 9 6 22 15 22 15 15 11 15 1 7 6 5 7 3 1 2 7 1 2 0 ...... 4 7 8 63 18 85 30 67 14 82 89 58 . . . EW H 1 . 1 1 1 . . 5 1 7 . 6 1 . 2 1 . 1 . . 20 20 10 22 072 . . . 1 2 . 1 . 0 0 0 0 0 0 0 0 ± α 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± 0 ± ± ± ± ± ± ± ± ± ± ± 2 3 9 8 8 7 3 0 8 0 9 1 ± ± ± ± ...... ± 1 6 4 . . . 7 9 1 3 145 40 51 33 07 40 84 17 75 61 47 81 . 16 . 42 . 21 . . . 12 15 . . 28 67 17 37 − − − − − 45 . . . 31 59 . 8 0 − − − − − − − − 400 2 9 − − − . 16 25 11 16 23 19 − − 13 12 − − 2 − − − − − − − − − - - - -0.24 1 - 2 2 2 2 3 3 3 4 3 3 4 2 10 10 12 09 28 27 19 14 38 13 − − − − − − − − − − − − ...... 031 048 076 077 020 012 046 009 023 010 025 006 087 076 025 019 009 015 025 020 044 048 032 053 036 022 m) (Å) line shape ( 0 0 0 0 0 0 0 0 ...... 0 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − + − + − + − + − 10 10 10 10 10 10 10 10 10 10 10 10 − + − + − − + − + − − + + − − + + + − + − + − + + − · · · · · · · · · · · · 22 22 16 51 34 19 ) ) ) ) ) ) ) ) ) ) ) ) . . . . . 0 0 0 0 0 − 054 138 089 284 167 060 258 050 136 147 085 085 200 58 35 09 13 12 16 28 29 26 14 70 78 08 16 73 79 16 36 56 69 21 21 26 21 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 . − + − − + + + − + − + − + − + − + − + − + − + − 29 22 53 83 41 12 26 23 27 80 16 27 ...... (4 (0 (0 (0 (1 (0 (7 (2 (3 (1 (0 2 2 3 2 2 3 2 2 2 08 22 10 06 28 24 − − − − − − − − − ...... 067 011 052 070 015 052 009 032 070 056 067 062 018 019 017 009 015 016 036 028 066 027 056 024 019 028 079 090 014 017 051 022 039 024 027 025 007 012 m) (3 . . 0 0 ...... 0 0 ...... 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − + − − + 10 10 10 10 10 10 10 10 10 + − + − + − + − + − + − − − − + + + − + − + − + − + − + − + − + − + + − + − 4 · · · · · · · · · . 10 45 11 ) ) ) ) ) ) ) ) ) 3 . . . 024 074 239 047 142 086 187 022 089 009 146 238 046 129 079 228 166 183 087 31 32 23 09 43 39 43 14 31 11 31 32 37 55 73 05 34 09 − ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − − + − + − + − + − + − + − + 24 . 05 09 42 33 46 70 23 44 96 ...... (1 HD 46060 - (0 HD 50083 (0 HD 50138 0 HD 56895B - - - - -0.16 1 HD 58647 0 HD 59319 (0 HD 68695 0 HD 76534 (2 HD 85567 0 HD 87403 (1 HD 87643 0 HD 94509 (1 HD 95881 0 HD 96042 (1 NameHD 31648 Near IR excess 0 Mid IR excess H HD 344261HD 34700HD 35187HD 35929 0 0 0 - - HD 319896 (0 HD 36112HD 36408HD 36982 (0 0 - - - - -0.13 5 HD 34282 0 HD 36917 0 HD 323771 0 HD 37258HD 37357HD 37371 0 0 - - HD 37490HD 37806HD 38087 (0 HD 38120 0 0 0 HD 39014 0 HD 41511 0 HD 45677 0 continued. Table 3.

Article number, page 25 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 0 5 5 1 5 3 0 3 5 4 7 8 1 1 6 9 1 3 3 1 7 2 2 ...... 45 02 54 04 50 01 47 03 40 10 34 11 5 26 38 51 37 11 11 2 44 92 53 04 23 16 34 15 26 35 15 1 − 8 1 1 0 8 3 1 1 8 0 . . . . 1 1 1 0 ...... 1 0 8 6 ...... 1 . 1 . . 1 0 . . 019 000 019 008 020 006 086 056 063 028 012 006 044 033 010 038 + − + − + − . . . . 0 0 0 0 0 0 0 0 0 0 0 0 . 0 0 . 0 0 ...... 0 0 0 0 0 0 . . 0 0 0 0 . . + − + − + − + − + − + − + − + − + − + − + − 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9 9 5 + − + − + − + − + − + − + − + − + − + − + − + − + − 10 9 4 6 0 5 2 6 4 4 3 1 + − + − + − + − + − + − + − + − ...... · 4 2 0 2 3 6 4 0 03 23 12 54 08 04 04 73 07 02 05 20 89 ) ...... 12 18 11 1 0 2 2 0 2 0 2 0 0 0 0 6 017 070 010 018 016 123 042 204 62 61 ...... 0 0 0 0 0 0 0 0 0 0 + − 75 . (2 ) (Myr) 6 3 0 0 6 0 7 3 6 5 1 4 2 2 3 2 9 2 0 0 8 3 2 0 7 9 ...... 15 9 25 11 60 15 11 12 53 31 47 14 59 37 11 15 15 10 13 61 20 13 13 13 37 31 30 09 2 14 19 4 12 15 15 11 7 12 19 12

6 6 . . 2 1 1 0 2 0 ...... 5 3 . . . . 6 4 1 1 9 6 1 1 4 3 5 3 6 3 2 1 ...... 090 090 076 076 094 094 + − + − + − + − 0 0 0 0 . . 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − 0 0 0 0 0 0 M + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 8 8 8 1 1 1 9 4 0 5 5 7 0 + − + − + − Mass Age ...... 9 00 56 03 56 82 08 70 82 59 80 02 12 84 25 17 43 ...... 6 2 513 885 810 . . . i - 1 - 12 - 18 ∗ ∗ ∗ -* 1 V 15 02 51 . . . 1.242.03 3 25 0.46 4 0.54 1 1.87 3 0.33 2 1.83 10 4.02 16 0.96 23 0.58 2 1.44 16 3.84 22 0.62 8 0.09 13 3.85 12 1.28 1 0.96 3 -0.18 1 -0.10 2 -0.13 2 0 0 59.72 9 23 − − - 3.19 1 - 2.00 11 7 3 2 7 2 8 2 7 6 7 2 2 2 3 4 α 21 21 26 25 26 26 31 22 22 10 37 21 31 s s s s s s s s d d d d d P P s s s d d d d d d d d d P - - 50.67 3 - - 0.37 3 - - - 8 - - - 32 - - -0.14 1 2 2 3 1 7 5 1 5 5 5 5 1 2 5 1 1 2 2 7 1 6 1 1 2 32 37 15 27 17 10 1 7 6 3 2 2 9 4 2 0 8 7 6 4 1 ...... 5 0 87 54 24 76 37 12 00 99 . . EW H . . 1 1 2 . 1 1 1 . 1 . 3 1 . 1 1 2 1 . 1 1 80 72 13 28 050 4 . . 2 0 0 0 0 0 0 . 0 ± α 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 0 ± ± ± ± ± ± ± ± ± 6 5 1 2 3 8 4 2 5 6 2 2 0 7 8 ± ± ...... ± 0 2 . . 371 62 83 16 97 29 09 00 291 186 55 32 . . . 30 41 59 . 31 31 29 . 38 18 24 . 32 . 39 64 49 − . 41 93 . − − 1 113 594 4 − − − − − − − − − − − − − 590 − − 22 14 21 . 61 − − − − 15 15 108 1 − − − − − − − − 7 9 2 3 2 2 2 3 2 3 3 2 . . 47 67 37 52 30 12 43 23 15 42 02 25 14 02 − − − − − − − − − − ...... 0 1 ...... 018 008 073 035 037 024 032 021 021 003 016 007 013 006 010 006 035 015 064 021 030 049 022 022 019 019 017 018 085 088 m) (Å) line shape ( . . . . 0 0 . . . . 0 . . 0 . . 0 . . 0 0 . . 0 . . 0 0 0 0 . . . . 0 0 ...... − + 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ − + − + − − + + − + − + + − 10 10 10 10 10 10 10 10 10 10 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − . · · · · · · · · · · 3 22 03 72 58 02 84 30 02 ) ) ) ) ) ) ) ) ) ) ...... 0 1 1 1 0 0 0 − 021 155 080 071 139 219 004 163 019 596 013 015 027 041 173 44 53 71 81 45 64 71 77 35 82 31 12 61 37 67 28 53 73 90 92 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . − − + − + + − + − + + − + − + − + − + − 26 63 79 73 79 07 99 28 96 87 ...... (1 (1 (1 (0 (0 (0 (0 (10 (12 2 2 2 2 2 2 2 2 2 15 08 34 25 50 18 11 09 13 06 10 04 16 01 − − − − − − − − − ...... 016 025 006 029 008 018 005 014 014 085 006 018 082 056 012 005 029 022 018 034 021 012 008 017 016 015 017 019 061 064 023 011 m) (3 . . . . 0 0 . 0 0 . 0 0 . . . . 0 0 0 0 0 0 . . 0 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − + − + − + − + − + − + − 10 10 10 10 10 10 10 10 10 + − + − − + + − + − − + − + + − − + − + − + − + − + − + + − + − 4 · · · · · · · · · . 25 24 72 17 32 10 01 ) ) ) ) ) ) ) ) ) 3 ...... 018 126 050 067 026 006 021 021 011 017 012 126 109 189 180 419 27 38 32 78 27 42 32 31 76 26 15 16 80 36 52 54 60 40 − ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − − + − + − + + − 24 . 50 26 63 23 84 17 27 83 27 ...... (1 MWC 657 0 MWC 878 0 MWC 953 (0 NSV 2968 (0 NV Ori 0 NameHD 9672 Near IR excess (0 Mid IR excess H Hen 3-938IL Cep 0 - - Hen 3-847IP Per 0 0 KK OphLKHa 260LKHa 338LkHa 215LkHa 257 0 1 0 0 0 HD 97048 0 Hen 3-823 (0 LkHa 208LkHa 259LkHa 324LkHa 339 0 (0 (1 0 Hen 3-1191 0 MQ Cas - - - - -* 1 HD 98922 0 HR 5999Hen 3-1121S 0 - (0 MWC 1021 0 Hen 3-1121 0 HT CMa 0 MWC 1080 0 HU CMa (6 MWC 137 (0 MWC 297 (1 MWC 342 0 MWC 593 - - MWC 655 0 continued. Table 3.

Article number, page 26 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 2 5 1 4 5 6 1 9 0 3 8 3 4 8 9 4 2 6 3 8 4 5 2 1 3 1 2 2 2 9 4 4 5 ...... 10 4 14 1 18 1 39 87 15 1 63 35 70 48 10 03 45 73 13 3 27 26 14 1 14 1 22 21 10 1 − . . . . 5 4 . . 1 0 5 4 3 0 1 1 8 0 5 0 2 2 . . 1 1 8 2 ...... 1 0 9 3 4 0 0 2 3 0 0 1 012 001 036 028 012 072 003 027 + − + − + − + − + − + − + − + − 0 0 0 0 0 0 . . 0 0 0 0 ...... 0 0 0 0 + − + − + − + − + − + − − + − + − + − + + − − + − + − + + − + − 0 0 0 0 0 0 0 0 3 6 3 2 5 7 1 + − − + − + − + − + − + − + 10 3 1 2 1 7 0 6 8 7 8 9 0 2 9 7 3 + − + − + + − − ...... 10 · 0 0 4 0 0 0 1 1 7 2 4 68 06 77 48 73 76 75 ) ...... 13 14 11 10 10 1 0 2 3 0 7 1 011 078 013 042 58 12 ...... 0 0 0 0 0 0 + − 12 . (1 ) (Myr) 5 4 5 5 3 2 5 0 3 4 0 1 0 7 6 2 3 3 7 3 4 2 ...... 30 27 12 31 36 33 16 09 27 15 15 06 55 15 73 51 17 23 32 09 51 45 21 16 12 12 73 45 44 24 10 12 23 17 52 38 10 10 29 28 36 23

. . 4 3 . . . . 3 2 . . 5 9 . . 1 0 ...... 2 1 ...... 5 7 8 . . 3 . . . . 4 5 . . . . 2 2 1 1 ...... 3 2 077 077 093 093 092 069 061 061 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . 0 0 . 0 0 0 0 . . . . 0 0 0 0 0 0 . . + − + − + − + + − − + − + − + − + + − − + − 0 0 0 0 0 0 0 0 M + − + − + − + − + − − + − + − + − + − + − + + − − + − + − + − + − + + − + − + − + − 7 9 3 7 3 0 1 0 2 3 2 + − + − + − + − Mass Age ...... 91 21 52 26 15 18 48 57 33 93 28 07 94 95 41 71 35 99 90 94 46 ...... 2 543 850 384 228 . . . . i - 1 - 2 - 2 ∗ V 08 . 5.79 1 4.84 10 0.711.31 2 8 0.78 5 0.18 2 0.12 11 0.51 2 0.05 1 0.41 28 2.35 1 1.04 2 1.15 31 3.28 5 0.18 1 1.01 3 -0.14 1 -0.09 2 -0.11 1 -0.00 2 -0.20 2 -0.13 30 -0.11 4 -0.07 3 -0.12 1 -0.12 4 15.40 12 10.10 5 17.20 1 13 - - 2 - 0.49 7 4 4 4 α 31 31 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 21 25 31 21 21 25 21 21 21 21 21 21 31 P P P s s s s s s s s s s d d d d d d d d d d d d d d P P P P - - 0.25 1 -P 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 15 15 15 15 15 15 15 15 15 15 13 15 15 15 15 4 2 5 0 3 0 2 5 2 6 ...... 6 0 5 2 8 73 51 85 57 61 70 90 83 . . . . . EW H 1 1 . 1 . . 4 1 . . 2 . 2 1 . 1 1 . 45 10 25 15 35 20 20 45 65 . 1 . 6 2 . . 1 . . . 1 . . 0 0 0 0 0 0 0 0 050 α 0 0 . 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 0 ± ± ± ± ± ± ± ± 9 2 4 9 6 5 7 9 9 0 ± ± ± ± ± ± ± ± ± ...... 5 2 1 2 8 . . . . . ± 1 9 3 76 32 03 02 20 57 60 77 . . 38 60 . 18 . 76 . . . 42 . 28 66 19 94 78 45 78 26 16 53 − − − . 43 . . . 36 . 50 . . 34 . . 9 110 − − − − − − 5 7 121 8 6 − − − − 18 12 77 30 30 96 190 − − 18 . 14 12 18 22 − − − − − 123 − − − − − − 2 − − − − − − − - - 3 4 3 3 2 35 33 14 06 10 04 16 03 11 11 − − − − − ...... 042 045 012 020 017 034 060 022 025 059 051 018 086 017 034 054 029 017 024 037 019 014 055 050 041 026 005 030 015 023 044 028 004 029 020 040 031 024 034 071 m) (Å) line shape ( . . . . 0 0 . . . . 0 0 ...... 0 0 ...... 0 0 . . 0 0 . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ − + + − + − + − + − 10 10 10 10 10 − + − + − + + − + − − + + − + − + − + − + − − − + + − + + − − + − + − + + − + − · · · · · 22 12 15 31 03 59 ) ) ) ) ) . . . . . 0 0 1 0 0 − 116 242 057 087 048 058 052 112 102 087 182 151 254 164 107 127 371 214 090 199 19 45 15 14 79 27 66 44 57 35 ...... 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . − + + − + − + − + − 48 27 67 50 14 . . . . . (0 (0 (4 3 3 2 2 2 2 39 40 19 03 12 12 − − − − − − ...... 017 047 020 046 021 022 008 015 046 012 047 021 091 056 007 012 049 027 014 026 020 009 035 022 019 010 049 017 019 035 019 014 013 024 011 012 060 036 086 056 065 056 m) (3 ...... 0 0 ...... 0 0 . . 0 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − + − + − 10 10 10 10 10 10 + − + − + − − + + − + − + − + − − + − + − − − + + + + − + − + − + − − + + − + − + − 4 · · · · · · . 35 04 60 ) ) ) ) ) ) 3 . . . 056 086 045 019 036 157 104 111 141 044 143 025 054 128 107 039 017 055 070 051 038 53 42 42 50 46 80 29 32 26 17 63 92 − ...... 0 0 0 0 0 0 0 0 0 0 0 0 + − − + − + − + + − + − 24 . 62 56 04 20 57 46 ...... (1 PDS 241 (2 PDS 27 0 PDS 24 0 PDS 277 0 PDS 286 (1 PDS 290 - (0 PDS 297 (0 PDS 324 - - PDS 415N 0 PDS 33PDS 389 0 PDS 431 0 (1 PDS 469PDS 520PDS 543 0 1 - (0 PDS 37PDS 477 0 0 PDS 34 0 PDS 69 - - PDS 361S 0 PDS 344 (1 NamePDS 002 Near IR excess Mid IR 0 excess H PDS 004PDS 021 0 0 PDS 022 0 PDS 130PDS 133 0 0 PDS 129PDS 134 0 - - PDS 123PDS 126 0 PDS 0 138 - - - d PDS 025PDS 124 (9 0 PDS 174PDS 211 0 0 continued. Table 3.

Article number, page 27 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 8 7 2 0 0 8 1 2 8 5 2 1 9 7 6 7 1 9 1 6 6 1 ...... 11 4 39 19 19 19 91 37 47 47 15 13 75 60 72 35 50 44 30 07 49 78 42 54 40 69 14 2 56 67 11 07 13 1 73 67 11 07 . . 8 0 . . 2 0 . . 8 2 1 0 ...... 1 0 0 1 ...... 7 4 0 1 . . . . 1 0 ...... 1 1 . . 0 1 095 065 026 012 091 046 073 047 082 049 043 000 + − + − + − 0 0 . . 0 0 0 0 . . . . 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − 0 0 0 0 0 0 0 0 0 0 0 0 9 6 6 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 2 5 0 8 8 5 2 4 5 2 4 + − + − + − + − + − + − ...... 4 2 2 9 3 7 8 1 7 96 89 29 15 22 16 10 37 84 13 33 41 63 68 10 98 ...... 11 12 0 1 4 4 0 1 0 0 0 1 9 0 1 0 0 1 023 010 188 075 081 127 ...... 0 0 0 0 0 0 ) (Myr) 0 9 0 9 5 3 0 0 1 8 6 6 2 1 5 9 8 1 3 1 ...... 35 26 48 64 19 10 69 48 51 39 59 38 49 48 14 14 59 18 09 23 10 10 32 22 46 19 15 08 22 19 14 11 80 76 4 15 72 55

7 5 ...... 1 0 3 1 . . 1 1 ...... 1 0 . . 1 0 4 2 . . 3 1 ...... 3 1 1 1 094 083 079 079 094 085 080 080 091 081 065 065 077 077 + − 0 0 0 0 . . 0 0 0 0 0 0 . . . . 0 0 0 0 . . 0 0 . . 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 − + − + − + − + − + − + − + − + − + − + 0 0 0 0 0 0 0 0 0 0 0 0 0 0 M + − + − + − + − + − + − + − + − + − + − + + − − + − + − + − + − + − + − 8 1 4 0 7 4 9 3 1 3 + − + − + − + − + − + − + − Mass Age ...... 89 00 88 03 06 89 11 85 89 50 98 31 18 82 94 10 82 54 ...... 1 3 3 2 605 612 300 659 572 706 550 ...... 1 1 1 1 1 i ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ -* 2 -* 2 -* 2 -* 1 -* 2 V 70 91 78 23 91 67 25 12 49 ...... 0.67 4 3.49 20 3.69 2 1.020.47 9 1 1.37 3 0.07 1 2.29 1 2.88 6 1.28 2 2.10 9 2.06 2 -0.21 5 -0.09 5 3 7 2 3 4 11.23 5 11 16 10 22 - -0.07 7 - 7.82 3 3 3 3 3 3 3 3 7 3 3 3 8 8 3 7 4 3 α 18 19 21 20 22 21 31 31 21 21 26 31 s d d d d d d d d d d d d d d d P s d d d d d d d d P P P - - 1.39 2 - - - 28 - - - 4 - - 11.43 8 - - -0.24 9 3 1 1 5 3 7 1 3 3 1 1 3 3 3 1 1 3 1 3 1 5 7 3 27 27 17 11 18 20 19 22 5 5 5 5 7 4 4 4 6 4 1 7 4 ...... 0 0 6 66 77 48 68 98 43 77 47 73 35 . . . EW H 1 . 1 1 1 1 . . 1 . . 4 . . 1 1 1 2 1 1 . . . 19 60 60 31 4 . 2 1 0 0 0 0 0 0 0 0 0 0 050 α 0 . ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 0 ± ± ± ± ± ± ± ± ± ± 7 2 3 6 2 9 6 6 9 4 5 7 2 ± ...... 7 6 7 . . . ± 6 4 10 40 60 80 60 12 70 69 70 70 132 195 126 . 12 32 69 63 95 49 . 11 . . . 12 11 . . . 17 . . 16 − − . 45 39 80 − − − 1 147 − − − − − − − − − − − − − 22 25 18 22 32 14 25 15 11 165 − − 21 . − − − − − − − − − − 10 − - 2 2 3 2 3 08 05 13 12 20 53 02 16 15 06 78 35 12 09 61 18 34 14 13 13 − − − − − ...... 037 045 063 019 040 019 010 011 048 035 009 062 030 011 080 059 066 058 072 069 063 059 040 031 012 010 013 038 039 033 049 054 040 051 037 033 m) (Å) line shape ( . . 0 0 ...... 0 0 0 0 0 0 0 0 . . . . . 0 . 0 0 0 0 0 ...... 0 . . . . 0 . 0 0 ...... 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − − − − + + − + − + + − + − + − + − + 10 10 10 10 10 + − − − + + + − − + + − + − + − − + − + − + − − − + + − + + + − + − + − · · · · · 22 16 12 32 09 22 29 30 09 54 02 ) ) ) ) ) ...... 0 1 0 1 0 0 0 0 0 0 − 227 158 180 065 102 197 169 157 119 399 159 216 309 348 287 145 142 147 29 94 15 12 80 52 12 15 51 32 ...... 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . + − + − + − + − + − 82 45 69 32 19 . . . . . (0 (1 (1 (12 2 2 2 09 20 05 15 13 11 36 37 11 11 21 03 − − − ...... 069 049 022 050 059 061 022 050 030 037 012 038 021 011 010 012 070 023 058 078 048 022 024 026 041 039 032 020 024 011 028 026 040 054 080 071 031 046 040 012 014 012 052 042 075 080 065 091 m) (3 ...... 0 0 ...... 0 0 ...... 0 0 . . 0 0 . . 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − − + − + − + + − − + 10 10 10 − + − + − + + − + − + − − + − + − + + − − + − + − + − + − + − + − + − + + − + − + − + − + − + − 4 · · · . 82 03 27 23 29 29 ) ) ) 3 ...... 136 177 216 066 027 367 042 160 086 219 239 054 151 040 108 124 181 266 036 167 155 151 207 397 23 31 44 67 28 15 − ...... 0 0 0 0 0 0 − + − + + − 24 . 49 29 12 . . . (1 V1977 Cyg 0 V2019 Cyg 0 V346 Ori 0 V350 Ori 0 V351 Ori 0 V361 Cep (0 V373 Cep 0 V374 Cep (1 V380 Ori 0 V388 Vel 0 V431 Sct 0 V718 ScoVV SerVX Cas 0 0 0 V594 CasV594 CygV669 CepV921 Sco 0 0 0 0 V599 Ori 0 NamePX Vul Near IR excess Mid IR excess 0 H RR Tau 0 RY Ori 0 SAO 185668 - (0 UX Ori 0 SAO 220669 (0 TY CrAV1012 OriV1295 Aql 0 0 0 V1478 Cyg 0 SV CepT Ori 0 V1493 0 CygV1685 CygV1686 Cyg 0 - - - - V1787 Ori 0 V1818 Ori 0 continued. Table 3.

Article number, page 28 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 2 1 . . 84 71 43 44 1 0 . . . . 0 0 0 0 + − + − + − 2 . 0 95 08 . . 1 5 ) (Myr) 9 5 . . 11 10 29 20

2 2 . . . . 0 0 0 0 + − M + − + − 7 Mass Age . 95 82 . . 1 i ∗ V 47 . 1.54 8 3.01 2 3 3 2 3 α s d d 2 3 3 0 . 57 29 EW H 1 . . 0 0 α ± ± ± 2 . 10 80 21 . . 9 − 19 − − 2 − 044 054 051 039 m) (Å) line shape ( . . . . 0 0 0 0 µ 10 + − + − · 22 ) − 225 165 80 30 . . . . 0 0 4 0 0 . + − 65 . (0 2 − 060 074 057 045 m) (3 . . . . 0 0 0 0 µ 10 + − + − 4 · . ) 3 267 216 29 81 − . . . . 0 0 − + 24 . 57 . (1 NameWRAY 15-1435 (0 Near IR excess Mid IR excess H WW Vul 0 XY Per A 0 line profile classification is as follows: single-peaked (s), double-peaked (d) and showing a P-Cygni profile (P), both regular or inverse. An asterisk together with the variability α The H continued. ˇ cerová et al. (2013); (11) Hernández et al. (2005); (12) Dunkin et al. (1997); (13) Pogodin et al. (2012); (14) Miroshnichenko et al. (1999); (15) Sartori et al. (2010); (16) Polster et al. (2012); (17) Table 3. References. indicator indicates that the source hadet been al. catalogued (2011a); as UXOR (4) type AbabakrKu in et the literature. al. References (2016); for (5)Manoj EW et values Hernández and al. et line (2006); al. shapes: (18) (1) (2004);et Miroshnichenko Fairlamb al. (6) et et (2012); Baines al. al. (24) (2017); et (2004); Spezzi (2) al.(1998); (19) Carmona et (31) (2006); Miroshnichenko et al. X-Shooter al. et (7) (2008); (2010); spectra, al. Wheelwright (25) (3) 2015, (2002); et Houet Mendigutía priv. (20) et al. al. comm, Borges al. (2016); (2010); from Fernandes (2016); (36) (8) ESO et (26) Miroshnichenko van observing al. Grinin et den program (2007); & al. 084.C-0952A; Ancker (21) Rostopchina (1998); (32) et Vieira (1996); (37) Ababakr et al. (27) Miroshnichenko al. et (1998); Vieira et (2003); al. et (9) al. (22) (2017); al. (2000). Oudmaijer Boehm (33) (2011); & & Zuckerman (28) Drew Catala et Acke (1999); et (1995); al. (10) al. (23) (2008); (2005); Nakano (34) (29) Oudmaijer, Herbig et & al. Bell (1998); (1988); (35) (30) Frasca Oudmaijer, et al.

Article number, page 29 of 39 A&A proofs: manuscript no. Vioque_et_al_2018

Table 4. IR excess at each bandpass (defined as Fobserved/FCK) for each Herbig Ae/Be star belonging to the high quality sample of 218 sources.

Name J H Ks W1 W2 W3 W4 1.24µm 1.66µm 2.16µm 3.4µm 4.6µm 12µm 22µm AB Aur 2.79 6.01 13.23 28.77 - 240.49 2244.34 AK Sco 1.34 1.88 3.14 7.78 11.61 49.83 400.83 AS 310 0.68 0.71 0.81 - - - - AS 470 1.21 1.28 1.71 2.15 2.77 3.21 5.36 AS 477 2.00 - - 27.30 59.04 104.62 475.86 BD+30 549 1.81 2.08 2.35 2.79 3.39 51.04 378.04 BD+41 3731 0.95 0.93 0.92 0.94 1.03 2.54 35.38 BF Ori 1.50 2.39 4.56 10.69 21.87 147.30 617.52 BH Cep 1.23 1.79 3.24 7.10 10.97 52.96 585.75 BO Cep 1.39 1.76 2.27 4.22 5.86 15.29 319.78 CO Ori 1.94 2.52 4.28 6.55 13.63 36.94 172.89 CPM 25 2.82 7.14 19.80 73.01 215.11 2106.36 13608.54 CQ Tau 1.06 1.92 4.31 8.63 19.75 150.12 1840.10 DG Cir 8.79 15.43 30.63 66.60 174.10 1164.79 11702.50 GSC 1876-0892 1.69 3.44 8.57 22.43 59.77 650.69 6589.47 GSC 3975-0579 1.54 2.90 7.03 21.27 36.92 79.38 933.24 GSC 6546-3156 1.57 1.91 2.37 6.11 10.31 717.95 5359.18 GSC 8143-1225 1.75 2.55 4.77 9.99 15.83 26.31 361.41 GSC 8581-2002 1.96 1.97 2.12 1.98 2.03 7.17 76.96 GSC 8645-1401 1.38 2.20 4.18 9.16 15.58 69.63 719.30 GSC 8994-3902 1.14 1.18 1.28 - - - - HBC 217 1.32 1.63 2.42 5.27 8.16 33.06 1021.50 HBC 222 1.30 1.64 2.74 6.10 9.02 18.50 323.71 HBC 334 2.04 2.72 3.52 10.86 20.90 788.83 11504.97 HBC 442 1.29 1.40 1.91 3.40 5.54 69.45 550.58 HBC 7 1.22 1.26 1.55 1.61 2.34 3.56 8.26 HBC 705 1.45 1.55 1.94 1.79 2.64 4.36 16.23 HD 100453 1.34 2.00 4.26 9.24 24.07 134.69 1857.99 HD 100546 1.34 2.09 3.64 10.62 21.19 490.87 7717.64 HD 101412 1.60 2.34 4.88 15.60 44.31 209.20 1006.20 HD 104237 1.50 2.40 4.60 8.98 21.58 84.07 367.56 HD 114981 1.23 1.24 1.34 1.41 1.39 1.32 4.53 HD 130437 2.19 2.56 3.24 4.23 5.85 8.52 32.43 HD 132947 1.30 1.42 1.54 - - - - HD 135344 0.56 0.47 0.49 0.49 0.52 0.56 46.31 HD 135344B 1.38 2.05 4.04 8.55 17.29 16.84 420.92 HD 139614 1.39 1.76 3.16 6.59 12.20 150.84 2119.75 HD 141569 1.29 1.32 1.45 2.77 2.69 6.48 102.87 HD 141926 1.59 1.80 2.31 3.01 4.38 5.53 18.99 HD 142527 2.45 4.16 8.24 15.15 36.02 158.75 1216.10 HD 142666 1.64 2.50 4.63 9.18 20.00 170.41 908.47 HD 143006 1.22 1.50 2.56 7.36 14.36 24.41 365.00 HD 144432 1.86 2.85 5.34 8.29 17.05 146.13 837.23 HD 149914 1.30 1.33 1.42 1.36 1.50 1.06 5.65 HD 150193 1.90 3.19 6.08 12.60 28.58 194.31 980.17 HD 155448 - 0.94 0.94 - - - - HD 158643 0.97 1.17 1.80 3.31 - 26.38 72.09 HD 163296 1.81 3.36 7.07 17.49 37.20 202.40 853.49

Article number, page 30 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars

Table 4. continued.

Name J H Ks W1 W2 W3 W4 1.24µm 1.66µm 2.16µm 3.4µm 4.6µm 12µm 22µm HD 169142 1.28 1.72 2.71 2.86 4.77 48.62 1212.98 HD 17081 1.08 1.10 1.11 1.06 1.27 0.69 1.56 HD 174571 0.99 1.02 1.08 0.67 0.72 0.76 4.03 HD 176386 1.73 1.77 2.04 3.89 3.54 - 1703.78 HD 179218 1.04 1.37 2.54 6.64 16.19 289.90 1956.36 HD 199603 0.98 0.95 1.03 0.97 1.14 0.61 0.76 HD 200775 2.25 3.98 8.36 25.21 - 149.52 2914.79 HD 235495 3.04 6.46 14.38 36.37 74.80 332.15 1821.12 HD 244314 1.70 2.81 5.47 13.82 22.29 143.73 778.93 HD 244604 1.78 3.20 7.17 11.50 23.27 137.15 839.31 HD 245185 1.91 3.16 6.58 17.93 36.27 705.27 4958.59 HD 249879 1.31 1.92 5.02 23.39 59.15 474.67 2392.06 HD 250550 3.11 7.68 18.39 44.57 97.12 614.45 6073.51 HD 259431 2.08 4.10 9.61 27.60 94.90 305.32 2145.16 HD 287823 1.29 1.99 4.48 12.55 24.20 76.66 1144.68 HD 288012 0.63 0.59 0.64 - - - - HD 290380 1.41 2.12 3.72 7.95 12.64 51.05 701.84 HD 290409 1.67 2.42 4.04 4.53 6.53 219.21 3060.97 HD 290500 2.10 2.82 4.14 12.80 20.09 100.80 2143.13 HD 290764 1.71 2.95 6.34 15.76 28.69 70.06 1819.14 HD 290770 2.15 4.07 9.14 19.26 34.92 200.05 1406.43 HD 305298 1.46 1.65 1.70 2.63 4.76 271.40 12010.35 HD 313571 1.86 2.18 2.65 4.29 4.69 10.16 152.52 HD 31648 1.60 2.68 5.49 9.48 22.53 120.18 611.63 HD 319896 1.43 1.61 1.98 5.75 9.81 22.37 414.03 HD 323771 2.54 5.90 13.58 38.02 76.64 312.52 2185.92 HD 34282 1.82 3.78 8.31 17.85 28.81 104.86 1196.27 HD 344261 1.15 1.24 1.24 1.23 1.20 0.81 - HD 34700 1.00 1.02 1.24 - - - - HD 35187 1.90 2.74 4.65 7.72 13.40 59.77 569.45 HD 35929 1.35 1.47 1.97 3.06 5.72 12.01 29.45 HD 36112 1.63 2.71 5.61 14.70 31.22 99.40 992.76 HD 36408 1.06 1.01 1.02 0.79 0.81 0.74 8.02 HD 36917 1.77 2.20 2.98 8.99 20.56 84.71 264.73 HD 36982 1.80 2.02 2.38 - - - - HD 37258 1.97 3.27 6.62 15.84 28.99 189.17 970.64 HD 37357 1.54 2.36 4.22 7.84 14.86 102.39 696.70 HD 37371 0.35 0.36 0.35 0.32 0.40 0.77 7.26 HD 37490 0.88 1.02 1.19 2.35 4.17 3.60 7.30 HD 37806 2.20 4.90 10.98 22.93 60.67 293.60 1185.58 HD 38087 2.01 2.45 2.77 2.94 2.84 2.04 275.15 HD 38120 1.80 3.11 6.10 10.25 20.61 658.25 5227.02 HD 39014 0.98 1.15 1.17 1.26 1.75 0.82 1.00 HD 41511 4.36 8.56 12.99 - - 89.68 289.86 HD 45677 3.19 7.44 32.20 - - 4700.04 28602.62 HD 46060 0.79 0.75 0.74 0.80 0.77 3.47 171.16 HD 50083 1.23 1.34 1.70 2.17 3.60 4.74 22.96 HD 50138 2.13 4.29 10.45 26.11 - 444.08 2375.59 HD 56895B 1.12 1.07 1.19 - - - -

Article number, page 31 of 39 A&A proofs: manuscript no. Vioque_et_al_2018

Table 4. continued.

Name J H Ks W1 W2 W3 W4 1.24µm 1.66µm 2.16µm 3.4µm 4.6µm 12µm 22µm HD 58647 1.29 1.73 3.27 8.58 27.12 39.00 105.20 HD 59319 1.15 1.11 1.16 1.24 1.21 0.88 129.98 HD 68695 1.53 2.61 5.20 15.95 26.36 96.34 1314.61 HD 76534 1.31 1.31 1.40 2.21 2.55 3.90 29.05 HD 85567 1.87 3.74 8.55 22.94 58.45 209.67 836.94 HD 87403 1.35 1.36 1.49 1.67 1.76 1.15 10.21 HD 87643 3.76 13.14 41.40 - - 2278.72 12716.58 HD 94509 1.28 1.40 1.53 2.17 2.55 3.48 17.17 HD 95881 2.28 4.48 10.93 33.15 86.18 295.79 1007.97 HD 96042 1.47 1.45 1.71 1.48 1.49 1.66 45.46 HD 9672 1.04 1.00 1.12 1.05 1.22 0.79 3.97 HD 97048 1.84 3.01 5.88 19.83 20.04 239.06 3219.89 HD 98922 2.15 4.46 11.09 - - 352.13 1148.15 HR 5999 1.59 2.82 6.25 - - 112.50 414.85 HT CMa 4.29 9.27 20.90 64.36 148.82 983.40 5396.58 HU CMa 1.58 3.04 7.04 20.04 42.48 535.48 3456.52 Hen 3-1121 1.73 4.13 10.01 26.17 62.07 699.79 6435.45 Hen 3-1121S 0.87 0.82 0.85 0.87 1.24 29.55 1962.39 Hen 3-1191 4.48 21.19 98.70 510.27 1995.88 10315.68 37787.17 Hen 3-823 1.35 1.45 1.73 2.15 2.63 3.62 9.41 Hen 3-847 1.38 2.77 12.52 127.89 761.32 9538.41 58854.42 Hen 3-938 2.93 8.65 27.33 69.17 213.51 823.51 3658.69 IL Cep 0.85 0.84 0.89 0.92 1.12 1.12 - IP Per 1.80 3.18 6.85 15.40 27.39 36.94 645.90 KK Oph 3.04 12.33 41.94 102.96 243.95 1149.69 4781.29 LKHa 260 1.49 3.11 - 22.31 53.84 388.18 312.56 LKHa 338 4.59 11.05 38.14 193.55 635.30 6577.01 38467.51 LkHa 208 1.73 2.31 4.08 12.66 43.63 1204.97 7023.91 LkHa 215 2.25 4.00 7.74 27.52 53.10 183.08 1801.94 LkHa 257 1.55 2.63 4.83 8.08 13.28 58.14 1381.77 LkHa 259 2.90 3.65 5.33 10.63 24.83 580.94 10317.88 LkHa 324 - 1.06 1.13 1.61 2.73 73.32 653.45 LkHa 339 0.75 1.07 2.15 5.96 13.88 243.18 2014.27 MQ Cas 5.39 19.10 58.45 - - - - MWC 1021 1.54 1.76 2.77 - - 37.54 74.27 MWC 1080 3.17 8.57 20.05 - - 579.77 2359.22 MWC 137 1.94 3.25 8.25 30.58 86.52 328.65 3052.44 MWC 297 1.55 3.71 8.52 - - 280.89 4022.76 MWC 342 2.99 6.55 15.48 - - 717.92 3665.94 MWC 593 1.61 1.83 2.24 - - - - MWC 655 3.26 4.01 6.05 5.07 6.92 41.93 1084.04 MWC 657 1.41 2.98 6.76 18.91 66.61 242.44 800.51 MWC 878 1.61 2.72 7.30 43.34 157.49 976.37 2629.98 MWC 953 1.26 1.42 1.65 1.99 2.20 3.13 102.09 NSV 2968 2.31 4.35 10.22 12.67 41.90 297.56 2246.08 NV Ori 1.19 1.64 2.92 5.05 8.88 88.23 553.15 PDS 002 1.08 1.24 1.76 3.39 4.95 51.00 718.19 PDS 004 1.43 2.17 3.79 7.57 15.55 160.54 1084.87 PDS 021 1.93 4.07 8.97 18.07 41.18 595.37 4130.11

Article number, page 32 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars

Table 4. continued.

Name J H Ks W1 W2 W3 W4 1.24µm 1.66µm 2.16µm 3.4µm 4.6µm 12µm 22µm PDS 022 1.21 1.48 2.47 6.19 13.77 340.58 4104.11 PDS 025 3.07 6.73 15.36 37.43 63.06 373.77 3752.05 PDS 123 2.58 6.31 14.05 32.32 65.17 326.20 2513.33 PDS 124 1.69 2.51 4.95 15.08 30.21 456.06 2857.71 PDS 126 1.51 2.16 4.09 8.08 14.23 79.11 390.36 PDS 129 1.38 1.58 2.17 3.86 6.50 78.22 630.91 PDS 130 2.13 4.13 8.84 22.69 45.11 500.19 3967.51 PDS 133 2.13 6.27 16.14 23.84 77.26 645.47 6061.49 PDS 134 1.79 1.90 2.09 - - - - PDS 138 1.20 1.20 1.24 - - - - PDS 174 2.24 3.43 5.38 12.27 26.19 537.78 29319.27 PDS 211 1.35 2.16 4.29 7.32 18.45 433.58 3154.52 PDS 24 2.15 3.93 8.37 18.41 38.89 852.57 9838.27 PDS 241 1.30 1.51 1.86 7.04 19.17 5641.92 83493.60 PDS 27 1.78 4.54 12.10 31.59 134.49 1091.96 7522.03 PDS 277 1.23 1.54 2.80 6.42 10.08 70.34 1354.70 PDS 286 1.73 1.85 2.39 0.99 1.33 1.82 21.31 PDS 290 0.64 0.67 0.76 0.63 0.61 1.88 111.26 PDS 297 1.12 1.12 1.15 1.13 1.16 1.11 18.69 PDS 324 1.34 1.53 1.79 - - - - PDS 33 1.80 2.83 4.94 11.66 21.41 458.11 4167.59 PDS 34 2.61 5.93 14.05 31.75 63.65 605.75 7602.75 PDS 344 1.78 2.17 2.65 4.71 8.69 390.97 3434.61 PDS 361S 2.87 3.39 3.74 5.80 6.63 4.67 - PDS 37 0.71 2.38 6.48 18.45 55.76 526.99 4328.48 PDS 389 1.55 2.02 3.23 4.82 8.86 38.90 328.32 PDS 415N 0.92 1.14 1.62 4.08 9.07 21.21 332.90 PDS 431 1.41 1.41 1.55 1.50 1.54 31.35 - PDS 469 1.44 1.85 2.61 4.02 6.10 214.67 1185.71 PDS 477 3.47 9.52 24.09 61.85 129.93 1201.43 9937.68 PDS 520 4.96 8.63 18.99 41.85 94.93 429.65 3731.82 PDS 543 0.71 0.62 0.65 0.76 1.94 19.23 572.04 PDS 69 1.15 2.19 4.42 - - - - PX Vul 1.75 2.68 4.57 10.61 18.18 46.64 341.19 RR Tau 1.23 3.42 8.44 24.63 48.23 156.34 996.01 RY Ori 1.35 1.63 2.70 4.49 7.27 58.82 358.89 SAO 185668 0.87 0.85 0.84 0.82 1.01 25.56 881.06 SAO 220669 1.43 1.43 1.46 1.65 1.56 3.78 251.64 SV Cep 2.62 5.05 10.71 23.65 49.17 743.44 4567.70 T Ori 2.65 5.88 14.38 35.54 52.41 192.59 1528.40 TY CrA 1.59 2.16 2.66 5.32 14.92 3250.79 6107.38 UX Ori 2.89 4.96 10.85 23.75 46.70 338.18 2472.27 V1012 Ori 1.42 2.75 5.85 10.93 20.05 94.49 1839.69 V1295 Aql 1.30 2.07 4.38 11.48 29.27 118.54 396.22 V1478 Cyg 2.70 4.71 13.20 - - 202.01 860.06 V1493 Cyg 2.39 3.92 7.04 16.28 35.08 197.28 900.73 V1685 Cyg 1.67 3.51 7.74 - - - 3124.43 V1686 Cyg 3.11 9.36 28.85 - - - - V1787 Ori 2.03 3.14 6.49 15.63 33.22 236.07 1369.96

Article number, page 33 of 39 A&A proofs: manuscript no. Vioque_et_al_2018

Table 4. continued.

Name J H Ks W1 W2 W3 W4 1.24µm 1.66µm 2.16µm 3.4µm 4.6µm 12µm 22µm V1818 Ori 0.36 1.94 7.77 28.86 86.74 453.29 2698.23 V1977 Cyg 2.09 4.14 10.33 25.91 66.33 153.92 794.03 V2019 Cyg 1.36 1.88 3.05 6.09 10.78 - 6223.77 V346 Ori 1.14 1.69 3.15 9.45 15.76 44.81 776.76 V350 Ori 1.28 2.45 5.47 11.91 20.34 133.53 1164.99 V351 Ori 1.73 2.44 4.68 11.46 22.68 44.22 691.22 V361 Cep 1.08 1.20 1.42 2.00 2.93 6.86 1741.54 V373 Cep 3.30 8.69 21.53 139.91 732.16 2133.92 42216.18 V374 Cep 1.57 1.83 2.42 1.98 2.76 4.93 27.19 V380 Ori 2.31 5.36 12.67 27.79 74.38 387.31 1548.31 V388 Vel 1.34 3.55 9.33 40.19 94.40 1665.31 21691.47 V431 Sct 1.34 2.96 15.27 129.70 - 4848.99 34190.80 V594 Cas 3.39 8.94 21.62 49.46 113.26 522.05 2764.98 V594 Cyg 26.23 29.73 43.29 25.55 34.83 54.10 169.97 V599 Ori 1.70 2.62 5.27 6.94 12.85 41.44 959.30 V669 Cep 1.83 3.48 10.37 44.40 174.02 2080.63 9413.59 V718 Sco 1.72 2.24 3.81 7.66 16.83 155.15 707.61 V921 Sco 3.58 8.44 24.92 - - 3378.91 57901.14 VV Ser 3.57 8.83 21.97 46.08 121.07 428.03 1766.30 VX Cas 2.00 4.62 11.25 27.13 45.17 276.95 2735.55 WRAY 15-1435 2.24 2.88 3.34 4.39 7.36 192.05 3602.51 WW Vul 2.19 4.60 10.53 17.32 32.02 231.68 1366.23 XY Per A 2.61 4.51 9.33 21.53 47.11 131.33 670.35

Article number, page 34 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 2 4 4 2 5 1 20 31 29 21 22 V Binary 7.36 Yes 9.50 Yes 9.80 Yes 9.63 Yes 9.25 Yes 11.5318.54 - - 13.91 Yes 15.52 - 15.12 - 16.76 - 14.94 - 17.00 - 13.55 - 10.55 Yes 12.79 - 10.92 Yes 11.51 - 12.79 - 13.07 Yes 11.97 - 13.46 - 15.55 - 13.41 - 13.54 - 16.60 - 11.75 - 17.46 - 11.85 - 12.79 - 12.60 Yes 13.10 - 15.17 Yes 12.97 - V 25 24 11 14 50 50 23 15 14 18 19 16 25 25 50 50 27 20 16 25 13 21 16 18 50 50 19 17 18 31 12 12 29 23 14 18 25 29 19 16 50 50 26 29 19 16 19 16 12 12 50 50 17 45 12 16 ...... 043 050 070 000 070 080 050 030 020 000 040 050 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 . . . . 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 83 12 60 18 24 05 51 93 82 85 51 97 50 77 72 13 39 16 44 53 90 63 12 13 51 87 68 37 ...... 0 8 1 3 5 1 2 1 2 0 0 2 4 3 8 2 1 5 1 1 2 2 5 2 4 4 5 3 051 000 570 870 110 030 ...... 2 0 0 1 1 1 ) (mag) (mag) 47 32 47 35 58 43 32 23 46 34 23 17 53 38 49 37 52 37 28 25 39 39 30 22 44 32 60 47 31 26 15 21 27 23 33 24 41 43 62 56 11 09 19 15 26 21 48 34 33 28 64 68 51 29 34 23 20 16 27 30 63 39 27 29

...... 057 057 072 059 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 . . 0 0 0 0 0 0 0 0 0 0 L 0 0 0 0 − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + + − + − 47 73 01 52 03 65 80 50 88 66 13 85 12 29 58 85 46 70 88 32 98 43 35 93 89 00 06 38 60 13 56 25 Log(L) A ...... 0 0 1 1 3 0 0 0 1 1 3 1 2 5 4 5 2 1 2 2 0 0 1 0 2 0 0 1 3 0 3 2 673 394 . . − − − − 0 0 − ff 80 70 e 180 160 160 220 180 160 150 150 160 160 250 250 250 250 250 250 250 250 120 250 180 160 180 160 180 160 250 250 500 500 130 130 + − 800 1000 3000 800 2000 3900 1000 1000 1000 1000 3000 800 1000 1000 1700 1400 5000 3000 1000 1000 2000 3900 3900 4500 4500 5000 1000 1000 2000 2000 1500 800 4000 4000 + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 5350 7250 8150 7830 8150 6400 7990 8750 7000 7750 7000 8150 8150 8150 9750 8500 12500 15000 16500 32900 29500 14000 16500 15800 11900 19500 12500 32900 29000 24500 12500 36900 11500 12000 2 6 4 8 5 7 ...... 30 17 40 22 58 42 12 11 13 9 18 16 47 28 22 19 4 4 3 3 2 2 520 270 160 100 440 190 470 220 370 150 390 220 400 170 430 170 390 170 550 370 570 390 650 420 630 380 470 220 610 360 660 420 470 260 350 150 170 100 410 180 580 330 150 50 110 50 + − + − + − + − + − + − + − + − − + − + − + + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − + − 6 4 5 . . . 37 243 39 350 86 980 11 850 52 760 47 900 43 670 12 1390 48 690 79 130 83 600 38 750 65 1670 12 149 53 880 27 1730 12 2430 21 101 57 176 17 1390 12 496 38 690 70 95 13 820 40 790 12 117 38 1470 89 230 ...... 034 605 042 2980 036 3280 094 2590 059 343 082 355 0 0 . 0 0 0 0 0 0 0 0 0 0 . . . 0 0 0 0 0 0 0 0 0 0 . 0 . 0 0 0 0 0 0 0 0 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 10 79 66 12 16 27 90 49 67 77 03 71 84 69 52 30 39 87 70 27 99 96 54 14 53 52 94 30 ...... 4 0 9 0 1 0 1 649 191 173 162 910 811 ...... − − − − − − − 0 − 65:38:38 0 65:39:53 42:14:02 0 44:17:16 0 26:19:55 0 14:52:57 0 31:06:20 0 08:08:21 0 04:04:35 1 29:32:50 0 67:57:39 2 09:48:02 1 50:14:21 0 28:19:15 8 10:26:05 + + + + + + + + + + + + + + + (h:m:s) (deg:m:s) (mas) (pc) (K) ( Be star belonging to the low quality sample of 34 sources. / DK ChaGSC 5360-1033 05:57:49.5 -14:05:34 12:53:17.1 1 -77:07:11 4 NameBP Psc 23:22:24.7 RA -02:13:42 DEC 2 Parallax Distance T GSC 5988-2257 07:41:41.1 -20:00:13 GSC 6542-2339 07:24:37.0 -24:34:47 1 HBC 1 00:07:02.6 HBC 324 00:07:30.7 HBC 694 20:24:29.5 HBC 717 20:52:06.0 HD 245906 05:39:30.5 HD 53367 07:04:25.5 -10:27:16 7 HD 72106B 08:29:34.9 -38:36:21 0 Hen 2-80 12:22:23.2 -63:17:17 0 MWC 314 19:21:34.0 MWC 623 19:56:31.5 MWC 930 18:26:25.2 -07:13:18 NX Pup 07:19:28.3 -44:35:11 PDS 144S 15:49:15.3 -26:00:55 6 PDS 229N 06:55:40.0 -03:09:50 PDS 322 10:52:08.7 -56:12:07 PDS 364 13:20:03.6 -62:23:54 PDS 371 13:47:31.4 -36:39:50 9 PDS 453 17:20:56.1 -26:03:31 5 PDS 530 18:41:34.4 PDS 551 18:55:23.0 PDS 581 19:36:18.9 PV Cep 20:45:54.0 R CrA 19:01:53.7 -36:57:09 10 UY Ori 05:32:00.3 -04:55:54 2 V590 Mon 06:40:44.6 V645 Cyg 21:39:58.3 V892 Tau 04:18:40.6 VY Mon 06:31:06.9 Z CMa 07:03:43.2 -11:33:06 4 Main parameters of each Herbig Ae Table 5.

Article number, page 35 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 , Mass and Age. The references for binarity are: (1) Baines et al. (2006); (2) Wheelwright et al. (2010); L and V taken from the following sources in order of choice: Fairlamb et al. 2015; Montesinos et al. 2009; Hernández et al. 2004; Mendigutía et al. et al. (2017); (19) Mayer et al. (2016); (20) Folsom et al. (2008); (21) Corporon & Lagrange (1999); (22) Doering & Meixner (2009); (23) Chelli et V ff ,A ff e Atmospheric parameters T 2012; Carmona et al. 2010;If Chen not et al. available 2016; they Alecian were et derived al. as 2013; described Sartori et in al. Sect. 2010; 2.2. Manoj See et Sect. al. 2006; 3 Hernández for et derivation al. of 2005; Vieira et al. 2003; APASS Data Release 9 and the SIMBAD database. References. (3) Leinert et al. (1997);(10) (4) Liu Maheswar et et al. al.(17) (2000); (2002); Morrell (11) (5) & Biller Wheelwright Levato et etal. (1991); al. al. (1995); (18) (2011); (2012); (24) Lazare (6) (12) Miroshnichenko Alecian et(2009); Schütz et (31) al. et al. Frasca (2002); al. (2013); et (25) (7) (2011); al. Friedemann Hamaguchi (13) (2016); et et (32) Boersma al. al. Marston et (1996); (2008); & al. (26) (8) McCollum (2009); Kraus (2008); Dunhill (14) (33) et et Zhang Malkov al. al. et et (2015); (2012); al. (9) (27) al. (2017). Coulson Torres (2006); et & (15) al. Walther Arellano (1995); (2000); Ferro (28) & Aspin Giridhar (1998); (29) (2003); Connelley (16) et Kubát al. et (2008); al. (30) (2010); Millour et al.

Article number, page 36 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars 2 1 3 4 4 9 8 3 6 0 8 6 7 3 7 8 6 5 3 ...... 49 07 83 59 14 1 20 13 3 17 75 75 52 50 − 6 7 . 1 1 2 3 . . . 8 1 . . 2 1 2 1 2 1 5 2 0 . 0 . . . 010 005 + − + − 0 0 0 0 0 0 . . 0 0 0 0 + − + − + − + + + − − − + − + − + − 0 0 6 5 + − + − + − + − + − 10 2 1 2 7 2 2 8 1 6 + − ...... · 5 1 1 2 2 2 0 08 80 88 28 45 ) . . . . . 13 17 0 0 0 0 1 015 80 50 . . . 0 0 0 + − 00 . (10 ) (Myr) 5 0 0 8 0 2 2 5 1 5 5 7 7 1 9 3 4 7 6 5 ...... 29 50 07 26 90 39 52 23 13 11 47 13 90 50

7 4 2 0 1 1 . . 1 0 1 0 4 3 . . 2 0 . . 1 0 2 1 ...... 1 0 068 068 0 0 0 0 0 0 0 0 0 0 0 0 . 0 0 . + − + − + − + − + − + − + − + − + − + − 0 0 M + − + + − − + − + − + − + − 8 8 4 6 4 7 2 4 0 4 + − Mass Age ...... 07 89 30 30 51 50 90 ...... 369 . i - 3 - 24 --- - 3 ---- 2 - 18 - 1 - 5 ---- 2 - 2 ------3 ---- 2 ------* - - -* 8 -* 5 V ------5 - - 1 ------1 - - 1 - - 3 8 2 8 α 21 21 31 21 36 28 21 21 16 21 21 31 31 26 21 21 28 26 21 31 21 21 d d P s s s s s s d d d d d d d d d d d P P P P P ----- 5 1 5 5 1 1 2 1 1 1 5 1 1 1 1 15 27 24 15 29 27 33 15 17 35 16 15 15 15 15 28 17 0 0 6 5 1 6 0 2 5 3 ...... 4 2 1 4 5 5 3 0 89 70 18 84 95 ...... EW H . 2 . 7 . 1 7 1 . 1 . 3 1 1 1 12 15 14 10 20 75 35 20 4 1 3 . . . 1 . 2 2 4 2 0 0 0 0 0 020 α 0 . 0 0 0 ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± 0 ± ± ± ± ± 3 0 8 5 7 3 0 0 2 2 ± ± ± ± ...... 2 2 1 7 0 3 0 2 ...... ± 5 2 16 47 59 80 00 125 130 134 201 . 88 . 26 25 . 69 . . 10 54 29 83 83 36 42 − − 37 . 51 59 94 43 95 28 40 . . . − − − − 155 − − − − − − − 9 1 − − − − − − − − 24 23 12 14 26 792 − . 19 14 − − − − − − − 3 − − 2 0 5 6 1 9 1 2 3 0 4 6 7 1 9 ...... 20 79 25 25 12 8 27 17 26 25 31 24 58 51 27 21 11 11 23 15 24 15 15 11 43 14 − − − . . . . 3 2 1 1 ...... 1 0 6 5 . . . . 1 0 3 4 . . 046 078 018 031 024 072 043 058 016 049 025 057 009 011 049 068 027 015 m) (Å) line shape ( + − + − 0 0 0 . . 0 0 0 0 0 0 0 ...... 0 0 0 0 0 . 0 . . 0 0 . . . 0 . . 0 0 . 0 . + − + − + − + − − + + − 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ − + + − + − + − − + − + + − − − + + + − + − 10 10 10 6 9 9 8 7 0 − + − + − + − − + + − − + + − + + − ...... 23 · · · 2 3 0 6 22 237 11 32 81 44 43 51 53 37 25 25 73 ) ) ) ...... 24 12 1 0 0 0 0 0 0 1 0 1 0 − 119 091 114 057 930 029 227 057 147 17 36 32 18 11 08 ...... 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 . + − + − − + 27 42 66 . . . (0 (0 (0 5 2 6 9 2 3 1 2 2 0 2 2 ...... 19 17 37 10 19 68 41 24 12 12 11 07 10 10 45 18 11 15 22 15 34 29 − − − − ...... 1 0 . . 7 5 . . 1 1 ...... 1 . 2 023 026 062 053 028 037 021 012 066 037 033 020 070 045 037 028 042 028 087 030 018 021 012 010 m) (3 + − . . 0 0 0 0 0 0 0 0 ...... 0 0 . . . . 0 0 0 0 . . . . 0 0 0 0 0 . . 0 + − + − + − − + 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 µ + − + − + − + − + − + − + − + − − + − + 10 10 10 10 0 7 6 7 + − + − + − + − + − + − + − + − − + − + + − + − . . . . 4 · · · · . 99 17 32 25 47 34 78 50 88 46 ) ) ) ) 3 ...... 017 145 077 146 131 084 061 231 044 104 124 077 16 33 39 28 27 21 34 23 − ...... 0 0 0 0 0 0 0 0 + − + − + − − + 24 . 97 79 23 53 . . . . (1 Be star belonging to the low quality sample of 34 sources. / HBC 324 0 GSC 6542-2339HBC 1 - 83 - HBC 694 0 VY Mon 0 HBC 717 0 DK ChaGSC 5360-1033GSC 5988-2257 0 8 - - Z CMa 0 NameBP Psc Near IR excess Mid IR excess 0 H V892 Tau 1 MWC 314 0 V645 Cyg (0 HD 72106B 0 V590 Mon 0 HD 245906HD 53367Hen 2-80 (0 0 0 PDS 530PDS 581R CrAUY Ori 0 0 0 16 MWC 623 0 PDS 551PV Cep 0 4 MWC 930NX PupPDS 322 0 (0 ------PDS 229NPDS 364PDS 453 0 0 0 PDS 144S 2 PDS 371 (0 Other parameters of each Herbig Ae Table 6.

Article number, page 37 of 39 A&A proofs: manuscript no. Vioque_et_al_2018 ) for these objects could not be derived as they are not astrometrically well i V line profile classification is as follows: single-peaked (s), double-peaked (d) and showing a P-Cygni profile (P), both regular or inverse. An asterisk together with the variability α ˇ cerová et al. (2013); (11) Hernández et al. (2005); (12) Dunkin et al. (1997); (13) Pogodin et al. (2012); (14) Miroshnichenko et al. (1999); (15) Sartori et al. (2010); (16) Polster The H References. indicator indicates that thebehaved. source Similarly, had many of been these cataloguedmasses sources as and fall UXOR ages outside type in the the in Pre-Main(3) cases Sequence the Mendigutía they tracks et literature. were al. and Variability computable (2011a); isochrones indicator(1999); but (4) in values (10) these Ababakr the ( Ku et values HR al. have diagram to (2016); and (5) be no Hernández taken masses et with or al. caution. ages (2004); References could (6) for be Baines EW derived et values for al. and them. (2006); line We (7) decided shapes: Wheelwright to (1) et present Fairlamb al. the et (2010); al. (8) (2017); van den (2) Ancker Carmona et et al. al. (1998); (2010); (9) Oudmaijer & Drew et al. (2012); (17)(1995); Manoj (23) et Nakano al. et (2006); al.(30) (18) (2012); Oudmaijer, Miroshnichenko (24) et et Spezzi al. al. et (1998);al. al. (2004); (31) (1998); (2008); (19) X-Shooter (35) (25) Miroshnichenko spectra, Frasca Hou et 2015, et et al. priv. al. al. comm, (2016); (2002); (2016); from (36) (20) (26) ESO Miroshnichenko Borges Grinin observing et Fernandes & al. program Rostopchina et (1998); 084.C-0952A; (1996); al. (37) (32) (27) Miroshnichenko (2007); Ababakr Vieira et (21) et et al. al. Vieira al. (2000). et (2017); (2011); al. (33) (28) Zuckerman (2003); Acke et et (22) al. al. Boehm (2008); (2005); & (29) (34) Catala Herbig Oudmaijer, et & Bell (1988);

Article number, page 38 of 39 M. Vioque et al.: Gaia DR2 study of Herbig Ae/Be stars

Table 7. IR excess at each bandpass (defined as Fobserved/FCK) for each Herbig Ae/Be star belonging to the low quality sample of 34 sources.

Name J H Ks W1 W2 W3 W4 1.24µm 1.66µm 2.16µm 3.4µm 4.6µm 12µm 22µm BP Psc 1.49 2.28 4.20 10.16 22.35 168.78 2200.40 DK Cha 14.73 49.13 170.17 325.14 - 7497.91 50750.83 GSC 5360-1033 2.65 6.78 12.87 - - - - GSC 5988-2257 2.62 6.21 16.62 47.97 117.47 1245.27 6594.26 GSC 6542-2339 2.19 3.64 5.68 - - - - HBC 1 88.74 484.43 2157.04 9200.72 27745.90 156330.65 1383272.17 HBC 324 0.93 1.57 3.60 15.52 54.31 943.29 13361.89 HBC 694 3.40 8.51 21.52 103.52 255.39 1995.47 31718.13 HBC 717 1.49 2.35 4.58 12.21 23.09 149.10 1086.81 HD 245906 1.65 2.56 4.15 6.60 9.67 58.51 465.93 HD 53367 1.27 1.32 1.39 3.04 4.84 5.37 116.70 HD 72106B 1.79 1.92 2.73 5.48 10.88 138.92 976.86 Hen 2-80 2.60 6.89 17.96 64.65 211.29 2910.69 15743.40 MWC 314 2.32 2.65 3.51 2.91 6.01 7.36 23.11 MWC 623 4.04 7.46 12.17 28.53 78.51 193.41 397.67 MWC 930 0.35 0.33 0.35 0.48 0.86 0.74 2.16 NX Pup 1.51 4.32 13.40 50.79 137.41 359.52 1521.35 PDS 144S 3.81 14.93 58.24 317.57 681.57 2883.63 21943.52 PDS 229N 1.36 - - 11.85 22.39 173.19 938.57 PDS 322 1.08 1.24 1.58 3.81 8.96 181.89 22064.34 PDS 364 1.95 4.35 10.24 25.65 54.14 737.90 10364.41 PDS 371 1.03 3.51 9.47 22.73 55.17 247.26 2992.30 PDS 453 2.35 3.96 7.79 22.88 44.93 279.44 4685.85 PDS 530 0.86 3.03 13.43 62.07 146.37 960.03 7986.35 PDS 551 32.07 111.81 327.37 1230.71 4108.60 29748.04 220447.30 PDS 581 1.71 9.18 41.39 207.15 1047.80 4498.40 47015.82 PV Cep 2.83 24.69 147.34 483.75 2383.40 21362.00 157219.18 R CrA 18.56 89.19 567.83 - - 12353.87 184275.03 UY Ori 1.70 3.77 11.11 36.23 84.99 1932.72 12539.63 V590 Mon 1.88 4.42 12.28 44.17 104.49 2075.86 21654.33 V645 Cyg 0.74 2.67 19.54 - - 12523.05 124226.58 V892 Tau 18.38 57.37 142.10 149.55 392.45 14054.94 136211.22 VY Mon 2.70 5.67 13.53 38.47 - 897.71 7399.24 Z CMa 1.21 2.85 9.21 - - 373.86 2204.64

Article number, page 39 of 39