Entropy

Correlation Functions in Open Quantum-Classical Systems

Chang‐Yu Hsieh, Raymond Kapral

Published Dec 27, 2013• 10 citations

Abstract

Quantum time correlation functions are often the principal objects of interest in experimental investigations of the dynamics of quantum systems. For instance, transport properties, such as diffusion and reaction rate coefficients, can be obtained by integrating these functions. The evaluation of such correlation functions entails sampling from quantum equilibrium density operators and quantum time evolution of operators. For condensed phase and complex systems, where quantum dynamics is difficult to carry out, approximations must often be made to compute these functions. We present a general scheme for the computation of correlation functions, which preserves the full quantum equilibrium structure of the system and approximates the time evolution with quantum-classical Liouville dynamics. Several aspects of the scheme are discussed, including a practical and general approach to sample the quantum equilibrium density, the properties of the quantum-classical Liouville equation in the context of correlation function computations, simulation schemes for the approximate dynamics and their interpretation and connections to other approximate quantum dynamical methods.

DOI: 10.3390/e16010200

Related papers

Dynamics of quantum–classical systems in nonequilibrium environments

Jeremy Schofield, Raymond Kapral

The dynamics of a quantum system coupled to a classical environment and subject to constraints that drive it out of equilibrium are described. The evolution of the system is governed by the quantum-classical Liouville equation. Rather than evaluating the evolution of the mixed quantum-classical density operator, we derive exact equations of motion for the nonequilibrium average values of a set of operators or variables, along with correlation function expressions for the dissipative coefficients that enter these equations. These equations are obtained by requiring that the exact nonequilibrium averages are equal to local nonequilibrium averages that depend on auxiliary fields whose values satisfy evolution equations obtained using projection operator methods. The results are illustrated by deriving reaction-diffusion equations coupled to fluid hydrodynamic equations for a solution of quantum particles that can exist in two metastable states. Nonequilibrium steady states are discussed along with the reaction rate and diffusion correlation functions that characterize such states.

arXiv:2411.18713The Journal of Chemical PhysicsFeb 24, 20250 citations

A Study of Quantum-classical Dynamics in the Mapping Basis

Ali M. Nassimi

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.

Aug 31, 20110 citations

Quantum time-dependent statistical mechanics

Mark E. Tuckerman

Abstract Chapter 14 covers the statistical mechanics of driven quantum systems and the non-equilibrium ensembles associated with them. The topic is presented both from the state-vector and density matrix perspectives. The chapter begins with a discussion of driving forces resulting from light-matter interactions and their treatment via time-dependent perturbation theory. The Fermi Golden Rule is derived, and it is shown how to incorporate the Golden Rule into the calculation of frequency spectra within the canonical ensemble, which then connects to canonical quantum time correlation functions. The formalism is derived again from the quantum Liouville equation. The discussion then turns to the challenge of computing quantum time correlation functions and approximate techniques based on the Feynman path integral of Chapter 12. These schemes include centroid and ring-polymer molecular dynamics as well as open-path integral methods

Aug 1, 20231 citation

Recent advances in the calculation of dynamical correlation functions

J. Florencio, O. F. de Alcantara Bonfim

We review various theoretical methods that have been used in recent years to calculate dynamical correlation functions of many-body systems. Time-dependent correlation functions and their associated frequency spectral densities are the quantities of interest, for they play a central role in both the theoretical and experimental understanding of dynamic properties. The calculation of the relaxation function is rather difficult in most cases of interest, except for a few examples where exact analytic expressions are allowed. For most of systems of interest approximation schemes must be used. The method of recurrence relation has, at its foundation, the solution of Heisenberg equation of motion of an operator in a many-body interacting system. Insights have been gained from theorems that were discovered with that method. For instance, the absence of pure exponential behavior for the relaxation functions of any Hamiltonian system. The method of recurrence relations was used in quantum systems such as dense electron gas, transverse Ising model, Heisenberg model, XY model, Heisenberg model with Dzyaloshinskii-Moriya interactions, as well as classical harmonic oscillator chains. Effects of disorder were considered in some of those systems. In the cases where analytical solutions were not feasible, approximation schemes were used, but are highly model-dependent. Another important approach is the numerically exact diagonalization method. It is used in finite-sized systems, which sometimes provides very reliable information of the dynamics at the infinite-size limit. In this work, we discuss the most relevant applications of the method of recurrence relations and numerical calculations based on exact diagonalizations.

arXiv:2005.14222May 28, 20201 citation

Approximating Time-Dependent Quantum Statistical Properties

Sara Bonella, Giovanni Ciccotti

Computing quantum dynamics in condensed matter systems is an open challenge due to the exponential scaling of exact algorithms with the number of degrees of freedom. Current methods try to reduce the cost of the calculation using classical dynamics as the key ingredient of approximations of the quantum time evolution. Two main approaches exist, quantum classical and semi-classical, but they suffer from various difficulties, in particular when trying to go beyond the classical approximation. It may then be useful to reconsider the problem focusing on statistical time-dependent averages rather than directly on the dynamics. In this paper, we discuss a recently developed scheme for calculating symmetrized correlation functions. In this scheme, the full (complex time) evolution is broken into segments alternating thermal and real-time propagation, and the latter is reduced to classical dynamics via a linearization approximation. Increasing the number of segments systematically improves the result with respect to full classical dynamics, but at a cost which is still prohibitive. If only one segment is considered, a cumulant expansion can be used to obtain a computationally efficient algorithm, which has proven accurate for condensed phase systems in moderately quantum regimes. This scheme is summarized in the second part of the paper. We conclude by outlining how the cumulant expansion formally provides a way to improve convergence also for more than one segment. Future work will focus on testing the numerical performance of this extension and, more importantly, on investigating the limit for the number of segments that goes to infinity of the approximate expression for the symmetrized correlation function to assess formally its convergence to the exact result.

EntropyDec 27, 20133 citations

A perturbation theory for multi-time correlation functions in open quantum systems

Piotr Szańkowski

Dynamical maps are the principal subject of the open system theory. Formally, the dynamical map of a given open quantum system is a density matrix transformation that takes any initial state and sends it to the state at a later time. Physically, it encapsulates the system’s evolution due to coupling with its environment. Hence, the theory provides a flexible and accurate framework for computing expectation values of open system observables. However, expectation values-or more generally, single-time correlation functions-capture only the simplest aspects of a quantum system’s dynamics. A complete characterization of the dynamics requires access to multi-time correlation functions as well: phenomena like detailed balance, linear and non-linear response, non-equilibrium transport in general, or even sequential measurements of system observables are all described in terms of multi-time correlations. For closed systems, such correlations are well-defined, even though knowledge of the system’s state alone is insufficient to determine them fully. In contrast, the standard dynamical map formalism for open systems does not account for multi-time correlations, as it is fundamentally limited to describing state evolution. Here, we extend the scope of open quantum system theory by developing a systematic perturbation theory for computing multi-time correlation functions.

arXiv:2502.19137SciPost PhysicsSep 3, 20251 citation