Next Article in Journal
Islands in Generalized Dilaton Theories
Next Article in Special Issue
Spin-Gravity Coupling in a Rotating Universe
Previous Article in Journal
Fixed-Point Estimation by Iterative Strategies and Stability Analysis with Applications
Previous Article in Special Issue
A Simple Cosmology in G4v
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian

1
School of Physics, Sun Yat-sen University, Guangzhou 510275, China
2
School of Science, Sun Yat-sen University, Shenzhen 518107, China
3
School of Physics and Astronomy, Sun Yat-sen University, Zhuhai 519000, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2023, 15(7), 1401; https://doi.org/10.3390/sym15071401
Submission received: 16 May 2023 / Revised: 6 July 2023 / Accepted: 7 July 2023 / Published: 11 July 2023
(This article belongs to the Special Issue Symmetry in Gravity Research)

Abstract

:
The equations of motion for a Lagrangian mainly refer to the acceleration equations, which can be obtained by the Euler–Lagrange equations. In the post-Newtonian Lagrangian form of general relativity, the Lagrangian systems can only maintain a certain post-Newtonian order and are incoherent Lagrangians since the higher-order terms are omitted. This truncation can cause some changes in the constant of motion. However, in celestial mechanics, Hamiltonians are more commonly used than Lagrangians. The conversion from Lagrangianto Hamiltonian can be achieved through the Legendre transformation. The coordinate momentum separable Hamiltonian can be computed by the symplectic algorithm, whereas the inseparable Hamiltonian can be used to compute the evolution of motion by the phase-space expansion method. Our recent work involves the design of a multi-factor correction map for the phase-space expansion method, known as the correction map method. In this paper, we compare the performance of the implicit algorithm in post-Newtonian Lagrangians and the correction map method in post-Newtonian Hamiltonians. Specifically, we investigate the extent to which both methods can uphold invariance of the motion’s constants, such as energy conservation and angular momentum preservation. Ultimately, the results of numerical simulations demonstrate the superior performance of the correction map method, particularly with respect to angular momentum conservation.

1. Introduction

Compact binary systems, composed of neutron stars or black holes, etc., are of immense interest to experimental and theoretical researchers as sources of gravitational waves for broadband laser interferometers. The temporal progression of binary systems encompassing compact objects can be elucidated through the implementation of Einstein’s equations of general relativity. Explicit symplectic integrators are supposed to be the ideal candidate with several benefits for the numerical simulations in these systems. They are designed to preserve the symplectic structure, which guarantees the precision and stability of numerical solutions over long time intervals. However, Einstein’s equations of general relativity describe the motion of strong gravitational systems for which exact solutions are very difficult to obtain, but there have been some efforts to address this issue. For example, Xin Wu [1] developed the explicit symplectic integrators for Hamiltonian systems in curved spacetimes, particularly for black hole spacetimes. The paper [2] discusses the motion of charged particles around the Schwarzschild black hole with an external magnetic field. Therefore, the orbital motion described in this paper is a two-body problem, involving the motion of a charged particle around a Schwarzschild black hole. Ying wang [3,4] constructs the explicit symplectic integrators in general relativity, specifically for the Hamiltonian of Schwarzschild spacetime geometry. The integrators are useful for the long-term integration of N-body Hamiltonian systems and modeling the chaotic motion of charged particles around a black hole with an external magnetic field. Xin Wu [5] discusses the construction of explicit symplectic integrators for Kerr black holes in general relativity. The authors introduce a time transformation function to the Hamiltonian of Kerr geometry to obtain a time-transformed Hamiltonian consisting of five splitting parts whose analytical solutions are explicit functions of the new coordinate time. Wei Sun [6] proposes an explicit symplectic integrator for the Kerr spacetime geometry to simulate the nonintegrable dynamics of charged particles moving around the Kerr black hole in an external magnetic field. The algorithm shows good numerical performance and is used to study the dynamics of order and chaos of charged particles.
Despite the efficacy of this method, another useful and well-developed approach is the form of the post-Newtonian (PN) Lagrangian or Hamiltonian [7] approximation [7,8,9]. Arun et al. [10] investigated inspiralling compact binaries in quasi-elliptical orbits and provided a comprehensive analysis of the third post-Newtonian energy flux. Their study focused on understanding the energy loss due to gravitational radiation and its implications for compact binary systems. Tessmer and Schäfer [11] studied the eccentric motion of spinning compact binaries. They examined the dynamics of these systems with non-circular orbits, considering the effects of spin and exploring the consequences of eccentricity on the gravitational wave signals emitted during inspiral. Hinder et al. [12] developed an eccentric binary black hole waveform model by combining numerical relativity simulations with post-Newtonian theory. Their work aimed to accurately describe the complete inspiral–merger–ringdown phase of eccentric binary black hole systems, providing insights into the gravitational waveforms emitted during these events. Chattaraj et al. [13] conducted high-accuracy comparisons between post-Newtonian theory and numerical relativity simulations, specifically focusing on eccentric binary black holes. They investigated the influence of higher modes on the waveforms and developed a model that incorporates eccentricity and accurately describes the inspiral, merger, and ringdown phases. Chowdhury and Khlopov [14] studied an eccentric binary black hole system within the framework of post-Newtonian theory. Their research aimed to understand the behavior of binary black holes with non-circular orbits, providing insights into the dynamics and gravitational wave emissions of eccentric binary systems. These approximations provide high-precision theoretical templates of gravitational waveforms, although their higher-order terms are truncated, which affects their equivalence [15,16,17].
The choice of approximation and the selection of the algorithm becomes crucial to ensuring an accurate and effective description of the trajectory evolution of compact binaries systems and matching corresponding gravitational waveforms.
The PN Lagrangian equations of motion are derived from the Euler–Lagrangian equations of a PN Lagrangian formulation, denoted as L ( r , v ) . By calculating the partial derivative of L with respect to velocity, we obtain the generalized momentum p = L / v . Similarly, the acceleration equations a = f ( r , v , a ) = L / r can be derived and we can obtain a coherent Lagrangian [18,19,20]. By limiting the inclusion of accelerations up to a certain PN order in the Lagrangian, the accelerations a in the function f will be modified to a ; it only has lower-order terms, i.e., a = f ( r , v , a ) .The acceleration equations become incoherent, due to the higher-order PN terms disappearing, leading to a loss of some values of the constants of motion during subsequent evolution. The same problem occurs with the post-Newtonian Hamiltonian form. The error of the constant of motion can be used as an indicator to test the performance of different algorithms in both approximate forms.
Various algorithms are available for the calculation of post-Newtonian Lagrangian quantities. For example, in optimizing the fifth-order Runge–Kutta method as a high-precision integrator, Zhong [21] employed corrections to all integrals within the conservative 3PN order Hamiltonian. Tsang [22] introduced an implicit symplectic integrator that accounts for 2.5PN gravitational radiation reaction terms in the Newtonian two-body problem. This approach effectively captures the effects of radiation reactions. Lubich [23] devised an explicit and implicit mixed symplectic integration technique that facilitates the splitting of orbital and spin contributions. By employing this approach, the dynamics of both orbital and spin variables can be accurately simulated. Zhong [24] proposed fourth-order canonical explicit and implicit mixed symplectic methods. These methods offer improved accuracy and stability in the computation of post-Newtonian quantities. Seyrich [25] developed Gauss Runge–Kutta implicit canonical symplectic schemes that preserve the structural properties of the system. These schemes ensure long-term numerical stability and accuracy. These algorithms, with their distinct methodologies, contribute to advancing the computation of post-Newtonian Lagrangian quantities, addressing specific aspects such as precision, radiation reaction, spin contributions, stability, and structural preservation.
Regarding the post-Newtonian Hamiltonian, the phase-space expansion method [26,27,28] is a usable algorithm. The Hamiltonian lacks separability and does not possess a coordinate momentum or multiple integrable splitting components. Pihajoki [26] extended the phase space variables by copying the coordinates and momenta. We achieved a Hamiltonian splitting form so that the explicit leapfrog algorithms become available. The permutation map of momentum was designed to suppress the interaction of the original and extended variables. Liu [29] devised a sequential mapping of coordinate and momentum permutations and constructed fourth-order phase-space expansion explicit method compositions of two triple products of the usual second-order leapfrog. These algorithms suffer a clear failure when calculating the chaotic orbits of celestial systems. The interactions between the original variables and the extended one become increasingly strong and show considerably different values, whereas they are supposed to be equivalent. Midpoint and correction maps [30,31] have been proposed to ensure the equivalence of the original variables and the copy one. Recently, we proposed a multi-factor correction map that yields a higher accuracy of the phase-space expansion method without significant computational resource increases [32]. This paper aims to design a multi-factor correction map for post-Newton Hamiltonian and examine its performance.
This article is divided into several sections. In Section 2, we revisit the Lagrangian and Hamiltonian equations of motion for compact binary systems within the post-Newtonian (PN) approximation. Section 3 presents the introduction of the phase-space expansion method and the development of a correction map for the post-Newtonian Hamiltonian. In Section 4, we conduct a comparative analysis of the accuracy of numerical solutions obtained using the implicit midpoint method in the computation of the post-Newtonian Lagrangian and the correction map method for the post-Newtonian Hamiltonian. Finally, in Section 5, we conclude.

2. PN Lagrangian and Hamiltonian in Compact Binary

Let us consider a compact binary system governed by a PN Lagrangian L ( r , v ) up to the m-th order, where r and v represent the position and velocity vectors, respectively. The Euler–Lagrangian equation is given by Equation (1), where the generalized momentum p is defined by Equation (2). The expression for p is given by Equation (2) and represents a nonlinear algebraic equation of v .
d p d t = L r .
Here
p = L v ;
d r d t = v .
We note that Equations (2) and (3) are differential equations, and that ( r , p ) are treated as the integration variables, whereas v is not. However, we can substitute Equation (2) into Equation (3) to obtain the corresponding acceleration equation, given by Equation (4). Here, a N , a 1 P N , a 2 P N , , a m P N correspond to the Newtonian term, the 1st, 2nd post-Newtonian-order term to the m-th post-Newtonian-order contributions for the accelerations.
d v d t = a N + a 1 P N + a 2 P N + a m P N .
When considering only the m-th PN order term in Equation (4), all terms higher than the m-th PN order are truncated. Consequently, Equation (4) does not align with the PN Lagrangian L, and Equations (3) and (4) are treated as incoherent PN equations of motion in the Lagrangian L. However, when utilizing Equations (3) and (4), the variables ( r , v ) can be used as a set of integration variables instead of the variables ( r , p ) . Nevertheless, this approach does not fully maintain constants of motion, such as the energy integral expressed by
E = v · p L .
In this paper, L is the dimensionless post-Newtonian (PN) Lagrangian formulation for compact binaries. The evolution of binaries can be given by the expression:
L = L N + L 1 P N .
L N and L 1 P N denote the non-relativistic and 1PN contributions to the Lagrangian, respectively. For simplicity, higher-order terms are not considered. The non-relativistic part is expressed as:
L N = r ˙ 2 2 + 1 r ,
whereas the 1PN part is given by [8]:
L 1 P N = 1 c 2 1 8 ( 1 3 η ) v 4 + 1 2 r [ ( 3 + η ) v 2 + η r 2 ( r · v ) 2 1 r ] .
Here η = μ / M is the dimensionless mass parameter. The reduced mass, μ , is defined as M 1 M 2 / M = β ( 1 + β ) 2 , and β = M 1 / M 2 is the mass ratio, where M 1 and M 2 represent the masses of the two bodies constituting the binary system and the total mass is denoted as M = M 1 + M 2 . Additionally, c is the speed of light and G represents the constant of gravity given in natural units with c = G = 1 . c is retained in some of the latter equations, and it can be ignored in the actual calculation.
The equations for the evolution of the system can be derived from the Lagrangian formulation. According to Equations (1) and (2), the equation for the evolution of the momentum can be written as:
d p d t = r r 3 1 + 1 c 2 [ 3 η 2 r 2 ( r · v ) 2 + 3 + η 2 v 2 1 r ] + η c 2 r 3 ( r · v ) v ,
where r is the separation vector between the two masses. The non-relativistic and 1PN contributions to this equation are, respectively, expressed in the first and second terms. The expression for the generalized momentum p till 1pN is given in terms of the velocity v as:
p = v + 1 c 2 v 2 2 ( 1 3 η ) v + 1 r [ η r 2 ( r · v ) r + ( 3 + η ) v ] .
The value of momentum can be obtained from Equation (10) once the velocity is known and vice versa. However, velocity needs to be solved iteratively since it cannot be obtained directly from the Lagrangian. The first post-Newtonian relative acceleration equation is used to determine the value of velocity, given by:
d v d t = a N + a 1 P N .
The sub-terms are
a N = r r 3 ,
a 1 P N = 1 r 2 c 2 r r [ ( 1 + 3 η ) v 2 2 r ( 2 + η ) 3 η 2 r 2 ( r · v ) 2 ] 2 r ( 2 η ) ( r · v ) v ,
Here, a N and a 1 P N describe the non-relativistic and 1PN contributions to the acceleration. Using Equations (6) and (10), the energy integral in Equation (5) can be expressed as.
E = v 2 2 1 r + 1 c 2 3 8 ( 1 3 η ) v 4 + 1 2 r [ ( 3 + η ) v 2 + η r 2 ( r · v ) 2 + 1 r ] .
In summary, the Lagrangian formulation of the evolution of binaries provides a mathematical framework for studying their motion. The momentum and velocity of the system can be determined from the equations derived from the Lagrangian, which include non-relativistic and relativistic contributions to the acceleration. The dimensionless PN Lagrangian offers valuable insights for the dynamics system, enabling the study of gravitational wave emission caused by binary systems.
With Equations (3) and (11), we can obtain the numerical solution ( r · v ) by using the fourth-order implicit midpoint method ( I M 4 ).
PN Hamiltonian form H can be derived through Lagrangian L using the Legendre transformation,
H = p · r ˙ L .
Then we obtain the 1PN Hamiltonian,
H = H N + H 1 P N .
In order to compare the effect of higher-order PN terms on the error in the constants of motion, we introduce the 2PN post-Newton Hamiltonian,
H = H N + H 1 P N + H 2 P N .
The expressions for the sub-terms in Hamiltonians (16) and (17) are, respectively, given by
H N = T ( p ) + V ( r ) = p 2 2 1 r ,
H 1 P N = 1 8 ( 3 η 1 ) p 4 1 2 [ ( 3 + η ) p 2 + η r ( r · p ) 2 ] 1 r + 1 2 r 2 ,
H 2 P N = 1 16 ( 1 5 η + 5 η 2 ) p 6 + 1 8 [ ( 5 20 η 3 η 2 ) p 4 2 η 2 r ( r · p ) 2 p 2 3 η 2 r ( r · p ) 4 ] 1 r + 1 2 [ ( 5 + 8 η ) p 2 + 3 η r ( r · p ) 2 ] 1 r 2 1 4 ( 1 + 3 η ) 1 r 3 ,
Due to the disappearance of higher-order terms, H and H are approximately equal to E, and not strictly equivalent. The integrators used in the Hamiltonian H and H will be described in the next section.

3. Phase-Space Expansion Method with a Multi-Factors Correction Map

Since neither the Hamiltonian H nor H can be separated into multiple integrable parts, the symplectic leapfrog method cannot be applied directly to these Hamiltonians unless they are suitably modified to a splitting form. An effective approach to solving this problem is the phase-space expansion method. Pihajoki [26] introduced a new pair of canonical and conjugate variables ( r ˜ , p ˜ ) from the original variables ( r , p ) . This doubles the phase-space variables, ( r , r ˜ , p , p ˜ ) and constructs a new Hamiltonian H ˜ ( r , r ˜ , p , p ˜ ) using two identical Hamiltonians H 1 and H 2 :
H ˜ ( r , r ˜ , p , p ˜ ) = H 1 ( r , p ˜ ) + H 2 ( r ˜ , p ) .
where both H 1 and H 2 should be equal to the original Hamiltonian H. The new Hamiltonian H ˜ already exhibits two integrable components. A conventional second-order leapfrog algorithm can be employed for its integration:
A 2 ( h ) = H 2 ( h 2 ) H 1 ( h ) H 2 ( h 2 ) ,
where h represents the time step, H 1 and H 2 are Hamiltonian operators. The code corresponding to integration (22) from nth to ( n + 1 ) th step is
r n + 1 2 = r n + h 2 p H 2 ( r ˜ n , p n ) p ˜ n + 1 2 = p ˜ n h 2 r ˜ H 2 ( r ˜ n , p n ) r ˜ n + 1 = r ˜ n + h p ˜ H 1 ( r n + 1 2 , p ˜ n + 1 2 ) p n + 1 = p n h r H 1 ( r n + 1 2 , p ˜ n + 1 2 ) r n + 1 = r n + 1 2 + h 2 p H 2 ( r ˜ n + 1 , p n + 1 ) p ˜ n + 1 = p ˜ n + 1 2 h 2 r ˜ H 2 ( r ˜ n + 1 , p n + 1 ) .
It can be seen that the calculation of the numerical solution ( r , p ˜ ) of H 2 requires the numerical solution ( r ˜ , p ) of H 1 , and vice versa, so there is an energy exchange between H 1 and H 2 , and even if the initial conditions are the same, H 1 and H 2 will become unequal in the later evolution unless the errors are constant equal to 0. To submit the accuracy, we construct a fourth-order algorithm using Yoshida’s triplet product
A 4 ( h ) = A 2 ( λ 3 h ) A 2 ( λ 2 h ) A 2 ( λ 1 h ) .
The time coefficients λ 1 , λ 2 , and λ 3 are identical to those presented in [33] and are set to λ 1 = λ 3 = 1 / ( 2 2 1 / 3 ) and λ 2 = 1 2 λ 1 . Algorithm A 4 is utilized to obtain a set of numerical solutions ( r , r ˜ , p , p ˜ ) . It is crucial to note that the original variables ( r , p ) and their counterparts ( r ˜ , p ˜ ) are intended to be identical at each integration step, but in reality, they exhibit discrepancies. The interaction between the solutions ( r , p ˜ ) of H 1 and ( r ˜ , p ) of H 2 leads to their divergence over time. To ensure the equivalence of the original variables and their copy one, Pihajoki [26] proposed a momentum permutation map, whereas Liu [29] proposed coordinate and momentum permutation maps. These applications were successful in several examples, but not in chaotic orbits. However, our previous work [31,32] proposed the manifold corrections map, which effectively overcame the challenge. Unlike that in the original paper [31,32], the correction map for the 1PN Hamiltonian is
M 1 P N = γ 2 , γ 2 , 0 , 0 γ 2 , γ 2 , 0 , 0 0 , 0 , α 2 , α 2 0 , 0 , α 2 , α 2 .
Here, the momentum scaling factor γ and the coordinate scaling factor α are incorporated into M 1 P N and can be obtained by solving the following formulations:
T ( α p + α p ˜ 2 ) = T ˜ ( p , p ˜ ) 2 = T 1 ( p ˜ ) + T 2 ( p ) 2 ,
V ( γ r + γ r ˜ 2 ) + H 1 P N ( γ r + γ r ˜ 2 , α p + α p ˜ 2 ) = V ˜ ( r , r ˜ ) + H ˜ 1 P N ( r , r ˜ , p , p ˜ ) 2 .
Equation (26) provides α = 2 ( p 2 + p ˜ 2 ) ( p + p ˜ ) 2 , and Newton’s method is used to obtain γ from Equation (27). Then, the fourth-order phase-space expansion method with multi-factor correction map for 1PN Hamiltonian H is established as
CM 1 P N ( h ) = M 1 P N A 4 ( h ) .
Similarly, for the 2PN Hamiltonian H , the aforementioned steps are applicable. The correction map M 2 P N for the Hamiltonian H follows the same structure as M 1 P N .
M 2 P N = γ 2 , γ 2 , 0 , 0 γ 2 , γ 2 , 0 , 0 0 , 0 , α 2 , α 2 0 , 0 , α 2 , α 2 .
However, the solution for γ is replaced by the following equation
V ( γ r + γ r ˜ 2 ) + H 1 P N ( γ r + γ r ˜ 2 , α p + α p ˜ 2 ) + H 2 P N ( γ r + γ r ˜ 2 , α p + α p ˜ 2 ) = V ˜ ( r , r ˜ ) + H ˜ 1 P N ( r , r ˜ , p , p ˜ ) + H ˜ 2 P N ( r , r ˜ , p , p ˜ ) 2 .
For convenience, such algorithms are referred to as the correction map method. The correction map method for 2PN Hamiltonian H is set up as
CM 2 P N ( h ) = M 2 P N A 4 ( h ) .
The A 4 algorithm is treated as an explicit symplectic method serving to the new Hamiltonian H ˜ , ensuring effective preservation of the energy of H ˜ , i.e., Δ H ˜ = Δ H 1 + Δ H 2 0 . The error evolution of H 1 and H 2 shows a clear time–axis symmetry, as demonstrated in Figure 1. Specifically, if H 1 calculates more energy than the initial energy, H 2 will calculate less, and vice versa. Taking advantage of this symmetry, we designed a manifold correction mapping approach to optimize the performance of the A 4 algorithm. The M 1 P N and M 2 P N corrections imposed on the solutions of A 4 serve three primary purposes. The A 4 algorithm serves multiple purposes in relation to the new Hamiltonian H ˜ . Firstly, it ensures that H 1 is equal to H 2 to prevent energy discrepancies that could impede the availability of numerical solutions. Secondly, it maintains the constancy of H ˜ after the correction, effectively suppressing the growth of energy errors. Thirdly, it reduces the energy deviation of each subterm of H from half of the corresponding subterm of H ˜ through the correction process. As the A 4 algorithm is an explicit symplectic method serving the new Hamiltonian H ˜ , it accurately calculates the total energy as well as the energy of each individual subterm in H ˜ ; thus, the algorithm altering these energies is not desirable, as it may weaken the algorithm’s stability and precision.
In Section 4, we will set initial values and perform numerical simulations of post-Newtonian Lagrangian and post-Newtonian Hamiltonian to compare the differences between the algorithms in terms of maintaining the constants of motion.

4. Numerical Simulation

This section showcases the outcomes of our numerical simulations, where we compare the post-Newtonian Lagrangian and post-Newtonian Hamiltonian algorithms in maintaining the constants of motion. To this end, we set initial values for a specific orbit, named orbit 1, with initial conditions ( β ; r , v ) = ( 5 4 ; 10 , 0 , 0 , 0 , 0.52 , 0 ) . The initial value of the momentum p in the post-Newtonian Hamiltonian is obtained from Equation (10). We use the fourth-order implicit midpoint method ( I M 4 ) to calculate the 1PN Lagrangian, whereas the algorithms CM 1 P N and CM 2 P N are used to calculate the Hamiltonian H and H , respectively. We take a fixed step size of h = 1 and plot the energy errors in Figure 2a,b. We observe that the CM 1 P N algorithm designed for the Hamiltonian H has significantly better accuracy in terms of energy error compared to I M 4 . However, the accuracy of CM 1 P N drops considerably in the H error behavior, as expected due to the vanishing of the 2PN term, whereas the error in CM 2 P N is at an order of 10 8 . Finally, Figure 2c shows that CM 2 P N performs the best in terms of accuracy and long-term stability compared to CM 1 P N and I M 4 .
Aside from ensuring energy conservation, we also track the preservation of orbital angular momentum Ł = r × p = [ 1 + 1 c 2 ( 1 3 η 2 v 2 + 3 + η r ) ] r × v , and examine its error as another performance metric for the algorithms. Figure 3 depicts the angular momentum errors Δ Ł = Ł Ł 0 with Ł = |Ł|, where Ł 0 represents the initial value. We deduce that the performance of the I M 4 algorithm in terms of angular momentum error is similar to its performance in energy error, with the worst accuracy but good long-term stability. Conversely, both CM 1 P N and CM 2 P N exhibit exceptional performance in terms of angular momentum error, with very little difference between them and significantly superior to I M 4 . However, there is a noticeable error growth in CM 1 P N and CM 2 P N . Our simulations show that the algorithm CM 1 P N in the 1PN Hamiltonian has a significant advantage over the I M 4 algorithm in maintaining the conservation of orbital angular momentum and a small accuracy advantage in maintaining energy integrals. The algorithm CM 2 P N in the 2PN Hamiltonian also has a significant advantage in maintaining angular momentum, while being comparable to CM 1 P N , with some improvement in the accuracy of the energy.
To validate the aforementioned conclusion, additional numerical simulations will be conducted in a different orbit, referred to as orbit 2. The initial conditions for orbit 2 are set as ( β ; r , v ) = ( 5 4 ; 10 , 0 , 0 , 0 , 0.52 , 0 ) . In Figure 4, three categories of energy error in orbit 2 will be depicted as follows: (a) Energy error analysis of orbit 2 for I M 4 , CM 1 P N , and CM 2 P N will be presented in Figure 4a. (b) Figure 4b will display the energy error analysis of orbit 2 specifically for I M 4 . (c) The energy error analysis of orbit 2, focusing on I M 4 , will be illustrated in Figure 4c.
It is evident from the figures that the performance of I M 4 , CM 1 P N , and CM 2 P N in orbit 2 closely resembles that of orbit 1. I M 4 continues to exhibit the highest error, CM 1 P N demonstrates a widening gap with I M 4 , and CM 2 P N remains the most advanced in terms of accuracy.
Turning to the angular momentum errors depicted in Figure 5, it is observed that there is no significant improvement for I M 4 , which still exhibits considerable deviation compared to the first-order post-Newtonian approximation, CM 1 P N . Furthermore, the inclusion of the second-order post-Newtonian term in CM 2 P N does not contribute significantly to reducing the angular momentum error.
Summarizing the findings from the numerical simulations conducted for both orbit 1 and orbit 2, we can conclude that CM 1 P N and CM 2 P N perform better for the calculation of the post-Newtonian approximation to the Hamiltonian. They exhibit a slight advantage in terms of energy error while demonstrating a notably superior accuracy in the calculation of angular momentum.

5. Summary

The exact equations of motion for a post-Newtonian Lagrangian formalism are the Euler–Lagrange equations, which consist of a coherent Lagrangian without any truncated terms. However, when the post-Newtonian Lagrangian form of general relativity maintains only a certain post-Newtonian order, it is referred to as the incoherent Lagrangian, with higher-order terms of the acceleration truncated. Incoherent Lagrangian can be numerically simulated using the Runge–Kutta method and implicit algorithms. Therefore, in the incoherent Lagrangian, motion constants such as energy integrals are only approximately conserved. The retention of the Hamiltonian in a certain post-Newtonian (PN) order leads to a high-order truncation problem. In addition, if the Hamiltonian is separable, symplectic algorithms can be used, which provide excellent performance. For the case where the Hamiltonian is inseparable, symplectic-like algorithms such as the phase-space expansion method with correction map, namely correction map method, can be used. The phase-space expansion method with correction map is referred to as the correction map method, which utilizes the symmetry of energy errors in H 1 and H 2 in the new Hamiltonian H ˜ to improve the accuracy and stability of the algorithm.
A comparison was made between the performance of the implicit midpoint method in the incoherent Lagrangian and the correction map method in the PN Hamiltonian. Under the 1PN Newtonian approximation, the correction map method performed better in terms of energy error, exhibiting higher accuracy and comparable stability. On the other hand, with regards to angular momentum error, the correction map method was significantly higher, reaching an order of 10 11 , whereas the implicit midpoint method was only 10 1 . Similarly, under the 2PN Newtonian Hamiltonian, the manifold correction mapping method further improved the accuracy of energy error, but there was no noticeable impact on the angular momentum error.
In conclusion, we compared the implicit midpoint methods for solving the equations of motion in post-Newtonian Lagrangians and the correction map method for PN Hamiltonians and investigate the extent to which both methods can uphold invariance of the motion’s constants, such as energy conservation and angular momentum preservation. Ultimately, the results of numerical simulations demonstrate the superior performance of the correction map method, particularly with respect to angular momentum conservation. Compared to incoherent Lagrangian, we recommend using the manifold correction map method for the Hamiltonian of compact binaries as a numerical tool.

Author Contributions

Conceptualization, W.L., H.-H.Z., J.L. and J.F.; methodology, J.L.; software, J.L.; validation, W.L., H.-H.Z., J.L. and J.F.; formal analysis, J.L.; investigation, H.-H.Z.; resources, W.L.; data curation, H.-H.Z.; writing—original draft preparation, J.L. and J.F.; writing—review and editing, W.L., H.-H.Z., J.L. and J.F.; visualization, J.L.; supervision, W.L.; project administration, H.-H.Z.; funding acquisition, W.L., H.-H.Z. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the National Natural Science Foundation of China (NSFC) through Grant Nos. 12203108, 12275367, and 12073089. Additional support was provided by the Fundamental Research Funds for the Central Universities and the Sun Yat-sen University Science Foundation. The Fundamental Research Funds for the Central Universities play a significant role in supporting fundamental research activities within the academic institutions, fostering a conducive environment for intellectual exploration and discovery. The authors acknowledge these funding sources for their valuable contributions to the completion of this study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
I M 4 Fourth-order implicit symplectic integrator
CM 1 P N Correction map method for 1PN Hamiltonian
CM 2 P N Correction map method for 2PN Hamiltonian

References

  1. Wu, X.; Wang, Y.; Sun, W.; Liu, F.Y.; Han, W.B. Explicit Symplectic Methods in Black Hole Spacetimes. Astrophys. J. 2022, 940, 166. [Google Scholar] [CrossRef]
  2. Zhou, N.Y.; Zhang, H.X.; Liu, W.F.; Wu, X. A Note on the Construction of Explicit Symplectic Integrators for Schwarzschild Spacetimes. Astrophys. J. 2022, 927, 160, Erratum in Astrophys. J. 2023, 947, 94. [Google Scholar] [CrossRef]
  3. Wang, Y.; Sun, W.; Liu, F.Y.; Wu, X. Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes. Astrophys. J. 2021, 907, 66. [Google Scholar] [CrossRef]
  4. Wang, Y.; Sun, W.; Liu, F.Y.; Wu, X. Construction of Explicit Symplectic Integrators in General Relativity. II. Reissner-Nordström Black Holes. Astrophys. J. 2021, 909, 22. [Google Scholar] [CrossRef]
  5. Wu, X.; Wang, Y.; Sun, W.; Liu, F.Y. Construction of Explicit Symplectic Integrators in General Relativity. IV. Kerr Black Holes. Astrophys. J. 2021, 914, 63. [Google Scholar] [CrossRef]
  6. Sun, W.; Wang, Y.; Liu, F.Y.; Wu, X. Applying explicit symplectic integrator to study chaos of charged particles around magnetized Kerr black hole. Eur. Phys. J. C 2021, 81, 785. [Google Scholar] [CrossRef]
  7. Pan, G.F.; Wu, X.; Liang, E.W. Extended phase-space symplectic-like integrators for coherent post-Newtonian Euler–Lagrange equations. Phys. Rev. D 2021, 104, 044055. [Google Scholar] [CrossRef]
  8. Blanchet, L.; Iyer, B.R. Third post-Newtonian dynamics of compact binaries: Equations of motion in the centre-of-mass frame. Class. Quantum Gravity 2003, 20, 755. [Google Scholar] [CrossRef]
  9. Tanay, S.; Stein, L.C.; Ghersi, J.T.G. Integrability of eccentric, spinning black hole binaries up to second post-Newtonian order. Phys. Rev. D 2021, 103, 064066. [Google Scholar] [CrossRef]
  10. Arun, K.G.; Blanchet, L.; Iyer, B.; Qusailah, M.S. Inspiralling compact binaries in quasi-elliptical orbits: The complete third post-Newtonian energy flux. Phys. Rev. D 2008, 77, 064035. [Google Scholar] [CrossRef] [Green Version]
  11. Tessmer, M.; Schäfer, G. Eccentric motion of spinning compact binaries. Phys. Rev. D 2014, 89, 104055. [Google Scholar] [CrossRef] [Green Version]
  12. Hinder, I.; Kidder, L.; Pfeiffer, P. Eccentric binary black hole inspiral-merger-ringdown gravitational waveform model from numerical relativity and post-Newtonian theory. Phys. Rev. D 2018, 98, 044015. [Google Scholar] [CrossRef] [Green Version]
  13. Chattaraj, A.; RoyChowdhury, T.; Divyajyoti; Mishra, C.; Gupta, A. High accuracy post-Newtonian and numerical relativity comparisons involving higher modes for eccentric binary black holes and a dominant mode eccentric inspiral-merger-ringdown model. Phys. Rev. D 2022, 106, 124008. [Google Scholar] [CrossRef]
  14. Chowdhury, S.R.; Khlopov, M. An eccentric binary black hole in post-Newtonian theory. Symmetry 2022, 14, 510. [Google Scholar] [CrossRef]
  15. Wu, X.; Mei, L.; Huang, G.; Liu, S. Analytical and numerical studies on differences between lagrangian and hamiltonian approaches at the same post-newtonian order. Phys. Rev. D 2015, 91, 024042. [Google Scholar] [CrossRef]
  16. Wu, X.; Huang, G. Ruling out chaos in comparable mass compact binary systems with one body spinning. Mon. Not. R. Astron. Soc 2015, 452, 3167–3178. [Google Scholar] [CrossRef] [Green Version]
  17. Huang, L.; Wu, X.; Ma, D.Z. Second post-Newtonian Lagrangian dynamics of spinning compact binaries. Eur. Phys. J. C 2016, 76, 488. [Google Scholar] [CrossRef] [Green Version]
  18. Li, D.; Wang, Y.; Deng, C.; Wu, X. Coherent post-Newtonian Lagrangian equations of motion. Eur. Phys. J. Plus 2021, 135, 390. [Google Scholar] [CrossRef]
  19. Li, D.; Wu, X.; Liang, E. Effect of the Quadrupole–Monopole Interaction on Chaos in Compact Binaries. Ann. Der Phys. 2019, 531, 1900136. [Google Scholar] [CrossRef]
  20. Li, D.; Wu, X.; Liang, E. Figure-eight orbits in three post-Newtonian formulations of triple black holes. Phys. Rev. D 2021, 104, 044039. [Google Scholar] [CrossRef]
  21. Zhong, S.Y.; Wu, X. Manifold corrections on spinning compact binaries. Phys. Rev. D 2010, 81, 104037. [Google Scholar] [CrossRef]
  22. Tsang, D.; Galley, C.R.; Stein, L.C. “Slimplectic” integrators: Variational integrators for general nonconservative systems. Astrophys. J. Lett. 2015, 809, L9. [Google Scholar] [CrossRef] [Green Version]
  23. Lubich, C.; Walther, B.; Bruegmann, B. Symplectic integration of post-Newtonian equations of motion with spin. Phys. Rev. D 2010, 81, 104025. [Google Scholar] [CrossRef] [Green Version]
  24. Zhong, S.Y.; Wu, X.; Liu, S.Q.; Deng, X.F. Global symplectic structure-preserving integrators for spinning compact binaries. Phys. Rev. D 2010, 82, 124040. [Google Scholar] [CrossRef]
  25. Seyrich, J. Gauss collocation methods for efficient structure preserving integration of post-Newtonian equations of motio. Phys. Rev. D 2013, 87, 084064. [Google Scholar] [CrossRef] [Green Version]
  26. Pihajoki, P. Explicit methods in extended phase space for inseparable Hamiltonian problems. Celest. Mech. Dyn. Astron. 2015, 121, 211–231. [Google Scholar] [CrossRef] [Green Version]
  27. Li, D.; Wu, X. Modification of logarithmic Hamiltonians and application of explicit symplectic-like integrators. Mon. Not. R. Astron. Soc 2017, 469, 3031–3041. [Google Scholar] [CrossRef]
  28. Li, D.; Wu, X. Chaotic motion of neutral and charged particles in a magnetized Ernst-Schwarzschild spacetime. Eur. Phys. J. Plus 2019, 134, 96. [Google Scholar] [CrossRef] [Green Version]
  29. Liu, L.; Wu, X.; Huang, G.; Liu, F.Y. Higher order explicit symmetric integrators for inseparable forms of coordinates and momenta. Mon. Not. R. Astron. Soc 2016, 459, 1968. [Google Scholar] [CrossRef] [Green Version]
  30. Luo, J.; Wu, X.; Huang, G.; Liu, F.Y. Explicit symplectic-like integrators with midpoint permutations for spinning compact binaries. Astrophys. J. 2017, 834, 64. [Google Scholar] [CrossRef]
  31. Luo, J.; Lin, W.P.; Yang, L. Explicit symplectic-like integration with corrected map for inseparable Hamiltonian. Mon. Not. R. Astron. Soc 2021, 501, 1511. [Google Scholar] [CrossRef]
  32. Luo, J.; Feng, J.; Zhang, H.H.; Lin, W.P. Performance of different correction maps in the extended phase-space method for spinning compact binaries. Mon. Not. R. Astron. Soc 2022, 518, 6132. [Google Scholar] [CrossRef]
  33. Yoshida, H. Construction of higher order symplectic integrators. Phys. Lett. A 1990, 150, 262. [Google Scholar] [CrossRef]
Figure 1. The energy errors Δ E of the Hamiltonian in Equation (16), as computed using the A 4 algorithm after the extended phase space, can be expressed as Δ H ˜ = H ˜ 2 H ( 0 ) = Δ H 1 + Δ H 2 , whereas Δ H i = H i ( t ) H ( 0 ) and H i ( t ) represent the value of the Hamiltonian H i at time t. H ( 0 ) denotes the initial value of the Hamiltonian in Equation (16). Time–axis symmetry exists between Δ H 1 (red dot) and Δ H 2 (blue dash).
Figure 1. The energy errors Δ E of the Hamiltonian in Equation (16), as computed using the A 4 algorithm after the extended phase space, can be expressed as Δ H ˜ = H ˜ 2 H ( 0 ) = Δ H 1 + Δ H 2 , whereas Δ H i = H i ( t ) H ( 0 ) and H i ( t ) represent the value of the Hamiltonian H i at time t. H ( 0 ) denotes the initial value of the Hamiltonian in Equation (16). Time–axis symmetry exists between Δ H 1 (red dot) and Δ H 2 (blue dash).
Symmetry 15 01401 g001
Figure 2. Different energy errors ( Δ H , Δ H , Δ E ) in orbit 1. (a) The energy error of H, denoted as Δ H = | H ( t ) H ( 0 ) | , where H ( t ) represents the value of the Hamiltonian H at time t, and H ( 0 ) is the initial value. (b) The energy error of H , denoted as Δ H = | H ( t ) H ( 0 ) | . (c) The energy error of E, denoted as Δ E = | E ( t ) E ( 0 ) | . The algorithm IM 4 is drawn with a black line, whereas CM 1 P N and CM 2 P N are drawn with red and blue lines, respectively.
Figure 2. Different energy errors ( Δ H , Δ H , Δ E ) in orbit 1. (a) The energy error of H, denoted as Δ H = | H ( t ) H ( 0 ) | , where H ( t ) represents the value of the Hamiltonian H at time t, and H ( 0 ) is the initial value. (b) The energy error of H , denoted as Δ H = | H ( t ) H ( 0 ) | . (c) The energy error of E, denoted as Δ E = | E ( t ) E ( 0 ) | . The algorithm IM 4 is drawn with a black line, whereas CM 1 P N and CM 2 P N are drawn with red and blue lines, respectively.
Symmetry 15 01401 g002
Figure 3. The angular momentum errors in orbit 2, denoted as Δ Ł = Ł Ł 0 , where Ł = |Ł| and Ł is calculated using three different methods: I M 4 (represented by black), CM 1 P N (represented by red), and CM 2 P N (represented by blue).
Figure 3. The angular momentum errors in orbit 2, denoted as Δ Ł = Ł Ł 0 , where Ł = |Ł| and Ł is calculated using three different methods: I M 4 (represented by black), CM 1 P N (represented by red), and CM 2 P N (represented by blue).
Symmetry 15 01401 g003
Figure 4. Different energy errors ( Δ H , Δ H , Δ E ) in orbit 2. (a) The energy error of H, represented as Δ H , is calculated as the absolute difference between the value of the Hamiltonian H at time t ( H ( t ) ) and its initial value ( H ( 0 ) ). (b) The energy error of H , denoted as Δ H , is determined as the absolute difference between the value of the Hamiltonian H at time t ( H ( t ) ) and its initial value ( H ( 0 ) ). (c) The energy error of E, denoted as Δ E , is computed as the absolute difference between the value of E at time t ( E ( t ) ) and its initial value ( E ( 0 ) ). The algorithm IM4 is represented by a solid black line, whereas C1PN and C2PN are indicated by dashed red and blue lines, respectively. The performance of each algorithm in orbit 2 is similar to that of orbit 1.
Figure 4. Different energy errors ( Δ H , Δ H , Δ E ) in orbit 2. (a) The energy error of H, represented as Δ H , is calculated as the absolute difference between the value of the Hamiltonian H at time t ( H ( t ) ) and its initial value ( H ( 0 ) ). (b) The energy error of H , denoted as Δ H , is determined as the absolute difference between the value of the Hamiltonian H at time t ( H ( t ) ) and its initial value ( H ( 0 ) ). (c) The energy error of E, denoted as Δ E , is computed as the absolute difference between the value of E at time t ( E ( t ) ) and its initial value ( E ( 0 ) ). The algorithm IM4 is represented by a solid black line, whereas C1PN and C2PN are indicated by dashed red and blue lines, respectively. The performance of each algorithm in orbit 2 is similar to that of orbit 1.
Symmetry 15 01401 g004
Figure 5. The errors of the angular momentum Δ Ł = Ł Ł 0 in orbit 2, here Ł = |Ł|, Ł are calculated by I M 4 (black), CM 1 P N (red), CM 2 P N (blue).
Figure 5. The errors of the angular momentum Δ Ł = Ł Ł 0 in orbit 2, here Ł = |Ł|, Ł are calculated by I M 4 (black), CM 1 P N (red), CM 2 P N (blue).
Symmetry 15 01401 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, J.; Feng, J.; Zhang, H.-H.; Lin, W. Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian. Symmetry 2023, 15, 1401. https://doi.org/10.3390/sym15071401

AMA Style

Luo J, Feng J, Zhang H-H, Lin W. Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian. Symmetry. 2023; 15(7):1401. https://doi.org/10.3390/sym15071401

Chicago/Turabian Style

Luo, Junjie, Jie Feng, Hong-Hao Zhang, and Weipeng Lin. 2023. "Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian" Symmetry 15, no. 7: 1401. https://doi.org/10.3390/sym15071401

APA Style

Luo, J., Feng, J., Zhang, H. -H., & Lin, W. (2023). Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian. Symmetry, 15(7), 1401. https://doi.org/10.3390/sym15071401

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop