Search for: All records

ABSTRACT The relationship between magnetic field strength B and gas density n in the interstellar medium is of fundamental importance. We present and compare Bayesian analyses of the B–n relation for two comprehensive observational data sets: a Zeeman data set and 700 observations using the Davis–Chandrasekhar–Fermi (DCF) method. Using a hierarchical Bayesian analysis we present a general, multiscale broken power-law relation, $$B=B_0(n/n_0)^{\alpha }$$, with $$\alpha =\alpha _1$$ for $$n< n_0$$ and $$\alpha _2$$ for $$n>n_0$$, and with $$B_0$$ the field strength at $$n_0$$. For the Zeeman data, we find: $$\alpha _1={0.15^{+0.06}_{-0.09}}$$ for diffuse gas and $$\alpha _2 = {0.53^{+0.09}_{-0.07}}$$ for dense gas with $$n_0 = 0.40^{+1.30}_{-0.30}\times 10^4$$ cm$$^{-3}$$. For the DCF data, we find: $$\alpha _1={0.26^{+0.01}_{-0.01}}$$ and $$\alpha _2={0.77_{-0.15}^{+0.14}}$$, with $$n_0=14.00^{+10.00}_{-7.00}\times 10^4$$ cm$$^{-3}$$, where the uncertainties give 68 per cent credible intervals. We perform a similar analysis on nineteen numerical magnetohydrodynamic simulations covering a wide range of physical conditions from protostellar discs to dwarf and Milky Way-like galaxies, computed with the arepo, flash, pencil, and ramses codes. The resulting exponents depend on several physical factors such as dynamo effects and their time-scales, turbulence, and initial seed field strength. We find that the dwarf and Milky Way-like galaxy simulations produce results closest to the observations. 

Abstract The Aeosproject introduces a series of high-resolution cosmological simulations that model star-by-star chemical enrichment and galaxy formation in the early Universe, achieving 1 pc resolution. These simulations capture the complexities of galaxy evolution within the first ~300 Myr by modeling individual stars and their feedback processes. By incorporating chemical yields from individual stars, Aeosgenerates galaxies with diverse stellar chemical abundances, linking them to hierarchical galaxy formation and early nucleosynthetic events. These simulations underscore the importance of chemical abundance patterns in ancient stars as vital probes of early nucleosynthesis, star formation histories, and galaxy formation. We examine the metallicity floors of various elements resulting from Population III enrichment, providing best-fit values for eight different metals (e.g., [O/H] = −4.0) to guide simulations without Population III models. Additionally, we identify galaxies that begin star formation with Population II after external enrichment and investigate the frequency of carbon-enhanced metal-poor stars at varying metallicities. The Aeossimulations offer detailed insights into the relationship between star formation, feedback, and chemical enrichment. Future work will extend these simulations to later epochs to interpret the diverse stellar populations of the Milky Way and its satellites. 

Abstract The first generations of stars left their chemical fingerprints on metal-poor stars in the Milky Way and its surrounding dwarf galaxies. While instantaneous and homogeneous enrichment implies that groups of conatal stars should have the same element abundances, small amplitudes of abundance scatter are seen at fixed [Fe/H]. Measurements of intrinsic abundance scatter have been made with small high-resolution spectroscopic data sets where measurement uncertainty is small compared to this scatter. In this work, we present a method to use mid-resolution survey data, which have larger errors, to make this measurement. Using APOGEE Data Release 17, we calculate the intrinsic scatter of Al, O, Mg, Si, Ti, Ni, and Mn relative to Fe for 333 metal-poor stars across six classical dwarf galaxies around the Milky Way, and 1604 stars across 19 globular clusters (GCs). We calibrate the reported abundance errors in bins of signal-to-noise ratio and [Fe/H] using a high-fidelity halo data set. Applying these calibrated errors to the APOGEE data, we find small amplitudes of average intrinsic abundance scatter in dwarf galaxies ranging from 0.03 to 0.09 dex, with a median value of 0.047 dex. For the GCs, we find intrinsic scatters ranging from 0.01 to 0.11 dex, with particularly high scatter for Al and O. Our measurements of intrinsic abundance scatter place important upper bounds, which are limited by our calibration, on the intrinsic scatter in these systems, as well as constraints on their underlying star formation history and mixing that we can look to simulations to interpret. 

Abstract Magnetic fields are now widely recognized as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modeled. ISM models exist of turbulent magnetic fields using small-scale dynamo (SSD). Others model the large-scale dynamo (LSD) organizing magnetic fields at the scale of the disk or spiral arms. Separately, neither can fully describe the galactic magnetic field dynamics nor topology. We model the LSD and SSD together at a sufficient resolution to use the low explicit Lagrangian resistivity required. The galactic SSD saturates within 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr versus 1–2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due toα-quenching near the midplane as the growing mean-field produces a magneticαthat opposes the kineticα. The magnetic energy in our models of the LSD shows a slightly sublinear response to increasing resolution, indicating that we are converging toward the physical solution at 1 pc resolution. Clustering supernovae in OB associations increases the growth rates for both the SSD and the LSD, compared to a horizontally uniform supernova distribution. 

ABSTRACT Protoplanetary discs spend their lives in the dense environment of a star-forming region. While there, they can be affected by nearby stars through external photoevaporation and dynamic truncations. We present simulations that use the amuse framework to couple the torch model for star cluster formation from a molecular cloud with a model for the evolution of protoplanetary discs under these two environmental processes. We compare simulations with and without extinction of photoevaporation-driving radiation. We find that the majority of discs in our simulations are considerably shielded from photoevaporation-driving radiation for at least 0.5 Myr after the formation of the first massive stars. Radiation shielding increases disc lifetimes by an order of magnitude and can let a disc retain more solid material for planet formation. The reduction in external photoevaporation leaves discs larger and more easily dynamically truncated, although external photoevaporation remains the dominant mass-loss process. Finally, we find that the correlation between disc mass and projected distance to the most massive nearby star (often interpreted as a sign of external photoevaporation) can be erased by the presence of less massive stars that dominate their local radiation field. Overall, we find that the presence and dynamics of gas in embedded clusters with massive stars is important for the evolution of protoplanetary discs. 

ABSTRACT The formation and evolution of galaxies have proved sensitive to the inclusion of stellar feedback, which is therefore crucial to any successful galaxy model. We present INFERNO, a new model for hydrodynamic simulations of galaxies, which incorporates resolved stellar objects with star-by-star calculations of when and where the injection of enriched material, momentum, and energy takes place. INFERNO treats early stellar kinematics to include phenomena such as walkaway and runaway stars. We employ this innovative model on simulations of a dwarf galaxy and demonstrate that our physically motivated stellar feedback model can drive vigorous galactic winds. This is quantified by mass and metal loading factors in the range of 10–100, and an energy loading factor close to unity. Outflows are established close to the disc, are highly multiphase, spanning almost 8 orders of magnitude in temperature, and with a clear dichotomy between mass ejected in cold, slow-moving (T ≲ 5 × 104 K, v < 100 km s−1) gas and energy ejected in hot, fast-moving (T > 106 K, v > 100 km s−1) gas. In contrast to massive disc galaxies, we find a surprisingly weak impact of the early stellar kinematics, with runaway stars having little to no effect on our results, despite exploding in diffuse gas outside the dense star-forming gas, as well as outside the galactic disc entirely. We demonstrate that this weak impact in dwarf galaxies stems from a combination of strong feedback and a porous interstellar medium, which obscure any unique signatures that runaway stars provide. 

ABSTRACT The cold neutral medium (CNM) is an important part of the galactic gas cycle and a precondition for the formation of molecular and star-forming gas, yet its distribution is still not fully understood. In this work, we present extremely high resolution simulations of spiral galaxies with time-dependent chemistry such that we can track the formation of the CNM, its distribution within the galaxy, and its correlation with star formation. We find no strong radial dependence between the CNM fraction and total neutral atomic hydrogen (H i) due to the decreasing interstellar radiation field counterbalancing the decreasing gas column density at larger galactic radii. However, the CNM fraction does increase in spiral arms where the CNM distribution is clumpy, rather than continuous, overlapping more closely with H2. The CNM does not extend out radially as far as H i, and the vertical scale height is smaller in the outer galaxy compared to H i with no flaring. The CNM column density scales with total mid-plane pressure and disappears from the gas phase below values of PT/kB = 1000 K cm−3. We find that the star formation rate density follows a similar scaling law with CNM column density to the total gas Kennicutt–Schmidt law. In the outer galaxy, we produce realistic vertical velocity dispersions in the H i purely from galactic dynamics, but our models do not predict CNM at the extremely large radii observed in H i absorption studies of the Milky Way. We suggest that extended spiral arms might produce isolated clumps of CNM at these radii. 

Abstract Magnetic fields grow quickly, even at early cosmological times, suggesting the action of a small-scale dynamo (SSD) in the interstellar medium (ISM) of galaxies. Many studies have focused on idealized, isotropic, homogeneous, turbulent driving of the SSD. Here we analyze more realistic simulations of supernova-driven turbulence to understand how it drives an SSD. We find that SSD growth rates are intermittently variable as a result of the evolving multiphase ISM structure. Rapid growth in the magnetic field typically occurs in hot gas, with the highest overall growth rates occurring when the fractional volume of hot gas is large. SSD growth rates correlate most strongly with vorticity and fluid Reynolds number, which also both correlate strongly with gas temperature. Rotational energy exceeds irrotational energy in all phases, but particularly in the hot phase while SSD growth is most rapid. Supernova rate does not significantly affect the ISM average kinetic energy density. Rather, higher temperatures associated with high supernova rates tend to increase SSD growth rates. SSD saturates with total magnetic energy density around 5% of equipartition to kinetic energy density, increasing slightly with magnetic Prandtl number. While magnetic energy density in the hot gas can exceed that of the other phases when SSD grows most rapidly, it saturates below 5% of equipartition with kinetic energy in the hot gas, while in the cold gas it attains 100%. Fast, intermittent growth of the magnetic field appears to be a characteristic behavior of supernova-driven, multiphase turbulence. 

Abstract Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suiteTorch, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer,N-body dynamics, and stellar feedback, to model four initially identical 10⁴M_⊙giant molecular clouds with a Gaussian density profile peaking at 521.5 cm⁻³. Using theTorchsoftware suite through theAMUSEframework, we modify three of the models, to ensure that the first star that forms is very massive (50, 70, and 100M_⊙). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star-forming region. The star formation rate is suppressed, reducing the total mass of the stars formed. Our fiducial control model, without an early massive star, has a larger star formation rate and total efficiency by up to a factor of 3, and a higher average star formation efficiency per freefall time by up to a factor of 7. Early-forming massive stars promote the buildup of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster. 

ABSTRACT We perform simulations of star cluster formation to investigate the morphological evolution of embedded star clusters in the earliest stages of their evolution. We conduct our simulations with Torch, which uses the Amuse framework to couple state-of-the-art stellar dynamics to star formation, radiation, stellar winds, and hydrodynamics in Flash. We simulate a suite of 104 M⊙ clouds at 0.0683 pc resolution for ∼2 Myr after the onset of star formation, with virial parameters α = 0.8, 2.0, 4.0 and different random samplings of the stellar initial mass function and prescriptions for primordial binaries. Our simulations result in a population of embedded clusters with realistic morphologies (sizes, densities, and ellipticities) that reproduce the known trend of clouds with higher initial α having lower star formation efficiencies. Our key results are as follows: (1) Cluster mass growth is not monotonic, and clusters can lose up to half of their mass while they are embedded. (2) Cluster morphology is not correlated with cluster mass and changes over ∼0.01 Myr time-scales. (3) The morphology of an embedded cluster is not indicative of its long-term evolution but only of its recent history: radius and ellipticity increase sharply when a cluster accretes stars. (4) The dynamical evolution of very young embedded clusters with masses ≲1000 M⊙ is dominated by the overall gravitational potential of the star-forming region rather than by internal dynamical processes such as two- or few-body relaxation.