Split-operator

A more powerfiil method for evaluating the time derivative of the wavefiinction is the split-operator method [39]. Flere we start by fomially solving ihd ild. = /7 / with the solution V fD = e Note that //is... 

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... 

An alternative to split operator methods is to use iterative approaches. In these metiiods, one notes that the wavefiinction is fomially "tt(0) = exp(-i/7oi " ), and the action of the exponential operator is obtained by repetitive application of //on a function (i.e. on the computer, by repetitive applications of the sparse matrix... 

This procedure is then repeated after each time step. Comparison with Eq. (2) shows that the result is the velocity Verlet integrator and we have thus derived it from a split-operator technique which is not the way that it was originally derived. A simple interchange of the Ly and L2 operators yields an entirely equivalent integrator. 

An alternative relax-and-drive procedure can be based on a strictly unitary treatment where the advance from Iq to t is done with a norm-conserving propagation such as provided by the split-operator propagation technique.(49, 50) This however is more laborious, and although it conserves the norm of the density matrix, it is not necessarily more accurate because of possible inaccuracies in the individual (complex) density matrix elements. It can however be used to advantage when the dimension of the density matrix is small and exponentiation of matrices can be easily done.(51, 52)... 

The following two sections present a brief overview of the split-operator method, as used in several recent applications [41, 42, 61, 62], and of the basis set expansion approach. 

The initial wavefunction is then expanded as in Eq. (2), and the wavepacket is propagated using the split operator method42 ... 

PES evaluations. In either case (coupled or split-operator with frozen Gaussian propagator)... 

Figure 1.18 shows energy levels for d orbitals in crystal fields of differing symmetry. The splitting operated by the octahedral field is much higher than that of the tetrahedral field (A = lower than the effect imposed by the square... 

Our main concern in this section is with the actual propagation forward in time of the wavepacket. The standard ways of solving the time-dependent Schrodinger equation are the Chebyshev expansion method proposed and popularised by Kossloff [16,18,20,37 0] and the split-operator method of Feit and Fleck [19,163,164]. I will not discuss these methods here as they have been amply reviewed in the references just quoted. Comparative studies [17-19] show conclusively that the Chebyshev expansion method is the most accurate and stable but the split-operator method allows for explicit time dependence in the Hamiltonian operator and is often faster when ultimate accuracy is not required. All methods for solving the time propagation of the wavepacket require the repeated operation of the Hamiltonian operator on the wavepacket. It is this aspect of the propagation that I will discuss in this section. 

In order to better understand the experimental results, we performed quantum mechanical calculations using the fast Fourier transform (FFT) split-operator technique, which was previously employed by Meier and Engel [38]... 

The time-dependence of the wavepacket evolving on any potential surface can be numerically determined by using the split operator technique of Feit and Fleck [10-15]. A good introductory overview of the method is given in Ref. [12]. We will discuss a potential in two coordinates because this example is relevant to the experimental spectra. The time-dependent Schrodinger equation in two coordinates Qx and Qy is... 

Z.G. Sun, S.Y. Lee, H. Guo, D.H. Zhang, Comparison of second-order split operator and Chebyshev propagator in wave packet based state-to-state reactive scattering calculations, J. Chem. Phys. 130 (2009) 174102. 

Extension of the equilibrium model to column or field conditions requires coupling the ion-exchange equations with the transport equations for the 5 aqueous species (Eq. 1). To accomplish this coupling, we have adopted the split-operator approach (e.g., Miller and Rabideau, 1993), which provides considerable flexibility in adjusting the sorption submodel. In addition to the above conceptual model, we are pursuing more complex formulations that couple cation exchange with pore diffusion, surface diffusion, or combined pore/surface diffusion (e.g., Robinson et al., 1994 DePaoli and Perona, 1996 Ma et al., 1996). However, the currently available data are inadequate to parameterize such models, and the need for a kinetic formulation for the low-flow conditions expected for sorbing barriers has not been established. These issues will be addressed in a future publication. 

Miller, C. T., and Rabideau, A. J. (1993). Development of split-operator Petrov-Galerkin methods for simulating transport and diffusion problems, Water Resources Research, 29(7), 2227-2240. 

---