Related papers: KRIOS: A new basis-expansion $N$-body code for col…
The dynamics of a large stellar (globular) cluster containing N=128,000 stars has been simulated by a direct summation (O(N$^2$)) method by using a heterogeneous platform. Preliminary simulations have been carried out on model systems with…
We present a new N-body and gas dynamics code, called Nyx, for large-scale cosmological simulations. Nyx follows the temporal evolution of a system of discrete dark matter particles gravitationally coupled to an inviscid ideal fluid in an…
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…
One of the fundamental aspects of statistical behaviour in many-body systems is exponential divergence of neighbouring orbits, which is often discussed in terms of Liapounov exponents. Here we study this topic for the classical…
We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm.…
We have developed a new numerical code to simulate the joint evolution of a massive black hole (MBH) and a surrounding stellar cluster at the center of a galaxy. The physics treated in the simulations include: 2-body relaxation, stellar…
Dark matter N-body simulations provide a powerful tool to model the clustering of galaxies and help interpret the results of galaxy redshift surveys. However, the galaxy properties predicted from N-body simulations are not necessarily…
Cosmological $N$-body simulations are the standard tool to study the emergence of the observed large-scale structure of the Universe. Such simulations usually solve for the gravitational dynamics of matter within the Newtonian…
We apply the joint lensing and dynamics code for the analysis of early-type galaxies, CAULDRON, to a rotating N-body stellar system with dark matter halo which significantly violates the two major assumptions of the method, i.e. axial…
We systematically study the effects of collisions on the overall dynamical evolution of dense star clusters using Monte Carlo simulations over many relaxation times. We derive many observable properties of these clusters, including their…
Our knowledge of the strengths of small bodies in the Solar System is limited by our poor understanding of their internal structures, and this, in turn, clouds our understanding of the formation and evolution of these bodies. Observations…
Planetary systems with several planets in compact orbital configurations such as TRAPPIST-1 are surely affected by tidal effects. Its study provides us with important insight about its evolution. We developed a second generation of a N-body…
The large dynamic range in some astrophysical N-body problems led to the use of adaptive multi-time-steps; however, the search for optimal strategies is still challenging. We numerically quantify the performance of the hierarchical…
For the Newtonian (gravitational) $n$-body problem in the Euclidean $d$-dimensional space, the simplest possible solutions are provided by those rigid motions (homographic solutions) in which each body moves along a Keplerian orbit and the…
We present a novel solver for an analogue to Poisson's equation in the framework of modified Newtonian dynamics (MOND). This equation is highly non-linear and hence standard codes based upon tree structures and/or FFT's in general are not…
Context. Simulating stellar dynamics in a molecular cloud environment is numerically challenging due to the strong coupling between young stars and their surrounding gas, and the large range of length and time scales. Aims. This paper is…
The distribution of galaxies, halo abundance, and peculiar velocities are influenced by non-linear gravitational interactions, making the study of non-linear evolution crucial for accurate cosmological predictions. We explore these aspects…
Direct gravitational simulations of n-body systems have a time complexity O(n^2), which gets computationally expensive as the number of bodies increases. Distributing this workload to multiple cores significantly speeds up the computation…
Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural…
We implement an efficient method to quantify time-dependent orbital complexity in gravitational $N$-body simulations. The technique, which we name DWaTIM, is based on a discrete wavelet transform of velocity orbital time series. The wavelet…