Researcher profile

Arka Banerjee

Arka Banerjee contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Sahyadri: A simulation suite for the cosmology dependence of the Cosmic Web

We present Sahyadri, a suite of cosmological $N$-body simulations designed to enable precision studies of the low-redshift Universe with next-generation spectroscopic surveys. Sahyadri includes systematic variations of six cosmological parameters around Planck 2018 constraints, with seed-matched initial conditions enabling cosmological parameter derivatives. Each simulation evolves $2048^3$ particles in a periodic box of side length $200$ $h^{-1}$ Mpc, yielding a particle mass of $m_{\rm{p}} = 8.1 \times 10^{7}\,h^{-1}\,M_{\odot}$ in the fiducial Planck 2018 cosmology. This resolution enables robust identification of dark matter halos down to $M_{\rm min} = 3.2 \times 10^{9}$ $h^{-1}$ $M_\odot$, which represents a factor of $\sim$25 improvement over the AbacusSummit suite, and is over two orders of magnitude better than the Quijote and Aemulus suites. We estimate that approximately 40% of DESI BGS galaxies at redshift $z < 0.15$ - roughly 1.6 million objects - reside in halos accessible to Sahyadri but beyond the reach of existing parameter-varying simulation suites. We demonstrate Sahyadri&#39;s capabilities through measurements of the matter power spectrum, halo mass function and power spectrum, and beyond 2-point statistics such as the Voronoi volume function and $k^{\rm th}$ nearest neighbour statistics, showing excellent agreement with theoretical predictions and significant sensitivity to $Ω_{\rm m}$ variations. We implement a custom compression scheme reducing storage requirements by a factor of $\sim$3 while maintaining sub-percent clustering accuracy. Key data products will be made publicly available.

preprint2024arXiv

On the computation of coarse cohomology

The purpose of this article is to relate coarse cohomology of metric spaces with a more computable cohomology. We introduce a notion of boundedly supported cohomology and prove that coarse cohomology of many spaces are isomorphic to the boundedly supported cohomology. Boundedly supported cohomology coincides with compactly supported Alexander--Spanier cohomology if the space is proper and contractible. Our work generalizes an earlier result of Roe which says that the coarse cohomology is isomorphic to the compactly supported Alexander-Spanier cohomology if the space is uniformly contractible. As an application of our main theorem, we obtain that coarse cohomology of the complement can be computed in terms of Alexander-Spanier cohomology for many spaces.

preprint2022arXiv

Beware of Fake $ν$s: The Effect of Massive Neutrinos on the Nonlinear Evolution of Cosmic Structure

Massive neutrinos suppress the growth of cosmic structure on small, non-linear, scales. It is thus often proposed that using statistics beyond the power spectrum can tighten constraints on the neutrino mass by extracting additional information from these non-linear scales. We study the information content regarding neutrino mass at the field level, quantifying how much of this information arises from the difference in non-linear evolution between a cosmology with 1 fluid (CDM) and 2 fluids (CDM + neutrinos). We do so by running two $N$-body simulations, one with and one without massive neutrinos; both with the same phases, and matching their linear power spectrum at a given, low, redshift. This effectively isolates the information encoded in the linear initial conditions from the non-linear cosmic evolution. We demonstrate that for $k \lesssim 1\,h/{\rm Mpc}$, and for a single redshift, there is negligible difference in the real-space CDM field between the two simulations. This suggests that all the information regarding neutrino mass is in the linear power spectrum set by the initial conditions. Thus any probe based on the CDM field alone will have negligible constraining power beyond that which exists at the linear level over the same range of scales. Consequently, any probe based on the halo field will contain little information beyond the linear power. We find similar results for the matter field responsible for weak lensing. We also demonstrate that there may be much information beyond the power spectrum in the 3d matter field, however, this is not observable in modern surveys via dark matter halos or weak lensing. Finally, we show that there is additional information to be found in redshift space.

preprint2022arXiv

Detection of spatial clustering in the 1000 richest SDSS DR8 redMaPPer clusters with Nearest Neighbor distributions

Distances to the $k$-nearest-neighbor ($k$NN) data points from volume-filling query points are a sensitive probe of spatial clustering. Here we present the first application of $k$NN summary statistics to observational clustering measurement, using the 1000 richest redMaPPer clusters ($0.1\leqslant z\leqslant 0.3$) from the SDSS DR8 catalog. A clustering signal is defined as a difference in the cumulative distribution functions (CDFs) of $k$NN distances from fixed query points to the observed clusters versus a set of unclustered random points. We find that the $k=1,2$-NN CDFs of redMaPPer deviate significantly from the randoms&#39; across scales of 35 to 155 Mpc, which is a robust signature of clustering. In addition to $k$NN, we also measure the two-point correlation function for the same set of redMaPPer clusters versus random points, which shows a noisier and less significant clustering signal within the same radial scales. Quantitatively, the $χ^2$ distribution for both the $k$NN-CDFs and the two-point correlation function measured on the randoms peak at $χ^2\sim 50$ (null hypothesis), whereas the $k$NN-CDFs ($χ^2\sim 300$, $p = 1.54\times 10^{-36}$) pick up a much more significant clustering signal than the two-point function ($χ^2\sim 100$, $p = 1.16\times 10^{-6}$) when measured on redMaPPer. Finally, the measured 3NN and 4NN CDFs deviate from the predicted $k=3, 4$-NN CDFs assuming an ideal Gaussian field, indicating a non-Gaussian clustering signal for redMaPPer clusters, although its origin might not be cosmological due to observational systematics. Therefore, $k$NN serves as a more sensitive probe of clustering complementary to the two point correlation function, providing a novel approach for constraining cosmology and galaxy-halo connection.

preprint2022arXiv

Modeling Nearest Neighbor distributions of biased tracers using Hybrid Effective Field Theory

We investigate the application of Hybrid Effective Field Theory (HEFT) -- which combines a Lagrangian bias expansion with subsequent particle dynamics from $N$-body simulations -- to the modeling of $k$-Nearest Neighbor Cumulative Distribution Functions ($k{\rm NN}$-${\rm CDF}$s) of biased tracers of the cosmological matter field. The $k{\rm NN}$-${\rm CDF}$s are sensitive to all higher order connected $N$-point functions in the data, but are computationally cheap to compute. We develop the formalism to predict the $k{\rm NN}$-${\rm CDF}$s of discrete tracers of a continuous field from the statistics of the continuous field itself. Using this formalism, we demonstrate how $k{\rm NN}$-${\rm CDF}$ statistics of a set of biased tracers, such as halos or galaxies, of the cosmological matter field can be modeled given a set of low-redshift HEFT component fields and bias parameter values. These are the same ingredients needed to predict the two-point clustering. For a specific sample of halos, we show that both the two-point clustering \textit{and} the $k{\rm NN}$-${\rm CDF}$s can be well-fit on quasi-linear scales ($\gtrsim 20 h^{-1}{\rm Mpc}$) by the second-order HEFT formalism with the \textit{same values} of the bias parameters, implying that joint modeling of the two is possible. Finally, using a Fisher matrix analysis, we show that including $k{\rm NN}$-${\rm CDF}$ measurements over the range of allowed scales in the HEFT framework can improve the constraints on $σ_8$ by roughly a factor of $3$, compared to the case where only two-point measurements are considered. Combining the statistical power of $k{\rm NN}$ measurements with the modeling power of HEFT, therefore, represents an exciting prospect for extracting greater information from small-scale cosmological clustering.

preprint2022arXiv

Signatures of Light Massive Relics on nonlinear structure formation

Cosmologies with Light Massive Relics (LiMRs) as a subdominant component of the dark sector are well-motivated from a particle physics perspective, and can also have implications for the $σ_8$ tension between early and late time probes of matter clustering. The effects of LiMRs on the Cosmic Microwave Background (CMB) and structure formation on large (linear) scales have been investigated extensively. In this paper, we initiate a systematic study of the effects of LiMRs on smaller, nonlinear scales using cosmological $N$-body simulations; focusing on quantities relevant for photometric galaxy surveys. For most of our study, we use a particular model of nonthermal LiMRs but the methods developed easily generalize to a large class of models of LiMRs -- we explicitly demonstrate this by considering the Dodelson-Widrow form of the velocity distribution. We find that, in general, the effects of LiMR on small scales are distinct from those of a $Λ$CDM universe, even when the value of $σ_8$ is matched between the models. We show that weak lensing measurements around massive clusters, between $\sim 0.1 h^{-1}$Mpc and $\sim 10 h^{-1}$Mpc, should have sufficient signal-to-noise in future surveys to distinguish between $Λ$CDM and LiMR models that are tuned to fit both CMB data and large (linear) scale structure data at late times. Furthermore, we find that different LiMR cosmologies which are indistinguishable by conventional linear probes can be distinguished by these probes if their velocity distributions are sufficiently different. LiMR models can, therefore, be best tested and constrained by jointly analyzing data from CMB and late-time structure formation on both large \textit{and} small scales.

preprint2022arXiv

Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the &#34;Cosmic Frontier&#34;, which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Cosmological Simulations for Dark Matter Physics

Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CMB-S4. The time is now to build a novel simulation program to interpret observations so that we can identify novel signatures of DM microphysics across a large dynamic range of length scales and cosmic time. This white paper identifies the key elements that are needed for such a simulation program. We identify areas of growth on both the particle theory side as well as the simulation algorithm and implementation side, so that we can robustly simulate the cosmic evolution of DM for well-motivated models. We recommend that simulations include a fully calibrated and well-tested treatment of baryonic physics, and that outputs should connect with observations in the space of observables. We identify the tools and methods currently available to make predictions and the path forward for building more of these tools. A strong cosmic DM simulation program is key to translating cosmological observations to robust constraints on DM fundamental physics, and provides a connection to lab-based probes of DM physics.

preprint2022arXiv

Tidal disruption of solitons in self-interacting ultralight axion dark matter

Ultralight axions (ULAs) are promising dark matter candidates that can have a distinct impact on the formation and evolution of structure on nonlinear scales relative to the cold, collisionless dark matter (CDM) paradigm. However, most studies of structure formation in ULA models do not include the effects of self-interactions, which are expected to arise generically. Here, we study how the tidal evolution of solitons is affected by ULA self-interaction strength and sign. Specifically, using the pseudospectral solver UltraDark.jl, we simulate the tidal disruption of self-interacting solitonic cores as they orbit a $10^{11}~M_{\mathrm{\odot}}$ Navarro-Frenk-White CDM host halo potential for a range of orbital parameters, assuming a fiducial ULA particle mass of $10^{-22}\mathrm{eV}$. We find that repulsive (attractive) self-interactions significantly accelerate (decelerate) soliton tidal disruption. We also identify a degeneracy between the self-interaction strength and soliton mass that determines the efficiency of tidal disruption, such that disruption timescales are affected at the $\sim 50\%$ level for variations in the dimensionless ULA self-coupling from $λ=-10^{-92}$ to $λ=10^{-92}$.

preprint2022arXiv

Tracer-Field Cross-Correlations with $k$-Nearest Neighbor Distributions

In astronomy and cosmology, significant effort is devoted to characterizing and understanding spatial cross-correlations between points - e.g. galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous field - e.g. weak lensing and Cosmic Microwave Background (CMB) maps. Recently, we introduced the $k$-nearest neighbor formalism to better characterize the clustering of discrete (point) datasets. Here we extend it to the point-field cross-correlation analysis. It combines $k$NN measurements of the point data set with measurements of the field smoothed on many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two datasets even when there are no linear (Gaussian) correlations. We further demonstrate that this framework is far more effective than the two-point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlations between halos and the underlying matter field in a cosmological simulation, between $10h^{-1}{\rm Mpc}$ and $30h^{-1}{\rm Mpc}$, is detected at $>5σ$ significance using the technique presented here, when the two-point cross-correlation significance is $\sim 1σ$. Finally, we show that the $k$NN cross-correlations of halos and the matter field can be well-modeled on quasilinear scales by the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for the two-point cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations with this method makes it a promising tool for various cosmological applications.

preprint2021arXiv

A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos

Quasi-N-body simulations, such as FastPM, provide a fast way to simulate cosmological structure formation, but have yet to adequately include the effects of massive neutrinos. We present a method to include neutrino particles in FastPM, enabling computation of the CDM and total matter power spectra to percent-level accuracy in the non-linear regime. The CDM-neutrino cross-power can also be computed at a sufficient accuracy to constrain cosmological observables. To avoid the shot noise that typically plagues neutrino particle simulations, we employ a quasi-random algorithm to sample the relevant Fermi-Dirac distribution when setting the initial neutrino thermal velocities. We additionally develop an effective distribution function to describe a set of non-degenerate neutrinos as a single particle to speed up non-degenerate simulations. The simulation is accurate for the full range of physical interest, $M_ν\lesssim 0.6$eV, and applicable to redshifts $z\lesssim2$. Such accuracy can be achieved by initializing particles with the two-fluid approximation transfer functions (using the REPS package). Convergence can be reached in $\sim 25$ steps, with a starting redshift of $z=99$. Probing progressively smaller scales only requires an increase in the number of CDM particles being simulated, while the number of neutrino particles can remain fixed at a value less than or similar to the number of CDM particles. In turn, the percentage increase in runtime-per-step due to neutrino particles is between $\sim 5-20\%$ for runs with $1024^3$ CDM particles, and decreases as the number of CDM particles is increased. The code has been made publicly available, providing an invaluable resource to produce fast predictions for cosmological surveys and studying reconstruction.

preprint2021arXiv

Nearest Neighbor distributions: new statistical measures for cosmological clustering

The use of summary statistics beyond the two-point correlation function to analyze the non-Gaussian clustering on small scales is an active field of research in cosmology. In this paper, we explore a set of new summary statistics -- the $k$-Nearest Neighbor Cumulative Distribution Functions ($k{\rm NN}$-${\rm CDF}$). This is the empirical cumulative distribution function of distances from a set of volume-filling, Poisson distributed random points to the $k$-nearest data points, and is sensitive to all connected $N$-point correlations in the data. The $k{\rm NN}$-${\rm CDF}$ can be used to measure counts in cell, void probability distributions and higher $N$-point correlation functions, all using the same formalism exploiting fast searches with spatial tree data structures. We demonstrate how it can be computed efficiently from various data sets - both discrete points, and the generalization for continuous fields. We use data from a large suite of $N$-body simulations to explore the sensitivity of this new statistic to various cosmological parameters, compared to the two-point correlation function, while using the same range of scales. We demonstrate that the use of $k{\rm NN}$-${\rm CDF}$ improves the constraints on the cosmological parameters by more than a factor of $2$ when applied to the clustering of dark matter in the range of scales between $10h^{-1}{\rm Mpc}$ and $40h^{-1}{\rm Mpc}$. We also show that relative improvement is even greater when applied on the same scales to the clustering of halos in the simulations at a fixed number density, both in real space, as well as in redshift space. Since the $k{\rm NN}$-${\rm CDF}$ are sensitive to all higher order connected correlation functions in the data, the gains over traditional two-point analyses are expected to grow as progressively smaller scales are included in the analysis of cosmological data.

preprint2020arXiv

A Cross-Correlation Study of High-energy Neutrinos and Tracers of Large-Scale Structure

The origin of the bulk of the astrophysical neutrinos detected by the IceCube Observatory remains a mystery. Previous source-finding analyses compare the directions of IceCube events and individual sources in astrophysical catalogs. The source association method is technically challenging when the number of source candidates is much larger than the number of the observed astrophysical neutrinos. We show that in this large source number regime, a cross-correlation analysis of neutrino data and source catalog can instead be used to constrain potential source populations for the high-energy astrophysical neutrinos, and provide spatial evidence for the existence of astrophysical neutrinos. We present an analysis of the cross-correlation of the IceCube 2010-2012 point-source data and a WISE-2MASS galaxy sample. While we find no significant detection of cross-correlation with the publicly available neutrino dataset, we show that, when applied to the full IceCube data, which has a longer observation time and higher astrophysical neutrino purity, our method has sufficient statistical power to detect a cross-correlation signal if the neutrino sources trace the Large Scale Structure of the Universe.

preprint2020arXiv

Investigating the use of field solvers for simulating classical systems

We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the non linear regime. We also investigate extensions of the Schrödinger-Poisson system that employ multiple stacked cold streams, and the von Neumann-Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. We present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.

preprint2020arXiv

Signatures of Velocity-Dependent Dark Matter Self-Interactions in Milky Way-mass Halos

We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, $w$, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of $w$ that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of $w$. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of $\sim 200\ \rm{km\ s}^{-1}$, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

preprint2019arXiv

Fisher for complements: Extracting cosmology and neutrino mass from the counts-in-cells PDF

We comprehensively analyse the cosmology dependence of counts-in-cell statistics. We focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly nonlinear scales. Based on large-deviation statistics, we parametrise the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the nonlinear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass $M_ν$ and can disentangle it from the clustering amplitude $σ_8$. Using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the $νΛ$CDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first $νΛ$CDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly nonlinear scales. The matter PDF is particularly powerful for constraining the matter density $Ω_m$, clustering amplitude $σ_8$ and the total neutrino mass $M_ν$. Adding the mildly nonlinear matter PDF to the mildly nonlinear matter power spectrum improves constraints on $Ω_m$ by a factor of 5 and $σ_8$ by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc/h)$^3$. We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

preprint2019arXiv

Signatures of Self-Interacting dark matter on cluster density profile and subhalo distributions

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($σ/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.

preprint2019arXiv

Weighing neutrinos with the halo environment

Nonlinear objects like halos and voids exhibit a scale-dependent bias on linear scales in massive neutrino cosmologies. The shape of this scale-dependent bias is a unique signature of the neutrino masses, but the amplitude of the signal is generally small, of the order of $f_ν$, the contribution of neutrinos to the total matter content ($\lesssim 1\%$). In this paper, we demonstrate for the first time how the strength of this signal can be substantially enhanced by using information about the halo environment at a range of scales. This enhancement is achieved by using certain combinations of the large scale Cold Dark Matter and total matter environments of halos, both of which are measurable from galaxy clustering and weak lensing surveys.