---
title: Driven Caldeira-Leggett Model Overview
url: https://www.emergentmind.com/topics/driven-caldeira-leggett-model
type: topic
---

# Driven Caldeira-Leggett Model Overview

The driven Caldeira-Leggett model denotes a family of open-quantum-system constructions built on the Caldeira-Leggett system-plus-bath framework, in which time dependence enters through direct forcing of the distinguished system coordinate, explicit driving of the bath coordinates, or non-equilibrium preparation of the reservoir. Taken together, the literature suggests that the term does not identify a single universally fixed formalism: in some works it means an externally driven system coupled to a harmonic bath, in others a bath that is itself driven by the external field, and in condensed-phase spectroscopy it can refer more loosely to the dissipative, noise-driven dynamics generated by the Caldeira-Leggett Hamiltonian and its generalized Langevin reduction [1807.01557] [1712.06397] [1601.04470].

## 1. Standard framework and the meaning of “driving”

The original Caldeira-Leggett construction describes a distinguished degree of freedom coupled linearly to many harmonic oscillators. In the formulation emphasized in the historical review of Caldeira’s work, the system-plus-environment Lagrangian is  
\[
L = \frac{1}{2}M\dot q^2 - V(q) + \frac{1}{2}\sum_j\left(m_j\dot x_j^2-m_j\omega_j^2x_j^2\right)-\sum_jF_j(q)x_j+\sum_j\frac{F_j^2(q)}{2m_j\omega_j^2},
\]
with bilinear coupling \(F_j(q)=c_j q\) in the standard model. The last term is the counter-term, whose stated role is to cancel the bath-induced renormalization of the bare system potential so that the dressed potential in the Hamiltonian is the intended physical one [2604.26524].

A Hamiltonian form used in driven-bath studies is  
\[
H_{\rm tot}(t)=H_S+H_R+H_{\rm int}+H_{\rm ext}(t),
\]
with
\[
H_S=\frac{p^2}{2M}+V(q), \qquad
H_R=\sum_n \left(\frac{p_n^2}{2m_n}+\frac{1}{2}m_n\omega_n^2 x_n^2\right),
\]
\[
H_{\rm int}=-q\sum_n c_n x_n+q^2\sum_n \frac{c_n^2}{2m_n\omega_n^2},
\]
and spectral density
\[
J(\omega)=\frac{\pi}{2}\sum_n \frac{c_n^2}{m_n\omega_n}\delta(\omega-\omega_n).
\]
Within this architecture, “driving” can be introduced in structurally distinct ways: by adding an explicit external term acting on the system coordinate, by forcing the bath coordinates themselves, or by engineering the bath state so that it acts as a source of deterministic or stochastic work [1807.01557] [2405.00215].

This distinction matters because the reduced system does not merely feel a local force. Once the bath is eliminated, the effective dynamics contains friction, fluctuations, and, in general, nonlocal memory kernels. A recurrent theme across the literature is that driving the environment is not equivalent to simply increasing the direct drive on the system; the bath has its own dynamics and mediates retarded response [1807.01557].

## 2. Explicitly driven baths and retarded effective forcing

A direct extension of the standard model is obtained by coupling the same external force \(F(t)\) both to the system coordinate and to each bath mode:
\[
H_{\rm ext}(t)= -\left(d_0 q+\sum_n d_n x_n\right)F(t).
\]
The bath equations of motion then read
\[
\dot x_n(t)=\frac{p_n(t)}{m_n},\qquad
\dot p_n(t)=-m_n\omega_n^2 x_n(t)+c_n q(t)+d_nF(t),
\]
so the drive enters the bath oscillator dynamics as a direct force term [1807.01557].

Eliminating the bath yields the exact operator equation
\[
M\ddot q(t)+M\int_0^t ds\,\gamma(t-s)\dot q(s)+\frac{\partial V(q(t))}{\partial q(t)}
= \xi(t)+d_0F(t)+\int_0^t ds\,\lambda(t-s)F(s),
\]
with the usual friction kernel
\[
\gamma(t)=\frac{1}{M}\sum_n \frac{c_n^2}{m_n\omega_n^2}\cos(\omega_n t)
\]
and a bath-driving delay kernel
\[
\lambda(t)=\sum_n \frac{c_n d_n}{m_n\omega_n}\sin(\omega_n t).
\]
The entire effect of bath driving is therefore captured by the effective force
\[
F_{\rm eff}(t)=d_0F(t)+\int_0^t ds\,\lambda(t-s)F(s).
\]
The bath-induced part is, in principle, non-Markovian, because it depends on the full history of the applied field [1807.01557].

For a harmonic central potential,
\[
V(q)=\frac{1}{2}M\omega_0^2 q^2,
\]
the long-time response to \(F(t)=F\cos(\omega t)\) is governed by the dynamic susceptibility
\[
\chi(\omega)=\frac{1}{M}\frac{d_0+\hat\lambda(-i\omega)}{\omega_0^2-i\omega\hat\gamma(-i\omega)-\omega^2}.
\]
The reported qualitative effect of bath driving is a low-frequency enhancement of the dispersive part, a maximum at zero frequency in \(\chi'(\omega)\), a shoulder-like feature in the absorptive part, and a slight shift of spectral weight toward lower frequencies [1807.01557].

The Rubin chain provides an explicit microscopic realization. A heavy central particle coupled to two semi-infinite harmonic chains can be mapped to a driven Caldeira-Leggett bath when the chain particles are also driven. In that mapping the Rubin spectral density is
\[
J(\omega)=m\omega\sqrt{\omega_D^2-\omega^2}, \qquad 0\le \omega\le \omega_D,
\]
and the damping kernel is
\[
\gamma(t)=\frac{m\omega_D}{M}\frac{J_1(\omega_D t)}{t},
\]
which decays oscillatory and algebraically, making the non-Markovian character explicit [1807.01557].

## 3. Driven spin-boson dynamics from coupled equilibrium

A second major realization of driven Caldeira-Leggett physics is the spin-boson model driven from an initially equilibrated coupled state. In the normal-mode representation the Hamiltonian is
\[
H_{\mathrm{tot}}(t)=\epsilon(t)\sigma_z + \Delta \sigma_x + \sum_\lambda \hbar\omega_\lambda b_\lambda^\dagger b_\lambda - \sigma_z \sum_\lambda c_\lambda (b_\lambda + b_\lambda^\dagger),
\]
with Ohmic spectral density
\[
I(\omega)=\alpha \omega \left[1+\left(\frac{\omega}{\omega_c}\right)^2\right]^{-2}.
\]
The central conceptual point is that the initial condition is not factorized; instead, the system and bath start in the thermal equilibrium state of the full interacting Hamiltonian [1712.06397].

The exact treatment is based on the Extended Stochastic Liouville Equation. Equilibrium preparation is performed in imaginary time through
\[
-\hbar \partial_\tau \overline{\rho}(\tau) = \left[ H_q(t_0)-\bar{\mu}(\tau)f(q) \right]\overline{\rho}(\tau),
\]
and the subsequent real-time evolution is
\[
i\hbar \partial_t \tilde{\rho}(t) = [H_q(t),\tilde{\rho}(t)]_- - \left\{ \eta(t)[f(q),\tilde{\rho}(t)]_- + \frac{\hbar}{2}\nu(t)[f(q),\tilde{\rho}(t)]_+ \right\}.
\]
The reduced density matrix is obtained by averaging over the stochastic noises, and the imaginary-time endpoint supplies the real-time initial condition [1712.06397].

Two regimes were studied. For constant bias, the exact partition-free dynamics remains stationary on average, confirming that the imaginary-time preparation reproduces the coupled thermal equilibrium. For a Landau-Zener sweep,
\[
\epsilon(t)=\kappa t,
\]
the transient dynamics and asymptotic state differ substantially from partitioned approximations. The reported effects include suppression of short-time oscillations by cross-time correlations between imaginary-time preparation and real-time evolution, movement of the asymptote toward the zero-temperature Landau-Zener limit as temperature is lowered, and suppression of coherent oscillations with increasing bath coupling \(\alpha\) [1712.06397].

This formulation shows that in driven Caldeira-Leggett dynamics the choice of initial condition is not a technical detail. Initial system-bath correlations can affect both transients and late-time observables.

## 4. Non-equilibrium Gaussian reservoirs as deterministic and stochastic drives

A more recent non-equilibrium version of the model prepares the bath modes in squeezed and displaced thermal states before system-bath coupling is switched on. The total Hamiltonian is
\[
\hat H(t)=\hat H_S+\sum_{\nu=1}^N \hat H_\nu+\sum_{\nu=1}^N \chi_\nu(t)\hat V_\nu,
\]
with
\[
\hat V_\nu=\sum_k V_{\nu,k}(\hat X)\otimes \hat X_{\nu,k}.
\]
The non-equilibrium ingredient lies in the reservoir preparation, not in an externally imposed field acting during the evolution [2405.00215].

Displacement and squeezing play different dynamical roles. For displacement,
\[
\hat D^\dagger(\alpha_{\nu,k}) \hat X_{\nu,k}\hat D(\alpha_{\nu,k}) =\hat X_{\nu,k}+L_{\nu,k},
\]
and, in the weak-coupling, strong-displacement limit, a single displaced mode yields approximately unitary reduced dynamics,
\[
\dot{\hat\rho}_S(t)= -\frac{i}{\hbar}\Big[\hat H_S+\chi_4 L_{4,1}V_{4,1}(\hat X)\cos(\omega_{4,1} t),\,\hat\rho_S(t)\Big].
\]
More generally,
\[
\hat H_{\rm eff}(t)=\hat H_S+f_4(t)V_4(\hat X),
\]
with
\[
f_4(t)=\chi_4\sum_k c_{4,k}A_{4,k}\cos(\omega_{4,k}t+\phi_k).
\]
The displaced reservoir therefore acts as a source of effective deterministic time dependence in the system Hamiltonian [2405.00215].

Squeezing modifies fluctuations instead of shifting the mean. In the weak-coupling, strong-squeezing limit, the system experiences stochastic driving,
\[
\dot{\hat\rho}_S^{(c)}(t) = -\frac{i}{\hbar}\big[\hat H_S+\xi \hat\Gamma(t),\,\hat\rho}_S^{(c)}(t)\big],
\]
with Gaussian \(\xi\). The same work shows that squeezing breaks the equilibrium fluctuation-dissipation relation, whereas displacement does not alter the fluctuation kernel in that way [2405.00215].

The thermodynamic interpretation is correspondingly modified. The work invested in preparing a non-equilibrium reservoir is
\[
\langle W_\nu\rangle= \operatorname{Tr}\!\left[\hat H_\nu\left(\hat\rho_\nu(0)-\frac{e^{-\beta_\nu \hat H_\nu}}{Z_\nu}\right)\right],
\]
and the total reservoir energy change, identified as heat, is
\[
\langle Q_\nu\rangle=\langle W_\nu\rangle+\langle \Delta E_\nu\rangle.
\]
The paper’s stated conclusion is that the apparent breakdown of equilibrium thermodynamic relations is reconciled with the second law once the energy spent to generate the squeezed or displaced initial state is included [2405.00215].

## 5. Spectroscopy, generalized Langevin descriptions, and the “noise-driven” interpretation

In condensed-phase spectroscopy, the expression “driven Caldeira-Leggett model” is used in a broader sense. One detailed study states that the term is not a separate formalism so much as the open-system, dissipative, noise-driven dynamics implied by the Caldeira-Leggett Hamiltonian and its generalized Langevin reduction. In that reduction the system momentum obeys
\[
\dot p(t)=F[x(t)]-\int_0^t \xi(t-\tau)p(\tau)\,d\tau+R(t),
\]
with memory kernel
\[
\xi(t)=\sum_i \frac{g_i^2}{m\omega_i^2}\cos(\omega_i t),
\]
and noise satisfying
\[
\langle R(0)R(t)\rangle = mkT\,\xi(t).
\]
The bath thus appears through deterministic forces, non-Markovian friction, and stochastic fluctuations consistent with thermal equilibrium [1601.04470].

For anharmonic system potentials, parameterization from molecular dynamics requires both the momentum autocorrelation function
\[
C_{pp}(t)=\langle p(t)p(0)\rangle
\]
and the momentum-force correlation function
\[
C_{pF}(t)=\langle F(t)p(0)\rangle,
\]
linked by the Volterra equation
\[
\dot C_{pp}(t)=C_{pF}(t)-\int_0^t \xi(t-\tau)C_{pp}(\tau)\,d\tau.
\]
The extracted spectral density is
\[
\hat{\xi}(\omega)=\frac{1+\hat C_{pF}(\omega)}{\hat C_{pp}(\omega)}-i\omega.
\]
This framework is technically attractive for vibrational spectroscopy, but its applicability is limited by the invertibility problem: in general there is no one-to-one mapping between a real many-body system and a unique Caldeira-Leggett or generalized Langevin representation [1601.04470] [1412.1688].

The sharpest statement of this limitation is that the mapping can be established self-consistently only when the system part is effectively harmonic. For a harmonic effective potential,
\[
V(x)=\frac{1}{2}m\tilde{\omega}^2 x^2,
\]
one has
\[
C_{pF}(t) = -\tilde{\omega}^2 \int_0^t C_{pp}(\tau)\, d\tau,
\]
so the two correlations are not independent. For anharmonic systems they remain genuinely independent, and compressing them into a single memory kernel loses information [1412.1688].

The same literature proposes empirical criteria for when the reduced driven/open-system description is likely to work: Gaussianity of the fluctuating force distribution, independence of the extracted spectral density on the choice of system potential, and, to a lesser extent, linearity of the system-bath coupling on the system side. A further stated implication is that an effectively harmonic mapping may reproduce linear spectra but cannot faithfully represent nonlinear spectroscopic dynamics, because higher-order response functions vanish for a purely harmonic system [1601.04470] [1412.1688].

## 6. Related nonequilibrium bath engineering: damped baths and thermal gradients

Some extensions modify the environment itself rather than applying an explicit periodic drive. One example replaces the ordinary bath oscillators by damped Caldirola-Kanai oscillators,
\[
\mathcal{L}_B=\sum_k e^{\mu_k t}\left(\frac12 m_k\dot R_k^2-\frac12 m_k\omega_{0_k}^2R_k^2\right),
\]
which leads, in the high-temperature Ohmic limit, to the master equation
\[
\frac{\partial \tilde{\rho}}{\partial t} = \frac{i\hbar}{2M}\left( \frac{\partial^2 \tilde{\rho}}{\partial x^2} - \frac{\partial^2 \tilde{\rho}}{\partial y^2} \right) -\frac{i}{\hbar}\big(V_F(x)-V_F(y)\big)\tilde{\rho} -\gamma (x-y)\left( \frac{\partial \tilde{\rho}}{\partial x} - \frac{\partial \tilde{\rho}}{\partial y} \right) -\frac{2M\gamma kT}{\hbar^2}(x-y)^2\tilde{\rho},
\]
with modified potential
\[
V_F(x)=V_R(x)+V_D(x), \qquad
V_D(x)=-\frac12 M\gamma\mu x^2.
\]
The new term is an inverted harmonic oscillator contribution induced by damping of the bath oscillators. In the double-well example reported in that work, the modification changes both short-time decoherence and long-time well-transfer probability, and a rescaling
\[
\gamma=\frac{\gamma_0}{a},\qquad \mu=\mu_0 a
\]
keeps \(\gamma\mu\) fixed and preserves the full potential while changing bath damping [2303.09516].

A different nonequilibrium extension addresses thermal gradients. In one construction, bath oscillators are driven by an external force proportional to the local gradient,
\[
m_k\ddot{q}_k = -m_k \omega_k^2 q_k + c_k x - \alpha_k T'(x),
\]
leading, in the Ohmic limit and for \(\alpha_k=\tilde{\alpha} c_k\), to
\[
M\ddot{x}=-V'(x) -\eta_{eff}[x] \dot{x}-\kappa T'(x) + F_L(t),
\]
with
\[
\eta_{eff}[x] = \eta(1-\tilde{\alpha}T''(x)).
\]
This model produces an explicit thermophoretic force \(-\kappa T'(x)\) and, for constant gradient, admits a clean Hamiltonian formulation. A second construction replaces the single bath by a continuum of local baths at temperatures \(T(X)\); in that case the gradient enters primarily through a space-dependent diffusion coefficient. Both models yield transport toward colder regions [2603.26980].

These variants are not “driven” in exactly the same sense as periodic forcing, but they show that within the Caldeira-Leggett logic environmental control parameters can enter through bath dynamics, bath statistics, or spatial bath structure.

## 7. Conceptual limits, non-Markovianity, and baseline equilibrium behavior

The driven Caldeira-Leggett model inherits the conceptual limitations of the undriven framework. The high-temperature Caldeira-Leggett master equation is historically central, but it is not exact and, as emphasized in the historical review, can violate complete positivity at low temperatures and short times. For a harmonic oscillator linearly coupled to a bosonic bath, the exact reduced dynamics is instead given by the Hu-Paz-Zhang master equation [2604.26524].

A stronger claim has been advanced for the microscopic position-position model: dissipation is a genuinely non-Markovian feature. In that analysis, the standard Caldeira-Leggett master equation
\[
\frac{d\rho}{dt} = -\frac{i}{\hbar}[H_S,\rho] -\frac{2m\gamma k_B T}{\hbar^2}[q,[q,\rho]] -\frac{i\gamma}{\hbar}[q,\{p,\rho\}]
\]
should be interpreted as an approximation to non-Markovian dynamics rather than as a true dissipative Markovian generator. In the strict Markov limit, a position-coupled bath yields only decoherence,
\[
\frac{d\rho}{dt} = -\frac{i}{\hbar}[H_S,\rho] -\frac{2m\gamma k_B T}{\hbar^2}[q,[q,\rho]],
\]
while the friction term originates from the imaginary part of the bath correlation function and therefore from bath memory [1701.05024].

For driven settings, this point is especially consequential. A plausible implication is that effective friction-like responses generated by bath driving, delayed forcing, or non-equilibrium reservoir engineering should be understood against a non-Markovian background rather than as purely local-in-time damping. The solvable equilibrium oscillator provides the reference case: an initially pure system coupled to a thermal bath approaches a thermal Gaussian state, its two-time correlation function tends to the thermal oscillator correlator, and late-time imaginary-time periodicity emerges only after thermalization [1212.3538].

Taken together, these results place the driven Caldeira-Leggett model in a precise niche. It is a flexible microscopic framework for driven dissipative quantum dynamics, response theory, decoherence, and engineered reservoir physics, but its predictive status depends on the bath model, the initial state, the presence or absence of memory, and the validity of the harmonic-bath reduction for the problem at hand.

Source: https://www.emergentmind.com/topics/driven-caldeira-leggett-model