---
title: 'Rainbow Chain: Exponential Hamiltonians'
url: https://www.emergentmind.com/topics/rainbow-chain
type: topic
---

# Rainbow Chain: Exponential Hamiltonians

Rainbow chain denotes a class of one-dimensional inhomogeneous local Hamiltonians whose nearest-neighbour couplings decay exponentially away from the centre. In the standard free-fermion realization, this deformation produces concentric singlet-type correlations across the central bond, a half-chain entanglement entropy that grows linearly with system size, and, in the weak-inhomogeneity regime, a continuum description in terms of a \(c=1\) conformal field theory on a curved background or, equivalently, a thermofield double at temperature \(T=h/(2\pi)\) [1503.02695, 1611.08559]. Recent work has used the same deformation to organize entanglement spectra into sharply resolved chiral conformal towers and to formulate a Wilsonian numerical renormalisation group for interacting exponentially decaying chains [2412.09685]. Closely related constructions include central-bond defect deformations, holographic descriptions based on AdS\(_3\) foliations, and tensor-network realizations of rainbow-entangled phases in Motzkin and Fredkin models [1912.10788, 1812.10023, 1811.11974].

## 1. Exponentially deformed lattice construction

The defining lattice feature is an exponential envelope of couplings away from the strongest central bond. In the free-fermion formulation emphasized by Szabó, one starts from a one-dimensional local Hamiltonian and deforms the nearest-neighbour couplings so that they fall off exponentially away from the centre,
\[
H=\sum_{n=1}^{L-1}J_n\,(c_n^\dagger c_{n+1}+c_{n+1}^\dagger c_n),\qquad
J_n=e^{-h(n-\tfrac12)}.
\]
Here \(h>0\) is the inhomogeneity, or rainbow, parameter. An equivalent notation is \(\lambda\equiv e^{-h}\), so that \(J_n=\lambda^{n-\tfrac12}\) [2412.09685].

An equivalent mirror-symmetric open-chain parametrization uses \(2L\) sites labelled by half-odd integers,
\[
H=-\frac12\,J_0\,(c_{+\tfrac12}^\dagger c_{-\tfrac12}+\mathrm{h.c.})
-\sum_{n=\tfrac12}^{L-\tfrac32}\frac{J_n}{2}\,
(c_n^\dagger c_{n+1}+c_{-n}^\dagger c_{-(n+1)}+\mathrm{h.c.}),
\]
with \(J_0=1\), \(J_n=\alpha^{2n}\), and \(\alpha=e^{-h/2}\) [1503.02695]. Both parametrizations implement the same physical idea: the entanglement cut is placed at the strongest bond, while all other bonds become exponentially weaker away from the centre.

In the strong-inhomogeneity limit, this geometry produces a “rainbow” of long bonds arching across the cut. In the language used in the earlier literature, the same state is the concentric singlet phase. This suggests that the term “rainbow chain” refers not only to a particular Hamiltonian profile but also to a characteristic entanglement pattern induced by exponential inhomogeneity.

## 2. Ground state and entanglement structure

For \(\alpha\to0\), or equivalently \(h\to\infty\), the ground state becomes a product of Bell-pair-like bonds connecting sites symmetrically placed about the centre,
\[
\bigl|\mathrm{rainbow}\bigr\rangle
=\prod_{k=1}^L\frac1{\sqrt2}
\bigl(c^\dagger_{-k+\tfrac12}\pm c^\dagger_{+k-\tfrac12}\bigr)\,|0\rangle.
\]
In this limit exactly \(L\) Bell pairs straddle the centre, and the half-chain entropy is
\[
S(A)=L\log2,
\]
which is a strict volume law in one dimension [1503.02695]. In the defect formulation of the model, the same strong-disorder picture is derived by a Dasgupta–Ma RG, yielding a product of concentric bonds for \(0<\gamma<1\) and area-law dimerized phases for \(\gamma<0\) or \(\gamma>1\) [1912.10788].

The weak-inhomogeneity regime \(0<h\ll1\) admits a continuum treatment. The inhomogeneous chain is mapped to a uniform CFT on the stretched coordinate
\[
\tilde x=\mathrm{sign}(x)\frac{e^{h|x|}-1}{h},
\qquad
\tilde L=\frac{e^{hL}-1}{h},
\]
and the half-chain von Neumann entropy becomes
\[
S(A)=\frac{c}{6}\ln(\tilde L)+\mathrm{const}
=\frac{c}{6}\ln\!\Bigl(\frac{e^{hL}-1}{h}\Bigr)+c'.
\]
For the rainbow chain \(c=1\). When \(hL\gg1\), this reduces to
\[
S(A)\simeq \frac{c}{6}\,hL,
\]
so the entropy again obeys a volume law, now derived from conformal methods rather than from the explicit Bell-pair limit [1503.02695, 1611.08559].

The same curved-CFT framework yields closed forms for more general bipartitions. For an interval anchored at the edge, \(A=[-L,x]\), the \(n\)-th Rényi entropy can be written as
\[
S_n(x)=\frac{n+1}{12n}\ln Y(x)+\frac{E_n}{2},
\]
with
\[
Y(x)=8\,e^{-h|x|}\,\frac{e^{hL}-1}{\pi h}\,
\cos\!\Bigl(\frac{\pi}{2}\frac{e^{h|x|}-1}{e^{hL}-1}\Bigr),
\]
and the von Neumann entropy follows from the limit \(n\to1\) [1611.08559]. These formulas show that the volume-law behavior is not restricted to the half-chain cut.

## 3. Curved-space, thermofield-double, and holographic descriptions

The continuum limit of the free rainbow chain is a massless Dirac fermion on a static curved background. After linearizing around the two Fermi points, the resulting Lorentzian metric is
\[
g_{00}=-e^{-2h|x|},\qquad g_{11}=+1,
\]
or, in Euclidean signature,
\[
ds^2=e^{2\sigma(x)}dt^2+dx^2,\qquad e^{\sigma(x)}=e^{-h|x|}.
\]
In the zweibein formalism the scalar curvature is
\[
R(x)=4h\,\delta(x)-h^2,
\]
so the background has constant negative curvature \(R=-h^2\) away from the origin [1611.08559].

This geometric reformulation leads directly to an effective temperature. The analytic continuation of the coordinate map extends the ground state to a strip of width \(\beta=2\pi/h\), and the state can be written in thermofield-double form,
\[
|\mathrm{rainbow}\rangle
=\sum_n e^{-\beta E_n/2}\,|n\rangle_L\otimes|n\rangle_R,
\qquad
T=\frac{1}{\beta}=\frac{h}{2\pi}.
\]
Equivalently, the reduced density matrix of one half-chain is approximately thermal,
\[
\rho_L\simeq e^{-\beta H_{CFT}},\qquad \beta=\frac{2\pi}{h}
\]
[1503.02695]. The appearance of a thermofield double and an effective temperature \(T=h/(2\pi)\) is described as strongly reminiscent of Hartle–Hawking and Unruh/Rindler constructions in the original discussion [1503.02695].

A holographic realization arrives at the same boundary geometry. By foliating AdS\(_3\) with AdS\(_2\) slices, MacCormack et al. obtained the induced boundary metric
\[
ds^2_{\rm bdy}=-e^{-2h|x|}dt^2+dx^2,
\]
with envelope \(f(x)=e^{-h|x|}\), in exact agreement with the rainbow-chain hopping amplitudes. The Ryu–Takayanagi prescription then reproduces the entanglement structure, including the half-chain result
\[
S_A^{\rm half}
=\frac{c}{3}\ln\!\Bigl[\frac{2}{h\epsilon}\sinh\!\Bigl(\frac{h\ell}{2}\Bigr)\Bigr],
\]
which crosses over to linear growth with slope \((c/6)h\) at \(\ell\gg1/h\), and the centre-cut result
\[
S_A^{\rm defect}=\frac{c}{3}\ln\!\frac{2}{h\epsilon},
\]
which is independent of subsystem size [1812.10023]. A plausible implication is that the same exponential deformation admits mutually consistent lattice, curved-CFT, and holographic descriptions.

## 4. Entanglement Hamiltonian and chiral conformal towers

A central development is the use of rainbow chains to read off detailed chiral CFT spectra from a single ground-state wavefunction. For a long uniform chain at a \((1+1)\)-dimensional conformal critical point, the reduced density matrix of the left half is
\[
\rho_A\equiv \mathrm{Tr}_B\,|\Psi\rangle\langle\Psi|
=e^{-2\pi H_{\mathrm{ent},A}},
\]
which defines an entanglement Hamiltonian. In free-fermion chains, \(H_{\mathrm{ent},A}\) is quadratic,
\[
H_{\mathrm{ent},A}=\sum_k \epsilon_k\,a_k^\dagger a_k,
\]
so the \(\epsilon_k\) form the single-particle entanglement spectrum [2412.09685].

Cardy, Läuchli, and Szabó’s key observation is that in continuum CFT the same \(\rho_A\) can be viewed as a thermal density matrix of a chiral half-chain of length \(\ell\) at inverse temperature \(\beta=2\pi\),
\[
\rho_A\propto e^{-2\pi H_{CFT}(\ell)}.
\]
On a segment with suitable boundary conditions, this \(H_{CFT}(\ell)\) has a purely chiral Virasoro spectrum, so the entanglement spectrum arranges into chiral conformal towers [2412.09685].

The conformal map underlying this construction is
\[
w=\log\tan z,\qquad z\in(-\pi/4,\pi/4)\subset \mathrm{Re}\,z,
\]
and, for a ring,
\[
w=\log\!\Bigl[i\,\frac{z+1}{z-1}\Bigr],\qquad |z|=1.
\]
Uniform discretization in \(w\)-space induces physical couplings in \(z\)-space proportional to \(|dw/dz|\), hence exponentially varying couplings and the rainbow Hamiltonian [2412.09685].

In the quadratic rainbow chain, the many-body entanglement energies
\[
\lambda_\alpha=\sum_{k\in\mathrm{occupied}}\epsilon_k
\]
form multiplets whose degeneracies match the chiral CFT degeneracies \(p(\Delta)\), and the zero-mode shift is the continuum primary dimension. Finite-size corrections scale empirically as \(O(1/w^2)\), in contrast to the \(O(1/\ln L)\) behavior of a uniform chain [2412.09685].

| Model | Entanglement-spectrum tower structure | Distinctive observation |
|---|---|---|
| Complex fermion, \(c=1\) | Degeneracies \(1,1,2,3,5,7,11,15,\dots\) in each charge sector \(q\); charge offsets \(\Delta_q^2/2\) | Multiplets resolved up to the 7th level |
| Majorana / critical Ising, \(c=\tfrac12\) | Even parity: identity tower; odd parity: \(\epsilon\)-tower; with SPT boundary conditions: \(\sigma\)-tower in both parities | Mirror eigenvalue on a periodic rainbow ring distinguishes even vs. odd Virasoro level |
| Three-state Potts, \(c=\tfrac45\) | Trivial \(S_3\) sector: \((0)\oplus(3)\); irrep \(E\): two copies of \((2/3)\) | Strong boundary fields induce \((1/8)\oplus(13/8)\) or \((1/40)\oplus(21/40)\) |

These results show that the rainbow deformation suppresses finite-size broadening enough to resolve eleven or more conformal levels in some examples, and that the entanglement spectrum can serve as a practical probe of chiral BCFT data even when the parent Hamiltonian is inaccessible [2412.09685].

## 5. Interacting rainbow chains and Wilsonian numerical renormalisation group

For non-free theories, the entanglement Hamiltonian cannot be reconstructed from a correlation matrix. Szabó therefore complements the rainbow-chain construction with a Wilsonian numerical renormalisation group tailored to exponentially decaying Hamiltonians. The algorithm folds the MPO of the rainbow Hamiltonian about the strongest central bond, turning sites \(i\) and \(L-i+1\) into a single supersite. The low-energy eigenstates of the first \(2i\) physical sites are represented by an MPS tensor \(M^{(i)}\) with two physical legs and one virtual leg. At step \(i\), one forms an effective Hamiltonian \(H_{\mathrm{eff}}^{(i+1)}\), diagonalizes it, retains only the \(\chi_{\mathrm{lowest}}\) eigenstates, stores them as \(M^{(i+1)}\), repeats until the flow converges to a fixed-point spectrum, and finally “unzips” the folded MPS by successive SVDs to recover the usual ordering and the entanglement spectrum [2412.09685].

Convergence is monitored by the stability of the low-energy levels as \(i\to L/2\) and by the decay of discarded weights. In the three-state Potts example,
\[
H=-\sum J_{i+\tfrac12}(X_iX_{i+1}^\dagger+\mathrm{h.c.})-\sum h_i(Z_i+Z_i^\dagger),
\]
the critical point \(J=h\) realizes the \(c=4/5\) minimal model. Using Wilson-NRG with \(\chi_{\mathrm{NRG}}=400\), the folded-chain fixed-point spectrum exhibits the towers \((0)\oplus(3)\) in the trivial \(S_3\) sector and two copies of \((2/3)\) in the two-dimensional irrep \(E\). After unzipping, the ground-state entanglement spectrum shows the same exact degeneracies. Strong boundary fields can be used to force the \((1/8)\oplus(13/8)\) or \((1/40)\oplus(21/40)\) towers, in agreement with BCFT predictions. The finite-size and fixed-\(\chi\) corrections scale as \(\sim1/w^2\) and \(\sim1/\chi^2\) [2412.09685].

This interacting extension is significant because it detaches the conformal-spectrum extraction problem from free-fermion integrability. A plausible implication is that rainbow deformation functions as a spectroscopy device for chiral CFT data across both quadratic and interacting critical chains.

## 6. Defects, symmetry variants, and related rainbow phases

A natural deformation introduces a tunable central defect \(\gamma\). In an open spin-\(\tfrac12\) chain of \(2L\) sites, the inhomogeneous XX Hamiltonian with bond-centred defect is
\[
H_L(h,\gamma)
=-\frac12\sum_{m=-L+1}^{L-1}
J_m\,c^\dagger_{m-\tfrac12}c_{m+\tfrac12}+\mathrm{h.c.},
\]
with
\[
J_m=
\begin{cases}
e^{-h|m|},&m\neq0,\\[4pt]
e^{-h\gamma},&m=0.
\end{cases}
\]
For \(\gamma=\tfrac12\) one recovers the usual rainbow couplings [1912.10788].

The strong-inhomogeneity phase diagram has three regimes. For \(0<\gamma<1\), the ground state remains a rainbow valence-bond state generated by the RG as a sequence of concentric bonds between sites \(\pm(i+\tfrac12)\). For \(\gamma<0\) or \(\gamma>1\), the ground state becomes a trivial dimer covering, with area-law entanglement bounded by \(\ln2\). At \(\gamma=0,1\), the three central couplings tie and the RG yields local multi-site modes and a valence-bond structure built from 4-site operators; cuts carry \(O(1)\) entanglement per block but do not grow with system size [1912.10788].

In the weak-inhomogeneity regime, the defect generates a scattering amplitude
\[
\tau=e^{-h(\gamma-\tfrac12)},
\]
and the half-chain entropy becomes
\[
S(L)\simeq \frac{c_{\mathrm{eff}}(\tau)}{6}\,hL+\cdots,
\]
with
\[
c_{\mathrm{eff}}(\tau)=\frac{6}{\pi^2}I[s(\tau)],
\qquad
s(\tau)=\frac{2\tau}{1+\tau^2}.
\]
The entanglement properties are invariant under \(\gamma\leftrightarrow1-\gamma\), and the self-dual point \(\gamma=\tfrac12\) is the defect-free rainbow chain. For \(\gamma\neq\tfrac12\), the single-particle entanglement spectrum is gapped, with the gap increasing with defect strength; at \(\gamma=\tfrac12\), one recovers the \(1/L\) spacing characteristic of CFT on a strip [1912.10788].

A site-centred version maps, after folding, to an inhomogeneous Su–Schrieffer–Heeger ladder. In the strong-coupling limit and for \(\gamma>1\), a trivial dimerized region forms around the defect, while an SSH SPT region persists toward the ends. The entanglement spectrum of a central block \(B_\ell\) shows three zero modes for \(\ell<\lfloor\gamma\rfloor\), then one topological zero mode for \(\ell\gg\gamma\) [1912.10788].

Related rainbow-entangled phases also arise outside the exponentially inhomogeneous free-fermion chain. In area-deformed Motzkin and Fredkin models, the deformation parameter \(t\) drives a transition from an area-law phase (\(t<1\)) to a rainbow or volume phase (\(t>1\)), where
\[
S_n(t>1)=n\ln c+O(1)
\]
for the left half-chain in the colorful case [1811.11974]. The associated “step-pyramid” tensor network exactly contracts to the ground state and interpolates between low-entanglement and rainbow-pairing regimes. A higher-dimensional free-fermion extension on a \(2L\times2L\) square lattice gives
\[
S(A)=A(\alpha)L^2+B(\alpha)L\ln L+C(\alpha)L,
\]
which is a genuine two-dimensional volume law for \(\alpha<1\) [1503.02695]. These related constructions indicate that “rainbow” has become a broader label for concentric, maximally nonlocal entanglement structures generated either by explicit inhomogeneous couplings or by weighted superpositions of extended paths.

Source: https://www.emergentmind.com/topics/rainbow-chain