skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert
ABSTRACT We analyse the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for $$\sim 9 \rm Myr$$, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at $$\sim 8{{\ \rm per\ cent}}$$ star formation efficiency (SFE), and the first supernova at $$8.3\, \rm Myr$$ comes too late to influence star formation significantly. The per-free-fall SFE is dynamic, accelerating from 0 per cent to $$\sim 18{{\ \rm per\ cent}}$$ before dropping quickly to <1 per cent, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief ($$\sim 1\, \mathrm{Myr}$$) compact ($$\sim 1\, \rm pc$$) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope α = −2 and a maximum mass of 55 M⊙. Stellar accretion takes $$\sim 400\, \rm kyr$$ on average, but $$\gtrsim 1\,\rm Myr$$ for >10 M⊙ stars, so massive stars finish growing latest. The fraction of stars in multiples increase as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to previous versions that neglected different feedback physics entirely. But more detailed comparison with synthetic observations will be needed to constrain the theoretical uncertainties.  more » « less
Award ID(s):
1748571 1652522
PAR ID:
10364043
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
512
Issue:
1
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 216-232
Size(s):
p. 216-232
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT It remains a major challenge to derive a theory of cloud-scale ($$\lesssim100$$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially resolved (∼100 pc) CO-to-H α flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically $$10\!-\!30\,{\rm Myr}$$, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $$\Sigma _{\rm H_2}\ge 8\,\rm {M_\odot}\,{{\rm pc}}^{-2}$$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $$\Sigma _{\rm H_2}\le 8\,\rm {M_\odot}\,{{\rm pc}}^{-2}$$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H α (75–90 per cent of the cloud lifetime), GMCs disperse within just $$1\!-\!5\,{\rm Myr}$$ once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4–10 per cent. These results show that galactic star formation is governed by cloud-scale, environmentally dependent, dynamical processes driving rapid evolutionary cycling. GMCs and H ii regions are the fundamental units undergoing these lifecycles, with mean separations of $$100\!-\!300\,{{\rm pc}}$$ in star-forming discs. Future work should characterize the multiscale physics and mass flows driving these lifecycles. 
    more » « less
  2. Abstract Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throughout collapse, star formation, and disruption. Low measured SFEs ($${\ll} 1\hbox{ per cent}$$) are ‘real’ but correspond to early stages; the true ‘per-freefall’ SFE where most stars actually form can be much larger. Very high ($${\gg} 10\hbox{ per cent}$$) values are often artificially enhanced by rapid gas dispersal. Simulations including stellar feedback reproduce observed GMC-scale SFEs, but simulations without feedback produce 20× larger SFEs. Radiative feedback dominates among mechanisms simulated. An anticorrelation of SFE with cloud mass is shown to be an observational artefact. We also explore individual dense ‘clumps’ within GMCs and show that (with feedback) their bulk properties agree well with observations. Predicted SFEs within the dense clumps are ∼2× larger than observed, possibly indicating physics other than feedback from massive (main-sequence) stars is needed to regulate their collapse. 
    more » « less
  3. ABSTRACT We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star-forming galaxy at z ∼ 2 under the presence of strong quasar winds. Our simulations include multiphase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disc winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-centre clump with stellar mass $${\rm M}_{\star }\sim 10^{7}\, {\rm M}_{\odot }$$, effective radius $${\rm R}_{\rm 1/2\, \rm Clump}\sim 20\, {\rm pc}$$, and surface density $$\Sigma _{\star } \sim 10^{4}\, {\rm M}_{\odot }\, {\rm pc}^{-2}$$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches $${\rm SFR} \sim 50\, {\rm M}_{\odot }\, {\rm yr}^{-1}$$ and $$\Sigma _{\rm SFR} \sim 10^{4}\, {\rm M}_{\odot }\, {\rm yr}^{-1}\, {\rm kpc}^{-2}$$, converting most of the progenitor gas cloud into stars in ∼2 Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time (∼35 Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters. 
    more » « less
  4. null (Ed.)
    ABSTRACT Feedback from massive stars plays a key role in molecular cloud evolution. After the onset of star formation, the young stellar population is exposed by photoionization, winds, supernovae, and radiation pressure from massive stars. Recent observations of nearby galaxies have provided the evolutionary timeline between molecular clouds and exposed young stars, but the duration of the embedded phase of massive star formation is still ill-constrained. We measure how long massive stellar populations remain embedded within their natal cloud, by applying a statistical method to six nearby galaxies at $$20{-}100~\mbox{$${\rm ~pc}$$}$$ resolution, using CO, Spitzer 24$$\rm \, \mu m$$, and H α emission as tracers of molecular clouds, embedded star formation, and exposed star formation, respectively. We find that the embedded phase (with CO and 24$$\rm \, \mu m$$ emission) lasts for 2−7 Myr and constitutes $$17{-}47{{\ \rm per\ cent}}$$ of the cloud lifetime. During approximately the first half of this phase, the region is invisible in H α, making it heavily obscured. For the second half of this phase, the region also emits in H α and is partially exposed. Once the cloud has been dispersed by feedback, 24$$\rm \, \mu m$$ emission no longer traces ongoing star formation, but remains detectable for another 2−9 Myr through the emission from ambient CO-dark gas, tracing star formation that recently ended. The short duration of massive star formation suggests that pre-supernova feedback (photoionization and winds) is important in disrupting molecular clouds. The measured time-scales do not show significant correlations with environmental properties (e.g. metallicity). Future JWST observations will enable these measurements routinely across the nearby galaxy population. 
    more » « less
  5. Abstract Molecular clouds are supported by turbulence and magnetic fields, but quantifying their influence on cloud life cycle and star formation efficiency (SFE) remains an open question. We perform radiation magnetohydrodynamic simulations of star-forming giant molecular clouds (GMCs) with UV radiation feedback, in which the propagation of UV radiation via ray tracing is coupled to hydrogen photochemistry. We consider 10 GMC models that vary in either initial virial parameter (1 ≤ α vir,0 ≤ 5) or dimensionless mass-to-magnetic flux ratio (0.5 ≤ μ Φ,0 ≤ 8 and ∞ ); the initial mass 10 5 M ⊙ and radius 20 pc are fixed. Each model is run with five different initial turbulence realizations. In most models, the duration of star formation and the timescale for molecular gas removal (primarily by photoevaporation) are 4–8 Myr. Both the final SFE ( ε * ) and time-averaged SFE per freefall time ( ε ff ) are reduced by strong turbulence and magnetic fields. The median ε * ranges between 2.1% and 9.5%. The median ε ff ranges between 1.0% and 8.0%, and anticorrelates with α vir,0 , in qualitative agreement with previous analytic theory and simulations. However, the time-dependent α vir ( t ) and ε ff,obs ( t ) based on instantaneous gas properties and cluster luminosity are positively correlated due to rapid evolution, making observational validation of star formation theory difficult. Our median ε ff,obs ( t ) ≈ 2% is similar to observed values. We show that the traditional virial parameter estimates the true gravitational boundedness within a factor of 2 on average, but neglect of magnetic support and velocity anisotropy can sometimes produce large departures from traditional virial parameter estimates. Magnetically subcritical GMCs are unlikely to represent sites of massive star formation given their unrealistic columnar outflows, prolonged lifetime, and low escape fraction of radiation. 
    more » « less