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

Physical Chemistry Chemical Physics

Density matrix and purity evolution in dissipative two-level systems: II. Relaxation

Sambarta Chatterjee, Nancy Makri

Published Jan 1, 2021• 13 citations
Publisher

Abstract

We investigate the time evolution of the reduced density matrix (RDM) and its purity in the dynamics of a two-level system coupled to a dissipative harmonic bath, when the system is initially placed in one of its eigenstates. We point out that the symmetry of the initial condition confines the motion of the RDM elements to a one-dimensional subspace and show that the purity always goes through its maximally mixed value at some time during relaxation, but subsequently recovers and (under low-temperature, weakly dissipative conditions) can rise to values that approach unity. These behaviors are quantified through accurate path integral calculations. Under low-temperature, weakly dissipative conditions, we observe unusual, nonmonotonic population dynamics when the two-level system is initially placed in its ground state. We also analyze the origin of the system-bath interactions responsible for the nonmonotonic behavior of purity during relaxation. Our results show that classical dephasing processes arising from site level fluctuations lead to a monotonic decay of purity, and that the quantum mechanical decoherence events associated with spontaneous phonon emission are responsible for the subsequent recovery of purity. Last, we show that coupling with a low-temperature bath can purify a mixed two-level system. In the case of the maximally mixed initial RDM, the purity increases monotonically even during short time.

DOI
10.1039/d0cp05528j

Related papers

Density matrix and purity evolution in dissipative two-level systems: I. Theory and path integral results for tunneling dynamics

Sambarta Chatterjee, Nancy Makri

The time evolution of the purity (the trace of the square of the reduced density matrix) and von Neumann entropy in a symmetric two-level system coupled to a dissipative harmonic bath is investigated through analytical arguments and accurate path integral calculations on simple models and the singly excited bacteriochlorophyll dimer. A simple theoretical analysis establishes bounds and limiting behaviors. The contributions to purity from a purely incoherent term obtained from the diagonal elements of the reduced density matrix, a term associated with the difference of the two eigenstate populations, and a third term related to the square of the time derivative of a site population, are discussed in various regimes. In the case of tunneling dynamics from a localized initial condition, the complex interplay among these contributions leads to the recovery of purity under low-temperature, weakly dissipative conditions. Memory effects from the bath are found to play a critical role to the dynamics of purity. It is shown that the strictly quantum mechanical decoherence process associated with spontaneous phonon emission is responsible for the long-time recovery of purity. These analytical and numerical results show clearly that the loss of quantum coherence during the evolution toward equilibrium does not necessarily imply the decay of purity, and that the time scales relevant to these two processes may be entirely different.

Physical Chemistry Chemical Physics•Jan 1, 2021•10 citations

Recovery of Purity in Dissipative Tunneling Dynamics

Sambarta Chatterjee, Nancy Makri

The time evolution of purity for an initially localized state of a symmetric two-level system coupled to a dissipative bath is investigated using numerically exact real-time path integral methods. With strong system-bath coupling and high temperature, the purity decays monotonically to its fully mixed value, with a short-time Gaussian behavior, which is subsequently followed by exponential evolution. However, under low-temperature and weak coupling conditions, a substantial recovery of purity is observed. A simple theoretical analysis reveals three contributions that correspond to a completely incoherent, eigenstate population difference and rate terms. The last two of these terms can counter the early drop of purity and are responsible for its rebound. These findings caution against using purity as a measure of decoherence in the dynamics of quantum dissipative systems.

The Journal of Physical Chemistry Letters•Sep 25, 2020•13 citations

Non-Markovian qubit dynamics in a thermal field bath: Relaxation, decoherence, and entanglement

Sanjiv Shresta, Charis Anastopoulos, Adrian A. Drǎgulescu +1 more

We study the non-Markovian dynamics of a qubit made up of a two-level atom interacting with an electromagnetic field (EMF) initially at finite temperature. Unlike most earlier studies where the bath is assumed to be fixed, we study the complete evolution of the combined qubit-EMF system, thus allowing for the coherent backaction from the bath on the qubit and the qubit on the bath in a self-consistent manner. In this way we can see the development of quantum correlations and entanglement between the system and its environment, and how that affects the decoherence and relaxation of the system. We find nonexponential decay for both the diagonal and nondiagonal matrix elements of the qubit's reduced density matrix in the pointer basis. The former shows the qubit relaxing to thermal equilibrium with the bath, while the latter shows the decoherence rate beginning at the usually predicted thermal rate, but changing to the zero-temperature value as the qubit and bath become entangled. The decoherence and relaxation rates are comparable, as in the zero-temperature case. On the entanglement of a qubit with the EMF we calculated the fidelity and the von Neumann entropy, which is a measure of the purity of the density matrix. The present more accurate non-Markovian calculations predict lower loss of fidelity and purity as compared with the Markovian results. Generally speaking, with the inclusion of quantum correlations between the qubit and its environment, the non-Markovian processes tend to slow down the drive of the system to equilibrium, prolonging the decoherence and better preserving the fidelity and purity of the system.

arXiv:quant-ph/0408084•Physical Review A•Feb 22, 2005•59 citations

Characterization of a qubit in presence of dissipation and external driving

Francesco Nesi

In this PhD thesis we have been concerned with dissipative and decoherence processes which take place specifically for quantum systems in a dissipative bistable potential. This potential describes a Hamiltonian system whose spectrum is composed by a ladder of levels, grouped as doublets under the potential barrier and being almost equally spaced above it. The system interacts with a thermal environment, fully characterized by the so-called bath spectral density. The dynamical quantity we looked at is the so-called reduced density matrix (RDM), obtained by tracing out the full density matrix the environmental degrees of freedom. The matrix elements for the RDM are exactly evaluated by means of the Feynman-Vernon real-time path-integral technique. Depending on which regime of parameters we are interested in, i.e. temperature and coupling strength with the reservoir, a few or several energy levels contribute to the dynamics. In particular, a two-level system (TLS) is first analyzed, corresponding to the regime where only the lowest doublet is energetically accessible. We discussed in the beginning the case of a TLS interacting with Ohmic and super-Ohmic baths and then the situation in which a symmetric TLS is assumed to be coupled with an external detector, like a dc-SQUID, the latter being coupled to an Ohmic environment. This model provides a better description of a real qubit, as frequently discussed in various experiments. In the former case, we developed a theory called Weakly-Interacting Blip Approximation (WIBA), which is nonperturbative in the coupling strength as well as in the tunneling frequency and gives reliable results over a wide range of temperatures and coupling strength to the thermal environment. In the latter case, the dynamics was evaluated upon starting from the Non-Interacting Blip Approximation (NIBA) and approximating the NIBA kernels up to first order in the detector-bath coupling strength. Analytical results, valid also at finite detuning, were obtained and we showed that our theory correctly reproduces the onset of multiple oscillation frequencies in the dynamics due to the entanglement between TLS and detector. Finally, we discussed the more general case of an M-level system coupled with an Ohmic environment. This situation describes the case of a bistable potential where the temperature is high enough that higher-lying energy levels are being populated and the system cannot be reduced to a simple TLS. Many are the experimental situations which can be described by our theory. We developed a novel approximation called Weakly-Interacting VR-Blip Approximation (VR-WIBA) for the dynamics in a driven dissipative bistable potential. In the regimes of low to moderate temperatures and coupling strengths, we exploited a separation of time scales between tunneling and vibrational relaxation events and we could derive analytical kernels for the intra- and inter-well dynamics.

•Jan 1, 2007•1 citation

Quantum dynamics of a two-level system in a structured environment: Numerical study beyond perturbation theory

Johanna Charlotte Escher, Joachim Ankerhold

The dissipative quantum dynamics of a two-level system interacting with a structured reservoir consisting of a damped harmonic mode is investigated by means of the numerically exact path integral Monte Carlo method. This approach provides benchmark results in a broad range of parameter space, in particular, in those domains which are not accessible by approximate methods and alternative numerical schemes, that is, strong coupling between system and harmonic mode and from low to high temperatures. For low temperatures the numerical data are quantitatively in agreement with the noninteracting blip approximation only in the regimes of very weak and very strong coupling. It turns out that the entangled dynamics of the two-level system and the harmonic mode is relatively robust so that its signatures are observable up to relatively strong friction and high temperatures. Nonequilibrium initial preparations of the reservoir with respect to the initial state of the system give rise, for strong interactions, to a stepwise decay of the population, thus displaying coherent dynamics of the bath. The impact of an additional Ohmic bath coupled directly to the two-level system is studied as well, including the case where both reservoirs carry different temperatures.

Physical Review A•Mar 30, 2011•12 citations

Time Evolution of Bath Properties in Spin-Boson Dynamics

Sohang Kundu, Nancy Makri

The dynamical behaviors of a two-level system (TLS) coupled to a harmonic dissipative bath has been studied extensively using a variety of analytical and numerical methods. The focus of the vast majority of these studies has been on the properties of the TLS, averaged with respect to the bath degrees of freedom. In this work, we use real-time path integral methods to probe the behavior of select bath degrees of freedom during the dynamics of a symmetric two-level system (TLS) coupled to a dissipative bath by calculating system-bath densities (SBD) and coordinate expectation values. Overall, the SBD motion on each diabatic state is simpler than the motion of the total density. In the weak coupling regime, which characterizes the parameters of oscillators that comprise such a bath, the SBD on each TLS state remains primarily compact and Gaussian-like, such that its peak is well characterized by the mode expectation value. In the absence of a dissipative environment, nonadiabatic density depletion leads to spikes in coordinate expectation values. The evolution of the SBD peak trajectory for two discrete modes exhibits Lissajous patterns with frequency-dependent shapes that strongly resemble classical trajectory motion on a torus. These patterns become more complex when the coupling of the mode to the TLS is increased outside of this regime, leading to persistent small amplitude oscillations in the TLS populations characterized by a very slow decay and SBD trajectories that exhibit behaviors reminiscent of chaotic classical systems. Indirect coupling to a dissipative bath has a stabilizing effect on the dynamics, eliminating spikes, synchronizing the SBD motion on the two diabatic states and regularizing the SBD trajectory to simple rectangular Lissajous-like shapes with a slowly shrinking boundary, regardless of the mode frequencies.

The Journal of Physical Chemistry B•Jul 20, 2021•11 citations
We investigate the time evolution of the reduced density matrix (RDM) and its purity in the dynamics of a two-level system coupled to a dissipative harmonic bath, when the system is initially placed in one of its eigenstates. We point out that the symmetry of the initial condition confines the motion of the RDM elements to a one-dimensional subspace and show that the purity always goes through its maximally mixed value at some time during relaxation, but subsequently recovers and (under low-temperature, weakly dissipative conditions) can rise to values that approach unity. These behaviors are quantified through accurate path integral calculations. Under low-temperature, weakly dissipative conditions, we observe unusual, nonmonotonic population dynamics when the two-level system is initially placed in its ground state. We also analyze the origin of the system-bath interactions responsible for the nonmonotonic behavior of purity during relaxation. Our results show that classical dephasing processes arising from site level fluctuations lead to a monotonic decay of purity, and that the quantum mechanical decoherence events associated with spontaneous phonon emission are responsible for the subsequent recovery of purity. Last, we show that coupling with a low-temperature bath can purify a mixed two-level system. In the case of the maximally mixed initial RDM, the purity increases monotonically even during short time.
The time evolution of the purity (the trace of the square of the reduced density matrix) and von Neumann entropy in a symmetric two-level system coupled to a dissipative harmonic bath is investigated through analytical arguments and accurate path integral calculations on simple models and the singly excited bacteriochlorophyll dimer. A simple theoretical analysis establishes bounds and limiting behaviors. The contributions to purity from a purely incoherent term obtained from the diagonal elements of the reduced density matrix, a term associated with the difference of the two eigenstate populations, and a third term related to the square of the time derivative of a site population, are discussed in various regimes. In the case of tunneling dynamics from a localized initial condition, the complex interplay among these contributions leads to the recovery of purity under low-temperature, weakly dissipative conditions. Memory effects from the bath are found to play a critical role to the dynamics of purity. It is shown that the strictly quantum mechanical decoherence process associated with spontaneous phonon emission is responsible for the long-time recovery of purity. These analytical and numerical results show clearly that the loss of quantum coherence during the evolution toward equilibrium does not necessarily imply the decay of purity, and that the time scales relevant to these two processes may be entirely different.
The time evolution of purity for an initially localized state of a symmetric two-level system coupled to a dissipative bath is investigated using numerically exact real-time path integral methods. With strong system-bath coupling and high temperature, the purity decays monotonically to its fully mixed value, with a short-time Gaussian behavior, which is subsequently followed by exponential evolution. However, under low-temperature and weak coupling conditions, a substantial recovery of purity is observed. A simple theoretical analysis reveals three contributions that correspond to a completely incoherent, eigenstate population difference and rate terms. The last two of these terms can counter the early drop of purity and are responsible for its rebound. These findings caution against using purity as a measure of decoherence in the dynamics of quantum dissipative systems.
We study the non-Markovian dynamics of a qubit made up of a two-level atom interacting with an electromagnetic field (EMF) initially at finite temperature. Unlike most earlier studies where the bath is assumed to be fixed, we study the complete evolution of the combined qubit-EMF system, thus allowing for the coherent backaction from the bath on the qubit and the qubit on the bath in a self-consistent manner. In this way we can see the development of quantum correlations and entanglement between the system and its environment, and how that affects the decoherence and relaxation of the system. We find nonexponential decay for both the diagonal and nondiagonal matrix elements of the qubit's reduced density matrix in the pointer basis. The former shows the qubit relaxing to thermal equilibrium with the bath, while the latter shows the decoherence rate beginning at the usually predicted thermal rate, but changing to the zero-temperature value as the qubit and bath become entangled. The decoherence and relaxation rates are comparable, as in the zero-temperature case. On the entanglement of a qubit with the EMF we calculated the fidelity and the von Neumann entropy, which is a measure of the purity of the density matrix. The present more accurate non-Markovian calculations predict lower loss of fidelity and purity as compared with the Markovian results. Generally speaking, with the inclusion of quantum correlations between the qubit and its environment, the non-Markovian processes tend to slow down the drive of the system to equilibrium, prolonging the decoherence and better preserving the fidelity and purity of the system.
In this PhD thesis we have been concerned with dissipative and decoherence processes which take place specifically for quantum systems in a dissipative bistable potential. This potential describes a Hamiltonian system whose spectrum is composed by a ladder of levels, grouped as doublets under the potential barrier and being almost equally spaced above it. The system interacts with a thermal environment, fully characterized by the so-called bath spectral density. The dynamical quantity we looked at is the so-called reduced density matrix (RDM), obtained by tracing out the full density matrix the environmental degrees of freedom. The matrix elements for the RDM are exactly evaluated by means of the Feynman-Vernon real-time path-integral technique. Depending on which regime of parameters we are interested in, i.e. temperature and coupling strength with the reservoir, a few or several energy levels contribute to the dynamics. In particular, a two-level system (TLS) is first analyzed, corresponding to the regime where only the lowest doublet is energetically accessible. We discussed in the beginning the case of a TLS interacting with Ohmic and super-Ohmic baths and then the situation in which a symmetric TLS is assumed to be coupled with an external detector, like a dc-SQUID, the latter being coupled to an Ohmic environment. This model provides a better description of a real qubit, as frequently discussed in various experiments. In the former case, we developed a theory called Weakly-Interacting Blip Approximation (WIBA), which is nonperturbative in the coupling strength as well as in the tunneling frequency and gives reliable results over a wide range of temperatures and coupling strength to the thermal environment. In the latter case, the dynamics was evaluated upon starting from the Non-Interacting Blip Approximation (NIBA) and approximating the NIBA kernels up to first order in the detector-bath coupling strength. Analytical results, valid also at finite detuning, were obtained and we showed that our theory correctly reproduces the onset of multiple oscillation frequencies in the dynamics due to the entanglement between TLS and detector. Finally, we discussed the more general case of an M-level system coupled with an Ohmic environment. This situation describes the case of a bistable potential where the temperature is high enough that higher-lying energy levels are being populated and the system cannot be reduced to a simple TLS. Many are the experimental situations which can be described by our theory. We developed a novel approximation called Weakly-Interacting VR-Blip Approximation (VR-WIBA) for the dynamics in a driven dissipative bistable potential. In the regimes of low to moderate temperatures and coupling strengths, we exploited a separation of time scales between tunneling and vibrational relaxation events and we could derive analytical kernels for the intra- and inter-well dynamics.
The dissipative quantum dynamics of a two-level system interacting with a structured reservoir consisting of a damped harmonic mode is investigated by means of the numerically exact path integral Monte Carlo method. This approach provides benchmark results in a broad range of parameter space, in particular, in those domains which are not accessible by approximate methods and alternative numerical schemes, that is, strong coupling between system and harmonic mode and from low to high temperatures. For low temperatures the numerical data are quantitatively in agreement with the noninteracting blip approximation only in the regimes of very weak and very strong coupling. It turns out that the entangled dynamics of the two-level system and the harmonic mode is relatively robust so that its signatures are observable up to relatively strong friction and high temperatures. Nonequilibrium initial preparations of the reservoir with respect to the initial state of the system give rise, for strong interactions, to a stepwise decay of the population, thus displaying coherent dynamics of the bath. The impact of an additional Ohmic bath coupled directly to the two-level system is studied as well, including the case where both reservoirs carry different temperatures.
The dynamical behaviors of a two-level system (TLS) coupled to a harmonic dissipative bath has been studied extensively using a variety of analytical and numerical methods. The focus of the vast majority of these studies has been on the properties of the TLS, averaged with respect to the bath degrees of freedom. In this work, we use real-time path integral methods to probe the behavior of select bath degrees of freedom during the dynamics of a symmetric two-level system (TLS) coupled to a dissipative bath by calculating system-bath densities (SBD) and coordinate expectation values. Overall, the SBD motion on each diabatic state is simpler than the motion of the total density. In the weak coupling regime, which characterizes the parameters of oscillators that comprise such a bath, the SBD on each TLS state remains primarily compact and Gaussian-like, such that its peak is well characterized by the mode expectation value. In the absence of a dissipative environment, nonadiabatic density depletion leads to spikes in coordinate expectation values. The evolution of the SBD peak trajectory for two discrete modes exhibits Lissajous patterns with frequency-dependent shapes that strongly resemble classical trajectory motion on a torus. These patterns become more complex when the coupling of the mode to the TLS is increased outside of this regime, leading to persistent small amplitude oscillations in the TLS populations characterized by a very slow decay and SBD trajectories that exhibit behaviors reminiscent of chaotic classical systems. Indirect coupling to a dissipative bath has a stabilizing effect on the dynamics, eliminating spikes, synchronizing the SBD motion on the two diabatic states and regularizing the SBD trajectory to simple rectangular Lissajous-like shapes with a slowly shrinking boundary, regardless of the mode frequencies.