Search

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., 20 pages, 21 figures, submitted to MNRAS

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., Comment: 24 pages, 20 figures, extra figures at https://github.com/guszejnovdavid/STARFORGE_IMF_paper_extra_plots

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 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., Accepted for publication in MNRAS. 17 pages, 10 figures. Comments welcome. Some data products and plotting code available at https://github.com/mikegrudic/StarforgeFullPhysics. Visit http://www.starforge.space to learn more about the STARFORGE Project. Simulation movie available at http://www.starforge.space/M2e4_fullphysics_flyaround.mp4

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 days 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 M >1 M 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. 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.

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., Comment: 18 pages, 15 figures, submitted to MNRAS

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., Comment: 25 pages, 22 figures. Updated to match published version

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., Comment: Submitted to MNRAS. 11 pages, 8 figures. See http://www.starforge.space/turbsphere_sigma_gas.mp4 for an animated version of Figure 2

This is an open access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited., Context. The infrared dark cloud (IRDC) SDC335.579-0.292 (hereafter, SDC335) is a massive (similar to 5000 M-circle dot) star-forming cloud which has been found to be globally collapsing towards one of the most massive star forming cores in the Galaxy, which is located at its centre. SDC335 is known to host three high-mass protostellar objects at early stages of their evolution and archival ALMA Cycle 0 data (at similar to 5 resolution) indicate the presence of at least one molecular outflow in the region detected in HNC. 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 during their formation.Aims. 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 young massive protostars in this cloud.Methods. We imaged the molecular outflows in SDC335 using new data from the Australia Telescope Compact Array of SiO and Class I CH3OH maser emission (at a resolution of similar to 3 '') alongside observations of four CO transitions made with the Atacama Pathfinder EXperiment and archival Atacama Large Millimeter/submillimeter Array (ALMA) CO, (CO)-C-13 (similar to 1 ''), and HNC data. We introduced a generalised argument to constrain outflow inclination angles based on observed outflow properties. We then used the properties of each outflow to infer the accretion rates on the protostellar sources driving them. These accretion properties allowed us to deduce the evolutionary characteristics of the sources. Shock-tracing SiO emission and CH3OH Class I maser emission allowed us to locate regions of interaction between the outflows and material infalling to the central region via the filamentary arms of SDC335.Results. We identify three molecular outflows in SDC335 - one associated with each of the known compact HII regions in the IRDC. These outflows have velocity ranges of similar to 10 km s(-1) and temperatures of similar to 60 K. The two most massive sources (separated by similar to 9000 AU) have outflows with axes which are, in projection, perpendicular. A well-collimated jet-like structure with a velocity gradient of similar to 155 km s(-1) pc(-1) is detected in the lobes of one of the outflows. 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-circle dot. The measured total accretion rate, inferred from the outflows, onto the protostars is 1.4(+/- 0.1) x 10(-3) M-circle dot yr(-1), which is comparable to the total mass infall rate toward the cloud centre on parsec scales of 2.5(+/- 1.0) x 10(-3) M-circle dot yr(-1), suggesting a near-continuous flow of material from cloud to core scales. Finally, we identify multiple regions where the outflows interact with the infalling material in the cloud's six filamentary arms, creating shocked regions and pumping Class I methanol maser emission. These regions provide useful case studies for future investigations of the disruptive effect of young massive stars on their natal clouds.© A. Avison et al. 2021., The Australia Telescope Compact Array is part of the Australia Telescope which is funded by the Commonwealth of Australia for operation as a National Facility managed by CSIRO. The authors would like to thank all ATNF staff past and present who helped during the ATCA observation used in this paper, particularly those who provided A.A. with curry. The authors would also like to thank the anonymous referee for their input into the paper after initial submission which has helped to improve the work. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00474.S and #2012.0.00781.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This publication is based on data acquired with the Atacama Pathfinder EXperiment (APEX). APEX is a collaboration between the Max-Planck-Institut fuer Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory. A.A. is funded by the STFC at the UK ARC Node. G.A.F acknowledges financial support from the State Agency for Research of the Spanish MCIU through the AYA2017-84390-C2-1-R grant (co-funded by FEDER) and through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofisica de Andalucia (SEV-2017-0709). N.P. wishes to acknowledge support under STFC consolidated grants ST/N000706/1 and ST/S00033X/1. A.D.C acknowledges the support from the UK STFC consolidated grant ST/N000706/1. A.L.R acknowledges support from NASA through Einstein Postdoctoral Fellowship grant number PF7-180166 awarded by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for NASA under contract NAS8-03060. This research made use of APLpy, an open-source plotting package for Python hosted at http://aplpy.github.com.This research made use of Astropy (http://www.astropy.org), a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018).

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 low densities. The simulation with outflows has an excess of diffuse gas compared to the simulations without outflows, exhibits an increased average sonic Mach number, and maintains a slower star formation rate (SFR) 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 SFRs and longer depletion times.

Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks., Comment: 57 pages, Space Science Reviews, topical collection Star Formation

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., Comment: 53 pages, 7 figures, To appear in Space Science Reviews (submitted January 31, 2020; accepted May 7, 2020), topical collection on Star formation

Low-energy electron diffraction and low-energy electron microscopy were used to study the (3 × 1) phase of Ag on Ge(111) for temperatures less than 540 °C. This phase was observed when depositing 0.05–0.1 ML Ag at 370 °C. The ( 3×3)R30o phase formed when depositing 0.3 ML Ag at 170 °C, and then annealing to 200−360 °C resulted in a nonreversible phase transition to [(4 × 4) + (3 × 1)] phases at 360 °C. The 3 phase appears to be metastable for temperatures

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 $\alpha = 0.03$ to $\alpha = 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., Comment: ApJ Accepted

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., Accepted to ApJ, 22 pages, 12 figures

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., 13 pages, 8 figures, ApJ, in press

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., 23 pages, 22 figures, 2 tables. Submitted to MNRAS. Simulation movies can be found here: https://anna-rosen.com/movies/

Intermediate and massive stars drive fast and powerful isotropic winds that interact with the winds of nearby stars in star clusters and the surrounding interstellar medium (ISM). Wind-ISM collisions generate astrospheres around these stars that contain hot $T\sim 10^7$ K gas that adiabatically expands. As individual bubbles expand and collide they become unstable, potentially driving turbulence in star clusters. In this paper we use hydrodynamic simulations to model a densely populated young star cluster within a homogeneous cloud to study stellar wind collisions with the surrounding ISM. We model a mass-segregated cluster of 20 B-type young main sequence stars with masses ranging from 3--17 $M_{\odot}$. We evolve the winds for $\sim$11 kyrs and show that wind-ISM collisions and over-lapping wind-blown bubbles around B-stars mixes the hot gas and ISM material generating Kolmogorov-like turbulence on small scales early in its evolution. We discuss how turbulence driven by stellar winds may impact the subsequent generation of star formation in the cluster, 12 pages, 5 figures, Accepted for publication in ApJL

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., Comment: 21 pages, 16 figures, accepted by ApJ

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., 27 pages, 8 figures, submitted to Journal of Computational Physics. Referee comments received and will be addressed in final version

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., 11 pages, 6 figures. Submitted to the Journal of Open Source Software. The source code is available at https://github.com/yt-project/unyt and the documentation is available at http://unyt.readthedocs.io/en/latest/. Accepted for publication in the Journal of Open Source Software

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., Comment: 17 pages, 12 figures, submitted to MNRAS

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes., 29 pages, 7 figures, Proceedings Joint Discussion 2. To be published in Hightlights of Astronomy, ed. T. Montmerle

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., 12 pages, 10 figures, Accepted for publication in ApJ