GPU accelerated manifold correction method for spinning compact binaries

The graphics processing unit (GPU) acceleration of the manifold correction algorithm based on the compute unified device architecture (CUDA) technology is designed to simulate the dynamic evolution of the Post-Newtonian (PN) Hamiltonian formulation of spinning compact binaries. The feasibility and the efficiency of parallel computation on GPU have been confirmed by various numerical experiments. The numerical comparisons show that the accuracy on GPU execution of manifold corrections method has a good agreement with the execution of codes on merely central processing unit (CPU-based) method. The acceleration ability when the codes are implemented on GPU can increase enormously through the use of shared memory and register optimization techniques without additional hardware costs, implying that the speedup is nearly 13 times as compared with the codes executed on CPU for phase space scan (including 314×314$314 \times 314$ orbits). In addition, GPU-accelerated manifold correction method is used to numerically study how dynamics are affected by the spin-induced quadrupole–monopole interaction for black hole binary system.

[1]  Xin Wu,et al.  Resurvey of order and chaos in spinning compact binaries , 2008 .

[2]  Coalescing binary systems of compact objects to (post)5/2-Newtonian order. V. Spin effects. , 1995, Physical review. D, Particles and fields.

[3]  Diego Rossinelli,et al.  GPU accelerated simulations of bluff body flows using vortex particle methods , 2010, J. Comput. Phys..

[4]  Xin Wu,et al.  Application of the logarithmic Hamiltonian algorithm to the circular restricted three-body problem with some post-Newtonian terms , 2015 .

[5]  Y. Zhu,et al.  Effects of manifold correction methods on chaos indicators , 2015 .

[6]  F. Rasio,et al.  Ruling out chaos in compact binary systems. , 2001, Physical review letters.

[7]  Xin Wu,et al.  Revisit on ``Ruling out chaos in compact binary systems'' , 2007 .

[8]  Second post-Newtonian gravitational wave polarizations for compact binaries in elliptical orbits , 2001, gr-qc/0110100.

[9]  Xin Wu,et al.  A Note on the Equivalence of Post-Newtonian Lagrangian and Hamiltonian Formulations , 2016 .

[10]  Xin Wu,et al.  EXPLICIT SYMPLECTIC-LIKE INTEGRATORS WITH MIDPOINT PERMUTATIONS FOR SPINNING COMPACT BINARIES , 2017 .

[11]  Xin Wu,et al.  Dynamics of spin effects of compact binaries , 2013 .

[12]  Xin Wu,et al.  New adaptive time step symplectic integrator: an application to the elliptic restricted three-body problem ⁄ , 2014 .

[13]  R. F. O’Connell,et al.  The gravitational interaction: Spin, rotation, and quantum effects-a review , 1979 .

[14]  Chaos and order in models of black hole pairs , 2006, gr-qc/0612003.

[15]  Xin Wu,et al.  Dynamics of the post-Newtonian circular restricted three-body problem with compact objects , 2014 .

[16]  P. Nacozy.  The use of integrals in numerical integrations of theN-body problem , 1971 .

[17]  Xinfa Deng,et al.  Global symplectic structure-preserving integrators for spinning compact binaries , 2010 .

[18]  D. G. Bettis,et al.  Recent developments of integrating the gravitational problem ofN-bodies , 1971 .

[19]  T. Fukushima.  Efficient Orbit Integration by Dual Scaling for Consistency of Kepler Energy and Laplace Integral , 2003 .

[20]  Xin Wu,et al.  Implementation of the velocity scaling method for elliptic restricted three-body problems , 2016 .

[21]  Xin Wu,et al.  Simulations of Dissipative Circular Restricted Three-body Problems Using the Velocity-scaling Correction Method , 2018 .

[22]  R. Ruth,et al.  Fourth-order symplectic integration , 1990 .

[23]  Amitabh Varshney,et al.  Parallel, stochastic measurement of molecular surface area. , 2008, Journal of molecular graphics & modelling.

[24]  Applying explicit symplectic-like methods to nonconservative nonseparable systems , 2017 .

[25]  Xuan Li,et al.  GPU‐accelerated finite‐difference time‐domain method for dielectric media based on CUDA , 2016 .

[26]  J. Baumgarte,et al.  Numerical stabilization of the differential equations of Keplerian motion , 1972 .

[27]  Transition from inspiral to plunge in precessing binaries of spinning black holes , 2005, gr-qc/0508067.

[28]  Thibault Damour,et al.  Coalescence of two spinning black holes: an effective one-body approach , 2001, gr-qc/0103018.

[29]  Fate of chaotic binaries , 2000, gr-qc/0010100.

[30]  Xin Wu,et al.  On preference of Yoshida construction over Forest–Ruth fourth-order symplectic algorithm , 2013 .

[31]  J. Wisdom,et al.  Symplectic maps for the N-body problem. , 1991 .

[32]  Eric Darve,et al.  Large calculation of the flow over a hypersonic vehicle using a GPU , 2008, J. Comput. Phys..

[33]  Xin Wu,et al.  Computation of Lyapunov exponents in general relativity , 2003 .

[34]  Rodney A. Kennedy,et al.  Parallel computation of mutual information on the GPU with application to real-time registration of 3D medical images , 2010, Comput. Methods Programs Biomed..

[35]  Xin Wu,et al.  Extending Nacozy’s Approach to Correct All Orbital Elements for Each of Multiple Bodies , 2008 .

[36]  Wen-biao Han,et al.  The adjustment-stabilization method for constrained systems , 2007, Comput. Phys. Commun..

[37]  Igor Omelyan,et al.  Optimized Forest–Ruth- and Suzuki-like algorithms for integration of motion in many-body systems , 2002 .

[38]  Xin Wu,et al.  Chaos in two black holes with next-to-leading order spin–spin interactions , 2014, 1403.0378.

[39]  Levin.  Gravity waves, chaos, and spinning compact binaries , 2000, Physical review letters.

[40]  S. McWilliams,et al.  Inspiral of generic black hole binaries: spin, precession and eccentricity , 2010, 1009.2533.

[41]  Xin Wu,et al.  Second post-Newtonian Lagrangian dynamics of spinning compact binaries , 2016, 1604.05810.

[42]  Yi Xie,et al.  Symplectic structure of post-Newtonian Hamiltonian for spinning compact binaries , 2010, 1004.4549.

[43]  Xin Wu,et al.  Manifold corrections on spinning compact binaries , 2010 .

[44]  Uri M. Ascher,et al.  Stabilization of invariants of discretized differential systems , 1997, Numerical Algorithms.

[45]  Alice C. Quillen,et al.  QYMSYM: A GPU-accelerated hybrid symplectic integrator that permits close encounters , 2010, 1007.3458.

[46]  Fei Zhang.  Energy corrections in Hamiltonian dynamics simulations , 1996 .

[47]  P. Nacozy.  The Use of Integrals in Numerical Integrations of the N-Body Problem (Papers appear in the Proceedings of IAU Colloquium No. 10 Gravitational N-Body Problem (ed. by Myron Lecar), R. Reidel Publ. Co. , Dordrecht-Holland.) , 1971 .

[48]  Xin Wu,et al.  Lyapunov indices with two nearby trajectories in a curved spacetime , 2006, 1006.5251.

[49]  Quan Zhang,et al.  A fast calculation method of optical transfer function using GPU parallel computation , 2015 .