Related papers: An adaptive algorithm for n-body field expansions
We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f_0. Our method uses importance sampling to find the sampling DF f_s that minimizes the…
This paper makes the first systematic attempt to determine using perturbation theory the positions of images by gravitational lensing due to arbitrary number of coplanar masses without any symmetry on a plane, as a function of lens and…
We propose a simple and computationally fast method for performing N-body simulations for a large class of modified gravity theories with a screening mechanism such as chameleons, symmetrons and galileons. By combining the linear…
We reexamine the recently introduced basis-set correction theory based on density-functional theory consisting in correcting the basis-set incompleteness error of wave-function methods using a density functional. We use a one-dimensional…
The envelope theory is a method to compute approximate eigensolutions of quantum $N$-body Hamiltonians with a quite general structure in $D$ dimensions. The advantages of the method are that it is easy to implement and that $N$ is treated…
Potential-density pair basis sets can be used for highly efficient N-body simulation codes, but they suffer from a lack of versatility, i.e. a basis set has to be constructed for each different class of stellar system. We present numerical…
We analyzed the performance of a perturbation theory for nonlinear cosmological dynamics, based on the Lagrangian description of hydrodynamics. In our previous paper, we solved hydrodynamic equations for a self-gravitating fluid with…
We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the…
In a Boson Sampling quantum optical experiment we send $n$ individual photons into an $m$-mode interferometer and we measure the occupation pattern on the output. The statistics of this process depending on the permanent of a matrix…
We present a simple and efficient method to set up spherical structure models for N-body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale…
The derivation of Debye shielding and Landau damping from the $N$-body description of plasmas requires many pages of heavy kinetic calculations in classical textbooks and is done in distinct, unrelated chapters. Using Newton's second law…
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…
We extensively develop a method of implementing mean-field calculations for deformed nuclei, using the Gaussian expansion method (GEM). This GEM algorithm has the following advantages: (i) it can efficiently describe the energy-dependent…
The problem of solving perturbatively the equations describing the evolution of self-gravitating collisionless matter in an expanding universe considerably simplifies when directly formulated in terms of the gravitational and velocity…
Astrophysical Challenges which demand the solution of the one million (or more) gravitating body problem are briefly discussed for the fields of cosmology, galactic nuclei and globular star clusters. Results from the classical three-body…
We describe and implement an adaptive particle-mesh algorithm to solve the Poisson equation for grid-based hydrodynamics codes with nested grids. The algorithm is implemented and extensively tested within the astrophysical code Enzo against…
We present a new scheme to compensate for the small-scales approximations resulting from Particle-Mesh (PM) schemes for cosmological N-body simulations. This kind of simulations are fast and low computational cost realizations of the large…
Tormen and Bertschinger have presented an algorithm which allows the dynamic range of N-body simulations to be extended by adding long-wavelength power to an evolved N-body simulation. This procedure is of considerable interest as it will…
We calculate the cosmological evolution of the 1-point probability distribution function (PDF), using an analytic approximation that combines gravitational perturbation theory with the Edgeworth expansion of the PDF. Our method applies…
The LCDM model is the most commonly admitted to describe our Universe. In spite of a great success with regard to the large scale structure formation, some problems are still unresolved at galactic scales. Alternative scenarios have to be…