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

arXiv:gr-qc/9401016

How to calculate a decoherence matrix numerically and the microscopic mechanism for decoherence in the Caldeira-Leggett model

Hans-Jürgen Pohle

Published Jan 17, 1994• 0 citations
arXivHTMLPDF

Abstract

A central object in the interpretation of quantum mechanics of closed systems with decoherent histories is the decoherence matrix. But only for a very small number of models one is able to give explicit expressions for its elements. So numerical methods are required. Unfortunately the dimensions of this matrices are usually very high, which makes also a direct numerical calculation impossible. A solution of this problem would be given by a method which only calculates the dominant matrix elements. This includes to make a decision about the dominance of an element before it will be calculated. In this paper I will develop an algorithm that combines the numerical calculation of the elements of the decoherence matrix with a permanent estimation, so that finally the dominant elements will be calculated only. As an example I apply this procedure to the Caldeira- Leggett-modell.

DOI
10.48550/arxiv.gr-qc/9401016

Related papers

Quantum mechanics of history: The decoherence functional in quantum mechanics

Fay Dowker B, J. J. Halliwell

We study a formulation of quantum mechanics in which the central notion is that of a quantum-mechanical history---a sequence of events at a succession of times. The primary aim is to identify sets of ``decoherent'' (or ``consistent'') histories for the system. These are quantum-mechanical histories suffering negligible interference with each other, and, therefore, to which probabilities may be assigned. These histories may be found for a given system using the so-called decoherence functional. When the decoherence functional is exactly diagonal, probabilities may be assigned to the histories, and all probability sum rules are satisfied exactly. We propose a condition for approximate decoherence, and argue that it implies that most probability sum rules will be satisfied to approximately the same degree. We also derive an inequality bounding the size of the off-diagonal terms of the decoherence functional. We calculate the decoherence functional for some simple one-dimensional systems, with a variety of initial states. For these systems, we explore the extent to which decoherence is produced using two different types of coarse graining. The first type of coarse graining involves imprecise specification of the particle's position. The second involves coupling the particle to a thermal bath of harmonic oscillators and ignoring the details of the bath (the Caldeira-Leggett model). We argue that both types of coarse graining are necessary in general. We explicitly exhibit the degree of decoherence as a function of the temperature of the bath, and of the width to within which the particle's position is specified. We study the diagonal elements of the decoherence functional, representing the probabilities for the possible histories of the system. To the extent that the histories decohere, we show that the probability distributions are peaked about the classical histories of the system, with the distribution of their initial positions and momenta given by a smeared version of the Wigner function. We discuss this result in connection with earlier uses of the Wigner function in this context. We find that there is a certain amount of tension between the demands of decoherence and peaking about classical paths.

Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Aug 15, 1992•205 citations

Quantum decoherence in the Caldeira-Leggett model by the real-time path integral on a computer

Jun Nishimura, Hiromasa Watanabe

A bstract We propose first-principle calculations of an open system based on the real-time path integral formalism treating the environment as well as the system of our interest together on a computer. The sign problem that occurs in applying Monte Carlo methods can be overcome in general by using the so-called Lefschetz thimble method, which has been developed over the past decade. Here we focus on the Caldeira-Leggett model, which is well known, in particular, as a model of quantum decoherence. In this case, the calculation simplifies drastically since the path integral becomes Gaussian for typical initial conditions. The relevant saddle point, which is unique and complex, can be determined by solving a linear equation with a huge but sparse coefficient matrix, and the integration over the Lefschetz thimble can be performed analytically. Thus we obtain, without assumptions or approximations, the reduced density matrix after a long-time evolution, tracing out a large number of harmonic oscillators in the environment. In particular, we confirm the dependence of the decoherence time on the coupling constant and the temperature that has been predicted from the master equation in a certain parameter regime.

arXiv:2503.20699•Journal of High Energy Physics•Sep 24, 2025•2 citations

Numerical studies on real time evolution in quantum systems - Quantum decoherence and the beginning of the universe

Naoyuki Yamamori

Real-time path integrals provide a natural formulation of quantum dynamics also in situations where a Hamiltonian description is unavailable or impractical, but their direct evaluation is obstructed by severe phase cancellation from the oscillatory weight.A guiding theme of this thesis is that, after complexifying the integration variables, real-time path integrals are naturally organized by complex saddle points.In favorable cases the original contour can be replaced by steepest-descent manifolds (Lefschetz thimbles) attached to the relevant saddles, providing a controlled saddle-based description of real-time observables.In Part I, we study quantum decoherence in a closed-system real-time path-integral formulation: the total system (subsystem plus environment) evolves unitarily, and the environment is traced out only at the final time to obtain the reduced density matrix.For the Caldeira-Leggett model with Gaussian initial conditions, the discretized effective action becomes a quadratic form, so that the reduced density matrix is characterized by a complex saddle point and the thimble integral reduces to a finite-dimensional Gaussian integral.The complex saddle is obtained by solving a large but sparse linear system, allowing an exact evaluation of the reduced density matrix without relying on the usual master-equation approximations.We then extend this methodology to a three-dimensional cubic lattice toy model by separating center-of-mass and relative coordinates, treating the center-of-mass mode as the subsystem and the relative modes as an effective environment.This enables controlled numerical studies of how the suppression of off-diagonal matrix elements depends on the coupling, inverse temperature, and system size, and allows us to extract the early-time scaling behavior of decoherence diagnostics from the complex-saddle evaluation.In Part II, we investigate spacetime dynamics in the Lorentzian IKKT matrix model.To define a Lorentzian formulation that is not continuously connected to the Euclidean model by the standard contour deformation, we include a Lorentz invariant mass term with coefficient γ, for which the classical equations admit real, smooth, expanding solutions when γ > 0. Since the Lorentz group is noncompact, the path integral can contain a divergent volume from the boost directions; we remove this volume nonperturbatively by implementing a boost "gauge-fixing" that selects the frame minimizing tr(A 0 ) 2 , leading to the condition tr(A 0 A i ) = 0 with the associated Faddeev-Popov determinant.We then apply the complex Langevin method to the resulting complex-weight matrix integral.Time and spatial observables are defined from matrix configurations, with holomorphic time ordering implemented by auxiliary variables, and spatial structure characterized through block observables such as the extent R 2 (t) and the moment-of-inertia tensor.In the bosonic model we find approximately real time and no indication of spontaneous breaking of the spatial symmetry in the i parameter range studied.Including fermions, we control the singular-drift problem induced by the Pfaffian by introducing a fermionic mass regulator m f and an anisotropic mass deformation, which enhances fermionic effects while maintaining numerical stability.In this setup we observe a real, expanding spatial extent and a dynamical dominance of three spatial directions at late times, consistent with a (3 + 1)-dimensional expanding spacetime.We finally discuss the intended limiting procedure-taking the large-N limit first and then approaching the undeformed theory-and comment on future directions, including a comparison with partially supersymmetric deformations that preserve 16 supercharges.

•Mar 24, 2026•0 citations

First Principles Numerical Demonstration of Emergent Decoherent Histories

Philipp Strasberg, Teresa E. Reinhard, Joseph Schindler

Within the histories formalism the decoherence functional is a formal tool to investigate the emergence of classicality in isolated quantum systems, yet an explicit evaluation of it from first principles has not been reported. We provide such an evaluation for up to five-time histories based on exact numerical diagonalization of the Schrödinger equation. We find a robust emergence of decoherence for slow and coarse observables of a generic random matrix model and extract a finite-size scaling law by varying the Hilbert space dimension over 4 orders of magnitude. Specifically, we conjecture and observe an exponential suppression of coherent effects as a function of the particle number of the system. This suggests a solution to the preferred basis problem of the many-worlds interpretation (or the set selection problem of the histories formalism) within a minimal theoretical framework without relying on environmentally induced decoherence, quantum Darwinism, Markov approximations, low-entropy initial states, or ensemble averages. Published by the American Physical Society 2024

arXiv:2304.10258•Physical Review X•Oct 30, 2024•18 citations

Adapted Caldeira-Leggett model

Andreas Albrecht, Rose Baunach, Andrew Arrasmith

We present a variant of the Caldeira-Leggett (CL) model of a harmonic oscillator coupled to an environment. The CL model is a standard tool for studying the physics of decoherence. Our ``adapted Caldeira-Leggett'' (ACL) model is built in a finite Hilbert space which makes it suitable for numerical studies. Taking a numerical approach allows us to avoid the limitations of standard approximation schemes used with the CL model. We are able to evolve the ACL model in a fully reversible unitary manner, without the built-in time asymmetry and other assumptions that come with the master equation methods typically used. We have used the ACL model to study topics in the field of decoherence and einselection where the full unitary evolution is essential to our work. Those results (reported in companion papers) include an examination of the relationship between einselection and the arrow of time, and studies of the very earliest stages of einselection. This paper provides details about the ACL model and our numerical methods. Our numerical approach makes it straightforward to explore and plot any property of the physical system. Thus we believe the examples and illustrations we present here may provide a helpful resource for those wishing to improve their familiarity with standard decoherence results, as well as those looking to probe the underpinnings of our companion papers. We expect the ACL model will be a useful tool for exploring additional phenomena that cannot be studied using traditional approximation schemes.

arXiv:2105.14040•Physical Review Research•Jun 23, 2023•6 citations

Decoherence: a numerical study

Chris Nagele, Oliver Janssen, Matthew Kleban

Abstract We study quantum decoherence numerically in a system consisting of a relativistic quantum field theory coupled to a measuring device that is itself coupled to an environment. The measuring device and environment are treated as quantum, non-relativistic particles. We solve the Schrödinger equation for the wave function of this tripartite system using exact diagonalization. Although computational limitations on the size of the Hilbert space prevent us from exploring the regime where the device and environment consist of a truly macroscopic number of degrees of freedom, we nevertheless see clear evidence of decoherence: after tracing out the environment, the density matrix describing the system and measuring device evolves quickly towards a matrix that is close to diagonal in a subspace of pointer states. We measure the speed with which decoherence spreads in the relativistic quantum field theory for a range of parameters. We find that it is less than the speed of light but faster than the speed of the massive charges in the initial state.

arXiv:2010.04803•Journal of Physics A Mathematical and Theoretical•Feb 6, 2023•3 citations
A central object in the interpretation of quantum mechanics of closed systems with decoherent histories is the decoherence matrix. But only for a very small number of models one is able to give explicit expressions for its elements. So numerical methods are required. Unfortunately the dimensions of this matrices are usually very high, which makes also a direct numerical calculation impossible. A solution of this problem would be given by a method which only calculates the dominant matrix elements. This includes to make a decision about the dominance of an element before it will be calculated. In this paper I will develop an algorithm that combines the numerical calculation of the elements of the decoherence matrix with a permanent estimation, so that finally the dominant elements will be calculated only. As an example I apply this procedure to the Caldeira- Leggett-modell.
We study a formulation of quantum mechanics in which the central notion is that of a quantum-mechanical history---a sequence of events at a succession of times. The primary aim is to identify sets of ``decoherent'' (or ``consistent'') histories for the system. These are quantum-mechanical histories suffering negligible interference with each other, and, therefore, to which probabilities may be assigned. These histories may be found for a given system using the so-called decoherence functional. When the decoherence functional is exactly diagonal, probabilities may be assigned to the histories, and all probability sum rules are satisfied exactly. We propose a condition for approximate decoherence, and argue that it implies that most probability sum rules will be satisfied to approximately the same degree. We also derive an inequality bounding the size of the off-diagonal terms of the decoherence functional. We calculate the decoherence functional for some simple one-dimensional systems, with a variety of initial states. For these systems, we explore the extent to which decoherence is produced using two different types of coarse graining. The first type of coarse graining involves imprecise specification of the particle's position. The second involves coupling the particle to a thermal bath of harmonic oscillators and ignoring the details of the bath (the Caldeira-Leggett model). We argue that both types of coarse graining are necessary in general. We explicitly exhibit the degree of decoherence as a function of the temperature of the bath, and of the width to within which the particle's position is specified. We study the diagonal elements of the decoherence functional, representing the probabilities for the possible histories of the system. To the extent that the histories decohere, we show that the probability distributions are peaked about the classical histories of the system, with the distribution of their initial positions and momenta given by a smeared version of the Wigner function. We discuss this result in connection with earlier uses of the Wigner function in this context. We find that there is a certain amount of tension between the demands of decoherence and peaking about classical paths.
A bstract We propose first-principle calculations of an open system based on the real-time path integral formalism treating the environment as well as the system of our interest together on a computer. The sign problem that occurs in applying Monte Carlo methods can be overcome in general by using the so-called Lefschetz thimble method, which has been developed over the past decade. Here we focus on the Caldeira-Leggett model, which is well known, in particular, as a model of quantum decoherence. In this case, the calculation simplifies drastically since the path integral becomes Gaussian for typical initial conditions. The relevant saddle point, which is unique and complex, can be determined by solving a linear equation with a huge but sparse coefficient matrix, and the integration over the Lefschetz thimble can be performed analytically. Thus we obtain, without assumptions or approximations, the reduced density matrix after a long-time evolution, tracing out a large number of harmonic oscillators in the environment. In particular, we confirm the dependence of the decoherence time on the coupling constant and the temperature that has been predicted from the master equation in a certain parameter regime.
Real-time path integrals provide a natural formulation of quantum dynamics also in situations where a Hamiltonian description is unavailable or impractical, but their direct evaluation is obstructed by severe phase cancellation from the oscillatory weight.A guiding theme of this thesis is that, after complexifying the integration variables, real-time path integrals are naturally organized by complex saddle points.In favorable cases the original contour can be replaced by steepest-descent manifolds (Lefschetz thimbles) attached to the relevant saddles, providing a controlled saddle-based description of real-time observables.In Part I, we study quantum decoherence in a closed-system real-time path-integral formulation: the total system (subsystem plus environment) evolves unitarily, and the environment is traced out only at the final time to obtain the reduced density matrix.For the Caldeira-Leggett model with Gaussian initial conditions, the discretized effective action becomes a quadratic form, so that the reduced density matrix is characterized by a complex saddle point and the thimble integral reduces to a finite-dimensional Gaussian integral.The complex saddle is obtained by solving a large but sparse linear system, allowing an exact evaluation of the reduced density matrix without relying on the usual master-equation approximations.We then extend this methodology to a three-dimensional cubic lattice toy model by separating center-of-mass and relative coordinates, treating the center-of-mass mode as the subsystem and the relative modes as an effective environment.This enables controlled numerical studies of how the suppression of off-diagonal matrix elements depends on the coupling, inverse temperature, and system size, and allows us to extract the early-time scaling behavior of decoherence diagnostics from the complex-saddle evaluation.In Part II, we investigate spacetime dynamics in the Lorentzian IKKT matrix model.To define a Lorentzian formulation that is not continuously connected to the Euclidean model by the standard contour deformation, we include a Lorentz invariant mass term with coefficient γ, for which the classical equations admit real, smooth, expanding solutions when γ > 0. Since the Lorentz group is noncompact, the path integral can contain a divergent volume from the boost directions; we remove this volume nonperturbatively by implementing a boost "gauge-fixing" that selects the frame minimizing tr(A 0 ) 2 , leading to the condition tr(A 0 A i ) = 0 with the associated Faddeev-Popov determinant.We then apply the complex Langevin method to the resulting complex-weight matrix integral.Time and spatial observables are defined from matrix configurations, with holomorphic time ordering implemented by auxiliary variables, and spatial structure characterized through block observables such as the extent R 2 (t) and the moment-of-inertia tensor.In the bosonic model we find approximately real time and no indication of spontaneous breaking of the spatial symmetry in the i parameter range studied.Including fermions, we control the singular-drift problem induced by the Pfaffian by introducing a fermionic mass regulator m f and an anisotropic mass deformation, which enhances fermionic effects while maintaining numerical stability.In this setup we observe a real, expanding spatial extent and a dynamical dominance of three spatial directions at late times, consistent with a (3 + 1)-dimensional expanding spacetime.We finally discuss the intended limiting procedure-taking the large-N limit first and then approaching the undeformed theory-and comment on future directions, including a comparison with partially supersymmetric deformations that preserve 16 supercharges.
Within the histories formalism the decoherence functional is a formal tool to investigate the emergence of classicality in isolated quantum systems, yet an explicit evaluation of it from first principles has not been reported. We provide such an evaluation for up to five-time histories based on exact numerical diagonalization of the Schrödinger equation. We find a robust emergence of decoherence for slow and coarse observables of a generic random matrix model and extract a finite-size scaling law by varying the Hilbert space dimension over 4 orders of magnitude. Specifically, we conjecture and observe an exponential suppression of coherent effects as a function of the particle number of the system. This suggests a solution to the preferred basis problem of the many-worlds interpretation (or the set selection problem of the histories formalism) within a minimal theoretical framework without relying on environmentally induced decoherence, quantum Darwinism, Markov approximations, low-entropy initial states, or ensemble averages. Published by the American Physical Society 2024
We present a variant of the Caldeira-Leggett (CL) model of a harmonic oscillator coupled to an environment. The CL model is a standard tool for studying the physics of decoherence. Our ``adapted Caldeira-Leggett'' (ACL) model is built in a finite Hilbert space which makes it suitable for numerical studies. Taking a numerical approach allows us to avoid the limitations of standard approximation schemes used with the CL model. We are able to evolve the ACL model in a fully reversible unitary manner, without the built-in time asymmetry and other assumptions that come with the master equation methods typically used. We have used the ACL model to study topics in the field of decoherence and einselection where the full unitary evolution is essential to our work. Those results (reported in companion papers) include an examination of the relationship between einselection and the arrow of time, and studies of the very earliest stages of einselection. This paper provides details about the ACL model and our numerical methods. Our numerical approach makes it straightforward to explore and plot any property of the physical system. Thus we believe the examples and illustrations we present here may provide a helpful resource for those wishing to improve their familiarity with standard decoherence results, as well as those looking to probe the underpinnings of our companion papers. We expect the ACL model will be a useful tool for exploring additional phenomena that cannot be studied using traditional approximation schemes.
Abstract We study quantum decoherence numerically in a system consisting of a relativistic quantum field theory coupled to a measuring device that is itself coupled to an environment. The measuring device and environment are treated as quantum, non-relativistic particles. We solve the Schrödinger equation for the wave function of this tripartite system using exact diagonalization. Although computational limitations on the size of the Hilbert space prevent us from exploring the regime where the device and environment consist of a truly macroscopic number of degrees of freedom, we nevertheless see clear evidence of decoherence: after tracing out the environment, the density matrix describing the system and measuring device evolves quickly towards a matrix that is close to diagonal in a subspace of pointer states. We measure the speed with which decoherence spreads in the relativistic quantum field theory for a range of parameters. We find that it is less than the speed of light but faster than the speed of the massive charges in the initial state.