Reaction-Coordinate Polaron-Transform Mapping

• Reaction-coordinate polaron-transform mapping is a framework that employs two exact transformations and a controlled truncation to derive effective Hamiltonians for strongly coupled open quantum systems.
• It reorganizes the system-bath interaction by isolating a collective environmental mode and then applying a polaron transformation to dress system parameters.
• This method provides analytic insights into phenomena such as energy renormalization, tunneling suppression, and transport turnover in nonequilibrium setups.

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 (Anto-Sztrikacs et al., 2022).

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 (Anto-Sztrikacs et al., 2022).

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 (Nazir et al., 2018). 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 (Strasberg et al., 2016).

The term should not be used indiscriminately. Several 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 (Shubrook et al., 2024). 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 (Barr et al., 13 Jan 2025). 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 (Lee et al., 2020). RCPT denotes the specific two-step construction in which both ingredients are present in sequence (Anto-Sztrikacs et al., 2022).

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, $\hat H_s$0 is the system-RC coupling, and $\hat H_s$1 are residual-bath operators. The RC is defined by

$\hat H_s$2

Thus the dominant collective displacement coordinate that carried the original system-bath coupling becomes explicit (Anto-Sztrikacs et al., 2022).

The RC parameters are fixed by moments of the original spectral density,

$\hat H_s$3

and the residual spectral density

$\hat H_s$4

obeys the exact relation

$\hat H_s$5

A notable scaling property follows: if $\hat H_s$6, then $\hat H_s$7, while $\hat H_s$8 does not change with $\hat H_s$9. This is the formal reason strong original coupling is absorbed primarily into the system-RC sector rather than into the residual bath (Anto-Sztrikacs et al., 2022).

For the Brownian spectral density

$\hat S$0

the mapped residual bath is Ohmic,

$\hat S$1

exactly in the $\hat S$2 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 (Nazir et al., 2018).

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 S$3

It shifts the RC operator according to

$\hat S$4

and defines the transformed system Hamiltonian

$\hat S$5

The full transformed Hamiltonian is

$\hat S$6

namely

$\hat S$7

At this stage the strong coupling has been moved into the transformed system sector and into dressed system operators (Anto-Sztrikacs et al., 2022).

The only approximation in the RCPT framework is a controlled truncation of the RC manifold. Keeping only the RC ground state $\hat S$8 yields

$\hat S$9

The effective system Hamiltonian is therefore

$\hat c_k^\dagger,\hat c_k$0

which can be written in closed form as

$\hat c_k^\dagger,\hat c_k$1

The residual coupling is encoded by an effective spectral density

$\hat c_k^\dagger,\hat c_k$2

This final Hamiltonian mirrors the initial one in operator structure, but with dressed system parameters and a weakened remaining system-environment coupling (Anto-Sztrikacs et al., 2022).

The truncation is justified when the RC frequency is the dominant energy scale, particularly

$\hat c_k^\dagger,\hat c_k$3

with $\hat c_k^\dagger,\hat c_k$4 a characteristic system scale and $\hat c_k^\dagger,\hat c_k$5 the bath temperature. The authors additionally assume $\hat c_k^\dagger,\hat c_k$6 in practice so that $\hat c_k^\dagger,\hat c_k$7 remains the dominant scale (Anto-Sztrikacs et al., 2022).

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 (Nazir et al., 2018). 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 (Lee et al., 2020). RCPT combines these two ideas sequentially: first the boundary is redrawn, then the enlarged system is displaced (Anto-Sztrikacs et al., 2022).

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 (Shubrook et al., 2024). 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 (Strasberg et al., 2016). A transport study beyond second order applied a small-polaron transformation after RC mapping as an interpretive tool, identifying an effective term

$\hat c_k^\dagger,\hat c_k$8

responsible for an inter-bath current that scales as $\hat c_k^\dagger,\hat c_k$9 (Anto-Sztrikacs et al., 2022). By contrast, multi-RC machine-learning work uses exact RC remapping followed by Born-Markov treatment of the residual baths, with no polaron step (Barr et al., 13 Jan 2025). 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 (Anto-Sztrikacs et al., 2022).

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,

$\nu_k$0

RCPT yields

$\nu_k$1

For $\nu_k$2, there is no renormalization because $\nu_k$3; for $\nu_k$4,

$\nu_k$5

RCPT therefore exposes both exponential quenching of the level splitting and the appearance of new bath-induced tunneling terms (Anto-Sztrikacs et al., 2022).

For nonequilibrium heat transport in the spin-boson model with one strongly coupled RC extracted from each bath, the effective qubit splitting becomes

$\nu_k$6

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 (Anto-Sztrikacs et al., 2022).

In the autonomous three-level absorption refrigerator, RCPT dresses the cold transition and modifies the cooling condition through the dressed levels

$\nu_k$7

so that the cooling window is shifted by strong coupling (Anto-Sztrikacs et al., 2022).

In phonon-assisted charge transport through a double quantum dot, the effective onsite energies become

$\nu_k$8

$\nu_k$9

while the lead couplings are renormalized as

$t_k$0

The paper interprets the resulting transport turnover as a competition between level renormalization toward degeneracy and exponential suppression of hybridization (Anto-Sztrikacs et al., 2022).

For dissipative spin chains with one RC extracted per site, the effective Hamiltonian suppresses local splittings and the $t_k$1 and $t_k$2 couplings while $t_k$3 survives. This places strong dissipation directly into an effective many-body Hamiltonian and provides analytic access to how dissipation reshapes the chain (Anto-Sztrikacs et al., 2022).

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 $t_k$4, 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 (Anto-Sztrikacs et al., 2022).

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 (Correa et al., 2019). 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 (Strasberg et al., 2016). 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” (Shubrook et al., 2024). 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 (Anto-Sztrikacs et al., 2022).