On the origin of magnetic fields in stars II: The effect of numerical resolution

This is the author accepted manuscript. The final version is available from Oxford University Press via the DOI in this recordAre the kG-strength magnetic fields observed in young stars a fossil field left over from their formation or are they generated by a dynamo? Our previous numerical study concluded that magnetic fields must originate by a dynamo process. Here, we continue that investigation by performing even higher numerical resolution calculations of the gravitational collapse of a 1 Msun rotating, magnetised molecular cloud core through the first and second collapse phases until stellar densities are reached. Each model includes Ohmic resistivity, ambipolar diffusion, and the Hall effect. We test six numerical resolutions, using between 10^5 and 3.0 x 10^7 particles to model the cloud. At all but the lowest resolutions, magnetic walls form in the outer parts of the first hydrostatic core, with the maximum magnetic field strength located within the wall rather than at the centre of the core. At high resolution, this magnetic wall is disrupted by the Hall effect, producing a magnetic field with a spiral-shaped distribution of intensity. As the second collapse occurs, this field is dragged inward and grows in strength, with the maximum field strength increasing with resolution. As the second core forms, the maximum field strength exceeds 1 kG in our highest resolution simulations, and the stellar core field strength exceeds this threshold at the highest resolution. Our resolution study suggests that kG-strength magnetic fields may be implanted in low-mass stars during their formation, and may persist over long timescales given that the diffusion timescale for the magnetic field exceeds the age of the Universe.European Union FP7University of St AndrewsAustralian Research Council (ARC

Infall of gas as the formation mechanism of stars up to 20 times more massive than the Sun

Theory predicts and observations confirm that low-mass stars (like the Sun) in their early life grow by accreting gas from the surrounding material. But for stars ~ 10 times more massive than the Sun (~10 M_sun), the powerful stellar radiation is expected to inhibit accretion and thus limit the growth of their mass. Clearly, stars with masses >10 M_sun exist, so there must be a way for them to form. The problem may be solved by non-spherical accretion, which allows some of the stellar photons to escape along the symmetry axis where the density is lower. The recent detection of rotating disks and toroids around very young massive stars has lent support to the idea that high-mass (> 8 M_sun) stars could form in this way. Here we report observations of an ammonia line towards a high-mass star forming region. We conclude from the data that the gas is falling inwards towards a very young star of ~20 M_sun, in line with theoretical predictions of non-spherical accretion.Comment: 11 pages, 2 figure

A Minimum Column Density of 1 g cm^-2 for Massive Star Formation

Massive stars are very rare, but their extreme luminosities make them both the only type of young star we can observe in distant galaxies and the dominant energy sources in the universe today. They form rarely because efficient radiative cooling keeps most star-forming gas clouds close to isothermal as they collapse, and this favors fragmentation into stars <~1 Msun. Heating of a cloud by accreting low-mass stars within it can prevent fragmentation and allow formation of massive stars, but what properties a cloud must have to form massive stars, and thus where massive stars form in a galaxy, has not yet been determined. Here we show that only clouds with column densities >~ 1 g cm^-2 can avoid fragmentation and form massive stars. This threshold, and the environmental variation of the stellar initial mass function (IMF) that it implies, naturally explain the characteristic column densities of massive star clusters and the difference between the radial profiles of Halpha and UV emission in galactic disks. The existence of a threshold also implies that there should be detectable variations in the IMF with environment within the Galaxy and in the characteristic column densities of massive star clusters between galaxies, and that star formation rates in some galactic environments may have been systematically underestimated.Comment: Accepted for publication in Nature; Nature manuscript style; main text: 14 pages, 3 figures; supplementary text: 8 pages, 1 figur

Radiation hydrodynamic simulations of massive star formation via gravitationally trapped H II regions - spherically symmetric ionized accretion flows

This paper investigates the gravitational trapping of H II regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealized spherically symmetric radiation hydrodynamical simulations of the early evolution of H II regions including the gravity of the central source. As with analytic steady-state solutions of spherically symmetric ionized Bondi accretion flows, we find gravitationally trapped H II regions with accretion through the ionization front on to the source. We found that, for a constant ionizing luminosity, fluctuations in the ionization front are unstable. This instability only occurs in this spherically symmetric accretion geometry. In the context of massive star formation, the ionizing luminosity increases with time as the source accretes mass. The maximum radius of the recurring H II region increases on the accretion time-scale until it reaches the sonic radius, where the infall velocity equals the sound speed of the ionized gas, after which it enters a pressure-driven expansion phase. This expansion prevents accretion of gas through the ionization front, the accretion rate on to the star decreases to zero, and it stops growing from accretion. Because of the time required for any significant change in stellar mass and luminosity through accretion our simulations keep both mass and luminosity constant and follow the evolution from trapped to expanding in a piecewise manner. Implications of this evolution of H II regions include a continuation of accretion of material on to forming stars for a period after the star starts to emit ionizing radiation, and an extension of the lifetime of ultracompact H II regions

ALMA observations of the Extended Green Object G19.01−0.03 – I. A Keplerian disc in a massive protostellar system

Using the Atacama Large Millimetre/submillimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA), we observed the Extended Green Object (EGO) G19.01-0.03 with sub-arcsec resolution from 1.05 mm to 5.01 cm wavelengths. Our 0.4 arcsec 1600 au angular resolution ALMA observations reveal a velocity gradient across the millimetre core MM1, oriented perpendicular to the previously known bipolar molecular outflow, which is consistently traced by 20 lines of 8 molecular species with a range of excitation temperatures, including complex organic molecules (COMs). Kinematic modelling shows the data are well described by models that include a disc in Keplerian rotation and infall, with an enclosed mass of 40-70 M (within a 2000 au outer radius) for a disc inclination angle of i = 40, of which 5.4-7.2 M is attributed to the disc. Our new VLA observations show that the 6.7 GHz Class II methanol masers associated with MM1 fo a partial ellipse, consistent with an inclined ring, with a velocity gradient consistent with that of the theal gas. The disc-to-star mass ratio suggests the disc is likely to be unstable and may be fragmenting into as-yet-undetected low-mass stellar companions. Modelling the centimetre-millimetre spectral energy distribution of MM1 shows the ALMA 1.05 mm continuum emission is dominated by dust, whilst a free-free component, interpreted as a hypercompact H ii region, is required to explain the VLA 5 cm emission. The high enclosed mass derived for a source with a moderate bolometric luminosity (104L) suggests that the MM1 disc may feed an unresolved high-mass binary system

A Substantial Population of Low Mass Stars in Luminous Elliptical Galaxies

The stellar initial mass function (IMF) describes the mass distribution of stars at the time of their formation and is of fundamental importance for many areas of astrophysics. The IMF is reasonably well constrained in the disk of the Milky Way but we have very little direct information on the form of the IMF in other galaxies and at earlier cosmic epochs. Here we investigate the stellar mass function in elliptical galaxies by measuring the strength of the Na I doublet and the Wing-Ford molecular FeH band in their spectra. These lines are strong in stars with masses <0.3 Msun and weak or absent in all other types of stars. We unambiguously detect both signatures, consistent with previous studies that were based on data of lower signal-to-noise ratio. The direct detection of the light of low mass stars implies that they are very abundant in elliptical galaxies, making up >80% of the total number of stars and contributing >60% of the total stellar mass. We infer that the IMF in massive star-forming galaxies in the early Universe produced many more low mass stars than the IMF in the Milky Way disk, and was probably slightly steeper than the Salpeter form in the mass range 0.1 - 1 Msun.Comment: To appear in Natur

A disk of dust and molecular gas around a high-mass protostar

The processes leading to the birth of low-mass stars such as our Sun have been well studied, but the formation of high-mass (> 8 x Sun's mass) stars has heretofore remained poorly understood. Recent observational studies suggest that high-mass stars may form in essentially the same way as low-mass stars, namely via an accretion process, instead of via merging of several low-mass (< 8 Msun) stars. However, there is as yet no conclusive evidence. Here, we report the discovery of a flattened disk-like structure observed at submillimeter wavelengths, centered on a massive 15 Msun protostar in the Cepheus-A region. The disk, with a radius of about 330 astronomical units (AU) and a mass of 1 to 8 Msun, is detected in dust continuum as well as in molecular line emission. Its perpendicular orientation to, and spatial coincidence with the central embedded powerful bipolar radio jet, provides the best evidence yet that massive stars form via disk accretion in direct analogy to the formation of low-mass stars

The stellar and sub-stellar IMF of simple and composite populations

The current knowledge on the stellar IMF is documented. It appears to become top-heavy when the star-formation rate density surpasses about 0.1Msun/(yr pc^3) on a pc scale and it may become increasingly bottom-heavy with increasing metallicity and in increasingly massive early-type galaxies. It declines quite steeply below about 0.07Msun with brown dwarfs (BDs) and very low mass stars having their own IMF. The most massive star of mass mmax formed in an embedded cluster with stellar mass Mecl correlates strongly with Mecl being a result of gravitation-driven but resource-limited growth and fragmentation induced starvation. There is no convincing evidence whatsoever that massive stars do form in isolation. Various methods of discretising a stellar population are introduced: optimal sampling leads to a mass distribution that perfectly represents the exact form of the desired IMF and the mmax-to-Mecl relation, while random sampling results in statistical variations of the shape of the IMF. The observed mmax-to-Mecl correlation and the small spread of IMF power-law indices together suggest that optimally sampling the IMF may be the more realistic description of star formation than random sampling from a universal IMF with a constant upper mass limit. Composite populations on galaxy scales, which are formed from many pc scale star formation events, need to be described by the integrated galactic IMF. This IGIMF varies systematically from top-light to top-heavy in dependence of galaxy type and star formation rate, with dramatic implications for theories of galaxy formation and evolution.Comment: 167 pages, 37 figures, 3 tables, published in Stellar Systems and Galactic Structure, Vol.5, Springer. This revised version is consistent with the published version and includes additional references and minor additions to the text as well as a recomputed Table 1. ISBN 978-90-481-8817-

The Role of Column Density in the Formation of Stars and Black Holes

The stellar mass in disk galaxies scales approximately with the fourth power of the rotation velocity, and the masses of the central black holes in galactic nuclei scale approximately with the fourth power of the bulge velocity dispersion. It is shown here that these relations can be accounted for if, in a forming galaxy with an isothermal mass distribution, gas with a column density above about 8 Msun/pc^2 goes into stars while gas with a column density above about 2 g/cm^2 (10^4 Msun/pc^2) goes into a central black hole. The lower critical value is close to the column density of about 10 Msun/pc^2 at which atomic gas becomes molecular, and the upper value agrees approximately with the column density of about 1 g/cm^2 at which the gas becomes optically thick to its cooling radiation. These results are plausible because molecule formation is evidently necessary for star formation, and because the onset of a high optical depth in a galactic nucleus may suppress continuing star formation and favour the growth of a central black hole.Comment: Accepted by Nature Physic

A Triple Protostar System Formed via Fragmentation of a Gravitationally Unstable Disk

Binary and multiple star systems are a frequent outcome of the star formation process, and as a result, almost half of all sun-like stars have at least one companion star. Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large scale fragmentation of turbulent gas cores and filaments or smaller scale fragmentation of a massive protostellar disk due to gravitational instability. Observational evidence for turbulent fragmentation on scales of $>$1000~AU has recently emerged. Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems. The triple protostar system L1448 IRS3B is an ideal candidate to search for evidence of disk fragmentation. L1448 IRS3B is in an early phase of the star formation process, likely less than 150,000 years in age, and all protostars in the system are separated by $<$200~AU. Here we report observations of dust and molecular gas emission that reveal a disk with spiral structure surrounding the three protostars. Two protostars near the center of the disk are separated by 61 AU, and a tertiary protostar is coincident with a spiral arm in the outer disk at a 183 AU separation. The inferred mass of the central pair of protostellar objects is $\sim$1 M$_{sun}$, while the disk surrounding the three protostars has a total mass of $\sim$0.30 M_{\sun}. The tertiary protostar itself has a minimum mass of $\sim$0.085 M$_{sun}$. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150~AU and 320~AU, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.Comment: Published in Nature on Oct. 27th. 24 pages, 8 figure