Related papers: Reducing phase error in long numerical binary blac…
Numerical effects are known to plague adaptive mesh refinement (AMR) codes when treating massive particles, e.g. representing massive black holes (MBHs). In an evolving background, they can experience strong, spurious perturbations and then…
We construct approximate initial data for non-spinning black hole binary systems by asymptotically matching the 4-metrics of two tidally perturbed Schwarzschild solutions in isotropic coordinates to a resummed post-Newtonian 4-metric in…
The success of the moving puncture method for the numerical simulation of black hole systems can be partially explained by the properties of stationary solutions of the 1+log coordinate condition. We compute stationary 1+log slices of the…
A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic…
We present a new iterative method to reduce eccentricity in black-hole-binary simulations. Given a good first estimate of low-eccentricity starting momenta, we evolve puncture initial data for ~4 orbits and construct improved initial…
In this article, we present a simple technique for boosting the order of accuracy of finite difference schemes for time dependent partial differential equations by optimally selecting the time step used to advance the numerical solution and…
We investigate the orbital decay of a massive BH embedded in a dark matter halo and a stellar bulge, using both analytical and numerical simulations with the aim of developing and validating a reliable dynamical friction (DF) correction…
We construct a closed-form, fully analytical 4-metric that approximately represents the spacetime evolution of non-precessing, spinning black hole binaries from infinite separations up to a few orbits prior to merger. We employ the…
Exponential time differencing methods is a power tool for high-performance numerical simulation of computationally challenging problems in condensed matter physics, fluid dynamics, chemical and biological physics, where mathematical models…
The characteristic approach to numerical relativity is a useful tool in evolving gravitational systems. In the past this has been implemented using two patches of stereographic angular coordinates. In other applications, a six-patch angular…
Tidal coupling between members of a compact binary system can have an interesting and important influence on that binary's dynamical inspiral. Tidal coupling also distorts the binary's members, changing them (at lowest order) from spheres…
The properties of modified Hayward black hole space-time can be investigated through analyzing the particle geodesics. By means of a detailed analysis of the corresponding effective potentials for a massive particle, we find all possible…
We introduce new numerical integration operators which compose the mass and stiffness matrices of a modified spectral element method for simulation of elastic wave propagation. While these operators use the same quadrature nodes as does the…
We present a novel approach to the numerical computation of quasi-normal modes, based on the first-order (in radial derivative) formulation of the equations of motion and using a matrix version of the continued fraction method. This…
Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein's equations with the harmonic slicing condition is presented.…
We show that puncture data for quasicircular binary black hole orbits allow a special gauge choice that realizes some of the necessary conditions for the existence of an approximate helical Killing vector field. Introducing free parameters…
Classical and new numerical schemes are generated using evolutionary computing. Differential Evolution is used to find the coefficients of finite difference approximations of function derivatives, and of single and multi-step integration…
We describe a grid generation procedure designed to construct new classes of orthogonal coordinate systems for binary black hole spacetimes. The computed coordinates offer an alternative approach to current methods, in addition to providing…
Stencil composition uses the idea of function composition, wherein two stencils with arbitrary orders of derivative are composed to obtain a stencil with a derivative order equal to sum of the orders of the composing stencils. In this…
The luminosity measurement at the projected International Linear $e^+e^-$ Collider ILC is planned to be performed with forward Bhabha scattering with an accuracy of the order of $10^{-4}$. A theoretical prediction of the differential…