-
The TEMPO Survey II: Science Cases Leveraged from a Proposed 30-Day Time Domain Survey of the Orion Nebula with the Nancy Grace Roman Space Telescope
Authors:
Melinda Soares-Furtado,
Mary Anne Limbach,
Andrew Vanderburg,
John Bally,
Juliette Becker,
Anna L. Rosen,
Luke G. Bouma,
Johanna M. Vos,
Steve B. Howell,
Thomas G. Beatty,
William M. J. Best,
Anne Marie Cody,
Adam Distler,
Elena D'Onghia,
René Heller,
Brandon S. Hensley,
Natalie R. Hinkel,
Brian Jackson,
Marina Kounkel,
Adam Kraus,
Andrew W. Mann,
Nicholas T. Marston,
Massimo Robberto,
Joseph E. Rodriguez,
Jason H. Steffen
, et al. (4 additional authors not shown)
Abstract:
The TEMPO (Transiting Exosatellites, Moons, and Planets in Orion) Survey is a proposed 30-day observational campaign using the Nancy Grace Roman Space Telescope. By providing deep, high-resolution, short-cadence infrared photometry of a dynamic star-forming region, TEMPO will investigate the demographics of exosatellites orbiting free-floating planets and brown dwarfs -- a largely unexplored disco…
▽ More
The TEMPO (Transiting Exosatellites, Moons, and Planets in Orion) Survey is a proposed 30-day observational campaign using the Nancy Grace Roman Space Telescope. By providing deep, high-resolution, short-cadence infrared photometry of a dynamic star-forming region, TEMPO will investigate the demographics of exosatellites orbiting free-floating planets and brown dwarfs -- a largely unexplored discovery space. Here, we present the simulated detection yields of three populations: extrasolar moon analogs orbiting free-floating planets, exosatellites orbiting brown dwarfs, and exoplanets orbiting young stars. Additionally, we outline a comprehensive range of anticipated scientific outcomes accompanying such a survey. These science drivers include: obtaining observational constraints to test prevailing theories of moon, planet, and star formation; directly detecting widely separated exoplanets orbiting young stars; investigating the variability of young stars and brown dwarfs; constraining the low-mass end of the stellar initial mass function; constructing the distribution of dust in the Orion Nebula and mapping evolution in the near-infrared extinction law; mapping emission features that trace the shocked gas in the region; constructing a dynamical map of Orion members using proper motions; and searching for extragalactic sources and transients via deep extragalactic observations reaching a limiting magnitude of $m_{AB}=29.7$\,mag (F146 filter).
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
Gamma-ray Emission from a Young Star Cluster in the Star-Forming Region RCW 38
Authors:
Paarmita Pandey,
Laura A. Lopez,
Anna L. Rosen,
Todd A. Thompson,
Tim Linden,
Ian Blackstone,
Lachlan Lancaster
Abstract:
We report the detection of gamma-ray emission near the young Milky Way star cluster ($\approx$ 0.5 Myr old) in the star-forming region RCW 38. Using 15 years of data from the Fermi-LAT, we find a significant ($ σ> 22$) detection coincident with the cluster, producing a total $γ$-ray luminosity of $L_γ = (2.66\pm 0.92) \times 10^{34}$ erg s$^{-1}$ adopting a power-law spectral model (…
▽ More
We report the detection of gamma-ray emission near the young Milky Way star cluster ($\approx$ 0.5 Myr old) in the star-forming region RCW 38. Using 15 years of data from the Fermi-LAT, we find a significant ($ σ> 22$) detection coincident with the cluster, producing a total $γ$-ray luminosity of $L_γ = (2.66\pm 0.92) \times 10^{34}$ erg s$^{-1}$ adopting a power-law spectral model ($Γ= 2.34\pm0.04$) in the 0.1 $-$ 500 GeV band. We estimate the total wind power to be $7 \times 10^{36}$ erg s$^{-1}$, corresponding to a CR acceleration efficiency of $η_{\rm CR} \simeq 0.4$ for a diffusion coefficient consistent with the local interstellar medium of $D = 10^{28}$ cm$^{2}$ s$^{-1}$. Alternatively, the $γ$-ray luminosity could also account for a lower acceleration efficiency of 0.1 if the diffusion coefficient in the star-forming region is smaller, $D\simeq 2.5\times10^{27}\,{\rm cm^2\,\,s^{-1}}$. We compare the hot-gas pressure from Chandra X-ray analysis to the CR pressure and find the former is four orders of magnitude greater, suggesting that the CR pressure is not dynamically important relative to the stellar wind feedback. As RCW 38 is too young for supernovae to have occurred, the high CR acceleration efficiency in RCW 38 demonstrates that stellar winds may be an important source of Galactic cosmic rays.
△ Less
Submitted 29 April, 2024;
originally announced April 2024.
-
FORGE'd in FIRE III: The IMF in Quasar Accretion Disks from STARFORGE
Authors:
Philip F. Hopkins,
Michael Y. Grudic,
Kyle Kremer,
Stella S. R. Offner,
David Guszejnov,
Anna L. Rosen
Abstract:
Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10…
▽ More
Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10 pc from the SMBH. We show it is critical to include physics often previously neglected, including magnetic fields, radiation, and (proto)stellar feedback. Closer to the SMBH, star formation is suppressed, but the (rare) stars that do form exhibit top-heavy IMFs. Stars can form only in special locations (e.g. magnetic field switches) in the outer QAD. Protostars accrete their natal cores rapidly but then dynamically decouple from the gas and wander, ceasing accretion on timescales ~100 yr. Their jets control initial core accretion, but the ejecta are swept up into the larger-scale QAD flow without much dynamical effect. The strong tidal environment strongly suppresses common-core multiplicity. The IMF shape depends sensitively on un-resolved dynamics of protostellar disks (PSDs), as the global dynamical times can become incredibly short ($\ll$ yr) and tidal fields are incredibly strong, so whether PSDs can efficiently transport angular momentum or fragment catastrophically at $\lesssim 10$ au scales requires novel PSD simulations to properly address. Most analytic IMF models and analogies with planet formation in PSDs fail qualitatively to explain the simulation IMFs, though we discuss a couple of viable models.
△ Less
Submitted 28 August, 2024; v1 submitted 11 April, 2024;
originally announced April 2024.
-
Detection of Diffuse Hot Gas Around the Young, Potential Superstar Cluster H72.97-69.39
Authors:
Trinity L. Webb,
Jennifer A. Rodriguez,
Laura A. Lopez,
Anna L. Rosen,
Lachlan Lancaster,
Omnarayani Nayak,
Anna F. McLeod,
Paarmita Pandey,
Grace M. Olivier
Abstract:
We present the first Chandra X-ray observations of H72.97-69.39, a highly-embedded, potential super-star cluster (SSC) in its infancy located in the star-forming complex N79 of the Large Magellanic Cloud. We detect particularly hard, diffuse X-ray emission that is coincident with the young stellar object (YSO) clusters identified with JWST, and the hot gas fills cavities in the dense gas mapped by…
▽ More
We present the first Chandra X-ray observations of H72.97-69.39, a highly-embedded, potential super-star cluster (SSC) in its infancy located in the star-forming complex N79 of the Large Magellanic Cloud. We detect particularly hard, diffuse X-ray emission that is coincident with the young stellar object (YSO) clusters identified with JWST, and the hot gas fills cavities in the dense gas mapped by ALMA. The X-ray spectra are best fit with either a thermal plasma or power-law model, and assuming the former, we show that the X-ray luminosity of L_X = (1.5 +- 0.3)e34 erg/s is a factor of ~20 below the expectation for a fully-confined wind bubble. Our results suggest that stellar wind feedback produces diffuse hot gas in the earliest stages of massive star cluster formation and that wind energy can be lost quickly via either turbulent mixing followed by radiative cooling or by physical leakage.
△ Less
Submitted 21 February, 2024;
originally announced February 2024.
-
Stellar Populations in STARFORGE: The Origin and Evolution of Star Clusters and Associations
Authors:
Juan P. Farias,
Stella S. R. Offner,
Michael Y. Grudić,
Dávid Guszejnov,
Anna L. Rosen
Abstract:
Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the STARFORGE simulation suite -- a set of radiation-magne…
▽ More
Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the STARFORGE simulation suite -- a set of radiation-magnetohydrodynamic simulations of star-forming turbulent clouds that include all key stellar feedback processes inherent to star formation. We use Nbody6++GPU to follow the evolution of the young stellar systems after gas removal. We use HDBSCAN to define stellar groups and analyze the stellar kinematics to identify the true bound star clusters. The conditions modeled by the simulations, i.e., global cloud surface densities below 0.15 g cm$^{-2}$,, star formation efficiencies below 15%, and gas expulsion timescales shorter than a free fall time, primarily produce expanding stellar associations and small clusters. The largest star clusters, which have $\sim$1000 bound members, form in the densest and lowest velocity dispersion clouds, representing $\sim$32 and 39% of the stars in the simulations, respectively. The cloud's early dynamical state plays a significant role in setting the classical star formation efficiency versus bound fraction relation. All stellar groups follow a narrow mass-velocity dispersion power law relation at 10 Myr with a power law index of 0.21. This correlation result in a distinct mass-size relationship for bound clusters. We also provide valuable constraints on the gas dispersal timescale during the star formation process and analyze the implications for the formation of bound systems.
△ Less
Submitted 20 September, 2023;
originally announced September 2023.
-
What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates
Authors:
Sabrina M. Appel,
Blakesley Burkhart,
Vadim A. Semenov,
Christoph Federrath,
Anna L. Rosen,
Jonathan C. Tan
Abstract:
We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycl…
▽ More
We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high and low density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.
△ Less
Submitted 18 January, 2023;
originally announced January 2023.
-
A Multiwavelength Study of the Massive Colliding Wind Binary WR 20a: A Possible Progenitor for Fast-Spinning LIGO Binary Black Hole Mergers
Authors:
Grace M. Olivier,
Laura A. Lopez,
Katie Auchettl,
Anna L. Rosen,
Aldo Batta,
Kathryn F. Neugent,
Enrico Ramirez-Ruiz,
Tharindu Jayasinghe,
Patrick J. Vallely,
Dominick M. Rowan
Abstract:
WR 20a is the most massive close-in binary known in our Galaxy. It is composed of two $\approx$80 M$_\odot$ Wolf-Rayet stars with a short period of $\approx$3.7 days in the open cluster Westerlund 2. As such, WR 20a presents us with a unique laboratory for studying the currently uncertain physics of binary evolution and compact object formation as well as for studying the wind collision region in…
▽ More
WR 20a is the most massive close-in binary known in our Galaxy. It is composed of two $\approx$80 M$_\odot$ Wolf-Rayet stars with a short period of $\approx$3.7 days in the open cluster Westerlund 2. As such, WR 20a presents us with a unique laboratory for studying the currently uncertain physics of binary evolution and compact object formation as well as for studying the wind collision region in an massive eclipsing binary system. We use deep Chandra observations of WR 20a to study the time variability of the wind collision region between the two Wolf-Rayet stars and are able to produce an X-ray light curve covering $\approx$2/3 of its orbital period. We find that the X-ray light curve is asymmetric because the flux of one peak is 2.5$σ$ larger than the flux of the other peak. This asymmetry could be caused by asymmetric mass-loss from the two stars or by the lopsidedness of the wind collision region due to the unusually fast rotation of the system. The X-ray light curve is also shifted in phase space when compared to the optical light curves measured by TESS and ASAS-SN. Additionally, we explore the ultimate fate of this system by modeling the resultant binary black hole merger expected at the end of the two stars' lives. We conclude that this system will evolve to be a representative of the sub-population of LIGO progenitors of fast-spinning binary black hole merger events.
△ Less
Submitted 5 December, 2022;
originally announced December 2022.
-
The TEMPO Survey I: Predicting Yields of the Transiting Exosatellites, Moons, and Planets from a 30-day Survey of Orion with the Nancy Grace Roman Space Telescope
Authors:
Mary Anne Limbach,
Melinda Soares-Furtado,
Andrew Vanderburg,
William M. J. Best,
Ann Marie Cody,
Elena D'Onghia,
René Heller,
Brandon S. Hensley,
Marina Kounkel,
Adam Kraus,
Andrew W. Mann,
Massimo Robberto,
Anna L. Rosen,
Richard Townsend,
Johanna M. Vos
Abstract:
We present design considerations for the Transiting Exosatellites, Moons, and Planets in Orion (TEMPO) Survey with the Nancy Grace Roman Space Telescope. This proposed 30-day survey is designed to detect a population of transiting extrasolar satellites, moons, and planets in the Orion Nebula Cluster (ONC). The young (1-3 Myr), densely-populated ONC harbors about a thousand bright brown dwarfs (BDs…
▽ More
We present design considerations for the Transiting Exosatellites, Moons, and Planets in Orion (TEMPO) Survey with the Nancy Grace Roman Space Telescope. This proposed 30-day survey is designed to detect a population of transiting extrasolar satellites, moons, and planets in the Orion Nebula Cluster (ONC). The young (1-3 Myr), densely-populated ONC harbors about a thousand bright brown dwarfs (BDs) and free-floating planetary-mass objects (FFPs). TEMPO offers sufficient photometric precision to monitor FFPs with ${\rm M}\geq1{\rm M}_{\rm J}$ for transiting satellites. The survey is also capable of detecting FFPs down to sub-Saturn masses via direct imaging, although follow-up confirmation will be challenging. TEMPO yield estimates include 14 (3-22) exomoons/satellites transiting FFPs and 54 (8-100) satellites transiting BDs. Of this population, approximately $50\%$ of companions would be "super-Titans" (Titan to Earth mass). Yield estimates also include approximately $150$ exoplanets transiting young Orion stars, of which $>50\%$ will orbit mid-to-late M dwarfs and approximately ten will be proto-habitable zone, terrestrial ($0.1{\rm M}_{\oplus} - 5{\rm M}_{\oplus}$) exoplanets. TEMPO would provide the first census demographics of small exosatellites orbiting FFPs and BDs, while simultaneously offering insights into exoplanet evolution at the earliest stages. This detected exosatellite population is likely to be markedly different from the current census of exoplanets with similar masses (e.g., Earth-mass exosatellites that still possess H/He envelopes). Although our yield estimates are highly uncertain, as there are no known exoplanets or exomoons analogous to these satellites, the TEMPO survey would test the prevailing theories of exosatellite formation and evolution, which limit the certainty surrounding detection yields.
△ Less
Submitted 26 September, 2022;
originally announced September 2022.
-
Effects of the environment on the multiplicity properties of stars in the STARFORGE simulations
Authors:
Dávid Guszejnov,
Aman N. Raju,
Stella S. R. Offner,
Michael Y. Grudić,
Claude-André Faucher-Giguère,
Philip F. Hopkins,
Anna L. Rosen
Abstract:
Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbul…
▽ More
Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation under-predicts these values. The discrepancy is likely due to the lack of disk fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of $>1\,\mathrm{M_\odot}$ stars, however, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anti-correlates with stellar density. This relationship also explains most multiplicity variations between runs, i.e., variations in the initial conditions that increase stellar density (increased surface density, reduced turbulence) decrease multiplicity. While other parameters, such as metallicity, interstellar radiation, and geometry significantly affect the star formation history or the IMF, varying them produces no clear trend in stellar multiplicity properties.
△ Less
Submitted 28 September, 2022; v1 submitted 4 August, 2022;
originally announced August 2022.
-
Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations
Authors:
Dávid Guszejnov,
Michael Y. Grudić,
Stella S. R. Offner,
Claude-André Faucher-Giguère,
Philip F. Hopkins,
Anna L. Rosen
Abstract:
One of the key mysteries of star formation is the origin of the stellar initial mass function (IMF). The IMF is observed to be nearly universal in the Milky Way and its satellites, and significant variations are only inferred in extreme environments, such as the cores of massive elliptical galaxies. In this work we present simulations from the STARFORGE project that are the first cloud-scale RMHD…
▽ More
One of the key mysteries of star formation is the origin of the stellar initial mass function (IMF). The IMF is observed to be nearly universal in the Milky Way and its satellites, and significant variations are only inferred in extreme environments, such as the cores of massive elliptical galaxies. In this work we present simulations from the STARFORGE project that are the first cloud-scale RMHD simulations that follow individual stars and include all relevant physical processes. The simulations include detailed gas thermodynamics, as well as stellar feedback in the form of protostellar jets, stellar radiation, winds and supernovae. In this work we focus on how stellar radiation, winds and supernovae impact star-forming clouds. Radiative feedback plays a major role in quenching star formation and disrupting the cloud, however the IMF peak is predominantly set by protostellar jet physics. We find the effect of stellar winds is minor, and supernovae occur too late}to affect the IMF or quench star formation. We also investigate the effects of initial conditions on the IMF. The IMF is insensitive to the initial turbulence, cloud mass and cloud surface density, even though these parameters significantly shape the star formation history of the cloud, including the final star formation efficiency. The characteristic stellar mass depends weakly on metallicity and the interstellar radiation field. Finally, while turbulent driving and the level of magnetization strongly influences the star formation history, they only influence the high-mass slope of the IMF.
△ Less
Submitted 2 August, 2022; v1 submitted 20 May, 2022;
originally announced May 2022.
-
A Massive Star is Born: How Feedback from Stellar Winds, Radiation Pressure, and Collimated Outflows Limits Accretion onto Massive Stars
Authors:
Anna L. Rosen
Abstract:
Massive protostars attain high luminosities as they are actively accreting and the radiation pressure exerted on the gas in the star's atmosphere may launch isotropic high-velocity winds. These winds will collide with the surrounding gas producing shock-heated ($T\sim 10^7$ K) tenuous gas that adiabatically expands and pushes on the dense gas that may otherwise be accreted. We present a suite of 3…
▽ More
Massive protostars attain high luminosities as they are actively accreting and the radiation pressure exerted on the gas in the star's atmosphere may launch isotropic high-velocity winds. These winds will collide with the surrounding gas producing shock-heated ($T\sim 10^7$ K) tenuous gas that adiabatically expands and pushes on the dense gas that may otherwise be accreted. We present a suite of 3D radiation-magnetohydrodynamic simulations of the collapse of massive prestellar cores and include radiative feedback from the stellar and dust-reprocessed radiation fields, collimated outflows, and, for the first time, isotropic stellar winds to model how these processes affect the formation of massive stars. We find that winds are initially launched when the massive protostar is still accreting and its wind properties evolve as the protostar contracts to the main-sequence. Wind feedback drives asymmetric adiabatic wind bubbles that have a bipolar morphology because the dense circumstellar material pinches the expansion of the hot shock-heated gas. We term this the "wind tunnel effect." If the core is magnetized, wind feedback is less efficient at driving adiabatic wind bubbles initially because magnetic tension delays their growth. We find that wind feedback eventually quenches accretion onto $\sim$30 $\rm{M_{\rm \odot}}$ protostars that form from the collapse of the isolated cores simulated here. Hence, our results suggest that $\gtrsim$30 $\rm{M_{\rm \odot}}$ stars likely require larger-scale dynamical inflows from their host cloud to overcome wind feedback. Additionally, we discuss the implications of observing adiabatic wind bubbles with \textit{Chandra} while the massive protostars are still highly embedded.
△ Less
Submitted 1 November, 2022; v1 submitted 20 April, 2022;
originally announced April 2022.
-
Cluster assembly and the origin of mass segregation in the STARFORGE simulations
Authors:
Dávid Guszejnov,
Carleen Markey,
Stella S. R. Offner,
Michael Y. Grudić,
Claude-André Faucher-Giguère,
Anna L. Rosen,
Philip F. Hopkins
Abstract:
Stars form in dense, clustered environments, where feedback from newly formed stars eventually ejects the gas, terminating star formation and leaving behind one or more star clusters. Using the STARFORGE simulations, it is possible to simulate this process in its entirety within a molecular cloud, while explicitly evolving the gas radiation and magnetic fields and following the formation of indivi…
▽ More
Stars form in dense, clustered environments, where feedback from newly formed stars eventually ejects the gas, terminating star formation and leaving behind one or more star clusters. Using the STARFORGE simulations, it is possible to simulate this process in its entirety within a molecular cloud, while explicitly evolving the gas radiation and magnetic fields and following the formation of individual, low-mass stars. We find that individual star-formation sites merge to form ever larger structures, while still accreting gas. Thus clusters are assembled through a series of mergers. During the cluster assembly process a small fraction of stars are ejected from their clusters; we find no significant difference between the mass distribution of the ejected stellar population and that of stars inside clusters. The star-formation sites that are the building blocks of clusters start out mass segregated with one or a few massive stars at their center. As they merge the newly formed clusters maintain this feature, causing them to have mass-segregated substructures without themselves being centrally condensed. The merged clusters relax to a centrally condensed mass segregated configuration through dynamical interactions between their members, but this process does not finish before feedback expels the remaining gas from the cluster. In the simulated runs the gas-free clusters then become unbound and break up. We find that turbulent driving and a periodic cloud geometry can significantly reduce clustering and prevent gas expulsion. Meanwhile, the initial surface density and level of turbulence have little qualitative effect on cluster evolution, despite the significantly different star formation histories.
△ Less
Submitted 19 May, 2022; v1 submitted 5 January, 2022;
originally announced January 2022.
-
The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert
Authors:
Michael Y. Grudić,
Dávid Guszejnov,
Stella S. R. Offner,
Anna L. Rosen,
Aman N. Raju,
Claude-André Faucher-Giguère,
Philip F. Hopkins
Abstract:
We analyze 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…
▽ More
We analyze 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\%$ star formation efficiency (SFE), and the first supernova at $8.3 \rm Myr$ comes too late to influence star formation significantly. The per-freefall SFE is dynamic, accelerating from 0 to $\sim 18\%$ before dropping quickly to <1%, 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_\odot$. Stellar accretion takes $\sim 400 \rm kyr$ on average, but $\gtrsim 1\rm Myr$ for $>10 M_\odot$ stars, so massive stars finish growing latest. The fraction of stars in multiples increases as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to variations which neglect different feedback physics entirely. But more detailed comparison with synthetic observations is necessary to constrain the theoretical uncertainties.
△ Less
Submitted 24 March, 2022; v1 submitted 3 January, 2022;
originally announced January 2022.
-
Less wrong: a more realistic initial condition for simulations of turbulent molecular clouds
Authors:
Henry B. Lane,
Michael Y. Grudić,
Dávid Guszejnov,
Stella S. R. Offner,
Claude-André Faucher-Giguère,
Anna L. Rosen
Abstract:
Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or "stirred" turbulence with periodic boundary con…
▽ More
Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or "stirred" turbulence with periodic boundary conditions (which may not model real GMC boundary conditions). We develop and test a new GMC simulation setup (called TURBSPHERE) that combines advantages of both approaches: we continuously stir an isolated cloud to model the energy cascade from larger scales, and use a static potential to confine the gas. The resulting cloud and surrounding envelope achieve a quasi-equilibrium state with the desired hallmarks of supersonic ISM turbulence (e.g. density PDF and a $\sim k^{-2}$ velocity power spectrum), whose bulk properties can be tuned as desired. We use the final stirred state as initial conditions for star formation simulations with self-gravity, both with and without continued driving and protostellar jet feedback, respectively. We then disentangle the respective effects of the turbulent cascade, simulation geometry, external driving, and gravity/MHD boundary conditions on the resulting star formation. Without external driving, the new setup obtains results similar to previous simple spherical cloud setups, but external driving can suppress star formation considerably in the new setup. Periodic box simulations with the same dimensions and turbulence parameters form stars significantly slower, highlighting the importance of boundary conditions and the presence or absence of a global collapse mode in the results of star formation calculations.
△ Less
Submitted 27 October, 2021;
originally announced October 2021.
-
The Effects of Magnetic Fields and Outflow Feedback on the Shape and Evolution of the Density PDF in Turbulent Star-Forming Clouds
Authors:
Sabrina M. Appel,
Blakesley Burkhart,
Vadim A. Semenov,
Christoph Federrath,
Anna L. Rosen
Abstract:
Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densitie…
▽ More
Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densities and the inclusion of stellar feedback from protostellar outflows and heating produces significant time-varying deviations from a lognormal distribution at the low densities. The simulation with outflows has an excess of diffuse gas compared to the simulations without outflows, exhibits increased average sonic Mach number, and maintains a slower star formation rate over the entire duration of the run. We study the mass transfer between the diffuse gas in the lognormal peak of the PDF, the collapsing gas in the power-law tail, and the stars. We find that the mass fraction in the power-law tail is constant, such that the stars form out of the power-law gas at the same rate at which the gas from the lognormal part replenishes the power-law. We find that turbulence does not provide significant support in the dense gas associated with the power-law tail. When including outflows and magnetic fields in addition to driven turbulence, the rate of mass transfer from the lognormal to the power-law, and then to the stars, becomes significantly slower, resulting in slower star formation rates and longer depletion times.
△ Less
Submitted 27 September, 2021;
originally announced September 2021.
-
Blowing Bubbles around Intermediate-Mass Stars: Feedback from Main-Sequence Winds is not Enough
Authors:
Anna L. Rosen,
Stella S. R. Offner,
Michael M. Foley,
Laura A. Lopez
Abstract:
Numerous spherical ``shells" have been observed in young star-forming environments that host low- and intermediate-mass stars. These observations suggest that these shells may be produced by isotropic stellar wind feedback from young main-sequence stars. However, the driving mechanism for these shells remains uncertain because the momentum injected by winds is too low to explain their sizes and dy…
▽ More
Numerous spherical ``shells" have been observed in young star-forming environments that host low- and intermediate-mass stars. These observations suggest that these shells may be produced by isotropic stellar wind feedback from young main-sequence stars. However, the driving mechanism for these shells remains uncertain because the momentum injected by winds is too low to explain their sizes and dynamics due to their low mass-loss rates. However, these studies neglect how the wind kinetic energy is transferred to the ISM and instead assume it is instantly lost via radiation, suggesting that these shells are momentum-driven. Intermediate-mass stars have fast ($v_w \gtrsim 1000$ km/s) stellar winds and therefore the energy injected by winds should produce energy-driven adiabatic wind bubbles that are larger than momentum-driven wind bubbles. Here, we explore if energy-driven wind feedback can produce the observed shells by performing a series of 3D magneto-hydrodynamic simulations of wind feedback from intermediate-mass and high-mass stars that are placed in a magnetized, turbulent molecular cloud. We find that, for the high-mass stars modeled, energy-driven wind feedback produces $\sim$pc scale wind bubbles in molecular clouds that agree with the observed shell sizes but winds from intermediate-mass stars can not produce similar shells because of their lower mass-loss rates and velocities. Therefore, such shells must be driven by other feedback processes inherent to low- and intermediate-mass star formation.
△ Less
Submitted 26 July, 2021;
originally announced July 2021.
-
Dust in the Wind with Resonant Drag Instabilities: I. The Dynamics of Dust-Driven Outflows in GMCs and HII Regions
Authors:
Philip F. Hopkins,
Anna L. Rosen,
Jonathan Squire,
Georgia V. Panopoulou,
Nadine H. Soliman,
Darryl Seligman,
Ulrich P. Steinwandel
Abstract:
Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly-coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of resonant drag instabilities (RDIs) which de-couple the…
▽ More
Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly-coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of resonant drag instabilities (RDIs) which de-couple the dust and gas dynamics and could qualitatively change the nonlinear outcome of these outflows. We present the first simulations of radiation-dust driven outflows in stratified, inhomogeneous media, including explicit grain dynamics and a realistic spectrum of grain sizes and charge, magnetic fields and Lorentz forces on grains (which dramatically enhance the RDIs), Coulomb and Epstein drag forces, and explicit radiation transport allowing for different grain absorption and scattering properties. In this paper we consider conditions resembling giant molecular clouds (GMCs), HII regions, and distributed starbursts, where optical depths are modest, single-scattering effects dominate radiation-dust coupling, Lorentz forces dominate over drag on grains, and the fastest-growing RDIs are similar, such as magnetosonic and fast-gyro RDIs. These RDIs generically produce strong size-dependent dust clustering, growing nonlinear on timescales that are much shorter than the characteristic times of the outflow. The instabilities produce filamentary and plume-like or 'horsehead' nebular morphologies that are remarkably similar to observed dust structures in GMCs and HII regions. Additionally, in some cases they strongly alter the magnetic field structure and topology relative to filaments. Despite driving strong micro-scale dust clumping which leaves some gas behind, an order-unity fraction of the gas is always efficiently entrained by dust.
△ Less
Submitted 8 November, 2022; v1 submitted 9 July, 2021;
originally announced July 2021.
-
Continuity of accretion from clumps to Class 0 high-mass protostars in SDC335
Authors:
A. Avison,
G. A. Fuller,
N. Peretto,
A. Duarte-Cabral,
A. L. Rosen,
A. Traficante,
J. E. Pineda,
R. Güsten,
N. Cunningham
Abstract:
The IRDC SDC335.579-0.292 (SDC335) is a massive star-forming cloud found to be globally collapsing towards one of the most massive star forming cores in the Galaxy. SDC335 hosts three high-mass protostellar objects at early stages of their evolution and archival ALMA Cycle 0 data indicate the presence of at least one molecular outflow in the region. Observations of molecular outflows from massive…
▽ More
The IRDC SDC335.579-0.292 (SDC335) is a massive star-forming cloud found to be globally collapsing towards one of the most massive star forming cores in the Galaxy. SDC335 hosts three high-mass protostellar objects at early stages of their evolution and archival ALMA Cycle 0 data indicate the presence of at least one molecular outflow in the region. Observations of molecular outflows from massive protostellar objects allow us to estimate the accretion rates of the protostars as well as to assess the disruptive impact that stars have on their natal clouds. The aim of this work is to identify and analyse the properties of the protostellar-driven molecular outflows within SDC335 and use these outflows to help refine the properties of the protostars. We imaged the molecular outflows in SDC335 using new data from the ATCA of SiO and Class I CH$_3$OH maser emission (~3 arcsec) alongside observations of four CO transitions made with APEX and archival ALMA CO, $^{13}$CO (~1 arcsec), and HNC data. We introduced a generalised argument to constrain outflow inclination angles based on observed outflow properties. We used the properties of each outflow to infer the accretion rates on the protostellar sources driving them and to deduce the evolutionary characteristics of the sources. We identify three molecular outflows in SDC335, one associated with each of the known compact HII regions. The outflow properties show that the SDC335 protostars are in the early stages (Class 0) of their evolution, with the potential to form stars in excess of 50 M$_{\odot}$. The measured total accretion rate onto the protostars is $1.4(\pm 0.1) \times 10^{-3}$M$_{\odot}$ yr$^{-1}$, comparable to the total mass infall rate toward the cloud centre on parsec scales of 2.5$(\pm 1.0) \times 10^{-3}$M$_{\odot}$ yr$^{-1}$, suggesting a near-continuous flow of material from cloud to core scales. [abridged].
△ Less
Submitted 16 December, 2020;
originally announced December 2020.
-
Evolution of Stellar Feedback in HII Regions
Authors:
Grace M. Olivier,
Laura A. Lopez,
Anna L. Rosen,
Omnarayani Nayak,
Megan Reiter,
Mark R. Krumholz,
Alberto D. Bolatto
Abstract:
Stellar feedback is needed to produce realistic giant molecular clouds (GMCs) and galaxies in simulations, but due to limited numerical resolution, feedback must be implemented using subgrid models. Observational work is an important means to test and anchor these models, but limited studies have assessed the relative dynamical role of multiple feedback modes, particularly at the earliest stages o…
▽ More
Stellar feedback is needed to produce realistic giant molecular clouds (GMCs) and galaxies in simulations, but due to limited numerical resolution, feedback must be implemented using subgrid models. Observational work is an important means to test and anchor these models, but limited studies have assessed the relative dynamical role of multiple feedback modes, particularly at the earliest stages of expansion when HII regions are still deeply embedded. In this paper, we use multiwavelength (radio, infrared, and X-ray) data to measure the pressures associated with direct radiation ($P_{\rm dir}$), dust-processed radiation ($P_{\rm IR}$), photoionization heating ($P_{\rm HII}$), and shock-heating from stellar winds ($P_{\rm X}$) in a sample of 106 young, resolved HII regions with radii $\lesssim$0.5 pc to determine how stellar feedback drives their expansion. We find that the $P_{\rm IR}$ dominates in 84% of the regions and that the median $P_{\rm dir}$ and $P_{\rm HII}$ are smaller than the median $P_{\rm IR}$ by factors of $\approx 6$ and $\approx 9$, respectively. Based on the radial dependences of the pressure terms, we show that HII regions transition from $P_{\rm IR}$-dominated to $P_{\rm HII}$-dominated at radii of $\sim$3 pc. We find a median trapping factor of $f_{\rm trap} \sim$ 8 without any radial dependence for the sample, suggesting this value can be adopted in sub-grid feedback models. Moreover, we show that the total pressure is greater than the gravitational pressure in the majority of our sample, indicating that the feedback is sufficient to expel gas from the regions.
△ Less
Submitted 8 December, 2020; v1 submitted 21 September, 2020;
originally announced September 2020.
-
The Role of Outflows, Radiation Pressure, and Magnetic Fields in Massive Star Formation
Authors:
Anna L. Rosen,
Mark R. Krumholz
Abstract:
Stellar feedback in the form of radiation pressure and magnetically-driven collimated outflows may limit the maximum mass that a star can achieve and affect the star-formation efficiency of massive pre-stellar cores. Here we present a series of 3D adaptive mesh refinement radiation-magnetohydrodynamic simulations of the collapse of initially turbulent, massive pre-stellar cores. Our simulations in…
▽ More
Stellar feedback in the form of radiation pressure and magnetically-driven collimated outflows may limit the maximum mass that a star can achieve and affect the star-formation efficiency of massive pre-stellar cores. Here we present a series of 3D adaptive mesh refinement radiation-magnetohydrodynamic simulations of the collapse of initially turbulent, massive pre-stellar cores. Our simulations include radiative feedback from both the direct stellar and dust-reprocessed radiation fields, and collimated outflow feedback from the accreting stars. We find that protostellar outflows punches holes in the dusty circumstellar gas along the star's polar directions, thereby increasing the size of optically thin regions through which radiation can escape. Precession of the outflows as the star's spin axis changes due to the turbulent accretion flow further broadens the outflow, and causes more material to be entrained. Additionally, the presence of magnetic fields in the entrained material leads to broader entrained outflows that escape the core. We compare the injected and entrained outflow properties and find that the entrained outflow mass is a factor of $\sim$3 larger than the injected mass and the momentum and energy contained in the entrained material are $\sim$25% and $\sim$5% of the injected momentum and energy, respectively. As a result, we find that, when one includes both outflows and radiation pressure, the former are a much more effective and important feedback mechanism, even for massive stars with significant radiative outputs.
△ Less
Submitted 8 June, 2020;
originally announced June 2020.
-
Zooming in on Individual Star Formation: Low- and High-mass Stars
Authors:
Anna L. Rosen,
Stella S. R. Offner,
Sarah I. Sadavoy,
Asmita Bhandare,
Enrique Vázquez-Semadeni,
Adam Ginsburg
Abstract:
Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and…
▽ More
Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and size scales. During the last decade a variety of technological and computing advances have transformed our understanding of star formation through the use of multi-wavelength observations, large scale observational surveys, and multi-physics multi-dimensional numerical simulations. Additionally, the use of synthetic observations of simulations have provided a useful tool to interpret observational data and evaluate the importance of various physical processes on different scales in star formation. Here, we review these recent advancements in both high- ($M \gtrsim 8 \, M_{\rm \odot}$) and low-mass star formation.
△ Less
Submitted 15 May, 2020;
originally announced May 2020.
-
Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity
Authors:
Paul C. Duffell,
Daniel D'Orazio,
Andrea Derdzinski,
Zoltan Haiman,
Andrew MacFadyen,
Anna L. Rosen,
Jonathan Zrake
Abstract:
Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sen…
▽ More
Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $α= 0.03$ to $α= 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.
△ Less
Submitted 31 July, 2020; v1 submitted 13 November, 2019;
originally announced November 2019.
-
Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores
Authors:
Anna L. Rosen,
Pak Shing Li,
Qizhou Zhang,
Blakesley Burkhart
Abstract:
Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subv…
▽ More
Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.
△ Less
Submitted 5 November, 2019; v1 submitted 26 February, 2019;
originally announced February 2019.
-
unyt: Handle, manipulate, and convert data with units in Python
Authors:
Nathan J. Goldbaum,
John A. ZuHone,
Matthew J. Turk,
Kacper Kowalik,
Anna L. Rosen
Abstract:
Software that processes real-world data or that models a physical system must have some way of managing units. While simple approaches like the understood convention that all data are in a unit system (such as the MKS SI unit system) do work in practice, they are fraught with possible sources of error both by developers and users of the software. In this paper we present unyt, a Python library bas…
▽ More
Software that processes real-world data or that models a physical system must have some way of managing units. While simple approaches like the understood convention that all data are in a unit system (such as the MKS SI unit system) do work in practice, they are fraught with possible sources of error both by developers and users of the software. In this paper we present unyt, a Python library based on NumPy and SymPy for handling data that has units. It is designed both to aid quick interactive calculations and to be tightly integrated into a larger Python application or library. We compare unyt with two other Python libraries for handling units, Pint and astropy.units, and find that unyt is faster, has higher test coverage, and has fewer lines of code.
△ Less
Submitted 1 October, 2018; v1 submitted 6 June, 2018;
originally announced June 2018.
-
An Unstable Truth: How Massive Stars get their Mass
Authors:
Anna L. Rosen,
Mark R. Krumholz,
Christopher F. McKee,
Richard I. Klein
Abstract:
The pressure exerted by massive stars' radiation fields is an important mechanism regulating their formation. Detailed simulation of massive star formation therefore requires an accurate treatment of radiation. However, all published simulations have either used a diffusion approximation of limited validity; have only been able to simulate a single star fixed in space, thereby suppressing potentia…
▽ More
The pressure exerted by massive stars' radiation fields is an important mechanism regulating their formation. Detailed simulation of massive star formation therefore requires an accurate treatment of radiation. However, all published simulations have either used a diffusion approximation of limited validity; have only been able to simulate a single star fixed in space, thereby suppressing potentially-important instabilities; or did not provide adequate resolution at locations where instabilities may develop. To remedy this we have developed a new, highly accurate radiation algorithm that properly treats the absorption of the direct radiation field from stars and the re-emission and processing by interstellar dust. We use our new tool to perform three-dimensional radiation-hydrodynamic simulations of the collapse of massive pre-stellar cores with laminar and turbulent initial conditions and properly resolve regions where we expect instabilities to grow. We find that mass is channeled to the stellar system via gravitational and Rayleigh-Taylor (RT) instabilities, in agreement with previous results using stars capable of moving, but in disagreement with methods where the star is held fixed or with simulations that do not adequately resolve the development of RT instabilities. For laminar initial conditions, proper treatment of the direct radiation field produces later onset of instability, but does not suppress it entirely provided the edges of radiation-dominated bubbles are adequately resolved. Instabilities arise immediately for turbulent pre-stellar cores because the initial turbulence seeds the instabilities. Our results suggest that RT features are significant and should be present around accreting massive stars throughout their formation.
△ Less
Submitted 11 July, 2016;
originally announced July 2016.
-
Hybrid Adaptive Ray-Moment Method (HARM$^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids
Authors:
Anna L. Rosen,
Mark R. Krumholz,
Jeffrey S. Oishi,
Aaron T. Lee,
Richard I. Klein
Abstract:
We present a highly-parallel multi-frequency hybrid radiation hydrodynamics algorithm that combines a spatially-adaptive long characteristics method for the radiation field from point sources with a moment method that handles the diffuse radiation field produced by a volume-filling fluid. Our Hybrid Adaptive Ray-Moment Method (HARM$^2$) operates on patch-based adaptive grids, is compatible with as…
▽ More
We present a highly-parallel multi-frequency hybrid radiation hydrodynamics algorithm that combines a spatially-adaptive long characteristics method for the radiation field from point sources with a moment method that handles the diffuse radiation field produced by a volume-filling fluid. Our Hybrid Adaptive Ray-Moment Method (HARM$^2$) operates on patch-based adaptive grids, is compatible with asynchronous time stepping, and works with any moment method. In comparison to previous long characteristics methods, we have greatly improved the parallel performance of the adaptive long-characteristics method by developing a new completely asynchronous and non-blocking communication algorithm. As a result of this improvement, our implementation achieves near-perfect scaling up to $\mathcal{O}(10^3)$ processors on distributed memory machines. We present a series of tests to demonstrate the accuracy and performance of the method.
△ Less
Submitted 6 July, 2016;
originally announced July 2016.
-
Gone With the Wind: Where is the Missing Stellar Wind Energy from Massive Star Clusters?
Authors:
Anna L. Rosen,
Laura A. Lopez,
Mark R. Krumholz,
Enrico Ramirez-Ruiz
Abstract:
Star clusters larger than $\sim 10^{3}$ $M_\odot$ contain multiple hot stars that launch fast stellar winds. The integrated kinetic energy carried by these winds is comparable to that delivered by supernova explosions, suggesting that at early times winds could be an important form of feedback on the surrounding cold material from which the star cluster formed. However, the interaction of these wi…
▽ More
Star clusters larger than $\sim 10^{3}$ $M_\odot$ contain multiple hot stars that launch fast stellar winds. The integrated kinetic energy carried by these winds is comparable to that delivered by supernova explosions, suggesting that at early times winds could be an important form of feedback on the surrounding cold material from which the star cluster formed. However, the interaction of these winds with the surrounding clumpy, turbulent, cold gas is complex and poorly understood. Here we investigate this problem via an accounting exercise: we use empirically determined properties of four well-studied massive star clusters to determine where the energy injected by stellar winds ultimately ends up. We consider a range of kinetic energy loss channels, including radiative cooling, mechanical work on the cold interstellar medium, thermal conduction, heating of dust via collisions by the hot gas, and bulk advection of thermal energy by the hot gas. We show that, for at least some of the clusters, none of these channels can account for more than a small fraction of the injected energy. We suggest that turbulent mixing at the hot-cold interface or physical leakage of the hot gas from the HII region can efficiently remove the kinetic energy injected by the massive stars in young star clusters. Even for the clusters where we are able to account for all the injected kinetic energy, we show that our accounting sets strong constraints on the importance of stellar winds as a mechanism for feedback on the cold interstellar medium.
△ Less
Submitted 6 May, 2014;
originally announced May 2014.
-
What Sets the Initial Rotation Rates of Massive Stars?
Authors:
Anna L. Rosen,
Mark R. Krumholz,
Enrico Ramirez-Ruiz
Abstract:
The physical mechanisms that set the initial rotation rates in massive stars are a crucial unknown in current star formation theory. Observations of young, massive stars provide evidence that they form in a similar fashion to their low-mass counterparts. The magnetic coupling between a star and its accretion disk may be sufficient to spin down low-mass pre-main sequence (PMS) stars to well below b…
▽ More
The physical mechanisms that set the initial rotation rates in massive stars are a crucial unknown in current star formation theory. Observations of young, massive stars provide evidence that they form in a similar fashion to their low-mass counterparts. The magnetic coupling between a star and its accretion disk may be sufficient to spin down low-mass pre-main sequence (PMS) stars to well below breakup at the end stage of their formation when the accretion rate is low. However, we show that these magnetic torques are insufficient to spin down massive PMS stars due to their short formation times and high accretion rates. We develop a model for the angular momentum evolution of stars over a wide range in mass, considering both magnetic and gravitational torques. We find that magnetic torques are unable to spin down either low or high mass stars during the main accretion phase, and that massive stars cannot be spun down significantly by magnetic torques during the end stage of their formation either. Spin-down occurs only if massive stars' disk lifetimes are substantially longer or their magnetic fields are much stronger than current observations suggest.
△ Less
Submitted 19 January, 2012;
originally announced January 2012.