• We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.
• ### Open clusters: III. Fundamental parameters of B stars in NGC 6087, NGC 6250, NGC 6383 and NGC 6530. B type stars with circumstellar envelopes(1711.01311)

Nov. 3, 2017 astro-ph.SR
Context. Stellar physical properties of star clusters are poorly known. Aims. Our goals are to perform a spectrophotometric study of the B star population in open clusters to derive accurate stellar parameters, search for the presence of circumstellar envelopes, and discuss the characteristics of these stars. Methods. The BCD spectrophotometric system is a powerful tool to obtain fundamental parameters and infer the main properties of open clusters: distance modulus, color excess, and age. We inspected the Balmer discontinuity to seek circumstellar disks and identify Be-star candidates. High resolution spectra in the H$\alpha$ region are used to confirm the Be nature. Results. We provide Teff , log g, Mv , Mbol and spectral types for a sample of 68 stars in the field of the open clusters NGC 6087, NGC 6250, NGC 6383, and NGC 6530, as well as the cluster distances, ages and reddening. Then, based on a sample of 230 B stars in the direction of the 11 open clusters studied along this series of three papers, we report 6 new Be stars, 4 blue straggler candidates, and 15 B-type stars (called Bdd) with a double Balmer discontinuity. We also find that the majority of the Be stars are dwarfs and present a maximum at the spectral type B2-B4 in young and intermediate-age open clusters. Another maximum of Be stars is observed at the spectral type B6-B8 in open clusters older than 40 Myr, where the population of Bdd stars also becomes relevant. Conclusions. Our results support previous statements that the Be phenomenon is present along the whole main sequence band and occurs in very different evolutionary states. We find clear evidence of an increase of stars with circumstellar envelopes with cluster age. The Be phenomenon reaches its maximum in clusters of intermediate age and the number of B stars with circumstellar envelopes (Be+Bdd stars) is also high for the older clusters.
• ### Gravity darkening in stars with surface differential rotation(1705.06485)

Aug. 18, 2017 astro-ph.SR
The interpretation of stellar apparent fundamental parameters must be treated consistently with the characteristics of their surface rotation law. We develop a model to determine the distribution of the effective temperature and gravity, which depend on the surface differential rotation law and on the stellar external geometry. The basic assumptions are: a) the atmosphric layers are in radiative equilibrium; b) the bolometric flux is anti-parallel with the effective gravity; c) the angular velocity in the surface is Omega=Omega(theta,alpha,k)] where (alpha,k) are free parameters. The effective temperature varies with co-latitude theta' and on the differential-rotation law through the effective gravity and the gravity-darkening function (GDF). Although the derived expressions can be treated numerically, for some low integer values of k, analytical forms of the integral ofcharacteristic curves', on which the determination of the GDF relies, are obtained. The effects of the quantities (eta,alpha,k) (eta=ratio between centrifugal and gravitational accelerations at the equator) on the determination of the Vsini parameter and on the `gravity-darkening exponent' are studied. Depending on the values of (eta,alpha,k) the velocity V in the derived Vsini may strongly deviate from the equatorial rotational velocity. It is shown that the von Zeipel's-like gravity-darkening exponent $\beta_1$ depends on all parameters $(\eta,\alpha,k)$ and that its value also depends on the viewing-angle i. Hence, there no unique interpretation of this exponent determined empirically in terms of (i,alpha). The data on rotating stars should be analyzed by taking into account the rotational effects through the GDF, by assuming k=2 as a first approximation. Instead of the classical pair (eta,beta1), it would be more useful to determine the quantities (eta,alpha,i) to characterize stellar rotation.
• ### Critical study of the distribution of rotational velocities of Be stars; II: Differential rotation and some hidden effects interfering with the interpretation of the Vsin i parameter(1702.07684)

Feb. 24, 2017 astro-ph.SR
We assume that stars may undergo surface differential rotation to study its impact on the interpretation of $V\!\sin i$ and on the observed distribution $\Phi(u)$ of ratios of true rotational velocities $u=V/V_\rm c$ ($V_\rm c$ is the equatorial critical velocity). We discuss some phenomena affecting the formation of spectral lines and their broadening, which can obliterate the information carried by $V\!\sin i$ concerning the actual stellar rotation. We studied the line broadening produced by several differential rotational laws, but adopted Maunder's expression $\Omega(\theta)=\Omega_o(1+\alpha\cos^2\theta)$ as an attempt to account for all of these laws with the lowest possible number of free parameters. We studied the effect of the differential rotation parameter $\alpha$ on the measured $V\!\sin i$ parameter and on the distribution $\Phi(u)$ of ratios $u=V/V_\rm c$. We conclude that the inferred $V\!\sin i$ is smaller than implied by the actual equatorial linear rotation velocity $V_\rm eq$ if the stars rotate with $\alpha<0$, but is larger if the stars have $\alpha>0$. For a given $|\alpha|$ the deviations of $V\!\sin i$ are larger when $\alpha<0$. If the studied Be stars have on average $\alpha<0$, the number of rotators with $V_\rm eq\simeq0.9V_\rm c$ is larger than expected from the observed distribution $\Phi(u)$; if these stars have on average $\alpha>0$, this number is lower than expected. We discuss seven phenomena that contribute either to narrow or broaden spectral lines, which blur the information on the rotation carried by $V\!\sin i$ and, in particular, to decide whether the Be phenomenon mostly rely on the critical rotation. We show that two-dimensional radiation transfer calculations are needed in rapid rotators to diagnose the stellar rotation more reliably.
• ### A test field for Gaia. Radial velocity catalogue of stars in the South Ecliptic Pole(1609.08044)

Sept. 26, 2016 astro-ph.SR
Gaia is a space mission currently measuring the five astrometric parameters as well as spectrophotometry of at least 1 billion stars to G = 20.7 mag with unprecedented precision. The sixth parameter in phase space (radial velocity) is also measured thanks to medium-resolution spectroscopy being obtained for the 150 million brightest stars. During the commissioning phase, two fields, one around each ecliptic pole, have been repeatedly observed to assess and to improve the overall satellite performances as well as the associated reduction and analysis software. A ground-based photometric and spectroscopic survey was therefore initiated in 2007, and is still running in order to gather as much information as possible about the stars in these fields. This work is of particular interest to the validation of the Radial Velocity Spectrometer (RVS) outputs. The paper presents the radial velocity measurements performed for the Southern targets in the 12 - 17 R magnitude range on high- to mid-resolution spectra obtained with the GIRAFFE and UVES spectrographs.
• ### Physical parameters of IPHAS-selected classical Be stars. (I. Determination procedure and evaluation of the results.)(1604.05935)

April 20, 2016 astro-ph.SR
We present a semi-automatic procedure to obtain fundamental physical parameters and distances of classical Be (CBe) stars, based on the Barbier-Chalonge-Divan (BCD) spectrophotometric system. Our aim is to apply this procedure to a large sample of CBe stars detected by the IPHAS photometric survey, to determine their fundamental physical parameters and to explore their suitability as galactic structure tracers. In this paper we describe the methodology used and the validation of the procedure by comparing our results with those obtained from different independent astrophysical techniques for subsamples of stars in common with other studies. We also present a test case study of the galactic structure in the direction of the Perseus Galactic Arm, in order to compare our results with others recently obtained with different techniques and the same sample of stars. We did not find any significant clustering of stars at the expected positions of the Perseus and Outer Galactic Arms, in agreement with previous studies in the same area that we used for verification.
• ### The Gaia astrophysical parameters inference system (Apsis). Pre-launch description(1309.2157)

The Gaia satellite will survey the entire celestial sphere down to 20th magnitude, obtaining astrometry, photometry, and low resolution spectrophotometry on one billion astronomical sources, plus radial velocities for over one hundred million stars. Its main objective is to take a census of the stellar content of our Galaxy, with the goal of revealing its formation and evolution. Gaia's unique feature is the measurement of parallaxes and proper motions with hitherto unparalleled accuracy for many objects. As a survey, the physical properties of most of these objects are unknown. Here we describe the data analysis system put together by the Gaia consortium to classify these objects and to infer their astrophysical properties using the satellite's data. This system covers single stars, (unresolved) binary stars, quasars, and galaxies, all covering a wide parameter space. Multiple methods are used for many types of stars, producing multiple results for the end user according to different models and assumptions. Prior to its application to real Gaia data the accuracy of these methods cannot be assessed definitively. But as an example of the current performance, we can attain internal accuracies (RMS residuals) on F,G,K,M dwarfs and giants at G=15 (V=15-17) for a wide range of metallicites and interstellar extinctions of around 100K in effective temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]), and 0.25dex in surface gravity (logg). The accuracy is a strong function of the parameters themselves, varying by a factor of more than two up or down over this parameter range. After its launch in November 2013, Gaia will nominally observe for five years, during which the system we describe will continue to evolve in light of experience with the real data.
• ### Rotational velocities of A-type stars IV. Evolution of rotational velocities(1201.2052)

Jan. 10, 2012 astro-ph.SR
In previous works of this series, we have shown that late B- and early A-type stars have genuine bimodal distributions of rotational velocities and that late A-type stars lack slow rotators. The distributions of the surface angular velocity ratio \Omega/\Omega_crit (\Omega_crit is the critical angular velocity) have peculiar shapes according to spectral type groups, which can be caused by evolutionary properties. We aim to review the properties of these rotational velocity distributions in some detail as a function of stellar mass and age. We have gathered v sin i for a sample of 2014 B6- to F2-type stars. We have determined the masses and ages for these objects with stellar evolution models. The (Teff, log L/Lsun)-parameters were determined from the uvby-\beta photometry and the HIPPARCOS parallaxes. The velocity distributions show two regimes that depend on the stellar mass. Stars less massive than 2.5 Msun have a unimodal equatorial velocity distribution and show a monotonical acceleration with age on the main sequence (MS). Stars more massive have a bimodal equatorial velocity distribution. Contrarily to theoretical predictions, the equatorial velocities of stars from about 1.7 Msun to 3.2 Msun undergo a strong acceleration in the first third of the MS evolutionary phase, while in the last third of the MS they evolve roughly as if there were no angular momentum redistribution in the external stellar layers. The studied stars might start in the ZAMS not necessarily as rigid rotators, but with a total angular momentum lower than the critical one of rigid rotators. The stars seem to evolve as differential rotators all the way of their MS life span and the variation of the observed rotational velocities proceeds with characteristic time scales \delta(t)\sim 0.2 t_MS, where t_MS is the time spent by a star in the MS.
• ### A 10-hour period revealed in optical spectra of the highly variable WN8 Wolf-Rayet star WR 123(1104.5182)

April 27, 2011 astro-ph.SR
Aims. What is the origin of the large-amplitude variability in Wolf-Rayet WN8 stars in general and WR123 in particular? A dedicated spectroscopic campaign targets the ten-hour period previously found in the high-precision photometric data obtained by the MOST satellite. Methods. In June-August 2003 we obtained a series of high signal-to-noise, mid-resolution spectra from several sites in the {\lambda}{\lambda} 4000 - 6940 A^{\circ} domain. We also followed the star with occasional broadband (Johnson V) photometry. The acquired spectroscopy allowed a detailed study of spectral variability on timescales from \sim 5 minutes to months. Results. We find that all observed spectral lines of a given chemical element tend to show similar variations and that there is a good correlation between the lines of different elements, without any significant time delays, save the strong absorption components of the Hei lines, which tend to vary differently from the emission parts. We find a single sustained periodicity, P \sim 9.8 h, which is likely related to the relatively stable pulsations found in MOST photometry obtained one year later. In addition, seemingly stochastic, large-amplitude variations are also seen in all spectral lines on timescales of several hours to several days.
• ### The WR 140 periastron passage 2009: first results from MONS and other optical sources(1101.1430)

Jan. 7, 2011 astro-ph.SR
We present the results from the spectroscopic follow-up of WR140 (WC7 + O4-5) during its last periastron passage in January 2009. This object is known as the archetype of colliding wind binaries and has a relatively large period (~ 8 years) and eccentricity (~ 0.89). We provide updated values for the orbital parameters, new estimates for the WR and O star masses and new constraints on the mass-loss rates.
• ### Differential rotation in rapidly rotating early-type stars. I. Motivations for combined spectroscopic and interferometric studies(1012.1707)

Dec. 8, 2010 astro-ph.SR
Since the external regions of the envelopes of rapidly rotating early-type stars are unstable to convection, a coupling may exist between the convection and the internal rotation. We explore what can be learned from spectroscopic and interferometric observations about the properties of the rotation law in the external layers of these objects. Using simple relations between the entropy and specific rotational quantities, some of which are found to be efficient at accounting for the solar differential rotation in the convective region, we derived analytical solutions that represent possible differential rotations in the envelope of early-type stars. A surface latitudinal differential rotation may not only be an external imprint of the inner rotation, but induces changes in the stellar geometry, the gravitational darkening, the aspect of spectral line profiles, and the emitted spectral energy distribution. By studying the equation of the surface of stars with non-conservative rotation laws, we conclude that objects undergo geometrical deformations that are a function of the latitudinal differential rotation able to be scrutinized both spectroscopically and by interferometry. The combination of Fourier analysis of spectral lines with model atmospheres provides independent estimates of the surface latitudinal differential rotation and the inclination angle. Models of stars at different evolutionary stages rotating with internal conservative rotation laws were calculated to show that the Roche approximation can be safely used to account for the gravitational potential. The surface temperature gradient in rapid rotators induce an acceleration to the surface angular velocity. A non-zero differential rotation parameter may indicate that the rotation is neither rigid nor shellular underneath the stellar surface.
• ### A semi-automatic procedure for abundance determination of A- and F-type stars(0904.0105)

April 1, 2009 astro-ph.SR, astro-ph.IM
A variety of physical processes leading to different types of pulsations and chemical compositions is observed among A- and F-type stars. To investigate the underlying mechanisms responsible for these processes in stars with similar locations in the H-R diagram, an accurate abundance determination is needed, among others. Here, we describe a semi-automatic procedure developed to determine chemical abundances of various elements ranging from helium to mercury for this type of stars. We test our procedure on synthetic spectra, demonstrating that our procedure provides abundances consistent with the input values, even when the stellar parameters are offset by reasonable observational errors. For a fast-rotating star such as Vega, our analysis is consistent with those carried out with other plane-parallel model atmospheres. Simulations show that the offsets from the input abundances increase for stars with low inclination angle of about 4 degrees. For this inclination angle, we also show that the distribution of the iron abundance found in different regions is bimodal. Furthermore, the effect of rapid rotation can be seen in the peculiar behaviour of the H_beta line.
• ### Fundamental parameters of B Supergiants from the BCD System(0903.5134)

March 30, 2009 astro-ph.SR
Effective temperatures of early-type supergiants are important to test stellar atmosphere- and internal structure-models of massive and intermediate mass objects at different evolutionary phases. However, these Teff values are more or less discrepant depending on the method used to determine them. We aim to obtain a new calibration of the Teff parameter for early-type supergiants as a function of observational quantities that are highly sensitive to the ionization balance in the photosphere and its gas pressure, independent of the interstellar extinction, and as much as possible model-independent. The observational quantities that best address our aims are the (lambda1, D) parameters of the BCD spectrophotometric system. They describe the energy distribution around the Balmer discontinuity, which is highly sensitive to Teff and logg. We perform a calibration of the (lambda1, D) parameters into Teff using effective temperatures derived with the bolometric-flux method for 217 program stars, whose individual absolute relative uncertainties are on average 0.05. The final comparison of calculated with obtained Teff values in the (lambda1,D) calibration further show that the latter have total relative uncertainties, which on average are of the order of 0.05 for all spectral types and luminosity classes. The effective temperatures of OB supergiants derived in this work agree on average within some 2000 K with other determinations found in the literature, except those issued from wind-free non-LTE plane-parallel models of stellar atmospheres, which produce effective temperatures that can be overestimated by up to more than 5000 K near Teff = 25000 K.
• ### New high resolution synthetic stellar libraries for the Gaia Mission(0812.0297)

Dec. 1, 2008 astro-ph
High resolution synthetic stellar libraries are of fundamental importance for the preparation of the Gaia Mission. We present new sets of spectral stellar libraries covering two spectral ranges: 300 --1100 nm at 0.1 nm resolution, and 840 -- 890 nm at 0.001 nm resolution. These libraries span a large range in atmospheric parameters, from super-metal-rich to very metal-poor (-5.0 $<$[Fe/H]$<$+1.0), from cool to hot (\teff=3000--50000 K) stars, including peculiar abundance variations. The spectral resolution, spectral type coverage and number of models represent a substantial improvement over previous libraries used in population synthesis models and in atmospheric analysis.
• ### Fundamental Parameters of He-Weak and He-Strong Stars(0705.0541)

May 3, 2007 astro-ph
We carried out low resolution spectroscopic observations in the wavelength range 3400-4700 A of 20 He-weak and 8 He-strong stars to determine their fundamental parameters by means of the Divan-Chalonge-Barbier (BCD) spectrophotometric system. For a few He-weak stars we also estimate the effective temperatures and the angular diameters by integrating absolute fluxes observed over a wide spectral range. Non-LTE model calculations are carried out to study the influence of the He/H abundance ratio on the emergent radiation of He-strong stars and on their Teff determination. We find that the effective temperatures, surface gravities and bolometric absolute magnitudes of He-weak stars estimated with the BCD system and the integrated flux method are in good agreement between each other, and they also agree with previous determinations based on several different methods. The mean discrepancy between the visual absolute magnitudes derived using the Hipparcos parallaxes and the BCD values is on average 0.3 mag for He-weak stars, while it is 0.5 mag for He-strong stars. For He-strong stars, we note that the BCD calibration, based on stars in the solar environment, leads to overestimated values of Teff. By means of model atmosphere calculations with enhanced He/H abundance ratios we show that larger He/H ratios produce smaller BD which naturally explains the Teff overestimation. We take advantage of these calculations to introduce a method to estimate the He/H abundance ratio in He-strong stars. The BD of HD 37479 suggests that the Teff of this star remains fairly constant as the star spectrum undergoes changes in the intensity of H and He absorption lines. Data for the He-strong star HD 66765 are reported for the first time.
• ### Thermo-statistical description of gas mixtures from space partitions(cond-mat/0609598)

The new mathematical framework based on the free energy of pure classical fluids presented in [R. D. Rohrmann, Physica A 347, 221 (2005)] is extended to multi-component systems to determine thermodynamic and structural properties of chemically complex fluids. Presently, the theory focuses on $D$-dimensional mixtures in the low-density limit (packing factor $\eta < 0.01$). The formalism combines the free-energy minimization technique with space partitions that assign an available volume $v$ to each particle. $v$ is related to the closeness of the nearest neighbor and provides an useful tool to evaluate the perturbations experimented by particles in a fluid. The theory shows a close relationship between statistical geometry and statistical mechanics. New, unconventional thermodynamic variables and mathematical identities are derived as a result of the space division. Thermodynamic potentials $\mu_{il}$, conjugate variable of the populations $N_{il}$ of particles class $i$ with the nearest neighbors of class $l$ are defined and their relationships with the usual chemical potentials $\mu_i$ are established. Systems of hard spheres are treated as illustrative examples and their thermodynamics functions are derived analytically. The low-density expressions obtained agree nicely with those of scaled-particle theory and Percus-Yevick approximation. Several pair distribution functions are introduced and evaluated. Analytical expressions are also presented for hard spheres with attractive forces due to K\^ac-tails and square-well potentials. Finally, we derive general chemical equilibrium conditions.
• ### Rotation in the ZAMS: Be and Bn stars(astro-ph/0601067)

Jan. 4, 2006 astro-ph
We show that Be stars belong to a high velocity tail of a single B-type star rotational velocity distribution in the MS. This implies that: 1) the number fraction N(Be)/N(Be+B) is independent of the mass; 2) Bn stars having ZAMS rotational velocities higher than a given limit might become Be stars.
• ### Differential rotation in early type stars(astro-ph/0601068)

Jan. 4, 2006 astro-ph
Using 2D models of rotating stars, the interferometric measurements of alpha Eri and its fundamental parameters corrected for gravitational darkening effects we infer that the star might have a core rotating 2.7 times faster than the surface. We explore the consequences on spectral lines produced by surface differential rotation combined with the effects due to a kind of internal differential rotation with rotational energies higher than allowed for rigid rotation which induce geometrical deformations that do not distinguish strongly from those carried by the rigid rotation.
• ### On the evolutionary status of Be stars. I. Field Be stars near the Sun(astro-ph/0509119)

Sept. 6, 2005 astro-ph
A sample of 97 galactic field Be stars were studied by taking into account the effects induced by the fast rotation on their fundamental parameters. All program stars were observed in the BCD spectrophotometric system in order to minimize the perturbations produced by the circumstellar environment on the spectral photospheric signatures. This is one of the first attempts at determining stellar masses and ages by simultaneously using model atmospheres and evolutionary tracks, both calculated for rotating objects. The stellar ages ($\tau$) normalized to the respective inferred time that each rotating star can spend in the main sequence phase ($\tau\_{\rm MS}$) reveal a mass-dependent trend. This trend shows that: a) there are Be stars spread over the whole interval $0 \la \tau/\tau\_{\rm MS} \la 1$ of the main sequence evolutionary phase; b) the distribution of points in the ($\tau/\tau\_{\rm MS},M/M\_{\odot}$) diagram indicates that in massive stars ($M \ga 12M\_{\odot}$) the Be phenomenon is present at smaller $\tau/\tau\_{\rm MS}$ age ratios than for less massive stars ($M \la 12M\_{\odot}$). This distribution can be due to: $i$) higher mass-loss rates in massive objets, which can act to reduce the surface fast rotation; $ii$) circulation time scales to transport angular momentum from the core to the surface, which are longer the lower the stellar mass.
• ### Rotational velocities of A-type stars II. Measurement of vsini in the northern hemisphere(astro-ph/0205255)

May 16, 2002 astro-ph
This work is the second part of the set of measurements of vsini for A-type stars, begun by Royer et al. (2002). Spectra of 249 B8 to F2-type stars brighter than V=7 have been collected at Observatoire de Haute-Provence (OHP). Fourier transforms of several line profiles in the range 4200--4600 A are used to derive vsini from the frequency of the first zero. Statistical analysis of the sample indicates that measurement error mainly depends on vsini and this relative error of the rotational velocity is found to be about 5% on average. The systematic shift with respect to standard values from Slettebak et al. (1975), previously found in the first paper, is here confirmed. Comparisons with data from the literature agree with our findings: vsini values from Slettebak et al. are underestimated and the relation between both scales follows a linear law: vsini(new) = 1.03 vsini(old) + 7.7. Finally, these data are combined with those from the previous paper (Royer et al. 2002), together with the catalogue of Abt & Morrell (1995). The resulting sample includes some 2150 stars with homogenized rotational velocities.