Search for: All records

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 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. 

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. 

Mass segregation is seen in many star clusters, but whether massive stars form in the center of a cluster or migrate there dynamically is still debated.N-body simulations show that early dynamical mass segregation is possible when sub-clusters merge to form a dense core with a small crossing time. However, the effect of gas dynamics on both the formation and dynamics of the stars could inhibit the formation of the dense core. We aim to study the dynamical mass segregation of star cluster models that include gas dynamics and selfconsistently form stars from the dense substructure in the gas. Our models use the TORCH framework, which is based on AMUSE and includes stellar and magnetized gas dynamics, as well as stellar evolution and feedback from radiation, stellar winds, and supernovae. Our models consist of three star clusters forming from initial turbulent spherical clouds of mass 10⁴, 10⁵, 10⁶M_⊙and radius 11.7 pc that have final stellar masses of 3.6 × 10³M_⊙, 6.5 × 10⁴M_⊙, and 8.9 × 10⁵M_⊙, respectively. There is no primordial mass segregation in the model by construction. All three clusters become dynamically mass segregated at early times via collapse confirming that this mechanism occurs within sub-clusters forming directly out of the dense substructure in the gas. The dynamics of the embedded gas and stellar feedback do not inhibit the collapse of the cluster. We find that each model cluster becomes mass segregated within 2 Myr of the onset of star formation, reaching the levels observed in young clusters in the Milky Way. However, we note that the exact values are highly time-variable during these early phases of evolution. Massive stars that segregate to the center during core collapse are likely to be dynamically ejected, a process that can decrease the overall level of mass segregation again. 

Two main mechanisms have classically been proposed for the formation of runaway stars. In the binary supernova scenario (BSS), a massive star in a binary explodes as a supernova, ejecting its companion. In the dynamical ejection scenario, a star is ejected during a strong dynamical encounter between multiple stars. We propose a third mechanism for the formation of runaway stars: the subcluster ejection scenario (SCES), where a subset of stars from an infalling subcluster is ejected out of the cluster via a tidal interaction with the contracting gravitational potential of the assembling cluster. We demonstrate the SCES in a star-by-star simulation of the formation of a young massive cluster from a 10⁶M_⊙gas cloud using theTORCHframework. This star cluster forms hierarchically through a sequence of subcluster mergers determined by the initial turbulent, spherical conditions of the gas. We find that these mergers drive the formation of runaway stars in our model. Late-forming subclusters fall into the central potential, where they are tidally disrupted, forming tidal tails of runaway stars that are distributed highly anisotropically. Runaways formed in the same SCES have similar ages, velocities, and ejection directions. Surveying observations, we identify several SCES candidate groups with anisotropic ejection directions. The SCES is capable of producing runaway binaries: two wide dynamical binaries in infalling subclusters were tightened through ejection. This allows for another velocity kick via subsequent via a subsequent BSS ejection. An SCES-BSS ejection is a possible avenue for the creation of hypervelocity stars unbound to the Galaxy. The SCES occurs when subcluster formation is resolved. We expect nonspherical initial gas distributions to increase the number of calculated runaway stars, bringing it closer to observed values. The observation of groups of runaway stars formed via the SCES can thus reveal the assembly history of their natal clusters. 

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. 

We report the first high-resolution, detailed abundances of 21 elements for giants in the Galactic bulge/bar within 1° of the Galactic plane, where high extinction has rendered such studies challenging. Our high-signal-to-noise-ratio and high-resolution, near-infrared spectra of seven M giants in the inner bulge, located at (l,b) = (0°, +1°), are observed using the IGRINS spectrograph. We report the first multichemical study of the inner Galactic bulge by investigating, relative to a robust new solar neighborhood sample, the abundance trends of 21 elements, including the relatively difficult to study heavy elements. The elements studied are: F, Mg, Si, S, Ca, Na, Al, K, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, Ce, Nd, and Yb. We investigate bulge membership of all seven stars using distances and orbital simulations, and we find that the most metal-poor star may be a halo interloper. Our investigation also shows that the inner bulge as close as 1° north of the Galactic Center displays a similarity to the inner disk sequence, following the high-[α/Fe] envelope of the solar vicinity metal-rich population, though no firm conclusions for a different enrichment history are evident from this sample. We find a small fraction of metal-poor stars ([Fe/H] > −0.5), but most of our stars are mainly of supersolar metallicity. Fluorine is found to be enhanced at high metallicity compared to the solar neighborhood, but confirmation with a larger sample is required. We will apply this approach to explore the populations of the nuclear stellar disk and the nuclear star cluster. 

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. 

Context.The star cluster initial mass function is observed to have an inverse power law exponent around 2, yet there is no consensus on what determines this distribution, and why some variation is observed in different galaxies. Furthermore, the cluster formation efficiency (CFE) covers a range of values, particularly when considering different environments. These clusters are often used to empirically constrain star formation and as fundamental units for stellar feedback models. Detailed galaxy models must therefore accurately capture the basic properties of observed clusters to be considered predictive. Aims.We study how feedback mechanisms acting on different timescales and with different energy budgets affect the star cluster mass function and CFE. Methods.We use hydrodynamical simulations of a dwarf galaxy as a laboratory to study star cluster formation. We test different combinations of stellar feedback mechanisms, including stellar winds, ionizing radiation, and supernovae (SNe). Results.Each feedback mechanism affects the CFE and cluster mass function. Increasing the feedback budget by combining the different types of feedback decreases the CFE by reducing the number of massive clusters. Ionizing radiation is found to be especially influential. This effect depends on the timing of feedback initiation, as shown by comparing early and late feedback. Early feedback occurs from ionizing radiation and stellar winds with onset immediately after a massive star is formed. Late feedback occurs when energy injection only starts after the main-sequence lifetime of the most massive SN progenitor, a timing that is further influenced by the choice of the most massive SN progenitor. Late feedback alone results in a broad, flat mass function, approaching a log-normal shape in the complete absence of feedback. Early feedback, on the other hand, produces a power-law cluster mass function with lower CFE, albeit with a steeper slope than that usually observed. 

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 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 doesn’t 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 midplane 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. We suggest that grand design spiral arms might produce isolated clumps of CNM at these radii.