---
title: Inhomogeneous Open XX Spin Chain
url: https://www.emergentmind.com/topics/inhomogeneous-open-xx-spin-chain
type: topic
---

# Inhomogeneous Open XX Spin Chain

An inhomogeneous open XX spin chain is a nearest-neighbor spin-\(\tfrac12\) chain on a finite interval, with site-dependent couplings and/or local magnetic fields and with no periodic closing bond. Its central mathematical feature is that, despite spatial inhomogeneity and open boundaries, it remains reducible by Jordan–Wigner fermionization to a quadratic free-fermion problem governed by a tridiagonal Jacobi matrix. This reduction underlies a wide range of exact results: explicit diagonalization in the one-excitation sector, orthogonal-polynomial constructions, finite-size nonlocal symmetries, perfect state transfer, exact nonequilibrium steady states under thermal boundary driving, and detailed entanglement and correlation analyses [2411.09487][2508.20184].

## 1. Definition and scope

A standard open-chain Hamiltonian is
\[
\mathcal{H} = -\frac{1}{2}\sum_{n=0}^{N-1}J_n\bigl(\sigma_n^x\sigma_{n+1}^x+\sigma_n^y\sigma_{n+1}^y\bigr) -\frac{1}{2}\sum_{n=0}^{N} B_n \sigma_n^z ,
\]
with real \(J_n\) and \(B_n\). Here “open” means that the interaction runs only over bonds \((n,n+1)\) for \(n=0,\dots,N-1\), with no periodic link between sites \(N\) and \(0\); “generic inhomogeneous open XX spin chain” means precisely that the couplings \(J_n\) and fields \(B_n\) are otherwise arbitrary real numbers [2508.20184]. An equivalent convention often used in the review literature is
\[
\mathcal{H} = -\frac{1}{2}\sum_{\ell=0}^{N-1} J_\ell\left(\sigma_\ell^x \sigma_{\ell+1}^x+\sigma_\ell^y \sigma_{\ell+1}^y\right) + \frac{1}{2}\sum_{\ell=0}^{N} B_\ell (1+\sigma_\ell^z),
\]
which differs only by an additive constant and sign conventions for the field term [2411.09487].

The model conserves total \(z\)-magnetization. In one convention this is
\[
S^z=\frac12\sum_{n=0}^{N}\sigma_n^z,\qquad [\mathcal H,S^z]=0,
\]
and in the review notation it is stated as
\[
S^z=\sum_{\ell=0}^{N}\sigma_\ell^z, \qquad [\mathcal{H},S^z]=0.
\]
Accordingly, the Hilbert space splits into fixed-excitation sectors [2508.20184][2411.09487].

Two recurrent terminological confusions are worth excluding. First, the homogeneous periodic XX0 chain is a different object: the paper on XX0 phase structure studies the homogeneous, isotropic, periodic XX0 chain at zero magnetic field rather than an inhomogeneous open chain, even though its free-fermion and matrix-integral techniques remain structurally relevant [1612.03463]. Second, in the integrability literature on open XXX/XXZ chains, “inhomogeneous” can denote either site/spectral inhomogeneities \(\{\theta_j\}\) or an additive inhomogeneous term in a \(T\)-\(Q\) equation; the open XXX paper is explicitly of the latter type and sets site inhomogeneities to zero [1307.5049].

The open XX chain also sits as the free-fermion specialization of the open XXZ chain. In the trigonometric parametrization summarized in the SOV literature, the XX chain is obtained at the free-fermion point \(\Delta=0\), corresponding to \(\eta=i\pi/2\) modulo \(i\pi\) [1307.3960].

## 2. Jacobi matrices, fermionization, and orthogonal polynomials

The Jordan–Wigner map sends the spin chain to free fermions. One explicit form is
\[
c_n^\dagger = \sigma_0^z\cdots \sigma_{n-1}^z \sigma_n^+,\qquad c_n = \sigma_0^z\cdots \sigma_{n-1}^z \sigma_n^-,
\]
which transforms the Hamiltonian into
\[
\mathcal{H} = \sum_{n=0}^{N-1} J_n (c_n^\dagger c_{n+1}+c_{n+1}^\dagger c_n) - \sum_{n=0}^{N} B_n c_n^\dagger c_n,
\]
up to an additive constant [2508.20184]. In the review notation,
\[
c_\ell^\dagger = \sigma_0^z \sigma_1^z \cdots \sigma_{\ell-1}^z \sigma_\ell^+, \qquad c_\ell = \sigma_0^z \sigma_1^z \cdots \sigma_{\ell-1}^z \sigma_\ell^-,
\]
and
\[
\mathcal{H} = \sum_{k,\ell=0}^{N}\Lambda_{k\ell} c_k^\dagger c_\ell,
\]
with \(\Lambda\) the tridiagonal Jacobi matrix determined by \(J_\ell\) and \(B_\ell\) [2411.09487].

In the one-excitation sector, the open chain is exactly a Jacobi-matrix problem. In the notation used for perfect state transfer,
\[
\mathcal H e_n = J_{n+1}e_{n+1}+b_n e_n+J_n e_{n-1}, \qquad J_0=J_{N+1}=0,
\]
while in the review notation the amplitudes satisfy
\[
\omega_k \phi_\ell(\omega_k) = J_\ell \phi_{\ell+1}(\omega_k) + B_\ell \phi_\ell(\omega_k) + J_{\ell-1}\phi_{\ell-1}(\omega_k).
\]
This is the defining three-term recurrence of an orthogonal-polynomial system; the one-particle eigenfunctions are therefore orthogonal-polynomial wavefunctions associated with the Jacobi data \((J_\ell,B_\ell)\) [2512.12419][2411.09487].

The diagonalization then proceeds by normal modes
\[
\tilde c_k=\sum_{\ell=0}^{N}\phi_\ell(\omega_k)c_\ell,
\qquad
\mathcal{H}=\sum_{k=0}^{N}\omega_k \tilde c_k^\dagger \tilde c_k,
\]
and the many-body eigenstates are occupation-number states in the \(\tilde c_k\) basis [2411.09487]. This reduction is the common technical substrate behind the symmetry, transport, and entanglement results reviewed below.

A distinguished solvable subclass is obtained by identifying the chain with a finite orthogonal polynomial system. The quasi-exactly solvable construction gives the correspondence
\[
P_{n+1}(E)=(E-b_n)P_n-a_n P_{n-1},
\]
with the associated open XX chain specified by \(J_n=\sqrt{a_{n+1}}\) and \(B_n=b_n\); it classifies six inequivalent canonical families, including two Lamé-related elliptic families [2007.00369]. The perfect-state-transfer literature supplies further explicit orthogonal-polynomial families: Krawtchouk, \(q\)-Racah, and para \(q\)-Racah chains, all realized as explicit inhomogeneous Jacobi matrices [2512.12419].

## 3. Exact \(su(2)\) symmetry and two-fold degeneracy

A particularly strong finite-size result is that a generic inhomogeneous open XX chain acquires an exact \(su(2)\) symmetry after a suitable uniform shift of the transverse field. If the one-particle Hamiltonian has eigenvalues \(\omega_k\), one chooses a particular eigenvalue \(\omega_K\) and shifts
\[
H \mapsto H-\omega_K I ,
\qquad
B_n \mapsto B_n+\omega_K \qquad \text{for all } n,
\]
so that
\[
\omega_K=0.
\]
After this tuning, the zero-mode fermions commute with the Hamiltonian,
\[
[\mathcal{H},\tilde c_K]=0,\qquad [\mathcal{H},\tilde c_K^\dagger]=0,
\]
and generate exact conserved operators [2508.20184].

In fermionic form the generators are
\[
\mathbb S^x=\frac12(\tilde c_K^\dagger+\tilde c_K),\qquad
\mathbb S^y=-\frac{i}{2}(\tilde c_K^\dagger-\tilde c_K),\qquad
\mathbb S^z=-i[\mathbb S^x,\mathbb S^y].
\]
In spin variables they become nonlocal Jordan–Wigner-string operators,
\[
\mathbb S^x=\frac12\sum_{n=0}^N a_n\, \sigma_0^z\cdots \sigma_{n-1}^z \sigma_n^x,
\qquad
\mathbb S^y=\frac12\sum_{n=0}^N a_n\, \sigma_0^z\cdots \sigma_{n-1}^z \sigma_n^y,
\]
with
\[
a_n=\phi_n(\omega_K)=\phi_n(0).
\]
The coefficients are therefore exactly the components of the normalized zero-energy eigenvector of the one-particle Jacobi matrix [2508.20184].

These operators satisfy
\[
[\mathbb S^x,\mathbb S^y]=i\mathbb S^z,\qquad
[\mathbb S^x,\mathbb S^z]=-i\mathbb S^y,\qquad
[\mathbb S^y,\mathbb S^z]=i\mathbb S^x,
\]
and commute with the Hamiltonian. Moreover,
\[
(\mathbb S^x)^2=(\mathbb S^y)^2=(\mathbb S^z)^2=\frac14 I,
\qquad
\vec{\mathbb S}^{\,2}=\frac34 I.
\]
The representation is thus a direct sum of spin-\(\tfrac12\) irreducibles, and every many-body energy level is at least doubly degenerate. The paper emphasizes that this is exact for finite \(N\), not an asymptotic statement [2508.20184].

Operationally, the coefficients \(a_n\) satisfy the zero-mode recurrence
\[
J_{n-1} a_{n-1} - B_n a_n + J_n a_{n+1}=0, \qquad n=1,\dots,N-1,
\]
with boundary equations
\[
-B_0 a_0 + J_0 a_1=0,\qquad J_{N-1}a_{N-1}-B_N a_N=0,
\]
and normalization \(\sum_{n=0}^N a_n^2=1\). This furnishes the explicit symmetry generator for any open chain once the uniform field has been adjusted to create the zero mode [2508.20184].

The orthogonal-polynomial viewpoint makes the symmetry concrete. In the homogeneous case the zero-mode profile is a sine wave; in the Krawtchouk chain it is
\[
a_n=\phi_n(0) = (-1)^n p^{n/2}(1-p)^{(N-n)/2}\sqrt{\binom{N}{n}},
\]
showing that the “inhomogeneous” symmetry generators are computable in closed form in algebraic families [2508.20184].

## 4. Boundary driving, nonequilibrium steady states, and transport

In open-system settings the same chain is coupled at its ends to thermal reservoirs. A global weak-coupling treatment couples boundary spins to bosonic baths and yields a Lindblad equation diagonal in the exact fermionic normal modes. In the review notation, the dissipator is
\[
\mathds D[\rho] = \lambda^2\sum_{k=0}^{N} \left[ d_k\left( \tilde c_k^\dagger \rho \tilde c_k -\frac12\{\rho,\tilde c_k \tilde c_k^\dagger\} \right) + \tilde d_k\left( \tilde c_k \rho \tilde c_k^\dagger -\frac12\{\rho,\tilde c_k^\dagger \tilde c_k\} \right) \right],
\]
with rates controlled by boundary amplitudes \(\phi_0(\omega_k)^2\) and \(\phi_N(\omega_k)^2\) of the normal modes [2411.09487].

The nonequilibrium steady state remains exactly solvable. It is diagonal in the many-body normal-mode basis,
\[
\rho_\infty = \sum_{\vec n}\lambda_{\vec n}\, |\vec n\rrangle\llangle \vec n|,
\]
with
\[
\lambda_{\vec n} = \frac{y_{\vec n}}{\prod_{k=0}^{N}(d_k+\tilde d_k)},
\qquad
y_{\vec n} = \prod_{k=0}^{N}\big(n_k\tilde d_k+(1-n_k)d_k\big).
\]
Thus each mode is populated by a balance of injection and extraction rates determined by the baths and by its boundary weights [2411.09487].

The exact steady heat current from the left bath is
\[
\mathfrak h_L = 2\pi h^2\lambda^2 \sum_k \frac{ \omega_k\, \phi_0(\omega_k)^2\phi_N(\omega_k)^2 \big(n_0(\omega_k)-n_N(\omega_k)\big) }{ \phi_0(\omega_k)^2\big(2n_0(\omega_k)+1\big) + \phi_N(\omega_k)^2\big(2n_N(\omega_k)+1\big) }.
\]
This formula makes the effect of inhomogeneity transparent: transport is controlled by the spectrum \(\omega_k\) and by the boundary values of the eigenfunctions [2411.09487].

Mirror symmetry is a special case. When the couplings and fields are mirror symmetric, the boundary amplitudes satisfy \(\phi_0(\omega_k)^2=\phi_N(\omega_k)^2\), and the current simplifies to the matrix expression
\[
\mathfrak h_L = \pi \lambda^2 \bra{0} \Lambda\, \frac{ \sinh\!\left(\frac{\beta_N-\beta_0}{2}\Lambda\right) }{ \sinh\!\left(\frac{\beta_0+\beta_N}{2}\Lambda\right) } \ket{0}.
\]
In this setting the high-temperature thermal conductivity scales as
\[
\kappa \sim \frac{\pi \lambda^2 h^2 N}{2T} B_0,
\]
so every mirror-symmetric inhomogeneous XX chain is ballistic at high temperature in this model [2411.09487].

Breaking mirror symmetry changes the picture qualitatively. Closed-form current formulas for boundary-driven inhomogeneous XX chains with bosonic baths show that inhomogeneities breaking mirror symmetry significantly reduce both heat and spin conductivities for small temperature differences, and the corresponding current bounds can decay exponentially with \(N\) under linear or random perturbations of mirror-symmetric chains [2408.04818]. The same work also stresses that perfect state transfer does not imply large nonequilibrium current: the mirror-symmetric Krawtchouk chain is a perfect-state-transfer chain, yet in low-temperature linear response its conductivity can decrease exponentially with system size [2408.04818].

A distinct but related global-Lindblad construction produces thermal rectification in asymmetric inhomogeneous XX chains. In that framework the mode occupations are weighted averages of left and right Fermi functions, the current is an exact mode sum, and rectification arises precisely when the effective edge couplings of the normal modes differ, \(g_{L,k}\neq g_{R,k}\), so that swapping hot and cold baths changes the weighting of the transmitting modes. The paper gives analytic examples where rectification remains finite in the thermodynamic limit and where ballistic transport and rectification coexist [2012.04811].

## 5. Spectral engineering and perfect state transfer

The one-excitation reduction makes the open XX chain a natural platform for perfect state transfer (PST). In the \(N+1\)-site formulation,
\[
H= \frac{1}{2}\sum_{\ell=0}^{N-1} J_{\ell+1} \left(\sigma^x_\ell \sigma^x_{\ell+1} + \sigma^y_\ell \sigma^y_{\ell+1}\right) +\frac12\sum_{\ell=0}^{N} b_\ell(\sigma^z_\ell+1),
\]
and the one-excitation Hamiltonian is the Jacobi matrix
\[
\mathcal H e_n = J_{n+1}e_{n+1}+b_n e_n+J_n e_{n-1},
\qquad
J_0=J_{N+1}=0.
\]
PST from site \(0\) to site \(N\) at time \(T\) is characterized by two conditions: persymmetry,
\[
b_n=b_{N-n},\qquad J_n=J_{N-n+1},
\]
and the odd-spacing condition
\[
\epsilon_{x+1}-\epsilon_x=\frac{\pi}{T}M_x,\qquad M_x\in 2\mathbb Z+1,\quad M_x>0.
\]
These are necessary and sufficient in the one-excitation sector [2512.12419].

Orthogonal-polynomial constructions provide explicit PST families. For the \(q\)-Racah chain, the couplings and fields are generated by recurrence coefficients,
\[
J_n=\mu\sqrt{A_{n-1}C_n},\qquad b_n=\mu B_n,
\]
with eigenvalues
\[
\epsilon_x=\mu( q^{-x} + \delta q^{x-N}),\qquad x=0,1,\dots,N.
\]
Under the persymmetric parameter choices
\[
\alpha\beta q^{N+1}=-1,\qquad \delta=\alpha^2 q^{N+1},
\]
or, in the second regime, \(\alpha\beta q^{N+1}=1\), one obtains explicit inhomogeneous XX chains of \(q\)-Racah or para \(q\)-Racah type with PST at time \(T\) once the odd-spacing constraints are satisfied [2512.12419].

The review literature places these constructions in a broader inverse-spectral program. The canonical exactly solvable example is the Krawtchouk chain
\[
J_n=\sqrt{p(1-p)}\sqrt{(n+1)(N-n)},\qquad
B_n=p(N-n)+(1-p)n,
\]
whose case \(p=\tfrac12\) is mirror symmetric and exhibits PST [2411.09487].

Mirror symmetry is therefore a common structural thread between PST and transport, but the two notions are not equivalent. The transport work makes this explicit: the Krawtchouk chain at \(p=\tfrac12\) has PST, yet low-temperature nonequilibrium conductivity can be much smaller than that of the homogeneous chain because NESS transport depends on boundary spectral weights of low-energy modes, not on coherent end-to-end transfer at one special time [2408.04818].

## 6. Entanglement, correlations, and continuum limits

Ground-state entanglement in inhomogeneous open XX chains is likewise controlled by the one-particle Jacobi matrix. For a subsystem \(A\), the restricted correlation matrix is
\[
C_{mn}=\llangle\Psi_0|c_m^\dagger c_n|\Psi_0\rrangle = \sum_{k=0}^{K}\phi_m(\omega_k)\phi_n(\omega_k),
\]
and the entanglement Hamiltonian matrix is
\[
h=\ln\big((1-C)/C\big).
\]
For families tied to Askey-scheme orthogonal polynomials, a commuting tridiagonal Heun operator can often be constructed and provides a good approximation to the entanglement Hamiltonian [2411.09487].

A broad algebraic class is defined by couplings whose squares are polynomial functions of the site index of degree at most four,
\[
J_n^2 = \text{polynomial in } n,\qquad \deg(J_n^2)\le 4.
\]
At half filling and in a constant magnetic field, these chains admit a continuum description as a massless Dirac fermion in a curved static background
\[
ds^2=J(x)^2dt^2-dx^2,
\qquad
\tilde x(x):=\int_0^x\frac{ds}{J(s)}.
\]
The Rényi entropy of an interval then takes the inhomogeneous conformal form
\[
S_\alpha= \frac1{12}(1+\alpha^{-1})\log f(N,\lambda)+c'_\alpha+\cdots,
\]
with the homogeneous lengths replaced by conformal lengths \((\widetilde{\Delta x},\tilde\ell,\tilde N)\). For \(\alpha\ge1\), the open-chain parity oscillations survive and are accurately reproduced by a conformally modified Fagotti–Calabrese term \(\mu_\alpha (-1)^L f(N,\lambda)^{-1/\alpha}\) in the standard half-filled, constant-field regime [2107.12200].

Outside that regime the open boundaries become especially visible. For arbitrary filling and/or inhomogeneous magnetic field, the numerical results show that a block of spins at each end of the chain becomes disentangled from the rest, so appreciable entanglement is concentrated in a central region rather than across the full interval [2107.12200]. The same paper identifies concrete sextic, Krawtchouk, and Lamé chains where the conformal-length picture remains quantitatively accurate to varying degrees.

The quasi-exactly solvable Lamé construction yields a particularly unusual asymptotic law. For the symmetric zero-field Lamé chain
\[
J_n=\sqrt{(n+1)(N-n-1)\Big(N-n-\tfrac32\Big)\Big(n+\tfrac12\Big)}, \qquad B_n=0,
\]
the half-filled boundary-block Rényi entropy satisfies
\[
S_{A,\alpha}\simeq \frac{1}{12}(1+\alpha^{-1})\log\!\left[ \frac{2}{\pi}\lambda(1-\lambda)N\log N \right]+\tilde\gamma_\alpha.
\]
The leading behavior is therefore the standard open-chain \(c=1\) logarithm plus an additional \(\log\log N\) correction, traced to the logarithmic growth of the conformal length in the corresponding curved-space Dirac problem [2007.00369].

Dynamical inhomogeneity can also be imposed through the initial state. For the XX limit of an XXZ chain prepared in a domain-wall profile by a spatially varying longitudinal field and then quenched to the homogeneous XX Hamiltonian, the exact magnetization profile is
\[
m^{z}(j,t) = -\frac{1}{2} + \sum_{m>0}J^{2}_{j+m}\!\left(\sqrt{(Jt)^{2} + \alpha^{2}}\right),
\]
with ballistic wall broadening
\[
W(t)=\sqrt{(Jt)^2+\alpha^2}.
\]
Near the center, \(m^z\) and \(C^{zz}\) relax to homogeneous values, but the transverse correlator retains an oscillatory memory of the initial wall height,
\[
C^{xx}(j,j+n,t\to\infty)\sim C^{xx}_{eq}(n)\cos\!\left(\frac{2\pi n}{\lambda}\right), \qquad \lambda=\frac{2}{m_0}.
\]
The paper states that, for an open XX chain, these formulas should be trusted as exact bulk or pre-reflection results; boundary reflections will eventually modify the late-time global behavior [1002.4446].

The inhomogeneous open XX spin chain is thus best understood as a unifying free-fermion/Jacobi-matrix framework rather than a single special model. Spatially varying couplings and fields alter the spectral data of the tridiagonal one-particle operator, and those altered eigenvalues and boundary amplitudes govern exact symmetry generation, state-transfer design, open-system transport, entanglement structure, and the persistence or suppression of boundary and bulk correlations.

Source: https://www.emergentmind.com/topics/inhomogeneous-open-xx-spin-chain