Related papers: The Outer Rim Simulation: A Path to Many-Core Supe…
We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same…
Cosmological N-body simulations rank among the most computationally intensive efforts today. A key challenge is the analysis of structure, substructure, and the merger history for many billions of compact particle clusters, called halos.…
We report scaling results on the world's largest supercomputer of our recently developed Billions-Body Molecular Dynamics (BBMD) package, which was especially designed for massively parallel simulations of the atomic dynamics in structural…
We introduce Project GIBLE (Gas Is Better resoLved around galaxiEs), a suite of cosmological zoom-in simulations where gas in the circumgalactic medium (CGM) is preferentially simulated at ultra-high numerical resolution. Our initial sample…
Constraining neutrino mass through cosmological observations relies on precise simulations to calibrate their effects on large scale structure, while these simulations must overcome computational challenges like dealing with large velocity…
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…
High performance computing numerical simulations are today one of the more effective instruments to implement and study new theoretical models, and they are mandatory during the preparatory phase and operational phase of any scientific…
We profile the impact of computation and inter-processor communication on the energy consumption and on the scaling of cortical simulations approaching the real-time regime on distributed computing platforms. Also, the speed and energy…
We report on improvements made over the past two decades to our adaptive treecode N-body method (HOT). A mathematical and computational approach to the cosmological N-body problem is described, with performance and scalability measured up…
Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and…
We present a scheme to extend the halo mass resolution of N-body simulations of the hierarchical clustering of dark matter. The method uses the density field of the simulation to predict the number of sub-resolution dark matter haloes…
Galaxy simulations have come a long way from the early days of simple N-body calculations, which considered only gravitational interactions, to the complex, multi-physics models used today. Beginning with initial conditions representative…
We present recent updates and improvements of the graphical processing unit (GPU) N-body code GENGA. Modern state-of-the-art simulations of planet formation require the use of a very high number of particles to accurately resolve planetary…
Upcoming galaxy surveys will bring a wealth of information about the clustering of matter, but modeling small-scale structure beyond $\Lambda$CDM remains computationally challenging. While accurate N-body emulators exist to model the matter…
The Millennium Run is the largest simulation of the formation of structure within the $\Lambda$CDM cosmogony so far carried out. It uses $10^{10}$ particles to follow the dark matter distribution in a cubic region 500$h^{-1}$Mpc on a side,…
We present a method to construct semi-numerical ``simulations'', which can efficiently generate realizations of halo distributions and ionization maps at high redshifts. Our procedure combines an excursion-set approach with first-order…
To help understand our universe better, researchers and scientists currently run extreme-scale cosmology simulations on leadership supercomputers. However, such simulations can generate large amounts of scientific data, which often result…
We present a new C++ code for collisional N-body simulations of star clusters. The code uses the Hermite fourth-order scheme with block time steps, for advancing the particles in time, while the forces and neighboring particles are computed…
Approximately half of the matter in the Universe is "unbound" at z = 0, according to N-body simulations such as the Millennium Run. Here, we use the milli-Millennium simulation to examine the distribution of unbound matter in relation to…
We point out that superconducting quantum computers are prospective for the simulation of the dynamics of spin models far from equilibrium, including nonadiabatic phenomena and quenches. The important advantage of these machines is that…