Quartic Hamiltonians, and higher Hamiltonians at next-to-leading order, for the affine Gaudin model
Abstract
In this work we will use a general procedure to construct higher local Hamiltonians for the affine Gaudin model. We focus on the first non-trivial example, the quartic Hamiltonians. We show by direct calculation that the quartic Hamiltonians commute amongst themselves and with the quadratic Hamiltonians which define the model.
We go on to introduce a certain next-to-leading-order semi-classical limit of the model. In this limit, we are able to write down the full hierarchy of higher local Hamiltonians and prove that they commute.
1 Introduction
Gaudin models represent interesting theories that find applications in several contexts in mathematical physics. In particular, they provide a general framework to study rich classes of classical and quantum integrable 1+1 dimensional integrable theories [Gau76, Gau14, Vic18, Lac18]. Recently, it has also been found that classical untwisted Gaudin models provide a dual description of such theories to the one given by 4D holomorphic Chern-Simons theories [Vic21, LV21].
In this paper we will focus on quantum Gaudin models. To define them, one needs a collection of distinct points on the Riemann sphere and a Kac-Moody algebra , which can be of finite or affine type. These models are defined by the quadratic Hamiltonians, namely
(1) |
in the (suitably completed, in the affine case) tensor product . Here are the generators of , defined at site , and is the invariant bilinear form on .
Finite-type quantum Gaudin models have been deeply studied. In particular, it is known that the quadratic Hamiltonians (1) are part of a larger family of mutually commuting operators called the Bethe or Gaudin subalgebra [Fre04, Tal04, Ryb06]. Moreover, these higher Hamiltonians can be diagonalized with an elegant form of Bethe Ansatz, where the eigenvector is the Schechtman-Varchenko vector and the eigenvalues are encoded as functions on a space of -opers, the Langlands dual algebra of [FFR94, MTV06, MTV07, MV02, MV04].
On the other hand, affine-type Gaudin models are still far from being fully understood.
One approach to affine Gaudin models was proposed in a recent paper [KL21]. In this work the authors propose an integrable model called generalized affine Gaudin model, which reproduces the usual model in a certain limit. Their construction is based on a new realization of the subalgebra through a new class of vertex operators, and fits affine Gaudin models into the general procedure first given in the seminal work [BLZ97].
Another approach, first proposed in the pioneering work [FF07], is to treat affine Gaudin models in closer analogy with their finite-type cousins. In particular, the spectrum of higher Hamiltonians should be described by suitable functions on a space of opers. Some further conjectures of how this might work, at least for the local Hamiltonians, were made in [LVY18], where it was conjectured that the eigenvalues of higher local Hamiltonians of the affine Gaudin models, as well as the Hamiltonians themselves, are given by hypergeometric-type integrals in the spectral plane, namely
(2) |
where lives is a (multi)set of indices given by the exponents of the algebra , is a certain multi-valued function defined by the data of the levels of the modules attached to the marked points , is a Pochhammer contour in the spectral plane around any two of these points (see e.g. fig. 1) and can be thought as the Hamiltonian density.

The key step in computing the higher Hamiltonians is to characterize these Hamiltonian densities, which are obtained by defining a suitable state for each given exponent . In order to do that, it is possible to exploit the general properties consistent Hamiltonians must obey: they have to commute with the generators of the algebra defining the model as well as amongst themselves (here denotes the underlying finite algebra). As we will see in section 4.7, this is equivalent to the following requirements
(3) |
where represents the diagonal action of the positive modes of the generators of the algebra and the zero mode is intended in the vertex algebra sense (see section 2.4). We give a precise definition of translation and twisted derivative of a state in sections 2.4 and 4.4.
The general expectation is that there exists a state for every exponent of the affine algebra , and that it takes the following form
(4) |
where and is a certain symmetric invariant tensor of . This particular structure is justified by the semi-classical counterpart of these models, which have been thoroughly studied [Eva+99, Eva01, LMV17]; in particular, the precise choice of symmetric invariant tensors needed to ensure that the Hamiltonians Poisson-commute is well understood and related to Drinfel’d-Sokolov reduction [Eva01, DS85]. In the very simplest cases, including the cubic Hamiltonian in type , there are no quantum corrections needed [LVY20]. The first case in which quantum corrections are present appeared in the recent paper [KLT22]: it is the example of the quartic Hamiltonian in type . In this paper, we consider this example in detail. In particular, we carry out the full computation to show that the resulting quartic Hamiltonians commute amongst themselves (see theorem 4.5.1 and proposition 4.7.2).
This direct calculation reveals the following fact: it turns out that any density obeying the first of the two conditions in eq. 3 (i.e. the one which ensures the Hamiltonian commutes with the generators of the affine algebra), automatically also obeys the second condition (which is needed for the Hamiltonians to commute amongst themselves) at least for , as shown in corollary 4.5.1 below. We should stress that this property is, for the moment, highly non-obvious and suggestive; it would be good to get a more systematic understanding of why it should be true.
Already in this case of the exponent , i.e. of quartic Hamiltonians, the direct computations needed are very lengthy. This is especially true of the computations needed to show the mutual commutativity of the quartic Hamilonians. For higher exponents , direct calculations become computationally difficult even with the aid of computer algebra, but we are able to prove a result for all exponents by truncating to the next-to-leading order in . To introduce the dependence from the formal parameter , we perform a re-scaling of the generators, namely , , in such a way that every time we perform a commutation, we introduce a factor of . This procedure allows us to identify different quantum corrections by their dependence. In this spirit, we will then work modulo terms of order . We show that, modulo such terms, the Hamiltonian density for each exponent , , of , takes the following form
and prove that the resulting Hamiltonians commute up to and including terms of order .
***
The paper is organized as follows.
In section 2 we recall the main ideas behind the theory of Kac-Moody algebras, their local completion and the concept of vacuum Verma module. We also recall the definition and the main features of vertex algebras.
In section 3 we define the algebra of observables of the Gaudin model and we recall the definition of invariant tensor of an algebra , focusing in particular on the case of .
In section 4 we define the meromorphic states , which represent the building blocks for the states we want to define. We also present their commutation relations and we describe the gradation they give rise, which will be of fundamental importance to prove the main results of the paper. After that, we characterize the space of states that vanish under the diagonal action of zero modes: this is the first step to define higher Gaudin Hamiltonians. At this point we will present all quantum corrections to the quartic state and we will explain why they appear in this specific example and not in the other known ones. In the following subsection, we restrict the number of possible states by asking they be singular up to twisted derivative under the diagonal action of positive modes. We focus on the quadratic and quartic states. Afterwards, we prove the main result of the paper, where we also ask that the 0th product (in the vertex algebra sense) of these states vanishes modulo twisted derivatives and translates. We will show that there is one quartic state, up to re-scaling and modulo the addition of translates and twisted derivatives, that satisfies this last requirement. At this point, having a precise definition of the quartic state, we describe the general construction of the higher Hamiltonians, in the spirit of [LVY18]. We will explain in detail why we ask for these properties and why they are important at the level of the quantum Hamiltonians.
In section 5 we try to solve the same problem from a different point of view. Instead of focusing on one specific Hamiltonian, trying to work out its explicit definition, we want to characterize all Hamiltonians at arbitrary , but working at sub-leading order. In order to do that we introduce the formal parameter by making a particular re-scaling of the generators of the algebra. We will then prove similar theorems to those of the previous sections, working modulo terms at order .
Acknowledgements. The authors would like to thank Sylvain Lacroix for helpful discussions and for providing many detailed comments. The research of Charles Young is supported by the Leverhulme Trust, Research Project Grant number RPG-2021-092.
2 Vacuum Verma modules for
2.1 Loop realization of
We define the loop algebra as the algebra of Laurent polynomials in a formal variable with coefficient in the finite-dimensional Lie algebra . The Lie brackets on this algebra are given by
(5) |
where and are arbitrary Laurent polynomials in .
Let be the canonically normalized symmetric invariant bilinear form on . It is given by taking the trace in the defining representation:
(6) |
It is possible to extend the loop algebra by a one-dimensional central element ,
(7) |
This extension is called the affine Lie algebra , whose commutation relations are
(8) | |||
(9) |
We shall use the notation
(10) |
The commutation relations can be equivalently written as
(11) |
We can add to this algebra a one-dimensional derivation , such that and , for all and . It is possible to show that this algebra is isomorphic to the Kac-Moody algebra over of type , see e.g. [Kac90, ch. 7],
(12) |
2.2 as a Kac-Moody algebra
The Cartan matrix for the Kac-Moody algebra of type is defined as . The Cartan decomposition is given by . The Chevalley-Serre generators are , while and are respectively a basis for the Cartan subalgebra of simple coroots of and a basis for the dual Cartan subalgebra of simple roots of . The latter are related via the canonical pairing between the Cartan algebra and its dual,
(13) |
The fundamental commutation relations in are
(14) |
where and , together with the Serre relations
(15) |
The Kac-Moody algebra has a central element , which spans a one-dimensional centre. A basis for the Cartan subalgebra is given by the coroots together with the derivation element , which by definition satisfies
(16) |
If we remove the 0th row and column from , we obtain the Cartan matrix corresponding to the finite dimensional Lie algebra . This subalgebra of is generated by , and .
2.3 Local completion and vacuum Verma modules
For any , let us define as the quotient of the universal enveloping algebra of by the two-sided ideal generated by . For each , let us introduce the left ideal . The inverse limit
(17) |
is a complete topological algebra, called the local completion of at level . With this definition, the elements of are possibly infinite sums of the type of elements which do truncate to finite sums when working modulo any .
A module over is said to be smooth if, for all and all , for sufficiently large . A module has level if acts as zero on . Any smooth module of level over is also a module over the completion .
We can identify the subalgebra of positive modes and introduce the one-dimensional representation defined by
(18) |
We define , the vacuum Verma module at level , as the induced smooth -module
(19) |
This vector space is spanned by monomials of the form , with and strictly negative mode numbers . We call these vectors states.
Let us denote by the derivation on defined by and . By setting for any , one can interpret as a translation operator : the reason for this identification will be clear in the next section.
2.4 Vertex algebra structure
As we now recall, the vacuum Verma module defined in the previous section has the structure of vertex algebra. Namely, we have the state-field map , which for every state associates a formal power series in the variable ,
(20) |
where is the th mode of . By definition, if for some , then for all , i.e.
(21) |
The fields for all other states can be obtained with the following properties
(22) |
where we have introduced the normal ordered product between fields
(23) |
In fact, any state can be written as , and using eqs. 22 and 23 we can always explicitly compute . Summarising what we have introduced so far, we have
-
•
a space of states ,
-
•
a vacuum vector ,
-
•
a translation operator ,
-
•
the state field map as in eq. 20.
This structure, together with some additional properties (see e.g. [FB01, §2.4.4]), defines a vertex algebra on .
3 Construction of higher Hamiltonians
3.1 The algebra of observables
Let us introduce a set of complex numbers , where and for all . Consider the following tensor product of vacuum Verma modules
(24) |
This space can be interpreted as a module over the direct sum of copies of . Let us denote by the copy of in the th direct summand.
Let us denote by the one-dimensional vacuum representation of the “positive modes” Lie subalgebra , defined by and for all , and . Therefore is the induced -module, namely
(25) |
Repeating similar arguments to those of the previous sections, we can define as the quotient of by the two-sided ideal generated by for all . We have the isomorphism
(26) |
Thanks to this fact, can be presented as
(27) |
where is acting as the th tensor factor. Again, we can introduce the left ideals generated by for all , and . Let be the inverse limit. This space is a complete topological algebra and
(28) |
where denotes the completed tensor product. This space is called the algebra of observables of the Gaudin model.
The tensor product is again a vertex algebra. The state-field map is just as in section 2.4 but decorated with the extra index (i).
Let us introduce the map , which is the diagonal embedding of into , defined as
(29) |
It extends to an embedding . It is easy to check that
(30) |
where is the usual Killing form as in eq. 6.
Therefore descends to an embedding of the quotients , where , and hence of their completions
(31) |
3.2 Invariant tensors on
Let be a basis of , and let be its dual basis with respect to the non-degenerate Killing form of eq. 6. Let denote the structure constants, so that
(32) |
Here and in what follows we employ the summation convention on Lie algebra indices. Thanks to the non-degeneracy of the bilinear form, we may suppose our basis is chosen in such a way that
(33) |
By doing this, we no longer have to distinguish between upper and lower indices. The structure constants are then
(34) |
where is the usual Levi-Civita symbol.
(Concretely, in the defining representation we have
(35) |
Recall that for any finite-dimensional Lie algebra , a tensor is invariant if
(36) |
or equivalently, if its components satisfy
(37) |
where the indices take values from to . In our case of , the ring of invariant tensors is generated by and . We shall need the following syzygy relations between them:
(38) |
Note in particular the last of these, which will play a crucial role in the explicit calculations of the next sections. It can also be generalized to higher rank tensors (see e.g. [Itō93, §369 F]).
4 Quartic Hamiltonian
4.1 Meromorphic states
Let us introduce a set of points in the complex plane, chosen to be pairwise distinct, whenever . For any element we introduce the -valued meromorphic functions
(39) |
We are allowed to take derivatives of such functions, which will be denoted by or, in general, for each ,
(40) |
Considering two of these functions with different spectral parameters, we get the following commutation relations
(41) |
By taking the limit , we get the commutation relations for the same spectral parameter, namely
(42) |
We see that these , for and , span a Lie algebra of -valued meromorphic functions of with poles at the marked points.
It is helpful to be able to treat this as an abstract Lie algebra. Thus, let denote the Lie algebra over with basis consisting of and , for , and with the non-vanishing Lie brackets given by
(43) |
Let denote subalgebra generated by for , and , and let
(44) |
denote the module over induced from the trivial one-dimensional module over .
We call the the space of meromorphic states. It is again a vertex algebra, with the state-field map as in section 2.4 but decorated with extra indices. For each , one has the homomorphism of Lie algebras given by evaluating at . It gives rise to a map of vertex algebras.
There is a bi-gradation of in which (for any ) has weight and has weight . This yields a bi-gradation of
(45) |
For each , let denote the subspace of grade . We call elements of homogeneous meromorphic states of degree .
4.2 Diagonal action of the zero modes of
There is an evident diagonal action of the Lie algebra on the -module , defined in the same way as the action on in eq. 29. In particular, for any , the zero modes stabilize each subspace , namely
(46) |
An important fact is that every state in properly contracted with an -invariant tensor vanishes under the diagonal action of the zero modes. This follows directly from the defining property of invariant tensors in eq. 37. Let denote with the invariant subspace. We can characterize this space for small :
-
•
for , .
-
•
for , . Indeed, elements of are of the form . Such an element is in if and only if are the components of an -invariant tensor of rank 1. But there are no nonzero such tensors.
-
•
for , has dimension 1 and it is spanned by the state
(47) -
•
for , has dimension 2 and it is spanned by the states
(48) and
(49) -
•
for , has dimension 14. Below, we will make use of the following explicit choice of basis:
(50)
Note that to write these terms we have to choose an ordering prescription. In this work we sort level first in ascending order from left to right and after that, for a given level, we sort derivatives in descending order from left to right. For example .
4.3 Top terms
We can see from the above construction that in the case and , there is a particular state, that we will call top term, which is the state in that contains exactly generators:
(51) |
We do not have such state for , because we can always use the commutation relations to reduce the number of generators, as shown in eq. 48. This pattern continues, as we now describe.
Notice that the universal enveloping algebra has an increasing filtration
(52) |
in which the generators count as and the generators count as , cf. the commutation relations of in eq. 43. For example , and as well. It gives rise to a corresponding filtration, , on the space of meromorphic states.
Observe that if then . We see that
(53) |
for some tensor , where the brackets around the indices denote the operation of symmetrization,
(54) |
(and we may symmetrize without loss of generality because the non-symmetric pieces fall into , as for example in eq. 48). Let us call the top term of the state .
If this state is -invariant, , then is a symmetric invariant tensor. Nonzero such tensors exist only in even degrees, and up to rescaling they are, explicitly,
In what follows, our interest is in meromorphic states that have nonzero top term (in other words states whose principal symbol has maximal degree) and that are -invariant.
4.4 Singular vectors up to twisted derivative
Let us define the twisted derivative operator of degree with respect to the spectral parameter ,
(55) |
Note that this operator sends in the bigradation we introduced above.
We will say that a vector is singular up to twisted derivatives if for all we have
(56) |
for all non-negative modes , . This defines a subspace
(57) |
of vectors singular up to twisted derivatives.
Proposition 4.4.1.
The space of singular vectors is spanned by the quadratic state defined in (47).
Proof.
We need to show that
(58) |
for some meromorphic states , for all and . For there is nothing to check since identically, by the definition of . It is enough to check the action of the first modes , since any higher modes can be expressed in terms of their brackets, i.e. etc. From direct calculations we get that
(59) |
where
(60) |
∎
More non-trivially, for we have the following result.
Proposition 4.4.2.
The space of singular vectors is of dimension 7. A choice of basis is given by the state
(61) |
together with the double translate state
(62) |
and the following twisted derivative states
(63) | |||
Proof.
Let . We may write it in our basis (50),
(64) |
for some coefficients with , and ask what conditions the requirement of being singular up to twisted derivatives, (56), places on these coefficients. It is enough to demand that
(65) |
for some meromorphic states , for all and . For zero modes there is nothing to check since exactly, by definition of . It is then enough to check the action of first modes, , since any higher modes can be expressed in terms of their brackets, etc. So we are to check under what conditions
(66) |
for some . By direct calculation, one finds that solutions exist precisely if the coefficient obey the relations
(67) |
When they do obey these relations, the required functions are given by
(68) |
where
The basis reported in the proposition can be obtained by the one defined by the restrictions (67) by a change of basis. ∎
The proposition above is in agreement with the calculation of the quartic Hamiltonian density (the analogue of our ), recently presented in [KLT22]. In the present conventions, the latter is given by
(69) |
and it does333To match conventions, note that for us (70) and in [KLT22] the tensor called is given by (71) We thank Sylvain Lacroix for clarifying discussions on this point. indeed lie in the space .
4.5 Hamiltonian densities
Now, to state the main result of the paper, we need two reintroduce rational functions of two different spectral parameters, and , cf. eq. 41 and eq. 42.
Recall the Lie algebra over from section 4.1. Let be the Lie algebra with generators , , and for , and and commutation relations
(72) |
together with eq. 43 for the generators with parameter and the analogue with parameter . This Lie algebra and its modules are defined over the ground ring of polynomials in . We have the vertex algebra defined analogously to eq. 44 and the two obvious embedding maps of vertex algebras , which we write as and .
Moreover, there is a natural notion of “expanding around ”. Namely, there is a homomorphism of Lie algebras over defined by
(73) |
which is motivated by considering the Taylor expansion of the function from eq. 39. This gives rise to a map . We say a state is regular at modulo translates if there exists such that the image of under this map has no singularities in .
Theorem 4.5.1.
The elements and obey the relations
(74a) | ||||
(74b) | ||||
(74c) | ||||
(74d) |
where and are elements of . Moreover, , , are regular at modulo translates.
Proof.
The two statements of the theorem follow from direct calculations. In particular, when , we get
(75) |
We have computed , , and explicitly, with the aid of the computer algebra system FORM [Ver13, Kui+13]. The expressions for and , are given in appendix A. The expressions for and are extremely lengthy (more that 500 terms in total), and we do not reproduce them here.
Once the expression of is known, i.e. the functions and are found, it can be shown that the theorem is automatically satisfied for the product . This comes from the property of the th product between two states of a vertex algebra, namely
(76) |
Therefore, by swapping two states in a 0th product, we obtain and a series of terms which are nothing but translates and therefore can be absorbed in the definition of .
To prove the second part of the theorem one can expand according to eq. 73 and the result follows from direct calculation. ∎
Having established this statement for the particular choice of quartic density , we automatically get the following property for any element of . It is a slightly weaker property, because the condition on the twisted derivative terms on the right hand side is less rigid. As we shall see in section 4.7 below, it is sufficient for defining consistent Hamiltonians.
Corollary 4.5.1.
For any element , one has
(77) | ||||
(78) | ||||
(79) |
where and are elements of . Moreover, , , are regular at modulo translates.
Proof.
We already know from theorem 4.5.1 that there exists an element, , satisfying these relations. But we saw in proposition 4.4.2 that every element of is proportional to up to the addition of certain translates and twisted derivatives.
It follows from the property (76) that if we add to any translate then the statement of theorem 4.5.1 still holds, the only difference being a re-definition of the states . And it is evident that, if we add to any linear combination of the twisted derivatives in eq. 63 then the resulting vector still obeys the weaker relations given above. (One might worry about introducing singularities at , but note that for any meromorphic states and , the product is regular at , as is manifest if we expand about in the spectral plane before taking the vertex-algebra product: .) ∎
4.6 Gaudin Hamiltonian
Let us define the following state at non-critical level, i.e. ,
(80) |
where is now the image in of the density defined in (47) and where
(81) |
the term in the brackets being the Segal-Sugawara vector at site .
It is possible to show (see [LVY20]), that the operator is the image in of
(82) |
where
(83) |
is the th copy of the quadratic Casimir operator of in and are the Hamiltonians in (1).
Theorem 4.6.1.
Given the images in of the densities , , we have
(84) |
with and being the images in of the meromorphic states in theorem 4.5.1.
Proof.
The result follows from direct calculations, using the definitions of , , , in (75) and appendix A, respectively. ∎
As we will see in the next section, this requirement is sufficient to ensure the commutativity of local Hamiltonians, arising from the densities and , with the usual quadratic Gaudin Hamiltonians which define the model.
4.7 Commuting Hamiltonians
In this section, we will simply recall the ideas presented in [LVY18]. Consider two states and their formal zero modes . It is possible to show the following vertex algebra identity (see e.g. [FB01])
(85) |
This means that if one is able to find a family of operators whose 0th product vanishes (or that can be expressed as a translation, since by definition), then their formal zero modes form a commutative subalgebra of the algebra of observables .
The meromorphic function which is obtained by acting with on the module , i.e. setting the central elements to numbers,
(86) |
has a special role and it is called the twist function of the model. Let us define also
(87) |
The function is multi-valued. It becomes single-valued on a certain multi-sheeted cover of . Let be any closed contour in this cover. For example, could be the lift to this cover of a Pochhammer contour in around any two of the marked points. Then , for any integer , is single-valued along , and one can introduce the integral
(88) |
which has the fundamental property that, for any meromorphic function which is non-singular along
(89) |
Let us now define the following object in ,
(90) |
for , where are now the images in of the densities we have defined in the previous section.
Proposition 4.7.1.
The operators commute amongst themselves, with the generators of , and with the quadratic Hamiltonians .
Proof.
We can use eq. 85 to compute
(91) |
but we know from theorem 4.5.1 that the 0th product between those states can be expressed as a sum of twisted derivatives and translations. It is now straightforward to check that the result of the commutator is zero: on one side because for every state , on the other because of the property (89).
To prove the second statement, we recall that and the general property that for any two states of a vertex algebra we have
(92) |
where
(93) |
Note that this represents a generalization of eq. 85, see e.g. [VY17, FB01]. At this point we can consider the generators with and we get
(94) |
It is now straightforward to check that the result is zero, by using property (56), described in the relevant cases in propositions 4.4.1 and 4.4.2, and the property (89). ∎
4.7.1 Fourier modes
Even though the operators we have just defined have all the right characteristics to be well-defined Hamiltonians as pointed out in proposition 4.7.1, there is one last subtlety about these objects, related to the fact that we want their action on highest weight modules to be diagonalisable. In fact, considering such that and setting , induces a -gradation on the product of vacuum Verma modules, called the homogeneous gradation. Therefore if with degree , the degree of its modes is . The objects we have constructed , by definition, have , therefore .
This shows that in the homogeneous gradation these operators have degree : this means that if we consider a module over which has a trivial subspace of grade for large , then the operator has a non-zero eigenvalue.
A way to overcome this issue is to consider the notion of Fourier mode of the state : they have the property that we are looking for, namely . Additionally, they satisfy a similar relation to (85),
(95) |
with . One has for and it is possible to show that the following recursive formula holds:
(96) |
where is the Taylor series in given by
(97) |
and where the coefficients are defined by the requirement that and are the expansions, for and respectively, of the function
(98) |
The first relevant examples are
(99a) | ||||
(99b) |
where . These formulae are the Fourier-analog of the normal ordered product formula eq. 23, and they allow one to compute by recursion the Fourier modes of a general state .
Property (95) means that if the vertex algebra 0th product of and vanishes their Fourier zero-modes generate a commutative subalgebra, in homogeneous degree zero, of . We let
(100) |
for . By the same logic as for proposition 4.7.1, we have the following.
Proposition 4.7.2.
The operators have homogeneous degree 0 and they commute amongst themselves, with the generators of , and with the quadratic Hamiltonians .
5 Higher local Hamiltonians to sub-leading order
In the previous section, we have shown that it is possible to define quartic local Hamiltonians which commute among themselves and with the quadratic ones, together with the generators of . Following the same steps, one could in principle try to construct the Hamiltonians for every exponent of . However, the direct calculation (already lengthy in the case of , as we noted above) becomes computationally very demanding. What we shall do instead is work to next-to-leading order in a certain semiclassical limit, which will at least give a strong consistency check on the existence of the expected Hamiltonian densities.
Thus, let us introduce a formal parameter and work over . So in particular, all vector spaces above are now to be regarded as modules over . We consider the following rescaled generators:
(101) |
With this re-scaling the commutation relations become
(102) |
At this point we can identify the various quantum corrections by their dependence and work grade by grade. We shall work at next-to-leading order, i.e. the next order beyond the usual semi-classical calculation of Poisson brackets. Thus, we consider the densities of Hamiltonians up to and including the leading quantum corrections at order , and we compute commutators up to and including terms of order .
Remark.
It is worth remarking that the classical limit, section 5, that we take is not quite the standard one which recovers the usual classical Gaudin model (cf. [KLT22] for a very complete discussion of that limit), because for us the central charges remain in the limit. From our present perspective this is simply for computational convenience – this limit produces the simplest possible non-trivial check, and had we rescaled the central charges there would be more potential quantum correction terms already at next-to-leading order. But it might be interesting to consider this classical limit in its own right.
Having introduced the formal parameter , there is a gradation on the enveloping algebras in which and have grade one and has grade zero. Recall that denotes the space of homogeneous meromorphic states of degree , eq. 45. Let now denote the subspace consisting of states that are also of grade in this new gradation (i.e. which are sums of terms having exactly factors of or ).
Proposition 5.0.1.
Modulo terms of order there is, up to rescaling, exactly one state , , such that, for all and ,
(103) |
Explicitly, modulo terms in ,
(104) |
Proof.
Given the basis for , we need to show that there exist a function such that
(105) |
For , this is always true thanks to the invariance of the tensor (37). For the same reasons explained in the previous section, the only relevant check that one needs to make is the one for . From direct calculation, we get
(106) |
Note that this result is in accordance with the exact ones obtained for the quadratic () and the leading order of quartic () states, cf. eqs. 60 and 68. ∎
(Observe that this is consistent with proposition 4.4.2 because the vectors in eqs. 62 and 63 all come with factors of in the limit.)
We can now state the following theorem
Theorem 5.0.1.
Let be as in eq. 104 above, for all odd . We have
(107) |
where are given by
(108) |
(109) |
where
(110) |
Moreover, is a regular function for modulo translates and modulo terms proportional to .
Proof.
The 0th mode of can be inferred from a purely combinatorial reasoning. Let us start with the top term of (104). We know that computing the 0th mode, the number of generators in any term we get does not change, but the result will be a state of total depth and therefore we know there must be at least one generator with a positive mode. We can also use the fact that we are working at leading order in , therefore we could get at least one , one or a term with two , every other term will be . The only thing to fix is the combinatorial factor describing the number of possible ways to write such terms. The result is
With similar arguments we can compute the 0th mode of the correction term of (104), the result reads
where . At this point, acting with what we have obtained on and using repeatedly the commutation relations (41), we obtain eq. 107. ∎
Appendix A Full expression for and
The explicit expressions for and in obtained by direct calculations are
References
- [BLZ97] V.. Bazhanov, S.. Lukyanov and Al.. Zamolodchikov “Integrable structure of conformal field theory. 2. Q operator and DDV equation” In Commun. Math. Phys. 190, 1997, pp. 247–278 DOI: 10.1007/s002200050240
- [DS85] V.. Drinfel’d and V.. Sokolov “Lie algebras and equations of Korteweg-de Vries type” In Journal of Soviet Mathematics 30.2, 1985, pp. 1975–2036 DOI: 10.1007/BF02105860
- [Eva+99] J.M. Evans, M. Hassan, N.J. MacKay and A.J. Mountain “Local conserved charges in principal chiral models” In Nuclear Physics B 561.3 Elsevier BV, 1999, pp. 385–412 DOI: 10.1016/s0550-3213(99)00489-7
- [Eva01] J.. Evans “Integrable sigma-models and Drinfel’d–Sokolov hierarchies” In Nuclear Physics B 608.3 Elsevier BV, 2001, pp. 591–609 DOI: 10.1016/s0550-3213(01)00157-2
- [FB01] E. Frenkel and D. Ben-Zvi “Vertex Algebras and Algebraic Curves”, Mathematical surveys and monographs American Mathematical Society, 2001 URL: https://books.google.co.uk/books?id=mRgqy1f1QkIC
- [FF07] B. Feigin and E. Frenkel “Quantization of soliton systems and Langlands duality” arXiv, 2007 DOI: 10.48550/ARXIV.0705.2486
- [FFR94] B. Feigin, E. Frenkel and N. Reshetikhin “Gaudin model, Bethe Ansatz and critical level” In Communications in Mathematical Physics 166.1, 1994, pp. 27–62 DOI: 10.1007/BF02099300
- [Fre04] E. Frenkel “Gaudin model and opers” In Workshop on Infinite Dimensional Algebras and Quantum Integrable Systems, 2004 arXiv:math/0407524
- [Gau14] M. Gaudin “The Bethe Wavefunction” Cambridge University Press, 2014 DOI: 10.1017/CBO9781107053885
- [Gau76] M. Gaudin “Diagonalisation d’une classe d’hamiltoniens de spin” In Journal De Physique 37, 1976, pp. 1087–1098
- [Itō93] Kiyosi Itō “Encyclopedic dictionary of mathematics” Cambridge, Mass: MIT Press, 1993
- [Kac90] V.. Kac “Infinite-Dimensional Lie Algebras” Cambridge University Press, 1990 DOI: 10.1017/CBO9780511626234
- [KL21] G.. Kotousov and S.. Lukyanov “ODE/IQFT correspondence for the generalized affine (2) Gaudin model” In JHEP 09, 2021, pp. 201 DOI: 10.1007/JHEP09(2021)201
- [KLT22] G.. Kotousov, S. Lacroix and J. Teschner “Integrable sigma models at RG fixed points: quantisation as affine Gaudin models” arXiv, 2022 DOI: 10.48550/arxiv.2204.06554
- [Kui+13] J. Kuipers, T. Ueda, J… Vermaseren and J. Vollinga “FORM version 4.0” In Comput. Phys. Commun. 184.5, 2013, pp. 1453–1467 DOI: 10.1016/j.cpc.2012.12.028
- [Lac18] S. Lacroix “Integrable models with twist function and affine Gaudin models”, 2018 arXiv:1809.06811 [hep-th]
- [LMV17] S. Lacroix, M. Magro and B. Vicedo “Local charges in involution and hierarchies in integrable sigma-models” In Journal of High Energy Physics 2017.9 Springer ScienceBusiness Media LLC, 2017 DOI: 10.1007/jhep09(2017)117
- [LV21] S. Lacroix and B. Vicedo “Integrable E-Models, 4d Chern-Simons Theory and Affine Gaudin Models. I. Lagrangian Aspects” In Symmetry, Integrability and Geometry: Methods and Applications SIGMA (Symmetry, IntegrabilityGeometry: MethodsApplication), 2021 DOI: 10.3842/sigma.2021.058
- [LVY18] S. Lacroix, B. Vicedo and C.A.S. Young “Affine Gaudin models and hypergeometric functions on affine opers” In Advances in Mathematics 350, 2018 DOI: 10.1016/j.aim.2019.04.032
- [LVY20] S. Lacroix, B. Vicedo and C.A.S. Young “Cubic hypergeometric integrals of motion in affine Gaudin models” In Adv. Theor. Math. Phys. 24.1, 2020, pp. 155–187 DOI: 10.4310/ATMP.2020.v24.n1.a5
- [MTV06] E. Mukhin, V. Tarasov and A. Varchenko “Bethe eigenvectors of higher transfer matrices” In Journal of Statistical Mechanics: Theory and Experiment 2006.8, 2006, pp. 08002 DOI: 10.1088/1742-5468/2006/08/P08002
- [MTV07] E. Mukhin, V. Tarasov and A. Varchenko “Schubert calculus and representations of general linear group” arXiv, 2007 DOI: 10.48550/ARXIV.0711.4079
- [MV02] E. Mukhin and A. Varchenko “Critical points of master functions and flag varieties” arXiv, 2002 DOI: 10.48550/ARXIV.MATH/0209017
- [MV04] E. Mukhin and A. Varchenko “Norm of a Bethe Vector and the Hessian of the Master Function” arXiv, 2004 DOI: 10.48550/ARXIV.MATH/0402349
- [Ryb06] L. Rybnikov “Uniqueness of higher Gaudin Hamiltonians” In Reports on Mathematical Physics 61, 2006 DOI: 10.1016/S0034-4877(08)80013-4
- [Tal04] D. Talalaev “Quantization of the Gaudin system”, 2004 arXiv:hep-th/0404153
- [Ver13] J. Vermaseren “FORM”, 2013 URL: https://www.nikhef.nl/~form/maindir/documentation/documentation.html
- [Vic18] B. Vicedo “On Integrable Field Theories as Dihedral Affine Gaudin Models” In International Mathematics Research Notices 2020.15, 2018, pp. 4513–4601 DOI: 10.1093/imrn/rny128
- [Vic21] B. Vicedo “4D Chern–Simons theory and affine Gaudin models” In Letters in Mathematical Physics 111, 2021 DOI: 10.1007/s11005-021-01354-9
- [VY17] B. Vicedo and C.A.S. Young “Vertex Lie algebras and cyclotomic coinvariants” In Communications in Contemporary Mathematics 19.02, 2017, pp. 1650015 DOI: 10.1142/S0219199716500152