Massive young stellar object (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally-balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (ALMA) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at ~a few 100 au from the protostar can be resolved by interferometric ALMA Cycle 7 C43-8 and C43-10 observations at band 6 (1.2 mm), using a maximal 0.015" beam angular resolution and at least 10-30 min exposure time for sources at distances of 1-2 kpc. Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability.

In recent years, the discovery of massive quasars at z~7 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 discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

To explain the observed population of supermassive black holes at z~7, very massive seed black holes or, alternatively, super-Eddington scenarios are needed to reach final masses of the order of 10^9 solar masses. A popular explanation for massive seeds has been the direct collapse model, which predicts the formation of a single massive object due to the direct collapse of a massive gas cloud. Simulations over the last years have however shown that such a scenario is very difficult to achieve. A realistic model of black hole formation should therefore take fragmentation into account, and consider the interaction between stellar-dynamical and gas-dynamical processes. We present here numerical simulations pursued with the AMUSE code, employing an approximate treatment of the gas. Based on these simulations, we show that very massive black holes of 10^4-10^5 solar masses may form depending on the gas supply and the accretion onto the protostars.

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 time-scale of 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.

Models aiming to explain the formation of massive black hole seeds, and in particular the direct collapse scenario, face substantial difficulties. These are rooted in rather ad hoc and fine-tuned initial conditions, such as the simultaneous requirements of extremely low metallicities and strong radiation backgrounds. Here we explore a modification of such scenarios where a massive primordial star cluster is initially produced. Subsequent stellar collisions give rise to the formation of massive (10^4 - 10^5 solar mass) objects. Our calculations demonstrate that the interplay between stellar dynamics, gas accretion and protostellar evolution is particularly relevant. Gas accretion onto the protostars enhances their radii, resulting in an enhanced collisional cross section. We show that the fraction of collisions can increase from 0.1-1% of the initial population to about 10% when compared to gas-free models or models of protostellar clusters in the local Universe. We conclude that very massive objects can form in spite of initial fragmentation, making the first massive protostellar clusters viable candidate birth places for observed supermassive black holes.

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 solar masses per year, 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 solar masses for accretion rates of 0.1-10 solar masses per year, 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.

We discuss the physics and the evolution of a typical massive star passing through an evolutionary stage similar to that of Betelgeuse. After a brief introduction recalling various observed parameters of Betelgeuse, we discuss the Pre-Main-Sequence phase (PMS), the Main-Sequence (MS) phase, the physics governing the duration of the first crossing of the HR diagram, the red supergiant stage (RSG), the post-red supergiant phases and the final fate of solar metallicity stars with masses between 9 and 25 M$_\odot$. We examine the impact of different initial rotation and of various prescriptions for the mass loss rates during the red supergiant phase. We show that, whatever the initial rotation rate (chosen between 0 and 0.7$\times\upsilon_{\rm crit}$, $\upsilon_{\rm crit}$ being the surface equatorial velocity producing a centrifugal acceleration balancing exactly the gravity) and the mass loss rates during the RSG stage (varied between a standard value and 25 times that value), a 15 M$_\odot$ star always ends its lifetime as a RSG and explodes as a type II-P or II-L supernova.