Related papers: Fast and reliable symplectic integration for plane…
Numerical integration methods are central to the study of self-gravitating systems, particularly those comprised of many bodies or otherwise beyond the reach of analytical methods. Predictor-corrector schemes, both multi-step methods and…
We report on the properties of eclipsing binaries from the Kepler mission with a newly developed photometric modeling code, which uses the light curve, spectral energy distribution of each binary, and stellar evolution models to infer…
A pure two-body problem has seven integrals including the Kepler energy, the Laplace vector, and the angular momentum vector. However, only five of them are independent. When the five independent integrals are preserved, the two other…
The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the…
This paper presents key enhancements to our previous work~\cite{naghmouchi2024mixed} on a hybrid Benders decomposition (HBD) framework for solving mixed integer linear programs (MILPs). In our approach, the master problem is reformulated as…
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…
We propose a novel non-negative spherical relaxation for optimization problems over binary matrices with injectivity constraints, which in particular has applications in multi-matching and clustering. We relax respective binary matrix…
We study symplectic numerical integration of mechanical systems with a Hamiltonian specified in non-canonical coordinates and its application to guiding-center motion of charged plasma particles in magnetic confinement devices. The…
We present a novel method for efficient direct integration of gravitational N-body systems with a large variation in characteristic time scales. The method is based on a recursive and adaptive partitioning of the system based on the…
In this work, we study the stellar-dynamical hardening of unequal mass supermassive black hole (SMBH) binaries in the central regions of merging galactic nuclei. We present a comprehensive set of direct $N$-body simulations of the problem,…
We present implementations of a fourth-order symplectic integrator on graphic processing units for three $N$-body models with long-range interactions of general interest: the Hamiltonian Mean Field, Ring and two-dimensional self-gravitating…
We introduce new methods for the numerical solution of general Hamiltonian boundary value problems. The main feature of the new formulae is to produce numerical solutions along which the energy is precisely conserved, as is the case with…
We implement advanced Riemann solvers HLLC and HLLD \cite{Mignone:2005ft,MUB:2009} together with an advanced constrained transport scheme \cite{Gardiner:2007nc} in a numerical-relativity neutrino-radiation magnetohydrodynamics code. We…
We introduce a recent symplectic integration scheme derived for solving physically motivated systems with non-separable Hamiltonians. We show its relevance to Riemannian manifold Hamiltonian Monte Carlo (RMHMC) and provide an alternative to…
Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of…
It is well known that symplectic integrators lose their near energy preservation properties when variable step sizes are used. The most common approach to combine adaptive step sizes and symplectic integrators involves the Poincar\'e…
The Hyperspherical Harmonics (HH) method is one of the most accurate techniques to solve the quantum mechanical problem for nuclear systems with $A\le 4$. In particular, by applying the Rayleigh-Ritz or Kohn variational principle, both…
Upcoming Large Scale Structure surveys aim to achieve an unprecedented level of precision in measuring galaxy clustering. However, accurately modeling these statistics may require theoretical templates that go beyond second-order…
The superintegrability of several Hamiltonian systems defined on three-dimensional configuration spaces of constant curvature is studied. We first analyze the properties of the Killing vector fields, Noether symmetries and Noether momenta.…
We report new mid-eclipse times of the short-period sdB/dM binary HW Vir, which differ substantially from the times predicted by a previous model. The proposed orbits of the two planets in that model are found to be unstable. We present a…