BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Shi-Fan Chen

Shi-Fan Chen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.COastro-ph.IMcond-mat.stat-mechhep-th

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
4topics
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

12 published item(s)

preprint2026arXiv

Effective Theories of Redshift-Space Galaxy Peculiar Velocities

We present predictions for redshift-space peculiar velocity statistics in the Lagrangian and Eulerian formulations of the effective field theory (EFT) of large-scale structure. We compute 2-point pairwise velocity statistics up to the second moment at next-to-leading (1-loop) order, showing that they can be modeled together with redshift-space galaxy densities with a consistent set of EFT coefficients. We show that peculiar velocity statistics have a distinct dependence on long-wavelength bulk flows that necessitates a variation on the usual infrared (IR) resummation procedure used to model baryon acoustic oscillations (BAO) in galaxy clustering. This can be implemented recursively in powers of the velocity in both the Lagrangian and Eulerian frameworks. We validate our analytic calculations against fully nonlinear N-body simulations, demonstrating that they can be used to recover the growth rate at better than percent level precision, well beyond the statistical requirements of upcoming peculiar velocity surveys and measurements of the kinetic Sunyaev-Zeldovich (kSZ) effect. As part of this work, we release $\href{https://github.com/sfschen/velocisaurus}{\texttt{velocisaurus}}$, a fast $\texttt{Python}$ code for computing EFT predictions of peculiar velocity statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

A Lagrangian theory for galaxy shape statistics

We formulate the Lagrangian perturbation theory of galaxy intrinsic alignments and compute the resulting auto and cross power spectra of galaxy shapes, densities and matter to 1-loop order. Our model represents a consistent effective-theory description of galaxy shape including the resummation of long-wavelength displacements which damp baryon acoustic oscillations, and includes one linear, three quadratic and two cubic dimensionless bias coefficients at this order, along with counterterms and stochastic contributions whose structure we derive. We compare this Lagrangian model against the three-dimensional helicity spectra of halo shapes measured in N-body simulations by Akitsu et al (2023) and find excellent agreement on perturbative scales while testing a number of more restrictive bias parametrizations. The calculations presented are immediately relevant to analyses of both cosmic shear surveys and spectroscopic shape measurements, and we make a fast FFTLog-based code spinosaurus publicly available with this publication.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A new analysis of galaxy 2-point functions in the BOSS survey, including full-shape information and post-reconstruction BAO

We present a new method for consistent, joint analysis of the pre- and post-reconstruction two-point functions of the BOSS survey. The post-reconstruction correlation function is used to accurately measure the distance-redshift relation and expansion history, while the pre-reconstruction power spectrum multipoles constrain the broad-band shape and the rate-of-growth of large-scale structure. Our technique uses Lagrangian perturbation theory to self-consistently work at the level of two-point functions, i.e.\ directly with the measured data, without approximating the constraints with summary statistics normalized by the drag scale. Combining galaxies across the full redshift range and both hemispheres we constrain $Ω_m=0.303 \pm 0.0082$, $H_0=69.23 \pm 0.77$ and $σ_8=0.733 \pm 0.047$ within the context of $Λ$CDM. These constraints are in good agreement both with the Planck primary CMB anisotropy data and recent cosmic shear surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmological Analysis of Three-Dimensional BOSS Galaxy Clustering and Planck CMB Lensing Cross Correlations via Lagrangian Perturbation Theory

We present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain $Ω_\mathrm{m} = 0.303\pm 0.008$, $H_0 = 69.21\pm 0.78$ and $σ_8 = 0.743\pm 0.043$. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to $Ω_m = 0.300 \pm 0.008$, $H_0 = 69.21 \pm 0.77$ and $σ_8 = 0.707 \pm 0.035$, in tension with CMB and cosmic shear constraints. The combination of $Ω_m$ and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $σ_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($z\simeq 0.38$) and at large angular scales, which show a $20\%$ deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmological constraints from the tomographic cross-correlation of DESI Luminous Red Galaxies and Planck CMB lensing

We use luminous red galaxies selected from the imaging surveys that are being used for targeting by the Dark Energy Spectroscopic Instrument (DESI) in combination with CMB lensing maps from the Planck collaboration to probe the amplitude of large-scale structure over $0.4\le z\le 1$. Our galaxy sample, with an angular number density of approximately $500\,\mathrm{deg}^{-2}$ over 18,000 sq.deg., is divided into 4 tomographic bins by photometric redshift and the redshift distributions are calibrated using spectroscopy from DESI. We fit the galaxy autospectra and galaxy-convergence cross-spectra using models based on cosmological perturbation theory, restricting to large scales that are expected to be well described by such models. Within the context of $Λ$CDM, combining all 4 samples and using priors on the background cosmology from supernova and baryon acoustic oscillation measurements, we find $S_8=σ_8(Ω_m/0.3)^{0.5}=0.73\pm 0.03$. This result is lower than the prediction of the $Λ$CDM model conditioned on the Planck data. Our data prefer a slower growth of structure at low redshift than the model predictions, though at only modest significance.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Engineered swift equilibration for arbitrary geometries

Engineered swift equilibration (ESE) is a class of driving protocols that enforce an equilibrium distribution with respect to external control parameters at the beginning and end of rapid state transformations of open, classical non-equilibrium systems. ESE protocols have previously been derived and experimentally realized for Brownian particles in simple, one-dimensional, time-varying trapping potentials; one recent study considered ESE in two-dimensional Euclidean configuration space. Here we extend the ESE framework to generic, overdamped Brownian systems in arbitrary curved configuration space and illustrate our results with specific examples not amenable to previous techniques. Our approach may be used to impose the necessary dynamics to control the full temporal configurational distribution in a wide variety of experimentally realizable settings.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Priors on red galaxy stochasticity from hybrid effective field theory

We investigate the stochastic properties of typical red galaxy samples in a controlled numerical environment. We use Halo Occupation Distribution (HOD) modelling to create mock realizations of three separate bright red galaxy samples consistent with datasets used for clustering and lensing analyses in modern galaxy surveys. Second-order Hybrid Effective Field Theory (HEFT) is used as a field-level forward model to describe the full statistical distribution of these tracer samples, and their stochastic power spectra are directly measured and compared to the Poisson shot-noise prediction. While all of the galaxy samples we consider are hosted within haloes with sub-Poisson stochasticity, we observe that the galaxy samples themselves possess stochasticities that range from sub-Poisson to super-Poisson, in agreement with predictions from the halo model. As an application of our methodology, we place priors on the expected degree of non-Poisson stochasticity in cosmological analyses using such samples. We expect these priors will be useful in reducing the complexity of the full parameter space for future analyses using second-order Lagrangian bias models. More generally, the techniques outlined here present the first application of hybrid EFT methods to characterize models of the galaxy--halo connection at the field level, revealing new connections between once-disparate modelling frameworks.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Neural Network Acceleration of Large-scale Structure Theory Calculations

We make use of neural networks to accelerate the calculation of power spectra required for the analysis of galaxy clustering and weak gravitational lensing data. For modern perturbation theory codes, evaluation time for a single cosmology and redshift can take on the order of two seconds. In combination with the comparable time required to compute linear predictions using a Boltzmann solver, these calculations are the bottleneck for many contemporary large-scale structure analyses. In this work, we construct neural network-based surrogate models for Lagrangian perturbation theory (LPT) predictions of matter power spectra, real and redshift space galaxy power spectra, and galaxy--matter cross power spectra that attain $\sim 0.1\%$ (at one sigma) accuracy over a broad range of scales in a $w$CDM parameter space. The neural network surrogates can be evaluated in approximately one millisecond, a factor of 1000 times faster than the full Boltzmann code and LPT computations. In a simulated full-shape redshift space galaxy power spectrum analysis, we demonstrate that the posteriors obtained using our surrogates are accurate compared to those obtained using the full LPT model. We make our surrogate models public at https://github.com/sfschen/EmulateLSS, so that others may take advantage of the speed gains they provide to enable rapid iteration on analysis settings, something that is essential in complex contemporary large-scale structure analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotic expansions for the Large Scale Structure

We explore the deep ultraviolet (that is, short-distance) limit of the power spectrum (PS) and of the correlation function of a cold dark matter dominated Universe. While for large scales the PS can be written as a double series expansion, in powers of the linear PS and of the wavenumber $k$, we show that, in the opposite limit, it can be expressed via an expansion in powers of the form $1/k^{d+2n}$, where $d$ is the number of spatial dimensions, and $n$ is a non negative integer. The coefficients of the terms of the expansion are nonperturbative in the linear PS, and can be interpreted in terms of the probability density function for the displacement field, evaluated around specific configurations of the latter, that we identify. In the case of the Zel'dovich dynamics, these coefficients can be determined analytically, whereas for the exact dynamics they can be treated as fit, or nuisance, parameters. We confirm our findings with numerical simulations and discuss the necessary steps to match our results to those obtained for larger scales and to actual measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Consistent Modeling of Velocity Statistics and Redshift-Space Distortions in One-Loop Perturbation Theory

The peculiar velocities of biased tracers of the cosmic density field contain important information about the growth of large scale structure and generate anisotropy in the observed clustering of galaxies. Using N-body data, we show that velocity expansions for halo redshift-space power spectra are converged at the percent-level at perturbative scales for most line-of-sight angles $μ$ when the first three pairwise velocity moments are included, and that the third moment is well-approximated by a counterterm-like contribution. We compute these pairwise-velocity statistics in Fourier space using both Eulerian and Lagrangian one-loop perturbation theory using a cubic bias scheme and a complete set of counterterms and stochastic contributions. We compare the models and show that our models fit both real-space velocity statistics and redshift-space power spectra for both halos and a mock sample of galaxies at sub-percent level on perturbative scales using consistent sets of parameters, making them appealing choices for the upcoming era of spectroscopic, peculiar-velocity and kSZ surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simulations and symmetries

We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be efficiently achieved by modern N-body methods with a symmetries-based bias expansion which describes the clustering of any tracer on large scales. We show that at low redshifts, and for moderately biased tracers, the substitution of N-body-determined dynamics improves over an equivalent model using perturbation theory by more than a factor of two in scale, while at high redshifts and for highly biased tracers the gains are more modest. This hybrid approach lends itself well to emulation. By removing the need to identify halos and subhalos, and by not requiring any galaxy-formation-related parameters to be included, the emulation task is significantly simplified at the cost of modeling a more limited range in scale.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Reconstructed Power Spectrum in the Zeldovich Approximation

Density-reconstruction sharpens the baryon acoustic oscillations signal by undoing some of the smoothing incurred by nonlinear structure formation. In this paper we present an analytical model for reconstruction based on the Zeldovich approximation, which for the first time includes a complete set of counterterms and bias terms up to quadratic order and can fit real and redshift-space data pre- and post-reconstruction data in both Fourier and configuration space over a wide range of scales. We compare our model to n-body data at $z = 0$ from the {\tt DarkSky} simulation \cite{skillman14}, finding sub-percent agreement in both real space and in the redshift-space power spectrum monopole out to $k = 0.4\ h$ Mpc$^{-1}$, and out to $k = 0.2\ h$ Mpc$^{-1}$ in the quadrupole, with comparable agreement in configuration space. We compare our model with several popular existing alternatives, updating existing theoretical results for exponential damping in wiggle/no-wiggle splits of the BAO signal and discuss the usually-ignored effect of higher bias contributions on the reconstructed signal. In the appendices, we re-derive the former within our formalism, present exploratory results on higher-order corrections due to nonlinearities inherent to reconstruction, and present numerical techniques with which to calculate the redshift-space power spectrum of biased tracers within the Zeldovich approximation.

0 citations0 reviewsNo code linkedOpen access