Related papers: sCOLA: The N-body COLA Method Extended to the Spat…
The gravitationally-driven evolution of cold dark matter dominates the formation of structure in the Universe over a wide range of length scales. While the longest scales can be treated by perturbation theory, a fully quantitative…
The determination of the resolution of cosmological N-body simulations, i.e., the range of scales in which quantities measured in them represent accurately the continuum limit, is an important open question. We address it here using…
We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes…
We present numerical $N$-body simulation studies of large-scale structure formation. The main purpose of these studies is to analyze the several models of dark matter and the role they played in the process of large-scale structure…
The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics.…
Vision-language models (VLMs) such as CLIP demonstrate strong generalization in zero-shot classification but remain highly vulnerable to adversarial perturbations. Existing methods primarily focus on adversarial fine-tuning or prompt…
The cosmic large-scale structure (LSS) provides a unique testing ground for connecting fundamental physics to astronomical observations. Modeling the LSS requires numerical $N$-body simulations or perturbative techniques that both come with…
We present a new hybrid Smoothed Particle Hydrodynamics (SPH)/N-body method for modelling the collisional stellar dynamics of young clusters in a live gas background. By deriving the equations of motion from Lagrangian mechanics we obtain a…
Accurate predictions of weak lensing observables are essential for understanding the large-scale structure of the Universe and probing the nature of gravity. In this work, we present a lightcone implementation to generate maps of the weak…
We present numerical N-body simulation studies of large-scale structure formation. The main purpose of these studies is to analyze the several models of dark matter and the role they played in the process of large-scale structure formation.…
This posting announces public availability of the GALAXY software package developed by the author over the past 40 years. It is a highly efficient code for the evolution of (almost) isolated, collisionless stellar systems, both disk-like…
Approximations to the exact solutions for gravitational instability in the expanding Universe are extremely useful for understanding the evolution of large--scale structure. We report on a series of tests of Newtonian Lagrangian…
We train a neural network model to predict the full phase space evolution of cosmological N-body simulations. Its success implies that the neural network model is accurately approximating the Green's function expansion that relates the…
We build a field level emulator for cosmic structure formation that is accurate in the nonlinear regime. Our emulator consists of two convolutional neural networks trained to output the nonlinear displacements and velocities of N-body…
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 discuss the performance of direct summation codes used in the simulation of astrophysical stellar systems on highly distributed architectures. These codes compute the gravitational interaction among stars in an exact way and have an…
Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This…
We present a novel $N$-body simulation method that compactifies the infinite spatial extent of the Universe into a finite sphere with isotropic boundary conditions to follow the evolution of the large-scale structure. Our approach…
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…
Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogs to validate data analysis pipelines and identify potential systematics. We…