Related papers: Accelerating $N$-body simulation of self-gravitati…
We developed a new direct-tree hybrid N-body algorithm for fully self-consistent N-body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because…
A collisional $N$-body simulation using NBODY5 on a single CRAY YMP processor is followed well into the post-collapse regime. This is presently one of the largest particle numbers of all such models published, but some data for an even…
Fast N-body PM simulations with a small number of time steps such as FastPM or COLA have been remarkably successful in modeling the galaxy statistics, but their lack of small scale force resolution and long time steps cannot give accurate…
We report on the first fully consistent conventional cluster simulation which includes terms up to post^{5/2} Newtonian in the potential of the massive body. Numerical problems for treating extremely energetic binaries orbiting a single…
We train convolutional neural networks to correct the output of fast and approximate N-body simulations at the field level. Our model, Neural Enhanced COLA --NECOLA--, takes as input a snapshot generated by the computationally efficient…
Oscillation probability calculations are becoming increasingly CPU intensive in modern neutrino oscillation analyses. The independency of reweighting individual events in a Monte Carlo sample lends itself to parallel implementation on a…
A small body moving in the field of a much larger black hole and subjected to its own gravity moves on an accelerated world line in the background spacetime of the large black hole. The acceleration is produced by the body's gravitational…
Reliable predictions of general relativity theory are extracted using approximation methods. Among these, the powerful post-Newtonian approximation provides us with our best insights into the problems of motion and gravitational radiation…
In this work we present an efficient implementation of Canonical Monte Carlo simulation for Coulomb many body systems on graphics processing units (GPU). Our method takes advantage of the GPU Single Instruction, Multiple Data (SIMD)…
The conservation of energy, linear momentum and angular momentum are important drivers for our physical understanding of the evolution of the Universe. These quantities are also conserved in Newton's laws of motion under gravity…
We report on the results of a study of the motion of a four particle non-relativistic one-dimensional self-gravitating system. We show that the system can be visualized in terms of a single particle moving within a potential whose…
The reproducibility of experiments is one of the main principles of the scientific method. However, numerical N-body experiments, especially those of planetary systems, are currently not reproducible. In the most optimistic scenario, they…
We develop a close-limit approximation to the head-on collision of two neutron stars similar to that used to treat the merger of black hole binaries . This approximation can serve as a useful benchmark test for future fully nonlinear…
This paper presents, to the author's knowledge, the first graphics processing unit (GPU) accelerated program that solves the evolution of interacting scalar fields in an expanding universe. We present the implementation in NVIDIA's Compute…
We outline how auxiliary-field quantum Monte Carlo (AFQMC) can leverage graphical processing units (GPUs) to accelerate the simulation of solid state sytems. By exploiting conservation of crystal momentum in the one- and two-electron…
We have preliminary results on the parallelization of a Tree-Code for evaluating gravitational forces in N-body astrophysical systems. For our Cray T3D/CRAFT implementation, we have obtained an encouraging speed-up behavior, which reaches a…
This paper describes some applications of GPU acceleration in ab initio nuclear structure calculations. Specifically, we discuss GPU acceleration of the software package MFDn, a parallel nuclear structure eigensolver. We modify the matrix…
Molecular dynamics (MD) simulation is a powerful computational tool to study the behavior of macromolecular systems. But many simulations of this field are limited in spatial or temporal scale by the available computational resource. In…
Simulation speed depends on code structures, hence it is crucial how to build a fast algorithm. We solve the Allen-Cahn equation by an explicit finite difference method, so it requires grid calculations implemented by many for-loops in the…
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…