We apply a full-polaron master equation and a weak-coupling non-Markovian master equation to describe the steady-state time-averaged properties of a driven two-level system, an electron coherently tunneling between double quantum dots (DQDs), interacting with a bosonic phonon bath. Comparing the results obtained using these two master equations with those from a recent DQD experiment and its corresponding weak-coupling theoretical method, we find that the original parameter set used in the experiment and theoretical method is not in the weak-coupling parameter regime. By using the full-polaron master equation with a slight adjustment on only the value of the interdot separation in the original experimental parameter set, we find that a reasonable fit to the experimentally measured time-averaged steady-state population data can be achieved. The adjusted interdot separation is within the possible values allowed by the geometry of the surface gates that define the DQD in the experiment. Our full-polaron equation approach does not require the special renormalization scheme employed in their weak-coupling theoretical method, and can still describe the experimental results of driving-induced phonon-enhanced steplike shoulder behaviors in the experiment. This demonstrates that the full-polaron master equation approach is a correct and efficient tool to describe the steady-state properties of a driven spin-boson model in the case of strong system-environment coupling.
We apply a full-polaron master equation and a weak-coupling non-Markovian master equation to describe the steady-state time-averaged properties of a driven two-level system, an electron coherently tunneling between double quantum dots (DQDs), interacting with a bosonic phonon bath. Comparing the results obtained using these two master equations with those from a recent DQD experiment and its corresponding weak-coupling theoretical method, we find that the original parameter set used in the experiment and theoretical method is not in the weak-coupling parameter regime. By using the full-polaron master equation with a slight adjustment on only the value of the interdot separation in the original experimental parameter set, we find that a reasonable fit to the experimentally measured time-averaged steady-state population data can be achieved. The adjusted interdot separation is within the possible values allowed by the geometry of the surface gates that define the DQD in the experiment. Our full-polaron equation approach does not require the special renormalization scheme employed in their weak-coupling theoretical method, and can still describe the experimental results of driving-induced phonon-enhanced steplike shoulder behaviors in the experiment. This demonstrates that the full-polaron master equation approach is a correct and efficient tool to describe the steady-state properties of a driven spin-boson model in the case of strong system-environment coupling.
We investigate the effects of initial states on the dynamics of a quantum spin-boson model using an efficient variational polaron master equation approach. The transformation of a factorized system-environment initial state into the variational polaron frame introduces inhomogeneous terms in the master equation, essential for accurately describing non-Markovian dynamics beyond the weak coupling regime. To efficiently account for these effects, we formulate the dynamics in an extended auxiliary Liouville space, converting the resulting integro-differential equations into coupled ordinary differential equations. We further introduce a simplified variational polaron technique for computing the time evolution of expectation values of general observables, including those that do not commute with the polaron transformation, overcoming a major limitation of previous approaches and enabling the unified treatment of coherence and population dynamics. The effects of transformed initial states generated from arbitrary factorized initial states are systematically characterized through expectation values and trace-distance analysis, providing insight into the associated bath-relaxation dynamics over a broad range of tunneling amplitudes, temperatures, and coupling strengths. Comparison with the weak-coupling master equation reveals a delay phenomenon associated with multi-bath-quanta processes, pronounced at intermediate coupling strengths and tunneling amplitudes. We further show that coherent superposition initial states induce stronger nonequilibrium effects than localized states, with inhomogeneous terms playing a crucial role in the short-time dynamics. The proposed formalism applies to both undriven and resonantly driven spin-boson models, providing a computationally efficient framework for open quantum-system dynamics with arbitrary factorized initial states and general physical observables.
We study two anisotropically interacting spins coupled to optical phonons; we restrict our analysis to the regime of strong coupling to the environment, to the antiadiabatic region, and to the subspace with zero value for STz (the z component of the total spin). In the case where each spin is coupled to a different phonon bath, we assume that the system and the environment are initially uncorrelated (and form a simply separable state) in the polaronic frame of reference. By analyzing the polaron dynamics through a non-Markovian quantum master equation, we find that the system manifests a small amount of decoherence that decreases both with increasing nonadiabaticity and with enhancing strength of coupling; whereas, under the Markovian approximation, the polaronic system exhibits a decoherence-free behavior. For the situation where both spins are coupled to the same phonon bath, we also show that the system is decoherence-free in the subspace where STz is fixed. To suppress decoherence through quantum control, we employ a train of $\ensuremath{\pi}$ pulses and demonstrate that unitary evolution of the system can be retained. We propose realization of a weakly decohering charge qubit from an electron in an oxide-based (tunnel-coupled) double-quantum-dot system.
We develop a versatile master equation approach to describe the nonequilibrium dynamics of a two-level system in contact with a bosonic environment, which allows for the exploration of a wide range of parameter regimes within a single formalism. As an experimentally relevant example, we apply this technique to the study of excitonic Rabi rotations in a driven quantum dot, and compare its predictions to the numerical Feynman integral approach. We find excellent agreement between the two methods across a generally difficult range of parameters. In particular, the variational master equation technique captures effects usually considered to be nonperturbative, such as multiphonon processes and bath-induced driving renormalization, and can give reliable results even in regimes in which previous master equation approaches fail.
We derive a Born-Markov master equation describing the dissipation induced by a bath of lossy but coherently driven two-level systems (TLSs) coupled to a bosonic system via Jaynes-Cummings interaction. We analytically derive all the master-equation rates. We characterize these rates for the particular case of a single-mode system coupled to identical TLSs. We study the steady state of the system and its exotic properties stemming from the nonthermal stationary state of the driven TLS bath. These properties include dissipative amplification, bath-induced linear instability, and both coherent and dissipative squeezing. The master equation is valid for arbitrarily strong TLS driving and it can be generalized to include multilevel systems or other system-bath interaction terms, among others. Our work provides a tool to study and characterize TLS-induced decoherence, a key limiting factor in quantum technological devices based on, for instance, superconducting circuits, magnonic systems, or quantum acoustics.
While the quantum regression theorem (QRT) is the standard tool for computing multi-time correlation functions in open quantum systems, it relies on system-bath separability and an environment that remains in equilibrium, assumptions that are violated once dynamical correlations develop. Using the projection operator formalism, we derive an extension to the QRT that explicitly incorporates these correlation-induced corrections, enabling us to capture multi-time correlation functions in non-Markovian regimes. We apply this framework to the variational polaron master equation for the spin-boson model in ohmic and super-ohmic regimes, where the polaron transformation mixes system-bath degrees of freedom to produce a non-thermal effective environment. Benchmarking against numerically exact tensor-network simulations, we demonstrate quantitative agreement for single- and two-time observables, including linear-response spectra, even at strong coupling. Our approach broadens the reach of analytic master equations to multi-time observables in strong-coupling and non-Markovian regimes, where environmental memory effects and system-bath correlations are crucial.
Understanding the dynamics of open quantum systems in strong coupling and non-Markovian regimes remains a formidable theoretical challenge. One popular and well-established method of approximation in these circumstances is provided by the polaron master equation (PME). In this work, we re-evaluate and extend the validity of the PME to capture the impact of non-Markovian polaron dressing, induced by non-equilibrium open system dynamics. By comparing with numerically exact techniques, we confirm that while the standard PME successfully predicts the dynamics of system observables that commute with the polaron transformation (e.g., populations in the Pauli z-basis), it can struggle to fully capture those that do not (e.g., coherences). This limitation stems from the mixing of system and environment degrees of freedom inherent to the polaron transformation, which affects the accuracy of calculated expectation values within the polaron frame. Employing the Nakajima-Zwanzig projection operator formalism, we introduce correction terms that provide an accurate description of observables that do not commute with the transformation. We demonstrate the significance of the correction terms in two cases, the canonical spin-boson model and a dissipative time-dependent Landau-Zener protocol, where they are shown to impact the system dynamics on both short and long timescales.