Open XXZ Spin Chain: Integrable Boundaries

• The open XXZ spin chain is an integrable quantum model with open boundaries, defined through reflection algebras and boundary magnetic fields.
• It employs advanced techniques such as the Bethe ansatz, generalized T-Q relations, and separation of variables to elucidate its spectral properties.
• The model facilitates studies of boundary thermodynamics, transport phenomena, and quantum-state preparation via driven nonequilibrium dynamics.

The open XXZ spin chain is the open-boundary analogue of the XXZ Heisenberg chain, defined on a finite interval or half-infinite line and equipped with boundary terms that may be diagonal, non-diagonal, or, in nonequilibrium settings, dissipative. For spin $\tfrac12$, a standard diagonal-boundary Hamiltonian is

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$

while the integrable structure is formulated through reflection equations and a double-row transfer matrix. In the spin-$\tfrac12$ diagonal case it is described as the paradigmatic example of a quantum integrable model with open boundary conditions, and the broader literature treats it as a laboratory for boundary algebraic Bethe ansatz, generalized $T$-$Q$ systems, separation of variables, qKZ constructions, determinant formulas, and exact nonequilibrium steady states (Dyke et al., 2021, Faldella et al., 2013).

1. Hamiltonians, anisotropy, and boundary data

For the open spin-$\tfrac12$ XXZ chain with diagonal boundary magnetic fields, one common convention is

$${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$

with $L$ spins, anisotropy $\Delta$, and real boundary magnetic fields $h,h'$ at the two ends. In the special boundary choice

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$0

the model has $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$1 quantum-group symmetry, leading to degeneracies analogous to an isotropic $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$2-symmetric chain (Dyke et al., 2021).

In the massive antiferromagnetic regime one frequently writes

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$3

and parametrizes the diagonal boundary fields by

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$4

This is the parametrization used in ground-state overlap and boundary-mode analyses, where the boundary fields are longitudinal and the model is treated within the standard integrable open-chain framework with diagonal boundary conditions (Abetian et al., 2024, Grijalva et al., 2019).

Special anisotropies generate additional exact structure. At the combinatorial point

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$5

with diagonal boundary fields

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$6

the Hamiltonian admits the simple exact eigenvalue

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$7

conjectured to be the non-degenerate ground-state energy for all $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$8 (Hagendorf et al., 2021).

Beyond spin $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$9, the open spin-$\tfrac12$0 XXZ chain is formulated through commuting transfer matrices

$\tfrac12$1

with general integrable boundary parameters

$\tfrac12$2

These parameters encode non-diagonal boundary fields that are generally not simultaneously diagonalizable with the bulk spin operators (Baiyasi et al., 2012). In the half-infinite chain, the thermodynamic-limit Hamiltonian is written with anisotropy $\tfrac12$3 and boundary parameters $\tfrac12$4, and the analysis distinguishes generic non-diagonal, generic diagonal, and one-sided upper or lower non-diagonal boundaries (Baseilhac et al., 2012).

2. Reflection algebra, transfer matrices, and hidden symmetries

Integrability for the open XXZ chain is implemented through Sklyanin’s reflection-algebra construction. For spin $\tfrac12$5, one introduces the bulk monodromy matrix $\tfrac12$6, the reflected monodromy, and boundary monodromies such as

$\tfrac12$7

with commuting transfer matrix

$\tfrac12$8

The Hamiltonian is recovered from the transfer matrix by a logarithmic derivative at a special point in the homogeneous limit (Faldella et al., 2013). In the most general scalar case, the boundary $\tfrac12$9-matrix is

$T$0

which is the origin of non-diagonal boundary terms (Faldella et al., 2013).

For higher spin, the fused transfer matrices satisfy a fusion hierarchy. At root-of-unity anisotropy

$T$1

the hierarchy truncates and yields a finite functional equation of order $T$2 for the fundamental transfer-matrix eigenvalues. This truncation underlies the generalized $T$3-$T$4 analysis of the open spin-$T$5 chain with non-diagonal boundaries (Baiyasi et al., 2012).

In the half-infinite setting, the transfer matrix can be rewritten in terms of a current algebra $T$6 with currents

$T$7

and the boundary condition determines which coideal algebra appears in the first modes. For generic non-diagonal boundaries the hidden symmetry is the q-Onsager algebra, whereas for generic diagonal boundaries it is the augmented q-Onsager algebra (Baseilhac et al., 2012).

A notable open-chain distinction is that only the even charges

$T$8

are true conserved quantities on a finite open chain, while odd charges are conserved only up to boundary terms. This distinction is central in the classification of integrable deformations and in the appearance of boundary reflection phases in asymptotic open-chain Bethe equations (Beisert et al., 2013).

3. Bethe ansatz, generalized $T$9-$Q$0 systems, and separation of variables

For diagonal boundaries, eigenstates can be constructed by boundary algebraic or coordinate Bethe ansatz. In the coordinate formulation, a Bethe state with $Q$1 down spins is expanded over ordered spin-down positions, and for the open chain the wave function involves both permutations of the Bethe momenta and sign flips $Q$2. Consequently, for fixed $Q$3 the open-chain Bethe wave function contains

$Q$4

terms, compared with only $Q$5 terms for the periodic chain. The Bethe roots obey

$Q$6

and the energy is

$Q$7

In the diagonal-boundary massive regime, the boundary algebraic Bethe ansatz uses states

$Q$8

with energy

$Q$9

These formulations are standard starting points for spectral and overlap problems (Dyke et al., 2021, Abetian et al., 2024).

Generic non-diagonal boundaries remove the conventional ferromagnetic reference state. This is precisely the regime where quantum separation of variables becomes decisive. In the inhomogeneous open spin-$\tfrac12$0 XXZ chain with the most general integrable boundary terms, the transfer-matrix spectrum is simple and every eigenvalue is characterized by a system of $\tfrac12$1 quadratic equations. The eigenstates are reconstructed in an SOV basis from discrete Baxter-like relations, and scalar products of separate states reduce to determinants (Faldella et al., 2013, Niccoli, 2012).

At root-of-unity anisotropy, the non-diagonal open spin-$\tfrac12$2 chain admits generalized $\tfrac12$3-$\tfrac12$4 relations involving two independent $\tfrac12$5-functions,

$\tfrac12$6

with

$\tfrac12$7

The construction works for special boundary families with at most two arbitrary boundary parameters, and numerical evidence for $\tfrac12$8 shows completeness for small chains (Baiyasi et al., 2012). A related arbitrary-spin analysis at

$\tfrac12$9

with $${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$0 restricted to odd integers, derives $${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$1-$${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$2 relations and Bethe equations for selected non-diagonal boundary choices and reproduces the complete spectra for representative spin-$${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$3 and spin-1 chains (Murgan et al., 2014).

A different constrained regime is governed by the phantom roots criterion,

$${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$4

under which the Hilbert space splits into two invariant subspaces $${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$5 and $${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$6, each described by homogeneous Bethe equations. The corresponding eigenstates are expanded in factorized chiral-shock vectors, and the resulting chiral coordinate Bethe ansatz produces bulk scattering matrices, boundary reflection matrices, and energies in a basis built from spin-helix-like local states (Zhang et al., 2021).

At the combinatorial point $${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$7, a boundary qKZ solution $${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$8 yields an additional exact route. The generalized overlap $${\cal H} = -\tfrac{1}{2}\sum_{n=0}^{L-2} \left( \sigma^{x}_{n}\sigma^{x}_{n+1} + \sigma^{y}_{n}\sigma^{y}_{n+1} + \Delta\, \sigma^{z}_{n}\sigma^{z}_{n+1} \right) -\tfrac{1}{2}\left( h\, \sigma^{z}_{0} + h'\, \sigma^{z}_{L-1} \right),$$9 satisfies symmetry, reflection, and reduction properties and becomes expressible in terms of symplectic characters associated to a double-staircase partition. In the homogeneous limit this produces explicit determinant formulas for overlaps and fidelity observables (Hagendorf et al., 2021).

4. Scalar products, overlaps, and fidelity observables

For non-diagonal boundaries, determinant technology survives in several complementary guises. In the vertex–face formulation with a Drinfeld twist or factorizing $L$0-matrix, the face-type boundary matrices become diagonal and the transformed pseudo-particle operators are polarization free. Because the open chain then has two sets of Bethe states, one obtains determinant formulas for four families of scalar products and, in the on-shell limit, Gaudin-type determinant formulas for the norms (Yang et al., 2010).

Within the SoV framework, scalar products of separate states are first written as dressed Vandermonde determinants with intricate inhomogeneity dependence and are then transformed into alternative determinant forms in which the homogeneous limit can be taken straightforwardly. When one state is on shell, the result becomes a generalized Slavnov-type determinant; for special boundary choices relevant to the half-infinite chain, the formulas simplify further and are intended for form-factor and correlation-function calculations (Kitanine et al., 2018). This agrees with the broader SOV program in which determinant scalar products are the basis for norms, matrix elements of local operators, and quasi-local form factors (Faldella et al., 2013, Niccoli, 2012).

Off-shell scalar products for diagonal-boundary open chains can also be characterized directly as solutions of a pair of functional equations obtained from the reflection algebra. The outcome is a multiple contour integral representation whose notable feature is that the homogeneous limit is trivial in that representation (Galleas, 2014).

Two overlap problems have received exact treatment. In the massive antiferromagnetic regime, the normalized overlap between ground states before and after changing one boundary field admits an exact thermodynamic-limit formula. When the two states have mismatched boundary roots, the overlap is exponentially small,

$L$1

whereas in the real-root and compatible boundary-root cases it reduces to a closed expression in double $L$2-Pochhammer symbols, up to exponentially small corrections (Abetian et al., 2024). At $L$3, the logarithmic bipartite fidelity

$L$4

is obtained exactly at finite size from a determinant formula for the overlap $L$5. Its large-$L$6 asymptotics contains the universal term

$L$7

and the calculation confirms the conformal-field-theory prediction with $L$8 and $L$9 for the $\Delta$0 correction (Hagendorf et al., 2021).

5. Boundary thermodynamics, surface terms, and boundary-localized modes

The thermodynamics of the open XXZ spin-$\Delta$1 chain differs from the periodic problem by an $\Delta$2 boundary contribution. For system size $\Delta$3,

$\Delta$4

and $\Delta$5 is the surface free energy. A finite-Trotter expression for the boundary free energy can be reinterpreted as the partition function of the six-vertex model with reflecting ends; Tsuchiya’s determinant formula then permits the infinite-Trotter limit. The resulting exact representation for $\Delta$6 is written in terms of the same nonlinear integral equation that governs the periodic XXZ chain, and the boundary magnetization follows from

$\Delta$7

In the massless regime the leading thermal correction to the boundary magnetization is $\Delta$8, not $\Delta$9 (Kozlowski et al., 2012).

At zero temperature in the massive antiferromagnetic regime, the ground state may contain a boundary-localized Bethe excitation, the boundary root,

$h,h'$0

For even chain length and equal boundary fields satisfying

$h,h'$1

the spectrum is gapped and the ground state is doubly degenerate up to exponentially small corrections in $h,h'$2. The two lowest states differ by whether the boundary root is localized near the left or right edge (Grijalva et al., 2019).

This boundary root controls both static and dynamical edge observables. The thermodynamic boundary magnetization decomposes as

$h,h'$3

and, strikingly, the magnetization at the left edge can depend on the right boundary field even in the half-infinite-chain limit because the ground-state presence of the boundary root depends on both boundaries. The same quantity determines the long-time plateau of the zero-temperature boundary autocorrelation: $h,h'$4 For odd $h,h'$5, by contrast, the ground-state structure and the location of the boundary-magnetization discontinuity are different, and the quasi-degenerate even-$h,h'$6 scenario does not occur in the same way (Grijalva et al., 2019).

6. Boundary driving, transport, and quantum-state preparation

A second major branch of the subject concerns boundary-driven open XXZ chains in Lindblad form. For a homogeneous nearest-neighbor Heisenberg XXZ spin-$h,h'$7 chain driven only at the edges, a weak-coupling expansion of the nonequilibrium steady state produces an explicit matrix-product operator $h,h'$8 satisfying

$h,h'$9

The Hermitian combination $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$00 becomes, in the thermodynamic limit, an exact pseudolocal conservation law with nonzero overlap with the current. Through Mazur’s inequality this yields a rigorous positive lower bound on the high-temperature spin Drude weight for

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$01

and the bound is a nonvanishing fractal function of the anisotropy (Prosen, 2011).

Under strong boundary dissipation in the easy-plane regime, the nonequilibrium steady state can instead approach a pure spin-helix state,

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$02

provided the anisotropy is tuned to the resonance condition

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$03

In that Zeno-limit regime the entropy vanishes, the energy current is zero, and the spin current becomes

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$04

so the transport crosses over from diffusive $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$05 at weak dissipation to ballistic $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$06 at strong dissipation (Popkov et al., 2017).

An exact nonequilibrium steady state is also known for a one-end driven geometry in which site $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$07 is coupled to a source bath $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$08 and the right end carries an arbitrary coherent field

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$09

The steady state has matrix-product form

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$10

with $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$11 an infinite-dimensional $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$12 Lax operator and $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$13 determined by a three-term recurrence fixed by the boundary field (Popkov et al., 17 Apr 2026).

Exact eigenstates of the open XXZ chain also appear in quantum-information settings. For real solutions of the Bethe equations, there is a probabilistic quantum algorithm that prepares Bethe states of the open chain on a quantum computer. A Bethe state of $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$14 spins with $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$15 down spins contains

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$16

terms, while the algorithm uses

$$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$17

qubits, namely $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$18 system qubits, $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$19 permutation-label qubits, and $$\mathbf{H}\equiv \mathbf{H}_{h^-,h^+}= \sum_{i=1}^{L-1} \left[ \sigma_i^x \sigma_{i+1}^x +\sigma_i^y \sigma_{i+1}^y +\Delta(\sigma_i^z \sigma_{i+1}^z-1) \right] +h^- \sigma_1^z+h^+\sigma_L^z,$$20 faucet qubits. The success probability decreases with the number of down spins, but amplitude amplification can be used to boost it (Dyke et al., 2021).

The open XXZ spin chain therefore encompasses several technically distinct but structurally connected objects: the finite and half-infinite Hamiltonian chain with integrable reflection boundaries, generalized higher-spin and non-diagonal boundary systems, boundary-driven Lindblad chains, and explicit overlap or fidelity problems. Across these settings, the recurring mathematical themes are reflection algebras, boundary scattering, determinant formulas, and boundary-sensitive spectral data. This suggests that the open chain is not merely the periodic XXZ chain with endpoints added, but a boundary-dominated integrable system whose most characteristic phenomena arise precisely because the ends are dynamical participants.