--- title: 'Open XXZ Spin Chain: Integrable Boundaries' url: https://www.emergentmind.com/topics/open-xxz-spin-chain type: topic --- # Open XXZ Spin Chain: Integrable Boundaries 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 [2109.05607][1307.3960]. ## 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 \[ h'=-h=\frac{1}{2}(q-q^{-1}), \qquad \Delta=\frac{1}{2}(q+q^{-1}), \] the model has \(U_q(su(2))\) quantum-group symmetry, leading to degeneracies analogous to an isotropic \(su(2)\)-symmetric chain [2109.05607]. In the massive antiferromagnetic regime one frequently writes \[ \Delta=\cosh\zeta>1,\qquad \zeta>0, \] and parametrizes the diagonal boundary fields by \[ h^\pm=-\sinh\zeta\,\coth\xi^\pm . \] 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 [2409.15194][1901.10932]. Special anisotropies generate additional exact structure. At the combinatorial point \[ \Delta=-\frac12, \] with diagonal boundary fields \[ p = \frac12\left(\frac12 - x\right), \qquad \bar p = \frac12\left(\frac12-\frac1x\right), \] the Hamiltonian admits the simple exact eigenvalue \[ E_0 = -\frac{3N-1}{4}-\frac{(1-x)^2}{2x}, \] conjectured to be the non-degenerate ground-state energy for all \(x>0\) [2111.15223]. Beyond spin \(\tfrac12\), the open spin-\(s\) XXZ chain is formulated through commuting transfer matrices \[ t^{(j,s)}(u) = \operatorname{tr}_a\!\left[ K^{+ (j)}_a(u)\, T^{(j,s)}_a(u)\, K^{-(j)}_a(u)\, \widehat T^{(j,s)}_a(u) \right], \] with general integrable boundary parameters \[ (\alpha_-,\beta_-,\theta_-), \qquad (\alpha_+,\beta_+,\theta_+). \] These parameters encode non-diagonal boundary fields that are generally not simultaneously diagonalizable with the bulk spin operators [1206.0814]. In the half-infinite chain, the thermodynamic-limit Hamiltonian is written with anisotropy \(\Delta=(q+q^{-1})/2\) and boundary parameters \(\epsilon_\pm,k_\pm\), and the analysis distinguishes generic non-diagonal, generic diagonal, and one-sided upper or lower non-diagonal boundaries [1211.6304]. ## 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\), one introduces the bulk monodromy matrix \(M_0(\lambda)\), the reflected monodromy, and boundary monodromies such as \[ U_-(\lambda)=M_0(\lambda)K_-(\lambda)\widehat M_0(\lambda), \] with commuting transfer matrix \[ T(\lambda)=\operatorname{tr}_0\!\left[K_+(\lambda)\,U_-(\lambda)\right]. \] The Hamiltonian is recovered from the transfer matrix by a logarithmic derivative at a special point in the homogeneous limit [1307.3960]. In the most general scalar case, the boundary \(K\)-matrix is \[ K(\lambda;\xi,\kappa,\tau)= \begin{pmatrix} \sinh(\lambda+\xi) & \kappa e^\tau \sinh(2\lambda)\\[1mm] \kappa e^{-\tau}\sinh(2\lambda) & \sinh(\xi-\lambda) \end{pmatrix}, \] which is the origin of non-diagonal boundary terms [1307.3960]. For higher spin, the fused transfer matrices satisfy a fusion hierarchy. At root-of-unity anisotropy \[ \eta = \frac{i\pi}{p+1}, \qquad p=1,2,3,\dots, \] the hierarchy truncates and yields a finite functional equation of order \(p+1\) for the fundamental transfer-matrix eigenvalues. This truncation underlies the generalized \(T\)-\(Q\) analysis of the open spin-\(s\) chain with non-diagonal boundaries [1206.0814]. In the half-infinite setting, the transfer matrix can be rewritten in terms of a current algebra \(\mathcal{O}_q(\widehat{\mathfrak{sl}_2})\) with currents \[ W_\pm(u),\qquad Z_\pm(u), \] 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 [1211.6304]. A notable open-chain distinction is that only the even charges \[ Q_{2r},\qquad r=1,2,\dots \] 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 [1308.1584]. ## 3. Bethe ansatz, generalized \(T\)-\(Q\) 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 \(M\) 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 \(k_j\mapsto -k_j\). Consequently, for fixed \((x_0,\ldots,x_{M-1})\) the open-chain Bethe wave function contains \[ 2^{M} M! \] terms, compared with only \(M!\) terms for the periodic chain. The Bethe roots obey \[ \frac{\alpha(k_{j})\, \beta(k_{j})}{\alpha(-k_{j})\, \beta(-k_{j})} = \prod_{l=0; l\ne j}^{M-1} \frac{B(-k_{j}, k_{l})}{B(k_{j}, k_{l})}, \] and the energy is \[ E(\{k_j\}) = -\tfrac{1}{2}\left[(L-1)\Delta + h + h' \right] + 2\sum_{j=0}^{M-1} (\Delta-\cos(k_{j})) . \] In the diagonal-boundary massive regime, the boundary algebraic Bethe ansatz uses states \[ \ket{\{\lambda\}}=\prod_{j=1}^N \mathcal B(\lambda_j)\ket{0}, \] with energy \[ E(\{\lambda\}) = h^++h^- -\sum_{j=1}^N \frac{4\sinh^2\zeta}{\cosh\zeta-\cos(2\lambda_j)} . \] These formulations are standard starting points for spectral and overlap problems [2109.05607][2409.15194]. 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\) XXZ chain with the most general integrable boundary terms, the transfer-matrix spectrum is simple and every eigenvalue is characterized by a system of \(N\) quadratic equations. The eigenstates are reconstructed in an SOV basis from discrete Baxter-like relations, and scalar products of separate states reduce to determinants [1307.3960][1206.0646]. At root-of-unity anisotropy, the non-diagonal open spin-\(s\) chain admits generalized \(T\)-\(Q\) relations involving two independent \(Q\)-functions, \[ \Lambda^{(1/2,s)}(u) = h^{(1)}(u)\frac{Q_1(u-\eta)}{Q_2(u)} + h^{(2)}(u)\frac{Q_2(u+\eta)}{Q_1(u)} , \] with \[ Q_j(u)=\prod_{k=1}^{M_j} \sinh(u-u^{(j)}_k)\,\sinh(u+u^{(j)}_k+\eta), \qquad j=1,2. \] The construction works for special boundary families with at most two arbitrary boundary parameters, and numerical evidence for \(s=1\) shows completeness for small chains [1206.0814]. A related arbitrary-spin analysis at \[ \eta = i\pi \frac{r}{q}, \] with \(q\) restricted to odd integers, derives \(T\)-\(Q\) relations and Bethe equations for selected non-diagonal boundary choices and reproduces the complete spectra for representative spin-\(\tfrac12\) and spin-1 chains [1403.4313]. A different constrained regime is governed by the phantom roots criterion, \[ (N-2M-1)\eta = \theta_-+\theta_+ +\alpha_- - \alpha_+ +\beta_- - \beta_+ \quad \mod 2\pi, \] under which the Hilbert space splits into two invariant subspaces \(\mathcal G_M^+\) and \(\mathcal G_M^-\), 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 [2107.13266]. At the combinatorial point \(q=e^{2\pi i/3}\), a boundary qKZ solution \(|\Psi_N(z_1,\dots,z_N)\rangle\) yields an additional exact route. The generalized overlap \(\Omega_{N_1,N_2}\) 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 [2111.15223]. ## 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 \(F\)-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 [1011.4719]. 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 [1807.05197]. 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 [1307.3960][1206.0646]. 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 [1412.5389]. 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, \[ S(\{\lambda\},\{\mu\})=O(L^{-\infty}), \] whereas in the real-root and compatible boundary-root cases it reduces to a closed expression in double \(q\)-Pochhammer symbols, up to exponentially small corrections [2409.15194]. At \(\Delta=-\tfrac12\), the logarithmic bipartite fidelity \[ \mathcal F_{N_1,N_2} = -\ln \left(\frac{\left|\langle \psi_{N}|\left(|\psi_{N_1}\rangle\otimes |\psi_{N_2}\rangle\right)\right|^2}{\|\psi_{N}\|^2 \|\psi_{N_1}\|^2\|\psi_{N_2}\|^2}\right) \] is obtained exactly at finite size from a determinant formula for the overlap \(\mathcal O_{N_1,N_2}\). Its large-\(N\) asymptotics contains the universal term \[ \mathcal F_{N_1,N_2}\sim \frac16\ln N + \cdots, \] and the calculation confirms the conformal-field-theory prediction with \(c=1\) and \(g(\xi)=0\) for the \(N^{-1}\ln N\) correction [2111.15223]. ## 5. Boundary thermodynamics, surface terms, and boundary-localized modes The thermodynamics of the open XXZ spin-\(\tfrac12\) chain differs from the periodic problem by an \(O(1)\) boundary contribution. For system size \(M\), \[ F(M,T) = M f_{\rm bulk}(T) + f_{\rm surf}(T) + o(1), \] and \(f_{\rm surf}\) 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 \(f_{\rm surf}\) is written in terms of the same nonlinear integral equation that governs the periodic XXZ chain, and the boundary magnetization follows from \[ \langle \sigma_1^z \rangle_T = - \frac{\sinh^2(\xi_-)}{J\sinh(\eta)} \, \frac{\partial f_{surf}}{\partial \xi_-}. \] In the massless regime the leading thermal correction to the boundary magnetization is \(T^2\), not \(T\) [1201.5884]. At zero temperature in the massive antiferromagnetic regime, the ground state may contain a boundary-localized Bethe excitation, the boundary root, \[ \alpha_{\text{BR}}^\sigma = -i\left(\zeta/2+\xi_\sigma+\epsilon_\sigma\right), \qquad \epsilon_\sigma=O(L^{-\infty}) . \] For even chain length and equal boundary fields satisfying \[ h_+=h_-=h,\qquad |h|