Anisotropic Heisenberg SSH Spin Chains

• Anisotropic Heisenberg SSH spin chains are one-dimensional quantum spin systems that combine bond alternation with directional anisotropy to yield dimerized and topologically non-trivial phases.
• These chains are modeled via SSH-XXZ Hamiltonians and related variants, where Berry phase diagnostics, solitonic defects, and gap openings characterize the system’s phase transitions.
• Experimental setups and computational frameworks, including Jordan–Wigner and tensor-network methods, demonstrate controlled quantum-state transfer and anomalous transport within these structures.

Anisotropic Heisenberg SSH spin chains are one-dimensional quantum spin systems in which the bond alternation of the Su–Schrieffer–Heeger construction is embedded into Heisenberg exchange with directional anisotropy. In the simplest spin-$\tfrac12$ realizations, nearest-neighbor couplings alternate between two values while the $xy$ and $z$ exchanges need not coincide; related variants include bond-dependent anisotropy, next-nearest-neighbor and chiral interactions, long-range dipolar terms, and higher-spin or single-ion-anisotropy extensions. Across these realizations, the subject links dimerization, Berry-phase topology, hidden order, solitonic domain walls, anomalous transport, and quantum-state-transfer functionality (Li et al., 2018, Paul et al., 2017, Moragues et al., 25 Aug 2025, Qiao et al., 2019).

1. Core Hamiltonians and model classes

A canonical spin-$\tfrac12$ SSH-XXZ chain is the dimerized bond-dependent Hamiltonian

$$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$

where $\eta$ controls the dimerization and $\Delta_1,\Delta_2$ control the XXZ anisotropy on odd and even bonds (Li et al., 2018). A widely used specialization is the bond-alternating anisotropic Heisenberg chain

$$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$

with $J_1>0$ fixed antiferromagnetic and $J_2$ either antiferromagnetic or ferromagnetic, and with bond alternation measured by $xy$0 (Paul et al., 2017).

A ferromagnetic convention used for quantum-channel studies writes the SSH-Heisenberg Hamiltonian as

$xy$1

with $xy$2 the SSH-XX limit and $xy$3 the isotropic Heisenberg limit (Moragues et al., 25 Aug 2025). Because different authors attach the unit cell to different bonds, the sign of $xy$4 labeling the topological regime is not universal. This suggests a convention change rather than a substantive discrepancy.

Beyond nearest-neighbor dimerized XXZ chains, anisotropic Heisenberg SSH-type physics also appears in integrable $xy$5-$xy$6 constructions with chirality and twisted boundaries. One exactly solved example has XXZ nearest-neighbor exchange, isotropic next-nearest-neighbor coupling, and staggered scalar chirality, with nearest-neighbor parameters

$xy$7

and an antiperiodic boundary condition implemented by conjugation with $xy$8 (Qiao et al., 2019).

2. Dimerization, Berry phases, and bond parity

The topological characterization used most directly in anisotropic SSH spin chains is the Berry phase under twisted boundary conditions. In the SSH-XXZ chain, the topologically non-trivial phase is identified by Berry phase $xy$9, while $z$0 denotes the trivial phase. The discretized expression is

$z$1

with $z$2 the twist angle threaded through the chain (Li et al., 2018). In the same work, the entanglement-based dimerization measure

$z$3

distinguishes even-bond from odd-bond dimerization, where $z$4 and $z$5 are von Neumann entropies of nearest-neighbor two-site reduced density matrices on even and odd bonds. The topologically non-trivial phase is primarily associated with even-bond dimerization, and in the anisotropy-only SSH-like model

$z$6

is the condition for a topological non-trivial phase.

A central structural result is that the topological phase is not determined by periodicity but by bond parity. In aperiodic SSH-like XX chains built from contiguous blocks of $z$7 and $z$8 bonds, odd-odd block parity preserves $z$9 for $\tfrac12$0, even-even block parity forces $\tfrac12$1 for all $\tfrac12$2, and mixed parity shifts the transition away from $\tfrac12$3 (Li et al., 2018). This replaces translational invariance with a real-space parity criterion.

Uniform XXZ anisotropy does not displace the conventional XXZ critical points at $\tfrac12$4; instead, dimerization changes the internal bond structure of the antiferromagnetic and Tomonaga–Luttinger-liquid regimes. Away from the fully polarized region $\tfrac12$5, the topological non-trivial phase is essentially independent of $\tfrac12$6, while bond-dependent anisotropy can enlarge or suppress it. Positive $\tfrac12$7 strengthens even-bond antiferromagnetism and expands the $\tfrac12$8 region, whereas large positive $\tfrac12$9 or strongly negative $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$0 suppresses it (Li et al., 2018).

The recent quantum-channel formulation supplements Berry-phase diagnostics with direct edge-state localization. In the one-excitation sector, the edge-localization weight

$$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$1

measures how strongly the boundary excitation overlaps with the two most relevant eigenstates. Strong localization appears mainly in the topological-like sector for moderate $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$2, while near $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$3 no edge states appear for any $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$4 (Moragues et al., 25 Aug 2025).

3. Gapped phases, hidden order, and solitonic defects

In the bond-alternating spin-$$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$5 anisotropic Heisenberg chain, a spin gap develops as soon as non-uniformity in alternating bond strength is introduced in the AFM-AFM case, and remains non-zero throughout the AFM-FM case. Exact diagonalization, bond-operator theory, and Jordan–Wigner analysis show that dimer order and string orders coexist in the ground state. The string correlator is non-zero across most of the anisotropic region $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$6, and the Haldane phase exists throughout most of that region, excluding the uniform AFM line $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$7 and the special FM XY point $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$8 (Paul et al., 2017). In this sense, anisotropic Heisenberg SSH chains interpolate between spin-Peierls-like dimerization and effective spin-1 Haldane physics.

The same work makes the SSH analogy explicit at the one-particle level. In the XY limit $$H = -\sum_{j=1}^N\Big[(1+\eta)\big(\sigma_{2j-1}^x\sigma_{2j}^x+\sigma_{2j-1}^y\sigma_{2j}^y+\Delta_1\sigma_{2j-1}^z\sigma_{2j}^z\big) +(1-\eta)\big(\sigma_{2j}^x\sigma_{2j+1}^x+\sigma_{2j}^y\sigma_{2j+1}^y+\Delta_2\sigma_{2j}^z\sigma_{2j+1}^z\big)\Big],$$9, the Jordan–Wigner mapping yields a two-band fermionic structure

$\eta$0

with alternating effective hoppings inherited from $\eta$1 and $\eta$2. At $\eta$3 the spectrum is gapless; once bond alternation is introduced, a gap opens, just as in the fermionic SSH chain (Paul et al., 2017).

A complementary realization is the anisotropic $\eta$4-$\eta$5 sawtooth chain, which at the symmetric point $\eta$6 has two exact dimerized ground states $\eta$7 and $\eta$8, with exact ground-state energy

$\eta$9

These two vacua are related by translation and support kink and antikink domain walls (Paul et al., 2019). The kink is dispersionless and gapless at the symmetric point, whereas the antikink is dispersive and defines the spin gap. For 1-, 5-, and 9-site antikink clusters, analytic dispersions $\Delta_1,\Delta_2$0 can be derived; a variational combination of these clusters yields a minimized dispersion $\Delta_1,\Delta_2$1, and the variational spin gap is

$\Delta_1,\Delta_2$2

At $\Delta_1,\Delta_2$3, the variational estimate is $\Delta_1,\Delta_2$4, close to the numerically extrapolated $\Delta_1,\Delta_2$5 (Paul et al., 2019).

This solitonic picture differs from the Haldane-type hidden order of the bond-alternating AFM-FM chain. The sawtooth chain realizes a frustration-induced dimerized phase with SSH/Majumdar–Ghosh-type domain walls, whereas the bond-alternating AFM-FM chain realizes a Haldane phase with coexisting string and dimer order. The shared structural element is the existence of two dimerized sectors and low-energy defects interpolating between them.

4. Dynamics, transport, and quantum-state transfer

In the easy-axis regime $\Delta_1,\Delta_2$6, dimerization acts as an integrability-breaking perturbation of the XXZ chain. The staggered Hamiltonian perturbation

$\Delta_1,\Delta_2$7

is precisely an SSH-type alternation of the exchange couplings (Nardis et al., 2021). For large anisotropy, spin transport is subdiffusive with dynamical exponent $\Delta_1,\Delta_2$8 up to a timescale parametrically long in $\Delta_1,\Delta_2$9; in the limit of infinite anisotropy, subdiffusion persists at all times, while for finite large $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$0 a late crossover to diffusion occurs with diffusion constant scaling as $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$1, independent of the strength of the integrability breaking. This establishes dimerized XXZ chains as a concrete setting in which SSH-type modulation reshapes hydrodynamics.

Ultracold-atom realizations of anisotropic Heisenberg chains clarify a different dynamical aspect: transverse spins need not decay by transport alone. For the XXZ Hamiltonian

$$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$2

anisotropy produces local transverse dephasing, while superexchange generates an effective field $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$3 whose inhomogeneity across chains, halving at edges, and fluctuations in the presence of mobile holes produce additional dephasing channels (Jepsen et al., 2021). This suggests that an SSH-dimerized extension with bond-dependent $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$4, $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$5, and $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$6 will inherit not only split magnon bands and edge-localized modes but also structured dephasing tied to strong and weak bonds.

Quantum communication studies expose the tension between topology and usable dynamics. In the anisotropic Heisenberg SSH chain used as a quantum channel, one-excitation transfer is quantified by

$$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$7

while two-excitation and two-qubit protocols use analogous amplitudes in the one- and two-excitation sectors (Moragues et al., 25 Aug 2025). Edge states in the topological-like regime are robust to static disorder, but their energy splitting decays exponentially with $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$8, so transfer becomes very slow. By contrast, trivial regimes can yield much faster and often practically superior transmission. The same study finds that dipolar long-range interactions break total magnetization conservation, strongly degrading transfer unless a sufficiently large Zeeman field suppresses leakage, and that optimal control can accelerate transfer dramatically in trivial regimes but becomes increasingly difficult in strongly localized topological regimes (Moragues et al., 25 Aug 2025). A related study of tailored Heisenberg and XXZ chains with site-dependent couplings reaches a complementary conclusion: although $$H=\sum_{i=1}^{N/2}\Big[ J_1\big( S_{2i-1}^xS_{2i}^x+S_{2i-1}^yS_{2i}^y+\Delta S_{2i-1}^zS_{2i}^z\big) + J_2\big( S_{2i}^xS_{2i+1}^x+S_{2i}^yS_{2i+1}^y+\Delta S_{2i}^zS_{2i+1}^z\big) \Big],$$9 chains can give high average transfer probability with moderate couplings, isotropic Heisenberg chains are the most robust option under static disorder when worst-case performance is emphasized (Serra et al., 2021).

5. Exact solutions and computational frameworks

Several analytic and tensor-network methods have become standard for anisotropic Heisenberg SSH chains because the topic combines dimerization, boundary sensitivity, and interaction effects. For the bond-alternating spin-$J_1>0$0 chain, the bond-operator formalism treats each strong bond as a dimer carrying one singlet and three triplets subject to the local constraint

$J_1>0$1

and yields triplon dispersions

$J_1>0$2

while the Jordan–Wigner transformation maps the XY limit to an SSH-like free-fermion problem (Paul et al., 2017). Exact diagonalization with finite-size extrapolation then resolves the spin gap and the string and dimer order parameters.

Twisted-boundary integrable models require a different toolkit. In the anisotropic $J_1>0$3-$J_1>0$4 chain with chirality and antiperiodic boundary conditions, the twist breaks $J_1>0$5 symmetry, so standard algebraic Bethe Ansatz does not apply. The spectrum is instead obtained by off-diagonal Bethe Ansatz through an inhomogeneous $J_1>0$6-$J_1>0$7 relation of the form

$J_1>0$8

with corresponding inhomogeneous Bethe equations (Qiao et al., 2019). This exact framework gives ground-state energies, twisted boundary energies, and root-string structures, and it turns boundary twists into quantitatively controlled observables rather than formal probes.

For higher-spin anisotropic Heisenberg chains, matrix-product-state methods with periodic boundary conditions become essential when the observables themselves are twist sensitive. In the spin-1 anisotropic Heisenberg chain with single-ion anisotropy and transverse field, the supersolid phase is diagnosed by coexistence of the structure factor

$J_1>0$9

and the spin stiffness

$J_2$0

Because $J_2$1 vanishes trivially with open boundaries, periodic boundary conditions are indispensable (Rossini et al., 2011). The same exposition proposes that in SSH-type anisotropic Heisenberg chains these diagnostics can be transferred directly to dimerized Hamiltonians, with $J_2$2 detecting density-wave or dimer order and $J_2$3 measuring off-diagonal order.

6. Microscopic derivations, realizations, and extensions

A microscopic route from correlated electrons to SSH-like anisotropic spin chains is provided by the half-filled asymmetric ionic Hubbard chain with alternating on-site interaction. In the strong-coupling limit, its effective Hamiltonian is a spin-$J_2$4 anisotropic XXZ chain. At second order one obtains

$J_2$5

with

$J_2$6

and at fourth order the model acquires alternating next-nearest-neighbor couplings, alternating three-spin terms, and both uniform and staggered fields (Grusha et al., 2015). This provides a concrete mechanism by which sublattice asymmetry and spin-dependent hopping generate anisotropic, effectively dimerized spin chains.

Several experimental platforms now supply partial realizations of these Hamiltonians. Ultracold atoms implement XXZ superexchange and effective magnetic fields (Jepsen et al., 2021). Magnetic atom chains on surfaces, nanographene spin chains, and transmon arrays have motivated direct studies of anisotropic Heisenberg SSH channels with dimerized nearest-neighbor exchange, dipolar couplings, and local control fields (Moragues et al., 25 Aug 2025). Tailored exchange profiles in isotropic and anisotropic Heisenberg chains up to $J_2$7 show that near-perfect transfer can be achieved without time-dependent global control, albeit with markedly different robustness properties across anisotropies (Serra et al., 2021).

The subject also extends naturally beyond spin-$J_2$8 nearest-neighbor SSH chains. One proposed spin-1 SSH-type generalization introduces alternating couplings on top of XXZ exchange, single-ion anisotropy, and transverse field,

$J_2$9

and retains the supersolid observables $xy$00 and $xy$01 as natural diagnostics (Rossini et al., 2011). A plausible implication is that anisotropic Heisenberg SSH chains form a common language across dimerized topological phases, hidden-order phases, and coexistence phases with both diagonal and off-diagonal order.

The modern picture is therefore deliberately non-reductive. Dimerization can create Berry-phase topology, but bond parity can be more fundamental than strict periodicity (Li et al., 2018). Hidden order can coexist with local dimer order in anisotropic bond-alternating chains (Paul et al., 2017). Solitonic domain walls can be the correct low-energy carriers in frustrated geometries (Paul et al., 2019). Transport can become subdiffusive once easy-axis XXZ dynamics is weakly dimerized (Nardis et al., 2021). And topological edge localization, although robust, need not optimize quantum communication performance (Moragues et al., 25 Aug 2025). Anisotropic Heisenberg SSH spin chains are best understood not as a single model but as a family of dimerized interacting spin systems in which topology, frustration, anisotropy, and control compete on equal footing.