Hopf Algebra Symmetries of an Integrable Hamiltonian for Anyonic Pairing

Since the advent of Drinfel’d’s double construction, Hopf algebraic structures have been a centrepiece for many developments in the theory and analysis of integrable quantum systems. An integrable anyonic pairing Hamiltonian will be shown to admit Hopf algebra symmetries for particular values of its coupling parameters. While the integrable structure of the model relates to the well-known six-vertex solution of the Yang–Baxter equation, the Hopf algebra symmetries are not in terms of the quantum algebra Uq(sl(2)). Rather, they are associated with the Drinfel’d doubles of dihedral group algebras D(Dn).


Introduction
Integrable quantum systems which admit exact solutions are central in advancing understanding of many-body systems.Classic examples are provided by the Heisenberg spin chain [1], the Bose [2] and Fermi [3] gases with delta-function interactions, the Bardeen-Cooper-Schrieffer pairing Hamiltonian with uniform scattering interactions [4], and the Hubbard model in one dimension [5].With the development of the Quantum Inverse Scattering Method [6] as a systematic prescription for constructing integrable quantum systems through the Yang-Baxter equation [3,7,8], and solving them through the algebraic Bethe ansatz, it subsequently emerged that Hopf algebraic structures are fundamental in quantum integrability.The works of Jimbo [9] and Drinfel'd [10] were instrumental in formulating the notion of quantum algebras U q (g), deformations of the universal enveloping algebras of a Lie algebra g, which have the structure of a quasi-triangular Hopf algebra.The significance of the quasi-triangular structure is that it affords an algebraic solution of the Yang-Baxter equation.Matrix solutions of the Yang-Baxter equation are then generated through representations of these algebras.The simplest example of the two-dimensional loop representation of the untwisted affine quantum algebra U q (sl(2) (1) ) leads to the six-vertex model solution of the Yang-Baxter equation, which establishes integrability of the anisotropic (XXZ) Heisenberg chain.The precise form of six-vertex solution obtained depends on the choice of gradation for U q (sl(2) (1) ).The principal gradation leads to the symmetric solution, while the homogeneous gradation leads to an asymmetric solution [11].Only in the latter case is the solution invariant with respect to the action of the non-affine subalgebra U q (sl(2)).
The work of Drinfel'd [10] also provides a means to construct a quasi-triangular Hopf algebra from any Hopf algebra and the dual algebra, through a procedure known as the double construction.The double construction applied to finite group algebras [12] yields a framework in which to develop anyonic models that lead to notions of topological quantum computation [13].In a series of works [14][15][16], solutions of the Yang-Baxter associated with Drinfel'd doubles of dihedral group algebras, denoted D(D n ), have been studied.In particular, it was found that two-dimensional representations of these algebras belong to the aforementioned six-vertex model solution in the symmetric case.The symmetric solution was employed in [17] to construct an integrable anyonic pairing Hamiltonian, which generalises the pairing Hamiltonian with uniform scattering interactions solved by Richardson [4].Below, this integrable anyonic pairing Hamiltonian will be shown to admit Hopf algebra symmetries given by D(D n ) for particular values of the coupling parameters.

The Integrable Hamiltonian for Anyonic Pairing
Consider a general anyonic pairing Hamiltonian of the reduced Bardeen-Cooper-Schrieffer form, which acts on a Hilbert space H of dimension 4 L , given by Above, {ε j : j = 1, ..., L} represent single-particle energy levels (two-fold denegerate labelled by ±) and G kl are the pairing interaction coupling parameters of the model.For q = exp(iβ), β ∈ R the operators {a j± , a † j± : j = 1, ..., L} satisfy the relations and those relations obtained by taking Hermitian conjugates.Throughout, I is used to denote an identity operator.These types of anyonic operators are considered as q-deformations of fermionic operators, with the usual fermionic commutation relations recovered in the limit q → 1.The anyonic creation and annihilation operators may be realised in terms of the canonical fermionic operators {c j± , c † j± : j = 1, ..., L} through a generalised Jordan-Wigner transformation where n jσ = c † jσ c jσ .As with the more familiar fermionic pairing Hamiltonians, one of the notable features of Equation ( 1) is the blocking effect.For any unpaired anyon at level j, the action of the pairing interaction is zero since only paired anyons interact.This means that the Hilbert space can be decoupled into a product of paired and unpaired anyonic states in which the action of the Hamiltonian on the space for the unpaired anyons is automatically diagonal in the natural basis.In view of this property, the pair number operator commutes with Equation ( 1) and thus provides a good quantum number.Below, M will be used to denote the eigenvalues of the pair number operator.
In [17] it was shown that, for a suitable restriction on the coupling parameters, the Hamiltonian is integrable in the sense of the Quantum Inverse Scattering Method and admits an exact solution derived through the algebraic Bethe ansatz.To characterise the integrable manifold of the coupling parameter space, the set of parameters {α} ∪ {z j : j = 1, ..., L} are introduced with the following constraints imposed: The conserved operators for this integrable model are obtained via the Quantum Inverse Scattering Method in a standard manner.Here, the key steps are noted.A transfer matrix t(u) ∈ End(H) is constructed as where T (x) is the monodromy matrix and tr a is the partial trace over an auxiliary space labelled by a.
The monodromy matrix is required to satisfy the relation which is an operator equation on V ⊗ V ⊗ H, with the two auxiliary spaces labelled by a and b.Above, is the six-vertex solution (it is convenient for our purposes to express the deformation parameter as q 2 rather than the more familiar q) of the Yang-Baxter equation [3,7,8] which acts on the three-fold space V ⊗ V ⊗ V .The subscripts above refer to the spaces on which the operators act, e.g., Two important properties of R(x), which will be called upon later, are where t 2 denotes partial transposition in the second space of the tensor product.The monodromy matrix is where and L(x) = R(q −1 x).Bearing in mind the earlier comments regarding the blocking effect, we may write where Note that U is defined in a sector-dependent manner in terms of the eigenvalues M of N , which is legitimate since N is conserved.
A consequence of Equation (10), and the diagonal form of U , is that the transfer matrices form a commutative family; i.e., The transfer matrices can be expanded in a Laurent series such that, because of Equation ( 12), the co-efficients commute Finally it can be verified that the Hamiltonian Equation ( 1), subject to the constraints of Equations (2,3), is expressible as (the corresponding expression in [17] contains typographical errors, which are corrected here) establishing that {t (j) : j = 1, ..., L} provides a set of Abelian conserved operators for the system.In this sense the system is said to be integrable.
In the remainder of this work it will be shown that for certain further restrictions on the coupling parameters there are additional Hopf algebraic symmetries of the system.These non-Abelian symmetries are not related to a quantum algebra U q (sl(2)) structure, but are realised through the Drinfel'd doubles of dihedral group algebras.

Drinfel'd Doubles of Dihedral Group Algebras
The dihedral group D n has two generators στ satisfying: where e denotes the group identity.Considering D n as a group algebra, the Drinfel'd double [10] where g are the group elements and g * are their dual elements.This gives an algebra of dimension 4n 2 .Multiplication of dual elements is defined by where δ is the Kronecker delta function.The products h * g are computed using The algebra D(D n ) becomes a Hopf algebra by imposing the following coproduct, antipode and counit respectively: An important property of D(D n ) which will be called upon later is Defining e = g∈Dn g * the universal R-matrix is given by This can be shown to satisfy the relations for a quasi-triangular Hopf algebra as defined in [10]: where ∆ T is the opposite coproduct When n is even, D(D n ) admits eight one-dimensional irreducible representations, (n 2 − 4)/2 two-dimensional irreducible representations, and eight n/2-dimensional irreducible representations.When n is odd, D(D n ) admits two one-dimensional irreducible representations, (n 2 − 1)/2 two-dimensional irreducible representations, and two n-dimensional irreducible representations.The explicit irreducible representations are given in [14].Our interest will be in the two-dimensional irreducible representations.To describe them, let ω = e 2πi/n .Then these representations have the form π 1(j) (σe) = ω j 0 0 ω −j , π 1(j) (τ e) = 0 1 1 0 , π 1(j) (eg * ) = δ(g, e) 1 0 0 1 for j = 1, ..., (n − 2)/2 if n is even and j = 1, ..., (n − 1)/2 if n is odd, for j = 1, ..., (n − 2)/2 if n is even, and for j = 1, ..., n and where k = 1, ..., (n − 2)/2 if n is even, and k = 1, ..., (n − 1)/2 if n is odd.For any of the above two-dimensional representations π µ , µ = 1(j), 2(j), (j, k) the tensor product representation applied to the universal R-matrix Equation (18) yields the general form for some l = 1, ..., n.Choosing q = ω l in Equation ( 6) we then find where is the permutation operator on the tensor product space.This shows that the Baxterisation of the D(D n ) R-matrix in two-dimensional representations leads to the symmetric six-vertex model at q a root of unity, which was previously reported in [14].Baxterisation of the D(D n ) R-matrix in higher-dimensional representations lead to the Fateev-Zamolodchikov solution of the Yang-Baxter equation, as discussed in [15,16].
Having identified the relationship Equation ( 21) between the solution Equation ( 6) of the Yang-Baxter equation and representations of the universal R-matrix Equation (18) for D(D n ), we can now proceed to determine when D(D n ) is a symmetry algebra of the transfer matrix associated to the Hamiltonian Equation (1) subject to the constraints Equations 2 and 3.

Symmetries of the Transfer Matrix and Hamiltonian
First we define We then define a modified monodromy matrix Through use of Equations (7,22,23) it can be shown that this monodromy matrix satisfies a generalised version of Equation ( 5): The transfer matrix is again defined by Equation ( 4).From the results of [18] it is known that Equation ( 12) still holds by use of Equation ( 9).The action of D(D n ) on an L-fold tensor product space is given through iterated use of the co-product action Equation ( 16): Below, for ease of notation, we will omit the representation symbols π µ when dealing with tensor product representations obtained through this action.Whenever we have the monodromy matrix Equation (24) commutes with the action of D(D n ) as a consequence of Equations (8)(9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20)(21).From the results of [19], the transfer matrix obtained from Equation (24) also commutes with the action of D(D n ) due to Equation (17).
Observing that we may write we may simplify Equation (24) as Comparing Equations 11 and 26 and taking note of Equation ( 25), these matrices are made equal by choosing meaning that the transfer matrices obtained from the monodromy matrices Equations 10 and 24 are equal.Thus we have established that the transfer matrix associated to the integrable Hamiltonian Equation (1) subject to the constraints of Equations 2 and 3 commutes with action of the quasi-triangular Hopf algebra D(D n ) whenever Equations 25 and 27 hold.
A crucial point to bear in mind is that the transfer matrices were defined in a sector-dependent manner, where each sector is associated with a fixed number of Cooper pairs.However the D(D n ) action does not preserve sectors, and specifically τ e acts as a particle-hole transformation: whereas These relations follow from the above two-dimensional matrix representations for which it is seen that representations of σe and eg * are always diagonal in the basis in which the action of N is diagonal.In the same basis, representations of τ e are orthogonal matrices with non-zero off-diagonal entries.
Recall that the Hamiltonian is defined through the transfer matrix by Equation (13).Consequently, while D(D n ) is a symmetry of the transfer matrix obtained from Equation (24) in the conventional sense, the interpretation of D(D n ) as a symmetry of the Hamiltonian is more subtle as the choice Equation ( 27) is sector-dependent and thus α needs to be treated as an operator-valued quantity.From Equation (28) we have for α given by Equation ( 27) that for each sector where N has eigenvalue M Using Equation ( 13) we obtain From the trigonometric identity this then leads to the following anti-symmetry relation for the integrable Hamiltonian Equations (1-3) whenever Equations 25 and 27 hold This relation shows how the spectrum of the Hamiltonian maps under a particle-hole transformation M → L − M induced by Equation (28).On the other hand, Thus as a result of Equations 29 and 30, the action of σe and eg * leaves the spectrum of the Hamiltonian invariant in each sector with fixed M .Finally, if the above procedure is followed using the asymmetric R-matrix a transfer matrix is obtained which commutes with the co-product action of U q (sl(2)) [20,21].However in this setting the corresponding conserved operator t (L−2) contains additional interaction terms.As a result, an expression analogous to Equation (13) does not yield an operator in the form of Equation (1).

Conclusions
An analysis of an integrable Hamiltonian for anyonic pairing, as given by Equation (1) subject to Equations 2 and 3, was undertaken.Values of the coupling parameters were identified for which the model admits Hopf algebraic symmetries.In Section 2 the construction of the integrable model was outlined in terms of the Quantum Inverse Scattering Method.This was achieved through the symmetric, six-vertex solution of the Yang-Baxter equation.The Hamiltonian was identified through a conserved operator associated to the corresponding transfer matrix.In Section 3 a description of the quasi-triangular Hopf algebra D(D n ) was presented, including explicit expressions for all irreducible, two-dimensional representations.Through these representations it was established that the symmetric, six-vertex solution of the Yang-Baxter equation is related to representations of the universal R-matrix for D(D n ).These results were utilised in Section 4 to construct a transfer matrix which preserved the D(D n ) symmetry.From this transfer matrix, values of the coupling parameters were identified for which the Hamiltonian Equation (1) subject to Equations 2 and 3 has D(D n ) as a symmetry algebra.However the interpretation of D(D n ) as a symmetry algebra for the Hamiltonian is somewhat unconventional in that both commuting and anti-commuting actions for the generators were found.The anti-commuting action is associated with a particular D(D n ) generator that induces a particle-hole transformation.