An Introduction to Common Numerical Integration Codes Used in Dynamical Astronomy

GANBISS: a new GPU accelerated N-body code for binary star systems

We present a GPU accelerated N-body integrator using the Bulirsch-Stoer method, called GANBISS (GPU accelerated n-body code for binary star systems). It is designed to simulate the dynamical evolution of planetesimal disks in binary star systems which contain some thousand disk objects. However, it can also be used for studies of non-interacting massless bodies where up to 50 million objects can be studied in a simulation. GANBISS shows the energy and angular momentum conservation behavior of non-symplectic integration methods. The code is written in CUDA C and can be run on NVIDIA GPUs of compute capability of at least 3.5. A comparison of GPU and CPU computations indicates a speed-up of the GPU performance of up to 100 times-depending on the number of disk objects.

Can Asteroid Belts Exist in the Luyten’s System?

The extra-solar planetary system Luyten is relatively close (12.3 light years) to our Sun. The Luyten’s red dwarf star is orbited by four planets, two of them Earth-like (in mass) and in 4:1 resonance. Extra-solar systems might contain asteroid belts such as ours. Therefore, it is important to investigate whether it is possible to have a stable population of minor bodies and compare them to those in our system. The study of extra-solar systems is crucial for understanding the evolution of planetary systems in general. Here, we investigate the stability of two possible asteroid populations in the Luyten’s system: the main asteroid belt between the two inner and two outer planets, and an outer asteroid belt, situated beyond the planets. We also explore the likelihood of observing an asteroid or a dwarf planet in this system. Our study suggests that the existence of asteroid belts is possible, notably the main belt at 0.09–0.53 au from the star and an outer belt (with the inner boundary at 0.85 au and the outer boundary at ∼66,000 au). The average Yarkovsky drift for the Luyten’s main asteroid belt is ∼0.5×10−4 au/Myr for km-size objects. The Luyten’s system might host extra-solar minor bodies, some of which could be capable of entering our own system. Presently, no asteroids can be detected in the Luyten’s system, not even a Ceres-sized body, because the detection signal using the radial velocity method is at least two orders of magnitude less than that required for discerning such objects. The detection probability of an asteroid in the Luyten belt similar to Ceres is about 1.3%, which is less than the probability of finding Luyten B (∼3%).

On the Possibility of Habitable Trojan Planets in Binary Star Systems

Approximately 60% of all stars in the solar neighbourhood (up to 80% in our Milky Way) are members of binary or multiple star systems. This fact led to the speculations that many more planets may exist in binary systems than are currently known. To estimate the habitability of exoplanetary systems, we have to define the so-called habitable zone (HZ). The HZ is defined as a region around a star where a planet would receive enough radiation to maintain liquid water on its surface and to be able to build a stable atmosphere. We search for new dynamical configurations-where planets may stay in stable orbits-to increase the probability to find a planet like the Earth.