*Project supported by the National Natural Science Foundation of China (Grant Nos. 11201169, 11271195, and 41231173) and the Project of Graduate Education Innovation of Jiangsu Province, China (Grant No. CXLX13 366).
A local energy conservation law is proposed for the Klein–Gordon–Schrödinger equations, which is held in any local time–space region. The local property is independent of the boundary condition and more essential than the global energy conservation law. To develop a numerical method preserving the intrinsic properties as much as possible, we propose a local energy-preserving (LEP) scheme for the equations. The merit of the proposed scheme is that the local energy conservation law can hold exactly in any time–space region. With the periodic boundary conditions, the scheme also possesses the discrete change and global energy conservation laws. A nonlinear analysis shows that the LEP scheme converges to the exact solutions with order 𝒪( τ2 + h2). The theoretical properties are verified by numerical experiments.
The Klein– Gordon– Schrö dinger (KGS) equations
describe the interaction between conservative complex neutron field and neutral meson Yukawa in quantum field theory, where ψ (x, t) represents a complex scalar nucleon field, ϕ (x, t) is a real scalar meson field, and
there exist a charge conservation law
and a global energy conservation law (GECL)
Massive works[1– 3] with respect to mathematical and theoretical analyses of the KGS equations have been done in the past decades. Studying its numerical approximation has recently become a hot topic. Zhang derived a charge- and energy-preserving finite difference (C-C) scheme for the KGS equations. During the last decade, many researchers have paid attention to geometrical numerical integrators. Kong et al. noted that the KGS equations have a multi-symplectic structure and then investigated the multi-symplectic midpoint (M-M) scheme. The explicit multi-symplectic schemes were proposed in Refs.  and ; while the implicit symplectic Fourier pseudo-spectral schemes were given in Refs.  and . Hong et al. also proposed three schemes, which are called M-C, T-C, and T-T schemes, for the KGS equations. Some other methods can be found in Refs. – .
As is well known, the GECL (4) depends on the boundary conditions, such as periodic or homogeneous boundary conditions. However, in practice, the complex boundary conditions lead to the invalidity of the traditional GECL algorithms. Actually, there exists a local energy conservation law, which is more essential than the global energy conservation law. Let ψ (x, t) = p(x, t) + iq(x, t), where p(x, t) and q(x, t) are real-valued functions. By further introducing some variables, the KGS system (1) can be reformed into a combination of ordinary differential equations (ODEs)
Proposition 1 The system (5) possesses a local energy conservation law (LECL)
Proof In Eq. (5a), multiplying the first and the second lines by pt and qt, respectively, we have
Next, multiplying the first line of Eq. (5b) by ϕ t gives
The LECL (6) is a local property. In other words, it is independent of the boundary conditions. Thus, the LECL is more essential than the GECL in physics. Naturally, when constructing algorithms for system (5), we expect that the LECL can still hold. Hereafter, we call these numerical methods the local energy-preserving (LEP) algorithms. For the concept of LEP algorithm and its applications, please refer to Refs. – . The aim of this paper is to construct a LEP algorithm for the KGS equations, and then analyze its convergence. Finally, we will check the performance of the proposed scheme.
As usual, we introduce some notations: xj = a + (j – 1)h, j = 1, 2, … , N + 1; tk = kτ , k = 0, 1, 2, … K, where h = (b – a)/N and τ are the spatial length and the temporal step span, respectively. Let
Taking f = g gives
Similarly, we can obtain a series of analogous discrete Leibnite rules in the time direction.
LEP scheme We discrete the ODEs (5) as follows:
Eliminating the auxiliary variables m, n, v, α , f, g, and w, we obtain a scheme
or an equivalent scheme
Now, we analyze the local and global conservative properties of the proposed LEP scheme.
Theorem 1 The LEP scheme is a local energy conservative algorithm, which admits the discrete LECL
Proof In Eq. (11a), multiplying the first and the second lines by
Equation (15) can be reduced as follows:
Multiplying the first line of Eq. (11b) by
Corollary 1 For the initial and periodic boundary conditions (2), the LEP-I scheme preserves the discrete GECL
To prove it, we just need to sum the LECL (14) over all space index j and utilize the periodic boundary conditions.
In general, as a quadratic invariant, the charge conservation law (3) plays an important role in quantum physics. Therefore, it is necessary to discuss whether it can be captured.
Theorem 2 For the initial and periodic boundary conditions (2), the LEP scheme conserves the charge conservation law exactly,
namely, the scheme is unconditionally stable with respect to the initial values.
Proof Multiplying Eq. (13a) by
The first term of Eq. (20) becomes
where Im stands for the imaginary part. The second term
is a real-valued function, where Re stands for the real part. The last term
By using the discrete Sobolev inequality and the Gronwall inequality,  we have the following convergence results.
Since the proof is very trivial, the details are not presented here. It should be pointed out that, in order to make the convergence analysis, we must prove ∥ ψ k ∥ ∞ ≤ C and ∥ ϕ k ∥ ∞ ≤ C, and analyze the truncation error of the LEP scheme first.
In this section, some numerical experiments are carried out to show the performance of the LEP scheme. The performance of the proposed method will be exhibited in the following aspects: (i) the accuracy of single soliton solutions; (ii) the performance in preserving the local energy, charge, and global energy conservation laws. The accuracy of migration of a soliton at tk = kτ is measured by
In order to show the variation of the local energy at the k-th time level against time, we define
Example 1 (soliton) The KGS equations have the following soliton solution:
Here, ν represents the propagating velocity of the wave and x0 is the initial phase. We take the initial conditions ψ 0(x) = ψ (x, 0, ν , 0), ϕ 0(x) = ϕ (x, 0, ν , 0), ϕ 1(x) = ϕ t(x, 0, ν , 0), and the periodic boundary conditions ψ (– 40, t) = ψ (40, t), ϕ (– 40, t) = ϕ (40, t), ϕ t(– 40, t) = ϕ t(40, t).
We choose ν = 0.1, τ = 0.01, and h = 0.1. From Fig. 1, one can see that the errors in solution ψ grow nearly linearly while the ones in solution ϕ grow oscillatory. One can also find that figure 1 is very similar to that in Ref. . The local energy and the errors in the global invariants are listed in Table 1. It is clear that the local energy 𝓔 local is always in the scale of 10− 14 throughout computations. Therefore, the local energy is conserved exactly. Since the errors in charge and global energy oscillate near zero in the scale of 10− 15, the two global invariants are also captured.
Next, we compare the numerical performance of the LEP scheme with that of the existing methods.[4, 5, 9] We implement all methods with ν = 0.3, various temporal and spatial step sizes. Table 2 reports the errors in the solutions, and the changes in charge and global energy at T = 1. From the table, we find that the solutions of LEP are better than those of C-C, M-M, T-C, and T-T schemes, while a little worse than the M-C scheme. The errors in invariants show that our method conserves the charge and the global energy much better than the others.
Example 2 (solitons collision) In the following simulations, we study head-on collisions of two solitary waves. The initial conditions are chosen as ψ 0 = ψ (x, 0, v1, x1) + ψ (x, 0, v2, x2), ϕ 0 = ϕ (x, 0, v1, x1) + ϕ (x, 0, v2, x2), and ϕ 1 = ϕ t(x, 0, v1, x1) + ϕ t(x, 0, v2, x2). We consider the simulation over the region [– 40, 40] × [0, 80] with parametersτ = 0.01 and h = 0.1.
First, we investigate the collision of symmetric solitons, i.e., ν 1 = 0.4, x1 = – 20 and ν 2 = – 0.4, x2 = 20. Figure 2 displays the evolution of ψ and ϕ . The two solitons collide at about T = 40, and then seem to be trapped. Table 3 displays the local energy, and the errors in charge and global energy against time. The results coincide with the theoretical analysis.
Second, we consider the non-symmetric solitions collision with ν 1 = 0.4, x1 = – 20 and ν 2 = – 0.2, x2 = 15. The simulation results are illustrated in Fig. 3. It is clear that the soliton with the larger amplitude becomes larger and the one with the smaller amplitude becomes smaller after the collision. The local energy 𝓔 local, and the errors in charge and global energy are presented in Table 4. The results are all within the roundoff error of the machine. Again the proposed LEP scheme conserves the invariants exactly.
In this work, we propose a LEP scheme to simulate the KGS equations. The local property is independent of the boundary conditions. Our method preserves the discrete LECL in any time– space region. If the boundary conditions are suitable, the LEP method will be the GEP method. Therefore, the method in the present paper extends the applying scope of the traditional GEP algorithm. Numerical results indicate that the present scheme can provide accurate solutions in long-term computations, and also show excellent performance in conserving the local energy, global energy, and charge invariants. From Eq. (5), we can also derive local momentum and global momentum conservation laws for the KGS equations. As we know, they are also important invariants in physics. Similarly, we could consider local momentum-preserving schemes for the KGS equations.