Inhomogeneous Open XX Spin Chain

• The inhomogeneous open XX spin chain is a finite spin-½ chain with site-dependent couplings and magnetic fields, reducible to a quadratic free-fermion system via Jordan–Wigner transformation.
• Its reduction to a tridiagonal Jacobi matrix allows explicit diagonalization, underpinning analyses of perfect state transfer, nonequilibrium steady states, and spectral properties.
• The model exhibits exact finite-size su(2) symmetry and provides a versatile framework for investigating transport phenomena, entanglement structure, and continuum limits.

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 (Bernard et al., 2024, Crampé et al., 27 Aug 2025).

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 (Crampé et al., 27 Aug 2025). An equivalent convention often used in the review literature 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 ,$$0

which differs only by an additive constant and sign conventions for the field term (Bernard et al., 2024).

The model conserves total $$\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 ,$$1-magnetization. In one convention this 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 ,$$2

and in the review notation it is stated as

$$\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 ,$$3

Accordingly, the Hilbert space splits into fixed-excitation sectors (Crampé et al., 27 Aug 2025, Bernard et al., 2024).

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 (Saeedian et al., 2016). Second, in the integrability literature on open XXX/XXZ chains, “inhomogeneous” can denote either site/spectral inhomogeneities $$\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 ,$$4 or an additive inhomogeneous term in a $$\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 ,$$5-$$\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 ,$$6 equation; the open XXX paper is explicitly of the latter type and sets site inhomogeneities to zero (Nepomechie, 2013).

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 $$\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 ,$$7, corresponding to $$\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 ,$$8 modulo $$\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 ,$$9 (Faldella et al., 2013).

2. Jacobi matrices, fermionization, and orthogonal polynomials

The Jordan–Wigner map sends the spin chain to free fermions. One explicit form is

$J_n$0

which transforms the Hamiltonian into

$J_n$1

up to an additive constant (Crampé et al., 27 Aug 2025). In the review notation,

$J_n$2

and

$J_n$3

with $J_n$4 the tridiagonal Jacobi matrix determined by $J_n$5 and $J_n$6 (Bernard et al., 2024).

In the one-excitation sector, the open chain is exactly a Jacobi-matrix problem. In the notation used for perfect state transfer,

$J_n$7

while in the review notation the amplitudes satisfy

$J_n$8

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_n$9 (Crampe et al., 13 Dec 2025, Bernard et al., 2024).

The diagonalization then proceeds by normal modes

$B_n$0

and the many-body eigenstates are occupation-number states in the $B_n$1 basis (Bernard et al., 2024). 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

$B_n$2

with the associated open XX chain specified by $B_n$3 and $B_n$4; it classifies six inequivalent canonical families, including two Lamé-related elliptic families (Finkel et al., 2020). The perfect-state-transfer literature supplies further explicit orthogonal-polynomial families: Krawtchouk, $B_n$5-Racah, and para $B_n$6-Racah chains, all realized as explicit inhomogeneous Jacobi matrices (Crampe et al., 13 Dec 2025).

3. Exact $B_n$7 symmetry and two-fold degeneracy

A particularly strong finite-size result is that a generic inhomogeneous open XX chain acquires an exact $B_n$8 symmetry after a suitable uniform shift of the transverse field. If the one-particle Hamiltonian has eigenvalues $B_n$9, one chooses a particular eigenvalue $(n,n+1)$0 and shifts

$(n,n+1)$1

so that

$(n,n+1)$2

After this tuning, the zero-mode fermions commute with the Hamiltonian,

$(n,n+1)$3

and generate exact conserved operators (Crampé et al., 27 Aug 2025).

In fermionic form the generators are

$(n,n+1)$4

In spin variables they become nonlocal Jordan–Wigner-string operators,

$(n,n+1)$5

with

$(n,n+1)$6

The coefficients are therefore exactly the components of the normalized zero-energy eigenvector of the one-particle Jacobi matrix (Crampé et al., 27 Aug 2025).

These operators satisfy

$(n,n+1)$7

and commute with the Hamiltonian. Moreover,

$(n,n+1)$8

The representation is thus a direct sum of spin-$(n,n+1)$9 irreducibles, and every many-body energy level is at least doubly degenerate. The paper emphasizes that this is exact for finite $n=0,\dots,N-1$0, not an asymptotic statement (Crampé et al., 27 Aug 2025).

Operationally, the coefficients $n=0,\dots,N-1$1 satisfy the zero-mode recurrence

$n=0,\dots,N-1$2

with boundary equations

$n=0,\dots,N-1$3

and normalization $n=0,\dots,N-1$4. This furnishes the explicit symmetry generator for any open chain once the uniform field has been adjusted to create the zero mode (Crampé et al., 27 Aug 2025).

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

$n=0,\dots,N-1$5

showing that the “inhomogeneous” symmetry generators are computable in closed form in algebraic families (Crampé et al., 27 Aug 2025).

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

$n=0,\dots,N-1$6

with rates controlled by boundary amplitudes $n=0,\dots,N-1$7 and $n=0,\dots,N-1$8 of the normal modes (Bernard et al., 2024).

The nonequilibrium steady state remains exactly solvable. It is diagonal in the many-body normal-mode basis,

$n=0,\dots,N-1$9

with

$N$0

Thus each mode is populated by a balance of injection and extraction rates determined by the baths and by its boundary weights (Bernard et al., 2024).

The exact steady heat current from the left bath is

$N$1

This formula makes the effect of inhomogeneity transparent: transport is controlled by the spectrum $N$2 and by the boundary values of the eigenfunctions (Bernard et al., 2024).

Mirror symmetry is a special case. When the couplings and fields are mirror symmetric, the boundary amplitudes satisfy $N$3, and the current simplifies to the matrix expression

$N$4

In this setting the high-temperature thermal conductivity scales as

$N$5

so every mirror-symmetric inhomogeneous XX chain is ballistic at high temperature in this model (Bernard et al., 2024).

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$6 under linear or random perturbations of mirror-symmetric chains (Bernard et al., 2024). 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 (Bernard et al., 2024).

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, $N$7, 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 (Silva et al., 2020).

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$8-site formulation,

$N$9

and the one-excitation Hamiltonian is the Jacobi matrix

$0$0

PST from site $0$1 to site $0$2 at time $0$3 is characterized by two conditions: persymmetry,

$0$4

and the odd-spacing condition

$0$5

These are necessary and sufficient in the one-excitation sector (Crampe et al., 13 Dec 2025).

Orthogonal-polynomial constructions provide explicit PST families. For the $0$6-Racah chain, the couplings and fields are generated by recurrence coefficients,

$0$7

with eigenvalues

$0$8

Under the persymmetric parameter choices

$0$9

or, in the second regime, $J_n$0, one obtains explicit inhomogeneous XX chains of $J_n$1-Racah or para $J_n$2-Racah type with PST at time $J_n$3 once the odd-spacing constraints are satisfied (Crampe et al., 13 Dec 2025).

The review literature places these constructions in a broader inverse-spectral program. The canonical exactly solvable example is the Krawtchouk chain

$J_n$4

whose case $J_n$5 is mirror symmetric and exhibits PST (Bernard et al., 2024).

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 $J_n$6 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 (Bernard et al., 2024).

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 $J_n$7, the restricted correlation matrix is

$J_n$8

and the entanglement Hamiltonian matrix is

$J_n$9

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 (Bernard et al., 2024).

A broad algebraic class is defined by couplings whose squares are polynomial functions of the site index of degree at most four,

$B_n$0

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

$B_n$1

The Rényi entropy of an interval then takes the inhomogeneous conformal form

$B_n$2

with the homogeneous lengths replaced by conformal lengths $B_n$3. For $B_n$4, the open-chain parity oscillations survive and are accurately reproduced by a conformally modified Fagotti–Calabrese term $B_n$5 in the standard half-filled, constant-field regime (Finkel et al., 2021).

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 (Finkel et al., 2021). 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

$B_n$6

the half-filled boundary-block Rényi entropy satisfies

$B_n$7

The leading behavior is therefore the standard open-chain $B_n$8 logarithm plus an additional $B_n$9 correction, traced to the logarithmic growth of the conformal length in the corresponding curved-space Dirac problem (Finkel et al., 2020).

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

$$\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 ,$$00

with ballistic wall broadening

$$\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 ,$$01

Near the center, $$\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 ,$$02 and $$\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 ,$$03 relax to homogeneous values, but the transverse correlator retains an oscillatory memory of the initial wall height,

$$\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 ,$$04

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 (Lancaster et al., 2010).

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.