Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Matching-pursuit/split-operator-Fourier-transform computations of thermal correlation functions

Xin Chen, Yinghua Wu, Víctor S. Batista

Published Feb 1, 2005• 11 citations
Publisher

Abstract

A rigorous and practical methodology for evaluating thermal-equilibrium density matrices, finite-temperature time-dependent expectation values, and time-correlation functions is described. The method involves an extension of the matching-pursuit/split-operator-Fourier-transform method to the solution of the Bloch equation via imaginary-time propagation of the density matrix and the evaluation of Heisenberg time-evolution operators through real-time propagation in dynamically adaptive coherent-state representations.

DOI
10.1063/1.1848513

Related papers

Path integral density matrix dynamics: A method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems

Scott Habershon

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.

The Journal of Chemical Physics•Sep 12, 2013•8 citations

Matrix Product Density Operators: Simulation of Finite-Temperature and Dissipative Systems

Frank Verstraete, Juan José García‐Ripoll, J. I. Cirac

We show how to simulate numerically the evolution of 1D quantum systems under dissipation as well as in thermal equilibrium. The method applies to both finite and inhomogeneous systems, and it is based on two ideas: (a) a representation for density operators which extends that of matrix product states to mixed states; (b) an algorithm to approximate the evolution (in real or imaginary time) of matrix product states which is variational.

arXiv:cond-mat/0406426•Physical Review Letters•Nov 12, 2004•1,093 citations

Approaches in evaluating two-time correlation function

Sintayehu Tesfa

Derivation of the procedures that can be applied in evaluating two-time correlation function in terms of coherent-state propagator and corresponding Q-function is presented. On the basis that the involved functions are generally exponential in nature, obtaining the two-time second-order correlation function is essentially claimed to be reduced to carrying out relatively simple integrations. Fundamentally, the time dependence of the operators is transferred to the density operator. Moreover, manipulation in reordering the operators is performed by applying the usual trace operation. With all details, it is basically observed that the two-time correlation can be readily determined once the pertinent coherent-state propagator or Q-function is known. Since working with c-number equation is far more handy than the associated operator equation, it is expected that the results derived in this contribution can aid in easing the otherwise involving mathematical rigor.

arXiv:0910.2939•Oct 15, 2009•0 citations

Quantum Imaginary Time Propagation algorithm for preparing thermal states

Francesco Turro

Calculations at finite temperatures are fundamental in different scientific fields, from nuclear physics to condensed matter. Evolution in imaginary time is a prominent classical technique for preparing thermal states of quantum systems. We propose a new quantum algorithm that prepares thermal states based on the quantum imaginary time propagation method, using a diluted operator with ancilla qubits to overcome the non-unitarity nature of the imaginary time operator. The presented method is the first that allows us to obtain the correct thermal density matrix on a general quantum processor for a generic Hamiltonian. We prove its reliability in the actual quantum hardware computing thermal properties for two and three neutron systems.

arXiv:2306.16580•Jun 28, 2023•0 citations

A new approach to the time evolution of characteristic function of the density operator obtained by virtue of thermal entangled state representation

Feng Chen, Hong-Yi Fan

Science China Physics Mechanics and Astronomy•Sep 18, 2012•8 citations

Quantum tunneling dynamics in multidimensional systems: A matching-pursuit description

Yinghua Wu, Víctor S. Batista

Rigorous simulations of quantum tunneling dynamics in model systems with up to 20 coupled degrees of freedom are reported. The simulations implement an extension of the recently developed matching-pursuit/split-operator Fourier-transform method to complex-valued coherent-state representations. The resulting method recursively applies the time-evolution operator, as defined by the Trotter expansion to second order accuracy, in dynamically adaptive coherent-state representations generated by an approach that combines the matching-pursuit algorithm with a gradient-based optimization method.

The Journal of Chemical Physics•Jul 8, 2004•60 citations
A rigorous and practical methodology for evaluating thermal-equilibrium density matrices, finite-temperature time-dependent expectation values, and time-correlation functions is described. The method involves an extension of the matching-pursuit/split-operator-Fourier-transform method to the solution of the Bloch equation via imaginary-time propagation of the density matrix and the evaluation of Heisenberg time-evolution operators through real-time propagation in dynamically adaptive coherent-state representations.
We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.
Derivation of the procedures that can be applied in evaluating two-time correlation function in terms of coherent-state propagator and corresponding Q-function is presented. On the basis that the involved functions are generally exponential in nature, obtaining the two-time second-order correlation function is essentially claimed to be reduced to carrying out relatively simple integrations. Fundamentally, the time dependence of the operators is transferred to the density operator. Moreover, manipulation in reordering the operators is performed by applying the usual trace operation. With all details, it is basically observed that the two-time correlation can be readily determined once the pertinent coherent-state propagator or Q-function is known. Since working with c-number equation is far more handy than the associated operator equation, it is expected that the results derived in this contribution can aid in easing the otherwise involving mathematical rigor.
Calculations at finite temperatures are fundamental in different scientific fields, from nuclear physics to condensed matter. Evolution in imaginary time is a prominent classical technique for preparing thermal states of quantum systems. We propose a new quantum algorithm that prepares thermal states based on the quantum imaginary time propagation method, using a diluted operator with ancilla qubits to overcome the non-unitarity nature of the imaginary time operator. The presented method is the first that allows us to obtain the correct thermal density matrix on a general quantum processor for a generic Hamiltonian. We prove its reliability in the actual quantum hardware computing thermal properties for two and three neutron systems.
Rigorous simulations of quantum tunneling dynamics in model systems with up to 20 coupled degrees of freedom are reported. The simulations implement an extension of the recently developed matching-pursuit/split-operator Fourier-transform method to complex-valued coherent-state representations. The resulting method recursively applies the time-evolution operator, as defined by the Trotter expansion to second order accuracy, in dynamically adaptive coherent-state representations generated by an approach that combines the matching-pursuit algorithm with a gradient-based optimization method.