Related papers: Speeding up N-body simulations of modified gravity…
N-body simulations are the most powerful method to study the non-linear evolution of large-scale structure. However, they require large amounts of computational resources, making unfeasible their direct adoption in scenarios that require…
Model-independent constraints on modified gravity models hitherto exist mainly on linear scales. A recently developed formalism presented a consistent parameterisation that is valid on all scales. Using this approach, we perform…
We use the Millennium Simulation to quantify the statistical accuracy and precision of the escape velocity technique for measuring cluster-sized halo masses at z~0.1. We show that in 3D, one can measure nearly unbiased (<4%) halo masses…
We consider the extent to which future imaging surveys of galaxies can distinguish between dark energy and modified gravity models for the origin of the cosmic acceleration. Dynamical dark energy models may have similar expansion rates as…
We describe the construction of an N-body simulation that mimics the true velocity and mass-density fields in a box of side 256\hmpc about the Local Group , and the production of mock catalogs that mimic in detail current catalogs of…
High-fidelity physics simulations are powerful tools in the design and optimization of charged particle accelerators. However, the computational burden of these simulations often limits their use in practice for design optimization and…
The effects of discreteness arising from the use of the N-body method on the accuracy of simulations of cosmological structure formation are not currently well understood. After a discussion of how the relevant discretisation parameters…
Within the landscape of modified theories of gravity, progress in understanding the behaviour of, and developing tests for, screening mechanisms has been hindered by the complexity of the field equations involved, which are nonlinear in…
Even though Nearest Neighbor Gaussian Processes (NNGP) alleviate considerably MCMC implementation of Bayesian space-time models, they do not solve the convergence problems caused by high model dimension. Frugal alternatives such as response…
We generalize a recently developed method for accelerated Monte Carlo calculation of path integrals to the physically relevant case of generic many-body systems. This is done by developing an analytic procedure for constructing a hierarchy…
(Abridged) We describe a new iterative approach for the realization of equilibrium N-body systems for given density distributions. Our method uses elements of Schwarzschild's technique and of the made-to-measure method, but is based on a…
In a mean-field approximation within Kinetic Field Theory, it is possible to derive an accurate analytic expression for the power spectrum of present-day non-linear cosmic density fluctuations. It depends on the theory of gravity and the…
The evolution of disk galaxies in modified gravity is studied by using high-resolution N-body simulations. More specifically, we use the weak field limit of two modified gravity theories, i.e., nonlocal gravity (NLG) and…
We apply a halo velocity bias model, $\gamma_{f}$, within the Aemulus simulation suite for General Relativity (GR) to investigate its efficacy in identifying the signature of assembly bias and Modified Gravity (MG). In the investigation of…
We present an algorithm which is designed to allow the efficient identification and preliminary dynamical analysis of thousands of structures and substructures in large N-body simulations. First we utilise a refined density gradient system…
Measuring the sum of the three active neutrino masses, $M_\nu$, is one of the most important challenges in modern cosmology. Massive neutrinos imprint characteristic signatures on several cosmological observables in particular on the…
We present two new algebraic multilevel hierarchical matrix algorithms to perform fast matrix-vector product (MVP) for $N$-body problems in $d$ dimensions, namely efficient $\mathcal{H}^2_{*}$ (fully nested algorithm, i.e., $\mathcal{H}^2$…
Discovering the mass of neutrinos is a principle goal in high energy physics and cosmology. In addition to cosmological measurements based on two-point statistics, the neutrino mass can also be estimated by observations of neutrino wakes…
We present an approach for the inclusion of non-spherical constituents in high-resolution N-body discrete element method (DEM) simulations. We use aggregates composed of bonded spheres to model non-spherical components. Though the method…
We introduce $\texttt{Hi-COLA}$, a code designed to run fast, approximate $\textit{N}$-body simulations of non-linear structure formation in reduced Horndeski gravity. Given an input Lagrangian, $\texttt{Hi-COLA}$ dynamically constructs the…