The understanding of the rotational evolution of early-type stars is deeply related to that of anisotropic mass and angular momentum loss. In this paper, we aim to clarify the rotational evolution of rapidly rotating early-type stars along the main sequence (MS). We have used the 2D ESTER code to compute and evolve isolated rapidly rotating early-type stellar models along the MS, with and without anisotropic mass loss. We show that stars with $Z=0.02$ and masses between $5$ and $7~M_\odot$ reach criticality during the main sequence provided their initial angular velocity is larger than 50% of the Keplerian one. More massive stars are subject to radiation-driven winds and to an associated loss of mass and angular momentum. We find that this angular momentum extraction from the outer layers can prevent massive stars from reaching critical rotation and greatly reduce the degree of criticality at the end of the MS. Our model includes the so-called bi-stability jump of the $\dot{M}-T_{\rm eff}$ relation of 1D-models. This discontinuity now shows up in the latitude variations of the mass-flux surface density, endowing rotating massive stars with either a single-wind regime (no discontinuity) or a two-wind regime (a discontinuity). In the two-winds-regime, mass loss and angular momentum loss are strongly increased at low latitudes inducing a faster slow-down of the rotation. However, predicting the rotational fate of a massive star is difficult, mainly because of the non-linearity of the phenomena involved and their strong dependence on uncertain prescriptions. Moreover, the very existence of the bi-stability jump in mass-loss rate remains to be substantiated by observations.

The understanding of the evolution of early-type stars is tightly related to that of the effects of rapid rotation. For massive stars, rapid rotation combines with their strong radiation-driven wind. The aim of this paper is to investigate two questions that are prerequisite to the study of the evolution of massive rapidly rotating stars: (i) What is the critical angular velocity of a star when radiative acceleration is significant in its atmosphere? (ii) How do mass and angular momentum loss depend on the rotation rate? To investigate fast rotation, which makes stars oblate, we used the 2D ESTER models and a simplified approach, the $\omega$-model, which gives the latitudinal dependence of the radiative flux in a centrifugally flattened radiative envelope. We find that radiative acceleration only mildly influences the critical angular velocity, at least for stars with masses lower than 40 Msun. We explain this mild reduction of the critical angular velocity compared to the classical Keplerian angular velocity by the combined effects of gravity darkening and a reduced equatorial opacity that is due to the centrifugal acceleration. To answer the second question, we first devised a model of the local surface mass flux, which we calibrated with previously developed 1D models. The discontinuity (the so-called bi-stability jump) included in the $\dot{M}-T_{\rm eff}$ relation of 1D models means that the mass flux of a fast-rotating star is controlled by either a single wind or a two-wind regime. Mass and angular momentum losses are strong around the equator if the star is in the two-wind regime. We also show that the difficulty of selecting massive stars that are viewed pole-on makes detecting the discontinuity in the relation between mass loss and effective temperature also quite challenging.

Gravity darkening (GD) and flattening are important consequences of stellar rotation. The precise characterization of these effects across the HRD is crucial to a deeper understanding of stellar structure and evolution. We seek to characterize such important effects on Sargas, an evolved, fast-rotating, intermediate-mass star, located in a region of the HRD where they have never been directly measured as far as we know. We use our numerical model CHARRON to analyze interferometric (VLTI/PIONIER) and spectroscopic (VLT/UVES) observations through a MCMC model-fitting procedure. The visibilities and closure phases from the PIONIER data are particularly sensitive to rotational flattening and GD. Adopting the Roche approximation, we investigate two GD models: (1) the beta-model (classical von Zeipel's law), and (2) the omega-model. Using this approach we measure several physical parameters of Sargas, namely, equatorial radius, mass, equatorial rotation velocity, mean Teff, inclination and position angle of the rotation axis, and beta. In particular, we show that the measured beta leads to a surface flux distribution equivalent to the one given by the omega-model. Thanks to our results, we also show that Sargas is most probably located in a rare and interesting region of the H-R diagram: within the Hertzsprung gap and over the hot edge of the instability strip. These results show once more the power of optical/IR long-baseline interferometry, combined with high-resolution spectroscopy, to directly measure fast-rotation effects and stellar parameters, in particular GD. As was the case for a few fast rotators previously studied by interferometry, the omega-model provides a physically more profound description of Sargas' GD, without the need of a beta exponent.

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.

Solar Orbiter is a joint ESA-NASA mission planed for launch in October 2018. The science payload includes remote-sensing and in-situ instrumentation designed with the primary goal of understanding how the Sun creates and controls the heliosphere. The spacecraft will follow an elliptical orbit around the Sun, with perihelion as close as 0.28 AU. During the late orbit phase the orbital plane will reach inclinations above 30 degrees, allowing direct observations of the solar polar regions. The Energetic Particle Detector (EPD) is an instrument suite consisting of several sensors measuring electrons, protons and ions over a broad energy interval (2 keV to 15 MeV for electrons, 3 keV to 100 MeV for protons and few tens of keV/nuc to 450 MeV/nuc for ions), providing composition, spectra, timing and anisotropy information. We present an overview of Solar Orbiter from the energetic particle perspective, summarizing the capabilities of EPD and the opportunities that these new observations will provide for understanding how energetic particles are accelerated during solar eruptions and how they propagate through the Heliosphere.

Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditionner of the Jacobian matrix.

Recent observations of rapidly rotating stars have revealed the presence of regular patterns in their pulsation spectra. This has raised the question as to their physical origin, and in particular, whether they can be explained by an asymptotic frequency formula for low-degree acoustic modes, as recently discovered through numerical calculations and theoretical considerations. In this context, a key question is whether compositional/density gradients can adversely affect such patterns to the point of hindering their identification. To answer this question, we calculate frequency spectra using two-dimensional ESTER stellar models. These models use a multi-domain spectral approach, allowing us to easily insert a compositional discontinuity while retaining a high numerical accuracy. We analyse the effects of such discontinuities on both the frequencies and eigenfunctions of pulsation modes in the asymptotic regime. We find that although there is more scatter around the asymptotic frequency formula, the semi-large frequency separation can still be clearly identified in a spectrum of low-degree acoustic modes.

Two-dimensional models of rapidly rotating stars are already unavoidable for the interpretation of interferometric or asteroseismic data of this kind of stars. When combined with time evolution, they will allow the including of a more accurate physics for the computation of element transport and the determination of surface abundances. In addition, modeling the evolution of rotation will improve gyrochronology. Presently, two-dimensional ESTER models predict the structure and the large-scale flows (differential rotation and meridional circulation) of stars with mass larger than 1.7Msun at any rotation rate. Main sequence evolution can be mimicked by varying the hydrogen content of the convective core. Models have been successfully tested on half a dozen of nearby fast rotating stars observed with optical or infra-red interferometers. They are now the right tool to investigate the oscillation spectrum of early-type fast rotators.

ESTER two-dimensional models solve the steady state structure of fast rotating early-type stars including the large scale flows associated with the baroclinicity of the radiative zones. Models are compared successfully to the fundamental parameters of the two main components of the triple system $\delta$ Velorum that have been derived from interferometric and orbit measurements. Testing the models on the Be star Achernar ($\alpha$ Eri), we cannot reproduce the data and conclude that this star has left the main sequence and is likely crossing the Herzsprung gap. Computing main sequence evolution of fast rotating stars at constant angular momentum shows that their criticality increases with time suggesting that the Be phenomenon and the ensuing mass ejections is the result of evolution.

Many early-type stars have been measured with high angular velocities. In such stars, mode identification is difficult as the effects of fast and differential rotation are not well known. Using fundamental parameters measured by interferometry, the ESTER structure code and the TOP oscillation code, we investigate the oscillation spectrum of Rasalhague (alpha Ophiuchi), for which observations by the MOST satellite found 57 oscillations frequencies. Results do not show a clear identification of the modes and highlight the difficulties of asteroseismology for such stars with a very complex oscillation spectrum.

This work aims at presenting the first two-dimensional models of an isolated rapidly rotating star that include the derivation of the differential rotation and meridional circulation in a self-consistent way.We use spectral methods in multidomains, together with a Newton algorithm to determine the steady state solutions including differential rotation and meridional circulation for an isolated non-magnetic, rapidly rotating early-type star. In particular we devise an asymptotic method for small Ekman numbers (small viscosities) that removes the Ekman boundary layer and lifts the degeneracy of the inviscid baroclinic solutions.For the first time, realistic two-dimensional models of fast-rotating stars are computed with the actual baroclinic flows that predict the differential rotation and the meridional circulation for intermediate-mass and massive stars. These models nicely compare with available data of some nearby fast-rotating early-type stars like Ras Alhague ($\alpha$ Oph), Regulus ($\alpha$ Leo), and Vega ($\alpha$ Lyr). It is shown that baroclinicity drives a differential rotation with a slow pole, a fast equator, a fast core, and a slow envelope. The differential rotation is found to increase with mass, with evolution (here measured by the hydrogen mass fraction in the core), and with metallicity. The core-envelope interface is found to be a place of strong shear where mixing will be efficient.Two-dimensional models offer a new view of fast-rotating stars, especially of their differential rotation, which turns out to be strong at the core-envelope interface. They also offer more accurate models for interpreting the interferometric and spectroscopic data of early-type stars.

Context.Interpretation of light curves of many types of binary stars requires the inclusion of the (cor)relation between surface brightness and local effective gravity. Until recently, this correlation has always been modeled by a power law relating the flux or the effective temperature and the effective gravity, namely T_eff {\alpha} g_eff^{\beta}. Aims. We look for a simple model that can describe the variations of the flux at the surface of stars belonging to a binary system. Methods. This model assumes that the energy flux is a divergence-free vector anti-parallel to the effective gravity. The effective gravity is computed from the Roche model. Results. After explaining in a simple manner the old result of Lucy (1967), which says that {\beta}=0.08 for solar type stars, we first argue that one-dimensional models should no longer be used to evaluate gravity darkening laws. We compute the correlation between log T_eff and log g_eff using a new approach that is valid for synchronous, weakly magnetized, weakly irradiated binaries. We show that this correlation is approximately linear, validating the use of a power law relation between effective temperature and effective gravity as a first approximation. We further show that the exponent {\beta} of this power law is a slowly varying function, which we tabulate, of the mass ratio of the binary star and the Roche lobe filling factor of the stars of the system. The exponent {\beta} remains mostly in the interval (0.20, 0.25) if extreme mass ratios are eliminated. Conclusions. For binary stars that are synchronous, weakly magnetized and weakly irradiated, the gravity darkening exponent is well constrained and may be removed from the free parameters of the models.

Modelling isolated rotating stars at any rotation rate is a challenge for the next generation of stellar models. These models will couple dynamical aspects of rotating stars, like angular momentum and chemicals transport, with classical chemical evolution, gravitational contraction or mass-loss. Such modelling needs to be achieved in two dimensions, combining the calculation of the structure of the star, its mean flows and the time-evolution of the whole. We present here a first step in this challenging programme. It leads to the first self-consistent two-dimensional models of rotating stars in a steady state generated by the ESTER code. In these models the structure (pressure, density and temperature) and the flow fields are computed in a self-consistent way allowing the prediction of the differential rotation and the associated meridian circulation of the stars. After a presentation of the physical properties of such models and the numerical methods at work, we give the first grid of such models describing massive and intermediate-mass stars for a selection of rotation rates up to 90% of the breakup angular velocity.

Interpretation of interferometric observations of rapidly rotating stars requires a good model of their surface effective temperature. Until now, laws of the form T_eff \propto g_eff^{\beta} have been used, but they are only valid for slowly rotating stars. We propose a simple model that can describe the latitudinal variations in the flux of rotating stars at any rotation rate. This model assumes that the energy flux is a divergence-free vector that is antiparallel to the effective gravity. When mass distribution can be described by a Roche model, the latitudinal variations in the effective temperature only depend on a single parameter, namely the ratio of the equatorial velocity to the Keplerian velocity. We validate this model by comparing its predictions to those of the most realistic two-dimensional models of rotating stars issued from the ESTER code. The agreement is very good, as it is with the observations of two rapidly rotating stars, {\alpha} Aql and {\alpha} Leo. We suggest that as long as a gray atmosphere can be accepted, the inversion of data on flux distribution coming from interferometric observations of rotating stars uses such a model, which has just one free parameter.

Classical Be stars are rapid rotators of spectral type late O to early A and luminosity class V-III, wich exhibit Balmer emission lines and often a near infrared excess originating in an equatorially concentrated circumstellar envelope, both produced by sporadic mass ejection episodes. The causes of the abnormal mass loss (the so-called Be phenomenon) are as yet unknown. For the first time, we can now study in detail Be stars outside the Earth's atmosphere with sufficient temporal resolution. We investigate the variability of the Be Star CoRoT-ID 102761769 observed with the CoRoT satellite in the exoplanet field during the initial run. One low-resolution spectrum of the star was obtained with the INT telescope at the Observatorio del Roque de los Muchachos. A time series analysis was performed using both cleanest and singular spectrum analysis algorithms to the CoRoT light curve. To identify the pulsation modes of the observed frequencies, we computed a set of models representative of CoRoT-ID 102761769 by varying its main physical parameters inside the uncertainties discussed. We found two close frequencies related to the star. They are 2.465 $\rm c\,d^{-1}$ (28.5 $\mathrm{\mu Hz}$) and 2.441 $\rm c\,d^{-1}$ (28.2 $\mathrm{\mu Hz}$). The precision to which those frequencies were found is 0.018 $\rm c\,d^{-1}$ (0.2 $\mathrm{\mu Hz}$). The projected stellar rotation was estimated to be 120 $\rm km\,s^{-1}$ from the Fourier transform of spectral lines. If CoRoT-ID 102761769 is a typical Galactic Be star it rotates near the critical velocity. The critical rotation frequency of a typical B5-6 star is about 3.5 $\rm c\,d^{-1}$(40.5 $\mathrm{\mu Hz}$), which implies that the above frequencies are really caused by stellar pulsations rather than star's rotation.

Recent results from interferometry and asteroseismology require models of rapidly rotating stars that are more and more precise. We describe the basic structure and the hydrodynamics of a fully radiative star as a preliminary step towards more realistic models of rotating stars. We consider a solar mass of perfect gas enclosed in a spherical container. The gas is self-gravitating and rotating, and is the seat of nuclear heating, and heat diffusion is due to radiative diffusion with Kramers type opacities. Equations are solved numerically with spectral methods in two dimensions with a radial Gauss-Lobatto grid and spherical harmonics. We computed the centrifugally flattened structure of such a star: the von Zeipel model, which says that the energy flux is proportional to the local effective gravity is tested. We show that it overestimates the ratio of the polar to the equatorial energy flux by almost a factor 2. We also determine the Brunt-Vaisala frequency distribution and show that outer equatorial regions in a radiative zone are convectively unstable when the rotation is fast enough. We compute the differential rotation and meridional circulation stemming from the baroclinicity of the star and show that, in such radiative zones, equatorial regions rotate faster than polar ones. The surface differential rotation is also shown to reach a universal profile when rotation is slow enough (less than 36% of the breakup one), as long as viscosity and Prandlt number remain small.