Related papers: Speeding up N-body simulations of modified gravity…
Overcoming the time scale limitations of atomistics can be achieved by switching from the state-space representation of Molecular Dynamics (MD) to a statistical-mechanics-based representation in phase space, where approximations such as…
We present a novel convex formulation that weakly couples the Material Point Method (MPM) with rigid body dynamics through frictional contact, optimized for efficient GPU parallelization. Our approach features an asynchronous time-splitting…
We present a framework for the simulation of rigid and deformable bodies in the presence of contact and friction. Our method is based on a non-smooth Newton iteration that solves the underlying nonlinear complementarity problems (NCPs)…
Gravity inversion is the problem of estimating subsurface density distributions from observed gravitational field data. We consider the two-dimensional (2D) case, in which recovering density models from one-dimensional (1D) measurements…
We describe the use of Graphics Processing Units (GPUs) for speeding up the code NBODY6 which is widely used for direct $N$-body simulations. Over the years, the $N^2$ nature of the direct force calculation has proved a barrier for…
We present the numerical methods and GPU-accelerated implementation underlying a Total Lagrangian finite element framework for finite-deformation flexible multibody dynamics, introduced in the companion paper [1]. The framework supports…
We present the implementation of a fast estimator for the full dark matter bispectrum of a three-dimensional particle distribution relying on a separable modal expansion of the bispectrum. The computational cost of accurate bispectrum…
We present $\texttt{MG-evolution}$, an $N$-body code simulating the cosmological structure formation for parametrised modifications of gravity. It is built from the combination of parametrised linear theory with a parametrisation of the…
Direct gravitational simulations of n-body systems have a time complexity O(n^2), which gets computationally expensive as the number of bodies increases. Distributing this workload to multiple cores significantly speeds up the computation…
Recent advances in 3D Gaussian Splatting (3DGS) have focused on accelerating optimization while preserving reconstruction quality. However, many proposed methods entangle implementation-level improvements with fundamental algorithmic…
We investigate void properties in $f(R)$ models using N-body simulations, focusing on their differences from General Relativity (GR) and their detectability. In the Hu-Sawicki $f(R)$ modified gravity (MG) models, the halo number density…
We have developed a new three-dimensional algorithm, based on the standard P$^3$M method, for computing deflections due to weak gravitational lensing. We compare the results of this method with those of the two-dimensional planar approach,…
This paper is concerned with theories of gravity that contain a scalar coupled both conformally and disformally to matter through the metric. By systematically deriving the non-relativistic limit, it is shown that no new non-linear…
Full ray-tracing maps of gravitational lensing, constructed from N-Body simulations, represent a fundamental tool to interpret present and future weak lensing data. However the limitation of computational resources and storage capabilities…
We present a new method for generating initial conditions for numerical cosmological simulations in which massive neutrinos are treated as an extra set of N-body (collisionless) particles. It allows us to accurately follow the density field…
Learning dynamical systems that respect physical symmetries and constraints remains a fundamental challenge in data-driven modeling. Integrating physical laws with graph neural networks facilitates principled modeling of complex N-body…
In 'modified' gravity the observed acceleration of the universe is explained by changing the gravitational force law or the number of degrees of freedom in the gravitational sector. Both possibilities can be tested by measurements of…
In the fields of quantum mechanics and quantum information science, the traces of reduced density matrix powers play a crucial role in the study of quantum systems and have numerous important applications. In this paper, we propose a…
We present the first suite of cosmological N-body simulations that simultaneously include the effects of two different and theoretically independent extensions of the standard $\Lambda$CDM cosmological scenario - namely an $f(R)$ theory of…
Cosmic voids are progressively emerging as a new viable cosmological probe. Their abundance and density profiles are sensitive to modifications of gravity, as well as to dark energy and neutrinos. The main goal of this work is to…