Related papers: Speeding up N-body simulations of modified gravity…
We introduce FastPM, a highly-scalable approximated particle mesh N-body solver, which implements the particle mesh (PM) scheme enforcing correct linear displacement (1LPT) evolution via modified kick and drift factors. Employing a…
This paper presents a fast, economical particle-multiple-mesh N-body code optimized for large-N modelling of collisionless dynamical processes, such as black-hole wandering or bar-halo interactions, occurring within isolated galaxies. The…
Newtonian N-body simulations have been employed successfully over the past decades for the simulation of the cosmological large-scale structure. Such simulations usually ignore radiation perturbations (photons and massless neutrinos) and…
We investigate a hybrid numerical algorithm aimed at the large-scale cosmological N-body simulation for the on-going and the future high precious sky surveys. It makes use of a truncated Fast Multiple Method (FMM) for short-range gravity,…
N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new…
We have developed a simulation code with the techniques which enhance both spatial and time resolution of the PM method for which the spatial resolution is restricted by the spacing of structured mesh. The adaptive mesh refinement (AMR)…
Fitting a theoretical model to experimental data in a Bayesian manner using Markov chain Monte Carlo typically requires one to evaluate the model thousands (or millions) of times. When the model is a slow-to-compute physics simulation,…
By using observations of the Hulse-Taylor pulsar we constrain the gravitational wave (GW) speed to the level of $10^{-2}$. We apply this result to scalar-tensor theories that generalize Galileon 4 and 5 models, which display anomalous…
The gravitational N-body simulation in the Solar system was performed using different parallel approaches with the comparisons in the computational times and speed-up values being carried out under different model sizes and the number of…
Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed…
We investigate nonlinear structure formation in Horndeski gravity with a luminal gravitational wave speed ($c_T = 1$) using the spherical collapse model incorporating Vainshtein screening. We compute the critical and virial overdensities…
Gravity inversion allows us to constrain the interior mass distribution of a planetary body using the observed shape, rotation, and gravity. Traditionally, techniques developed for gravity inversion can be divided into Monte Carlo methods,…
We use N-body simulations to study the statistics of massive halos and redshift space distortions for theories with a standard \Lambda CDM expansion history and a galileon-type scalar field. The extra scalar field increases the…
The N-body problem is a classic problem involving a system of N discrete bodies mutually interacting in a dynamical system. At any moment in time there are N*(N - 1)/2 such interactions occurring. This scaling as N^2 leads to computational…
The statistics of dark matter halos is an essential component of understanding the nonlinear evolution in modified gravity cosmology. Based on a series of modified gravity N-body simulations, we investigate the halo mass function,…
We compare void size and clustering statistics for nDGP and $f(R)$ gravity models and GR using N-body simulations. We show how it is critical to consider the statistics derived from mock galaxy catalogs rather than the dark matter halos…
Extending 3D Gaussian Splatting (3DGS) to 4D physical simulation remains challenging. Based on the Material Point Method (MPM), existing methods either rely on manual parameter tuning or distill dynamics from video diffusion models,…
A major aim of cosmological surveys is to test deviations from the standard $\Lambda$CDM model, but the full scientific value of these surveys will only be realised through efficient simulation methods that keep up with the increasing…
Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic…
The Variational Gaussian wavepacket (VGW) method is an alternative to Path Integral Monte-Carlo (PIMC) for the computation of thermodynamic properties of many-body systems at thermal equilibrium. It provides a direct access to the thermal…