Search for: All records

Context.Submillimeter galaxies (SMGs) constitute a key population of bright star-forming galaxies at high-redshift. These galaxies challenge galaxy formation models, particularly regarding the reproduction of their observed number counts and redshift distributions. Furthermore, although SMGs contribute significantly to the cosmic star formation rate density (SFRD), their precise role remains uncertain. Upcoming surveys, such as the Ultra Deep Survey with the TolTEC camera, are expected to offer valuable insights into SMG properties and their broader impact in the Universe. Aims.Robust modeling of SMGs in a cosmological representative volume is necessary to investigate their nature in preparation for next-generation submillimeter surveys. Here, we test different parametric models for SMGs in large-volume hydrodynamical simulations, assess their contribution to the SFRD, and build expectations for future submillimeter surveys. Methods.We implement and test parametric relations derived from radiative transfer calculations across three cosmological simulation suites: EAGLE, IllustrisTNG, and FLAMINGO. We place particular emphasis on the FLAMINGO simulations due to their large volume and robust statistical sampling of SMGs. Based on the model that best reproduces observational number counts, we forecast submillimeter fluxes within the simulations, analyze the properties of SMGs, and evaluate their evolution over cosmic time. Results.Our results show that the FLAMINGO simulation reproduces the observed redshift distribution and source number counts of SMGs without requiring a top-heavy initial mass function. On the other hand, the EAGLE and IllustrisTNG simulations show a deficit of bright SMGs. We find that SMGs with S₈₅₀ > 1 mJy contribute up to ∼27% of the cosmic SFRD atz ∼ 2.6 in the FLAMINGO simulation, which is consistent with recent observations. Flux density functions reveal a rise in SMG abundance fromz = 6 toz = 2.5 that is followed by a sharp decline in the number of brighter SMGs fromz = 2.5 toz = 0. Leveraging the SMG population in FLAMINGO, we forecast that the TolTEC UDS will detect ∼80 000 sources over 0.8 deg²at 1.1 mm (at the 4σdetection limit), capturing about 50% of the cosmic SFRD atz ∼ 2.5. 

Aims.We investigate the physical properties and redshift evolution of simulated galaxies residing in unvirialized cosmic structures (i.e., protoclusters) at cosmic noon, to understand the influence of the environment on galaxy formation. This work is intended to build clear expectations for the ongoing ODIN (One-hundred-deg²DECam Imaging in Narrowbands) survey, which is mapping large-scale structures atz= 2.4,3.1, and 4.5 using Lyα-emitting galaxies (LAEs) as tracers. Methods.From the IllustrisTNG simulations, we define subregions centered on the most massive clusters ranked by total stellar mass atz= 0 and study the properties of galaxies within, including those of LAEs. To model the LAE population, we take a semi-analytical approach that assigns Lyαluminosity and equivalent width based on the UV luminosities to galaxies in a probabilistic manner. We investigate stellar mass, star formation rate (SFR), major merger events, and specific star formation rate of the population of star-forming galaxies and LAEs in the field- and protocluster environment and trace their evolution across cosmic time betweenz= 0−4. Results.We find that the overall shape of the UV luminosity function in simulated protocluster environments is characterized by a substantially shallower faint-end slope and a large excess on the bright end, signaling different formation histories for galaxies therein. The difference is milder for the Lyαluminosity function. While protocluster galaxies follow the same SFR-M_★scaling relation as average field galaxies, a larger fraction appears to have experienced major mergers in the last 200 Myr and as a result shows enhanced star formation at a ≈60% level, leading to a flatter distribution in both SFR and M_★relative to galaxies in the average field. We find that protocluster galaxies, including LAEs, begin to quench much earlier (z∼0.8−1.6) than field galaxies (z∼0.5−0.9); our result is in qualitative agreement with recent observational results and highlights the importance of large-scale environment on the overall formation history of galaxies. 

ABSTRACT The association of GRB170817A with a binary neutron star (BNS) merger has revealed that BNSs produce at least a fraction of short gamma-ray bursts (SGRBs). As gravitational wave (GW) detectors push their horizons, it is important to assess coupled electromagnetic (EM)/GW probabilities and maximize observational prospects. Here, we perform BNS population synthesis calculations with the code mobse, seeding the binaries in galaxies at three representative redshifts, $$z$$ = 0.01, 0.1, and 1 of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the microphysics parameters of GRB170817A, we numerically compute the broad-band light curves of jets from BNS mergers, with the afterglow brightness dependent on the local medium density at the merger site. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find a gamma-ray detection probability of $$\sim\!2{{\rm per\ cent}},10{{\rm per\ cent}},\mathrm{and}\ 40{{\rm per\ cent}}$$ for BNSs at $$z$$ = 1, 0.1, and 0.01, respectively, for the random case, rising to $$\sim\!75{{\rm per\ cent}}$$ for the $$z$$ = 0.01, GW-triggered aligned case. Afterglow detection probabilities of GW-triggered BNS mergers vary in the range of $$\sim \! 0.3 \!-\! 0.5{{\rm per\ cent}}$$, with higher values for aligned jets, and are comparable across the high- and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses. 

ABSTRACT We use the improved IllustrisTNG300 magnetohydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 h−1 Mpc and significant improvements on the subgrid model with respect to previous Illustris simulations, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over a large cosmological volume and halo mass range. We show at high statistical significance that the halo assembly bias signal (i.e. the secondary dependence of halo bias on halo formation redshift) manifests itself on the clustering of the galaxy population when this is split by stellar mass, colour, specific star formation rate, and surface density. A significant signal is also found for galaxy size: at fixed halo mass, larger galaxies are more tightly clustered than smaller galaxies. This effect, in contrast to the rest of the dependencies, seems to be uncorrelated with halo formation time, with some small correlation only detected for halo spin. We also explore the transmission of the spin bias signal, i.e. the secondary dependence of halo bias on halo spin. Although galaxy spin retains little information about the total halo spin, the correlation is enough to produce a significant galaxy spin bias signal. We discuss possible ways to probe this effect with observations.