Lecture 14

The quantum inverse scattering method, and its classical limit

What this lecture does

Course map: G2 and G5, P2 quantised. Lecture 10 solved the Heisenberg chain by a guess and left the question of why it is solvable. This lecture answers with Lecture 9’s three moves, quantised. The chain’s Lax operator satisfies an RLLRLL relation with Yang’s R\Rmat-matrix of Lecture 13, the monodromy satisfies RTTRTT, and the transfer matrices commute at all values of the spectral parameter. Then the lecture restores ℏ\hbar in the quantum Toda chain and recovers Lecture 9’s r\rmat-matrix as ℏ→0\hbar\to0, sign for sign.

Reading and tutorial. Material marked Reading is studied on your own; material marked Tutorial is worked in the tutorial. Everything else is the lectured core.

14.1 The question

Lecture 9 made the Toda ring’s integrals commute with one identity for one site; Lecture 10 solved the Heisenberg chain by a guess. So the question is

why do the Heisenberg chain’s charges commute, and what classical structure is this the quantisation of?

The answer: Lecture 9’s structure, with commutators for brackets and Yang’s R\Rmat for r\rmat (Theorems 14.1 to 14.3); and the classical limit of the quantum Toda chain proves that the dictionary is right (Corollary 14.6). The worked example is the Heisenberg chain’s Lax operator.

14.2 From brackets to commutators

Lecture 9’s site matrices obey the quadratic bracket (9.9); it composes along an ultralocal chain (Theorem 9.7), and its trace makes the transfer matrices Poisson-commute (Corollary 9.8). Quantisation makes the entries of LL operators and brackets commutators, {⋅,⋅}→[⋅,⋅]/iℏ\pb{\cdot}{\cdot}\to\comm{\cdot}{\cdot}/\ii\hbar, and the legs of Definition 9.2 survive. The quantum relation for the monodromy is the RTTRTT relation R12(λ−μ) T2(μ) T1(λ)=T1(λ) T2(μ) R12(λ−μ),R=1+iℏ r+O(ℏ2). \Rmat_{12}(\lambda-\mu)\,\leg{\mono}{2}(\mu)\,\leg{\mono}{1}(\lambda) =\leg{\mono}{1}(\lambda)\,\leg{\mono}{2}(\mu)\,\Rmat_{12}(\lambda-\mu), \qquad \Rmat=\Id+\ii\hbar\,\rmat+O(\hbar^{2}). (14.1) At order ℏ\hbar it returns (9.9); with the other order of the monodromies it would return r→−r\rmat\to-\rmat, so the order is forced (§14.6). Other books use the other order, each with its own R\Rmat and Lax operator: convert all three together, never one alone. Until §14.6, ℏ=1\hbar=1, and the auxiliary copies of C2\C^{2} are aa and bb, apart from the sites n=1,…,Nn=1,\dots,N.

14.3 The Lax operator of the Heisenberg chain

The quantum space is the Heisenberg chain of Lecture 10, H=(C2)⊗N\mathcal H=(\C^{2})^{\otimes N}, one spin 12\tfrac12 per site. The auxiliary space aa is one more C2\C^{2}. The Lax operator at site nn is Lan(λ)=λ+iPan=(λ+i2)+i σa⋅Sn=(λ+i2+iSnziSn−iSn+λ+i2−iSnz) ⁣a, \Lop{an}(\lambda)=\lambda+\ii P_{an} =\big(\lambda+\tfrac{\ii}{2}\big)+\ii\,\boldsymbol\sigma_a\cdot\mathbf S_n =\begin{pmatrix}\lambda+\tfrac{\ii}2+\ii S^{z}_n&\ii S^{-}_n\\[2pt] \ii S^{+}_n&\lambda+\tfrac{\ii}2-\ii S^{z}_n\end{pmatrix}_{\!a}, (14.2) with PanP_{an} the permutation of the auxiliary space with site nn (Definition 9.4), S=σ/2\mathbf S=\boldsymbol\sigma/2 and S±=Sx±iSyS^{\pm}=S^{x}\pm\ii S^{y}. The middle form uses P=12(1+σa⋅σn)P=\tfrac12(\Id+\boldsymbol\sigma_a\cdot\boldsymbol\sigma_n), the identity behind the Hamiltonian (10.2); the matrix is in the auxiliary space, with entries operators on site nn. Like the Toda site matrix it is linear in the spectral parameter and acts on one site.

Figure 14.1 draws it. Each space is a line: the auxiliary space horizontal, the site vertical. An operator on both is a vertex where they cross, and (14.2) says what happens there: with weight λ\lambda the two states go straight through, and with weight i\ii they swap lines, the P\perm of Definition 9.4.

Figure 14.1. The Lax operator as a vertex. The auxiliary line (indigo) crosses the line of site nn (grey). With weight λ\lambda both states go straight on; with weight i\ii the auxiliary state leaves on the site’s line and the site’s state on the auxiliary line. (Where lines run side by side, a swap is drawn as a crossing; for lines that meet at right angles the crossing is the identity and the swap is two turns.)

Two auxiliary spaces aa, bb and one site nn make three copies of C2\C^{2}, and on them the only operators needed are the three transpositions PabP_{ab}, PanP_{an}, PbnP_{bn}. Their products are permutations of three objects, and a product of two different transpositions is a 33-cycle, of which there are two: C=PabPan=PanPbn=PbnPab,C−1=PanPab=PabPbn=PbnPan. C=P_{ab}P_{an}=P_{an}P_{bn}=P_{bn}P_{ab},\qquad C^{-1}=P_{an}P_{ab}=P_{ab}P_{bn}=P_{bn}P_{an}. (14.3) (Follow a vector x⊗y⊗zx\otimes y\otimes z through each product; this is the computation of Proposition 13.4.)

Theorem 14.1 (The RLL relation). With R(λ)=λ−iP\Rmat(\lambda)=\lambda-\ii P, Rab(λ−μ) Lbn(μ) Lan(λ)=Lan(λ) Lbn(μ) Rab(λ−μ) \Rmat_{ab}(\lambda-\mu)\,\Lop{bn}(\mu)\,\Lop{an}(\lambda) =\Lop{an}(\lambda)\,\Lop{bn}(\mu)\,\Rmat_{ab}(\lambda-\mu) (14.4) for all λ,μ\lambda,\mu. With λ+iP\lambda+\ii P in place of R\Rmat it fails.

Proof. Take R=(λ−μ)+siPab\Rmat=(\lambda-\mu)+s\ii P_{ab} with a sign s=±1s=\pm1 and compute the difference Δ\Delta of the two sides. It has two parts. The first is (λ−μ)(LbnLan−LanLbn)(\lambda-\mu)\big(\Lop{bn}\Lop{an}-\Lop{an}\Lop{bn}\big), and in LbnLan−LanLbn\Lop{bn}\Lop{an}-\Lop{an}\Lop{bn} everything cancels except the products of two PP’s: PanPbn−PbnPan=C−C−1P_{an}P_{bn}-P_{bn}P_{an}=C-C^{-1}. The second is si(PabLbnLan−LanLbnPab)=si(iμ(PabPan−PanPab)+iλ(PabPbn−PbnPab)−(PabPbnPan−PanPbnPab)),s\ii\big(P_{ab}\Lop{bn}\Lop{an}-\Lop{an}\Lop{bn}P_{ab}\big) =s\ii\Big(\ii\mu\big(P_{ab}P_{an}-P_{an}P_{ab}\big) +\ii\lambda\big(P_{ab}P_{bn}-P_{bn}P_{ab}\big) -\big(P_{ab}P_{bn}P_{an}-P_{an}P_{bn}P_{ab}\big)\Big), the terms λμPab\lambda\mu P_{ab} having cancelled. By (14.3) the first bracket is C−C−1C-C^{-1}, the second C−1−CC^{-1}-C, and the triple products are both PbnP_{bn}, so they cancel. The second part is therefore s(λ−μ)(C−C−1)s(\lambda-\mu)(C-C^{-1}), and Δ=(1+s)(λ−μ)(C−C−1).\Delta=(1+s)(\lambda-\mu)\big(C-C^{-1}\big). This vanishes for s=−1s=-1 and is not zero for s=+1s=+1, since C≠C−1C\ne C^{-1}. 

The whole relation is a statement about the permutations of three objects (the check script confirms the residual symbolically).

14.4 From one site to the chain

The monodromy matrix

The monodromy matrix of the chain is the ordered product of the Lax operators along it, Ta(λ)=LaN(λ)⋯La2(λ)La1(λ)=(A(λ)B(λ)C(λ)D(λ)) ⁣a, \mono_a(\lambda)=\Lop{aN}(\lambda)\cdots\Lop{a2}(\lambda)\Lop{a1}(\lambda) =\begin{pmatrix}A(\lambda)&B(\lambda)\\ C(\lambda)&D(\lambda)\end{pmatrix}_{\!a}, (14.5) an operator on Ca2⊗H\C^{2}_a\otimes\mathcal H: in the auxiliary space a 2×22\times2 matrix whose entries A,B,C,DA,B,C,D are operators on the whole chain. As a picture it is a row: one auxiliary line crossing every site in turn (Figure 14.2), Figure 9.1 with operators at the vertices. Operators on different sites commute, so Lan\Lop{an} commutes with Lbm\Lop{bm} for n≠mn\ne m: they act on disjoint sets of spaces. That is the quantum form of ultralocality, and it is all the next proof uses.

Figure 14.2. The monodromy matrix as a row. The auxiliary line meets site 11 first, so La1\Lop{a1} is the rightmost factor. The two free ends of the auxiliary line are the two indices of the 2×22\times2 matrix (14.5); the ends of the site lines are the quantum space.

Theorem 14.2 (RLL implies RTT). The monodromy (14.5) satisfies Rab(λ−μ) Tb(μ) Ta(λ)=Ta(λ) Tb(μ) Rab(λ−μ). \Rmat_{ab}(\lambda-\mu)\,\mono_b(\mu)\,\mono_a(\lambda) =\mono_a(\lambda)\,\mono_b(\mu)\,\Rmat_{ab}(\lambda-\mu). (14.6) The same holds for any chain whose Lax operators obey (14.4) with a common R\Rmat and act on different sites.

Proof. In TbTa=LbN⋯Lb1 LaN⋯La1\mono_b\mono_a=\Lop{bN}\cdots\Lop{b1}\,\Lop{aN}\cdots\Lop{a1}, move each Lan\Lop{an} to the left past Lb,n−1,…,Lb1\Lop{b,n-1},\dots,\Lop{b1}, with which it commutes. The factors pair up by site: Tb(μ)Ta(λ)=(LbNLaN)(Lb,N−1La,N−1)⋯(Lb1La1).\mono_b(\mu)\mono_a(\lambda) =\big(\Lop{bN}\Lop{aN}\big)\big(\Lop{b,N-1}\Lop{a,N-1}\big)\cdots \big(\Lop{b1}\Lop{a1}\big). Now Rab\Rmat_{ab} enters on the left. By (14.4) at site NN, RabLbNLaN=LaNLbNRab\Rmat_{ab}\Lop{bN}\Lop{aN}=\Lop{aN}\Lop{bN}\Rmat_{ab}: it has passed the first pair and reversed it. The same step at each site in turn carries Rab\Rmat_{ab} to the right end, leaving (LaNLbN)⋯(La1Lb1) Rab(\Lop{aN}\Lop{bN})\cdots(\Lop{a1}\Lop{b1})\,\Rmat_{ab}, and undoing the pairing gives TaTbRab\mono_a\mono_b\Rmat_{ab}. 

This is the train argument (Figure 14.3): the crossing of the two auxiliary lines is pushed through the row one site at a time, as the bracket was pushed through the product by Leibniz’s rule in Theorem 9.7.

Figure 14.3. The train argument. Operators act left to right along the auxiliary lines, so the leftmost thing drawn is the rightmost factor. Bottom: the rows, then the crossing Rab\Rmat_{ab}. By (14.4) the crossing slides through site NN, then N−1N-1, …, until it has passed the row (top): (14.6).

The Lax operator is an R\Rmat-matrix, and RLL is the Yang–Baxter equation

Compare (14.2) with R\Rmat: Lan(λ)=λ+iPan=−Ran(−λ). \Lop{an}(\lambda)=\lambda+\ii P_{an}=-\Rmat_{an}(-\lambda). (14.7) The Lax operator is the R\Rmat-matrix itself, with the site as its second space. It is not Ran(λ)\Rmat_{an}(\lambda): the two differ by λ→−λ\lambda\to-\lambda. Substitute (14.7) into (14.4); the two minus signs cancel, Rab(λ−μ) Rbn(−μ) Ran(−λ)=Ran(−λ) Rbn(−μ) Rab(λ−μ).\Rmat_{ab}(\lambda-\mu)\,\Rmat_{bn}(-\mu)\,\Rmat_{an}(-\lambda) =\Rmat_{an}(-\lambda)\,\Rmat_{bn}(-\mu)\,\Rmat_{ab}(\lambda-\mu). Rename b→1b\to1, a→2a\to2, n→3n\to3, and put u=−μu=-\mu, v=−λv=-\lambda, so that λ−μ=u−v\lambda-\mu=u-v. Since P12=P21P_{12}=P_{21} this is R12(u−v) R13(u) R23(v)=R23(v) R13(u) R12(u−v), \Rmat_{12}(u-v)\,\Rmat_{13}(u)\,\Rmat_{23}(v) =\Rmat_{23}(v)\,\Rmat_{13}(u)\,\Rmat_{12}(u-v), (14.8) the Yang–Baxter equation in its standard form. With u=u13u=u_{13}, v=u23v=u_{23} and u−v=u12u-v=u_{12} it is the equation (13.6) that Theorem 13.3 derived from three-body consistency, for the matrix XB(u)=R(u)/(u+i)X_{\mathrm B}(u)=\Rmat(u)/(u+\ii) of the delta bosons (13.3). So Theorem 14.1 and Proposition 13.4 are one fact, read twice. In Lecture 13 the three lines were three particles and R\Rmat a two-body collision; in the Heisenberg chain one of the three lines is a site. The fermions of Lecture 13 scatter with XF=−L(u)/(u+i)X_{\mathrm F}=-\Lop{}(u)/(u+\ii), the Lax operator itself.

Figure 14.5 shows the equation as a picture. It is Figure 13.1, and Figure 14.4 is its special case with one line vertical.

Figure 14.4. The RLLRLL relation (14.4). Left: the site line nn meets aa (lower) and then bb, and afterwards the auxiliary lines cross. Right: they cross first, and nn meets bb and then aa. The crossing of aa and bb may be slid through the site line.
Figure 14.5. The Yang–Baxter equation (14.8) as a slide: line 22 passes either above or below the crossing of lines 11 and 33. Lines carry spectral parameters and a crossing of lines i,ji,j carries Rij\Rmat_{ij} at the difference of their parameters, which is why u13=u12+u23u_{13}=u_{12}+u_{23}. Reading left to right, the first crossing is the rightmost factor. Make line 33 a site and this is Figure 14.4.

14.5 Commuting transfer matrices

The transfer matrix is the trace of the monodromy over the auxiliary space, t(λ)=tra⁡Ta(λ)=A(λ)+D(λ), \transfer(\lambda)=\tra\mono_a(\lambda)=A(\lambda)+D(\lambda), (14.9) an operator on the chain alone: in the picture the auxiliary line is closed into a loop (Figure 14.6). It is a polynomial of degree NN in λ\lambda with operator coefficients, and leading term 2λN2\lambda^{N}.

Figure 14.6. The transfer matrix: the auxiliary line of Figure 14.2, closed into a loop by the trace. The loop has no ends, so it acts on the sites alone; this is the loop of Figure 9.1 with operators at the vertices.

Theorem 14.3 (Commuting transfer matrices). [t(λ),t(μ)]=0\comm{\transfer(\lambda)}{\transfer(\mu)}=0 for all λ,μ∈C\lambda,\mu\in\C.

Proof. R(u)R(−u)=(u−iP)(−u−iP)=−(u2+1)\Rmat(u)\Rmat(-u)=(u-\ii P)(-u-\ii P)=-(u^{2}+1), so R(u)\Rmat(u) is invertible for u≠±iu\ne\pm\ii. For such λ−μ\lambda-\mu, multiply (14.6) on the left by Rab−1\Rmat_{ab}^{-1}: Tb(μ)Ta(λ)=Rab−1 Ta(λ)Tb(μ) Rab\mono_b(\mu)\mono_a(\lambda)=\Rmat_{ab}^{-1}\,\mono_a(\lambda)\mono_b(\mu)\,\Rmat_{ab}. Take the trace over both auxiliary spaces. On the left, tr⁡abTbTa=t(μ)t(λ)\operatorname{tr}_{ab}\mono_b\mono_a=\transfer(\mu)\transfer(\lambda), because Ta\mono_a and Tb\mono_b have their auxiliary legs in different spaces. On the right, Rab\Rmat_{ab} acts on the auxiliary spaces only, as a numerical matrix, and the partial trace is cyclic for such a factor; so tr⁡ab(Rab−1TaTbRab)=tr⁡ab(TaTb)=t(λ)t(μ)\operatorname{tr}_{ab}\big(\Rmat_{ab}^{-1}\mono_a\mono_b\Rmat_{ab}\big) =\operatorname{tr}_{ab}\big(\mono_a\mono_b\big)=\transfer(\lambda)\transfer(\mu). Both sides are polynomials in λ,μ\lambda,\mu, so the identity extends to λ−μ=±i\lambda-\mu=\pm\ii. 

This is Corollary 9.8: there the trace of a commutator vanished, here the trace of a conjugate is the trace (Figure 14.7). The coefficients of t(λ)\transfer(\lambda) commute with one another. Whether the family contains a Hamiltonian worth having, and whether its members are local, is the business of Lecture 15.

Figure 14.7. Why transfer matrices commute. Insert RR−1=1\Rmat\Rmat^{-1}=\Id between the loops; R\Rmat slides through every site (Figure 14.3), exchanging the rows, and meets R−1\Rmat^{-1} round the back. The loops have passed through each other.

What RTTRTT gives, and what it does not

The proof used (14.6) and the invertibility of R\Rmat, and not the Yang–Baxter equation. The equation does two other jobs. It supplies Lax operators: by (14.7) the rational R\Rmat, read with a site as its second space, is a Lax operator for its own RLLRLL relation, which is why the Heisenberg chain has one. And it is sufficient, not forced, for the consistency of the RTTRTT relations taken as the definition of an algebra (the readings at the end). Finally, integrability does not require RTTRTT: the free-fermion charges of the Ising chain (Lecture 18) are constructed directly.

14.6 The classical limit, verified: the quantum Toda chain

The Heisenberg chain has no classical limit worth the name — a spin 12\tfrac12 has no small parameter. To test the correspondence of §14.2, quantise Lecture 9’s chain. Restore ℏ\hbar. The periodic Toda chain has canonical pairs with [qm,pn]=iℏ δmn\comm{q_m}{p_n}=\ii\hbar\,\delta_{mn}, represented by pn=−iℏ ∂/∂qnp_n=-\ii\hbar\,\partial/\partial q_n on functions of qq. Keep the site matrix of (9.2) as it stands, Ln(λ)=(λ−pne−qn−eqn0), \Lop{n}(\lambda)=\begin{pmatrix}\lambda-p_n&e^{-q_n}\\ -e^{q_n}&0\end{pmatrix}, (14.10) now with operator entries. No ordering question arises: each entry involves qnq_n or pnp_n, never both.

Theorem 14.4 (Quantum Toda: RLL). The site matrix (14.10) satisfies Rab(λ−μ) Lb(μ) La(λ)=La(λ) Lb(μ) Rab(λ−μ),R(λ)=λ−iℏP,\Rmat_{ab}(\lambda-\mu)\,\Lop{b}(\mu)\,\Lop{a}(\lambda) =\Lop{a}(\lambda)\,\Lop{b}(\mu)\,\Rmat_{ab}(\lambda-\mu),\qquad \Rmat(\lambda)=\lambda-\ii\hbar P, exactly, to all orders in ℏ\hbar, and fails with λ+iℏP\lambda+\ii\hbar P. Consequently the monodromy T=LN⋯L1\mono=\Lop{N}\cdots\Lop{1} satisfies RTTRTT, and [t(λ),t(μ)]=0\comm{\transfer(\lambda)}{\transfer(\mu)}=0.

Proof. Write L(λ)=λE11+L0\Lop{}(\lambda)=\lambda E_{11}+\Lop{}^{0} with L0=L(0)\Lop{}^{0}=\Lop{}(0). In the leg notation the (ij),(kl)(ij),(kl) entry of La(λ)Lb(μ)\Lop{a}(\lambda)\Lop{b}(\mu) is Lik(λ)Ljl(μ)L_{ik}(\lambda)L_{jl}(\mu), and of Lb(μ)La(λ)\Lop{b}(\mu)\Lop{a}(\lambda) it is the same product in the other order. So their difference is the matrix of commutators −[Lik,Ljl]-\comm{L_{ik}}{L_{jl}}, which do not depend on λ,μ\lambda,\mu. The only non-zero ones are [p,e∓q]=±iℏ e∓q\comm{p}{e^{\mp q}}=\pm\ii\hbar\,e^{\mp q}, exactly iℏ\ii\hbar times the Poisson brackets used in the proof of Proposition 9.9; that proof computed the whole matrix of brackets, so Lb(μ)La(λ)−La(λ)Lb(μ)=−iℏ (E11⊗L0−L0⊗E11)P.\Lop{b}(\mu)\Lop{a}(\lambda)-\Lop{a}(\lambda)\Lop{b}(\mu) =-\ii\hbar\,\big(E_{11}\otimes\Lop{}^{0}-\Lop{}^{0}\otimes E_{11}\big)P . Next, P Lb(μ)La(λ) P=La(μ)Lb(λ)P\,\Lop{b}(\mu)\Lop{a}(\lambda)\,P=\Lop{a}(\mu)\Lop{b}(\lambda) (conjugating by PP exchanges the legs), and expanding L=λE11+L0\Lop{}=\lambda E_{11}+\Lop{}^{0} gives La(μ)Lb(λ)−La(λ)Lb(μ)=(μ−λ)(E11⊗L0−L0⊗E11)\Lop{a}(\mu)\Lop{b}(\lambda)-\Lop{a}(\lambda)\Lop{b}(\mu) =(\mu-\lambda)\big(E_{11}\otimes\Lop{}^{0}-\Lop{}^{0}\otimes E_{11}\big); the L0⊗L0\Lop{}^{0}\otimes\Lop{}^{0} terms, in the same order on both sides, cancel. With R=(λ−μ)+siℏP\Rmat=(\lambda-\mu)+s\ii\hbar P the difference of the two sides of the relation is (λ−μ)(LbLa−LaLb)+siℏ(PLbLa−LaLbP)=−(1+s) iℏ (λ−μ)(E11⊗L0−L0⊗E11)P,(\lambda-\mu)\big(\Lop{b}\Lop{a}-\Lop{a}\Lop{b}\big) +s\ii\hbar\big(P\Lop{b}\Lop{a}-\Lop{a}\Lop{b}P\big) =-(1+s)\,\ii\hbar\,(\lambda-\mu)\big(E_{11}\otimes\Lop{}^{0}-\Lop{}^{0}\otimes E_{11}\big)P, using PLbLa=La(μ)Lb(λ)PP\Lop{b}\Lop{a}=\Lop{a}(\mu)\Lop{b}(\lambda)P. It vanishes for s=−1s=-1 only. Sites commute, so Theorem 14.2 and Theorem 14.3 apply verbatim, now with R(u)R(−u)=−(u2+ℏ2)\Rmat(u)\Rmat(-u)=-(u^{2}+\hbar^{2}). 

The coefficients of t(λ)=∑kckλN−k\transfer(\lambda)=\sum_kc_k\lambda^{N-k} are the quantum integrals. They are the expressions of (9.4), now as operators, with no ordering ambiguity since each term has at most one factor per site: c1=−∑pnc_1=-\sum p_n, and H=12c12−c2=12∑pn2+∑eqn−qn+1H=\tfrac12c_1^{2}-c_2=\tfrac12\sum p_n^{2}+\sum e^{q_n-q_{n+1}} is the quantum periodic Toda Hamiltonian, a function of the commuting ckc_k. Its spectral problem was reduced to one-variable quantisation conditions, for small NN by Gutzwiller (1981) and for all NN by Sklyanin’s separation of variables (1985) cited.

ℏ→0\hbar\to0

To compare with Lecture 9, divide by λ−μ\lambda-\mu: R(λ−μ)λ−μ=1+iℏ r(λ−μ),r(λ)=−Pλ,\frac{\Rmat(\lambda-\mu)}{\lambda-\mu}=\Id+\ii\hbar\,\rmat(\lambda-\mu), \qquad \rmat(\lambda)=-\frac{\perm}{\lambda}, exactly the r\rmat of Proposition 9.9, and (14.1)’s expansion with no O(ℏ2)O(\hbar^{2}) term. The limit needs a rule for turning operators into functions.

Lemma 14.5 (Symbols). Write an operator polynomial in the pp’s with all pp’s on the right, F=∑αfα(q) pαF=\sum_\alpha f_\alpha(q)\,p^{\alpha}, and call σ(F)=∑αfα(q) pα\sigma(F)=\sum_\alpha f_\alpha(q)\,p^{\alpha}, a function of (q,p)(q,p), its symbol. Then σ(FG)=∑α(−iℏ)∣α∣α! ∂pασ(F) ∂qασ(G),soσ([F,G])=iℏ {σ(F),σ(G)}+O(ℏ2).\sigma(FG)=\sum_{\alpha}\frac{(-\ii\hbar)^{\abs\alpha}}{\alpha!}\, \partial_p^{\alpha}\sigma(F)\,\partial_q^{\alpha}\sigma(G), \qquad\text{so}\qquad \sigma\big(\comm FG\big)=\ii\hbar\,\pb{\sigma(F)}{\sigma(G)}+O(\hbar^{2}).

Proof (★). Move each pn=−iℏ ∂qnp_n=-\ii\hbar\,\partial_{q_n} in FF to the right through GG’s functions of qq by Leibniz’s rule, pnkg=∑j(kj)(−iℏ)j(∂qnjg) pnk−jp_n^{k}g=\sum_j\binom kj(-\ii\hbar)^{j}(\partial_{q_n}^{j}g)\,p_n^{k-j}. The first-order terms of σ(FG)−σ(GF)\sigma(FG)-\sigma(GF) are −iℏ∑n(∂pnf ∂qng−∂png ∂qnf)=iℏ{f,g}-\ii\hbar\sum_n\big(\partial_{p_n}f\,\partial_{q_n}g-\partial_{p_n}g\,\partial_{q_n}f\big) =\ii\hbar\pb fg in the course’s convention. 

Corollary 14.6 (The classical limit of RTTRTT). The symbols of the entries of the quantum monodromy are the entries of the classical monodromy (9.3). At order ℏ\hbar the RTTRTT relation of Theorem 14.4 becomes {T1(λ),T2(μ)}=[r(λ−μ),T1(λ)T2(μ)],r(λ)=−Pλ:\pb{\leg{\mono}{1}(\lambda)}{\leg{\mono}{2}(\mu)} =\comm{\rmat(\lambda-\mu)}{\leg{\mono}{1}(\lambda)\leg{\mono}{2}(\mu)}, \qquad\rmat(\lambda)=-\frac{\perm}{\lambda}: Proposition 9.9 and Theorem 9.7, sign for sign.

Proof. Each entry of T\mono is a sum of products of entries of different site matrices, which commute, so its symbol is the classical product. Write RTTRTT as TaTb−TbTa=iℏ(r TbTa−TaTb r)\mono_a\mono_b-\mono_b\mono_a=\ii\hbar\big(\rmat\,\mono_b\mono_a-\mono_a\mono_b\,\rmat\big) and take symbols. On the left the (ij),(kl)(ij),(kl) entry is the commutator [Tik(λ),Tjl(μ)]\comm{T_{ik}(\lambda)}{T_{jl}(\mu)}, whose symbol is iℏ{Tik,Tjl}+O(ℏ2)\ii\hbar\pb{T_{ik}}{T_{jl}}+O(\hbar^{2}) by Lemma 14.5: the entry of iℏ{T1,T2}\ii\hbar\pb{\leg{\mono}{1}}{\leg{\mono}{2}} in Definition 9.2. On the right the symbols of TbTa\mono_b\mono_a and TaTb\mono_a\mono_b are both T1T2+O(ℏ)\leg{\mono}{1}\leg{\mono}{2}+O(\hbar). Divide by iℏ\ii\hbar and let ℏ→0\hbar\to0. 

Tracing over both legs, {t(λ),t(μ)}=0\pb{\transfer(\lambda)}{\transfer(\mu)}=0 is the order-ℏ\hbar term of Theorem 14.3. With λ+iℏP\lambda+\ii\hbar P, false for this chain, the same steps would give r=+P/λ\rmat=+\perm/\lambda: the classical side decides the quantum sign. The check script repeats this computation for two and three sites.

14.7 What it buys

  • A commuting family for the Heisenberg chain, from one relation for one site (Theorems 14.1 to 14.3). Lecture 15 obtains the chain’s Hamiltonian from its logarithmic derivative and diagonalises the family.

  • The Yang–Baxter equation, read a second time: the Lax operator is Yang’s R\Rmat with a site as one of its spaces, and RLLRLL is the equation of Lecture 13 (14.8).

  • G5, closed: the quantum Toda chain’s RTTRTT relation becomes Lecture 9’s quadratic bracket as ℏ→0\hbar\to0, sign for sign (Corollary 14.6).

14.8 Caveats, in one place

  • Status. Theorems 14.1 to 14.4 and Corollary 14.6 are proved; the proof of Lemma 14.5 is a reading. The quantum Toda chain’s spectral problem is cited.

  • What RTTRTT does not give. It gives commuting operators, not a Hamiltonian worth having or locality (Lecture 15), and not completeness, which stays Lecture 10’s cited question. The independence of the quantum Toda integrals is not shown.

  • The classical limit is tested on Toda; the Heisenberg chain has no small parameter.

Tutorial

14.9 Two sites by hand

For N=2N=2, t(λ)=tra⁡[(λ+iPa2)(λ+iPa1)]=2λ2+2iλ−tra⁡(Pa2Pa1)\transfer(\lambda)=\tra\big[(\lambda+\ii P_{a2})(\lambda+\ii P_{a1})\big] =2\lambda^{2}+2\ii\lambda-\tra\big(P_{a2}P_{a1}\big), using tra⁡1=2\tra\Id=2 and tra⁡Pan=1\tra P_{an}=\Id. By (14.3), Pa2Pa1=P21Pa2P_{a2}P_{a1}=P_{21}P_{a2}, so tra⁡(Pa2Pa1)=P12\tra(P_{a2}P_{a1})=P_{12} and t(λ)=2λ2+2iλ−P12,\transfer(\lambda)=2\lambda^{2}+2\ii\lambda-P_{12}, which commutes with itself at any two values, as it must: it involves only P12P_{12}. On the triplet P12=1P_{12}=1, on the singlet −1-1. Then t(0)=−P12\transfer(0)=-P_{12}, the swap of the two sites, and t(0)−1t′(0)=(−P12)(2i)=−2iP12\transfer(0)^{-1}\transfer'(0)=(-P_{12})(2\ii)=-2\ii P_{12}: Lecture 15 finds the chain’s Hamiltonian in this logarithmic derivative. Repeat for N=3N=3 and check that t(0)=−i P32P31\transfer(0)=-\ii\,P_{32}P_{31}.

Reading

14.10 The classical Yang–Baxter equation

The linear Sklyanin bracket (9.11) of Lecture 9’s readings is an assumption about a Poisson bracket, so it must be consistent with the Jacobi identity. Working out what that requires produces the classical form of Lecture 13’s equation.

Put three legs in play, V⊗V⊗VV\otimes V\otimes V, write Li\leg{L}{i} for LL in leg ii at spectral parameter λi\lambda_i, and abbreviate rij=rij(λi−λj)\leg{\rmat}{ij}=\leg{\rmat}{ij}(\lambda_i-\lambda_j). Two facts do the work: [rij,Lk]=0\comm{\leg{\rmat}{ij}}{\leg{L}{k}}=0 when k∉{i,j}k\notin\{i,j\}, because they act on disjoint legs; and r21=−r12\leg{\rmat}{21}=-\leg{\rmat}{12}, which Remark 9.13 showed to be sufficient for antisymmetry of the Poisson bracket but which the bracket does not by itself force. It is therefore a hypothesis of Theorem 14.7, not a consequence of Definition 9.12, and is assumed from here on.

Theorem 14.7 (Jacobi and the CYBE). For a Lax matrix satisfying (9.11) with an r\rmat obeying the skew condition (9.12), the Jacobi identity for L1,L2,L3\leg{L}{1},\leg{L}{2},\leg{L}{3} is equivalent to [  [r12,r13]+[r12,r23]+[r13,r23]  ,  L1+L2+L3  ]=0. \comm{\;\comm{\leg{\rmat}{12}}{\leg{\rmat}{13}} +\comm{\leg{\rmat}{12}}{\leg{\rmat}{23}} +\comm{\leg{\rmat}{13}}{\leg{\rmat}{23}}\;} {\;\leg{L}{1}+\leg{L}{2}+\leg{L}{3}\;}=0 . (14.11) In particular it holds whenever r\rmat satisfies the classical Yang–Baxter equation [r12,r13]+[r12,r23]+[r13,r23]=0. \comm{\leg{\rmat}{12}}{\leg{\rmat}{13}} +\comm{\leg{\rmat}{12}}{\leg{\rmat}{23}} +\comm{\leg{\rmat}{13}}{\leg{\rmat}{23}}=0 . (14.12)

Proof. Since r\rmat is numerical it passes through a Poisson bracket, so {{L1,L2},L3}=[r12,{L1,L3}+{L2,L3}]=[r12,[r13,L1+L3]+[r23,L2+L3]].\pb{\pb{\leg{L}{1}}{\leg{L}{2}}}{\leg{L}{3}} =\comm{\leg{\rmat}{12}}{\pb{\leg{L}{1}}{\leg{L}{3}}+\pb{\leg{L}{2}}{\leg{L}{3}}} =\comm{\leg{\rmat}{12}}{\comm{\leg{\rmat}{13}}{\leg{L}{1}+\leg{L}{3}} +\comm{\leg{\rmat}{23}}{\leg{L}{2}+\leg{L}{3}}} . Use [r12,L3]=0\comm{\leg{\rmat}{12}}{\leg{L}{3}}=0 and the matrix Jacobi identity in the form [A,[B,C]]=[[A,B],C]+[B,[A,C]]\comm{A}{\comm{B}{C}}=\comm{\comm{A}{B}}{C}+\comm{B}{\comm{A}{C}} on each of the two terms. The first gives [r12,[r13,L1+L3]]=[[r12,r13],L1+L3]+[r13,[r12,L1]],\comm{\leg{\rmat}{12}}{\comm{\leg{\rmat}{13}}{\leg{L}{1}+\leg{L}{3}}} =\comm{\comm{\leg{\rmat}{12}}{\leg{\rmat}{13}}}{\leg{L}{1}+\leg{L}{3}} +\comm{\leg{\rmat}{13}}{\comm{\leg{\rmat}{12}}{\leg{L}{1}}}, and the second the same with (13)→(23)(13)\to(23), L1→L2\leg{L}{1}\to\leg{L}{2}. Now add the two cyclic permutations of this, replacing r21,r31,r32\leg{\rmat}{21},\leg{\rmat}{31},\leg{\rmat}{32} by −r12,−r13,−r23-\leg{\rmat}{12},-\leg{\rmat}{13},-\leg{\rmat}{23} as (9.12) allows. The six leftover double commutators pair up, again by matrix Jacobi: for instance [r13,[r12,L1]]−[r12,[r13,L1]]=−[[r12,r13],L1].\comm{\leg{\rmat}{13}}{\comm{\leg{\rmat}{12}}{\leg{L}{1}}} -\comm{\leg{\rmat}{12}}{\comm{\leg{\rmat}{13}}{\leg{L}{1}}} =-\comm{\comm{\leg{\rmat}{12}}{\leg{\rmat}{13}}}{\leg{L}{1}} . Now collect. The first group contributes each double commutator of r\rmat’s twice, once from each of the two cyclic terms that produce it, and the two contributions carry different pairs of Li\leg{L}{i}: [r12,r13]: L1+L2 and L1+L3;[r12,r23]: L1+L2 and L2+L3;[r13,r23]: L1+L3 and L2+L3,\comm{\leg{\rmat}{12}}{\leg{\rmat}{13}}:\ \leg{L}{1}+\leg{L}{2}\ \text{and}\ \leg{L}{1}+\leg{L}{3};\quad \comm{\leg{\rmat}{12}}{\leg{\rmat}{23}}:\ \leg{L}{1}+\leg{L}{2}\ \text{and}\ \leg{L}{2}+\leg{L}{3};\quad \comm{\leg{\rmat}{13}}{\leg{\rmat}{23}}:\ \leg{L}{1}+\leg{L}{3}\ \text{and}\ \leg{L}{2}+\leg{L}{3}, so their totals are 2L1+L2+L32\leg{L}{1}+\leg{L}{2}+\leg{L}{3}, L1+2L2+L3\leg{L}{1}+2\leg{L}{2}+\leg{L}{3} and L1+L2+2L3\leg{L}{1}+\leg{L}{2}+2\leg{L}{3} respectively. The three surviving pairs of double commutators, computed as displayed above, contribute exactly −[[r12,r13],L1]-\comm{\comm{\leg{\rmat}{12}}{\leg{\rmat}{13}}}{\leg{L}{1}}, −[[r12,r23],L2]-\comm{\comm{\leg{\rmat}{12}}{\leg{\rmat}{23}}}{\leg{L}{2}} and −[[r13,r23],L3]-\comm{\comm{\leg{\rmat}{13}}{\leg{\rmat}{23}}}{\leg{L}{3}}, which remove the one doubled term in each line. All three are left with L1+L2+L3\leg{L}{1}+\leg{L}{2}+\leg{L}{3}, and that is (14.11). 

Remark 14.8 (Machine-checked). The collection above is exactly the kind of bookkeeping in which a sign is lost. It is verified in scripts/check-l8-algebra.py, for d=2d=2, with three explicit legs and symbolic entries throughout — sixty free parameters, so the result is an identity and not a sample. The Jacobi sum equals [X,L1+L2+L3]\comm{X}{\leg{L}{1}+\leg{L}{2}+\leg{L}{3}}, and — the point of Remark 14.9 — is not identically zero on its own.

Remark 14.9 (Jacobi does not force the CYBE, and saying so matters). Theorem 14.7 gives (14.11), and (14.12) is sufficient for it, not implied by it: the Jacobi identity only requires the triple sum to commute with L1+L2+L3\leg{L}{1}+\leg{L}{2}+\leg{L}{3}, and for a particular LL that is a weaker demand. What makes (14.12) the right thing to impose is that it mentions no LL at all. One solution r\rmat then serves every model whose Lax matrix it fits, which is why the classification of r\rmat-matrices is a subject and the classification of integrable models is not.

One variant deserves its name in passing. Replacing the zero on the right of (14.12) by a specific r\rmat-independent term gives the modified classical Yang–Baxter equation, whose solutions generate the same structures over a slightly larger class of Lie algebras. Its theory belongs to Lie-bialgebra and Poisson–Lie theory, which this course does not enter; it is named here so that meeting the phrase elsewhere is not a surprise.

Figure 14.8. The classical Yang–Baxter equation. Three lines cross pairwise, and the two ways of resolving the triple crossing — sliding line 22 past the crossing of 11 and 33 either before or after — give the same answer. This is the classical form of Figure 13.1.

Reading

The rational solution

Proposition 14.10 (The rational r\rmat-matrix). r(λ)=c P/λ\rmat(\lambda)=c\,\perm/\lambda satisfies the CYBE (14.12) for every constant cc.

Proof. Write λij=λi−λj\lambda_{ij}=\lambda_i-\lambda_j. The permutation operators satisfy P12P13=P13P23=P23P12\leg{\perm}{12}\leg{\perm}{13}=\leg{\perm}{13}\leg{\perm}{23} =\leg{\perm}{23}\leg{\perm}{12} and P13P12=P12P23=P23P13\leg{\perm}{13}\leg{\perm}{12}=\leg{\perm}{12}\leg{\perm}{23} =\leg{\perm}{23}\leg{\perm}{13} — both are the two ways of cycling three legs. Calling the difference of those two products KK, all three commutators are ±K\pm K: [P12,P13]=K,[P12,P23]=−K,[P13,P23]=K.\comm{\leg{\perm}{12}}{\leg{\perm}{13}}=K,\quad \comm{\leg{\perm}{12}}{\leg{\perm}{23}}=-K,\quad \comm{\leg{\perm}{13}}{\leg{\perm}{23}}=K . So the left of (14.12) is c2K(1λ12λ13−1λ12λ23+1λ13λ23)=c2K λ23−λ13+λ12λ12λ13λ23=0,c^{2}K\left(\frac{1}{\lambda_{12}\lambda_{13}}-\frac{1}{\lambda_{12}\lambda_{23}} +\frac{1}{\lambda_{13}\lambda_{23}}\right) =c^{2}K\,\frac{\lambda_{23}-\lambda_{13}+\lambda_{12}} {\lambda_{12}\lambda_{13}\lambda_{23}}=0, because λ12+λ23=λ13\lambda_{12}+\lambda_{23}=\lambda_{13} identically. 

The cancellation is a statement about the three spectral parameters and nothing else, which is why cc is free: the normalisation of r\rmat is fixed by the Lax matrix it is paired with, not by the CYBE. Lecture 9 fixes it for the periodic Toda chain, c=−1c=-1 in (9.10); the open chain of Lecture 9’s readings needs a different, constant r\rmat.

Two further families solve (14.12): the trigonometric r\rmat-matrices, in which 1/λ1/\lambda is replaced by cot⁡λ\cot\lambda or coth⁡λ\coth\lambda and the permutation is deformed, and the elliptic ones, built from Weierstrass functions. Belavin and Drinfeld classified the solutions under hypotheses this course has not set up: for r\rmat valued in g⊗g\glie\otimes\glie with g\glie a finite-dimensional simple Lie algebra, in difference form, and non-degenerate in their technical sense (the residue at the simple pole is proportional to the Casimir of an invariant form), every solution is equivalent to a rational, trigonometric or elliptic one cited, and the three degenerate into one another in limits. Those hypotheses are doing real work — dropped, the statement is false, and already fails for d=1d=1, where every commutator in (14.12) vanishes and any function of λ\lambda with a simple pole at the origin solves it. The same trichotomy reappears in the Calogero–Moser potentials (optional lecture B) and in the quantum chains XXX, XXZ and XYZ — it is the same list every time, which is not a coincidence.

Reading

The open chain’s r\rmat is not a solution

The open chain’s constant r\rmat of (9.22) (Lecture 9’s readings) is skew, PrP=−r\perm\rmat\perm=-\rmat, which is (9.12) for a constant r\rmat; so Theorem 14.7 applies to it. Does it solve the classical Yang–Baxter equation (14.12)? No. Its left-hand side is [r12,r13]+[r12,r23]+[r13,r23]=−14[P12,P13]  ≠0. \comm{\leg{\rmat}{12}}{\leg{\rmat}{13}}+\comm{\leg{\rmat}{12}}{\leg{\rmat}{23}} +\comm{\leg{\rmat}{13}}{\leg{\rmat}{23}} =-\tfrac14\comm{\leg{\perm}{12}}{\leg{\perm}{13}}\;\ne0 . (14.13) This identity is cited, not proved here: Babelon–Bernard–Talon, eq. (4.32), state it in operator form, and scripts/check-l9-algebra.py checks it symbolically for N=2,3N=2,3. An equation of this shape, with a fixed non-zero right-hand side in place of the zero, is the modified CYBE named above.

And yet the Jacobi identity holds, as it must: the bracket is the canonical one on R2N\R^{2N}. Theorem 14.7 asks only that the left of (14.13) commute with L1+L2+L3\leg{L}{1}+\leg{L}{2}+\leg{L}{3}, and it does, for every matrix LL: P12\leg{\perm}{12} commutes with L1+L2\leg{L}{1}+\leg{L}{2} (it only swaps them) and with L3\leg{L}{3} (another leg), so it commutes with the sum, and so does P13\leg{\perm}{13}, and so does their commutator. This is Remark 14.9 in a real model: the CYBE is sufficient for Jacobi, and the Toda chain satisfies Jacobi without it.

How unique is r\rmat? Adding x1+yPx\Id+y\perm changes nothing in (9.21), since both commute with L1+L2\leg{L}{1}+\leg{L}{2} (P\perm only swaps the two legs). Babelon–Bernard–Talon, §4.7, derive (9.22) rather than check it, and show it is the only skew constant solution; scripts/check-rmatrix.py (V8) finds, at sampled points for N=3N=3, no constant solutions beyond r+x1+yP\rmat+x\Id+y\perm. Of these, the skew one is r\rmat itself, which fails the CYBE, while r±12P\rmat\pm\tfrac12\perm do solve the CYBE (Babelon–Bernard–Talon, (4.33); checked for N=2,3N=2,3) and are not skew, so Theorem 14.7 as proved does not cover them. So no constant r\rmat for the open chain is both skew and a solution of the CYBE: the route “skew plus CYBE, hence Jacobi” is closed here, and Jacobi holds by the weaker condition above. That is not a defect of the model; it is why Lie-bialgebra theory has a modified equation.

Reading

14.11 Reordering three monodromies

Three auxiliary spaces 1,2,31,2,3 at parameters λ1,λ2,λ3\lambda_1,\lambda_2,\lambda_3 give three relations (14.6), one for each pair, with Rij=Rij(λi−λj)\Rmat_{ij}=\Rmat_{ij}(\lambda_i-\lambda_j).

Proposition 14.11 (Reordering three monodromies). Let T\mono satisfy (14.6) with an invertible R\Rmat. Reversing T1T2T3\mono_1\mono_2\mono_3 one neighbouring pair at a time gives, along the two possible routes, T1T2T3=WA (T3T2T1) WA−1=WB (T3T2T1) WB−1,WA=R12R13R23,WB=R23R13R12.\mono_1\mono_2\mono_3=W_A\,\big(\mono_3\mono_2\mono_1\big)\,W_A^{-1} =W_B\,\big(\mono_3\mono_2\mono_1\big)\,W_B^{-1},\qquad W_A=\Rmat_{12}\Rmat_{13}\Rmat_{23},\quad W_B=\Rmat_{23}\Rmat_{13}\Rmat_{12}. Hence WB−1WAW_B^{-1}W_A commutes with every T3T2T1\mono_3\mono_2\mono_1.

Proof. Write (14.6) as TiTj=Rij TjTi Rij−1\mono_i\mono_j=\Rmat_{ij}\,\mono_j\mono_i\,\Rmat_{ij}^{-1} for i<ji<j. Route A exchanges the pairs 1212, then 1313, then 2323: T1T2T3=R12T2T1T3R12−1=R12R13T2T3T1R13−1R12−1=WA T3T2T1 WA−1\mono_1\mono_2\mono_3=\Rmat_{12}\mono_2\mono_1\mono_3\Rmat_{12}^{-1} =\Rmat_{12}\Rmat_{13}\mono_2\mono_3\mono_1\Rmat_{13}^{-1}\Rmat_{12}^{-1} =W_A\,\mono_3\mono_2\mono_1\,W_A^{-1}, each Rij\Rmat_{ij} commuting with the monodromy it does not touch. Route B exchanges 2323, 1313, 1212 and gives WBW_B in the same way. Equating the two results gives the statement. 

For a model that satisfies RTTRTT this is automatic — the two routes rewrite one operator. It matters when RTTRTT is taken as the definition of an algebra generated by the entries of T\mono, as a Sklyanin bracket can be posited (Lecture 9’s readings): then the two routes yield the commutation condition of the proposition among the generators, and WA=WBW_A=W_B makes it automatic. That is sufficient, not necessary: for R=σx⊗σz\Rmat=\sigma^{x}\otimes \sigma^{z}, which fails the Yang–Baxter equation, WB=−WAW_B=-W_A and the condition is again automatic (the check script). It is Remark 14.9 quantised: there the Jacobi identity required a combination of r\rmat’s only to commute with L1+L2+L3\leg{L}{1}+\leg{L}{2}+\leg{L}{3}, here the associativity of the products requires WB−1WAW_B^{-1}W_A only to commute with the triple products, and for a particular model that is a weaker demand. What makes the Yang–Baxter equation the right thing to impose is, again, that it mentions no model. (The recipe Lan(λ)=−Ran(−λ)\Lop{an}(\lambda)=-\Rmat_{an}(-\lambda) of (14.7) used R12=R21\Rmat_{12}=\Rmat_{21}; for a solution without that symmetry it is Lan(λ)=−Rna(−λ)\Lop{an}(\lambda)=-\Rmat_{na}(-\lambda).) For the Heisenberg chain nothing is left to choose: the check script finds that (14.4) fixes R\Rmat up to a scalar for one and two sites, at generic λ,μ\lambda,\mu.

Reading

14.12 The Yang–Baxter equation at order ℏ2\hbar^{2}, and more

YBE →\to CYBE. Give each line ii a parameter λi\lambda_i and write Rij=Rij(λi−λj)\Rmat_{ij}=\Rmat_{ij}(\lambda_i-\lambda_j), so (14.8) is R12R13R23=R23R13R12\Rmat_{12}\Rmat_{13}\Rmat_{23}=\Rmat_{23}\Rmat_{13}\Rmat_{12}. Suppose R=1+iℏr+ℏ2s+O(ℏ3)\Rmat=\Id+\ii\hbar\rmat+\hbar^{2}s+O(\hbar^{3}), in difference form. The terms of order ℏ\hbar are r12+r13+r23\rmat_{12}+\rmat_{13}+\rmat_{23} on both sides. Those of order ℏ2\hbar^{2} containing ss are also the same on both sides; the rest are (iℏ)2(\ii\hbar)^{2} times the products of two r\rmat’s in the order of the factors, so the difference of the two sides at order ℏ2\hbar^{2} is −ℏ2([r12,r13]+[r12,r23]+[r13,r23]):-\hbar^{2}\Big(\comm{\rmat_{12}}{\rmat_{13}}+\comm{\rmat_{12}}{\rmat_{23}} +\comm{\rmat_{13}}{\rmat_{23}}\Big): a solution of the Yang–Baxter equation has a classical r\rmat that solves the CYBE (14.12), in the per-leg form of §14.10. For R=λ−iℏP\Rmat=\lambda-\ii\hbar P, normalised by 1/λ1/\lambda (which does not change (14.8), the scalars being the same three on both sides), r=−P/λ\rmat=-\perm/\lambda is Proposition 14.10 at c=−1c=-1. This is the picture Figure 14.8 promised: the classical slide is the first correction to the quantum one proved; the check script repeats it with generic matrices.

The quantum determinant. At λ−μ=i\lambda-\mu=\ii, R(i)=i(1−P)\Rmat(\ii)=\ii(1-P) is 2i2\ii times the projector onto the antisymmetric part of C2⊗C2\C^{2}\otimes\C^{2}, and RTTRTT says that Ta(μ+i)Tb(μ)\mono_a(\mu+\ii)\mono_b(\mu) maps the antisymmetric line to itself. The operator by which it multiplies that line (central; for XXX a number) is the quantum determinant, qdet⁡T(μ)=A(μ+i)D(μ)−B(μ+i)C(μ),\operatorname{qdet}\mono(\mu)=A(\mu+\ii)D(\mu)-B(\mu+\ii)C(\mu), the counterpart of det⁡T=1\det\mono=1 for Toda (9.3). It commutes with every entry of T\mono (Korepin, Bogoliubov and Izergin, §VII.8, in their conventions). For the universal algebra defined by RTTRTT on a normalised formal series T(u)=1+T(1)u−1+⋯\mono(u)=\Id+\mono^{(1)}u^{-1}+\cdots, the Yangian Y(gl2)Y(\gllie_2), the coefficients of the quantum determinant generate the centre (Molev, ch. 1) cited; a spin chain gives a representation of it. For the Heisenberg chain it is the number [μ(μ+2i)]N[\mu(\mu+2\ii)]^{N} (the check script, N≤3N\le3). Lecture 15 meets its two factors as the eigenvalues of AA and DD on the reference state.

Other spins. For a spin ss at each site, Lan(λ)=λ+i2+i σa⋅Sn\Lop{an}(\lambda)=\lambda+\tfrac{\ii}{2}+\ii\,\boldsymbol\sigma_a\cdot\mathbf S_n satisfies (14.4) with the same R\Rmat, by a verification that uses only the commutation relations of S\mathbf S cited: Korepin, Bogoliubov and Izergin, §VIII.4, (4.6)–(4.8), in their conventions; the check script verifies s=1s=1. The transfer matrix still commutes with itself; but the regular point is lost, since L(0)\Lop{}(0) is no longer a permutation, and a local Hamiltonian needs a transfer matrix with a spin-ss auxiliary space (fusion: Kulish, Reshetikhin and Sklyanin 1981; not treated here).

Reading

14.13 Why the zero-curvature equation is not what is quantised

Lecture 8’s and Lecture 9’s readings derived the conservation of t(λ)\transfer(\lambda) from the zero-curvature equation (8.23), L˙n=Vn+1Ln−LnVn\dot L_n=V_{n+1}L_n-L_nV_n, and Lecture 9 proved, separately, the Sklyanin bracket. This lecture quantised only the second. The zero-curvature equation is the equation of motion under one Hamiltonian, and in quantum mechanics equations of motion are not postulated: they follow from [⋅,H]\comm{\cdot}{H} once HH is known. What has to be quantised is the kinematics, the brackets, and from them the commuting family and its Hamiltonians follow. The site-by-site transport survives as the monodromy (14.5), and the bracket as RTTRTT. The Heisenberg equation for T(λ)\mono(\lambda) under HH is then a consequence, and [H,t(λ)]=0\comm{H}{\transfer(\lambda)}=0 is Theorem 14.3, not an input. That is the reason the classical limit of §14.6 was taken in the bracket, not in the equation of motion.

Reading

14.14 The Toda pair, quantised

Lecture 12 set the Toda displacement beside Lecture 11’s θ′\theta' as a counterpart. For the Toda pair itself the link can be checked. In Proposition 12.1 the relative coordinate x=q1−q2x=q_1-q_2 and its momentum p=12(p1−p2)p=\tfrac12(p_1-p_2) have Hrel=p2+exH_{\mathrm{rel}}=p^{2}+e^{x}, and the relative speed is v=2pv=2p. Quantised, −ℏ2ψ′′+exψ=p2ψ.-\hbar^{2}\psi''+e^{x}\psi=p^{2}\psi . With ℏ=1\hbar=1 and k=pk=p, the solution that decays under the wall is ψ=K2ik(2ex/2)\psi=K_{2\ii k}(2e^{x/2}), the modified Bessel function: the chain rule turns this equation into Bessel’s. Far from the wall, x→−∞x\to-\infty, the small-argument form of KνK_\nu gives cited: DLMF §§10.25–10.27 ψ≈12[Γ(−2ik) eikx+Γ(2ik) e−ikx]:\psi\approx\tfrac12\big[\Gamma(-2\ii k)\,e^{\ii kx}+\Gamma(2\ii k)\,e^{-\ii kx}\big]: an incoming wave eikxe^{\ii kx}, reflected with the phase Γ(2ik)/Γ(−2ik)\Gamma(2\ii k)/\Gamma(-2\ii k). Since Γ(−2ik)=Γ(2ik)‾\Gamma(-2\ii k)=\overline{\Gamma(2\ii k)}, the stationary-phase argument of §11.3 puts the incoming packet on x=2kt+xwx=2kt+x_w and the outgoing one on x=−2kt+xwx=-2kt+x_w, with xw=2Re⁡ψ0(2ik)x_w=2\operatorname{Re}\psi_0(2\ii k) and ψ0=Γ′/Γ\psi_0=\Gamma'/\Gamma the digamma function: the packet moves as if reflected by a hard wall at xwx_w. In the language of §11.3, with k=p1−p2=2pk=p_1-p_2=2p the momentum difference, the reflected wave’s phase is θ(k)=2arg⁡Γ(ik)\theta(k)=2\arg\Gamma(\ii k) up to a constant, and θ′(k)=2Re⁡ψ0(ik)=xw\theta'(k)=2\operatorname{Re}\psi_0(\ii k)=x_w: the wall is θ′\theta'. Restoring ℏ\hbar shifts xx by 2ln⁡ℏ2\ln\hbar and replaces kk by p/ℏp/\hbar, so xw=2ln⁡ℏ+2Re⁡ψ0(2ip/ℏ) ⟶ 2ln⁡(2p)=2ln⁡v(ℏ→0),x_w=2\ln\hbar+2\operatorname{Re}\psi_0(2\ii p/\hbar)\ \longrightarrow\ 2\ln(2p)=2\ln v \qquad(\hbar\to0), by Stirling’s formula, Re⁡ψ0(iy)=ln⁡y+O(y−2)\operatorname{Re}\psi_0(\ii y)=\ln y+O(y^{-2}). The classical orbit of Proposition 12.1, x=ln⁡v2−v∣t∣+o(1)x=\ln v^{2}-v\abs t+o(1), has its effective wall exactly there, and that wall is what displaces each momentum line by ∓2ln⁡v\mp2\ln v. So the classical displacement is the ℏ→0\hbar\to0 limit of the quantum one (scripts/check-l16-algebra.py checks the Bessel equation, the asymptotic expansion and the effective-wall limit; the stationary-phase step is argued, not checked; Sutherland §4.2 does the same for the whole hyperbolic family) sketched.

Summary

  • The Lax operator Lan(λ)=λ+iPan=−Ran(−λ)\Lop{an}(\lambda)=\lambda+\ii P_{an} =-\Rmat_{an}(-\lambda) satisfies RLLRLL with R(λ)=λ−iP\Rmat(\lambda)=\lambda-\ii P (Theorem 14.1), which is Lecture 13’s Yang–Baxter equation.
  • The train (Theorem 14.2): the monodromy satisfies RTTRTT, and [t(λ),t(μ)]=0\comm{\transfer(\lambda)}{\transfer(\mu)}=0 (Theorem 14.3).
  • The classical limit (Theorem 14.4 and Corollary 14.6): the quantum Toda chain satisfies RLLRLL exactly with R=λ−iℏP\Rmat=\lambda-\ii\hbar P, and at order ℏ\hbar its RTTRTT relation is Lecture 9’s quadratic bracket with r=−P/λ\rmat=-\perm/\lambda.
  • Next: Lecture 15 finds the Heisenberg Hamiltonian in t(λ)\transfer(\lambda) and diagonalises the family by the algebraic Bethe ansatz.

Worked in the tutorial: Exercise Sheet 14.

References

E. K. Sklyanin, L. A. Takhtajan and L. D. Faddeev, Theor. Math. Phys. 40 (1979) 688 (the method); L. D. Faddeev, How algebraic Bethe ansatz works for integrable model, Les Houches lectures 1995, arXiv:hep-th/9605187, §§1–3; V. E. Korepin, N. M. Bogoliubov and A. G. Izergin, Quantum Inverse Scattering Method and Correlation Functions, Cambridge University Press (1993), §§VI.1–VI.5, §VII.8 (the quantum determinant), §VIII.4 (higher spin); M. C. Gutzwiller, Ann. Phys. 133 (1981) 304 and E. K. Sklyanin, Lect. Notes Phys. 226 (1985) 196 (the quantum Toda chain); A. Molev, Yangians and Classical Lie Algebras, American Mathematical Society (2007), Chapter 1 (the Yangian and its centre); O. Babelon, D. Bernard and M. Talon, Introduction to Classical Integrable Systems, Cambridge University Press (2003), Ch. 2–3 (the classical Yang–Baxter equation) and §4.7 (the open chain’s r\rmat). The conventions for the ordering of RTTRTT, of the monodromy and for the sign in R\Rmat differ between these sources; see §14.2.

The whole bibliography, lecture by lecture →