Milky Way Research Articles

Validating full-spectrum fitting with a synthetic integral-field spectroscopic observation of the milky way

Abstract Ongoing deep IFS observations of disk galaxies provide opportunities for comparison with the Milky Way (MW) to understand galaxy evolution. However, such comparisons are marred by many challenges such as selection effects, differences in observations and methodology, and proper validation of full-spectrum fitting methods. In this study, we present a novel code GalCraft to address these challenges by generating mock IFS data cubes of the MW using simple stellar population models and a mock MW stellar catalog derived from E-Galaxia. We use the widely adopted full-spectrum fitting code pPXF to investigate the ability to recover kinematics and stellar populations for an edge-on mock MW IFS observation. We confirm that differences in kinematics, mean age, [M/H], and [α/Fe] between thin and thick disks can be distinguished. However, the age distribution is overestimated in the ranges between 2 − 4 and 12 − 14 Gyr compared to the expected values. This is likely due to the age spacing and degeneracy of SSP templates. We find systematic offsets in the recovered kinematics due to insufficient spectral resolution and the variation of line-of-sight velocity distribution with age and [M/H]. With future higher resolution and multi-[α/Fe] SSP templates, GalCraft will be useful to validate key signatures such as [α/Fe]-[M/H] distribution at different R and |z| and potentially infer radial migration and kinematic heating efficiency to study detailed chemodynamical evolution of MW-like galaxies.

Milky Way and nearby galaxy science with the SALTUS space observatory

Milky Way and nearby galaxy science with the SALTUS space observatory

The Interstellar Medium in Dwarf Irregular Galaxies

Dwarf irregular (dIrr) galaxies are among the most common type of galaxy in the Universe. They typically have gas-rich, low-surface-brightness, metal-poor, and relatively thick disks. Here, we summarize the current state of our knowledge of the interstellar medium (ISM), including atomic, molecular, and ionized gas, along with their dust properties and metals. We also discuss star-formation feedback, gas accretion, and mergers with other dwarfs that connect the ISM to the circumgalactic and intergalactic media. We highlight one of the most persistent mysteries: the nature of pervasive gas that is yet undetected as either molecular or cold hydrogen, the “dark gas.” Some highlights include the following: ▪Significant quantities of Hi are in far-outer gas disks.▪Cold Hi in dIrrs would be molecular in the Milky Way, making the chemical properties of star-forming clouds significantly different.▪Stellar feedback has a much larger impact in dIrrs than in spiral galaxies.▪The escape fraction of ionizing photons is significant, making dIrrs a plausible source for reionization in the early Universe.▪Observations suggest a significantly higher abundance of hydrogen (H2 or cold Hi) associated with CO in star-forming regions than that traced by the CO alone.

Hot gas accretion fuels star formation faster than cold accretion in high-redshift galaxies

ABSTRACT We use high-resolution ($\simeq$35pc) hydrodynamical simulations of galaxy formation to investigate the relation between gas accretion and star formation in galaxies hosted by dark matter haloes of mass $10^{12}$${{\mathrm{M}}_{\odot }}$ at $z = 2$. At high-redshift, cold-accreted gas is expected to be readily available for star formation, while gas accreted in a hot mode is expected to require a longer time to cool down before being able to form stars. Contrary to these expectations, we find that the majority of cold-accreted gas takes several hundred Myr longer to form stars than hot-accreted gas after it reaches the inner circumgalactic medium (CGM). Approximately 10 per cent of the cold-accreted gas flows rapidly through the inner CGM on to the galactic disc. The remaining 90 per cent is trapped in a turbulent accretion region that extends up to $\sim 50$ per cent of the virial radius, from which it takes several hundred Myr for the gas to be transported to the star-forming disc. In contrast, most hot shock-heated gas avoids this ‘slow track’, and accretes directly from the CGM on to the disc where stars can form. We find that shock-heating of cold gas after accretion in the inner CGM and supernova-driven outflows contribute to, but do not fully explain, the delay in star formation. These processes combined slow down the delivery of cold-accreted gas to the galactic disc and consequently limit the rate of star formation in Milky Way mass galaxies at $z \gt 2$.

Frank Drake is Alive! (Rethinking the Drake Equation for the Search For Biological Life)

In 1961, astronomer Frank Drake formulated the Drake Equation as a cornerstone for scientific discourse regarding the prevalence of communicative extraterrestrial civilizations within the Milky Way galaxy. This equation, often referred to as the “Classic Drake Equation”, outlines the key factors influencing the number of potential civilizations with which we might establish communication. This article submerges into the Drake Equation and proposes a simplified version focused on the broader detection of extraterrestrial non-intelligent life. The established terms of the equation, such as the rate of stellar formation, the fraction of stars harboring planetary systems, and the probability of such systems containing habitable planets, are re-examined and discussed. Additionally, a reevaluation of other factors is presented. Based on this revised framework, various scenarios are explored. As our technological capabilities continue to advance, the detection of biosignatures on exoplanets (incorporated into the suggested new version of the equation) is anticipated to offer insights into the search for life beyond Earth. Keywords: Drake Equation, Extraterrestrial Life, Exoplanetology, Habitability, Biomarkers, Communicative Civilizations

Regulating star formation in a magnetized disc galaxy

ABSTRACT We use high-resolution magnetohydrodynamic simulations of isolated disc galaxies to investigate the co-evolution of magnetic fields with a self-regulated, star-forming interstellar medium (ISM). The simulations are conducted using the ramses adaptive mesh refinement code on the standard agora initial condition, with gas cooling, star formation, and feedback. We run galaxies with a variety of initial magnetic field strengths. The fields evolve and achieve approximate saturation within 500 Myr, but at different levels. The galaxies reach a quasi-steady state, with slowly declining star formation due to both gas consumption and increase in the field strength at intermediate ISM densities. We connect this behaviour to differences in the gas properties and overall structure of the galaxies. Stronger magnetic fields limit supernova bubble sizes. Different cases support the ISM using varying combinations of magnetic pressure, turbulence, and thermal energy. Initially, $\gtrsim\!\! 1\ \mu \mathrm{ G}$ magnetic fields evolve modestly and dominate support at all radii. Conversely, initially weaker fields grow through feedback and turbulence but never dominate the support. This is reflected in the stability of the gas disc. This interplay determines the overall distribution of star formation in each case. We conclude that an initially weak field can grow to produce a realistic model of a local disc galaxy, but starting with typically assumed field strengths ($\gtrsim\!\! 1\ \mu \mathrm{ G}$) will not.

The JWST EXCELS survey: too much, too young, too fast? Ultra-massive quiescent galaxies at 3 < z < 5

ABSTRACT We report ultra-deep, medium-resolution spectroscopic observations for four quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot })\gt 11$ at $3 \lt z \lt 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first two galaxies are newly selected from PRIMER UDS imaging, both at $z=4.62$ and separated by 860 pkpc on the sky, within a larger structure for which we confirm several other members. Both formed at $z\simeq 8-10$. These systems could plausibly merge by the present day to produce a local massive elliptical galaxy. The other two ultra-massive quiescent galaxies are previously known at $z=3.99$ and 3.19, with the latter (ZF-UDS-7329) having been the subject of debate as potentially too old and too massive to be accommodated by the $\Lambda$-CDM halo-mass function. Both exhibit high stellar metallicities, and for ZF-UDS-7329 we are able to measure the $\alpha -$enhancement, obtaining [Mg/Fe] = $0.42^{+0.19}_{-0.17}$. We finally evaluate whether these four galaxies are consistent with the $\Lambda$-CDM halo-mass function using an extreme value statistics approach. We find that the $z=4.62$ objects and the $z=3.19$ object are unlikely within our area under the assumption of standard stellar fractions ($f_*\simeq 0.1-0.2$). However, these objects roughly align with the most massive galaxies expected under the assumption of 100 per cent conversion of baryons to stars ($f_*$=1). Our results suggest extreme galaxy formation physics during the first billion years, but no conflict with $\Lambda$-CDM cosmology.

A global view on star formation: The GLOSTAR Galactic plane survey

Context. The GLObal view on STAR formation in the Milky Way (GLOSTAR) survey studies star formation with the Very Large Array (VLA) and the Effelsberg 100 meter radio telescope in the Galactic plane between −2° < ℓ < 60° and |b| < 1°, and the Cygnus X region (76° < ℓ < 83° and −1° < b < 2°), with unprecedented sensitivity in both flux density (1σ ~50 μJy beam−1) and the capability of detecting emission with angular scales in the range from 1.″0 to the largest radio structures in the Galaxy on the order of a few degrees in size. Aims. Here, we provide a complete GLOSTAR-VLA D-configuration radio source catalog for the part of the Galactic disk covered. A catalog for the “pilot region” between 28° < ℓ < 36° has been published in a previous paper and here we present the complementary catalog for the area within 2° < ℓ < 28°, 36° < ℓ < 60° and |b| < 1°. Methods. Observations were taken with the VLA in a 4–8 GHz band to image 100 square degrees of the inner Galactic disk at a reference frequency of 5.8 GHz, using a total of 260 h of telescope time. We determined spectral indices (α; Sν ∝ να) inside the observed band and in the frequency range of 1.4–5.8 GHz by complementing our results with those from The HI/OH/Recombination line survey of the inner Milky Way (THOR), which covers 1–2 GHz. Results. The final images have an angular resolution of 18″ and an average sensitivity of 123 μJy beam−1. The sensitivity is better (~60 μJy beam−1) in areas free of extended emission. The complementary Galactic disk catalog presented in this work consists of 11 211 radio sources. Of these, 1965 are known large-scale structure sources such as star-forming region complexes, well-known supernova remnants (SNRs), SNR candidates, or parts thereof. The remaining 9227 are discrete individual sources. Source parameters–namely flux densities, sizes, spectral indices, and classifications-are reported. We identify 769 H II region candidates, 359 of which have been newly classified as such. The mean value of spectral indices of 225 H II regions is +0.14 ± 0.02, consistent with most of them emitting optically thin thermal radio emission. Combining our results with the previously published catalog of the pilot region, the final GLOSTAR-VLA D-configuration catalog contains 12 981 radio sources.

Reconciling M/L Ratios Across Cosmic Time: a Concordance IMF for Massive Galaxies

The stellar initial mass function (IMF) is thought to be bottom heavy in the cores of the most massive galaxies, with an excess of low-mass stars compared to the Milky Way. However, studies of the kinematics of quiescent galaxies at 2 < z < 5 find M/L ratios that indicate lighter IMFs. Light IMFs have also been proposed for the unexpected populations of luminous galaxies that JWST has uncovered at z > 7 to reduce tensions with galaxy formation models. Here we explore “ski slope” IMFs that are simultaneously bottom heavy, with a steep slope at low stellar masses, and top heavy, with a shallow slope at high masses. We derive a form of the IMF for massive galaxies that is consistent with measurements in the local universe and yet produces relatively low M/L ratios at high redshift. This concordance IMF has slopes γ 1 = 2.40 ± 0.09, γ 2 = 2.00 ± 0.14, and γ 3 = 1.85 ± 0.11 in the regimes 0.08 M ⊙ – 0.5 M ⊙, 0.5 M⊙ – 1 M ⊙, and >1 M ⊙, respectively. The IMF parameter α, the mass excess compared to a Milky Way IMF, ranges from log(α)≈+0.3 for present-day galaxies to log(α)≈−0.1 for their star-forming progenitors. The concordance IMF applies only to the central regions of the most massive galaxies, with velocity dispersions σ ∼ 300 km s−1, and their progenitors. However, it can be generalized using a previously measured relation between α and σ. We arrive at the following modification to the Kroupa (2001) IMF for galaxies with σ ≳ 160 km s−1: γ1≈1.3+4.3logσ160 , γ2≈2.3−1.2logσ160 , and γ3≈2.3−1.7logσ160 , with σ 160 = σ/160 km s−1. If galaxies grow primarily inside out, so that velocity dispersions are relatively stable, these relations should also hold at high redshift.

Close Encounters of Wide Binaries Induced by the Galactic Tide: Implications for Stellar Mergers and Gravitational-wave Sources

A substantial fraction of stars can be found in wide binaries with projected separations between ∼102 and 105 au. In the standard lore of binary physics, these would evolve as effectively single stars that remotely orbit one another on stationary Keplerian ellipses. However, embedded in their Galactic environment, the low binding energy of wide binaries makes them exceptionally prone to perturbations from the gravitational potential of the Milky Way and encounters with passing stars. Employing a fully relativistic N-body integration scheme, we study the impact of these perturbations on the orbital evolution of wide binaries along their trajectory through the Milky Way. Our analysis reveals that the torques exerted by the Galaxy can cause large-amplitude oscillations of the binary eccentricity to 1 − e ≲ 10−8. As a consequence, the wide binary members pass close to each other at periapsis, which, depending on the type of binary, potentially leads to a mass transfer or collision of stars or to an inspiral and subsequent merger of compact remnants due to gravitational-wave radiation. Based on a simulation of 105 wide binaries across the Galactic field, we find that this mechanism could significantly contribute to the rate of stellar collisions and binary black hole mergers as inferred from observations of luminous red novae and gravitational-wave events by LIGO/Virgo/Kagra. We conclude that the dynamics of wide binaries, despite their large mean separation, can give rise to extreme interactions between stars and compact remnants.

The Formation of Milky Way “Bones”: Ubiquitous HI Narrow Self-absorption Associated with CO Emission

Long and skinny molecular filaments running along Galactic spiral arms are known as “bones,” since they make up the skeleton of the Milky Way. However, their origin is still an open question. Here, we compare spectral images of HI taken by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) with archival CO and Herschel dust emission to investigate the conversion from HI to H2 in two typical Galactic bones, CFG028.68-0.28 and CFG047.06+0.26. Sensitive FAST HI images and an improved methodology enabled us to extract HI narrow self-absorption (HINSA) features associated with CO line emission on and off the filaments, revealing the ubiquity of HINSA toward distant clouds for the first time. The derived cold HI abundances, [HI]/[H2], of the two bones range from ∼(0.5 to 44.7) × 10−3, which reveal different degrees of HI–H2 conversion, and are similar to those of nearby, low-mass star-forming clouds, Planck Galactic cold clumps, and a nearby active high-mass star-forming region G176.51+00.20. The HI–H2 conversion has been ongoing for 2.2–13.2 Myr in the bones, a timescale comparable to that of massive star formation therein. Therefore, we are witnessing young giant molecular clouds (GMCs) with rapid massive star formation. Our study paves the way of using HINSA to study cloud formation in Galactic bones and, more generally, in distant GMCs in the FAST era.

Galactic Chemical Evolution Models Favor an Extended Type Ia Supernova Delay-time Distribution

Type Ia supernovae (SNe Ia) produce most of the Fe-peak elements in the Universe and therefore are a crucial ingredient in galactic chemical evolution models. SNe Ia do not explode immediately after star formation, and the delay-time distribution (DTD) has not been definitively determined by supernova surveys or theoretical models. Because the DTD also affects the relationship among age, [Fe/H], and [α/Fe] in chemical evolution models, comparison with observations of stars in the Milky Way is an important consistency check for any proposed DTD. We implement several popular forms of the DTD in combination with multiple star formation histories for the Milky Way in multizone chemical evolution models that include radial stellar migration. We compare our predicted interstellar medium abundance tracks, stellar abundance distributions, and stellar age distributions to the final data release of the Apache Point Observatory Galactic Evolution Experiment. We find that the DTD has the largest effect on the [α/Fe] distribution: a DTD with more prompt SNe Ia produces a stellar abundance distribution that is skewed toward a lower [α/Fe] ratio. While the DTD alone cannot explain the observed bimodality in the [α/Fe] distribution, in combination with an appropriate star formation history it affects the goodness of fit between the predicted and observed high-α sequence. Our model results favor an extended DTD with fewer prompt SNe Ia than the fiducial t −1 power law.

The Most Sensitive Radio Recombination Line Measurements Ever Made of the Galactic Warm Ionized Medium

Diffuse ionized gas pervades the disk of the Milky Way. We detect extremely faint emission from this Galactic warm ionized medium (WIM) using the Green Bank Telescope to make radio recombination line (RRL) observations toward two Milky Way sight lines: G20, (ℓ, b) = (20°, 0°), and G45, (ℓ, b) = (45°, 0°). We stack 18 consecutive Hnα transitions between 4.3 and 7.1 GHz to derive 〈Hnα〉 spectra that are sensitive to RRL emission from plasmas with emission measures EM ≳ 10 cm−6 pc. Each sight line has two Gaussian-shaped spectral components with emission measures that range between ∼100 and ∼300 cm−6 pc. Because there is no detectable RRL emission at negative LSR velocities, the emitting plasma must be located interior to the solar orbit. The G20 and G45 emission measures imply rms densities of 0.15 and 0.18 cm−3, respectively, if these sight lines are filled with homogeneous plasma. The observed 〈Hnβ〉/〈Hnα〉 line ratios are consistent with LTE excitation for the strongest components. The high-velocity component of G20 has a narrow line width, 13.5 km s−1, that sets an upper limit of ≲4000 K for the plasma electron temperature. This is inconsistent with the ansatz of a canonically pervasive, low-density, ∼10,000 K WIM plasma.

On the Kinematic Nature of Apparent Disks at High Redshifts: Local Counterparts are Not Dominated by Ordered Rotation but by Tangentially Anisotropic Random Motion

It is not straightforward to physically interpret the apparent morphology of galaxies. Recent observations by the James Webb Space Telescope (JWST) revealed a dominant galaxy population at high redshifts (z > 2) that were visually classified as disks for their flattened shapes and/or exponential light profiles. The extensively accepted interpretation is that they are dynamically cold disks supported by bulk rotation. However, it is long known that flattened shapes and exponential profiles are not exclusive for rotating disk structure. To break degeneracy and assess the rotational support of typical high-z galaxies in the JWST samples, those with active star formation and stellar masses lg(M⋆/M⊙)∼9 , we study the kinematics of their equal-mass counterparts at z = 0. While these local star-forming low-mass galaxies are photometrically similar to real dynamically cold disks, they are not supported by ordered rotation but primarily by random motion, and their flattened shapes result largely from tangential orbital anisotropy. Given the empirical and theoretical evidence that young galaxies are dynamically hotter at higher redshifts, our results suggest that the high-z JWST galaxies may not be cold disks but are dynamically warm/hot galaxies with flattened shapes driven by anisotropy. While both have low rotational support, local low-mass galaxies possess oblate shapes, contrasting the prolate shapes (i.e., cigar like) of low-mass systems at high redshifts. Such shape transition (prolate ⇒ oblate) indicates an associated change in orbital anisotropy (radial ⇒ tangential), with roots likely in the assembly of their host dark matter halos.

Abundances of Planetary Nebulae and Evolved Stars: Iron and Sulfur Depletion, and Carbon and Nitrogen Enrichment, in Low- and Intermediate-mass Stellar Populations in the Milky Way

We research the elemental abundances in Galactic planetary nebulae (PNe) compared with those of their stellar progenitors (red giant branch and asymptotic giant branch, AGB, stars), to explore and quantify the expected—i.e., due to AGB evolution or condensation onto grains—differences. We gleaned the current literature for the nebular abundances while we used the APOGEE DR 17 survey data for the stellar sample. We examined the elements in common between the nebular and stellar samples, namely, C, N, O, Fe, and S. We confirm that iron in PNe is mostly entrapped in grains, with an average depletion 〈D[Fe/H]〉 = 1.741 ± 0.486 dex, and we disclose a weak correlation between iron depletion and the [O/H] abundance, D[Fe/H] = (6.6003 ± 2.443) × [O/H] + (1.972 ± 0.199). Sulfur may also be mildly depleted in PNe, with 〈D[S/H]〉 = 0.179 ± 0.291 dex. We also found an indication of nitrogen enrichment for PNe 〈E[N/H]〉 = 0.393 ± 0.421 dex, with maximum enrichment (0.980 ± 0.243) occurring for the PNe whose progenitors have gone through the hot-bottom burning. The carbon enrichment is 〈E[C/H]〉 = 0.337 ± 0.463 dex when measured for the general PN populations. Our results will be relevant for future Galactic and extragalactic studies comparing nebular and stellar samples.

A global view on star formation: The GLOSTAR Galactic plane survey

Studies of Galactic H II regions are of crucial importance for studying star formation and the evolution of the interstellar medium. Gaining an insight into their physical characteristics contributes to a more comprehensive understanding of these phenomena. The GLOSTAR project aims to provide a GLObal view on STAR formation in the Milky Way by performing an unbiased and sensitive survey. This is achieved by using the extremely wideband (4–8 GHz) C-band receiver of the Karl G. Jansky Very Large Array and the Effelsberg 100 m telescope. Using radio recombination lines observed in the GLOSTAR survey with the VLA in D-configuration with a typical line sensitivity of 1 σ ~ 3.0 mJy beam−1 at ~5 km s−1 and an angular resolution of 25″, we cataloged 244 individual Galactic H II regions (−2° ≤ ℓ ≤ 60° and |b| ≤ 1°, and 76° ≤ ℓ ≤ 83° and −1° ≤ b ≤ 2°) and derived their physical properties. We examined the mid-infrared (MIR) morphology of these H II regions and find that a significant portion of them exhibit a bubble-like morphology in the GLIMPSE 8 μm emission. We also searched for associations with the dust continuum and sources of methanol maser emission, other tracers of young stellar objects, and find that 48% and 14% of our H II regions, respectively, are coextensive with those. We measured the electron temperature for a large sample of H II regions within Galactocentric distances spanning from 1.6 to 13.1 kpc and derived the Galactic electron temperature gradient as ~372 ± 28 K kpc−1 with an intercept of 4248 ± 161 K, which is consistent with previous studies.

Landscapes under the stars

Realizing he’d never seen the Milky Way with his own eyes, Colin White joined photographer Shaun Davey on a nighttime expedition to capture the heart of the galaxy.

Study of X-ray emission from the S147 nebula by SRG/eROSITA: Supernova-in-the-cavity scenario

The Simeis 147 nebula (S147) is particularly well known for a spectacular net of Hα-emitting filaments. It is often considered one of the largest and oldest (∼105 yr) cataloged supernova remnants in the Milky Way, although the kinematics of the pulsar PSR J0538+2817 suggests that this supernova remnant might be a factor of three younger. The former case is considered in a companion paper, while here we pursue the latter. Both studies are based on the data of SRG/eROSITA All-Sky Survey observations. Here, we confront the inferred properties of the X-ray emitting gas data with the scenario of a supernova explosion in a low-density cavity, such as a wind-blown-bubble. This scenario assumes that a ∼20 M⊙ progenitor star has had a low velocity with respect to the ambient interstellar medium, and so stayed close to the center of a dense shell created during its main-sequence evolution till the moment of the core-collapse explosion. The ejecta first propagate through the low-density cavity until they collide with the dense shell, and only then does the reverse shock go deeper into the ejecta and power the observed X-ray emission of the nebula. The part of the remnant inside the dense shell remains non-radiative till this point, plausibly in a state with Te &lt; Ti and nonequilibrium ionization. On the contrary, the forward shock becomes radiative immediately after entering the dense shell, and, being subject to instabilities, gives the nebula its characteristic “foamy” appearance in Hα and radio emission.

Gamma rays from dark matter spikes in EAGLE simulations

Intermediate Mass Black Holes (IMBHs) with a mass range between 100 M⊙ and 106 M⊙ are expected to be surrounded by high dark matter densities, so-called dark matter spikes.The high density of self-annihilating Weakly Interacting Massive Particles (WIMPs) in these spikes leads to copious gamma-ray production. Sufficiently nearby IMBHs could therefore appear as unidentified gamma-ray sources.However, the number of IMBHs and their distribution within our own Milky Way is currently unknown. In this work, we provide a mock catalogue of IMBHs and their dark matter spikes obtained from the EAGLE simulations, in which black holes with a mass of 105 M⊙/h are seeded into the centre of halos greater than 1010 M⊙/h to model black hole feedback influencing the formation of galaxies.The catalogue contains the coordinates and dark matter spike parameters for about 2500 IMBHs present in about 150 Milky Way-like galaxies. We expect about 15+9 -6 IMBHs within our own galaxy, mainly distributed in the Galactic Centre and the Galactic Plane. In the most optimistic scenario, we find that current and future gamma-ray observatories, such as Fermi-LAT, H.E.S.S. and CTAO, would be sensitive enough to probe the cross section of dark matter self-annihilation around IMBHs down to many orders of magnitude below the thermal relic cross section for dark matter particles with masses from GeV to TeV.We have made the IMBH mock catalogue and the source code for our analysis publicly available, providing the resources to study dark matter self-annihilation around IMBHs with current and upcoming gamma-ray observatories.

Suppressed Cosmic-Ray Energy Densities in Molecular Clouds from Streaming Instability-regulated Transport

Cosmic rays (CRs) are the primary driver of ionization in star-forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic-ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a 2000M ⊙ MC and the subsequent star formation including CRT using the STAR FORmation in Gaseous Environments framework implemented in the GIZMO code. We show that when CRT is streaming-dominated, the CR energy in the cloud is strongly attenuated due to energy losses from the streaming instability. Consequently, in a Milky Way–like environment the median CR ionization rate in the cloud is low (ζ ≲ 2 × 10−19 s−1) during the main star-forming epoch of the calculation and the impact of CRs on the star formation in the cloud is limited. However, in high-CR environments, the CR distribution in the cloud is elevated (ζ ≲ 6 × 10−18), and the relatively higher CR pressure outside the cloud causes slightly earlier cloud collapse and increases the star formation efficiency by 50% to ∼13%. The initial mass function is similar in all cases except with possible variations in a high-CR environment. Further studies are needed to explain the range of ionization rates observed in MCs and explore star formation in extreme CR environments.