Linearized density matrix propagation

The chapter is organized as follows The quantum-classical Liouville dynamics scheme is first outlined and a rigorous surface hopping trajectory algorithm for its implementation is presented. The iterative linearized density matrix propagation approach is then described and an approach for its implementation is presented. In the Model Simulations section the comparable performance of the two methods is documented for the generalized spin-boson model and numerical convergence issues are mentioned. In the Conclusions we review the perspectives of this study. 

In this section we present results using the two approaches described in the previous sections the Trotter factorized QCL (TQCL), and iterative linearized density matrix (ILDM) propagation schemes, to study the spin-boson model consisting of a two level system that is bi-linearly coupled to a bath with Mh harmonic modes. This popular model of a quantum system embedded in an environment is described by the following general hamiltonian ... 

Ultrasonic methods have been considered for non-destructive evaluation of compositional changes and mixture uniformity in filled polymers. In principle, it is feasible to determine the size, shape and distribution of filler particles, to detect agglomeration and to assess the extent of filler-matrix interaction, through appropriate application of ultrasonic procedures [48]. The method involves determination of the elastic behaviour of solids by measurement of ultrasonic wave velocity. Propagation of a plane wave in a linear elastic material can be related to its elastic modulus (E) and density (p) according to ... 

This approximation is better known as the time-dependent Hartree—Fock approximation (TDHF) (McLachlan and Ball, 1964) (see Section 11.1) or random phase approximation (RPA) (Rowe, 1968) and can also be derived as the linear response of an SCF wavefunction, as described in Section 11.2. Furthermore, the structure of the equations is the same as in time-dependent density functional theory (TD-DFT), although they differ in the expressions for the elements of the Hessian matrix E22. The polarization propagator in the RPA is then given as... 

---