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 Physical Chemistry Letters

Essential vs Removable Bath Anharmonicity: Path Integral Results with Model Electronic-Vibrational Baths

Nancy Makri

Published Sep 18, 2024• 8 citations
Publisher

Abstract

Fully quantum mechanical real-time path integral calculations on a two-level system (TLS) interacting with a simple or composite anharmonic bath are presented. The bath consists of a few or many units, each comprising two or three electronic states coupled to secondary harmonic vibrational baths. The primitive harmonic bath approximation leads to qualitatively wrong system dynamics at high temperatures, while the effective harmonic bath mapping becomes exact in the macroscopic bath limit. Significant deviations are observed when the number of bath units is small, including skewing and blue shifts in the population oscillations and enhanced or suppressed coherence. These effects reflect the essential anharmonic character of the bath and cannot be captured by harmonic bath treatments.

DOI
10.1021/acs.jpclett.4c02213

Related papers

Two-Tier Modular Anharmonic Small Matrix Path Integral with Composite Spin-Boson Baths

Nancy Makri

The anharmonic small matrix path integral (anh-SMatPI) algorithm is adapted to composite environments, where each anharmonic bath unit couples to the system and to its own harmonic bath. The propagation matrices are obtained from numerical evaluation of the composite influence functional, which is based on separate iterative path integral calculations for each bath module. The numerically exact nature of the algorithm is demonstrated on a two-level system (TLS) coupled to dissipative TLS baths of one, two, and 50 units.

Journal of Chemical Theory and Computation•Sep 20, 2024•6 citations

Iterative path integral calculations of non-linear spectra and entanglement dynamics

Mohammad Mehdi Sahrapour

We review an iterative path integral method that allows the computation of exact, long-time dynamics of small systems interacting with a dissipative bath. The method takes advantage of the finite memory of large environments at finite temperature to allow an iterative evaluation of the dynamics, thereby replacing an exponential scaling in simulation time with a linear scaling. This method is applied to calculate the dynamics of two model systems. The first consists of two-qubits interacting with a common bath. In this case we observe a variety of entanglement effects. If the qubits are initially separable, through an indirect coupling, the bath can create steady state entanglement between the qubits. This effect is due to the existence of a decoherence-free subspace as a result of the form of the qubit Hamiltonian and system-bath coupling. Entanglement created by the bath is shown to decrease with increasing temperature and system-bath coupling strength. However large system-bath coupling causes a faster increase in the entanglement. Initially entangled qubits lose their entanglement as a result of interactions with the bath, an effect that is heightened at higher temperatures. Direct coupling between the qubits is shown to slow the decay of entanglement and preserve some entanglement at long times; however at high temperatures this steady state entanglement becomes negligible. The second system we consider is vibrational degree of freedom coupled to a bath of harmonic oscillators or two-level systems. We compute four-time correlation functions which are used to calculate response functions relevant to third order infrared or seventh order Raman experiments for harmonic, Morse, and quadratic-quartic potentials. Our calculations reveal the role of potential features (anharmonicity and eigenvalue spectrum), both on short and long time scales, on the response function. Further, thermal excitation causes dramatic changes in the appearance of the response function, introducing symmetry with respect to the main diagonal. Finally, coupling to harmonic dissipative baths leads to decay of the response function (primarily along the t3 direction) and a broadening of the peaks in its Fourier transform. At high temperatures two-level-system baths are less efficient in destroying coherence than harmonic baths of similar parameters.

•May 24, 2013•0 citations

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

Path Integral over Equivalence Classes for Quantum Dynamics with Static Disorder

Nancy Makri

An efficient, fully quantum mechanical, real-time path integral method for including the effects of static disorder in the dynamics of systems coupled to common or local harmonic baths is presented. Rather than performing a large number of demanding calculations for different realizations of the system Hamiltonian, the influence of the bath is captured through a single evaluation of the path sum by grouping the system paths into equivalence classes of fixed system amplitudes. The method is illustrated with several analytical and numerical examples that show a variety of nontrivial effects arising from the interplay among coherence, dissipation, thermal fluctuations, and geometric phases.

The Journal of Physical Chemistry Letters•Jan 31, 2024•7 citations

Parsing the Influence Functional: Harmonic Bath Mapping and Anharmonic Small Matrix Path Integral

Nancy Makri

The influence functional (IF) encodes all of the information required for calculating the dynamical properties of a system in contact with its environment. A direct and simple procedure is introduced for extracting from a few numerical evaluations of the IF, without computing time correlation functions or evaluating integrals, the parameters required for path integral calculations, either within or beyond the harmonic mapping, and for assessing the accuracy of the harmonic bath approximation. In addition, the small matrix decomposition of the path integral (SMatPI) is extended to anharmonic environments and the required matrices are constructed directly from the IF.

The Journal of Physical Chemistry Letters•Apr 19, 2024•10 citations

Dissipative quantum dynamics with the surrogate Hamiltonian approach. A comparison between spin and harmonic baths

David Gelman, Christiane P. Koch, Ronnie Kosloff

The dissipative quantum dynamics of an anharmonic oscillator coupled to a bath is studied with the purpose of elucidating the differences between the relaxation to a spin bath and to a harmonic bath. Converged results are obtained for the spin bath by the surrogate Hamiltonian approach. This method is based on constructing a system-bath Hamiltonian, with a finite but large number of spin bath modes, that mimics exactly a bath with an infinite number of modes for a finite time interval. Convergence with respect to the number of simultaneous excitations of bath modes can be checked. The results are compared to calculations that include a finite number of harmonic modes carried out by using the multiconfiguration time-dependent Hartree method of Nest and Meyer [J. Chem. Phys. 119, 24 (2003)]. In the weak coupling regime, at zero temperature and for small excitations of the primary system, both methods converge to the Markovian limit. When initially the primary system is significantly excited, the spin bath can saturate restricting the energy acceptance. An interaction term between bath modes that spreads the excitation eliminates the saturation. The loss of phase between two cat states has been analyzed and the results for the spin and harmonic baths are almost identical. For stronger couplings, the dynamics induced by the two types of baths deviate. The accumulation and degree of entanglement between the bath modes have been characterized. Only in the spin bath the dynamics generate entanglement between the bath modes.

arXiv:quant-ph/0402144•The Journal of Chemical Physics•Jun 23, 2004•51 citations
Fully quantum mechanical real-time path integral calculations on a two-level system (TLS) interacting with a simple or composite anharmonic bath are presented. The bath consists of a few or many units, each comprising two or three electronic states coupled to secondary harmonic vibrational baths. The primitive harmonic bath approximation leads to qualitatively wrong system dynamics at high temperatures, while the effective harmonic bath mapping becomes exact in the macroscopic bath limit. Significant deviations are observed when the number of bath units is small, including skewing and blue shifts in the population oscillations and enhanced or suppressed coherence. These effects reflect the essential anharmonic character of the bath and cannot be captured by harmonic bath treatments.
The anharmonic small matrix path integral (anh-SMatPI) algorithm is adapted to composite environments, where each anharmonic bath unit couples to the system and to its own harmonic bath. The propagation matrices are obtained from numerical evaluation of the composite influence functional, which is based on separate iterative path integral calculations for each bath module. The numerically exact nature of the algorithm is demonstrated on a two-level system (TLS) coupled to dissipative TLS baths of one, two, and 50 units.
We review an iterative path integral method that allows the computation of exact, long-time dynamics of small systems interacting with a dissipative bath. The method takes advantage of the finite memory of large environments at finite temperature to allow an iterative evaluation of the dynamics, thereby replacing an exponential scaling in simulation time with a linear scaling. This method is applied to calculate the dynamics of two model systems. The first consists of two-qubits interacting with a common bath. In this case we observe a variety of entanglement effects. If the qubits are initially separable, through an indirect coupling, the bath can create steady state entanglement between the qubits. This effect is due to the existence of a decoherence-free subspace as a result of the form of the qubit Hamiltonian and system-bath coupling. Entanglement created by the bath is shown to decrease with increasing temperature and system-bath coupling strength. However large system-bath coupling causes a faster increase in the entanglement. Initially entangled qubits lose their entanglement as a result of interactions with the bath, an effect that is heightened at higher temperatures. Direct coupling between the qubits is shown to slow the decay of entanglement and preserve some entanglement at long times; however at high temperatures this steady state entanglement becomes negligible. The second system we consider is vibrational degree of freedom coupled to a bath of harmonic oscillators or two-level systems. We compute four-time correlation functions which are used to calculate response functions relevant to third order infrared or seventh order Raman experiments for harmonic, Morse, and quadratic-quartic potentials. Our calculations reveal the role of potential features (anharmonicity and eigenvalue spectrum), both on short and long time scales, on the response function. Further, thermal excitation causes dramatic changes in the appearance of the response function, introducing symmetry with respect to the main diagonal. Finally, coupling to harmonic dissipative baths leads to decay of the response function (primarily along the t3 direction) and a broadening of the peaks in its Fourier transform. At high temperatures two-level-system baths are less efficient in destroying coherence than harmonic baths of similar parameters.
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.
An efficient, fully quantum mechanical, real-time path integral method for including the effects of static disorder in the dynamics of systems coupled to common or local harmonic baths is presented. Rather than performing a large number of demanding calculations for different realizations of the system Hamiltonian, the influence of the bath is captured through a single evaluation of the path sum by grouping the system paths into equivalence classes of fixed system amplitudes. The method is illustrated with several analytical and numerical examples that show a variety of nontrivial effects arising from the interplay among coherence, dissipation, thermal fluctuations, and geometric phases.
The influence functional (IF) encodes all of the information required for calculating the dynamical properties of a system in contact with its environment. A direct and simple procedure is introduced for extracting from a few numerical evaluations of the IF, without computing time correlation functions or evaluating integrals, the parameters required for path integral calculations, either within or beyond the harmonic mapping, and for assessing the accuracy of the harmonic bath approximation. In addition, the small matrix decomposition of the path integral (SMatPI) is extended to anharmonic environments and the required matrices are constructed directly from the IF.
The dissipative quantum dynamics of an anharmonic oscillator coupled to a bath is studied with the purpose of elucidating the differences between the relaxation to a spin bath and to a harmonic bath. Converged results are obtained for the spin bath by the surrogate Hamiltonian approach. This method is based on constructing a system-bath Hamiltonian, with a finite but large number of spin bath modes, that mimics exactly a bath with an infinite number of modes for a finite time interval. Convergence with respect to the number of simultaneous excitations of bath modes can be checked. The results are compared to calculations that include a finite number of harmonic modes carried out by using the multiconfiguration time-dependent Hartree method of Nest and Meyer [J. Chem. Phys. 119, 24 (2003)]. In the weak coupling regime, at zero temperature and for small excitations of the primary system, both methods converge to the Markovian limit. When initially the primary system is significantly excited, the spin bath can saturate restricting the energy acceptance. An interaction term between bath modes that spreads the excitation eliminates the saturation. The loss of phase between two cat states has been analyzed and the results for the spin and harmonic baths are almost identical. For stronger couplings, the dynamics induced by the two types of baths deviate. The accumulation and degree of entanglement between the bath modes have been characterized. Only in the spin bath the dynamics generate entanglement between the bath modes.