Related papers: A Lagrangian Integrator for Planetary Accretion an…
We study the collisional evolution of km-sized planetesimals in tight binary star systems to investigate whether accretion towards protoplanets can proceed despite the strong gravitational perturbations from the secondary star. The orbits…
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…
We present a new method to calculate formation of cosmological structure in the Newtonian limit. The method is based on Lagrangian perturbation theory plus two key theoretical extensions. One advance involves identifying and fixing a…
We present a code allowing to evolve three-dimensional self-gravitating collisionless systems with a large number of particles N >= 10^7. FLY (Fast Level-based N-bodY code) is a fully parallel code based on a tree algorithm. It adopts…
The particle proper orthogonal decomposition (PPOD) is demonstrated on cases of particle flows in decaying homogeneous isotropic turbulence. Data is generated through one-way coupled simulations, where particle positions and velocities are…
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…
When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model…
It is well-known that a Lagrangian induces a compatible presymplectic form on the equation manifold (stationary surface, understood as a submanifold of the respective jet-space). Given an equation manifold and a compatible presymplectic…
Infinitesimal contraction analysis provides exponential convergence rates between arbitrary pairs of trajectories of a system by studying the system's linearization. An essentially equivalent viewpoint arises through stability analysis of a…
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…
In this work, we adapt a module for planetary formation within the hydrodynamic code FARGO3D. Planetary formation is modeled by a solid core accretion scenario, with the core growing in oligarchic regime. The initial superficial density of…
We report experimental results on the three dimensional Lagrangian acceleration in highly turbulent flows. Tracer particles are tracked optically using four silicon strip detectors from high energy physics that provide high temporal and…
We describe an updated version of our hybrid N-body-coagulation code for planet formation. In addition to the features of our 2006-2008 code, our treatment now includes algorithms for the 1D evolution of the viscous disk, the accretion of…
The gravitational $N$-body problem is a nearly universal problem in astrophysics which, despite its deceptive simplicity, still presents a significant computational challenge. For collisional systems such as dense star clusters, the need to…
The shearing sheet is a model dynamical system that is used to study the small-scale dynamics of astrophysical disks. Numerical simulations of particle trajectories in the shearing sheet usually employ the leapfrog integrator, but this…
Path integral molecular dynamics (PIMD), which maps a quantum particle onto a fictitious classical system of ring polymers and propagates the "beads" of this extended classical system using molecular dynamics, is widely used to capture…
We developed a new direct-tree hybrid N-body algorithm for fully self-consistent N-body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because…
Modern N-body techniques for planetary dynamics are generally based on symplectic algorithms specially adapted to the Kepler problem. These methods have proven very useful in studying planet formation, but typically require the timestep for…
We examine the structure of the Clifford algebra associated with a Hermitian bilinear form and apply the result to a dynamical model of the relativistic point particle. The dynamics of the particle is described by a Dirac spinor with…
We present a new method of identifying protostellar disc fragments in a simulation based on density derivatives, and analyse our data using this and the existing CLUMPFIND method, which is based on an ordered search over all particles in…