Related papers: Fast and reliable symplectic integration for plane…
We explore the probable chemical signature of planet formation in the remarkable binary system HD 106515. The A star hosts a massive long-period planet with 9 MJup detected by radial velocity. We also refine stellar and planetary parameters…
A Bayesian re-analysis of published radial velocity data sets is providing evidence for additional planetary candidates. The nonlinear model fitting is accomplished with a new hybrid Markov chain Monte Carlo (HMCMC) algorithm which…
This paper uses the HCT finite element method and mesh adaptation technology to solve the nonlinear plate bending problem and conducts error analysis on the iterative method, including a priori and a posteriori error estimates. Our…
We implement coordinates suitable for studying wide binary systems in TRACE, a hybrid integrator in the widely used open-source N-body integration package REBOUND. This is a regime in which traditional hybrid integrators perform poorly. The…
We present a comprehensive exploration of hierarchical triple black hole (BH) systems to address the "initial separation" problem in gravitational wave astrophysics. This problem arises because isolated BH binaries must have extremely small…
In the first part of this paper, inspired by the geometric method of Jean-Pierre Marec, we consider the two-impulse Hohmann transfer problem between two coplanar circular orbits as a constrained nonlinear programming problem. By using the…
Binary embedding is the problem of mapping points from a high-dimensional space to a Hamming cube in lower dimension while preserving pairwise distances. An efficient way to accomplish this is to make use of fast embedding techniques…
In this paper we consider the numerical solution of the Hamiltonian wave equation in two spatial dimension. We use the Mimetic Finite Difference (MFD) method to approximate the continuous problem combined with a symplectic integration in…
The partial case of the planar $N+1$ body problem, $N\ge2$, of the type of planetary system with satellites is studied. One of the bodies (the Sun) is assumed to be much heavier than the other bodies ("planets" and "satellites"), moreover…
Optimization problems is one of the most challenging applications of quantum computers, as well as one of the most relevants. As a consequence, it has attracted huge efforts to obtain a speedup over classical algorithms using quantum…
The general consensus in the N-body community is that statistical results of an ensemble of collisional N-body simulations are accurate, even though individual simulations are not. A way to test this hypothesis is to make a direct…
We show that, when applied to any non-canonical Hamiltonian system, any integrator that is symplectic for canonical Hamiltonian problems is actually conjugate symplectic for the non-canonical structure. This result is useful because it…
In previous papers, explicit symplectic integrators were designed for nonrotating black holes, such as a Schwarzschild black hole. However, they fail to work in the Kerr spacetime because not all variables can be separable, or not all…
Mounting discoveries of debris discs orbiting newly-formed stars and white dwarfs (WDs) showcase the importance of modeling the long-term evolution of small bodies in exosystems. WD debris discs are in particular thought to form from very…
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…
We present a practical algorithm based on symplectic splitting methods to integrate numerically in time the Schr\"odinger equation. When discretized in space, the Schr\"odinger equation can be recast as a classical Hamiltonian system…
Dynamic simulation of elastic bodies is a longstanding task in engineering and computer graphics. In graphics, numerical integrators like implicit Euler and BDF2 are preferred due to their stability at large time steps, but they tend to…
High-intensity focused ultrasound (HIFU) is a promising treatment modality for the non-invasive ablation of pathogenic tissue in many organs. Optimal treatment planning strategies based on high-performance computing methods are expected to…
In this work, we systematically examine the application of spatio-temporal splitting heuristics to the Multi-Robot Motion Planning (MRMP) problem in a graph-theoretic setting: a problem known to be NP-hard to optimally solve. Following the…
For simulations that deal only with dark matter or stellar systems, the conventional N-body technique is fast, memory efficient, and relatively simple to implement. However when including the effects of gas physics, mesh codes are at a…