SBSownak Bose

Sownak Bose

*This is an algorithmically generated author record based on publicly available information.

Current Organizations / Institutes

The organizations/institutes from the scholar's latest article


Sownak Bose has a total of 884 co-authors. Here are the top co-authors by number of shared publications:


Latest Publications

  • Open Access

    Halo mass functions in mixed cold and fuzzy dark matter models

    by Sarah C Johnston, Simon May, Tibor Dome, Sownak Bose, Alastair Basden, Carlton Baugh, Anastasia Fialkov, Alex Tocher

    We investigate the impact of mixed cold and fuzzy dark matter (referred to as MDM) cosmologies on the halo mass function (HMF) using numerical simulations performed with the AxiREPO framework. We consider models in which an ultralight axion-like component with mass $m =10^{-24.5}\, \mathrm{eV}$ constitutes a fraction $f \le 0.3$ of the total dark matter. To enable consistent halo identification in mixed-species scenarios, we develop a grid-based halo-finding pipeline that combines the particle-based cold dark matter (CDM) and wave-like fuzzy dark matter (FDM) components into a unified density field. We find that FDM traces the large-scale CDM distribution while suppressing small-scale structure through wave interference effects, leading to a reduction in the abundance of low-mass haloes and modifying the HMF in a manner dependent on redshift and FDM fraction. Increasing the FDM fraction produces a systematic downward shift in the HMF and modifies its high-mass slope. Motivated by these trends, we introduce a phenomenological model that maps CDM HMFs to their MDM counterparts using a suppression function with parameters dependent on redshift and FDM fraction. This model reproduces the simulated HMFs within approximately 0.1–0.2 dex across the parameter space explored ($1 \le z \le 4$, $f \le 0.3$). Our results provide a computationally efficient method for predicting structure formation in MDM cosmologies without requiring dedicated simulations for each parameter choice, and establish a framework for exploring the impact of MDM on cosmological structure formation.

  • Open Access

    The galaxy ultraviolet luminosity function from z = 7 to 20 in the COLIBRE simulations

    by Shengdong Lu, Carlos S Frenk, Cedric G Lacey, Andrea Gebek, Joop Schaye, Shaun Cole, Sownak Bose, Anna Durrant, Nick Andreadis, Maarten Baes, et al.

    James Webb Space Telescope (JWST) has enabled the detection of galaxies in the earliest stages of cosmic history. We compare the ultraviolet luminosity functions (UVLFs) at redshifts $z=7\mathrm{-}20$ predicted by a set of new cosmological hydrodynamical simulations, colibre (COLd ISM and Better REsolution), with observations, including those from JWST. The UV luminosities of colibre galaxies are derived using the radiative transfer code skirt, which tracks stellar emission and its processing through the multiphase interstellar medium and dust distribution predicted by colibre. We find that although colibre is consistent with the observed evolution of the stellar mass function up to $z=12$, its dust-attenuated UVLFs fall systematically below the observations at the bright end: at the number density of $10^{-6}\, \mathrm{Mpc^{-3}\, mag^{-1}}$, the brightest galaxies are underluminous by $\approx 1\, \rm mag$ at $z=7$, increasing to $\approx 2.5\, \rm mag$ at $z=15$. Accounting for observational uncertainties brings the colibre UVLFs closer to the observational data, but does not fully resolve the discrepancy. Ignoring dust attenuation allows colibre to produce sufficiently bright galaxies at $7\lesssim z \lesssim 12$, but at $z=15$, colibre still underpredicts the luminosities of the brightest galaxies, indicating the need for additional physical mechanisms to boost the UV luminosities at the earliest cosmic epochs, such as a ‘top-heavy’ stellar initial mass function. We fit the colibre UVLFs with Schechter functions and calculate the evolution of the best-fitting parameters. We find that the galaxy number density decreases, the characteristic luminosity becomes fainter and the faint-end slope becomes steeper towards higher redshifts. The UV luminosity density decreases by a factor of several hundred from $z = 7$ to $15$.


  • Open Access

    A modern halo streaming model for redshift space distortions

    by Cheng-Zong Ruan, Baojiu Li, Carlton M Baugh, Sownak Bose, Alexander Eggemeier, David F Mota

    Accurate modelling of redshift-space distortions (RSD) in galaxy clustering is essential for extracting cosmological information from current and forthcoming large-scale structure surveys. While perturbation theory is reliable on large scales, much of the constraining power lies at intermediate and small separations, where non-linear dynamics within and between dark matter haloes dominate. We present a halo streaming model for non-linear galaxy clustering in redshift space that is accurate and physically interpretable. Our framework combines the streaming model for RSD with a halo-model decomposition of the galaxy clustering into central/satellite and one-/two-halo contributions. We build dedicated emulators for the key physical ingredients, trained on a suite of N-body simulations: halo mass functions, real-space halo two-point correlation functions, and pairwise velocity moments. By emulating these modular building blocks rather than the final redshift-space observable, this approach preserves physical transparency, enables targeted optimization for each ingredient, and remains flexible to changes in tracer populations and galaxy–halo connection models. The resulting halo streaming model reproduces the simulated non-linear anisotropic clustering signal down to highly non-linear scales, while achieving the computational efficiency required for cosmological parameter inference. This framework is designed to support full-shape RSD analyses for surveys such as DESI and Euclid, facilitating precision measurements of structure growth and tests of gravity. All codes and trained emulators are publicly available in the freyja repository (https://github.com/chzruan/freyja).

  • Open Access

    The BUFFALO survey: the subhalo mass function for Abell 2744 with strong + weak gravitational lensing

    by Nency R Patel, Mathilde Jauzac, Anna Niemiec, Guillaume Mahler, David Lagattuta, David Harvey, Johannes Schwinn, Andrew Robertson, Lukas J Furtak, Sownak Bose, et al.

    We present a combined strong and weak gravitational lensing mass model of the massive merging galaxy cluster Abell 2744 ($z = 0.308$). The mass reconstruction combines strong-lensing constraints in the cluster core with weak-lensing measurements derived from deep imaging with the Hubble Space Telescope and the James Webb Space Telescope, enabling a reconstruction of the projected mass distribution from the core to the cluster outskirts. Using the resulting mass map, we identify cluster substructures and construct the subhalo mass function for Abell 2744. To interpret the observed substructure population, we apply identical detection methods to projected mass maps of massive galaxy clusters from the BAryons and HAloes of MAssive Systems (BAHAMAS) simulations. We observe a systematic excess of detected substructures in Abell 2744 compared to BAHAMAS. We interpret this excess as likely reflecting Abell 2744’s exceptionally disturbed dynamical state, characterized by complex ongoing mergers, rather than tension with Lambda cold dark matter or self-interacting dark matter predictions. However, confirming this interpretation will require extending the analysis to a larger sample of clusters spanning a range of dynamical states.

  • Open Access

    Galaxy luminosity functions from far-UV to submillimetre at z = 0 in the COLIBRE simulations

    by Shengdong Lu, Carlos S Frenk, Cedric G Lacey, Andrea Gebek, Joop Schaye, Shaun Cole, Sownak Bose, Nick Andreadis, Maarten Baes, Alejandro Benítez-Llambay, et al.

    We present predictions from the colibre cosmological hydrodynamical simulations of galaxy formation for the present-day galaxy luminosity functions (LFs) at wavelengths ranging from the far-ultraviolet (FUV) to the submillimetre. The simulations are post-processed with the radiative transfer code skirt, accounting for dust attenuation and emission using the distribution and properties of dust grains predicted directly by colibre. Results from simulations varying in mass resolution by a factor of $\sim 10^2$ ($\sim 10^5\text {-} 10^7\, \mathrm{{\rm M}_{\odot }}$) show very good convergence over most luminosity ranges. The colibre-skirt LFs match the data remarkably well from the FUV to $3.4\, \rm{\mu m}$ and also in the far-infrared and submillimetre wavelength range ($70\text {-}850\, \rm{\mu m}$). In the mid-infrared (MIR; $8\text {-}24\, \rm{\mu m}$), colibre-skirt matches the data well at low luminosities but significantly underpredicts the luminosities of MIR-bright galaxies, with the discrepancy increasing towards longer wavelengths. The total infrared LF, obtained by integrating the spectral energy distributions over $8\text {-}1000\, \rm{\mu m}$, also matches observations well at the faint end but underpredicts the number of very bright galaxies. The unprecedented agreement at all other wavelengths indicates that colibre, coupled with this calibration-free skirt post-processing framework, successfully predicts the properties of stellar populations at the present day and the amount and distribution of interstellar dust.

  • Open Access

    Modelling the expulsion of baryons from haloes: the role of feedback and of the cosmological constant

    by Oscar Veenema, Daniele Sorini, Sownak Bose

    The extent to which galactic-scale astrophysical processes conspire with the underlying cosmological model to expel baryons from haloes remains a central question in galaxy formation. We present an analytical model for the gas distribution within and beyond haloes, based on the balance between gravitational collapse, hydrostatic pressure, and cosmic expansion. Our model predicts the halo-centric distance enclosing a baryon mass fraction equal to the cosmic value $f_{\rm b} = \Omega _{\rm b}/\Omega _{\rm m}$ (‘closure radius’) in an arbitrary $\Lambda$CDM cosmology. We compare the predictions with the results of six variants of the EAGLE cosmological, hydrodynamical simulation, encompassing values of the cosmological constant ranging from 0 to 100 times its observed value in our Universe, $\Lambda _0$. Despite its simplicity, our model exhibits excellent agreement with the simulations for haloes with mass $M_{\rm 200c} \gt 10^{11} \, \mathrm{M}_{\odot }$ in the redshift range $0\lt z\lt 3$, suggesting that it captures the key astrophysical processes and highlighting its robustness to the cosmological parameters. Thus, it provides the first physical explanation for the empirical closure radius–halo mass relation previously observed in simulations. Furthermore, we find that dark energy plays a non-negligible role in baryon evacuation: the simulations reveal that in the fiducial cosmological model, the closure radius at $z\lt 2$ is $\sim 30 {{\ \rm per\ cent}}$ larger than in an Einstein–de Sitter universe. In cosmologies with $\Lambda \ge 10 \Lambda _0$, dark energy emerges as the dominant factor in this process – suggesting that, as our Universe transitions towards $\Lambda$-domination, dark energy eventually becomes the primary driver of baryon evacuation from massive haloes.

  • Open Access

    Secondary dependence of baryonic effects on the density profile of dark matter halos

    by Yikun Wang, Idit Zehavi, Sergio Contreras, Giovanni Aricò, Sownak Bose, Lars Hernquist

    Baryonic physics is anticipated to be a major source of systematic uncertainty in current and future large-scale cosmological surveys. We investigate how baryonic effects on halo density profiles vary with secondary halo properties at fixed halo mass, using the large-volume MillenniumTNG hydrodynamical simulation and its dark matter-only counterpart. We focus on the impact of halo concentration and large-scale environment on the ratio of density profiles of matched halos in the hydrodynamical and dark matter-only simulations. At redshift z = 0.0, we find a strong dependence on halo concentration, especially at lower halo mass (12.5 h / h -1 M ⊙ ) < 13.0), where more concentrated halos exhibit weaker inner enhancement and stronger intermediate-radius suppression at fixed halo mass, with variations reaching ∼ 15% at small scales and decreasing toward larger scales. This trend weakens and reverses at higher halo mass. In contrast, the secondary dependence on large-scale environment is weaker (∼ 2%) and largely scale-independent, with halos in denser regions exhibiting slightly weaker intermediate suppression. By separating internal profile redistribution from total mass suppression, we show that concentration impacts both components, whereas the environmental dependence is primarily associated with an overall mass shift. These secondary dependencies persist at z = 0.5 and correlate with variations in internal baryonic properties. We examine additional halo properties, including halo spin and velocity dispersion, and find significant secondary dependence. Overall, our results highlight the important role of secondary halo properties in modulating baryonic effects on halo density profiles, with potential implications for future modeling efforts.

  • Top Cited Papers
    Open Access

    Little red dots host black hole stars: a unified family of gas-reddened AGN revealed by JWST /NIRSpec spectroscopy

    by Anna de Graaff, Raphael E Hviding, Rohan P Naidu, Jenny E Greene, Tim B Miller, Joel Leja, Jorryt Matthee, Gabriel Brammer, Harley Katz, Rachel Bezanson, et al.

    We construct and characterize a sample of 146 little red dots (LRDs) across ${2.0{\lt} z{\lt} 9.3}$, selecting all sources with v-shaped ultraviolet (UV)–optical continua from James Webb Space Telescope (JWST) NIRSpec/PRISM spectra and compact morphologies in NIRCam/F444W imaging. We show that LRD continuum spectra are well described by modified blackbodies across $\sim 0.4-1.0\, \rm{\mu m}$, with typical $T\sim 5000\,$ K ($\lambda _{\rm peak}\sim 0.65\, \rm{\mu m}$) across 2 dex in luminosity, and a tail toward $T\sim 2000\,$ K. LRDs therefore trace a locus in the Hertzsprung–Russell diagram that is analogous to stars on the Hayashi track, strongly supporting the picture that LRDs are active galactic nucleus (AGN) embedded in optically thick dense gas envelopes. Hotter LRDs with $\lambda _{\rm peak}\lt 0.65\, \rm{\mu m}$ typically have strong Balmer breaks, redder UV slopes, and high optical luminosities; other LRDs show weak or no Balmer breaks, and wide variety in $\beta _{\rm UV}$ and $L_{5100}$. Crucially, we demonstrate that the UV–optical continua are strongly linked to the H$\alpha$, H$\beta$, [O iii], and O i line properties. There is a tight linear relation between the H$\alpha$ and optical continuum luminosities, as well as H$\alpha$ and $\rm {O\, {\small I} }_{\ \lambda 8446}$, indicating that Balmer, O i, and optical emission must primarily be powered by the same source. The Balmer decrement increases strongly toward higher $L_{\rm H\alpha }$, $L_{5100}$, and Balmer break strength, providing key evidence for luminosity-dependent effects of collisional (de-)excitation and resonant scattering in the gaseous envelopes. In contrast, we show that [O iii] emission likely originates from star-forming host galaxies, with strong variation in the AGN-to-host ratio among the LRD population. Our work presents an empirical description of the nature and structure of LRDs, defining a new benchmark for ongoing LRD model developments.


  • Open Access

    Galaxy formation in modified gravity – II. Galaxy halo connection and assembly bias

    by Michael Collier, Sownak Bose, Baojiu Li

    Modern surveys such as Dark Energy Spectroscopic Instrument and Euclid, which collect data for hundreds of millions of galaxies to map the large-scale structure (LSS) of the Universe, hold the key to determining the cosmological parameters and testing new physics. This ambition, however, is limited by uncertainties in the galaxy–halo connection: the link between observed galaxies and the underlying, unobservable matter field, by accounting for effects such as galaxy bias and assembly bias (AB). These are particularly poorly understood for modified gravity (MG) models, which are popular alternatives to the cosmological constant to explain accelerated expansion. We approach this problem using mock emission-line galaxy (ELG) and luminous red galaxy (LRG) catalogues in $f(R)$ gravity matching the specifications of ongoing Stage-IV galaxy surveys, generated from state-of-the-art MG hydrodynamical simulations. While the interplay between MG – especially the chameleon screening mechanism – and galaxy formation leaves complicated imprints in the galaxy–halo connection, a simple physical picture emerges in which halo and galaxy formation are enhanced for progressively more massive haloes over time. We confirm that the basic galaxy–halo connection model, the halo occupation distribution (HOD), in which galaxy occupation is determined solely by halo mass, underestimates galaxy clustering strength in Lmbda-cold dark matter by 10–20 per cent at $z\lesssim 1$ when neglecting AB, and demonstrate that MG introduces further complexity. Extending this model with a suitably chosen environment density as a secondary HOD variable reduces the AB effect in all models to 2–3 per cent for $z\lesssim 0.5$. This provides a well-motivated starting point for further works on minimizing the impact of AB when testing non-standard cosmological models with LSS.