The Bethe Ansatz

The planar Heisenberg chain is obtained from the unified treatment by taking

\[ \varphi (\lambda) = \sinh \lambda, \hspace{5mm} \eta = -i \zeta \]

with \(0 < \zeta \leq \pi\).

R-matrix

\begin{equation*} R (\lambda) = \left( \begin{array}{cccc} 1 & 0 & 0 & 0 \\ 0 & b (\lambda) & c (\lambda) & 0 \\ 0 & c (\lambda) & b (\lambda) & 0 \\ 0 & 0 & 0 & 1 \end{array} \right) \end{equation*} \begin{align*} a(\lambda) &= 1, & b(\lambda) &= \frac{\sinh(\lambda)}{\sinh(\lambda - i\zeta)}, \\ c(\lambda) &= \frac{-i\sin(\zeta)}{\sinh(\lambda -i\zeta)}, & d(\lambda) &= b(\lambda +i\zeta/2)^N = \left(\frac{\sinh (\lambda +i\zeta/2)}{\sinh(\lambda -i \zeta/2)}\right)^N \end{align*}

Transfer matrix

\begin{equation*} L_j (\lambda) = R_{aj} (\lambda +\zeta/2) \end{equation*}

Transfer matrix eigenvalue

\begin{align*} \tau (\lambda | \{ \lambda_j \}) &= \frac{\sinh(\bar{\lambda} - \lambda -i\zeta)}{\sinh (\bar{\lambda} - \lambda)} + \left[ \frac{\sinh(\lambda +i\zeta/2)}{\sinh (\lambda -i\zeta/2)} \right]^N \frac{\sinh(\bar{\lambda} - \lambda +i\zeta)}{\sinh (\bar{\lambda} - \lambda)} \end{align*}

Hamiltonian

\begin{equation*} H_{XXZ} = \frac{\sin \zeta}{2i} \frac{d}{d\lambda} \ln \tau (\lambda) |_{\lambda = -i\zeta/2} = \sum_{j=1}^N \left[S^x_j S^x_{j+1} + S^y_j S^y_{j+1} + \cos \zeta ~(S^z_j S^z_{j+1} - \frac{1}{4}) \right] \end{equation*}

Bethe equations

\[ \frac{a(\lambda_j)}{d(\lambda_j)} \frac{b(\lambda_j, \bar{\lambda}_j)}{b(\bar{\lambda}_j, \lambda_j)} = \left(\frac{\sinh (\lambda_j -i\zeta/2)}{\sinh (\lambda_j +i\zeta/2)}\right)^N \frac{\sinh (\lambda_j - \bar{\lambda}_j +i\zeta)}{\sinh(\lambda_j - \bar{\lambda}_j -i\zeta)} = 1 \]

or equivalently (using \(\sinh (-\lambda) = -\sinh (\lambda)\) and the fact that the term \(j\) is now included in the product)

\begin{align*} &r (\lambda_j) + s(\lambda_j | \{ \lambda \}) = 0, \\ &r (\lambda_j) \equiv \frac{a(\lambda_j)}{d(\lambda_j)} = \left(\frac{\sinh (\lambda_j -i\zeta/2)}{\sinh (\lambda_j +i \zeta/2)}\right)^N, \\ &s (\lambda_j | \{ \lambda \}) \equiv \frac{b(\bar{\lambda}, \lambda_j)}{b(\lambda_j, \bar{\lambda})} = \frac{\sinh (\lambda_j - \bar{\lambda} -i\zeta)}{\sinh(\lambda_j - \bar{\lambda} +i\zeta)} \end{align*} \begin{align*} &{\cal E} (\lambda_j | \{ \lambda \}) + 1 = 0, \\ &{\cal E} (\lambda | \{ \lambda \}) \equiv \frac{a(\lambda_j)}{d(\lambda_j)} \frac{b(\lambda_j, \bar{\lambda})}{b(\bar{\lambda}, \lambda_j)} = \left(\frac{\sinh (\lambda_j -i\eta/2)}{\sinh (\lambda_j +i\zeta/2)}\right)^N \frac{\sinh (\lambda_j - \bar{\lambda} +i\zeta)}{\sinh(\lambda_j - \bar{\lambda} -i\zeta)} \end{align*}

Energy

\[ E(\{ \lambda \}) = \sum_j \frac{-\sin^2 \zeta}{\cosh 2\lambda_j - \cos \zeta} \]

Momentum

\[ e^{i k(\lambda_j)} = \frac{\sinh(\lambda_j +i\zeta/2)}{\sinh(\lambda_j - i\zeta/2)} \]

\[ P(\{ \lambda \}) = i \ln \tau (-i\zeta/2 | \{ \lambda \}) = \frac{i}{N} \ln \bar{r} = i \ln \frac{\sinh(\bar{\lambda} -i\zeta/2)}{\sinh(\bar{\lambda} +i\zeta/2)} \]

\[ e^{-i \hat{P}} {\cal O}_j e^{i \hat{P}} = {\cal O}_{j+1} \]

Scalar products

\[ S_M (\{ \mu \}, \{ \lambda \}) \equiv \langle 0 | C(\bar{\mu}) B(\bar{\lambda}) | 0 \rangle \]

\begin{align*} &S_M (\{ \mu \}, \{ \lambda \}) = \frac{\mbox{det}~T(\{ \mu \}, \{ \lambda \})}{\mbox{det}~V(\{ \mu \}, \{ \lambda \})}, \\ &T_{jk} = \frac{\partial}{\partial \lambda_k} \tau (\mu_j | \{ \lambda \}), \hspace{5mm} V_{jk} = \frac{1}{\sinh (\mu_k - \lambda_j)} \end{align*}

\[ \mathbb{B}(\lambda) \equiv \frac{B(\lambda)}{d(\lambda)}, \hspace{5mm} \mathbb{C}(\lambda) \equiv \frac{C(\lambda)}{d(\lambda)} \]

Gaudin matrix and state norm

\begin{equation*} {\mathbb N}_M = \langle 0 | C (\bar{\lambda}) B (\bar{\lambda}) | 0 \rangle = \sinh^M(\eta) \prod_{j \neq k} \frac{\sinh(\lambda_j - \lambda_k + \eta)}{\sinh(\lambda_j - \lambda_k)} \det {\cal G} (\{ \lambda \}) \end{equation*} \begin{align*} {\cal G}_{jk} (\{ \lambda \}) &= -\frac{\partial}{\partial \lambda_k} \ln \left[ \frac{a(\lambda_j)}{d(\lambda_j)} \frac{b(\lambda_j, \bar{\lambda}_j)}{b(\bar{\lambda}_j, \lambda_j)}\right] = -\frac{\partial}{\partial \lambda_k} \ln \bigl(-{\cal E} (\lambda_j | \{ \lambda \})\bigr) \\ &= -\frac{\partial}{\partial \lambda_k} \ln r(\lambda_j) + \frac{\partial}{\partial \lambda_k} \ln s (\lambda_j | \{ \lambda \}) \\ &= -\delta_{jk} ~N ~\partial_\lambda \ln \left.\frac{\sinh(\lambda + \eta/2)}{\sinh(\lambda - \eta/2)} \right|_{\lambda_j} + \partial_{\lambda_k} \ln \frac{\sinh(\lambda_j - \bar{\lambda} + \eta)}{\sinh(\lambda_j - \bar{\lambda} - \eta)} \end{align*}

Matrix element for \(\sigma^-\)

\begin{equation*} F_j^- (\{ \mu \}_{M+1}, \{ \lambda \}_M) = \langle 0 | C(\bar{\mu}) ~\sigma^-_j ~B(\bar{\lambda}) | 0 \rangle \end{equation*} \begin{align*} F_j^- (\{ \mu \}_{M+1}, \{ \lambda \}_M) &= \frac{\phi_{j-1} (\{ \mu \})}{\phi_{j-1} (\{ \lambda \})} \frac{\sinh(\bar{\mu} + \eta/2)}{\sinh(\bar{\lambda} + \eta/2)} \frac{\det_{M+1} H^- (\{ \mu \}, \{ \lambda \})}{\prod_{M+1 \geq l > m \geq 1} \sinh(\mu_l - \mu_m) \prod_{1 \leq l < m \leq M} \sinh(\lambda_l - \lambda_m)} \end{align*}

\[ \phi_j (\{ \lambda \}) = \left(b(\bar{\lambda} - \eta/2)\right)^{-j} = e^{iP(\{ \lambda \}) j} \]

\begin{align*} \left. H^{-}_{jk} \right|_{k < M+1} &= \frac{\sinh(\eta)}{\sinh(\mu_j - \lambda_k)} \left( a(\lambda_k) \sinh(\bar{\mu}_j - \lambda_k + \eta) - d(\lambda_k) \sinh (\bar{\mu}_j - \lambda_k -\eta) \right), \nonumber \\ H^{-}_{j M+1} &= \frac{\sinh(\eta)}{\sinh(\mu_j + \eta/2) \sinh(\mu_j - \eta/2)}. \end{align*}

Matrix element for \(\sigma^z\)

\begin{equation*} F_j^z (\{ \mu \}_{M}, \{ \lambda \}_M) = \langle 0 | C(\bar{\mu}) ~\sigma^z_j ~B(\bar{\lambda}) | 0 \rangle \end{equation*} \begin{align*} F_j^z (\{ \mu \}_{M}, \{ \lambda \}_M) = \frac{\phi_{j-1} (\{ \mu \})}{\phi_{j-1} (\{ \lambda \})} \frac{\sinh(\bar{\mu} + \eta/2)}{\sinh(\bar{\lambda} + \eta/2)} \frac{\det_{M} \bigl(H (\{ \mu \}, \{ \lambda \}) - 2P (\{ \mu \}, \{ \lambda \}) \bigr)}{\prod_{M \geq l > m \geq 1} \sinh(\mu_l - \mu_m) \prod_{1 \leq l < m \leq M} \sinh(\lambda_l - \lambda_m)} \end{align*} \begin{align*} H_{jk} &= \frac{\sinh(\eta)}{\sinh(\mu_j - \lambda_k)} \left( a(\lambda_k) \sinh(\bar{\mu}_j - \lambda_k + \eta) - d(\lambda_k) \sinh (\bar{\mu}_j - \lambda_k -\eta) \right), \\ P_{jk} &= \sinh(\eta) \frac{\sinh(\bar{\lambda} - \lambda_k + \eta)}{\sinh(\mu_j - \eta/2) \sinh(\mu_j + \eta/2)} \end{align*}