Related papers: MLAPM - a C code for cosmological simulations
Execution-based evaluation of LLM-generated code implicitly treats successful execution as a proxy for correctness. In scientific simulation, this proxy is insufficient: a generated input file can run, mesh, and converge while encoding…
I present a large set of high resolution simulations, called CosmicGrowth Simulations, which were generated with either 8.6 billion or 29 billion particles. As the nominal cosmological model that can match nearly all observations on…
We present and test a code for two-fluid simulations of galaxy formation, one of the fluids being collision-less. The hydrodynamical evolution is solved through the SPH method while gravitational forces are calculated using a tree method.…
Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not…
We present a new Particle-Mesh cosmological N-body code for accurately solving the modified Poisson equation of the Quasi Linear formulation of MOND. We generate initial conditions for the Angus (2009) cosmological model, which is identical…
We establish, through coarse-grained computation, a connection between traditional, continuum numerical algorithms (initial value problems as well as fixed point algorithms) and atomistic simulations of the Larson model of micelle…
In the paper, gridless particle techniques are presented in order to solve problems involving electrostatic, collisionless plasmas. The method makes use of computational particles having the shape of spherical shells or of rings, and can be…
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…
In a previous paper we introduced a new method for simulating collisional gravitational $N$-body systems with linear time scaling on $N$, based on the Multi-Particle Collision (MPC) approach. This allows us to simulate globular clusters…
Many models beyond the Standard Model predict light and feebly interacting particles that are often long-lived. These long-lived particles (LLPs) in many cases can be produced from meson decays. In this work, we propose a simple and quick…
We present a new massively parallel code for N-body and cosmological hydrodynamical simulations of modified gravity models. The code employs a multigrid-accelerated Newton-Gauss-Seidel relaxation solver on an adaptive mesh to efficiently…
Deep cosmic microwave background polarization experiments allow a very precise internal reconstruction of the gravitational lensing signal in pricinple. For this aim, likelihood-based or Bayesian methods are typically necessary, where very…
High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on…
We present version 1.0 of the cosmological simulation code $\scriptstyle{\rm CO}N{\rm CEPT}$, designed for simulations of large-scale structure formation. $\scriptstyle{\rm CO}N{\rm CEPT}\, 1.0$ contains a P$^3$M gravity solver, with the…
SMILEI is a collaborative, open-source, object-oriented (C++) particle-in-cell code. To benefit from the latest advances in high-performance computing (HPC), SMILEI is co-developed by both physicists and HPC experts. The code's structures,…
Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input…
We present an interpretable emulator for the linear matter power spectrum (MPS) in the standard cosmological model $\Lambda$CDM, constructed via a physics-informed symbolic regression framework. By combining domain knowledge with a machine…
Precision measurements of the large scale structure of the Universe require large numbers of high fidelity mock catalogs to accurately assess, and account for, the presence of systematic effects. We introduce and test a scheme for…
The basic idea of fast Monte Carlo (MC) simulations is to perform particle-based MC simulations with the excluded-volume interactions modeled by "soft" repulsive potentials that allow particle overlapping. This gives much faster system…
A novel adaptive technique for electromagnetic Particle In Cell (PIC) plasma simulations is presented here. Two main issues are identified in designing adaptive techniques for PIC simulation: first, the choice of the size of the particle…