Star Formation Research Articles

A Study of the Abundance of Aluminium in Massive Stars

Abstract: This study presents the results of the calculation of the mass fractions of isotopes from hydrogen to nickel during the entire core burning phases of a 25 M⊙ star. A simple stellar model was used, considering four main parameters: temperature, density, initial composition, and time, without accounting for mass loss or mixing. The mass fractions of the isotopes were calculated using the open-source package NucNet Tools and the updated reaction rates from the JINA Reaclib database. A comparison of our results with existing similar data is performed and acceptable agreement is conspicuous. Hydrodynamic conditions in massive stars favor the production of 26Al, therefore particular attention was given to the isotopes 26Al and 27Al and the ratio of their mass fractions R_Al for comparison with the literature. The averaged value of the controversial R_Al during the star's lifetime is found to be 1.68×10-4. Keywords: Elemental abundances in stars, Evolution, Stellar, Nuclear astrophysics, Hydrostatic stellar nucleosynthesis, Nucleosynthesis: stellar, Stars formation. PACS: 97.10.Tk, 87.23.-n, 97.10.Cv, 26.20.-f, 7.10.Cv, 97.10.Bt.

The cosmic rate of pair-instability supernovae

Abstract Pair-instability supernovae (PISNe) have crucial implications for many astrophysical topics, including the search for very massive stars, the black hole mass spectrum, and galaxy chemical enrichment. To this end, we need to understand where PISNe are across cosmic time, and what are their favourable galactic environments. We present a new determination of the PISN rate as a function of redshift, obtained by combining up-to-date stellar evolution tracks from the PARSEC and FRANEC codes, with an up-to-date semi-empirical determination of the star formation rate and metallicity evolution of star-forming galaxies throughout cosmic history. We find the PISN rate to exhibit a huge dependence on the model assumptions, including the criterion to identify stars unstable to pair production, and the upper limit of the stellar initial mass function. Remarkably, the interplay between the maximum metallicity at which stars explode as PISNe, and the dispersion of the galaxy metallicity distribution, dominates the uncertainties, causing a ∼ seven-orders-of-magnitude PISN rate range. Furthermore, we show a comparison with the core-collapse supernova rate, and study the properties of the favourable PISN host galaxies. According to our results, the main contribution to the PISN rate comes from metallicities between ∼10−3 and 10−2, against the common assumption that views very-low-metallicity, Population III stars as exclusive or dominant PISN progenitors. The strong dependencies we find offer the opportunity to constrain stellar and galaxy evolution models based on possible future (or the lack of) PISN observations.

Extreme-value modelling of the brightest galaxies at z ≳ 9

Abstract Data from the James Webb Space Telescope have revealed an intriguing population of bright galaxies at high redshifts. In this work, we use extreme-value statistics to calculate the distribution (in UV magnitude) of the brightest galaxies in the redshift range 9 ≲ z ≲ 16. We combine the Generalised Extreme Value (GEV) approach with modelling of the galaxy luminosity function. We obtain predictions of the brightest galaxies for a suite of luminosity functions, including the Schechter and double power law functions, as well as a model parametrised by the stellar formation efficiency f*. We find that the JWST data is broadly consistent with f* of $5~{{\ \rm per\ cent}}-10~{{\ \rm per\ cent}}$, and that the brightest galaxy at z ∼ 16 will have $M_{\rm UV}\approx -23.5^{0.8}_{0.4}$. If f* is dependent on halo mass, we predict $M_{\rm UV}\approx -22.5^{0.5}_{1.5}$ for such an object. We show that extreme-value statistics not only predicts the magnitude of the brightest galaxies at high redshifts, but may also be able to distinguish between models of star formation in high-redshift galaxies.

Gas-phase Fe/O and Fe/N abundances in star-forming regions. Relations between nucleosynthesis, metallicity, and dust

Gas-phase Fe/O and Fe/N abundances in star-forming regions. Relations between nucleosynthesis, metallicity, and dust

Merging Gas-rich Galaxies That Harbor Low-luminosity Twin Quasars at z = 6.05: A Promising Progenitor of the Most Luminous Quasars

We present Atacama Large Millimeter/submillimeter Array [C ii] 158 μm line and underlying far-IR continuum emission observations (0.″57 × 0.″46 resolution) toward a quasar–quasar pair system recently discovered at z = 6.05. The quasar nuclei (C1 and C2) are faint (M 1450 ≳ −23 mag), but we detect very bright [C ii] emission bridging the 12 kpc between the two objects and extending beyond them (total luminosity L [C ii] ≃ 6 × 109 L ⊙). The [C ii]-based total star formation rate of the system is ∼550 M ⊙ yr−1 (the IR-based dust-obscured star formation is ∼100 M ⊙ yr−1), with a [C ii]-based total gas mass of ∼1011 M ⊙. The dynamical masses of the two galaxies are large (∼9 × 1010 M ⊙ for C1 and ∼5 × 1010 M ⊙ for C2). There is a smooth velocity gradient in [C ii], indicating that these quasars are a tidally interacting system. We identified a dynamically distinct, fast-[C ii] component around C1: detailed inspection of the line spectrum there reveals the presence of a broad-wing component, which we interpret as the indication of fast outflows with a velocity of ∼600 km s−1. The expected mass-loading factor of the outflows, after accounting for multiphase gas, is ≳2 − 3, which is intermediate between AGN-driven and starburst-driven outflows. Hydrodynamic simulations in the literature predict that this pair will evolve to a luminous (M 1450 ≲ −26 mag), starbursting (≳1000 M ⊙ yr−1) quasar after coalescence, one of the most extreme populations in the early Universe.

FORGE’d in FIRE III: The IMF in Quasar Accretion Disks from STARFORGE

Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10 pc from the SMBH. We show it is critical to include physics often previously neglected, including magnetic fields, radiation, and (proto)stellar feedback. Closer to the SMBH, star formation is suppressed, but the (rare) stars that do form exhibit top-heavy IMFs. Stars can form only in special locations (e.g. magnetic field switches) in the outer QAD. Protostars accrete their natal cores rapidly but then dynamically decouple from the gas and ‘wander,’ ceasing accretion on timescales ~100 yr. Their jets control initial core accretion, but the ejecta are ‘swept up’ into the larger-scale QAD flow without much dynamical effect. The strong tidal environment strongly suppresses common-core multiplicity. The IMF shape depends sensitively on un-resolved dynamics of protostellar disks (PSDs), as the global dynamical times can become incredibly short (< yr) and tidal fields are incredibly strong, so whether PSDs can efficiently transport angular momentum or fragment catastrophically at <10 au scales requires novel PSD simulations to properly address. Most analytic IMF models and analogies with planet formation in PSDs fail qualitatively to explain the simulation IMFs, though we discuss a couple of viable models.

Detection of Two Totally Eclipsing B-type Binaries with Extremely Low Mass Ratios

The detection of O- and B-type stars with extremely low-mass companions is very important for understanding the formation and evolution of binary stars. However, finding them remains a challenge because the low-mass components in such systems contribute such small flux to the total. During our search for pulsations among O- and B-type stars using the Transiting Exoplanet Survey Satellite (TESS) data, we found two short-period and B-type (B9) eclipsing binaries with orbital periods of 1.61613 and 2.37857 days. Photometric solutions of the two close binaries were derived by analyzing the TESS light curves with the Wilson–Devinney method. It is discovered that both of them are detached binaries with extremely low mass ratios of 0.067(2) for TIC 260342097 and 0.140(3) for TIC 209148631. The determined mass ratio indicates that TIC 260342097 is one of the lowest mass ratios among known B-type binary systems. We showed that the two systems have total eclipses with a broad and flat secondary minimum, suggesting that the photometric parameters could be derived reliably. The absolute parameters of the two binaries are estimated and it is found that the secondary components in the two systems are overluminous and oversize when compared with the normal low-mass and cool main-sequence (MS) stars. These findings may imply that the two systems are composed of a B-type MS primary and a cool pre-MS secondary with orbital periods shorter than 2.5 days. They are valuable targets to test theories of binary star formation and evolution.

Constraints on Primordial Magnetic Fields from High Redshift Stellar Mass Density

Primordial magnetic fields (PMFs) play a pivotal role in influencing small-scale fluctuations within the primordial density field, thereby enhancing the matter power spectrum within the context of the ΛCDM model at small scales. These amplified fluctuations accelerate the early formation of galactic halos and stars, which can be observed through advanced high-redshift observational techniques. Therefore, stellar mass density (SMD) observations, which provide significant opportunities for detailed studies of galaxies at small scales and high redshifts, offer a novel perspective on small-scale cosmic phenomena and constrain the characteristics of PMFs. In this study, we compile 14 SMD data points at redshifts z > 6 and derive stringent constraints on the parameters of PMFs, which include the amplitude of the magnetic field at a characteristic scale of λ = 1 Mpc, denoted as B 0, and the spectral index of the magnetic field power spectrum, n B . At 95% confidence level, we establish upper limits of B 0 < 4.44 nG and n B < −2.24, along with a star formation efficiency of approximately f*0∼0.1 . If we fix n B at specific values, such as −2.85, −2.9, and −2.95, the 95% upper limits for the amplitude of the magnetic field can be constrained to 1.33, 2.21, and 3.90 nG, respectively. Finally, we attempt to interpret recent early observations provided by the James Webb Space Telescope using the theory of PMFs and find that by selecting appropriate PMF parameters, it is possible to explain these results without significantly increasing the star formation efficiency.

The Lyα Nondetection by JWST NIRSpec of a Strong Lyα Emitter at z = 5.66 Confirmed by MUSE

The detections of Lyα emission in galaxies with redshifts above 5 are of utmost importance for constraining the cosmic reionization timeline, yet such detections are usually based on slit spectroscopy. Here we investigate the significant bias induced by slit placement on the estimate of Lyα escape fraction ( fescLyα ) by presenting a galaxy (dubbed A2744-z6Lya) at z = 5.66, where its deep JWST/NIRSpec prism spectroscopy completely misses the strong Lyα emission detected in the MUSE data. A2744-z6Lya exhibits a pronounced UV continuum with an extremely steep spectral slope of β=−2.640−0.008+0.008 , and it has a stellar mass of ∼108.75 M ⊙, a star formation rate of ∼5.95 M ⊙ yr−1, and a gas-phase metallicity of 12+log(O/H)∼7.90 . The observed flux and rest-frame equivalent width of its Lyα from MUSE spectroscopy are 1.2 × 10−16 erg s−1 cm−2 and 75 Å, equivalent to fescLyα=42%±1% . However, its Lyα nondetection from JWST NIRSpec gives a 5σ upper limit of <5%, in stark contrast to that derived from MUSE. To explore the reasons for this bias, we perform spatially resolved stellar population analysis of A2744-z6Lya using the JWST NIRCam and HST imaging data to construct 2D maps of star formation rate, dust extinction, and neutral hydrogen column density. We find that the absence of Lyα in the slit regions probably stems from both the resonance scattering effect of neutral hydrogen and dust extinction. Through analyzing an extreme case in detail, this work highlights the important caveat of inferring fescLyα from slit spectroscopy, particularly when using the JWST multiplexed NIRSpec microshutter assembly.

Nuclear acoustic envelope soliton in a relativistically degenerate magneto-rotating stellar plasma

The study explores a model called Relativistic Degenerate Magneto-Rotating Quantum Plasma (RDMRQP), which comprises a static heavy nucleus, an inertial non-degenerate light nucleus, and warm non-relativistic or ultra-relativistic electrons. The focus is on observing the emergence of Nucleus-Acoustic Envelope Solitons (NAESs). Using the reductive perturbation method, a Nonlinear Schrödinger equation (NLSE) is derived to characterize the properties of NAESs. We have analysed the Peregrine breather soliton solution. The investigation reveals that the temperature of warm degenerate species, plasma system’s rotational speed, and the presence of heavy nucleus species can alter the fundamental features (height and width) of NAESs in the WDMRQP system. The study emphasizes the existence of only positive NA wave potential. Additionally, a phase plane analysis is conducted to gain a deeper understanding of the parametric dependencies. Through detailed mathematical and numerical analysis, the study avoids overloading with complex mathematics while demonstrating parametric dependence via phase portrait analysis. The research augments the envelop soliton model with breather mode solutions and discusses modulation instability using the Benjamin-Feir Index, highlighting the significance of solitons in star formation. Envelop solitons, stable waves within stars and proto-stars, influence energy transport and stellar evolution, playing a crucial role in the accretion process and formation of stable structures. Key findings include the effects of various parameters on NAESs’ generation and propagation in which non-relativistic and ultra-relativistic electrons support NAESs, with amplitude and width influenced by temperature, rotational frequency, inclination angle, and the presence of a static heavy nucleus. The research is applicable to hot white dwarfs and neutron stars, suggesting further exploration of quantum effects and non-planar or arbitrary amplitude NAESs.

MagMaR III—Resisting the Pressure, Is the Magnetic Field Overwhelmed in NGC6334I?

We report on Atacama Large Millimeter/submillimeter Array observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from northwest to southeast with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS (J = 5 → 4) emission yields an outflow energy of approximately 3.5 × 1045 erg, aligning with previous interferometric studies. Utilizing the Davis–Chandrasekhar–Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 ± 1 mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary HII region is also larger than the magnetic energy, suggesting that protostellar feedback may be the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions.

From Halos to Galaxies. IX. Estimate of Halo Assembly History for SDSS Galaxy Groups

The properties of the galaxies are tightly connected to their host halo mass and halo assembly history. Accurate measurement of the halo assembly history in observation is challenging but crucial to the understanding of galaxy formation and evolution. The stellar-to-halo mass ratio (M */M h) for the centrals has often been used to indicate the halo assembly time t h,50 of the group, where t h,50 is the lookback time at which a halo has assembled half of its present-day virial mass. Using mock data from the semi-analytic models, we find that M */M h shows a significant scatter with t h,50, with a strong systematic difference between the group with a star-forming central (blue group) and passive central (red group). To improve the accuracy, we develop machine learning models to estimate t h,50 for galaxy groups using only observable quantities in the mocks. Since star formation quenching will decouple the co-growth of the dark matter and baryon, we train our models separately for blue and red groups. Our models have successfully recovered t h,50, within an accuracy of ∼1.09 Gyr. With careful calibrations of individual observable quantities in the mocks with Sloan Digital Sky Survey (SDSS) observations, we apply the trained models to the SDSS Yang et al. groups and derive the t h,50 for each group for the first time. The derived SDSS t h,50 distributions are in good agreement with that in the mocks, in particular for blue groups. The derived halo assembly history, together with the halo mass, make an important step forward in studying the halo–galaxy connections in observation.

Magnetohydrodynamical modeling of star-disk formation: from isolated spherical collapse towards incorporation of external dynamics

The formation of protostars and their disks has been understood as the result of the gravitational collapse phase of an accumulation of dense gas that determines the mass reservoir of the star-disk system. Against this background, the broadly applied scenario of considering the formation of disks has been to model the collapse of a dense core assuming spherical symmetry. Our understanding of the formation of star-disk systems is currently undergoing a reformation though. The picture evolves from interpreting disks as the sole outcome of the collapse of an isolated prestellar core to a more dynamic picture where disks are affected by the molecular cloud environment in which they form. In this review, we provide a status report of the state-of-the-art of spherical collapse models that are highly advanced in terms of the incorporated physics together with constraints from models that account for the possibility of infall onto star-disk systems in simplified test setups, as well as in multi-scale simulations that cover a dynamical range from the Giant Molecular Cloud environment down to the disk. Considering the observational constraints that favor a more dynamical picture of star formation, we finally discuss the challenges and prospects in linking the efforts of tackle the problem of star-disk formation in combined multi-scale, multi-physics simulations.

Giant molecular clouds and their type classification in M 74: Toward understanding star formation and cloud evolution

Abstract We investigated the giant molecular clouds (GMCs) in M 74 (NGC 628), using data obtained from the PHANGS (Physics at High Angular resolution in Nearby GalaxieS) project. We applied GMC types according to the activity of star formation: Type I without star formation, Type II with H$\alpha$ luminosity ($L_{\mathrm{H\alpha }}$) less than $10^{37.5}\ \rm{erg\ s ^{-1}}$, and Type III with $L_{\mathrm{H\alpha }}$ greater than $10^{37.5}\ \rm{erg\ s^{-1}}$. A total of 432 GMCs were identified, with 59, 201, and 172 GMCs, for Types I, II, and III, respectively. The size and mass of the GMCs range from 23 to 238 pc and $10^{4.9}$ to $10^{7.1}\, M_{\odot }$, indicating that the mass and radius increase from Types I to III. Clusters younger than 4 Myr and H ii regions are concentrated within 150 pc of a GMC, indicating a tight association between these young objects and GMCs. The virial ratio decreases from Type I to Type III, indicating that Type III GMCs are the most gravitationally relaxed among the three. We interpret that the GMCs evolve from Type I to Type III, as previously observed in the Large Magellanic Cloud. Based on a steady-state assumption, the estimated evolutionary timescales of Types I, II, and III are 1, 5, and 4 Myr, respectively. We assume that the timescale of Type III is equal to the age of the associated clusters, indicating a GMC lifetime of 10 Myr or longer. Although Chevance et al. (2020, MNRAS, 493, 2872) investigated GMCs using the same PHANGS dataset of M 74, they did not define a GMC, reaching an evolutionary picture with a 20 Myr duration of the non-star-forming phase, which is five times longer than 4 Myr. We compare the present results with those of Chevance et al. (2020, MNRAS, 493, 2872) and argue that defining individual GMCs is essential for understanding GMC evolution.

Probing the filamentary nature of star formation in the California giant molecular cloud

Context. Recent studies suggest that filamentary structures are representative of the initial conditions of star formation in molecular clouds and support a filament paradigm for star formation, potentially accounting for the origin of the stellar initial mass function (IMF). The detailed, local physical properties of molecular filaments remain poorly characterized, however. Aims. Using Herschel imaging observations of the California giant molecular cloud, we aim to further investigate the filament paradigm for low- to intermediate-mass star formation and to better understand the exact role of filaments in the origin of stellar masses. Methods. Using the multiscale, multiwavelength extraction method getsf, we identify starless cores, protostars, and filaments in the Herschel data set and separate these components from the background cloud contribution to determine accurate core and filament properties. Results. We find that filamentary structures contribute approximately 20% of the overall mass of the California cloud, while compact sources such as dense cores contribute a mere 2% of the total mass. Considering only dense gas (defined as gas with AV,bg &gt; 4.5–7), filaments and cores contribute ~66–73% and 10–14% of the dense gas mass, respectively. The transverse half-power diameter measured for California molecular filaments has a median undeconvolved value of 0.18 pc, consistent within a factor of 2 with the typical ~0.1 pc width of nearby filaments from the Herschel Gould Belt survey. A vast majority of identified prestellar cores (~82–90%) are located within ~0.1 pc of the spines of supercritical filamentary structures. Both the prestellar core mass function (CMF) and the distribution of filament masses per unit length or filament line mass function (FLMF) are consistent with power-law distributions at the high-mass end, ΔN/ΔlogM ∝ M−1.4±0.2 at M &gt; 1 M⊙ for the CMF and ΔN/Δlog Mline ∝ Mline−1.5±0.2 for the FLMF at Mline &gt; 10 M⊙ pc−1, which are both consistent with the Salpeter power-law IMF. Based on these results, we propose a revised model for the origin of the CMF in filaments, whereby the global prestellar CMF in a molecular cloud arises from the integration of the CMFs generated by individual thermally supercritical filaments within the cloud. Conclusions. Our findings support the existence a tight connection between the FLMF and the CMF/IMF and suggests that filamentary structures represent a critical evolutionary step in establishing a Salpeter-like mass function.

JWST MIRI and NIRCam observations of NGC 891 and its circumgalactic medium

We present new JWST observations of the nearby, prototypical edge-on, spiral galaxy NGC 891. The northern half of the disk was observed with NIRCam in its F150W and F277W filters. Absorption is clearly visible in the mid-plane of the F150W image, along with vertical dusty plumes that closely resemble the ones seen in the optical. A $ kpc ^2$ area of the lower circumgalactic medium (CGM) was mapped with MIRI F770W at 12 pc scales. Thanks to the sensitivity and resolution of JWST, we detect dust emission out to $ 4$ kpc from the disk, in the form of filaments, arcs, and super-bubbles. Some of these filaments can be traced back to regions with recent star formation activity, suggesting that feedback-driven galactic winds play an important role in regulating baryonic cycling. The presence of dust at these altitudes raises questions about the transport mechanisms at play and suggests that small dust grains are able to survive for several tens of million years after having been ejected by galactic winds in the disk-halo interface. We lay out several scenarios that could explain this emission: dust grains may be shielded in the outer layers of cool dense clouds expelled from the galaxy disk, and/or the emission comes from the mixing layers around these cool clumps where material from the hot gas is able to cool down and mix with these cool cloudlets. This first set of data and upcoming spectroscopy will be very helpful to understand the survival of dust grains in energetic environments, and their contribution to recycling baryonic material in the mid-plane of galaxies.

Multivariate Predictors of Lyman Continuum Escape. I. A Survival Analysis of the Low-redshift Lyman Continuum Survey* * Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with programs GO-15626, GO-13744,

To understand how galaxies reionized the Universe, we must determine how the escape fraction of Lyman continuum (LyC) photons (f esc) depends on galaxy properties. Using the z ∼ 0.3 Low-redshift Lyman Continuum Survey (LzLCS), we develop and analyze new multivariate predictors of f esc. These predictions use the Cox proportional hazards model, a survival analysis technique that incorporates both detections and upper limits. Our best model predicts the LzLCS f esc detections with an rms scatter of 0.31 dex, better than single-variable correlations. According to ranking techniques, the most important predictors of f esc are the equivalent width (EW) of Lyman-series absorption lines and the UV dust attenuation, which track line-of-sight absorption due to H i and dust. The H i absorption EW is uniquely crucial for predicting f esc for the strongest LyC emitters, which show properties similar to weaker LyC emitters and whose high f esc may therefore result from favorable orientation. In the absence of H i information, star formation rate surface density (ΣSFR) and [O iii]/[O ii] ratio are the most predictive variables and highlight the connection between feedback and f esc. We generate a model suitable for z > 6, which uses only the UV slope, ΣSFR, and [O iii]/[O ii]. We find that ΣSFR is more important in predicting f esc at higher stellar masses, whereas [O iii]/[O ii] plays a greater role at lower masses. We also analyze predictions for other parameters, such as the ionizing-to-nonionizing flux ratio and Lyα escape fraction. These multivariate models represent a promising tool for predicting f esc at high redshift.

Metallicity Dependence of Pressure-regulated Feedback-modulated Star Formation in the TIGRESS-NCR Simulation Suite

We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over Z′≡Z/Z⊙∼0.1-3 , gas surface density Σgas ∼ 5–150 M ⊙ pc−2, and stellar surface density Σstar ∼ 1–50 M ⊙ pc−2. The range of emergent star formation rate surface density is ΣSFR ∼ 10−4–0.5 M ⊙ kpc−2 yr−1, and ISM total midplane pressure is P tot/k B = 103–106 cm−3 K, with P tot equal to the ISM weight W . For given Σgas and Σstar, we find ΣSFR∝Z′0.3 . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as ϒth∝W−0.46Z′−0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒturb+mag∝W−0.22Z′−0.18 . The reduction in ΣSFR at low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved.

Evidence for a Shallow Evolution in the Volume Densities of Massive Galaxies at z = 4–8 from CEERS

We analyze the evolution of massive (log10[M ⋆/M ⊙] > 10) galaxies at z ∼ 4–8 selected from JWST Cosmic Evolution Early Release Survey (CEERS). We infer the physical properties of all galaxies in the CEERS NIRCam imaging through spectral energy distribution (SED) fitting with dense basis to select a sample of high-redshift massive galaxies. Where available we include constraints from additional CEERS observing modes, including 18 sources with MIRI photometric coverage, and 28 sources with spectroscopic confirmations from NIRSpec or NIRCam WFSS. We sample the recovered posteriors in stellar mass from SED fitting to infer the volume densities of massive galaxies across cosmic time, taking into consideration the potential for sample contamination by active galactic nuclei. We find that the evolving abundance of massive galaxies tracks expectations based on a constant baryon conversion efficiency in dark matter halos for z ∼ 4–8. At higher redshifts, we observe an excess abundance of massive galaxies relative to this simple model, resulting in a shallower decline of observed volume densities of massive galaxies. These higher abundances can be explained by modest changes to star formation physics and/or the efficiencies with which star formation occurs in massive dark matter halos, and are not in tension with modern cosmology.

Formation of Galactic Disks. II. The Physical Drivers of Disk Spin-up

Using a representative sample of Milky Way (MW)–like galaxies from the TNG50 cosmological simulation, we investigate physical processes driving the formation of galactic disks. A disk forms as a result of the interplay between inflow and outflow carrying angular momentum in and out of the galaxy. Interestingly, the inflow and outflow have remarkably similar distributions of angular momentum, suggesting an exchange of angular momentum and/or outflow recycling, leading to continuous feeding of prealigned material from the corotating circumgalactic medium. We show that the disk formation in TNG50 is correlated with stellar bulge formation, in qualitative agreement with a recent theoretical model of disk formation facilitated by steep gravitational potentials. Disk formation is also correlated with the formation of a hot circumgalactic halo with around half of the inflow occurring at subsonic and transonic velocities corresponding to Mach numbers of ≲2. In the context of recent theoretical works connecting disk settling and hot halo formation, our results imply that the subsonic part of the inflow may settle into a disk while the remaining supersonic inflow will perturb this disk via the chaotic cold accretion. We find that disks tend to form when the host halos become more massive than ∼(1–2) × 1011 M ⊙, consistent with previous theoretical findings and observational estimates of the predisk protogalaxy remnant in the MW. Our results do not prove that either corotating outflow recycling, gravitational potential steepening, or hot halo formation cause disk formation, but they show that all these processes occur concurrently and may play an important role in disk growth.