Related papers: Gurzadyan's Problem 5 and improvement of softening…
Pairwise forces between particles in cosmological N-body simulations are generally softened to avoid hard collisions. Physically, this softening corresponds to treating the particles as diffuse clouds rather than point masses. For particles…
Cosmological simulations of structure formation follow the collisionless evolution of dark matter starting from a nearly homogeneous field at early times down to the highly clustered configuration at redshift zero. The density field is…
In self-consistent N-body simulations of collisionless systems, gravitational interactions are modified on small scales to remove singularities and simplify the task of numerically integrating the equations of motion. This `gravitational…
Gravitational softening length is one of the key parameters to properly set up a cosmological $N$-body simulation. In this paper, we perform a large suit of high-resolution $N$-body simulations to revise the optimal softening scheme…
The gravitational potential is a key function involved in many astrophysical problems. Its evaluation inside continuous media from Newton's law is known to be challenging because of the diverging kernel 1/|r-r'|. This difficulty is…
We study the accumulation of errors in cosmological N-body algorithms that are caused by representing the continuous distribution of matter by massive particles, comparing the PPPM and Adaptive Multigrid codes. We use for this a new measure…
In N-body simulations of collisionless stellar systems, the forces are softened to reduce the shot noise. Softening modifies gravity at r=|x-y| smaller than softening length epsilon and the softened forces are increasingly biased for ever…
Modeling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve…
We show that the holographic Complexity = Volume proposal satisfies a very general notion of Momentum/Complexity correspondence (PC), based on the Momentum Constraint of General Relativity. It relates the rate of complexity variation with…
We review various approaches to approximating general relativistic effects in hydrodynamic simulations of stellar core collapse and post-bounce evolution. Different formulations of a modified Newtonian gravitational potential are presented.…
We present a two-dimensional (2D) Particle-Particle-Particle-Mesh (P$^3$M) algorithm with an optimized Green function and adaptive softening length for gravitational lensing studies in N-Body simulations. The analytical form of the…
We propose a symmetrized form of the softened gravitational potential which is a natural extension of the Plummer potential. The gravitational potential at the position of particle i (x_i,y_i,z_i), induced by particle j at (x_j,y_j,z_j), is…
The post-Newtonian (PN) perturbative framework has been successful in understanding the slow-motion, weak field limit of Einstein's theory of gravity on solar system scales, and for isolated astrophysical systems. The parameterized…
The Parameterised Post-Newtonian (PPN) approach is the default framework for performing precision tests of gravity in nearby astrophysical systems. In recent works we have extended this approach for cosmological applications, and in this…
Einstein's theory of gravity has been extensively tested on solar system scales, and for isolated astrophysical systems, using the perturbative framework known as the parameterized post-Newtonian (PPN) formalism. This framework is designed…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…
We develop an implementable stochastic proximal point (SPP) method for a class of weakly convex, composite optimization problems. The proposed stochastic proximal point algorithm incorporates a variance reduction mechanism and the resulting…
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,…
Many recent studies have highlighted certain failures of the standard Eulerian-space cosmological perturbation theory (SPT). Its problems include (1) not capturing large-scale bulk flows [leading to an O(1) error in the 1-loop SPT…
We compare the cosmological first-order post-Newtonian (1PN) approximation with the relativistic cosmological linear perturbation theory in a zero-pressure medium with the cosmological constant. We compare equations and solutions in several…