Aims. We study the impact of a subsolar metallicity on various properties of non-rotating and rotating stars, such as surface velocities and abundances, lifetimes, evolutionary tracks, and evolutionary scenarios. Methods. We provide a grid of single star models covering a mass range of 0.8 to 120 M 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). Results. 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/L ) > 5.5, is very sensitive to rotational mixing. Strong mixing significantly reduces the MS width. Here for the first time over the whole mass range, we confirm 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/L ) 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, which is 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.

Supermassive primordial stars forming in atomically-cooled halos at z ∼ 15 − 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 yr −1 until the general relativistic instability triggers its collapse to a black hole at masses of ∼ 10 5 M . 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 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 yr −1 the stars evolve as red, cool supergiants with surface temperatures below 10 4 K towards masses > 10 5 M , and become blue and hot, with surface temperatures above 10 5 K, only for rates 0.001 M 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.

Supermassive primordial stars are suspected to be the progenitors of the most massive quasars at z ∼ 6. Previous studies of such stars were either unable to resolve hydrodynamical timescales or considered stars in isolation, not in the extreme accretion flows in which they actually form. Therefore, they could not self-consistently predict their final masses at collapse, or those of the resulting supermassive black hole seeds, but rather invoked comparison to simple polytropic models. Here, we systematically examine the birth, evolution and collapse of accreting non-rotating supermassive stars under accretion rates of 0.01 − 10 M ⊙ yr −1 using the stellar evolution code Kepler. Our approach includes post-Newtonian corrections to the stellar structure and an adaptive nuclear network, and can transition to following the hydrodynamic evolution of supermassive stars after they encounter the general relativistic instability. We find that this instability triggers the collapse of the star at masses of 150, 000 − 330, 000 M ⊙ for accretion rates of 0.1 − 10 M ⊙ yr −1 , and that the final mass of the star scales roughly logarithmically with the rate. The structure of the star, and thus its stability against collapse, is sensitive to the treatment of convection, and the heat content of the outer accreted envelope. Comparison with other codes suggests differences here may lead to small deviations in the evolutionary state of the star as a function of time, that worsen with accretion rate. Since the general relativistic instability leads to the immediate death of these stars, our models place an upper limit on the masses of the first quasars at birth.

In recent years, the discovery of massive quasars at $z\sim7$ has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on the discussions held at the Monash Prato Centre from November 20 to 24, 2017, during the workshop ‘Titans of the Early Universe: The Origin of the First Supermassive Black Holes’.

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.

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 ΩΓ-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.

Context. The origin of the fossil magnetic fields detected in 5 to 10% of intermediate-mass main sequence stars is still highly debated. Aims. We want to bring observational constraints to a large population of intermediate-mass pre-main sequence (PMS) stars in order to test the theory that convective-dynamo fields generated during the PMS phases of stellar evolution can occasionally relax into fossil fields on the main sequence. Methods. Using distance estimations, photometric measurements, and spectropolarimetric data from HARPSpol and ESPaDOnS of 38 intermediate-mass PMS stars, we determined fundamental stellar parameters (Teff, L and v sin i) and measured surface magnetic field characteristics (including detection limits for non-detections, and longitudinal fields and basic topologies for positive detections). Using PMS evolutionary models, we determined the mass, radius, and internal structure of these stars. We compared different PMS models to check that our determinations were not model-dependant. We then compared the magnetic characteristics of our sample accounting for their stellar parameters and internal structures. Results. We detect magnetic fields in about half of our sample. About 90% of the magnetic stars have outer convective envelopes larger than ∼25% of the stellar radii, and heavier than ∼2% of the stellar mass. Going to higher mass, we find that the magnetic incidence in intermediate-mass stars drops very quickly, within a timescale on the order of few times 0.1 Myr. Finally, we propose that intermediate-mass T Tauri stars with large convective envelopes, close to the fully convective limit, have complex fields and that their dipole component strengths may decrease as the sizes of their convective envelopes decrease, similar to lower-mass T Tauri stars.

Context. 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. Aims. 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. Methods. We used the Geneva stellar evolution code to compute the evolution of stars with initial masses of 1 and 1.5 M⊙ 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. Results. 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. In the plots of initial orbital distance versus planetary mass, we show the regions that lead to engulfment or any significant changes in the orbit. As a result of orbital evolution, a region near the star can become devoid of planets after the PMS phase. We call this zone the planet desert, and its extent depends sensitively on stellar rotation. An examination of the planet distribution as a function of distance to the host star and mass can provide constraints on current computations.