The apparent lack of massive O-type stars near the zero-age main sequence (at ages < 2 Myr) is a topic widely discussed. Different explanations for this elusive detection have been proposed, but no firm conclusions have been reached yet. We reassess this empirical result benefiting from the high-quality spectroscopic observations of >400 Galactic O-type stars gathered by the IACOB and OWN surveys. We used temperatures and gravities from a iacob-gbat/fastwind spectroscopic analysis to locate our sample in the Kiel and spectroscopic HR diagrams. We evaluated the completeness of our sample of stars, observational biases using information from the Galactic O star catalog (GOSC), systematics of our methodology, and compare with other recent studies using smaller samples of Galactic O-type stars. We base our discussion on the spectroscopic HR diagram to avoid the use of uncertain distances. We performed a detailed study of the young cluster Trumpler-14 as an example of how Gaia cluster distances can help to construct the associated classical HR diagram. The apparent lack of massive O-type stars near the ZAMS with masses between 30 and 70 Msol persist even when spectroscopic results from a large, non-biased sample of stars are used. We do not find correlation between the dearth of stars and observational biases, limitations of our methodology, or the use of spectroscopic HR diagram instead of the classical one. Investigating the efficiency of mass accretion during the formation process we conclude that an adjustment of the accretion rate towards lower values could reconcile the hotter boundary of detected O-type stars and the theoretical birthline. Last, we discuss that the presence of a small sample of O2-O3.5 stars found closer to the ZAMS might be explained taking into account non-standard star evolution (e.g. binary interaction, mergers, or homogeneous evolution).

We address the origin of the observed bimodal rotational distribution of stars in massive young and intermediate age stellar clusters. This bimodality is seen as split main sequences at young ages and also has been recently directly observed in the $Vsini$ distribution of stars within massive young and intermediate age clusters. Previous models have invoked binary interactions as the origin of this bimodality, although these models are unable to reproduce all of the observational constraints on the problem. Here we suggest that such a bimodal rotational distribution is set up early within a cluster's life, i.e., within the first few Myr. Observations show that the period distribution of low-mass ($\la 2 M_\odot$) pre-main sequence (PMS) stars is bimodal in many young open clusters and we present a series of models to show that if such a bimodality exists for stars on the PMS that it is expected to manifest as a bimodal rotational velocity (at fixed mass/luminosity) on the main sequence for stars with masses in excess of $\sim1.5$~\msun. Such a bimodal period distribution of PMS stars may be caused by whether stars have lost (rapid rotators) or been able to retain (slow rotators) their circumstellar discs throughout their PMS lifetimes. We conclude with a series of predictions for observables based on our model.

Supermassive stars have been proposed as the progenitors of the massive ($\sim 10^{9}\,\rm{M}_{\odot}$) quasars observed at $z\sim7$. Prospects for directly detecting supermassive stars with next-generation facilities depend critically on their intrinsic lifetimes, as well as their formation rates. We use the 1D stellar evolution code Kepler to explore the theoretical limiting case of zero-metallicity, non-rotating stars, formed monolithically with initial masses between $10\,\rm{kM}_{\odot}$ and $190\,\rm{kM}_{\odot}$. We find that stars born with masses between $\sim60\,\rm{kM}_{\odot}$ and $\sim150\,\rm{kM}_{\odot}$ collapse at the end of the main sequence, burning stably for $\sim1.5\,\rm{Myr}$. More massive stars collapse directly through the general relativistic instability after only a thermal timescale of $\sim3\,\rm{kyr}$--$4\,\rm{kyr}$. The expected difficulty in producing such massive, thermally-relaxed objects, together with recent results for currently preferred rapidly-accreting formation models, suggests that such ``truly direct'' or ``dark'' collapses may not be typical for supermassive objects in the early Universe. We close by discussing the evolution of supermassive stars in the broader context of massive primordial stellar evolution and the possibility of supermassive stellar explosions.

We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.

Major mergers of gas-rich galaxies provide promising conditions for the formation of supermassive black holes (SMBHs; $\gtrsim10^5$ M$_\odot$) by direct collapse because they can trigger mass inflows as high as $10^4-10^5$ M$_\odot$ yr$^{-1}$ on sub-parsec scales. However, the channel of SMBH formation in this case, either dark collapse (direct collapse without prior stellar phase) or supermassive star (SMS; $\gtrsim10^4$ M$_\odot$), remains unknown. Here, we investigate the limit in accretion rate up to which stars can maintain hydrostatic equilibrium. We compute hydrostatic models of SMSs accreting at $1-1000$ M$_\odot$ yr$^{-1}$, and estimate the departures from equilibrium a posteriori by taking into account the finite speed of sound. We find that stars accreting above the atomic cooling limit ($\gtrsim10$ M$_\odot$ yr$^{-1}$) can only maintain hydrostatic equilibrium once they are supermassive. In this case, they evolve adiabatically with a hylotropic structure, that is, entropy is locally conserved and scales with the square root of the mass coordinate. Our results imply that stars can only become supermassive by accretion at the rates of atomically cooled haloes ($\sim0.1-10$ M$_\odot$ yr$^{-1}$). Once they are supermassive, larger rates are possible.

Globular clusters contain multiple stellar populations with peculiar chemical compositions. Pollution of the intracluster gas by an early population of fast-evolving stars is the most common scenario for explaining the observations. Stars with masses in excess of 1000 Msun have recently been suggested as potential polluters. We investigate the spectral properties of proto-GCs that would host a supermassive star. Our main goal is to quantify how such a star would affect the integrated light of the cluster, and to study the detectability of such objects. We computed non-LTE atmosphere models for SMS with various combinations of stellar parameters (luminosity, effective temperature, and mass) and metallicities appropriate for GCs, and we predict their emergent spectra. Using these spectra, we calculated the total emission of young proto-GCs with SMS as predicted in a previously reported scenario, and we computed synthetic photometry in UV, optical, and near-IR bands, in particular for the JWST. At an effective temperature of 10000 K, the spectrum of SMSs shows a Balmer break in emission. This feature is due to strong non-LTE effects (implied by the high luminosity) and is not observed in normal stars. The hydrogen lines also show a peculiar behavior, with Balmer lines in emission while higher series lines are in absorption. At 7000 K, the Balmer break shows a strong absorption. At high effective temperatures, the Lyman break is found in emission. Cool and luminous SMSs are found to dominate the integrated spectrum of the cluster, except for the UV range. The predicted magnitudes of these proto-GCs are m_AB~28-30 between 0.7 and 8 um and for redshifts z~4-10, which is detectable with the JWST. The peculiar observational features of cool SMSs imply that they might in principle be detected in color-color diagrams that probe the spectral energy distribution below and above the Balmer break.

Grids of stellar models are useful tools to derive the properties of stellar clusters, in particular young clusters hosting massive stars, and to provide information on the star formation process in various mass ranges. Because of their short evolutionary timescale, massive stars end their life while their low-mass siblings are still on the pre-main sequence (pre-MS) phase. Thus the study of young clusters requires consistent consideration of all the phases of stellar evolution. But despite the large number of grids that are available in the literature, a grid accounting for the evolution from the pre-MS accretion phase to the post-MS phase in the whole stellar mass range is still lacking. We build a grid of stellar models at solar metallicity with masses from 0.8 $M_\odot$ to 120 $M_\odot$, including pre-MS phase with accretion. We use the {\sc genec} code to run stellar models on this mass range. The accretion law is chosen to match the observations of pre-MS objects on the Hertzsprung-Russell diagram. We describe the evolutionary tracks and isochrones of our models. The grid is connected to previous MS and post-MS grids computed with the same numerical method and physical assumptions, which provides the widest grid in mass and age to date. Numerical tables of our models and corresponding isochrones are available online.

Supermassive stars (SMSs) are candidates for being progenitors of supermassive quasars at high redshifts. However, their formation process requires strong mechanisms that would be able to extract the angular momentum of the gas that the SMSs accrete. We investigate under which conditions the magnetic coupling between an accreting SMS and its winds can remove enough angular momentum for accretion to proceed from a Keplerian disc. We numerically computed the rotational properties of accreting SMSs that rotate at the $\Omega\Gamma$-limit and estimated the magnetic field that is required to maintain the rotation velocity at this limit using prescriptions from magnetohydrodynamical simulations of stellar winds. We find that a magnetic field of 10 kG at the stellar surface is required to satisfy the constraints on stellar rotation from the $\Omega\Gamma$-limit. Magnetic coupling between the envelope of SMSs and their winds could allow for SMS formation by accretion from a Keplerian disc, provided the magnetic field is at the upper end of present-day observed stellar fields. Such fields are consistent with primordial origins.

Massive protostars grow and evolve under the effect of rapid accretion of circumstellar gas and dust, falling at high rates ($\ge 10^{-4}$-$10^{-3}\, \rm M_{\odot}\, \rm yr^{-1}$). This mass infall has been shown, both numerically and observationally, to be episodically interspersed by accretion of dense gaseous clumps migrating through the circumstellar disc to the protostellar surface, causing sudden accretion and luminous bursts. Using numerical gravito-radiation-hydrodynamics and stellar evolution calculations, we demonstrate that, in addition to the known bloating of massive protostars, variable episodic accretion further influences their evolutionary tracks of massive young stellar objects (MYSOs). For each accretion-driven flare, they experience rapid excursions toward more luminous, but colder regions of the Hertzsprung-Russell diagram. During these excursions, which can occur up to the end of the pre-main-sequence evolution, the photosphere of massive protostars can episodically release much less energetic photons and MYSOs surreptitiously adopt the same spectral type as evolved massive (supergiants) stars. Each of these evolutionary loop brings the young high-mass stars close to the forbidden Hayashi region and might make their surrounding H II regions occasionally fainter, before they recover their quiescent, pre-burst surface properties. We interpret such cold, intermittent pre-main-sequence stellar evolutionary excursions and the dipping variability of HII regions as the signature of the presence of a fragmenting circumstellar accretion disc surrounding the MYSOs. We conjecture that this mechanism might equivalently affect young stars in the intermediate-mass regime.

The collapse of supermassive primordial stars in hot, atomically-cooled halos may have given birth to the first quasars at $z \sim$ 15 - 20. Recent numerical simulations of these rapidly accreting stars reveal that they are cool, red hypergiants shrouded by dense envelopes of pristine atomically-cooled gas at 6,000 - 8,000 K, with luminosities $L$ $\gtrsim$ 10$^{10}$ L$_{\odot}$. Could such luminous but cool objects be detected as the first stage of quasar formation in future near infrared (NIR) surveys? We have now calculated the spectra of supermassive primordial stars in their birth envelopes with the Cloudy code. We find that some of these stars will be visible to JWST at $z \lesssim$ 20 and that with modest gravitational lensing Euclid and WFIRST could detect them out to $z \sim$ 10 - 12. Rather than obscuring the star, its accretion envelope enhances its visibility in the NIR today by reprocessing its short-wavelength flux into photons that are just redward of the Lyman limit in the rest frame of the star.

We explore the initial conditions that lead to detached binary evolution of massive pre-MS binaries and ask how large a fraction of the observed binary systems may have been initially formed as low-mass protobinaries and later undergone a significant accretion phase while remaining detached. We develop a family of analytic models to describe the orbital separation, $a$, and mass ratio, $q$, evolution. For a given mass accretion rate onto the binary system, we define a recipe for distributing this mass between the two components. For this we introduce a parameter $\eta$, such that $\frac{\dot{M}_2}{\dot{M}_1}=q^\eta$ at any time, to determine the binary mass ratio evolution. Depending on the choice of $\eta$, any type of mass ratio evolution is possible. Furthermore, we use MESA, a detailed stellar structure code, to calculate an extensive grid of binary sequences where a protobinary undergoes accretion, and we identify the initial conditions that separates detached from non-detached pre-MS binary evolution. A value of $\eta$ around 2 allows accretion growth in detached systems to form close massive binaries on the Zero Age Main-Sequence with minimum orbital periods down to about 1.2 days for $M_{\rm 1, ZAMS}=20-30M_{\odot}$ twin-binaries. $\eta=2$ can also reproduce the observed population of binary systems with primary stars above $6M_{\odot}$. The whole observed range of massive close binaries can form via accretion growth in detached systems, making the binary formation channel of accretion growth a strong contender to explain the formation of massive close binaries, including progenitors of coalescing binary black holes.

When planets are formed from the protoplanetary disk and after the disk has dissipated, the evolution of their orbits is governed by tidal interactions, friction, and gravitational drag, and also by changes in the mass of the star and planet. These interactions may change the initial distribution of the distances between the planets and their host star by expanding the original orbit, by contracting it (which may cause an engulfment of the planet by the star), or by destroying the planet. We study the evolution of the orbit of a planet orbiting its host star under the effects of equilibrium tides, dynamical tides, drag (frictional and gravitational), and stellar mass loss. We used the Geneva stellar evolution code to compute the evolution of stars with initial masses of 1 and 1.5 solar mass with different rotation rates at solar metallicity. The star is evolved from the pre-main-sequence (PMS) up to the tip of the red giant branch. We used these models as input for computing the evolution of the planetary orbits. We explored the effects of changing the planet masses (of 1 Earth mass up to 20 Jupiter masses), the distance between the planet and the star (of 0.015 and more than 3 au), the mass, and the spin of the star. We present results when only the equilibrium tide was accounted for and when both equilibrium and dynamical tides were accounted for. The expression for the dynamical tide is a frequency-averaged dissipation of tidally excited inertial waves, obtained from a piecewise homogeneous two-layer stellar model. Gravity wave damping was neglected. Dynamical tides in convective zones have a significant effect on planetary orbits only during the PMS phase and only for fast-rotating stars. They have no significant effects during the PMS phase for initially slow-rotating stars and during the red giant branch phase, regardless of the initial rotation (abridged).

Supermassive stars born from pristine gas in atomically-cooled haloes are thought to be the progenitors of supermassive black holes at high redshifts. However, the way they accrete their mass is still an unsolved problem. In particular, for accretion to proceed, a large amount of angular momentum has to be extracted from the collapsing gas. Here, we investigate the constraints stellar evolution imposes on this angular momentum problem. We present an evolution model of a supermassive Population III star including simultaneously accretion and rotation. We find that, for supermassive stars to form by accretion, the accreted angular momentum has to be about 1% of the Keplerian angular momentum. This tight constraint comes from the $\Omega\Gamma$-limit, at which the combination of radiation pressure and centrifugal force cancels gravity. It implies that supermassive stars are slow rotators, with a surface velocity less than 10-20% of their first critical velocity, at which the centrifugal force alone cancels gravity. At such low velocities, the deformation of the star due to rotation is negligible.

Supermassive primordial stars forming in atomically-cooled halos at $z \sim15-20$ are currently thought to be the progenitors of the earliest quasars in the Universe. In this picture, the star evolves under accretion rates of $0.1 - 1$ $M_\odot$ yr$^{-1}$ until the general relativistic instability triggers its collapse to a black hole at masses of $\sim10^5$ $M_\odot$. However, the ability of the accretion flow to sustain such high rates depends crucially on the photospheric properties of the accreting star, because its ionising radiation could reduce or even halt accretion. Here we present new models of supermassive Population III protostars accreting at rates $0.001 - 10$ $M_\odot$ yr$^{-1}$, computed with the GENEVA stellar evolution code including general relativistic corrections to the internal structure. We use the polytropic stability criterion to estimate the mass at which the collapse occurs, which has been shown to give a lower limit of the actual mass at collapse in recent hydrodynamic simulations. We find that at accretion rates higher than $0.001$ $M_\odot$ yr$^{-1}$ the stars evolve as red, cool supergiants with surface temperatures below $10^4$ K towards masses $>10^5$ $M_\odot$, and become blue and hot, with surface temperatures above $10^5$ K, only for rates $\lesssim0.001$ $M_\odot$ yr$^{-1}$. Compared to previous studies, our results extend the range of masses and accretion rates at which the ionising feedback remains weak, reinforcing the case for direct collapse as the origin of the first quasars.

Rotation plays a key role in the star-formation process, from pre-stellar cores to pre-main-sequence (PMS) objects. Understanding the formation of massive stars requires taking into account the accretion of angular momentum during their PMS phase. We study the PMS evolution of objects destined to become massive stars by accretion, focusing on the links between the physical conditions of the environment and the rotational properties of young stars. In particular, we look at the physical conditions that allow the production of massive stars by accretion. We present PMS models computed with a new version of the Geneva Stellar Evolution code self-consistently including accretion and rotation according to various accretion scenarios for mass and angular momentum. We describe the internal distribution of angular momentum in PMS stars accreting at high rates and we show how the various physical conditions impact their internal structures, evolutionary tracks, and rotation velocities during the PMS and the early main sequence. We find that the smooth angular momentum accretion considered in previous studies leads to an angular momentum barrier and does not allow the formation of massive stars by accretion. A braking mechanism is needed in order to circumvent this angular momentum barrier. This mechanism has to be efficient enough to remove more than 2/3 of the angular momentum from the inner accretion disc. Due to the weak efficiency of angular momentum transport by shear instability and meridional circulation during the accretion phase, the internal rotation profiles of accreting stars reflect essentially the angular momentum accretion history. As a consequence, careful choice of the angular momentum accretion history allows circumvention of any limitation in mass and velocity, and production of stars of any mass and velocity compatible with structure equations.

Massive star formation requires the accretion of gas at high rate while the star is already bright. Its actual luminosity depends sensitively on the stellar structure. We compute pre-main-sequence tracks for massive and intermediate-mass stars with variable accretion rates and study the evolution of stellar radius, effective temperature and ionizing luminosity, starting at $2\,M_\odot$ with convective or radiative structures. The radiative case shows a much stronger swelling of the protostar for high accretion rates than the convective case. For radiative structures, the star is very sensitive to the accretion rate and reacts quickly to accretion bursts, leading to considerable changes in photospheric properties on timescales as short as 100 - 1000 yr. The evolution for convective structures is much less influenced by the instantaneous accretion rate, and produces a monotonically increasing ionizing flux that can be many orders of magnitude smaller than in the radiative case. For massive stars, it results in a delay of the HII region expansion by up to 10,000 yr. In the radiative case, the HII region can potentially be engulfed by the star during the swelling, which never happens in the convective case. We conclude that the early stellar structure has a large impact on the radiative feedback during the pre-main-sequence evolution of massive protostars and introduces an important uncertainty that should be taken into account. Because of their lower effective temperatures, our convective models may hint at a solution to an observed discrepancy between the luminosity distribution functions of massive young stellar objects and compact HII regions.

Massive stars likely form by accretion and the evolutionary track of an accreting forming star corresponds to what is called the birthline in the HR diagram. The shape of this birthline is quite sensitive to the evolution of the entropy in the accreting star. We first study the reasons why some birthlines published in past years present different behaviours for a given accretion rate. We then revisit the question of the accretion rate, which allows us to understand the distribution of the observed pre-main-sequence (pre-MS) stars in the Hertzsprung-Russell (HR) diagram. Finally, we identify the conditions needed to obtain a large inflation of the star along its pre-MS evolution that may push the birthline towards the Hayashi line in the upper part of the HR diagram. We present new pre-MS models including accretion at various rates and for different initial structures of the accreting core. From the observed upper envelope of pre-MS stars in the HR diagram, we deduce the accretion law that best matches the accretion history of most of the intermediate-mass stars. In the case of cold disc accretion, the existence of a significant swelling during the accretion phase, which leads to radii $\gtrsim100\,R_\odot$ and brings the star back to the red part of the HR diagram, depends sensitively on the initial conditions. For an accretion rate of $10^{-3}\,M_\odot\rm\,yr^{-1}$, only models starting from a core with a significant radiative region evolve back to the red part of the HR diagram. We also obtain that, in order to reproduce the observed upper envelope of pre-MS stars in the HR diagram with an accretion law deduced from the observed mass outflows in ultra-compact HII regions, the mass effectively accreted onto the star with respect to the total in falling matter decreases when the mass of the star increases.

The database of intermediate mass rotating stellar models presented in the past years by the Geneva Stellar Evolution Group can be used to build synthetic stellar populations that fully account for the effects of stellar rotation. However, up to now we still lacked stellar evolutionary tracks that rotate close to the critical limit during the whole MS phase. This occurs because the flat internal profile of rotation imposed at the Zero-Age MS (ZAMS) is modified by the action of meridional currents immediately after the ZAMS, causing the surface rotational velocity to decrease abruptly until it reaches a quasi-stationary state. We compute stellar models with non-solid rotation at the ZAMS to obtain stellar evolutionary tracks with a larger content of angular momentum, that attain rotational equatorial velocities close to the critical limit throughout their MS phase. These models have a longer MS lifetime and a higher surface chemical enrichment already at the end of the MS, particularly at Z=0.002. Stellar models with solid rotation at the ZAMS adequately represent the overall characteristics and evolution of differentially rotating models of identical angular momentum content, but with a lower initial surface rotational velocity. For these models we recommend to use as the initial rotational rate the values derived once the quasi-stationary state is reached, after the abrupt decrease in surface velocity. Interestingly, the initial equatorial rotational velocities are virtually metallicity independent for the stellar models we computed with the same mass and angular momentum content at the ZAMS. If, as some observational evidence indicates, B-type stars at Z=0.002 rotate with a higher equatorial velocity at the ZAMS than stars with Z=0.014, our finding would indicate that the angular momentum content of B-type stars in the SMC is higher than their Galactic counterparts.

(shortened) We provide a grid of single star models covering a mass range from 0.8 to 120 Msun with an initial metallicity Z = 0.002 with and without rotation. We discuss the impact of a change in the metallicity by comparing the current tracks with models computed with exactly the same physical ingredients but with a metallicity Z = 0.014 (solar). We show that the width of the main-sequence (MS) band in the upper part of the Hertzsprung-Russell diagram (HRD), for luminosity above log(L/Lsun) > 5.5, is very sensitive to rotational mixing. Strong mixing significantly reduces the MS width. We confirm, but here for the first time on the whole mass range, that surface enrichments are stronger at low metallicity provided that comparisons are made for equivalent initial mass, rotation and evolutionary stage. We show that the enhancement factor due to a lowering of the metallicity (all other factors kept constant) increases when the initial mass decreases. Present models predict an upper luminosity for the red supergiants (RSG) of log (L/Lsun) around 5.5 at Z = 0.002 in agreement with the observed upper limit of RSG in the Small Magellanic Cloud. We show that models using shear diffusion coefficient calibrated to reproduce the surface enrichments observed for MS B-type stars at Z = 0.014 can also reproduce the stronger enrichments observed at low metallicity. In the framework of the present models, we discuss the factors governing the timescale of the first crossing of the Hertzsprung gap after the MS phase. We show that any process favouring a deep localisation of the H-burning shell (steep gradient at the border of the H-burning convective core, low CNO content) and/or the low opacity of the H-rich envelope favour a blue position in the HRD for the whole or at least a significant fraction of the core He-burning phase.

The way angular momentum is built up in stars during their formation process may have an impact on their further evolution. In the frame of the cold disc accretion scenario, we study for the first time how angular momentum builds up inside the star during its formation and what are the consequences for its evolution on the main sequence (MS). Computation begins from a hydrostatic core on the Hayashi line of 0.7 Msol at solar metallicity (Z=0.014) rotating as a solid body. Accretion rates depending on the luminosity of the accreting object are considered varying between 1.5e-5 and 1.7e-3 Msol/yr. The accreted matter is assumed to have an angular velocity equal to that of the outer layer of the accreting star. Models are computed for a mass-range on the zero-age main sequence (ZAMS) between 2 and 22 Msol. We study how the internal and surface velocities vary as a function of time during the accretion phase and the evolution towards the ZAMS. Stellar models, whose evolution has been followed along the pre-MS phase, are found to exhibit a shallow gradient of angular velocity on the ZAMS. Interestingly, for masses on the ZAMS larger than 8 Msol, there exists a maximum surface velocity that can be reached through the present scenario of formation. Typically, for 14 Msol models, only stars with surface velocities on the ZAMS lower than about 45% of the critical velocity can be formed. To reach higher velocities would require to start from cores rotating above the critical limit. We find that this upper velocity limit is smaller for higher masses. In contrast, below 8 Msol, there is no restriction and the whole domain of velocities, up to the critical one, can be reached.

The rotation of stars has many interesting and important consequences for the photometric and chemical evolution of galaxies. Many of the predictions of models of stellar rotation are now compared with observations of surface abundances and velocities, with interferometric studies of fast rotating stars, with internal rotation profiles as they can be deduced by asteroseismology, to cite just a few observational constraints. In this paper, we investigate how the outputs of models depend on the prescriptions used for the diffusion coefficients included in the shellular rotating models. After recalling the various prescriptions found in the literature, we discuss their impacts on the evolutionary tracks and lifetimes of the Main-Sequence (MS) phase, the changes of the surface composition and velocities during the MS phase, the distribution of the core helium lifetime in the blue and the red part of the HR diagram, the extensions of the blue loops, the evolution of the angular momentum of the core, and the synthesis of primary nitrogen in fast-rotating metal-poor massive stars. While some of these outputs depend only slightly on the prescriptions used (for instance, the evolution of the surface velocities), most of them show a significant dependence. The models which best fit the changes of the surface abundances are those computed with the vertical shear diffusion coefficient of Maeder (1997) and the horizontal shear diffusion coefficient by Zahn (1992).

Aims: We study the influence of rotation and disc lifetime on lithium depletion of pre-main sequence (PMS) solar-type stars. Methods: The impact of rotational mixing and of the hydrostatic effects of rotation on lithium abundances are investigated by computing non-rotating and rotating PMS models that include a comprehensive treatment of shellular rotation. The influence of the disc lifetime is then studied by comparing the lithium content of PMS rotating models experiencing different durations of the disc-locking phase between 3 and 9 Myr. Results: The surface lithium abundance at the end of the PMS is decreased when rotational effects are included. During the beginning of the lithium depletion phase, only hydrostatic effects of rotation are at work. This results in a decrease in the lithium depletion rate for rotating models compared to non-rotating ones. When the convective envelope recedes from the stellar centre, rotational mixing begins to play an important role due to differential rotation near the bottom of the convective envelope. This mixing results in a decrease in the surface lithium abundance with a limited contribution from hydrostatic effects of rotation, which favours lithium depletion during the second part of the PMS evolution. The impact of rotation on PMS lithium depletion is also found to be sensitive to the duration of the disc-locking phase. When the disc lifetime increases, the PMS lithium abundance of a solar-type star decreases owing to the higher efficiency of rotational mixing in the radiative zone. A relationship between the surface rotation and lithium abundance at the end of the PMS is then obtained: slow rotators on the zero-age main sequence are predicted to be more lithium-depleted than fast rotators due to the increase in the disc lifetime.