Related papers: GLISSE: A GPU-optimized planetary system integrato…
We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage…
The rapid advancement of large-scale cosmological simulations has opened new avenues for cosmological and astrophysical research. However, the increasing diversity among cosmological simulation models presents a challenge to the robustness.…
In this paper, we propose a $\mu$-mode integrator for computing the solution of stiff evolution equations. The integrator is based on a $d$-dimensional splitting approach and uses exact (usually precomputed) one-dimensional matrix…
We introduce version two of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). The first version (Alexander & Gieles) assumed that cluster evolution is balanced for the majority of the life-cycle, meaning that the…
We present GRIP, a graph neural network accelerator architecture designed for low-latency inference. AcceleratingGNNs is challenging because they combine two distinct types of computation: arithmetic-intensive vertex-centric operations and…
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.…
We introduce GYRE, a new open-source stellar oscillation code which solves the adiabatic/non-adiabatic pulsation equations using a novel Magnus Multiple Shooting (MMS) numerical scheme. The code has a global error scaling of up to 6th order…
The rapid analysis of ongoing gravitational microlensing events has been integral to the successful detection and characterisation of cool planets orbiting low mass stars in the Galaxy. In this paper we present an implementation of search…
We study the gravitational stability of gaseous streams in the complex environment of a galaxy merger, because mergers are known to be places of ongoing massive cluster formation and bursts of star formation. We find an analytic stability…
We describe the use of Graphics Processing Units (GPUs) for speeding up the code NBODY6 which is widely used for direct $N$-body simulations. Over the years, the $N^2$ nature of the direct force calculation has proved a barrier for…
Bayesian inference is the workhorse of gravitational-wave astronomy, for example, determining the mass and spins of merging black holes, revealing the neutron star equation of state, and unveiling the population properties of compact…
Long-term integrations of asteroid orbits with high-accuracy numerical integrators are essential for understanding dynamical evolution and ejection from the Solar System, but are computationally expensive. Here, we investigate the dynamical…
Graphics Processing Units (GPUs) are now powerful and flexible systems adapted and used for other purposes than graphics calculations (General Purpose computation on GPU -- GPGPU). We present here a prototype to be integrated into…
In this work we revisit the problem of the dynamical stability of hierarchical triple systems with applications to circumbinary planetary orbits. We carry out more than 3 10^8 numerical simulations of planets between the size of Mercury and…
We present a scalable dissipative particle dynamics simulation code, fully implemented on the Graphics Processing Units (GPUs) using a hybrid CUDA/MPI programming model, which achieves 10-30 times speedup on a single GPU over 16 CPU cores…
We describe the implementation and performance of the ${\rm P^3T}$ (Particle-Particle Particle-Tree) scheme for simulating dense stellar systems. In ${\rm P^3T}$, the force experienced by a particle is split into short-range and long-range…
Investigating uncertainties in computer simulations can be prohibitive in terms of computational costs, since the simulator needs to be run over a large number of input values. Building an emulator, i.e. a statistical surrogate model of the…
Transitional disks are protoplanetary disks with large and deep central holes in the gas, possibly carved by young planets. Dong, R., & Dawson, R. 2016, ApJ, 825, 7 simulated systems with multiple giant planets that were capable of carving…
Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations…
The implementation of the GRAVITY+ Adaptive Optics (GPAO) system at VLTI enables unprecedented sensitivity and stability in optical interferometry. This allows high-precision characterization of directly imaged exoplanets at medium spectral…