---
title: Reaction-Coordinate Polaron-Transform Mapping
url: https://www.emergentmind.com/topics/reaction-coordinate-polaron-transform-mapping
type: topic
---

# Reaction-Coordinate Polaron-Transform Mapping

Reaction-coordinate polaron-transform mapping, usually abbreviated RCPT, is a framework for generating effective Hamiltonian models to treat nonequilibrium open quantum systems at strong coupling with their surroundings. It is based on two exact transformations of the Hamiltonian followed by a controlled truncation, and it ends with a new Hamiltonian with a weakened coupling to the environment. The resulting effective Hamiltonian mirrors the initial one, except that its parameters are dressed by the system-bath couplings [2211.05701].

## 1. Definition and conceptual setting

In the RCPT formulation, the starting point is a generic impurity-type Hamiltonian linearly coupled to a bosonic bath through bath displacements,
\[
\hat{H} = \hat{H}_s + \sum_k \nu_k \left(\hat{c}_k^{\dagger} + \frac{t_k}{\nu_k}\hat{S} \right) \left(\hat{c}_k + \frac{t_k}{\nu_k}\hat{S} \right).
\]
Here \(\hat H_s\) is the bare system Hamiltonian, \(\hat S\) is the system operator that couples to the bath, \(\hat c_k^\dagger,\hat c_k\) are bath operators, \(\nu_k\) are bath-mode frequencies, and \(t_k\) are coupling constants. The bath influence is encoded by the spectral density
\[
J(\omega)=\sum_k t_k^2 \delta(\omega-\nu_k).
\]
The RCPT construction combines two operations that, in the literature, often appear separately: a reaction-coordinate mapping and a polaron transformation [2211.05701].

Earlier reaction-coordinate work had already established the central physical idea behind the first step: one isolates a collective environmental mode—the reaction coordinate—and incorporates it into an enlarged “supersystem,” leaving a residual bath to be treated perturbatively [1805.08307]. In thermodynamic applications this is a redefinition of the system-environment boundary, not merely a frame change, and it is precisely this repartitioning that allows strong-coupling and non-Markovian effects to be retained explicitly inside the enlarged Hamiltonian [1602.01340].

The term should not be used indiscriminately. Several recent papers use a reaction-coordinate mapping without any polaron transformation. In the analysis of non-Markovian quantum heat statistics, for example, the method is explicitly “purely a reaction coordinate mapping”; no polaron unitary is introduced [2408.08829]. Likewise, machine-learning reconstruction of structured spectral densities uses an exact multi-reaction-coordinate remapping followed by a Born-Markov treatment of the residual bath only, and does not apply a polaron transformation [2501.07485]. Conversely, canonical small-polaron theory in the Wannier basis displaces all phonon modes directly through a canonical transformation but does not define a single collective phonon coordinate or a residual bath in the reaction-coordinate sense [2011.03620]. RCPT denotes the specific two-step construction in which both ingredients are present in sequence [2211.05701].

## 2. Exact reaction-coordinate stage

The first exact step is the extraction of a collective bosonic mode, the reaction coordinate (RC), into an enlarged system. The mapped Hamiltonian is
\[
\hat{H}_{RC} = \hat{H}_s +  \Omega \left(\hat{a}^{\dagger} + \frac{\lambda}{\Omega}\hat{S} \right) \left(\hat{a} + \frac{\lambda}{\Omega}\hat{S} \right) + \sum_k \omega_k \left(\hat{b}_k^{\dagger} + \frac{f_k}{\omega_k} (\hat{a}^{\dagger} + \hat{a}) \right) \left(\hat{b}_k + \frac{f_k}{\omega_k} (\hat{a}^{\dagger} + \hat{a}) \right),
\]
where \(\hat a^\dagger,\hat a\) are RC operators, \(\Omega\) is the RC frequency, \(\lambda\) is the system-RC coupling, and \(\hat b_k^\dagger,\hat b_k\) are residual-bath operators. The RC is defined by
\[
\lambda(\hat a^\dagger+\hat a)=\sum_k t_k(\hat c_k^\dagger+\hat c_k).
\]
Thus the dominant collective displacement coordinate that carried the original system-bath coupling becomes explicit [2211.05701].

The RC parameters are fixed by moments of the original spectral density,
\[
\lambda^2 = \frac{1}{\Omega} \int_0^\infty \omega J(\omega)\, d\omega,
\qquad
\Omega^2 = \frac{\int_0^\infty \omega^3 J(\omega)\, d\omega}{\int_0^\infty \omega J(\omega)\, d\omega},
\]
and the residual spectral density
\[
J_{RC}(\omega)=\sum_k f_k^2\delta(\omega-\omega_k)
\]
obeys the exact relation
\[
J_{RC}(\omega) = \frac{2\pi \lambda^2 J(\omega)} {\left[ P \int \frac{J(\omega')\, d\omega'}{\omega' - \omega} \right]^2 + \pi^2 J(\omega)^2}.
\]
A notable scaling property follows: if \(J(\omega)\to \alpha J(\omega)\), then \(\lambda \to \sqrt{\alpha}\,\lambda\), while \(J_{RC}(\omega)\) does not change with \(\alpha\). This is the formal reason strong original coupling is absorbed primarily into the system-RC sector rather than into the residual bath [2211.05701].

For the Brownian spectral density
\[
J(\omega)=\frac{4\gamma \Omega^2 \lambda^2 \omega}{(\omega^2-\Omega^2)^2 + (2\pi\gamma\Omega\omega)^2},
\]
the mapped residual bath is Ohmic,
\[
J_{RC}(\omega)=\gamma \omega e^{-|\omega|/\Lambda},
\]
exactly in the \(\Lambda\to\infty\) limit. This is a paradigmatic instance of the reaction-coordinate principle already emphasized in earlier RC thermodynamics: structured non-Markovian physics is moved into an enlarged system, while the remaining reservoir is rendered simpler [1805.08307].

## 3. Polaron transformation and effective Hamiltonian

After the exact RC mapping, RCPT applies a second exact transformation: a polaron transformation acting on the enlarged system,
\[
\hat U_P = e^{\frac{\lambda}{\Omega}(\hat a^\dagger-\hat a)\hat S}.
\]
It shifts the RC operator according to
\[
\hat U_P \hat a \hat U_P^\dagger = \hat a - \frac{\lambda}{\Omega}\hat S,
\]
and defines the transformed system Hamiltonian
\[
\hat{\tilde H}_s \equiv \hat U_P \hat H_s \hat U_P^\dagger.
\]
The full transformed Hamiltonian is
\[
\hat H_{RC-P}\equiv \hat U_P \hat H_{RC}\hat U_P^\dagger,
\]
namely
\[
\hat{H}_{RC-P} = \hat{\tilde{H}_{s} + \Omega \hat{a}^{\dagger}\hat{a} +\sum_k \omega_k \bigg\{\left[ \hat{b}_k^{\dagger} + \frac{f_k}{\omega_k} \left(\hat{a}^{\dagger} + \hat{a} -\frac{2\lambda}{\Omega}\hat{S}\right) \right] \left[ \hat{b}_k + \frac{f_k}{\omega_k} \left(\hat{a}^{\dagger} + \hat{a} -\frac{2\lambda}{\Omega}\hat{S}\right) \right]\bigg\}.
\]
At this stage the strong coupling has been moved into the transformed system sector and into dressed system operators [2211.05701].

The only approximation in the RCPT framework is a controlled truncation of the RC manifold. Keeping only the RC ground state \(|0\rangle\) yields
\[
\hat{H}^{eff}(\lambda) = \bra{0}\hat H_{RC-P}\ket{0}
= \bra{0}\hat{\tilde{H}_{s}\ket{0} + \sum_k \omega_k \left( \hat{b}_k^{\dagger} - \frac{2 \lambda f_k}{\Omega \omega_k}\hat{S} \right) \left( \hat{b}_k - \frac{2 \lambda f_k}{\Omega \omega_k}\hat{S} \right).
\]
The effective system Hamiltonian is therefore
\[
\hat H_s^{eff}(\lambda)=\bra{0}\hat{\tilde H}_s\ket{0},
\]
which can be written in closed form as
\[
\hat{H}_s^{eff}(\lambda) =
e^{-\frac{\lambda^2}{2\Omega^2}\hat S^2}
\left(
\sum_{n=0}^\infty \frac{\lambda^{2n}}{\Omega^{2n} n!}\, \hat S^n \hat H_s \hat S^n
\right)
e^{-\frac{\lambda^2}{2\Omega^2}\hat S^2}.
\]
The residual coupling is encoded by an effective spectral density
\[
J^{eff}(\omega)=\frac{4\lambda^2}{\Omega^2}J_{RC}(\omega).
\]
This final Hamiltonian mirrors the initial one in operator structure, but with dressed system parameters and a weakened remaining system-environment coupling [2211.05701].

The truncation is justified when the RC frequency is the dominant energy scale, particularly
\[
\Omega \gg \Delta,\lambda,T,
\]
with \(\Delta\) a characteristic system scale and \(T\) the bath temperature. The authors additionally assume \(\lambda/\Omega\ll1\) in practice so that \(\Omega\) remains the dominant scale [2211.05701].

## 4. Relation to reaction-coordinate and polaron methods

Reaction-coordinate mapping and the polaron transformation reorganize strong system-bath coupling in different ways. In the RC picture, a specific collective environmental mode is promoted into the system Hamiltonian, and the residual bath is coupled to that mode rather than directly to the original subsystem [1805.08307]. In the standard polaron picture, bath modes are displaced conditional on the system state, so that strong dressing is absorbed into transformed operators and renormalized parameters without changing the same system-bath partition in the same way [2011.03620]. RCPT combines these two ideas sequentially: first the boundary is redrawn, then the enlarged system is displaced [2211.05701].

The thermodynamic distinction is explicit in RC-only work. In strong-coupling heat statistics for the spin-boson model, the reaction coordinate is absorbed into the enlarged system and energy stored in the RC is therefore not automatically classified as heat; this is one reason the residual-environment definition of heat behaves more like a dissipative thermal reservoir in the highly non-Markovian regime considered [2408.08829]. This suggests that RCPT should be interpreted as more than a convenient frame change: the polaron step acts after a genuine repartitioning of the Hamiltonian.

Several adjacent constructions stop short of full RCPT. One thermodynamic study argued that the RC framework allows one to justify master equations derived in a polaron transformed reference frame for a single electron transistor coupled to vibrations; there the polaron master equation appears as the fast-RC limit of the RC description rather than as a separate exact mapping [1602.01340]. A transport study beyond second order applied a small-polaron transformation after RC mapping as an interpretive tool, identifying an effective term
\[
-i\frac{ \lambda_R\lambda_L}{\Omega_L} (\hat{a}_R^{\dagger} + \hat{a}_R)\hat{\sigma}_y (\hat{a}_L^{\dagger} - \hat{a}_L)
\]
responsible for an inter-bath current that scales as \(j_q \propto \lambda^4\) [2203.06165]. By contrast, multi-RC machine-learning work uses exact RC remapping followed by Born-Markov treatment of the residual baths, with no polaron step [2501.07485]. The named RCPT framework is therefore a specific, later synthesis rather than a generic label for any study that mentions both RC and polaron ideas [2211.05701].

## 5. Dressed physics and representative applications

The most distinctive output of RCPT is an effective Hamiltonian in which strong-coupling physics is encoded analytically in dressed parameters. In the generalized spin-boson model,
\[
\hat{H} = \Delta \hat{\sigma}_z + \sum_k \nu_k \left( \hat{c}_k^{\dagger} + \frac{t_k}{\nu_k}\hat{\sigma}_\theta \right) \left( \hat{c}_k + \frac{t_k}{\nu_k}\hat{\sigma}_\theta \right),
\qquad
\hat{\sigma}_\theta = \cos\theta\,\hat{\sigma}_z + \sin\theta\,\hat{\sigma}_x,
\]
RCPT yields
\[
\hat{H}_s^{eff}(\lambda) = \frac{\Delta}{2}\left[(1+e^{-2\lambda^2/\Omega^2}) + (1-e^{-2\lambda^2/\Omega^2})\cos(2\theta)\right]\hat{\sigma}_z + \frac{\Delta}{2}(1-e^{-2\lambda^2/\Omega^2})\sin(2\theta)\hat{\sigma}_x.
\]
For \(\theta=0\), there is no renormalization because \([\hat H_s,\hat S]=0\); for \(\theta=\pi/2\),
\[
\hat H_s^{eff}=\Delta e^{-2\lambda^2/\Omega^2}\hat\sigma_z.
\]
RCPT therefore exposes both exponential quenching of the level splitting and the appearance of new bath-induced tunneling terms [2211.05701].

For nonequilibrium heat transport in the spin-boson model with one strongly coupled RC extracted from each bath, the effective qubit splitting becomes
\[
\hat H_s^{eff} = \Delta e^{-\sum_{\alpha=L,R}2\lambda_\alpha^2/\Omega_\alpha^2}\hat\sigma_z.
\]
The paper uses this structure to explain turnover of the heat current: stronger coupling initially promotes transitions, but eventually suppresses the effective energy gap that is being transported [2211.05701].

In the autonomous three-level absorption refrigerator, RCPT dresses the cold transition and modifies the cooling condition through the dressed levels
\[
\epsilon_1(\lambda_c)=\frac{\Delta}{2}\left(1-e^{-2\lambda_c^2/\Omega_c^2}\right),
\qquad
\epsilon_2(\lambda_c)=\frac{\Delta}{2}\left(1+e^{-2\lambda_c^2/\Omega_c^2}\right),
\qquad
\epsilon_3(\lambda_c)=1,
\]
so that the cooling window is shifted by strong coupling [2211.05701].

In phonon-assisted charge transport through a double quantum dot, the effective onsite energies become
\[
\epsilon_L(\lambda)= \left[ \epsilon_L\cosh\!\left(\frac{\lambda^2}{\Omega^2}\right) +\epsilon_R\sinh\!\left(\frac{\lambda^2}{\Omega^2}\right) \right]e^{-\lambda^2/\Omega^2},
\]
\[
\epsilon_R(\lambda)= \left[ \epsilon_R\cosh\!\left(\frac{\lambda^2}{\Omega^2}\right) +\epsilon_L\sinh\!\left(\frac{\lambda^2}{\Omega^2}\right) \right]e^{-\lambda^2/\Omega^2},
\]
while the lead couplings are renormalized as
\[
h_{k,L}(\lambda)=h_{k,L}e^{-\lambda^2/(2\Omega^2)},
\qquad
h_{k,R}(\lambda)=h_{k,R}e^{-\lambda^2/(2\Omega^2)}.
\]
The paper interprets the resulting transport turnover as a competition between level renormalization toward degeneracy and exponential suppression of hybridization [2211.05701].

For dissipative spin chains with one RC extracted per site, the effective Hamiltonian suppresses local splittings and the \(J_y\) and \(J_z\) couplings while \(J_x\) survives. This places strong dissipation directly into an effective many-body Hamiltonian and provides analytic access to how dissipation reshapes the chain [2211.05701].

## 6. Validity, limitations, and thermodynamic interpretation

The approximation in RCPT is not the RC mapping or the polaron transformation themselves, both of which are exact, but the projection onto the RC ground state. Because the RC is frozen into \(|0\rangle\), transient non-Markovian dynamical features associated with explicit RC excitations can be lost; oscillations and memory effects tied to RC-system information exchange are therefore not fully captured. The method is correspondingly strongest for steady-state and low-energy properties, and less complete for full transient dynamics [2211.05701].

The framework also inherits the usual delicacy of residual-bath treatments. A benchmark study of RC mapping without the polaron step showed that the reduced stationary state of the original problem can remain accurate even when residual dissipation on the augmented system is very strong, but that the same weak-coupling master equation can grossly overestimate the stationary heat current across the wire [1905.13673]. This suggests that RC-based embeddings can be much more reliable for reduced states than for transport observables whenever the residual-bath approximation is stressed.

For thermodynamics, RC literature has consistently emphasized that strong coupling changes bookkeeping. In the supersystem description, internal energy includes the original system, the RC, and their interaction, and the stationary reduced state is linked to the Hamiltonian of mean force rather than to a Gibbs state of the bare subsystem alone [1602.01340]. In heat-counting problems, whether RC energy is counted as environmental energy or retained on the system side can qualitatively change the interpretation of “heat” [2408.08829]. A plausible implication is that RCPT should be interpreted thermodynamically at the mapped supersystem-residual-bath boundary, where the weakened residual coupling is introduced and where standard weak-coupling reasoning is intended to apply.

A final structural limitation concerns multi-bath and multi-RC settings. When several polaron transformations do not commute, the order of transformations may matter, and the construction is not unique. The approximation can be systematically improved by retaining higher RC manifolds, but the principal simplicity of RCPT comes precisely from avoiding that enlarged Hilbert-space growth [2211.05701].

Source: https://www.emergentmind.com/topics/reaction-coordinate-polaron-transform-mapping