© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1310.7455 · The Journal of Chemical Physics

The rigorous stochastic matrix multiplication scheme for the calculations of reduced equilibrium density matrices of open multilevel quantum systems

Xin Chen

Published Apr 15, 2014• 1 citation

Abstract

Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, i.e., how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernon influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems.

DOI: 10.1063/1.4870496

Related papers

The Rigorous Stochastic Matrix Product Method for the Calculations of Reduced Equilibrium Density Matrices of Multilevel Open Quantum Systems

Xin Chen

•Oct 28, 2013•0 citations

A rigorous approach to the Feynman-Vernon influence functional and its applications. I

Sergio A Albeverio, Laura Cattaneo, Sonia Mazzucchi +1 more

A rigorous representation of the Feynman-Vernon influence functional used to describe open quantum systems is given, based on the theory of infinite dimensional oscillatory integrals. An application to the case of the density matrices describing the Caldeira-Leggett model of two quantum systems with a quadratic interaction is treated.

Journal of Mathematical Physics•Oct 1, 2007•9 citations

Tensor propagator for iterative quantum time evolution of reduced density matrices. I. Theory

Nancy Makri, Dmitrii E. Makarov

For common condensed phase problems described by a low-dimensional system coupled to a harmonic bath, Feynman’s path integral formulation of time-dependent quantum mechanics leads to expressions for the reduced density matrix of the system where the effects of the harmonic environment enter through an influence functional that is nonlocal in time. In a recent Letter [Chem. Phys. Lett. 221, 482 (1994)], we demonstrated that the range of the nonlocal interactions is finite even at zero temperature, such that the nonlocal kernel extends over only a few time steps if the path integral is expressed in terms of accurate quasiadiabatic propagators. This feature arises from disruption of phase coherence in macroscopic environments and leads to Markovian dynamics for an augmented reduced density tensor, permitting iterative time evolution schemes. In the present paper we analyze the structure and properties of the relevant tensor propagator. Specifically, we show that the tensor multiplication scheme rigorously conserves the trace of the reduced density matrix, and that in cases of short-range nonlocality it leads to Redfield-type equations which are correct to all orders in perturbation theory and which take into account memory effects. We also argue that a simple eigenvector analysis reveals (without actual iteration) the nature of the dynamics and of the equilibrium state, and directly yields quantum reaction or relaxation rates.

The Journal of Chemical Physics•Mar 15, 1995•636 citations

An extension of stochastic hierarchy equations of motion for the equilibrium correlation functions

Yaling Ke, Yi Zhao

A traditional stochastic hierarchy equations of motion method is extended into the correlated real-time and imaginary-time propagations, in this paper, for its applications in calculating the equilibrium correlation functions. The central idea is based on a combined employment of stochastic unravelling and hierarchical techniques for the temperature-dependent and temperature-free parts of the influence functional, respectively, in the path integral formalism of the open quantum systems coupled to a harmonic bath. The feasibility and validity of the proposed method are justified in the emission spectra of homodimer compared to those obtained through the deterministic hierarchy equations of motion. Besides, it is interesting to find that the complex noises generated from a small portion of real-time and imaginary-time cross terms can be safely dropped to produce the stable and accurate position and flux correlation functions in a broad parameter regime.

The Journal of Chemical Physics•Jun 1, 2017•18 citations

Influence functionals: General methodology for subsystem calculations

Kurt Allinger, Mark A. Ratner

Feynman path integrals offer a convenient and elegant tool for dealing with the statistical physics of quantum-mechanical systems. By using the Trotter product formula, one can evaluate directly the propagation of many-body quantum systems in imaginary time. The idea of an influence functional arises when one wishes to consider division of the physical system into a primary system and a bath. We present a general approach to the problem of constructing influence functionals, one that is capable of dealing both with Boltzmann baths (classical oscillators) and quantum baths of bosons or fermions. The fermion bath of special interest is generally a submanifold of the electronic state such as the bridge structure connecting chromophores or electron localization subunits which are common in problems of mixed valency and superexchange. The fermion bath is treated by a general rewriting of the influence phase S as a sum of an eigenvalue part ${S}_{\mathrm{EV}}$ and a remainder, ${S}^{0}$, that describes transitions. Reduction of the original Hamiltonian leads to a reduced Hamiltonian with effective off-diagonal matrix elements and influence functionals which can include memory effects. We present a general formalism for construction of influence functionals, and discuss particular applications to systems of interest, especially in electronic-structure problems. The following paper presents a variational technique for finding effective Hamiltonians.

Physical Review A•Jan 1, 1989•8 citations

Quantum Stochastic Trajectories: The Fokker–Planck–Bohm Equation Driven by the Reduced Density Matrix

Francesco Avanzini, Giorgio J. Moro

The quantum molecular trajectory is the deterministic trajectory, arising from the Bohm theory, that describes the instantaneous positions of the nuclei of molecules by assuring the agreement with the predictions of quantum mechanics. Therefore, it provides the suitable framework for representing the geometry and the motions of molecules without neglecting their quantum nature. However, the quantum molecular trajectory is extremely demanding from the computational point of view, and this strongly limits its applications. To overcome such a drawback, we derive a stochastic representation of the quantum molecular trajectory, through projection operator techniques, for the degrees of freedom of an open quantum system. The resulting Fokker-Planck operator is parametrically dependent upon the reduced density matrix of the open system. Because of the pilot role played by the reduced density matrix, this stochastic approach is able to represent accurately the main features of the open system motions both at equilibrium and out of equilibrium with the environment. To verify this procedure, the predictions of the stochastic and deterministic representation are compared for a model system of six interacting harmonic oscillators, where one oscillator is taken as the open quantum system of interest. The undeniable advantage of the stochastic approach is that of providing a simplified and self-contained representation of the dynamics of the open system coordinates. Furthermore, it can be employed to study the out of equilibrium dynamics and the relaxation of quantum molecular motions during photoinduced processes, like photoinduced conformational changes and proton transfers.

The Journal of Physical Chemistry A•Feb 21, 2018•1 citation

Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, i.e., how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernon influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems.

For common condensed phase problems described by a low-dimensional system coupled to a harmonic bath, Feynman’s path integral formulation of time-dependent quantum mechanics leads to expressions for the reduced density matrix of the system where the effects of the harmonic environment enter through an influence functional that is nonlocal in time. In a recent Letter [Chem. Phys. Lett. 221, 482 (1994)], we demonstrated that the range of the nonlocal interactions is finite even at zero temperature, such that the nonlocal kernel extends over only a few time steps if the path integral is expressed in terms of accurate quasiadiabatic propagators. This feature arises from disruption of phase coherence in macroscopic environments and leads to Markovian dynamics for an augmented reduced density tensor, permitting iterative time evolution schemes. In the present paper we analyze the structure and properties of the relevant tensor propagator. Specifically, we show that the tensor multiplication scheme rigorously conserves the trace of the reduced density matrix, and that in cases of short-range nonlocality it leads to Redfield-type equations which are correct to all orders in perturbation theory and which take into account memory effects. We also argue that a simple eigenvector analysis reveals (without actual iteration) the nature of the dynamics and of the equilibrium state, and directly yields quantum reaction or relaxation rates.

A traditional stochastic hierarchy equations of motion method is extended into the correlated real-time and imaginary-time propagations, in this paper, for its applications in calculating the equilibrium correlation functions. The central idea is based on a combined employment of stochastic unravelling and hierarchical techniques for the temperature-dependent and temperature-free parts of the influence functional, respectively, in the path integral formalism of the open quantum systems coupled to a harmonic bath. The feasibility and validity of the proposed method are justified in the emission spectra of homodimer compared to those obtained through the deterministic hierarchy equations of motion. Besides, it is interesting to find that the complex noises generated from a small portion of real-time and imaginary-time cross terms can be safely dropped to produce the stable and accurate position and flux correlation functions in a broad parameter regime.

Feynman path integrals offer a convenient and elegant tool for dealing with the statistical physics of quantum-mechanical systems. By using the Trotter product formula, one can evaluate directly the propagation of many-body quantum systems in imaginary time. The idea of an influence functional arises when one wishes to consider division of the physical system into a primary system and a bath. We present a general approach to the problem of constructing influence functionals, one that is capable of dealing both with Boltzmann baths (classical oscillators) and quantum baths of bosons or fermions. The fermion bath of special interest is generally a submanifold of the electronic state such as the bridge structure connecting chromophores or electron localization subunits which are common in problems of mixed valency and superexchange. The fermion bath is treated by a general rewriting of the influence phase S as a sum of an eigenvalue part

{S}_{\mathrm{EV}}

and a remainder,

{S}^{0}

, that describes transitions. Reduction of the original Hamiltonian leads to a reduced Hamiltonian with effective off-diagonal matrix elements and influence functionals which can include memory effects. We present a general formalism for construction of influence functionals, and discuss particular applications to systems of interest, especially in electronic-structure problems. The following paper presents a variational technique for finding effective Hamiltonians.

The quantum molecular trajectory is the deterministic trajectory, arising from the Bohm theory, that describes the instantaneous positions of the nuclei of molecules by assuring the agreement with the predictions of quantum mechanics. Therefore, it provides the suitable framework for representing the geometry and the motions of molecules without neglecting their quantum nature. However, the quantum molecular trajectory is extremely demanding from the computational point of view, and this strongly limits its applications. To overcome such a drawback, we derive a stochastic representation of the quantum molecular trajectory, through projection operator techniques, for the degrees of freedom of an open quantum system. The resulting Fokker-Planck operator is parametrically dependent upon the reduced density matrix of the open system. Because of the pilot role played by the reduced density matrix, this stochastic approach is able to represent accurately the main features of the open system motions both at equilibrium and out of equilibrium with the environment. To verify this procedure, the predictions of the stochastic and deterministic representation are compared for a model system of six interacting harmonic oscillators, where one oscillator is taken as the open quantum system of interest. The undeniable advantage of the stochastic approach is that of providing a simplified and self-contained representation of the dynamics of the open system coordinates. Furthermore, it can be employed to study the out of equilibrium dynamics and the relaxation of quantum molecular motions during photoinduced processes, like photoinduced conformational changes and proton transfers.