Source author record

Anna L. Rosen

Anna L. Rosen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.HE astro-ph.IM astro-ph.EP physics.comp-ph

Catalog footprint

What is connected

16works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

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.

preprint2022arXiv

A Massive Star is Born: How Feedback from Stellar Winds, Radiation Pressure, and Collimated Outflows Limits Accretion onto Massive Stars

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.

preprint2022arXiv

Cluster assembly and the origin of mass segregation in the STARFORGE simulations

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.

preprint2022arXiv

Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations

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.

preprint2022arXiv

The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert

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.

preprint2021arXiv

Less wrong: a more realistic initial condition for simulations of turbulent molecular clouds

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.

preprint2021arXiv

The Effects of Magnetic Fields and Outflow Feedback on the Shape and Evolution of the Density PDF in Turbulent Star-Forming Clouds

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.

preprint2020arXiv

Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity

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.

preprint2020arXiv

Evolution of Stellar Feedback in HII Regions

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.

preprint2020arXiv

The Role of Outflows, Radiation Pressure, and Magnetic Fields in Massive Star Formation

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.

preprint2020arXiv

Zooming in on Individual Star Formation: Low- and High-mass Stars

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.

preprint2019arXiv

Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores

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.

preprint2016arXiv

An Unstable Truth: How Massive Stars get their Mass

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.

preprint2016arXiv

Hybrid Adaptive Ray-Moment Method (HARM$^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids

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.

preprint2014arXiv

Gone With the Wind: Where is the Missing Stellar Wind Energy from Massive Star Clusters?

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.

preprint2012arXiv

What Sets the Initial Rotation Rates of Massive Stars?

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.

Anna L. Rosen

What is connected

Connect this record

See the researcher in context

Building this map preview

16 published item(s)

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

A Massive Star is Born: How Feedback from Stellar Winds, Radiation Pressure, and Collimated Outflows Limits Accretion onto Massive Stars

Cluster assembly and the origin of mass segregation in the STARFORGE simulations

Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations

The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert

Less wrong: a more realistic initial condition for simulations of turbulent molecular clouds

The Effects of Magnetic Fields and Outflow Feedback on the Shape and Evolution of the Density PDF in Turbulent Star-Forming Clouds

Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity

Evolution of Stellar Feedback in HII Regions

The Role of Outflows, Radiation Pressure, and Magnetic Fields in Massive Star Formation

Zooming in on Individual Star Formation: Low- and High-mass Stars

Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores

An Unstable Truth: How Massive Stars get their Mass

Hybrid Adaptive Ray-Moment Method (HARM$^2$): A Highly Parallel Method for Radiation Hydrodynamics on Adaptive Grids

Gone With the Wind: Where is the Missing Stellar Wind Energy from Massive Star Clusters?

What Sets the Initial Rotation Rates of Massive Stars?