Galaxy Formation Research Articles

Anatomy of the Star Formation in a Tidally Disturbed Disk Galaxy: NGC 3718

We present a UV, optical, and near-IR study of the star-forming complexes in the nearby peculiar galaxy NGC 3718, using Ultraviolet Imaging telescope, Galaxy Evolution Explorer, Spitzer, and DECam Legacy Survey imaging data. The galaxy has a disturbed optical morphology owing to the multiple tidal arms, the warped disk, and the prominent curved dust lanes, but in the near-IR it appears to be a bulge-dominated galaxy. Its disturbed morphology makes it an excellent case to study star formation in a tidally disturbed galaxy that may have undergone a recent minor merger. To study the distribution and properties of the star-forming clumps (SFCs), we divided the galaxy within the R 25 (B band) radius into three parts—the upper, central, and lower regions. Using the UV band images, we investigated the warped star-forming disk, the extended tidal arms, and the distribution and sizes of the 182 SFCs. Their distribution is 49, 60, and 73 in the galaxy’s upper, central, and lower regions, respectively. We determined the UV color, star formation rates (SFRs), star formation density (ΣSFR), and ages of the SFCs. The central disk of the galaxy has a larger mean ΣSFR that is ∼3.3 and ∼1.6 times higher than the upper and lower regions, respectively. We also find that the SFCs in the central disk are older than those in the tidal arms. Our study thus shows that minor mergers can trigger the inside-out growth of galaxy disks, where the younger SFCs are in outer tidal arms and not in the inner disk.

Identification of High-redshift Galaxy Overdensities in GOODS-N and GOODS-S

We conduct a systematic search for high-redshift galaxy overdensities at 4.9 < z spec < 8.9 in both the Great Observatories Origins Deep Survey (GOODS)-N and GOODS-S fields using James Webb Space Telescope/Near-Infrared Camera (JWST/NIRCam) imaging from the JWST Advanced Deep Extragalactic Survey and JWST Extragalactic Medium-band Survey in addition to JWST/NIRCam wide field slitless spectroscopy from the First Reionization Epoch Spectroscopic Complete Survey. High-redshift galaxy candidates are identified using Hubble Space Telescope + JWST photometry spanning λ = 0.4–5.0 μm. We confirmed the redshifts for roughly a third of these galaxies using JWST spectroscopy over λ = 3.9–5.0 μm through identification of either Hα or OIIIλ5008 around the best-fit photometric redshift. The rest-ultraviolet magnitudes and continuum slopes of these galaxies were inferred from the photometry: the brightest and reddest objects appear in more dense environments and thus are surrounded by more galaxy neighbors than their fainter and bluer counterparts, suggesting accelerated galaxy evolution within overdense environments. We find 17 significant (δ gal ≥ 3.04, N gal ≥ 4) galaxy overdensities across both fields (seven in GOODS-N and 10 in GOODS-S), including the two highest redshift spectroscopically confirmed galaxy overdensities to date at zspec=7.954 and zspec=8.222 (representing densities around ∼6 and ∼12 times that of a random volume). We estimate the total halo mass of these large-scale structures to be 11.5≤log10Mhalo/M⊙≤13.4 using an empirical stellar mass-to-halo mass relation, which are likely underestimates as a result of incompleteness. These protocluster candidates are expected to evolve into massive galaxy clusters with log10Mhalo/M⊙≳14 by z = 0.

A Break in the Size–Stellar Mass Relation: Evidence for Quenching and Feedback in Dwarf Galaxies

Mapping stars and gas in nearby galaxies is fundamental for understanding their growth and the impact of their environment. This issue is addressed by comparing the stellar “edges” of galaxies D stellar, defined as the outermost diameter where in situ star formation significantly drops, with the gaseous distribution parameterized by the neutral atomic hydrogen diameter measured at 1 M ⊙ pc−2, D HI. By sampling a broad H i mass range 105 M ⊙ < M HI < 1011 M ⊙, we find several dwarf galaxies with M HI < 109 M ⊙ from the field and Fornax Cluster that are distinguished by D stellar ≫ D HI. For the cluster dwarfs, the average H i surface density near D stellar is ∼0.3 M ⊙ pc−2, reflecting the impact of quenching and outside-in gas removal from ram pressure and tidal interactions. In comparison, D stellar/D HI ranges between 0.5 and 2 in dwarf field galaxies, consistent with the expectations from stellar feedback. Only more massive disk galaxies in the field can thus be characterized by the common assumption that D stellar ≲ D HI. We discover a break in the D stellar–M ⋆ relation at m break ∼ 4 × 108 M ⊙ that potentially differentiates the low-mass regime, where the influence of stellar feedback and environmental processes more prominently regulates the sizes of nearby galaxies. Our results highlight the importance of combining deep optical and H i imaging for understanding galaxy evolution.

Evidence for Intrinsic X-Ray Weakness among Red Quasars at Cosmic Noon

Quasar feedback is a key ingredient in shaping galaxy evolution. A rare population of extremely red quasars (ERQs) at z = 2−3 are often associated with high-velocity [O iii]λ5008 outflows and may represent sites of strong feedback. In this paper, we present an X-ray study of 50 ERQs to investigate the link between the X-ray and outflow properties of these intriguing objects. Using hardness ratio analysis, we confirm that the ERQs are heavily obscured systems with gas column density reaching N H = 1023−24 cm−2. We identify 20 X-ray-nondetected ERQs at high mid-infrared (MIR) luminosities of ν L ν,6 μm ≳ 3 × 1046 erg s−1. By stacking the X-ray observations, we find that the nondetected ERQs are on average underluminous in X-rays by a factor of ∼10 for their MIR luminosities. We consider such X-ray weakness to be due to both heavy gas absorption and intrinsic factors. Moreover, we find that the X-ray-weak sources also display higher-velocity outflows. One option to explain this trend is that weaker X-rays facilitate more vigorous line-driven winds, which then accelerate the [O iii]-emitting gas to kiloparsec scales. Alternatively, super-Eddington accretion could also lead to intrinsic X-ray weakness and more powerful continuum-driven outflow.

Illuminating the Incidence of Extraplanar Dust Using Ultraviolet Reflection Nebulae with GALEX

Circumgalactic dust grains trace the circulation of mass and metals between star-forming regions and gaseous galactic halos, giving insight into feedback and tidal stripping processes. We perform a search for ultraviolet (UV) reflection nebulae produced by extraplanar dust around 551 nearby (D < 100 Mpc), edge-on disk galaxies using archival near-ultraviolet and far-ultraviolet images from the Galaxy Evolution Explorer (GALEX), accounting for the point-spread function (FWHM = 4″–5″). We detect extraplanar emission ubiquitously in stacks of galaxies binned by morphology and star formation rate, with scale heights of h h = 1–2.3 kpc and ≈10% of the total (reddened) flux in the galaxy found beyond the B-band isophotal level of μ B = 25 mag arcsec−2. This emission is detected in 7% of the individual galaxies, and an additional one-third have at least 5% of their total flux found beyond μ B = 25 mag arcsec−2 in a disk component. The extraplanar luminosities and colors are consistent with reflection nebulae rather than stellar halos and indicate that, on average, disk galaxies have an extraplanar dust mass of 5%–15% of that in their interstellar medium. This suggests that recycled material composes at least a third of the inner circumgalactic medium (R < 10 kpc) in ∼L* galaxies.

High-quality Extragalactic Legacy-field Monitoring (HELM) with DECam: Project Overview and First Data Release

High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4 m Blanco telescope. Since 2019 February, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the (u)gri(z) bands, over a collective sky area of ∼38 deg2. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependencies on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first ∼3.5 yr of light curves during 2019A–2022A.

Active Galactic Nucleus Jet-inflated Bubbles as Possible Origin of Odd Radio Circles

Odd radio circles (ORCs) are newly discovered extragalactic radio objects with an unknown origin. In this work, we carry out three-dimensional cosmic-ray (CR) magnetohydrodynamic simulations using the FLASH code and predict the radio morphology of end-on active galactic nucleus (AGN) jet-inflated bubbles considering hadronic emission. We consider CR proton (CRp)-dominated jets as they tend to inflate oblate bubbles, promising to reproduce the large inferred sizes of the ORCs when viewed end-on. We find that powerful and long-duration CRp-dominated jets can create bubbles with similar sizes (∼300–600 kpc) and radio morphology (circular and edge-brightened) to the observed ORCs in low-mass (M vir ∼ 8 × 1012 − 8 × 1013 M ⊙) halos. Given the same amount of input jet energy, longer-duration (thus lower-power) jets tend to create larger bubbles since high-power jets generate strong shocks that carry away a significant portion of the jet energy. The edge-brightened feature of the observed ORCs is naturally reproduced due to efficient hadronic collisions at the interface between the bubbles and the ambient medium. We further discuss the radio luminosity, X-ray detectability, and the possible origin of such strong AGN jets in the context of galaxy evolution. We conclude that end-on CR-dominated AGN bubbles could be a plausible scenario for the formation of ORCs.

The Three Hundred: The existence of massive dark matter-deficient satellite galaxies in cosmological simulations

The observation of a massive galaxy with an extremely low dark matter content (i.e. NGC 1277) has posed questions about how such objects form and evolve in a hierarchical universe. We here report on the finding of several massive, dark matter-deficient galaxies in a set of 324 galaxy clusters theoretically modelled by means of full-physics hydrodynamical simulations. We first focus on two example galaxies selected amongst the most massive and dark matter-deficient ones. By tracing the evolution of these galaxies, we find that their lack of dark matter is a result of multiple pericentre passages. While orbiting their host halo, tidal interactions gradually strip away dark matter while preserving the stellar component. A statistical analysis of all massive satellite galaxies in the simulated clusters shows that the stellar-to-total mass ratio today is strongly influenced by the number of orbits and the distance at pericentres. Galaxies with more orbits and closer pericentres are more dark matter-deficient. Additionally, we find that massive, dark matter-deficient galaxies at the present day are either the remnants of very massive galaxies at infall or former central galaxies of infalling groups. We conclude that such massive yet dark matter-deficient galaxies exist and are natural by-products of typical cluster galaxy evolution, with no specific requirement for an exotic formation scenario.

Using dust to constrain dark matter models

Abstract In this paper, we use hydrodynamic zoom-in simulations of Milky Way-type haloes to explore using dust as an observational tracer to discriminate between cold and warm dark matter universes. Comparing a cold and 3.5keV warm dark matter particle model, we tune the efficiency of galaxy formation in our simulations using a variable supernova rate to create Milky Way systems with similar satellite galaxy populations while keeping all other simulation parameters the same. Cold dark matter, having more substructure, requires a higher supernova efficiency than warm dark matter to achieve the same satellite galaxy number. These different supernova efficiencies create different dust distributions around their host galaxies, which we generate by post-processing the simulation output with the powderday codebase. Analysing the resulting dust in each simulation, we find ∼4.5 times more dust in our cold dark matter Milky Way halos compared with warm dark matter. The distribution of dust out to R200c is then explored, revealing that the warm dark matter simulations are noticeably less concentrated than their cold dark matter counterparts, although differences in substructure complicate the comparison. Our results indicate that dust is a possible unique probe to test theories of dark matter.

Effects of Bursty Star Formation on [C ii] Line Intensity Mapping of High-redshift Galaxies

Bursty star formation—a key prediction for high-redshift galaxies from cosmological simulations explicitly resolving stellar feedback in the interstellar medium—has recently been observed to prevail among galaxies at redshift z ≳ 6. Line intensity mapping (LIM) of the 158 μm [C ii] line as a star formation rate (SFR) indicator offers unique opportunities to tomographically constrain cosmic star formation at high redshift, in a way complementary to observations of individually detected galaxies. To understand the effects of bursty star formation on [C ii] LIM, which have remained unexplored in previous studies, we present an analytic modeling framework for high-z galaxy formation and [C ii] LIM signals that accounts for bursty star formation histories induced by delayed supernova feedback. We use it to explore and characterize how bursty star formation can impact and thus complicate the interpretation of the [C ii] luminosity function and power spectrum. Our simple analytic model indicates that bursty star formation mainly affects low-mass galaxies by boosting their average SFR and [C ii] luminosity, and in the [C ii] power spectrum it can create a substantial excess in the large-scale clustering term. This distortion results in a power spectrum shape that cannot be explained by invoking a mass-independent logarithmic scatter. We conclude that burstiness must be accounted for when modeling and analyzing [C ii] data sets from the early Universe, and that in the extreme, the signature of burstiness may be detectable with first-generation experiments such as TIME, CONCERTO, and CCAT-DSS.

Percolation Statistics in the MillenniumTNG Simulations

The statistical analysis of cosmic large-scale structure is most often based on simple two-point summary statistics, like the power spectrum or the two-point correlation function of a sample of galaxies or other types of tracers. In contrast, topological measures of clustering are also sensitive to higher-order correlations and thus offer the prospect to access additional information that may harbor important constraining power. We here revisit one such geometric measure of the cosmic web in the form of the so-called percolation analysis, using the recent MillenniumTNG simulation suite of the ΛCDM paradigm. We analyze continuum percolation statistics both for high-resolution dark matter particle distributions and for galaxy mock catalogs from a semianalytic galaxy formation model within a periodic simulation volume of 3000 Mpc on a side. For comparison, we also investigate the percolation statistics of random particle sets and neutrino distributions with two different summed particle masses. We find that the percolation statistics of the dark matter distribution evolves strongly with redshift and thus clustering strength, yielding a progressively lower percolation threshold toward later times. However, there is a sizable residual dependence on numerical resolution, which we interpret as a residual influence of different levels of shot noise. This is corroborated by our analysis of galaxy mock catalogs, whose results depend on sampling density more strongly than on galaxy selection criteria. While this limits the discriminative power of percolation statistics, our results suggest that it still remains useful as a complementary cosmological test when controlled for sampling density.

Galaxy Assembly Bias in the Stellar-to-halo Mass Relation for Red Central Galaxies from SDSS

We report evidence of galaxy assembly bias—the correlation between galaxy properties and biased secondary halo properties at fixed halo mass (M H)—in the stellar-to-halo mass relation for red central galaxies from the Sloan Digital Sky Survey. In the M H = 1011.5–1013.5 h −1 M ⊙ range, central galaxy stellar mass (M *) is correlated with the number density of galaxies within 10 h −1 Mpc (δ 10), a common proxy for halo formation time. This galaxy assembly bias signal is also present when M H, M *, and δ 10 are substituted with group luminosity, galaxy luminosity, and metrics of the large-scale density field. To associate differences in δ 10 with variations in halo formation time, we fitted a model that accounts for (1) errors in the M H measured by the J. L. Tinker group catalog and (2) the level of correlation between halo formation time and M * at fixed M H. Fitting of this model yields that (1) errors in M H are ∼0.15 dex and (2) halo formation time and M * are strongly correlated (Spearman’s rank correlation coefficient ∼0.85). At fixed M H, variations of ∼0.4 dex in M * are associated with ∼1–3 Gyr variations in halo formation time and galaxy formation time (from stellar population fitting). These results are indicative that halo properties other than M H can impact central galaxy assembly.

Baryonification extended to thermal Sunyaev Zel’dovich

Baryonification algorithms model the impact of galaxy formation and feedback on the matter field in gravity-only simulations by adopting physically motivated parametric prescriptions. In this paper, we extend these models to describe gas temperature and pressure, allowing for a self-consistent modelling of the thermal Sunyaev-Zel’dovich effect, weak gravitational lensing, and their cross-correlation, down to small scales. We validate our approach by showing that it can simultaneously reproduce the electron pressure, gas, stellar, and dark matter power spectra as measured in all BAHAMAS hydrodynamical simulations. Specifically, with only two additional free parameters, we can fit the electron pressure auto- and cross-power spectra at 10% while reproducing the suppression in the matter power spectrum induced by baryons at the per cent level, for different active galactic nuclei (AGN) feedback strengths in BAHAMAS. Furthermore, we reproduce BAHAMAS convergence and thermal Sunyaev Zel’dovich angular power spectra within 1% and 10% accuracy, respectively, down to ℓ = 5000. When used jointly with cosmological rescaling algorithms, the baryonification presented here allows for a fast and accurate exploration of cosmological and astrophysical scenarios. Therefore, it can be employed to create mock catalogues, lightcones, and large training sets for emulators aimed at interpreting forthcoming multi-wavelength observations of the large-scale structure of the Universe.

The Brown Dwarf Kinematics Project (BDKP). VI. Ultracool Dwarf Radial and Rotational Velocities from SDSS/APOGEE High-resolution Spectroscopy

We present precise measurements of radial (RV) and projected rotational ( vsini ) velocities of a sample of 258 M6 to L2 dwarfs with multiepoch, high-resolution (λ/Δλ = 22,500), near-infrared (1.514–1.696 μm) spectroscopic observations reported in the Apache Point Observatory Galactic Evolution Experiment Data Release 17. The spectra were modeled using a Markov Chain Monte Carlo forward-modeling method, which achieved median precisions of σ RV = 0.4 km s−1 and σvsini = 1.1 km s−1. One-half of our sample (138 sources) are previously known members of nearby young clusters and moving groups, and we identified three new kinematic members of the Argus or Carina Near moving groups, 2MASS J05402570+2448090, 2MASS J14093200+4138080, and 2MASS J21272531+5553150. Excluding these sources, we find that the majority of our sample have kinematics consistent with the Galactic thin disk, and 11 sources are associated with the intermediate thin/thick disk. The field sample has a velocity dispersion of 38.2 ± 0.3 km s−1, equivalent to an age of 3.30 ± 0.19 Gyr based on empirical age–velocity dispersion relations, and a median vsini of 17 km s−1. For 172 sources with multiepoch observations, we identified 37 as having significant RV variations, and determined preliminary orbit parameters for 26 sources with four or more epochs, nine of which are short-period binary candidates. For 40 sources with photometric variability periods from the literature less than 5 days and vsini > 20 km s−1, we find a decline in projected radii ( Rsini ) with age congruent with evolutionary models. Finally, we also present multiepoch RV and vsini measurements for additional 444 candidate ultracool dwarfs.

Data-driven Stellar Intrinsic Colors and Dust Reddenings for Spectrophotometric Data: From the Blue-edge Method to a Machine Learning Approach

Intrinsic colors (ICs) of stars are essential for studies on both stellar physics and dust reddening. In this work, we developed an XGBoost model to predict ICs with the atmospheric parameters T eff, logg , and [M/H]. The model was trained and tested for three colors at Gaia and Two Micron All Sky Survey bands with 1,040,446 low-reddening sources. The atmospheric parameters were determined by the Gaia DR3 General Stellar Parameterizer from Photometry (GSP-phot) module and were validated by comparing with the Apache Point Observatory Galactic Evolution Experiment and Large Sky Area Multi-Object fiber Spectroscopic Telescope. We further confirmed that the biases in GSP-phot parameters, especially for [M/H], do not present a significant impact on the IC prediction. The generalization error of the model estimated by the test set is 0.014 mag for (GBP−GRP)0 , 0.050 mag for (GBP−KS)0 , and 0.040 mag for (J−KS)0 . The model was applied to a sample containing 5,714,528 reddened stars with stellar parameters from R. Andrae et al. to calculate ICs and reddenings. The high consistency in the comparison of E(J − K S) between our results and literature values further validates the accuracy of the XGBoost model. The variation of E(G BP − K S)/E(G BP − G RP), a representation of the extinction law, with Galactic longitude is found on large scales. This work preliminarily presents the feasibility and accuracy of the machine learning approach for IC and dust reddening calculation, whose products could be widely applied to spectrophotometric data. The data sets and trained model can be accessed via Zenodo, doi:10.5281/zenodo.12787594. Models for more bands will be completed in following work.

Testing the sources of the peculiar abundances in globular clusters

This work aims to analyze some of the polluters proposed in the self-enrichment scenarios put forward to explain the multiple populations in globular clusters (GCs), extending previous studies. Three scenarios with different polluter stars were tested: asymptotic giant branch stars (AGBs), high-mass interacting binaries (IBs), and fast rotating massive stars (FRMSs). With abundance data available from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey and ΔY estimates from precise Hubble Space Telescope (HST) photometry, twenty-six clusters were studied, increasing the number of clusters in previous studies by more than a factor of three. We also included the study of the abundances of N, C, Mg, and Al, extending previous studies that mainly focused on the abundances of He, O, and Na. In addition, we constructed an empirical model to test whether one could explain the chemical signatures of the “enriched” population of GC stars with a fixed source and dilution process based on empirical data. In agreement with work by other authors, we found that the proposed polluters can generally predict the qualitative abundance patterns in GC stars and in some cases quantitatively predict some elements, but in most cases when we compare the model yields with the observations, we find that they cannot explain the entire set of observed abundance patterns. The empirical model succeeds in reproducing the abundances of Al for a given ΔY (and vice versa), showing that there is a direct relationship between Al and He, with one increasing proportionally to the other. However, the empirical model fails to reproduce the observed abundances of Na and N, in agreement with the results of previous works. The observed decoupling between the maximum abundances of CNO-cycle elements such as N and Na with those of Al and He provides new information and constraints for future models and could take us a step closer to understanding the origin of the peculiar abundance variations of GC stars.

Mapping Multiphase Metals in Star-forming Galaxies: A Spatially Resolved UV+Optical Study of NGC 5253

We present a pioneering, spatially resolved, multiphase gas abundance study on the blue compact dwarf galaxy NGC 5253, targeting 10 star-forming (SF) clusters inside six far-UV (FUV) Hubble Space Telescope/Cosmic Origins Spectrograph (COS) pointings with cospatial optical Very Large Telescope/MUSE observations throughout the galaxy. The SF regions span a wide range of ages (1–15 Myr) and are distributed at different radii (50–230 pc). We performed a robust absorption-line profile fitting on the COS spectra, covering 1065–1430 Å in the FUV, allowing an accurate computation of neutral-gas abundances for 13 different ions sampling eight elements. These values were then compared with the ionized-gas abundances, measured using the direct method on MUSE integrated spectra inside analog COS apertures. Our multiphase, spatially resolved comparisons find abundances, which are lower in the neutral gas than the ionized gas by 0.22, 0.80, and 0.58 dex for log(O/H), log(N/H), and log(N/O), respectively. We modeled the chemical abundance distributions and evaluated correlations as a function of the radius and age. It was found that, while N, O, and N/O abundances decrease as a function of age in the ionized gas, they increase with age in the neutral gas. No strong correlations for N, O, or N/O were observed as a function of the radius. The N/O and N/H offsets between the phases were found to decrease with age, providing evidence that chemical enrichment happens differentially, first in the ionized-gas phase around 2–5 Myrs (due to N-rich Wolf–Rayet stars) and then mixing out into the cold neutral gas on longer timescales of 10–15 Myr.

On the Cosmic Variance of the Merger Rate Density of Binary Neutron Stars

The cosmic variance on the star formation history may lead to bias in the merger rate density estimation of binary neutron star (BNS) mergers by compact binary population synthesis. In this paper, we take advantage of the large box size of the Millennium Simulation combined with the semianalytic galaxy formation model GABE and the parameterized population binary star evolution model to examine how much effect the cosmic variance will introduce on the estimation of the merger rate density of BNS mergers. We find that for subbox sizes of 100 and 200 Mpc, the variance of merger rate density σ R /R at different redshifts is about 23%–35% and 13%–20%, respectively. On the one hand, as for the variance of the detection rate on BNS mergers with the current LIGO–Virgo–KAGRA (LVK) detector network, this value is very small at ≲10%, which indicates ignoring the cosmic variance is reasonable for estimating the merger rate density from current LVK observations. On the other hand, with next-generation gravitational wave detectors, it is possible to localize BNS mergers within subboxes possessing a length of 40 Mpc for a source redshift z s < 0.2. In such a small box, the cosmic variance of the merger rate density is significant, i.e., the value of σ R /R is about ∼55%. This hints that estimating the merger rate density of BNS in different sky areas may provide useful information on the cosmic variance.

Constraining the Initial Mass Function via Stellar Transients

The stellar initial mass function (IMF) represents a fundamental quantity in astrophysics and cosmology describing the mass distribution of stars from low mass all the way up to massive and very massive stars. It is intimately linked to a wide variety of topics, including stellar and binary evolution, galaxy evolution, chemical enrichment, and cosmological reionization. Nonetheless, the IMF still remains highly uncertain. In this work, we aim to determine the IMF with a novel approach based on the observed rates of transients of stellar origin. We parametrize the IMF with a simple but flexible Larson shape, and insert it into a parametric model for the cosmic UV luminosity density, local stellar mass density, type Ia supernova (SN Ia), core-collapse supernova (CCSN), and long gamma-ray burst (LGRB) rates as a function of redshift. We constrain our free parameters by matching the model predictions to a set of empirical determinations for the corresponding quantities via a Bayesian Markov Chain Monte Carlo method. Remarkably, we are able to provide an independent IMF determination with a characteristic mass mc=0.10−0.08+0.24M⊙ and high-mass slope ξ=−2.53−0.27+0.24 that are in accordance with the widely used IMF parameterizations (e.g., Salpeter, Kroupa, Chabrier). Moreover, the adoption of an up-to-date recipe for the cosmic metallicity evolution allows us to constrain the maximum metallicity of LGRB progenitors to Zmax=0.12−0.05+0.29Z⊙. We also find which progenitor fraction actually leads to SN Ia or LGRB emission (e.g., due to binary interaction or jet-launching conditions), put constraints on the CCSN and LGRB progenitor mass ranges, and test the IMF universality. These results show the potential of this kind of approach for studying the IMF, its putative evolution with the galactic environment and cosmic history, and the properties of SN Ia, CCSN, and LGRB progenitors, especially considering the wealth of data incoming in the future.

The FLAMINGO project: a comparison of galaxy cluster samples selected on mass, X-ray luminosity, Compton-Y parameter, or galaxy richness

ABSTRACT Galaxy clusters provide an avenue to expand our knowledge of cosmology and galaxy evolution. Because it is difficult to accurately measure the total mass of a large number of individual clusters, cluster samples are typically selected using an observable proxy for mass. Selection effects are therefore a key problem in understanding galaxy cluster statistics. We make use of the $(2.8~\rm {Gpc})^3$ FLAMINGO hydrodynamical simulation to investigate how selection based on X-ray luminosity, thermal Sunyaev–Zeldovich effect or galaxy richness influences the halo mass distribution. We define our selection cuts based on the median value of the observable at a fixed mass and compare the resulting samples to a mass-selected sample. We find that all samples are skewed towards lower mass haloes. For X-ray luminosity and richness cuts below a critical value, scatter dominates over the trend with mass and the median mass becomes biased increasingly low with respect to a mass-selected sample. At $z\le 0.5$, observable cuts corresponding to median halo masses between $M_\text{500c}=10^{14}$ and $10^{15}~\rm {{\rm M}_{\odot }}$ give nearly unbiased median masses for all selection methods, but X-ray selection results in biased medians for higher masses. For cuts corresponding to median masses $\lt 10^{14}$ at $z\le 0.5$ and for all masses at $z\ge 1$, only Compton-Y selection yields nearly unbiased median masses. Importantly, even when the median mass is unbiased, the scatter is not because for each selection the sample is skewed towards lower masses than a mass-selected sample. Each selection leads to a different bias in secondary quantities like cool-core fraction, temperature, and gas fraction.