2024-10-30 | | Total: 36
The innermost portions of the Milky Way's stellar halo have avoided scrutiny until recently. The lack of wide-area survey data, made it difficult to reconstruct an uninterrupted view of the density distribution of the metal-poor stars inside the Solar radius. In this study, we utilize red giant branch (RGB) stars from Gaia, with metallicities estimated using spectro-photometry from Gaia Data Release 3. Accounting for Gaia's selection function, we examine the spatial distribution of metal-poor ([M/H]<-1.3) RGB stars, from the Galactic centre (r~1 kpc) out to beyond the Solar radius (r~18 kpc). Our best-fitting single-component cored power-law model shows a vertical flattening of ~0.5 and a slope -3.4, consistent with previous studies. Motivated by the mounting evidence for two distinct stellar populations in the inner halo, we additionally test a range of two-component models. One of the components models the tidal debris from the Gaia Sausage/Enceladus merger, while the other captures the Aurora population -- stars that predate the Galactic disk formation. Our best-fit two-component model suggests that both populations contribute equally around the Solar radius, but Aurora dominates the inner halo with a steeper power-law index of -4.5, in agreement with the nitrogen-rich star distribution measured by Horta et al. (2021).
After over three decades of unsuccessful attempts, we report the first detection of molecular gas emission in Malin 1, the largest spiral galaxy observed to date, and one of the most iconic giant low surface brightness galaxies. Using ALMA, we detect significant $^{12}$CO(J=1-0) emission in the galaxy's central region and tentatively identify CO emission across three regions on the disc. These observations allow for a better estimate of the H$_2$ mass and molecular gas mass surface density, both of which are remarkably low given the galaxy's scale. By integrating data on its HI mass, we derive a very low molecular-to-atomic gas mass ratio. Overall, our results highlight the minimal presence of molecular gas in Malin 1, contrasting sharply with its extensive, homogeneous atomic gas reservoir. For the first time, we position Malin 1 on the Kennicutt-Schmidt (K-S) diagram, where it falls below the main sequence for normal spirals, consistent with previous upper limits but now with more accurate figures. These findings are crucial for constraining our understanding of star formation processes in environments characterized by extremely low molecular gas densities and for refining models of galaxy formation, thereby improving predictions concerning the formation, evolution, and distribution of these giant, elusive galaxies.
We present the findings of a comprehensive and detailed analysis of merger tree data from ultra-high-resolution cosmological $N$-body simulations. The analysis, conducted with a particle mass resolution of $5 \times 10^3 h^{-1} M_{\odot}$ and a halo mass resolution of $10^7 h^{-1} M_{\odot}$, provides sufficient accuracy to suppress numerical artefacts. This study elucidates the dynamical evolution of subhaloes associated with the Milky Way-like host haloes. Unlike more massive dark matter haloes, which have been extensively studied, these subhaloes follow a distinct mass evolution pattern: an initial accretion phase, followed by a tidal stripping phase where mass is lost due to the tidal forces of the host halo. The transition from accretion to stripping, where subhaloes reach their maximum mass, occurs around a redshift of $z\simeq1$. Smaller subhaloes reach this point earlier, while larger ones do so later. Our analysis reveals that over 80 per cent of subhaloes have experienced mass loss, underscoring the universality of tidal stripping in subhalo evolution. Additionally, we derived the eccentricities and pericentre distances of subhalo orbits from the simulations and compare them with those of nearby satellite galaxies observed by the Gaia satellite. The results demonstrate a significant alignment between the orbital elements predicted by the cold dark matter model and the observed data, providing robust support for the model as a credible candidate for dark matter.
In the present-day universe, the most massive galaxies are ellipticals located in the cores of galaxy clusters, harboring the heaviest super-massive black holes (SMBHs). However the mechanisms that drive the early growth phase and subsequent transformation of these morphology and kinematics of galaxies remain elusive. Here we report (sub)kiloparsec scale observations of stars, gas, and dust in ADF22.A1, a bright dusty starburst galaxy at z=3.1, hosting a heavily obscured active galactic nucleus and residing in a proto-cluster core. ADF22.A1 is a giant spiral galaxy with the kinematics of a rotating disk with rotation velocity Vrot=530+/-10km/s and diameter larger than 30 kpc. The high specific stellar angular momentum of this system, j*=3400+/-600 kpc km/s, requires a mechanism to effectively spin-up ADF22.A1, indicating the importance of accretion from the cosmic web to supply both gas and angular momentum to galaxies in their early gas-rich starburst phase. In its inner region, gas flows along dust lanes in a bar connected with the bright dusty core and the estimated mass ratio of a bulge to SMBH matches the local relation, suggesting that bars are a key mechanism to shape the early co-evolution of these components. Comparison with cosmological simulations shows that ADF22.A1 will likely evolve into a massive elliptical at the present day, experiencing a significant reduction in angular momentum associated with subsequent galaxy mergers.
Early JWST observations have revealed the ubiquitous presence in the early Universe, up to z about16, of extreme baryon concentrations, namely forming globular clusters, extremely dense galaxies that may or may not be UV bright, and supermassive black holes in relatively low-mass galaxies. This paper is trying to pinpoint which physical conditions may have favored the formation of such concentrations, that appear to be very common at high redshifts while their formation being progressively more and more rare at lower redshifts. Building on local globular cluster evidence, it is argued that such conditions can consist in a combination of a 10 Myr extended feedback free time, coupled to low angular-momentum densities in deep local minima of the ISM vorticity field, where baryon concentrations are more likely to form. It is argued that the former condition would follow from more massive stars failing to explode as supernovae, and the latter one from low vorticity prevailing in the early Universe, in contrast to later times with their secular increase of the angular momentum density due to the cumulative effect of tidal interactions.
Many nucleosynthetic channels create the elements, but two-parameter models characterized by $\alpha$ and Fe nonetheless predict stellar abundances in the Galactic disk to accuracies of 0.02 to 0.05 dex for most measured elements, near the level of current abundance uncertainties. It is difficult to make individual measurements more precise than this to investigate lower-amplitude nucleosynthetic effects, but population studies of mean abundance patterns can reveal more subtle abundance differences. Here we look at the detailed abundances for 67315 stars from APOGEE DR17, but in the form of abundance residuals away from a best-fit two-parameter, data-driven nucleosynthetic model. We find that these residuals show complex structures with respect to age, guiding radius, and vertical action that are not random and are also not strongly correlated with sources of systematic error such as surface gravity, effective temperature, and radial velocity. The residual patterns, especially in Na, C+N, Ni, Mn, and Ce, trace kinematic structures in the Milky Way, such as the inner disk, thick disk, and flared outer disk. A principal component analysis suggests that most of the observed structure is low-dimensional and can be explained by a few eigenvectors. We find that some, but not all, of the effects in the low-$\alpha$ disk can be explained by dilution with fresh gas, so that abundance ratios resemble those of stars with higher metallicity. The patterns and maps we provide could be combined with accurate forward models of nucleosynthesis, star formation, and gas infall to provide a more detailed picture of star and element formation in different Milky Way components.
We present an exploration of the Milky Way's structural parameters using an all-sky sample of RC giants to map the stellar density from the inner to the outer parts of the Galactic disc. These evolved giants are considered to be standard candles due to their low intrinsic variance in their absolute luminosities, allowing us to estimate their distances with reasonable confidence. We exploit all-sky photometry from the AllWISE mid-infrared survey and the Gaia survey, along with astrometry from Gaia Data Release 3 and recent 3D extinction maps, to develop a probabilistic scheme in order to select with high confidence \rc{}-like stars. Our curated catalogue contains about 10 million sources, for which we estimate photometric distances based on the WISE $W1$ photometry. We then derive the selection function for our sample, which is the combined selection function of sources with both \gaia{} and \allwise{} photometry. Using the distances and accounting for the full selection function of our observables, we are able to fit a two-disc, multi-parameter model to constrain the scale height (\hz{}), scale-length (\rd{}), flaring, and the relative mass ratios of the two disc components. We illustrate and verify our methodology using mock catalogues of \rc{} stars. We find that the \rc{} population is best described by a flared thin disc with scale length \rd{}=$3.56\pm0.32$ kpc and scale height at the Sun of \hzsun{}=$0.17\pm0.01$ kpc, and a shorter and thicker disc with \rd{}=$2.59\pm0.11$ kpc, \hzsun{}=$0.45\pm0.11$ kpc, with no flare. The thicker disc constitutes 64\% of the \rc{} stellar mass beyond 3 kpc, while the thin disk shows evidence of being warped beyond 9 kpc from the Galactic center. The residuals between the predicted number density of RC stars from our axisymmetric model and the measured counts show possible evidence of a two-armed spiral perturbation in the disc of the Milky Way.
We present capivara, a fast and scalable multi-decomposition package designed to study astrophysical properties within distinct structural components of galaxies. Our spectro-decomposition code for analyzing integral field unit (IFU) data enables a more holistic approach, moving beyond conventional radial gradients and the bulge-plus-disk dichotomy. It facilitates comprehensive comparisons of integrated stellar ages and metallicities across various galactic structures. Our classification method naturally identifies outliers and organizes the different pixels based on their dominant spectral features. The algorithm leverages the scalability and GPU acceleration of Torch, outputting both a one-dimensional spectrum and a full data cube for each galaxy component, without relying on Voronoi binning. We demonstrate the capabilities of our approach using a sample of galaxies from the MaNGA survey, processing the resulting data cubes with the starlight spectral fitting code to derive both stellar population and ionized gas properties of the galaxy components. Our method effectively groups regions with similar spectral properties in both the continuum and emission lines. By aggregating the spectra of these regions, we enhance the signal-to-noise ratio of our analysis while significantly speeding up computations by reducing the number of spectra processed simultaneously. capivara will be freely available on GitHub.
Recent JWST observations have revealed a large population of intermediate/low-luminosity AGN at early times with peculiar properties, different from local AGN or luminous quasars. To better understand the physical conditions in the BLRs of these early AGN, we used the optical FeII (4434--4684 \AA) and the broad $\rm H \beta$ emission, and the ratio between their equivalent widths $R_{Fe}$, as a probe on a purposefully assembled sample. Specifically, we gathered a sample of 26 high redshift ($\langle z \rangle$=6.4) AGN, observed by JWST, with broad $\rm H\beta$ detection both in the high and low luminosity regimes (respectively 14 faint AGN and 12 quasars), to investigate their optical FeII emission properties. In addition, we carefully selected control samples at lower $z$. We found that the population of faint AGN ($\rm \log(L_{H \beta} / (erg \, s^{-1}))\lesssim 44$) exhibits a significantly lower FeII emission than their local counterparts ($R_{Fe}<$0.24 versus $R_{Fe}\simeq$0.85 in the control sample), while the quasars at the epoch of reionisation observed by JWST present a FeII emission profile that closely resembles that observed at $z<3$. We argue that the weakness of the FeII bump in the faint JWST AGN might be due to the reduced metallicity of their broad line region ($\lesssim 0.5~Z_{\odot}$), while luminous quasars have already reached chemical maturity ($\sim Z{_\odot}$ or higher). Lastly, we highlight an intriguing similarity between the spectral properties of the high redshift population of faint AGN with those harboured in local metal poor dwarf galaxies.
Blind spectroscopy of massive lensing galaxy clusters with MUSE has revealed large numbers of gravitationally-lensed Lyman-$ \alpha $ emitters exhibiting asymmetric profiles at $ 2.9 \leq z \leq 6.7 $, suggesting abundant outflows from low-mass star-forming galaxies in the early universe. Are these primaeval galaxies experiencing their first bursts of star formation, or established galaxies experiencing rejuvenation? With JWST rest-frame optical/NIR continuum imaging now available for many of these objects, we can search for older stellar populations. Here, we search for spectroscopic confirmation of outflows from these galaxies, finding a few high-signal-to-noise cases in which blueshifted interstellar absorption lines are detected. Next, we analyse the star formation histories with combined HST + JWST photometry. We find most them to be well characterised by very young, low metallicity stellar populations. However, despite the rest-frame optical/NIR coverage of JWST, we cannot place strict upper bounds on the mass in old stars (age $ > 100\,\text{Myr} $).
The evolution equation of two-point correlation function $\xi$ of galaxies can analytically describe the large scale structure of the galaxy distribution, and the solution depends also upon the initial condition. The primeval spectrum of the baryon acoustic oscillations (BAO) contains multi peaks that survived the Silk damping, and, as a relevant portion, two peaks of the primeval BAO spectrum fall into the range of current galaxy surveys. Incorporating this portion, we use a twin-peak initial power spectrum of the galaxies, and obtain the evolution solution from a redshift $z=8$ to $z=0$ in the Gaussian approximation. The outcome $\xi(r)$ at $z=0.6$ still exhibits the 100 Mpc periodic bumps as observed by the WiggleZ survey, a feature largely determined by the Jeans length $\lambda_J$ in the equation. In particular, due to the superposition of the twin peaks in the initial condition, $\xi(r)$ shows a shallow trough at $\sim 70 h^{-1}$Mpc and a deep trough at $\sim 140 h^{-1}$Mpc, agreeing with the observational data, much better than our previous work that used a simple one-peak initial spectrum.
The Milky Way hosts at least two modes in its present day distribution of Fe and alpha-elements. The exact cause of this bimodality is disputed, but one class of explanations involves the merger between the Milky Way and a relatively massive satellite (Gaia-Sausage-Enceladus) at z~2. However, reproducing this bimodality in simulations is not straightforward, with conflicting results on the prevalance, morphology, and mechanism behind multimodality. We present a case study of a galaxy in the Illustris TNG50 simulation which undergoes sequential phases of starburst, brief quiescence, and then rejuvenation. This scenario results in a pronounced abundance bimodality after a post-processing adjustment of the [alpha/Fe] of old stars designed to mimic a higher star formation efficiency in dense gas. The high- and low-alpha sequences are separated in time by the brief quiescent period, which is not associated with a merger but by the formation of a bar followed by AGN activity. This galaxy indicates a novel scenario in which the alpha-bimodality in the Milky Way is caused by the formation of the bar via AGN-induced quenching. In addition to a stellar age gap in the Milky Way, we predict that abundance bimodalities should be more common in barred as opposed to unbarred galaxies.
We examined the double red clump (RC) observed in the Galactic bulge, interpreted as a difference in distance ("X-shaped bulge scenario") or in chemical composition ("multiple population scenario"). To verify chemical differences between the RC groups, we performed low-resolution spectroscopy for RC and red giant branch (RGB) stars using Gemini-South/GMOS in three fields of the bulge, and collected diverse data from literature. We divided our sample stars not only into bright and faint RC groups, but also into bluer ([Fe/H] < -0.1) and redder ([Fe/H] > -0.1) groups following the recent u-band photometric studies. For the metal-poor stars, no statistically significant difference in CN index was detected between the bright and faint RC groups for all observed fields. However, we found, from cross-matching with high-resolution spectroscopic data, a sign of Na enhancement in the "metal-poor and bright" RC group compared to the "metal-poor and faint" group at (l,b)=(-1 deg,-8.5 deg). When the contributions of the RGB stars on the RC regimes are taken into account, the Na abundance difference between genuine RCs would correspond to approximately 0.23 dex, similar to globular cluster (GCs) with multiple populations. In contrast, the metal-rich stars do not show chemical differences between the bright and faint RCs. It implies that the double RC observed in the metal-poor component of the bulge might be linked to the multiple populations originated from GC-like subsystem, whereas that of the metal-rich component would have produced by the X-shaped structure. Our results support the previous studies suggesting composite nature of the Milky Way bulge.
Context. With the increasing number of discoveries of globular clusters in the inner Milky Way, the need for spectroscopic confirmation and further investigation of their stellar populations and chemodynamical properties has become crucial. Aims. Gran 5 is a newly reported low-mass globular cluster located close to the Galactic center, and it is thought to be an accreted object associated with the Gaia-Enceladus structure. This study aims to investigate the stellar populations of Gran 5 and their detailed chemical properties. Methods. We performed high-resolution near-infrared spectroscopy on seven stars in the field of Gran 5 using IGRINS on the Gemini-South telescope. Results. We identified six stars as cluster members and reveal that they are divided into two stellar populations with different metallicities, with mean [Fe/H] values of -0.76 dex and -0.55 dex, respectively. In addition, the chemodynamical properties of Gran 5 agree with those of in situ globular clusters. Conclusions. Our findings represent the first detection of two stellar populations with different metallicities in a low-mass globular cluster. This suggests that the metallicity variation in Gran 5 may have arisen from processes different from those in other globular clusters with metallicity variation, or that it may have lost a substantial amount of its initial mass during its evolution.
We present a new analytical model for the attenuation to Epoch of Reionization (EoR) galaxies by proximate neutral hydrogen gas. Many galaxy spectra in the EoR taken by JWST have shown a flux deficit at wavelengths just redward of the Lyman break, and this has been regarded as resulting from Ly$\alpha$ damping wing absorption by the increasing amount of neutral hydrogen in the line-of-sight. However, previous attenuation models for the intergalactic medium (IGM) commonly used in photometric redshift template-fitting codes assume that the Lyman break is rather sharp, which leads to systematic overestimation of photometric redshifts at $z>7$. In this letter, we build and empirically calibrate a new attenuation model that takes the increased Ly$\alpha$ damping wing absorption into account. Our model consists of the canonical IGM absorption and an additional absorption component due to dense neutral hydrogen gas clouds proximate to the galaxy, and we derive the redshift evolution of HI column density of the proximate clouds by calibrating the model using CANUCS JWST observations. The resulting total transmission curve resolves the photometric redshift bias at $z>7$, an improvement that is robust to choice of template-fitting code, template set, and photometric catalog used. The new attenuation model can be easily implemented in existing template-fitting codes, and significantly improves the photometric redshift performance in the EoR.
We propose a novel dust battery mechanism for generating seed magnetic fields in the early universe, in which charged dust grains are radiatively accelerated, inducing strong electric currents that subsequently generate magnetic fields. Our analysis demonstrates that this process is effective even at very low metallicities (approximately $ \sim 10^{-5} Z_\odot$), and capable of producing seed fields with significant amplitudes of $B \sim \rm \mu G$ around luminous sources over timescales of years to Myr and across spatial scales ranging from au to kpc. Crucially, we find that this mechanism is generically $\sim10^8$ times more effective than the radiatively-driven electron battery or Biermann battery in relatively cool gas ($\ll 10^{5}\,$K), including both neutral and ionized gas. Furthermore, our results suggest that, to first order, dissipation effects do not appear to significantly impede this process, and that it can feasibly generate coherent seed fields on macroscopically large ISM scales (much larger than turbulent dissipation scales or electron mean-free-paths in the ISM). These seed fields could then be amplified by subsequent dynamo actions to the observed magnetic fields in galaxies. Additionally, we propose a sub-grid model for integration into cosmological simulations, and the required electric-field expressions for magnetohydrodynamic-particle-in-a-cell (MHD-PIC) simulations that explicitly model dust dynamics. Finally, we explore the broad applicability of this mechanism across different scales and conditions, emphasizing its robustness compared to other known battery mechanisms.
We study the importance of several processes that influence the evolution of dust and its grain size distribution on spatially resolved scales in nearby galaxies. Here, we compiled several multi-wavelength observations for the nearby galaxies NGC628(M74), NGC5457(M101), NGC598(M33), and NGC300. We applied spatially resolved spectral energy distribution fitting to the latest iteration of infrared data to get constraints on the galaxy dust masses and the small-to-large grain abundance ratio. For comparison, we took the radial profiles of the stellar mass and gas mass surface density for NGC628 combined with its metallicity gradient in the literature to calibrate a single-galaxy simulation using the GADGET4-OSAKA code. The simulations include a parametrization to separate the dense and diffuse phases of the ISM where different dust-evolution mechanisms are in action. We find that our simulation can reproduce the radial profile of dust mass surface density but overestimates the SLR in NGC628. Changing the dust-accretion timescale has little impact on the dust mass or SLR, as most of the available metals are accreted onto dust grains at early times (< 3Gyr), except in the outer regions of the galaxy. This suggests we can only constrain the accretion timescale of galaxies at extremely low metallicities where accretion still competes with other mechanisms controlling the dust budget. The overestimation of the SLR likely results from (i) overly efficient shattering processes in the diffuse interstellar medium, which were calibrated to reproduce Milky Way-type galaxies and/or (ii) our use of a diffuse and dense gas density subgrid model that does not entirely capture the intricacies of the small-scale structure present in NGC628.
The observed prevalence of galaxies exhibiting bursty star formation histories (SFHs) at $z\gtrsim6$ has created new challenges and opportunities for understanding their formation pathways. The degenerate effects of the efficiency and burstiness of star formation on the observed UV luminosity function are separable by galaxy clustering. However, quantifying the timescales of burstiness requires more than just the continuum UV measurements. Here we develop a flexible semi-analytic framework for modeling both the amplitude of star formation rate (SFR) variations and their temporal correlation, from which the luminosity function and clustering can be derived for SFR indicators tracing different characteristic timescales (e.g., UV continuum and H$\alpha$ luminosities). Based on this framework, we study the prospect of using galaxy summary statistics to distinguish models where SFR fluctuations are prescribed by different power spectral density (PSD) forms. Using the Fisher matrix approach, we forecast the constraints on parameters in our PSD-based model that can be extracted from mock JWST observations of the UV and H$\alpha$ luminosity functions and clustering bias factors at $z\sim6$. These constraints demonstrate the feasibility of constraining the burstiness of high-$z$ galaxies solely from their one-point and two-point statistics and underscore the importance of combining tracers of both long-term and short-term SFR variations. Our flexible framework can be readily extended to characterize the SFH of high-redshift galaxies with a wider range of observational diagnostics.
Over the past decade, several millimeter interferometer programs have mapped the nearby star-forming galaxy M51 at a spatial resolution of ${\le}170$ pc. This study combines observations from three major programs: the PdBI Arcsecond Whirlpool Survey (PAWS), the SMA M51 large program (SMA-PAWS), and the Surveying the Whirlpool at Arcseconds with NOEMA (SWAN). The dataset includes the (1-0) and (2-1) rotational transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O isotopologues. The observations cover the $r{<}\rm 3\,kpc$ region including center and part of the disk, thereby ensuring strong detections of the weaker $^{13}$CO and C$^{18}$O lines. All observations are convolved in this analysis to an angular resolution of 4$''$, corresponding to a physical scale of ${\sim}$170 pc. We investigate empirical line ratio relations and quantitatively evaluate molecular gas conditions such as temperature, density, and the CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$). We employ two approaches to study the molecular gas conditions: (i) assuming local thermal equilibrium (LTE) to analytically determine the CO column density and $\alpha_{\rm CO}$, and (ii) using non-LTE modeling with RADEX to fit physical conditions to observed CO isotopologue intensities. We find that the $\alpha_{\rm CO}$ values {in the center and along the inner spiral arm} are $\sim$0.5 dex (LTE) and ${\sim}$0.1 dex (non-LTE) below the Milky Way inner disk value. The average non-LTE $\alpha_{\rm CO}$ is $2.4{\pm}0.5$ M$_\odot$ pc$^{-2}$ (K km s$^{-1}$)$^{-1}$. While both methods show dispersion due to underlying assumptions, the scatter is larger for LTE-derived values. This study underscores the necessity for robust CO line modeling to accurately constrain the molecular ISM's physical and chemical conditions in nearby galaxies.
The goal of this project is to construct an estimator for the masses of supermassive black holes in active galactic nuclei (AGNs) based on the broad Halpha emission line. We make use of published reverberation mapping data. We remeasure all Halpha time lags from the original data as we find that often the reverberation measurements are improved by detrending the light curves. We produce mass estimators that require only the Halpha luminosity and the width of the Halpha emission line as characterized by either the FWHM or the line dispersion. It is possible, on the basis of a single spectrum covering the Halpha emission line, to estimate the mass of the central supermassive black hole in AGNs, taking into account all three parameters believed to affect mass measurement: luminosity, line width, and Eddington ratio. The typical formal accuracy in such estimates is of order 0.2-0.3 dex.
Within the Milky Way (MW), younger stellar populations exhibit steeper (more negative) metallicity radial gradients; the origin of this trend remains debated. The FIRE-2 cosmological simulations of MW-mass galaxies show the same trend as the MW, which in FIRE-2 arises because the metallicity gradient of the interstellar medium (ISM), and thus of stars at birth, became steeper over time. We seek to understand this evolution in the context of inside-out radial growth of galaxies. Most FIRE-2 galaxies grew radially inside-out in both gas and stars; specifically, their surface density profiles, $\Sigma(R)$, became shallower over time. Combined with a realized superlinear (Kennicutt-Schmidt-like) relation between star formation and total gas density, the profile of the ratio $\Sigma_{\rm star}(R)/\Sigma_{\rm gas}(R)$ became shallower (flatter) over time. Thus, if metals stayed where they were injected into the ISM from stars, the metallicity gradient would become shallower over time, as some models predict. However, metallicity gradients in FIRE-2 became steeper over time, because of the additional effects of (radial) mixing of metals in the ISM. Specifically, the velocity dispersion and net radial advection of gas declined over time, as ISM turbulence decreased and the disk settled, leading to upside-down vertical growth. In FIRE-2, this evolution in metal mixing of gas associated with upside-down growth dominates over inside-out radial growth, causing the metallicity radial gradient of the ISM and of stars at birth to become steeper over time. We argue that this reflects the ISM history of the MW and of typical MW-mass galaxies.
We present new stellar population models, $\alpha$-MC, self-consistently taking into account non-solar $\rm [\alpha/Fe]$ abundances for both isochrones and stellar spectra. The $\alpha$-MC models are based on $\alpha$-enhanced MIST isochrones and C3K spectral libraries, which are publicly available in FSPS. Our new models cover a wide range of ages ($\rm \log (age/yr) = 5.0 - 10.3$), metallicities ($\rm [Fe/H]=[-2.5,+0.5]$ in steps of 0.25, $\rm [\alpha/Fe]=-0.2,+0.0,+0.2,+0.4,+0.6$), and wavelengths ($0.1-2.5\,\rm \mu m$). We investigate the separate and combined effects of $\alpha$-enhanced isochrones and stellar spectral libraries on simple stellar populations (SSPs), including their broadband colors, spectral indices, and full spectra. We find that the primary effect of $\alpha$-enhancement in isochrones is to lower the overall continuum levels and redden the continuum shapes, while $\alpha$-enhancement in stellar spectra mainly affects individual spectral lines. At constant $\rm [Fe/H]$, $\alpha$-enhancement has significant impacts on the broadband colors by $\rm \sim 0.1-0.4\,mag$ across all ages ($\rm 0.01 - 10\,Gyr$). The effects of $\alpha$-enhancement on colors at fixed $\rm [Z/H]$ are smaller, by $\rm \sim 0.1-0.2\,mag$. The spectral indices involving $\alpha$-elements, Ca4227 and Mg b, increase with $\rm [\alpha/Fe]$ (both at fixed $\rm [Fe/H]$ and fixed $\rm [Z/H]$) due to enhanced $\alpha$-abundances. At constant $\rm [Fe/H]$, $\alpha$-enhancement weakens most Fe-sensitive and Hydrogen Balmer lines. Our new self-consistent $\alpha$-enhanced models will be essential in deriving accurate physical properties of high-redshift galaxies, where $\alpha$-enhancement is expected to be common.
Despite the ubiquity of clumpy star-forming galaxies at high-redshift, the origin of clumps are still largely unconstrained due to the limited observations that can validate the mechanisms for clump formation. We postulate that if clumps form due to the accretion of metal-poor gas that leads to violent disk instability, clumpy galaxies should have lower gas-phase metallicities compared to non-clumpy galaxies. In this work, we obtain the near-infrared spectrum for 42 clumpy and non-clumpy star-forming galaxies of similar masses, SFRs, and colors at $z\approx0.7$ using the Gemini Near-Infrared Spectrograph (GNIRS) and infer their gas-phase metallicity from the {\nii} and {\halpha} line ratio. We find that clumpy galaxies have lower metallicities compared to non-clumpy galaxies, with an offset in the weighted average metallicity of $0.07\pm0.02$ dex. We also find an offset of $0.06\pm0.02$ dex between clumpy and non-clumpy galaxies in a comparable sample of 23 star-forming galaxies at $z\approx1.5$ using existing data from the FMOS-COSMOS survey. Similarly, lower {\nii}/{\halpha} ratio are typically found in galaxies that have more of their $\mathrm{UV_{rest}}$ luminosity originating from clumps, suggesting that \enquote{clumpier} galaxies are more metal poor. We also derive the intrinsic velocity dispersion and line-of-sight rotational velocity for galaxies from the GNIRS sample. The majority of galaxies have $\sigma_0/v_c \approx 0.2$, with no significant difference between clumpy and non-clumpy galaxies. Our result indicates that clump formation may be related to the inflow of metal-poor gas; however, the process that forms them does not necessarily require significant, long-term kinematic instability in the disk.
The James Webb Space Telescope (JWST) has revealed unexpectedly rapid galaxy assembly in the early universe, in tension with models of star and galaxy formation. In the gas conditions typical of early galaxies, particularly their low abundances of heavy elements (metals) and dust, the star-formation process is poorly understood. Some models predict that stars form in atomic gas at low metallicity, in contrast to forming in molecular gas as observed in higher-metallicity galaxies. To understand the very high star-formation rates at early epochs, it is necessary to determine whether molecular gas formation represents a bottleneck to star formation, or if it is plentiful even at extremely low metallicity. Despite repeated searches, star-forming molecular gas has not yet been observed in any galaxy below 7% of the Solar metallicity, leaving the question of how stars form at lower metallicities unresolved. Here, we report the detection of rotationally excited emission from molecular hydrogen in the star-forming region of the nearby, 3% Solar metallicity galaxy Leo P with the MIRI-MRS instrument onboard JWST. These observations place a lower limit on the molecular gas content of Leo P and, combined with our upper limit on carbon monoxide emission from a deep search of this galaxy, demonstrate that MIRI-MRS is sensitive to much smaller molecular gas masses at extremely low metallicity compared to the traditional observational tracer. This discovery pushes the maximum metallicity at which purely atomic gas may fuel star formation a factor of two lower, providing crucial empirical guidance for models of star formation in the early universe.
Active galactic nuclei (AGN) feedback from supermassive black holes (SMBHs) at the centers of galaxy clusters plays a key role in regulating star formation and shaping the intracluster medium (ICM), often manifesting through prominent X-ray cavities embedded in the cluster's hot atmosphere. Here we show that X-ray cavities arise naturally due to AGN feedback in TNG-Cluster. This is a new suite of magnetohydrodynamic cosmological simulations of galaxy formation and evolution, and hence of galaxy clusters, whereby cold dark matter, baryon dynamics, galactic astrophysics, and magnetic fields are evolved together consistently. We construct mock Chandra X-ray observations of the central regions of the 352 simulated clusters at $z=0$ and find that $\sim$39 per cent contain X-ray cavities. Identified X-ray cavities vary in configuration (single, pairs, or multiples), with some still attached to SMBHs, while others have buoyantly risen. Their size ranges from a few to several tens of kpc. In terms of gas physical properties, TNG-Cluster X-ray cavities are underdense compared to the surrounding halo and filled with hot gas ($\sim$10$^8$K); 25 per cent of them are surrounded by an X-ray bright and compressed rim associated with a weak shock (Mach number $\sim 1.5$). Clusters exhibiting X-ray cavities are preferentially strong or weak cool-cores, are dynamically relaxed, and host SMBHs accreting at low Eddington rates. We show that TNG-Cluster X-ray cavities originate from episodic, wind-like energy injections from central AGN. Our results illustrate the existence and diversity of X-ray cavities simulated in state-of-the-art models within realistic cosmological environments and show that these can form without necessarily invoking bipolar, collimated, or relativistic jets.