Related papers: A particle-based hybrid code for planet formation
The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that…
One of the current challenges of planet formation theory is to explain the enrichment of observed exoplanetary atmospheres. Past studies have focused on scenarios where either pebbles or planetesimals were the heavy element enrichment's…
We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N-body integrator, namely the Bulirsch--Stoer from the SWIFT package. We can then derive various observables:…
The final sizes, composition, and angular momenta of solid planetary bodies depend on the outcomes of collisions between planetary embryos. The most common numerical method for simulating embryo collisions is to combine a gravity solver…
Many exoplanetary systems are multiplanet configurations whose long-term dynamics are governed by N-body gravitational interactions. Consequently, their detection signatures cannot be adequately described by Keplerian orbits. Accurately…
We study the accumulation of errors in cosmological N-body algorithms that are caused by representing the continuous distribution of matter by massive particles, comparing the PPPM and Adaptive Multigrid codes. We use for this a new measure…
The formation of planetesimals requires the growth of dust particles through collisions. Micron-sized particles must grow by many orders of magnitude in mass. In order to understand and model the processes during this growth, the mechanical…
Planetary bodies form by accretion of smaller bodies. It has been suggested that a very efficient way to grow protoplanets is by accreting particles of size <<km (e.g., chondrules, boulders, or fragments of larger bodies) as they can be…
We present a novel implementation of a soft sphere, discrete elements code to simulate the dynamics of self-gravitating granular materials. The code is used to study the outcome of sub-sonic collisions between self-gravitating rubble piles…
At least 30\% of main sequence stars host planets with sizes of between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk--planet tidal…
We developed a Keplerian-based Hamiltonian splitting for solving the gravitational $N$-body problem. This splitting allows us to approximate the solution of a general $N$-body problem by a composition of multiple, independently evolved…
In this work, we provide a solution to the problem of computing collision stress in particle-tracking simulations. First, a formulation for the collision stress between particles is derived as an extension of the virial stress formula to…
The study of the interaction between solid objects and magnetohydrodynamic (MHD) fluids is of great importance in physics as consequence of the significant phenomena generated, such as planets interacting with stellar wind produced by their…
We review the basic dynamics and accretion of planetesimals by showing N-body simulations. The orbits of planetesimals evolve through two-body gravitational relaxation: viscous stirring increases the random velocity and dynamical friction…
Potential-density pair basis sets can be used for highly efficient N-body simulation codes, but they suffer from a lack of versatility, i.e. a basis set has to be constructed for each different class of stellar system. We present numerical…
Detailed modeling and simulation of biochemical systems is complicated by the problem of combinatorial complexity, an explosion in the number of species and reactions due to myriad protein-protein interactions and post-translational…
Robustly handling collisions between individual particles in a large particle-based simulation has been a challenging problem. We introduce particle merging-and-splitting, a simple scheme for robustly handling collisions between particles…
The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in…
Planets grow via the collisional accretion of small bodies in a protoplanetary disk. Such small bodies feel strong gas drag and their orbits are significantly affected by the gas flow and atmospheric structure around the planet. We…
We investigate planetesimal accretion via a direct N-body simulation of an annulus at 1 AU orbiting a 1 $M_{\odot}$ star. The planetesimal ring, which initially contains N = $10^6$ bodies is evolved into the oligarchic growth phase. Unlike…