Reaction-Coordinate Mapping

Reaction-coordinate mapping is a way to redraw the boundary between a quantum system and its environment. A collective environmental coordinate that is strongly coupled to the system, slowly relaxing, or spectrally structured is promoted into the system. The remaining residual bath is then treated as weaker, broader, and more nearly Markovian.

The guiding picture is:

original system + structured bath
    = enlarged system + residual bath

The original reduced system may show memory. The enlarged system can often be described by a time-local master equation. The memory has not disappeared; it has been carried by the reaction coordinate, which is now explicit.

For a closely related pole-based construction, see Pseudomode Methods. Pseudomodes are often tied to Lorentzian response poles or leaky resonances. Reaction-coordinate mappings emphasize collective coordinates of a bosonic environment and are especially useful for strong-coupling and thermodynamic questions.

Starting Hamiltonian

A common bosonic system–bath model is

$$H = H_S + A\otimes B + H_B,$$

with

$$B = \sum_k g_k \left( b_k+b_k^\dagger \right), \qquad H_B = \sum_k \hbar\omega_k b_k^\dagger b_k.$$

The bath is often summarized by a spectral density,

$$J(\omega) = \sum_k g_k^2\delta(\omega-\omega_k), \qquad \omega>0.$$

If $J(\omega)$ is broad and smooth near the system transition frequencies, a weak-coupling Markovian master equation may be adequate. If $J(\omega)$ has a narrow peak, a strong low-frequency component, or an underdamped vibrational structure, the system may retain correlations with a collective environmental degree of freedom. The reaction-coordinate mapping makes that degree of freedom explicit.

Mapped Hamiltonian

After the mapping, one writes an enlarged Hamiltonian of the form

$$H = H_{S+\mathrm{RC}} + H_{\mathrm{res}} + H_{\mathrm{RC-res}},$$

where

$$H_{S+\mathrm{RC}} = H_S + \hbar\Omega a^\dagger a + \lambda A \left( a+a^\dagger \right) + H_{\mathrm{ct}}.$$

Here $a$ is the reaction-coordinate mode, $\Omega$ is its frequency, $\lambda$ is its coupling to the system operator $A$, and $H_{\mathrm{ct}}$ denotes any counterterm or renormalization term required by the convention.

The residual bath has modes $c_q$:

$$H_{\mathrm{res}} = \sum_q \hbar\nu_q c_q^\dagger c_q,$$

and the reaction coordinate couples to it through

$$H_{\mathrm{RC-res}} = \left( a+a^\dagger \right) \sum_q h_q \left( c_q+c_q^\dagger \right).$$

The mapping is chosen so that tracing out the reaction coordinate and residual bath reproduces the original bath influence on $S$. The residual spectral density is not arbitrary; it is determined by the original $J(\omega)$ and the mapping convention.

Why This Helps

The original system may be strongly coupled to the structured part of the environment. Perturbation theory in $A\otimes B$ can then fail.

After the mapping, the strong coupling $\lambda A(a+a^\dagger)$ is included in $H_{S+\mathrm{RC}}$ and treated nonperturbatively. The remaining coupling to the residual bath may be weak enough for controlled Markov, secular, or coarse-grained approximations.

Thus the approximation is not:

system weakly coupled to original bath

but rather:

enlarged system weakly coupled to residual bath

This distinction is the whole point. The residual bath should be simpler than the original structured bath, but that still has to be checked.

Brownian-Oscillator Example

A standard structured spectral density has an underdamped Brownian-oscillator form,

$$J_{\mathrm{BO}}(\omega) = \frac{ 2\lambda^2\gamma\Omega^2\omega } { (\Omega^2-\omega^2)^2+\gamma^2\omega^2 }, \qquad \omega>0.$$

The parameters suggest the mapping directly:

• $\Omega$ is the reaction-coordinate frequency;
• $\lambda$ sets the system-reaction-coordinate coupling strength;
• $\gamma$ sets the damping of the reaction coordinate by the residual bath.

In the mapped model, the system couples to a single oscillator of frequency $\Omega$, and that oscillator is damped by a smoother residual bath. If $\gamma$ is small, the reaction coordinate stores excitation for a long time and memory is prominent. If $\gamma$ is large compared with the system-reaction-coordinate exchange scale, the reaction coordinate can often be eliminated again, leading back toward a Markovian effective model.

The exact prefactors depend on spectral-density conventions. The structural point is robust: a peaked spectral density can be represented as a distinguished oscillator plus a residual damping mechanism.

Master Equation on the Enlarged System

Once the reaction coordinate is included in the system, the density matrix $\varrho$ for $S+\mathrm{RC}$ obeys an equation schematically like

$$\frac{d\varrho}{dt} = - \frac{i}{\hbar} [H_{S+\mathrm{RC}},\varrho] + \mathcal L_{\mathrm{res}}(\varrho).$$

The residual dissipator $\mathcal L_{\mathrm{res}}$ should be derived from the coupling between the reaction coordinate and the residual bath. In a careful weak-coupling treatment, one diagonalizes $H_{S+\mathrm{RC}}$ and constructs jump operators at the Bohr frequencies of the enlarged Hamiltonian.

This is important. A local dissipator such as

$$\kappa\mathcal D[a]\varrho$$

may be a useful approximation in some regimes, but it is not automatically thermodynamically consistent when $a$ is strongly hybridized with $S$. The residual bath sees the dressed eigenstates of the enlarged system, not the bare reaction coordinate in isolation.

Reduced Memory

The original system state is obtained by tracing out the reaction coordinate:

$$\rho_S(t) = \operatorname{Tr}_{\mathrm{RC}}\varrho(t).$$

Even if $\varrho(t)$ follows a Markovian master equation, $\rho_S(t)$ can show memory because information can move from $S$ into the reaction coordinate and later return. The reaction coordinate acts as an explicit memory register.

This is the same structural lesson as Pseudomode Methods: Markovianity depends on where the system boundary is drawn. The enlarged description may be CP-divisible while the smaller reduced description is not.

Strong-Coupling Equilibrium

Reaction-coordinate mappings are especially useful for equilibrium and thermodynamic questions at strong coupling.

If the enlarged system thermalizes with a residual bath at inverse temperature $\beta$, its equilibrium state may be close to

$$\varrho_{\mathrm{eq}} = \frac{ e^{-\beta H_{S+\mathrm{RC}}} } { \operatorname{Tr}_{S,\mathrm{RC}} e^{-\beta H_{S+\mathrm{RC}}} }.$$

The reduced system equilibrium is then

$$\rho_S^{\mathrm{eq}} = \operatorname{Tr}_{\mathrm{RC}} \varrho_{\mathrm{eq}}.$$

This is generally not the bare Gibbs state

$$\frac{e^{-\beta H_S}}{\operatorname{Tr}_S e^{-\beta H_S}}.$$

It is instead related to a Hamiltonian of mean force. This distinction explains why demanding a bare-Gibbs steady state can be wrong outside the weak-coupling limit. See Thermal Master Equations for the weak-coupling thermal case.

Relation to Chain Mappings

Reaction-coordinate mapping can be viewed as the first step of a star-to-chain transformation. The original bath modes form a star around the system. Orthogonal-polynomial or Lanczos methods transform the bath into a chain:

system - site 0 - site 1 - site 2 - ...

The first site is the reaction coordinate. Keeping only that site in the enlarged system and treating the rest as a residual bath gives the basic reaction-coordinate picture. Keeping more sites gives a longer explicit memory chain and can improve accuracy for strongly structured environments.

This perspective is common in tensor-network simulations and numerically exact open-system methods. The tradeoff is transparent: more explicit bath sites mean less memory left in the residual bath but a larger system Hilbert space.

Practical Workflow

A reaction-coordinate calculation usually follows this sequence.

1. Specify the original $H_S$, coupling operator $A$, bath state, and spectral density $J(\omega)$.
2. Choose a reaction coordinate that captures the dominant structured part of the bath.
3. Derive or fit $\Omega$, $\lambda$, counterterms, and the residual spectral density.
4. Build $H_{S+\mathrm{RC}}$ and choose a reaction-coordinate Hilbert-space cutoff.
5. Derive $\mathcal L_{\mathrm{res}}$ for the residual bath, preferably in the dressed eigenbasis when coupling is strong.
6. Evolve $\varrho(t)$ or compute its steady state.
7. Trace out the reaction coordinate to obtain $\rho_S(t)$ and system observables.
8. Test convergence against cutoff, mapping depth, residual-bath approximation, and known limiting cases.

The method is constructive, but not automatic. The residual bath assumptions and numerical truncations are part of the model.

When to Use It

Reaction-coordinate mapping is useful when:

• one environmental feature is strongly coupled to the system;
• the bath spectrum has a pronounced peak or underdamped mode;
• weak-coupling rates built from the original $J(\omega)$ fail;
• equilibrium at strong coupling matters;
• the environment is better viewed as a slow coordinate plus a faster background;
• a finite explicit memory mode makes the model easier to simulate or interpret.

It is less natural when the environment has broad algebraic memory, many comparable peaks, a strongly non-Gaussian spin structure, or no clear coordinate to isolate. In those cases, hierarchical methods, tensor-network chain simulations, collision models, or direct microscopic modeling may be more appropriate.

For Gaussian baths whose correlation functions can be expanded into exponentials, Hierarchical Equations of Motion provide a systematic auxiliary-density-operator alternative to choosing a small number of explicit coordinates.

Common Mistakes

Calling the mapping an approximation by itself

The exact transformation of a harmonic bath can be formal. The approximations usually enter later: truncating the reaction-coordinate Hilbert space, keeping only one chain site, fitting a residual spectrum, or applying a Markovian residual-bath master equation.

Using local damping by habit

A dissipator $\kappa\mathcal D[a]$ is not automatically correct for a strongly coupled reaction coordinate. If $S$ and the reaction coordinate are hybridized, the residual bath induces transitions between dressed eigenstates.

Comparing to the wrong Gibbs state

At strong coupling, the reduced equilibrium state need not be the Gibbs state of $H_S$. Compare to the reduced enlarged equilibrium or Hamiltonian of mean force, not automatically to $e^{-\beta H_S}$.

Losing the counterterm

Linear coordinate couplings can shift potentials and frequencies. Depending on convention, a counterterm prevents double-counting or unintended renormalization. Dropping it can change the physical model.

Truncating the oscillator too aggressively

The reaction coordinate is a bosonic mode. Strong coupling, finite temperature, and transient excitation may require more Fock states than intuition from weak coupling suggests.

Forgetting the initial state

If the original bath is thermal, the mapped reaction coordinate may initially be correlated with the residual bath in the exact thermal state. Product initial states are convenient but should be justified for the intended protocol.

Exercises

Reduced Equilibrium Is Not Bare Gibbs

Assume the enlarged equilibrium state is

$$\varrho_{\mathrm{eq}} = \frac{ e^{-\beta H_{S+\mathrm{RC}}} } { \operatorname{Tr}_{S,\mathrm{RC}} e^{-\beta H_{S+\mathrm{RC}}} },$$

with

$$H_{S+\mathrm{RC}} = H_S+\hbar\Omega a^\dagger a+\lambda A(a+a^\dagger).$$

Why is $\rho_S^{\mathrm{eq}}=\operatorname{Tr}_{\mathrm{RC}}\varrho_{\mathrm{eq}}$ generally not equal to the bare Gibbs state of $H_S$?

Solution

If $\lambda=0$, the exponential factorizes:

$$e^{-\beta H_{S+\mathrm{RC}}} = e^{-\beta H_S} e^{-\beta\hbar\Omega a^\dagger a},$$

so tracing out the reaction coordinate leaves the bare Gibbs state of $H_S$.

For $\lambda\ne0$, the interaction term is inside the exponential. In general,

$$[H_S+\hbar\Omega a^\dagger a,\lambda A(a+a^\dagger)]\ne0,$$

and the exponential does not factorize. The reaction coordinate remains correlated with the system in equilibrium, so the reduced state contains coupling-dependent corrections. These corrections are summarized by the Hamiltonian of mean force.

Memory from an Explicit Coordinate

Explain how a Markovian master equation for $S+\mathrm{RC}$ can lead to non-Markovian reduced dynamics for $S$.

Solution

The enlarged state $\varrho(t)$ contains correlations between $S$ and the reaction coordinate. At an intermediate time $s$, the reduced state

$$\rho_S(s)=\operatorname{Tr}_{\mathrm{RC}}\varrho(s)$$

does not determine the full $\varrho(s)$. The future of $S$ can depend on excitation or phase information stored in the reaction coordinate. Therefore there may be no completely positive intermediate map acting only on $\rho_S(s)$, even though $S+\mathrm{RC}$ evolves Markovianly.

Local Versus Dressed Dissipation

Why can a local dissipator $\kappa\mathcal D[a]$ be unreliable when $\lambda A(a+a^\dagger)$ is large?

Solution

The residual bath couples to the reaction-coordinate coordinate, but the energy eigenstates of the enlarged system are eigenstates of

$$H_{S+\mathrm{RC}} = H_S+\hbar\Omega a^\dagger a+\lambda A(a+a^\dagger).$$

When $\lambda$ is large, these eigenstates are dressed mixtures of system and reaction-coordinate states. Thermal transitions induced by the residual bath should occur at Bohr frequencies of $H_{S+\mathrm{RC}}$. A bare local dissipator may impose rates at the wrong frequencies and can give an incorrect steady state or heat current.

One More Chain Site

In a chain mapping, what is the expected effect of keeping two explicit bath sites instead of only the first reaction coordinate?

Solution

Keeping a second site moves more environmental memory into the enlarged system. The residual bath seen by the second site is often smoother and shorter-correlated than the residual bath seen by the first site. The price is a larger Hilbert space and more numerical cost. If results change significantly when the second site is added, the one-coordinate mapping plus residual-bath approximation was not yet converged.

Cross-Links

• System–Bath Hamiltonians
• Memory Kernels
• Non-Markovian Dynamics
• Pseudomode Methods
• Hierarchical Equations of Motion
• Collision Models
• Strong Coupling
• Noise Spectra
• Caldeira–Leggett Model
• Thermal Master Equations
• Detailed Balance
• Approximation Checklist

References

• A. Garg, J. N. Onuchic, and V. Ambegaokar, “Effect of friction on electron transfer in biomolecules,” Journal of Chemical Physics 83, 4491-4503 (1985).
• A. W. Chin, A. Rivas, S. F. Huelga, and M. B. Plenio, “Exact mapping between system-reservoir quantum models and semi-infinite discrete chains using orthogonal polynomials,” Journal of Mathematical Physics 51, 092109 (2010).
• R. Martinazzo, B. Vacchini, K. H. Hughes, and I. Burghardt, “Communication: Universal Markovian reduction of Brownian particle dynamics,” Journal of Chemical Physics 134, 011101 (2011).
• J. Iles-Smith, N. Lambert, and A. Nazir, “Environmental dynamics, correlations and the emergence of noncanonical equilibrium states in open quantum systems,” Physical Review A 90, 032114 (2014).
• P. Strasberg, G. Schaller, N. Lambert, and T. Brandes, “Nonequilibrium thermodynamics in the strong coupling and non-Markovian regime based on a reaction coordinate mapping,” New Journal of Physics 18, 073007 (2016).
• A. Nazir and G. Schaller, “The reaction coordinate mapping in quantum thermodynamics,” in Thermodynamics in the Quantum Regime, Springer (2018).
• I. de Vega and D. Alonso, “Dynamics of non-Markovian open quantum systems,” Reviews of Modern Physics 89, 015001 (2017).