Split-operator short-time propagator method

We apply the split-operator method Eq. (3.7) to separately handle the potential and kinetic energy terms of the Hamiltonian. We further apply the split-operator scheme to separate the kinetic energy term into the two exponentially noncommutative parts Tr + Tr and Tg to obtain a numerical short-time propagation method. 

Once the grid (or two grids) are prepared, there are two similar types of approaches to propagate the initial wavefiinction forward with time. One approach is split-operator methods, [59] where the short-time propagator is divided into a kinetic and potential parts so that... 

The split-operator (SP) method is extremely popular and has been widely used in many practical applications. It approximates the short-time propagator by the equation... 

Within the split-operator scheme, for example, numerical short-time propagation would involve diagonalizing the interaction matrix (either together with the potential energy matrix or as a separate term) at each time step. We will see several other time-propagation methods in Sec. 5.2.1. 

The Hamiltonian and the coordinates are discretized by means of the generalized pseudospectral (GPS) method in prolate spheroidal coordinates [44-47], allowing optimal and nonuniform spatial grid distribution and accurate solution of the wave functions. The time-dependent Kohn-Sham Equation 3.5 can be solved accurately and efficiently by means of the split-operator method in the energy representation with spectral expansion of the propagator matrices [44-46,48]. We employ the following split operator, second-order short-time propagation formula [40] ... 

We may apply the split-operator method (Sec. 3.2.1) to the three matrices lijV) V/), and Vo for a short-time numerical propagation scheme. Ionization is now described by population of the neutral state Xn R,t) transferring to the ionized state partial-wave components Xc,kjix R,t) over time through the interaction represented by the matrix Vo(i ,f)-... 

Usually, the propagator (7(r, to) is approximated by various schemes [55,60,137], and there are plenty of wonderful articles that have explained each in detail, such as the split operator method and higher order split operator methods [11, 36, 130], Chebyshev polynomial expansion [131], Faber polynomial expansion [51, 146], short iterative Lanczos propagation method [95], Crank-Nicholson second-order differencing [10,56,57], symplectic method [14,45], recently proposed real Chebyshev method [24,44,125], and Multi-configuration Time-Dependent Hartree (MCTDH) Method [ 12,73,81-83]. For details, one may refer to the corresponding references. 

---