Acquiring the Symplectic Operator Based on Pure Mathematical Derivation Then Verifying It in the Intrinsic Problem of Nanodevices

: The symplectic algorithm can maintain the symplectic structure and intrinsic properties of the system, its cumulative error is small and suitable for multi-step calculation. At present, the widely accepted symplectic operators are obtained by solving the Hamilton equation based on artiﬁcial deﬁnitions and assumptions in advance. There are inevitable dispersion errors. We solve the equation by pure mathematical derivation without any artiﬁcial limitations and assumptions. The way to accurately obtain high-precision symplectic operators greatly reduces the dispersion error from the beginning. The numerical solution of the one-dimensional Schrödinger equation for describing the intrinsic problem of nanodevices is used as an application environment to compare the total energy distribution of the particle wave function in the box, thus verifying the properties of the Symplectic Operator based on Pure Mathematical Derivation by comparing with Finite-Di ﬀ erence Time-Domain (FDTD) and the widely accepted symplectic operator.


Introduction
The Hamilton system includes pharmacology, structural biology, semiconductors, superconductivity, plasma, celestial mechanics, materials, and partial differential equations. In the study of the Hamiltonian system, there are three efficient methods namely: Finite-Difference Time-Domain (FDTD), Runge-Kutta algorithm, and symplectic algorithms. The Runge-Kutta algorithm is the most widely used method, but its order has certain limitations, and the stability decreases with the increase of the order [1]. Although FDTD is also widely used, the calculation accuracy is low [2] and the time-space stability factor must satisfy a specific condition. Moreover, its numerical stability must be affected by the changes of dielectric constant in different media [3]. Compared with the traditional non-symplectic algorithm, the symplectic algorithm not only has higher precision, but also maintains the symplectic structure of the system without the infinite accumulation of energy errors. Those characteristics are greatly beneficial for solving equations.
There are many ways to achieve the symplectic algorithm. The FDTD method starts from the Maxwell equations that summarize the basic laws of macroscopic electromagnetic fields, and makes full use of the accuracy of finite difference to achieve accurate solution of electromagnetic problems. The FDTD method is very convenient for electromagnetic modeling of complex targets and has been widely used in solving electromagnetic problems such as electromagnetic propagation, radiation, and scattering. We start from Maxwell's equations and get it in Cartesian coordinates.
These six coupled partial differential equations are the basis of the FDTD algorithm, where E and H represent the electric field strength and the magnetic field strength, respectively. In the FDTD method, the time step ∆t and the spatial step size ∆x ∆y ∆z are not independent of each other, and their values must conform to a certain relationship.
It can be inferred that in order to ensure numerical stability, each discrete parameter should satisfy This is the numerical stability condition of the FDTD algorithm, which gives the relationship that should be satisfied between space and time interval.
Ruth [4] proposed a method for calculating the third-order symplectic algorithm equation in 1983. Ruth considered a system of differential equations governed by the Hamiltonian, He also developed an explicit third-order symplectic map and to indicate the method for higher order maps. He let the approximate nth order symplectic map be denoted by (x, p) = m n (t) x 0 , p 0 , where t is the time step (assumed small) and n is the order of the map Then he demonstrated a method for finding M n (t). This method is the foundation for the development of the symplectic algorithm. The dispersion characteristics are not suitable to higher-order computational equations. Iwatsu [5] and Cui et al. [6] considered the deficiencies of the Runge-Kutta method by limiting the range of a particular symplectic operator to determine other parameters, this method still cannot efficiently reduce the error to get accurate results. In order to achieve the symplectic method for stochastic Hamiltonian systems with multiplicative noise, Milstein et al. [7] proposed a new class of completely implicit methods for stochastic systems. The increment of the Wiener process in these completely implicit schemes is replaced by some truncated random variables. Mcmahon et al. [8] constructed different integrators based on the selection of discrete Lagrangian for approximation actions, which is beneficial for construction by combining discrete Lagrangian quantities in the explicit symplectic Runge-Kutta method. Xiao et al. [9] developed an explicit high-order, non-classical symplectic algorithm for an ideal two-fluid system. In this method, the fluid is discretized into particles, and the electromagnetic field and internal energy are regarded as discrete differential form fields on a fixed grid. It is suitable for large-scale simulation of large-scale physical problems and requires long-term fidelity and accuracy. Wen [10] defined a hidden symplectic format by using the finite term of the series to approximate the original generation functional in the Hamiltonian system, but it can only restore the physical process more completely and does not explore the huge application of the symplectic algorithm in application research. Qing et al. [11] applied the symplectic algorithm to the analysis of cracks without considering the stable hybrid scheme used in non-traditional solutions or classical hybrid models, nor extending the method to the 2D domain.
In the study of the intrinsic problem of nanodevices, Liang et al. [12] used the variational method to solve the Schrödinger equation, Feng [13] studied the characteristics of He~+ high-order harmonic radiation under butterfly metal nanostructures by numerically solving Schrödinger equation theory. These methods have strict requirements for the accuracy of the solution. Cheng et al. [14] considered the surface polarization effect for nanocrystal quantum dots in discrete background media and applied the perturbation method to solve the Schrödinger equation of excitons. This requires a smaller error range for the analytical solution. Liu et al. [15] conducted a theoretical study on the relationship between the semiconductor surface barrier and the surface energy level and derived the eigenequation of the surface electron wave function energy using the Schrödinger equation and matrix theory. Given the nonlinear relationship between the surface barrier and the surface energy level, it is necessary to use the exact analytical solution to give the simulation curve of the surface energy level solution set under different surface barriers.
How to achieve accurate high-order symplectic operators with minimum error is a difficult problem. Based on the above factors, we propose a method to acquire the symplectic operator based on pure mathematical derivation. We propose a method for accurately calculating symplectic operators by using the correspondence between coefficients of different series and order. This method establishes the corresponding equations and solves them by purely mathematical derivation, using the same unknown symplectic operator to represent other coefficients. The value of the error is minimized to obtain the exact value of the symplectic operator, and the symplectic operator with the best numerical dispersion characteristics is calculated.

Calculate Symplectic Operators
The general expression of the Hamilton equation is [16] where J = 0 n I n −I n 0 n is the standard symplectic matrix, 0 n , I n are n order 0 matrix and n order unit matrix, H = H(p 1 , p 2 , · · · , p n , q 1 , q 2 , · · · , q n ) is the Hamilton function independent of time t, Z is the transposed matrix of H, and the differential operator is introduced.
In the symplectic propagator technique, we use the non-dissipative m-level p-order explicit symplectic integral to approximate the time evolution matrix of the Maxwell's curl equation. where d l and c l are symplectic operator or symplectic propagator coefficients, and each step of time exp(c l ∆ t U) and exp(d l ∆ t V) is a basic symplectic transformation.
In order to determine the symplectic operator, we perform Taylor series expansion on the 0 point, then there is Comparing the coefficients in the LHS and RHS equations, we can get different order symplectic operators, where the second-order RHS coefficients are expressed as: Define the error function as error = 2 e 2 3,1 + e 2 3,2 + e 2 3,3 + . . . + e 2 3,5 + e 2 3,6 e 2 m,n represents the coefficient term of the corresponding symplectic operator.
In Figure 1, c 1 , c 2 , and d 2 are represented with d 1 . It can be seen that c 1 , c 2 , and d 2 vary linearly with the change of d 1 . Figure 2 shows the variation of the error function with d 1 . It can be seen that the minimum value of the error function is in the range of d 1 .    The third-order RHS coefficients are After exact calculation, the relationship between each symplectic operator and d 1 is The relationship between each symplectic operator and d 1 is shown in Figure 3.

4,14
The calculation method is similar to the previous one, and it is omitted here for reasons of space. From the above results, we can get the exact symplectic propagation coefficient through numerical calculation, and at this time, the unique symplectic operator of the system is optimal, and the error function can obtain the minimum value. The relationship between each symplectic operator and d1 is shown in Figure 3. In the fourth-order, each symplectic operator constitutes 14 equations, which are

Comparison with Different Symplectic Operators
For the ordinary numerical analysis method, as the calculation progresses, the numerical solution obtained and the cumulative error of the Hamilton system will become larger and larger. SFDTD represents the Symplectic Finite-Difference Time-Domain. The symplectic algorithm based on the decomposition symplectic operator can not only obtain the exact value of the symplectic propagation coefficient, but also effectively reduce the dispersion error. Table 1 gives the values of different symplectic operators. By comparing the second-order with different order symplectic operators, Figure 4 clearly shows that the dispersion error gradually increases with the increase of the order while keeping the series constant. In the same order, by comparing the symplectic operators of different series in Figure 5, it is also found that the higher the number of symplectic operators, the better the dispersion characteristics.    Comparing the dispersion error of different methods in Figure 6, we can clearly see that the imaginary operator that we have solved by pure mathematical derivation can achieve good dispersion characteristics.   Comparing the dispersion error of different methods in Figure 6, we can clearly see that the imaginary operator that we have solved by pure mathematical derivation can achieve good dispersion characteristics. Comparing the dispersion error of different methods in Figure 6, we can clearly see that the imaginary operator that we have solved by pure mathematical derivation can achieve good dispersion characteristics. . By Comparing the symplectic algorithm of the same order number obtained by the current research with the symplectic algorithm obtained by Ruth [4], we can see that our method is significantly better than the traditional method.

Application in the Intrinsic Problem of Nanodevices
In quantum mechanics, an infinitely deep square well is a potential well with a potential in the well of 0 and an infinite potential outside the well. The expression is [17] ( ) = { 0, 0 ≤ ≤ ∞, , a represents the length of the box, and x represents the position of the particle.
For the cube region, the eigen energy of the quantum well is = ℎ 2 2 2 2 2 , = 1,2, … , the resonance frequency is = ℎ = 49.3928 2 , = 1,2, …, the box size is ∈ [0,34], the space step is = 1, and the stability is = ℎ ∆ 2 = 0.1, = 41,000. The energy distribution of the spectrum of the particle wave function in the box was shown in see Figure 7. We calculate the resonance frequency by using different symplectic operators respectively. For analysis and comparison, it can be clearly seen from Figure 8 that the error between the symplectic algorithm and the analytical solution obtained by theoretical calculation is smaller. Figure 6. By Comparing the symplectic algorithm of the same order number obtained by the current research with the symplectic algorithm obtained by Ruth [4], we can see that our method is significantly better than the traditional method.

Application in the Intrinsic Problem of Nanodevices
In quantum mechanics, an infinitely deep square well is a potential well with a potential in the well of 0 and an infinite potential outside the well. The expression is [17].
, a represents the length of the box, and x represents the position of the particle.
For the cube region, the eigen energy of the quantum well is E n = h 2 π 2 2ma 2 n 2 , n = 1, 2, . . ., the resonance frequency is ω n = E n h = 49.3928n 2 , n = 1, 2, . . ., the box size is x ∈ [0, 34], the space step is τ = 1, and the stability is S δ = τh ∆ 2 x m = 0.1, N max = 41, 000. The energy distribution of the spectrum of the particle wave function in the box was shown in see Figure 7. We calculate the resonance frequency by using different symplectic operators respectively. For analysis and comparison, it can be clearly seen from Figure 8 that the error between the symplectic algorithm and the analytical solution obtained by theoretical calculation is smaller. ℎ = 1, and the stability is = ℎ ∆ 2 = 0.1, = 41,000. The energy distribution of the spectrum of the particle wave function in the box was shown in see Figure 7. We calculate the resonance frequency by using different symplectic operators respectively. For analysis and comparison, it can be clearly seen from Figure 8 that the error between the symplectic algorithm and the analytical solution obtained by theoretical calculation is smaller.

Conclusions
In view of the insufficiency of the dispersion error caused by simplification and many artificial assumptions such as symmetry and time reversibility in the process of solving traditional symplectic operators, we calculate different orders by solving the Hamilton equation with purely mathematical derivation. The new symplectic operators correspond to different series of numbers without making any assumptions and artificial simplifications, and reduce the dispersion error from the source. On the basis of comparing different symplectic operators such as (2,4), (4,2), (4,3), and (4,4), we have selected the fourth-order as the main research object. The fourth-order new symplectic operator is applied to the intrinsic problem of nanodevices to solve the one-dimensional Schrödinger equation. Comparing the symplectic operator method with the traditional FDTD method and Ruth's method, the comparison of the energy distribution of the wave function spectrum calculated by these three methods and the real error value of the wave function at the center frequency fully proves that our method is feasible because of its higher precision and lower dispersion error. This method can also be widely used in other fields, and its dispersion characteristics require further testing.
Author Contributions: These authors contributed equally to this work.
Funding: This work was supported by the National Natural Science Foundation of China 61271351 and 41827807.

Conflicts of Interest:
The authors declare no conflict of interest.

Conclusions
In view of the insufficiency of the dispersion error caused by simplification and many artificial assumptions such as symmetry and time reversibility in the process of solving traditional symplectic operators, we calculate different orders by solving the Hamilton equation with purely mathematical derivation. The new symplectic operators correspond to different series of numbers without making any assumptions and artificial simplifications, and reduce the dispersion error from the source. On the basis of comparing different symplectic operators such as (2,4), (4,2), (4,3), and (4,4), we have selected the fourth-order as the main research object. The fourth-order new symplectic operator is applied to the intrinsic problem of nanodevices to solve the one-dimensional Schrödinger equation. Comparing the symplectic operator method with the traditional FDTD method and Ruth's method, the comparison of the energy distribution of the wave function spectrum calculated by these three methods and the real error value of the wave function at the center frequency fully proves that our method is feasible because of its higher precision and lower dispersion error. This method can also be widely used in other fields, and its dispersion characteristics require further testing.