In molecular dynamics applications there is a growing interest in mixed quantum-classical models. The quantum-classical Liouville equation (QCL) describes most atoms of the molecular system under consideration by means of a classical phase space density but an important, small portion of the system by means of quantum mechanics. The QCL is derived from the full quantum dynamical (QD) description by applying the Wigner transform to the \\classical part" of the system only. We discuss the conditions under which the QCL model approximates the full QD evolution of the system. First, analysis of the asymptotic properties of the Wigner transform shows that solving the QCL yields a rst order approximation of full quantum dynamics. Second, we discuss the adiabatic limit of the QCL. This discussion shows that the QCL solutions may be interpretated as classical phase space densities, at least near the adiabatic limit. Third, it is demonstrated that the QCL yields good approximations of non-adia...
In molecular dynamics applications there is a growing interest in mixed quantum-classical models. The quantum-classical Liouville equation (QCL) describes most atoms of the molecular system under consideration by means of a classical phase space density but an important, small portion of the system by means of quantum mechanics. The QCL is derived from the full quantum dynamical (QD) description by applying the Wigner transform to the \classical part" of the system only. We discuss the conditions under which the QCL model approximates the full QD evolution of the system. First, analysis of the asymptotic properties of the Wigner transform shows that solving the QCL yields a rst order approximation of full quantum dynamics. Second, we discuss the adiabatic limit of the QCL. This discussion shows that the QCL solutions may be interpretated as classical phase space densities, at least near the adiabatic limit. Third, it is demonstrated that the QCL yields good approximations of non-adia...
A numerical technique for the solution of the Quantum Liouville equation is discussed which is expected to be useful for mixed states describing highly excited many body systems.;The equation is represented in phase-space, where the density operator becomes the Wigner function. The Liouville operator which determines its time evolution, separates into a classical and quantum term. When the Wigner function varies slowly in momentum, the quasi-classical approximation is valid in which one keeps only the first order quantum term in the Liouville equation.;The Wigner function is represented by a sample set of gaussian smeared points. In an iterative time development scheme, the points are propagated first along a classical phase-space trajectory, which is then followed by a quantum development. The latter is represented by a stochastic process incorporated into the numerical procedure as a generalized Markovian process.;This method is numerically tested for an anharmonic quartic potential. The contours of the Wigner function at different time intervals are compared with the purely classical and exact results. It is seen that it provides a definite improvement over the purely classical approximation. Various averages are also compared and sources of error discussed.;In addition it was discovered that intermediate mixed states show a 'quantum focusing' effect, wherein the function peaks to values greater than the original maximum. This quantum-mechanical effect did not appear in the pure case, and disappears because of Liouville's theorem for highly mixed states which are essentially classical. There are other numerical schemes, which treat the points as moving in an effective potential, wherein the points retain their identity. These again by Liouville's theorem cannot show this focusing effect.
Mixed quantum-classical equations of motion are derived for a quantum subsystem of light (mass m) particles coupled to a classical bath of massive (mass M) particles. The equation of motion follows from a partial Wigner transform over the bath degrees of freedom of the Liouville equation for the full quantum system, followed by an expansion in the small parameter μ=(m/M)1/2 in analogy with the theory of Brownian motion. The resulting mixed quantum-classical Liouville equation accounts for the coupled evolution of the subsystem and bath. The quantum subsystem is represented in an adiabatic (or other) basis and the series solution of the Liouville equation leads to a representation of the dynamics in an ensemble of surface-hopping trajectories. A generalized Pauli master equation for the evolution of the diagonal elements of the density matrix is derived by projection operator methods and its structure is analyzed in terms of surface-hopping trajectories.
The mixed quantum-classical Liouville equation (QCLE) provides an approximate perturbative framework for describing the dynamics of systems with coupled quantum and classical degrees of freedom of disparate thermal wavelengths. The evolution governed by the Liouville operator preserves many properties of full quantum dynamics, including the conservation of total population, energy, and purity, and has shown quantitative agreement with exact quantum results for the expectation values of many observables where direct comparisons are feasible. However, since the QCLE density matrix operator is obtained from the partial Wigner transform of the full quantum density matrix, its matrix elements can have negative values, implying that the diagonal matrix elements behave as pseudo-densities rather than densities of classical phase space. Here, we compare phase-space distributions generated by exact quantum dynamics with those produced by QCLE evolution from pure quantum initial states. We show that resonance effects in the off-diagonal matrix elements differ qualitatively, particularly for low-energy states. Furthermore, numerical and analytical results for low-dimensional models reveal that the QCLE can violate the positivity of marginal phase-space densities, a property that should hold at all times for any physical system. A perturbative analysis of a model system confirms that such violations arise generically. We also show that the violations of positivity of the marginal densities vanish as the initial energy of the system increases relative to the energy gap between subsystem states. These findings suggest that a negativity index, quantifying deviations from positivity, may provide a useful metric for assessing the validity of mixed quantum-classical descriptions.
The mixed quantum-classical Liouville equation (QCLE) provides an approximate perturbative framework for describing the dynamics of systems with coupled quantum and classical degrees of freedom of disparate thermal wavelengths. The evolution governed by the Liouville operator preserves many properties of full quantum dynamics, including the conservation of total population, energy, and purity, and has shown quantitative agreement with exact quantum results for the expectation values of many observables where direct comparisons are feasible. However, since the QCLE density matrix operator is obtained from the partial Wigner transform of the full quantum density matrix, its matrix elements can have negative values, implying that the diagonal matrix elements behave as pseudo-densities rather than densities of classical phase space. Here, we compare phase-space distributions generated by exact quantum dynamics with those produced by QCLE evolution from pure quantum initial states. We show that resonance effects in the off-diagonal matrix elements differ qualitatively, particularly for low-energy states. Furthermore, numerical and analytical results for low-dimensional models reveal that the QCLE can violate the positivity of marginal phase-space densities, a property that should hold at all times for any physical system. A perturbative analysis of a model system confirms that such violations arise generically. We also show that the violations of positivity of the marginal densities vanish as the initial energy of the system increases relative to the energy gap between subsystem states. These findings suggest that a negativity index, quantifying deviations from positivity, may provide a useful metric for assessing the validity of mixed quantum\textendash{}classical descriptions.
This paper presents a mathematical derivation of a model for quantum-classical molecular dynamics (QCMD) as a partial classical limit of the full Schrödinger equation. This limit is achieved in two steps: separation of the full wave function and short wave asymptotics for its ‘‘classical’’ part. Both steps can be rigorously justified under the same smallness assumptions. This throws some light on the time-dependent self-consistent-field method and on mixed quantum-semiclassical models, which also depend on the separation step. On the other hand, the theory leads to a characterization of the critical situations in which the QCMD model is in danger of largely deviating from the solution of full Schrödinger equation. These critical situations are exemplified in an illustrative numerical simulation: the collinear collision of a classical particle with a harmonic quantum oscillator.
Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.