---
title: 'Trigonometric Gaudin Models: Integrability Insights'
url: https://www.emergentmind.com/topics/trigonometric-gaudin-models
type: topic
---

# Trigonometric Gaudin Models: Integrability Insights

Searching arXiv for recent and foundational papers on trigonometric Gaudin models.
Trigonometric Gaudin models are integrable many-body systems whose interaction kernels are governed by trigonometric or hyperbolic dependence on the spectral or inhomogeneity parameters, rather than the rational dependence of the XXX-type Gaudin model. In the standard quantum-inverse-scattering perspective, they arise as quasi-classical limits of XXZ-type spin chains, while in classical form they are encoded by Lax matrices satisfying linear or boundary-modified \(r\)-matrix brackets with trigonometric \(r\)-matrices. Across the literature, the term includes bulk \(sl(2)\) and \(\mathfrak{gl}_N\) Gaudin magnets, open-boundary models obtained from reflection equations, higher-rank Richardson–Gaudin systems, quasi-trigonometric non-skew-symmetric deformations, and more geometric or dual formulations parameterized by moduli spaces or related to Calogero–Sutherland and Ruijsenaars–Schneider systems [1303.2481].

## 1. Algebraic definition and basic structures

The standard trigonometric Gaudin framework is characterized by a trigonometric classical \(r\)-matrix and a Lax matrix whose spectral invariants generate commuting Hamiltonians. For the classical \(sl(2)\) trigonometric Gaudin magnet, the Lax matrix can be written as
\[
L(\lambda)=
\begin{pmatrix}
A(\lambda) & B(\lambda)\\
C(\lambda) & -A(\lambda)
\end{pmatrix},
\]
with
\[
A(\lambda)=\sum_{j=1}^{N}\cot(\lambda-\lambda_j)\,s_j^3,\qquad
B(\lambda)=\sum_{j=1}^{N}\frac{s_j^-}{\sin(\lambda-\lambda_j)},\qquad
C(\lambda)=\sum_{j=1}^{N}\frac{s_j^+}{\sin(\lambda-\lambda_j)},
\]
where the site variables satisfy the direct-sum \(sl(2)\) Poisson algebra and the Casimirs \((s_j^3)^2+s_j^+s_j^- = s_j^2\) fix the dynamics on a product of two-spheres [1003.0389]. The corresponding trigonometric \(r\)-matrix is
\[
r_t(\lambda)=\frac{i}{\sin\lambda}
\begin{pmatrix}
\cos\lambda & 0 & 0 & 0\\
0 & 0 & 1 & 0\\
0 & 1 & 0 & 0\\
0 & 0 & 0 & \cos\lambda
\end{pmatrix},
\]
and the Lax matrix obeys the linear bracket
\[
\{L(\lambda)\otimes \mathbf 1,\mathbf 1\otimes L(\mu)\}
=
[r_t(\lambda-\mu),L(\lambda)\otimes\mathbf 1+\mathbf 1\otimes L(\mu)].
\]
The generating function of commuting integrals is
\[
-\det L(\lambda)=A^2(\lambda)+B(\lambda)C(\lambda),
\]
whose residues yield the Gaudin Hamiltonians [1003.0389].

For \(\mathfrak{gl}_N\), the trigonometric model can be encoded in the current algebra \(\mathfrak g^+=\mathfrak b^+\oplus t^{-1}\mathfrak{gl}_N[t^{-1}]\) with current matrix \(\mathcal L^+(u)\) satisfying
\[
[\mathcal L^+_1(u),\mathcal L^+_2(v)]
=
[\mathcal L^+_1(u)+\mathcal L^+_2(v),r_{12}(u/v)],
\]
where
\[
r(x)=\sum_{i,j=1}^N\left(\frac{1+x}{1-x}+\operatorname{sgn}(j-i)\right)e_{ij}\otimes e_{ji}.
\]
Under tensor products of vector evaluation representations with parameters \(a_1,\dots,a_\ell\), one has
\[
\mathcal L^+(u)\mapsto r_{01}(u/a_1)+\cdots+r_{0\ell}(u/a_\ell),
\]
and the residues of \(\operatorname{tr}\mathcal L^+(u)^2\) reproduce the standard quadratic trigonometric Gaudin Hamiltonians [1802.06499].

A recurrent structural theme is that trigonometric models are multiplicative rather than additive in spectral dependence. This appears in the \(\mathfrak{gl}_N\) current relation through \(u/v\), in multiplicative trigonometric \(r\)-matrices, and in geometric parameter spaces such as \(M_{n+2}=(\mathbb C^\times)^n\setminus \Delta/\mathbb C^\times\) for trigonometric Gaudin subalgebras [2407.06424]. This suggests a precise sense in which trigonometric Gaudin theory interpolates between additive rational-current formulations and multiplicative XXZ-type constructions.

## 2. Classical and quantum integrability

A standard conceptual statement is that rational Gaudin models are quasi-classical limits of XXX-type chains, whereas trigonometric Gaudin models are quasi-classical limits of XXZ-type chains [1303.2481]. In the periodic case, evaluating the transfer matrix at inhomogeneities and expanding in the quasi-classical parameter produces commuting Gaudin Hamiltonians. For the open XXZ chain, the same mechanism is implemented through the Sklyanin double-row transfer matrix
\[
t(\lambda)=\operatorname{tr}_0\bigl(K^+_0(\lambda)\,\mathcal T_0(\lambda)\bigr),
\]
whose quasi-classical expansion at \(\lambda=\alpha_k\) has the form
\[
t(\alpha_k)=\mathbb 1+\eta H_k+O(\eta^2),
\]
yielding commuting open-boundary trigonometric Gaudin Hamiltonians [1303.2481].

In the boundary-free \(sl(2)\) case, the Hamiltonians extracted from
\[
-\det L(\lambda)
=
\sum_{i=1}^N\left(\frac{s_i^2}{\sin^2(\lambda-\lambda_i)}+H_i\cot(\lambda-\lambda_i)\right)-H_0^2
\]
are
\[
H_i=
\sum_{k\neq i}
\frac{
2\cos(\lambda_i-\lambda_k)s_i^3 s_k^3+s_i^+s_k^-+s_i^-s_k^+
}{
\sin(\lambda_i-\lambda_k)
},
\qquad
H_0=\sum_j s_j^3,
\]
with \(\{H_i,H_j\}=0\) and \(\sum_i H_i=0\) [1003.0389]. The model also admits a rationalized description via \(z=e^{i\lambda}\), under which
\[
L(z)=\sigma_3L(-z)\sigma_3,
\]
so the trigonometric Gaudin magnet becomes a rational \(2N\)-pole object with a \(\mathbb Z_2\) reflection constraint [1003.0389].

In the quantum \(\mathfrak{gl}_N\) setting, higher commuting Hamiltonians arise from the \(q\to 1\) limit of XXZ Bethe subalgebra elements in \(Y_q(\mathfrak{gl}_N)\). The generating series
\[
\sum_{m=1}^\infty \theta_m\,y^m
=
\sum_{s=1}^\infty y^s\, \operatorname{tr}_{1,\dots,s}
\Bigl(
\mathcal T_{s-1,s}(y)\cdots \mathcal T_{12}(y)\,
\mathcal L_1\cdots \mathcal L_s
\Bigr),
\qquad
\mathcal L=2u\partial_u-\mathcal L^+(u),
\]
defines explicit higher Hamiltonians whose coefficients generate a commutative subalgebra of \(U(\mathfrak g^+)\) [1802.06499]. This is the trigonometric analogue of Talalaev-type higher Gaudin Hamiltonians and shows that the trigonometric theory is not restricted to quadratic Hamiltonians.

## 3. Boundary terms and reflection-algebra formulations

Open trigonometric Gaudin models are constructed through reflection equations. For the non-symmetric XXZ \(R\)-matrix, the left boundary matrix \(K^-(\lambda)\) satisfies
\[
R_{12}(\lambda-\mu)\,K^-_1(\lambda)\,R_{21}(\lambda+\mu)\,K^-_2(\mu)
=
K^-_2(\mu)\,R_{12}(\lambda+\mu)\,K^-_1(\lambda)\,R_{21}(\lambda-\mu),
\]
while the right boundary matrix \(K^+(\lambda)\) is related by
\[
K^+(\lambda)=K^-(-\lambda-\eta)\,M.
\]
With suitable normalization, the quasi-classical expansion of the open transfer matrix yields
\[
H_k=
\Gamma_k(\alpha_k)
+\sum_{l\neq k} r_{kl}(\alpha_k-\alpha_l)
+\sum_{l\neq k} K^-_k(\alpha_k)\,r_{lk}(\alpha_k+\alpha_l)\,K^-_k(-\alpha_k),
\]
where the three terms are, respectively, a local boundary term, the usual bulk interaction, and a reflected interaction depending on \(\alpha_k+\alpha_l\) [1303.2481].

This reflected term is the hallmark of open trigonometric Gaudin systems. Compared with the bulk Hamiltonian, open boundaries add “image-spin” interactions dressed by the boundary matrix. The same paper also shows that a classical-reflection-equation construction based on
\[
\mathcal L(\lambda)=L(\lambda)-K(\lambda)L(-\lambda)K(-\lambda)
\]
does not reproduce the same Hamiltonians as the quasi-classical limit of the open XXZ transfer matrix, making the notion of “open trigonometric Gaudin model” construction-dependent [1303.2481].

For triangular boundary conditions, the trigonometric \(sl(2)\) model admits a full algebraic Bethe ansatz. The left and right reflection matrices are specialized to upper-triangular form, and in the Gaudin limit the condition
\[
\lim_{\eta\to0}\bigl(K^+(\lambda)K^-(\lambda)\bigr)
=
\kappa^2 \sinh(\xi-\lambda)\sinh(\xi+\lambda)\,\mathbf 1
\]
ties the two boundaries together [1705.02235]. The corresponding generating function is
\[
\tau(\lambda)=\operatorname{tr}_0\mathcal L_0^2(\lambda),
\qquad
\mathcal L_0(\lambda)=L_0(\lambda)-K_0(\lambda)L_0(-\lambda)K_0^{-1}(\lambda),
\]
and the residues at \(\lambda=\pm \alpha_m\) give the trigonometric Gaudin Hamiltonians with boundary terms [1709.06419]. Appropriate Bethe vectors yield a simple off-shell action of \(\tau(\lambda)\), with Bethe equations
\[
2\rho(\mu_i)
+\frac{1}{\sinh(\xi+\mu_i)\sinh(\xi-\mu_i)}
-
\sum_{j\neq i}^M
\frac{2}{\sinh(\mu_i+\mu_j)\sinh(\mu_i-\mu_j)}
=0,
\]
which are the boundary trigonometric analogue of the rational Gaudin equations [1709.06419].

A broader generalization is the classical \(N\)-reflection equation. Given a solution \(k(\nu)\) of that equation, one defines
\[
\bar r_{ab}(\lambda,\nu)
=
\sum_{j=0}^{N-1}
g^{(j)}(\nu)\,
k_b^{(j)}(\nu)\,
r_{ab}\bigl(\lambda,\tau^j(\nu)\bigr)\,
\bigl(k_b^{(j)}(\nu)\bigr)^{-1},
\]
which again satisfies the classical Yang–Baxter equation, though it is generally non-skew-symmetric [1803.09931]. This produces new boundary, folded, or cyclotomic Gaudin-type models and includes explicit trigonometric \(N=2\) and \(N=3\) solutions with Möbius transformations of the spectral parameter. The paper does not write explicit trigonometric Gaudin Hamiltonians in full detail, but it gives the exact framework for constructing them through residues of \(\operatorname{tr}B(\lambda)^2\) [1803.09931].

## 4. Separation of variables, higher structures, and deformations

A notable extension of the standard trigonometric framework is the quasi-trigonometric \(\mathfrak{gl}(2)\) Gaudin model governed by a non-skew-symmetric, non-dynamical classical \(r\)-matrix
\[
r(u,v)=\left(\frac{1}{2}\frac{u+v}{v-u}+c_1\right)X_{11}\otimes X_{11}
+\left(\frac{1}{2}\frac{u+v}{v-u}+c_2\right)X_{22}\otimes X_{22}
+\frac{v}{v-u}X_{12}\otimes X_{21}
+\frac{u}{v-u}X_{21}\otimes X_{12},
\]
with \(c_2=-c_1\) the distinguished one-parameter subfamily supporting two full separation-of-variables constructions [2107.08376]. At \(c_1=c_2=0\) one recovers the usual skew-symmetric trigonometric \(r\)-matrix, but for other values the model is genuinely non-skew-symmetric and satisfies the generalized classical Yang–Baxter equation [2107.08376].

For \(N\) sites, the quasi-trigonometric Gaudin Lax matrix is
\[
L(v)=\sum_{i=1}^N(\cdots)+c(v),
\qquad
c(v)=2c_{11}X_{11}+2c_{22}X_{22},
\]
and the commuting Hamiltonians are the residues
\[
H_j=\frac{1}{2}\operatorname*{res}_{v=\nu_j}\operatorname{tr}L^2(v).
\]
Relative to the standard trigonometric Hamiltonians, the deformation adds the non-skew-symmetric \(r\)-matrix parameters \(c_1,c_2\), on-site quadratic terms, and external magnetic-field terms \(2c_{11}S_{11}^{(j)}+2c_{22}S_{22}^{(j)}\) [2107.08376].

The same model possesses an additional geometric integral
\[
M=\frac12\sum_{j=1}^N\bigl(S_{22}^{(j)}-S_{11}^{(j)}\bigr),
\]
central in the separation-of-variables construction. For \(c_2=-c_1\), two one-parameter families of separating functions \((A(u),B(u))\) generate complete canonical separated variables when \(k\neq 0\), and for generic parameters the separation curves are not the ordinary spectral curve
\[
\det(L(u)-\lambda\,{\rm Id})=0,
\]
but rather \(M\)-shifted spectral curves [2107.08376]. This is one of the clearest demonstrations that trigonometric Gaudin theory naturally extends beyond standard Sklyanin-type separated coordinates.

The classical trigonometric Gaudin magnet also admits explicit Bäcklund transformations, interpreted as exact integrable time-discretizations. The map is defined by
\[
\tilde L(z)=D(z)L(z)D^{-1}(z),
\]
with a Darboux matrix preserving the reflection symmetry \(L(z)=\sigma_3L(-z)\sigma_3\). In trigonometric variables the Darboux matrix takes the XXZ-type form
\[
D(\lambda)=\alpha
\begin{pmatrix}
\sin(\lambda-\lambda_0-\mu)+PQ\,\cos(\lambda-\lambda_0) & P\,\cos\mu\\
Q\,\sin(2\mu)-PQ^2\cos\mu & \sin(\lambda-\lambda_0+\mu)-PQ\,\cos(\lambda-\lambda_0)
\end{pmatrix},
\]
and the interpolating Hamiltonian is
\[
\mathcal H(\lambda_0)=\gamma(\lambda_0)=\sqrt{A^2(\lambda_0)+B(\lambda_0)C(\lambda_0)}.
\]
The rational Gaudin Bäcklund transformation is recovered by a small-angle limit [1003.0389].

A more indirect structural development is the bi-Hamiltonian approach to Gaudin algebras via compatible Lie brackets on quotient current algebras \( \mathfrak q[t]/(p)\). This framework explicitly captures rational Gaudin algebras and certain limits, but it does not explicitly construct standard trigonometric Gaudin Hamiltonians; its relevance to trigonometric models is therefore structural and suggestive rather than direct [2105.01020].

## 5. Dualities, moduli spaces, and generalized parameter spaces

The trigonometric quantum-classical duality relates the twisted inhomogeneous XXZ chain to the classical trigonometric Ruijsenaars–Schneider model, and in its quasi-classical/nonrelativistic limit it produces a duality between the quantum trigonometric Gaudin model and the classical Calogero–Sutherland system [1510.07509]. On the quantum side, the XXZ transfer matrix
\[
\hat T^{\rm XXZ}(z)=\operatorname{tr}_0\, V_0\, R_{01}(z-q_1)\cdots R_{0N}(z-q_N)
\]
is expanded with
\[
V^\varepsilon = 1_{n\times n}+\varepsilon\,\operatorname{diag}(v_1,\dots,v_n)+O(\varepsilon^2),
\]
and the trigonometric Gaudin Hamiltonians appear as residues
\[
\hat H_i^{\rm G} = \mathop{\rm Res}_{z=q_i}\hat T_2(z;\{q_i\}),
\]
with explicit formula
\[
\hat H_i^{\rm G}
=
\sum_{a=1}^n v_a e_{aa}^{(i)}
+
\sum_{j\neq i}
\left(
\frac{h}{\sinh(q_i-q_j)}\sum_{a\neq b} e_{ab}^{(i)}e_{ba}^{(j)}
+
h\,\coth(q_i-q_j)\sum_{a=1}^n e_{aa}^{(i)}e_{aa}^{(j)}
\right)
\]
[1510.07509]. The dual classical spectrum forms strings centered at the twist parameters \(v_a\), a trigonometric splitting phenomenon absent in the rational case.

A different geometric picture is provided by the family of trigonometric Gaudin subalgebras
\[
A^{\mathrm{trig}_\theta(z_1,\dots,z_n)\subset (U\mathfrak g^{\otimes n})^{\mathfrak h},
\]
defined for arbitrary complex simple \(\mathfrak g\) as images of homogeneous Gaudin algebras under a quantum Hamiltonian reduction map \(\psi_\theta\) [2407.06424]. Their quadratic Hamiltonians are
\[
H_{i,\theta}^{\mathrm{trig}}
=
\sum_{j\ne i}\frac{\Omega^{(ij)}}{z_i-z_j}
+\frac{\theta^{(i)}}{z_i}
-\sum_j\frac{\Omega_-^{(ij)}}{z_i},
\]
and the family is parameterized by \(M_{n+2}\), compactifying over \(\overline M_{n+2}\) [2407.06424]. A central degeneration formula is
\[
\lim_{\varepsilon\to 0}
A^{\mathrm{trig}_{\varepsilon^{-1}\chi}(1-\varepsilon z_1,\dots,1-\varepsilon z_n)
=
A_\chi(z_1,\dots,z_n),
\]
showing that rational inhomogeneous Gaudin subalgebras arise as limits of trigonometric ones [2407.06424]. The same work further identifies real loci with simple-spectrum regimes and interprets the resulting monodromy as an action of affine cactus-type groups.

For relativistic \({\rm GL}(NM)\) Gaudin-type systems, the elliptic theory admits a universal \(R\)-matrix formulation
\[
L^{ij}(z)=\sum_{a=1}^{n}\operatorname{tr}_2\!\Big(R^{\,z-z_a}_{12}(q_{ij}+\eta)\,P_{12}\,S^{ij,a}_2\Big),
\]
and the authors state explicitly that this description also provides trigonometric and rational analogues whenever the chosen \(R\)-matrix satisfies the associative Yang–Baxter equation and the required additional properties [2204.06137]. The paper does not work out a separate trigonometric function-theoretic model in full detail, but it places trigonometric relativistic Gaudin analogues inside a uniform \(R\)-matrix classification.

Some duality results with irregular singularities are likewise structurally relevant. They concern rational Gaudin algebras with Takiff data and do not directly construct trigonometric Hamiltonians, but they supply determinant/Berezinian duality mechanisms, oscillator realizations, and classical spectral identities that are closely comparable to the kinds of spectral correspondences appearing in trigonometric Gaudin theory [2507.00730]. This suggests a broader landscape in which trigonometric, cyclotomic, and irregular-rational models should be viewed as adjacent rather than isolated constructions.

## 6. Physical realizations and dynamical phenomena

Trigonometric Richardson–Gaudin models appear in a concrete open-system realization: a system of \(N\)-level atoms coupled to a Markovian reservoir can be mapped, after vectorization, to a Liouvillian that is a linear combination of commuting trigonometric \(SU(N)\) Richardson–Gaudin charges [2005.05155]. The trigonometric kernels are
\[
X_{mm'}=\frac{e^{i(z_{m'}-z_m)}}{\sin(z_{m'}-z_m)},
\qquad
Z_{m'm}=\cot(z_{m'}-z_m),
\]
and in the two-copy realization the Liouvillian spectrum is determined by nested Bethe equations of trigonometric Gaudin type [2005.05155]. For \(SU(2)\) this recovers the exact solution of the collective dissipative spin model, while for \(SU(3)\) the paper studies the steady state and dissipative gaps for finite systems and in the thermodynamic limit [2005.05155]. This places trigonometric Gaudin theory directly inside non-Hermitian and dissipative quantum dynamics.

On the classical side, the \(\mathrm{su}(2)\) rational and trigonometric Gaudin models can be viewed as generalized coupled angular momenta systems on
\[
\mathcal M=\mathbb S^2\times \mathbb S^2,
\]
with momentum map
\[
F=(J,H_{w,\mathbf t}),
\qquad
J=R_1z_1+R_2z_2.
\]
The trigonometric subfamily is \(w=0\), with Hamiltonian
\[
H_{0,\mathbf t}=t_0(z_1+z_2)^2+t_3(x_1x_2+y_1y_2)+t_4z_1z_2
\]
[2410.01372]. The rank-zero singularities \(m_0\) and \(m_2\) undergo Hamiltonian Hopf bifurcations at
\[
t_{4,m_0^\pm}=t_{4,m_2^\pm}
=
\pm \frac{2t_3\sqrt{R_1R_2}}{R_1+R_2},
\]
and the normal form computed up to sixth order determines when the bifurcation is non-degenerate, supercritical, or subcritical [2410.01372]. In the subcritical case a flap appears in the image of the momentum map, and the figures in that work show additional global bifurcations involving pleats and cusp collisions [2410.01372]. This suggests that even the low-rank trigonometric Gaudin model has a significantly richer singularity theory than its integrability alone would indicate.

Taken together, these developments show that “trigonometric Gaudin model” now denotes not a single formula but a coherent domain of integrable structures: XXZ quasi-classical limits, multiplicative current algebras, boundary reflection systems, higher Hamiltonian algebras, non-skew-symmetric quasi-trigonometric deformations, moduli-space families, quantum-classical dualities, and concrete realizations in dissipative and symplectic dynamics. A plausible implication is that future work will continue to treat trigonometric Gaudin systems less as isolated deformations of rational models and more as organizing centers connecting boundary integrability, geometric representation theory, and generalized \(R\)-matrix frameworks [2407.06424].

Source: https://www.emergentmind.com/topics/trigonometric-gaudin-models