Computer simulation of quantum dynamics in a classical spin environment

Quantum-Classical Hybrid Systems and Ehrenfest's Theorem

The conceptual analysis of quantum mechanics brings to light that a theory inherently consistent with observations should be able to describe both quantum and classical systems, i.e., quantum-classical hybrids. For example, the orthodox interpretation of measurements requires the transient creation of quantum-classical hybrids. Despite its limitations in defining the classical limit, Ehrenfest's theorem makes the simplest contact between quantum and classical mechanics. Here, we generalized the Ehrenfest theorem to bipartite quantum systems. To study quantum-classical hybrids, we employed a formalism based on operator-valued Wigner functions and quantum-classical brackets. We used this approach to derive the form of the Ehrenfest theorem for quantum-classical hybrids. We found that the time variation of the average energy of each component of the bipartite system is equal to the average of the symmetrized quantum dissipated power in both the quantum and the quantum-classical case. We expect that these theoretical results will be useful both to analyze quantum-classical hybrids and to develop self-consistent numerical algorithms for Ehrenfest-type simulations.

DECIDE: A Deterministic Mixed Quantum-Classical Dynamics Approach

Mixed quantum-classical dynamics provides an efficient way of simulating the dynamics of quantum subsystems coupled to many-body environments. Many processes, including proton-transfer reactions, electron-transfer reactions, and vibrational energy transport, for example, take place in such open systems. The most accurate algorithms for performing mixed quantum-classical simulations require very large ensembles of trajectories to obtain converged expectation values, which is computationally prohibitive for quantum subsystems containing even a few degrees of freedom. The recently developed “Deterministic evolution of coordinates with initial decoupled equations” (DECIDE) method has demonstrated high accuracy and low computational cost for a host of model systems; however, these applications relied on representing the equations of motion in subsystem and adiabatic energy bases. While these representations are convenient for certain systems, the position representation is convenient for many other systems, including systems undergoing proton- and electron-transfer reactions. Thus, in this review, we provide a step-by-step derivation of the DECIDE approach and demonstrate how to cast the DECIDE equations in a quantum harmonic oscillator position basis for a simple one-dimensional (1D) hydrogen bond model. After integrating the DECIDE equations of motion on this basis, we show that the total energy of the system is conserved for this model and calculate various quantities of interest. Limitations of casting the equations in an incomplete basis are also discussed.

Symmetry-Induced Emergence of a Pseudo-Qutrit in the Dipolar Coupling of Two Qubits

We investigate a system of two identical and distinguishable spins 1/2, with a direct magnetic dipole-dipole interaction, in an external magnetic field. Constraining the hyperfine tensor to exhibit axial symmetry generates the notable symmetry properties of the corresponding Hamiltonian model. In fact, we show that the reduction of the anisotropy induces the invariance of the Hamiltonian in the 3×3 subspace of the Hilbert space of the two spins in which S^2 invariably assumes its highest eigenvalue of 2. By means of appropriate mapping, it is then possible to choose initial density matrices of the two-spin system that evolve in such a way as to exactly simulate the time evolution of a pseudo-qutrit, in the sense that the the actual two-spin system nests the subdynamics of a qutrit regardless of the strength of the magnetic field. The occurrence of this dynamic similitude is investigated using two types of representation for the initial density matrix of the two spins. We show that the qutrit state emerges when the initial polarizations and probability vectors of the two spins are equal to each other. Further restrictions on the components of the probability vectors are reported and discussed.

Quasi-Lie Brackets and the Breaking of Time-Translation Symmetry for Quantum Systems Embedded in Classical Baths

Many open quantum systems encountered in both natural and synthetic situations are embedded in classical-like baths. Often, the bath degrees of freedom may be represented in terms of canonically conjugate coordinates, but in some cases they may require a non-canonical or non-Hamiltonian representation. Herein, we review an approach to the dynamics and statistical mechanics of quantum subsystems embedded in either non-canonical or non-Hamiltonian classical-like baths which is based on operator-valued quasi-probability functions. These functions typically evolve through the action of quasi-Lie brackets and their associated Quantum-Classical Liouville Equations, or through quasi-Lie brackets augmented by dissipative terms. Quasi-Lie brackets possess the unique feature that, while conserving the energy (which the Noether theorem links to time-translation symmetry), they violate the time-translation symmetry of their algebra. This fact can be heuristically understood in terms of the dynamics of the open quantum subsystem. We then describe an example in which a quantum subsystem is embedded in a bath of classical spins, which are described by non-canonical coordinates. In this case, it has been shown that an off-diagonal open-bath geometric phase enters into the propagation of the quantum-classical dynamics. Next, we discuss how non-Hamiltonian dynamics may be employed to generate the constant-temperature evolution of phase space degrees of freedom coupled to the quantum subsystem. Constant-temperature dynamics may be generated by either a classical Langevin stochastic process or a Nosé–Hoover deterministic thermostat. These two approaches are not equivalent but have different advantages and drawbacks. In all cases, the calculation of the operator-valued quasi-probability function allows one to compute time-dependent statistical averages of observables. This may be accomplished in practice using a hybrid Molecular Dynamics/Monte Carlo algorithms, which we outline herein.