Quick search Find article
Quick search
Find article
J. Phys. A: Math. Theor. 41 No 13 (4 April 2008) 1-9
doi:10.1088/1751-8113/41/13/132001
PII: S1751-8113(08)69572-4

FAST TRACK COMMUNICATION

On the ODE/IM correspondence for minimal models

Patrick Dorey1,2, Clare Dunning3, Ferdinando Gliozzi4 and Roberto Tateo4

1 Department of Mathematical Sciences, Durham University, Durham DH1 3LE, UK
2 Isaac Newton Institute, 20 Clarkson Road, Cambridge CB3 0EH, UK
3 IMSAS, University of Kent, Canterbury CT2 7NF, UK
4 Dip. di Fisica Teorica and INFN, Università di Torino, Via P Giuria 1, 10125 Torino, Italy

E-mail: p.e.dorey@durham.ac.ukt.c.dunning@kent.ac.ukgliozzi@to.infn.it and tateo@to.infn.it

Received 4 January 2008
Published 14 March 2008

Abstract. Within the framework of the ODE/IM correspondence, we show that the minimal conformal field theories with c < 1 emerge naturally from the monodromy properties of certain families of ordinary differential equations.

PACS numbers: 02.30.Hq, 02.30.Ik, 11.25.Hf

1. Introduction

This note is about a unifying programme, the ODE/IM correspondence, which links two-dimensional conformal field theories and integrable models to the spectral theory of ordinary differential equations. The first instance of this correspondence [1] was based on an identity between the transfer matrix eigenvalues of certain integrable models in their conformal limits [2, 3], and the spectral determinants [4, 5] of second-order ordinary differential equations. Since the initial results of [1] and then [6], the ODE/IM correspondence has been used in various branches of physics ranging from condensed matter [7] to PT-symmetric quantum mechanics [8], and from boundary conformal field theory [9] to the study of non-compact sigma models [10]. It has also been linked with the geometric Langlands correspondence [11]. A recent review containing many more references is [12].

Early examples of the correspondence concerned the ground states of the integrable lattice models, albeit with possibly twisted boundary conditions. In the conformal field theory setting this gave access to primary fields, but not to their descendants. However, in [13] Bazhanov, Lukyanov and Zamolodchikov conjectured that the descendant fields could be found through relatively-simple generalizations of the initial differential equation. The new equations were obtained by modifying the initial potential, which in general has a regular singularity at zero and an irregular singularity at infinity, by introducing a level-dependent number of additional regular singularities in the complex plane, subject to a zero-monodromy condition around these extra singularities, though not about the origin. In this short note we will show that the minimal c < 1 conformal field theories can equivalently be associated with ODEs governed by a trivial monodromy property in the whole complex plane, including the origin. A natural quantization condition on the coefficient of the regular (Fuchsian) singularity at the origin emerges, and since this coefficient is related to the Virasoro vacuum parameter p of [2, 3] this restricts the resulting conformal weights, so as to match precisely the Kac tables of the minimal models {\cal M}_{\p \q} .

2. The c < 1 minimal models

We begin with the basic observation of [1, 6], that the Schrödinger equation

Equation (1)

is related to conformal field theory. The Stokes relations associated with (1) imply constraints on its eigenvalues E  in {Ei}, given suitable boundary conditions, which coincide with the Bethe Ansatz equations (BAEs) for the twisted six-vertex model in its conformal (c  =  1) limit (see, for example, [12]). The same BAEs emerge from the study of c ≤ 1 CFTs in the framework developed by Bazhanov, Lukyanov and Zamolodchikov in [2, 3]. In the notation used in [13], equation (1) encodes the primary field of a Virasoro module with central charge c, vacuum parameter p and highest weight Δ, where

Equation (2)

In [14] it was observed that for 2M rational and suitable values of l, the solutions to (1) will all lie on a finite cover of the complex plane, and that this translates into a truncation of the fusion hierarchy [15] of the associated integrable model. This is of particular interest because, for such truncations, the central charges and field contents map to those of the minimal models with c < 1 (see, for example, the discussion in section 3 of [16]). However only the simplest case of 2M integer and l(l  +  1)  =  0 was discussed explicitly in [14], in part because the general monodromies of solutions are hard to unravel in the presentation (1). One of our aims in this note is to show that a much simpler treatment is possible, from which the Kac tables of minimal model primary fields emerge in a very natural fashion.

Working backwards, we first note that, for any two coprime integers \p<\q , the ground state of the minimal model {\cal M}_{\p \q} is found by setting

Equation (3)

in (1). This corresponds to the central charge c_{\p\q}=1-\frac{6}{\p\q}(\q{-}\p)^2 and (lowest-possible) conformal weight \Delta=\frac{1}{4\p\q}(1-(\q{-}\p)^2) .

We now observe, as in section 6 of [17], that the l(l  +  1)/x2 term in (1) can be eliminated, for this value of l, by the following transformation:

Equation (4)

With a further rescaling z\to (2/\p)^{2/\q}z , (1) becomes

Equation (5)

where

Equation (6)

Note that the change of variable has replaced the singular generalized potential P(x)=x^{2\q/\p-2}-E+(1/\p^2-1/4)x^{-2} , defined on a multi-sheeted Riemann surface, by a simple polynomial W(z)=z^{\p-2} (z^{\q-\p}- \skew3\tilde{E}) . In particular, any solution to (5) is automatically single-valued around z  =  0, and the truncation of the fusion hierarchy as explained in section 4 of [14] is made much more transparent.

To see which other primary states in the original model might have similarly-trivial monodromy, we keep l real with l  +  1/2 > 0, but otherwise arbitrary, and again perform the change of variable (4). The result is now

Equation (7)

where

Equation (8)

The Fuchsian singularity in (7) at z  =  0 means that the equation admits a pair of solutions which generally have the power series expansions

Equation (9)

where \lambda_1=\tilde{l} +1 and \lambda_2=-\tilde{l} are the two roots of the indicial equation

Equation (10)

A general solution to (7) can be expressed as y(z,\skew3\tilde{E},\tilde l)=\sigma\chi_1(z)+\tau \chi_2(z) , and we shall demand that the transformed ODE should be such that, for arbitrary \skew3\tilde{E} , the monodromy of y(z,\skew3\tilde{E},\tilde l) around z  =  0 is projectively trivial, which is to say that

Equation (11)

This condition ensures that the eigenvalues obtained by imposing the simultaneous decay of solutions in a pair of asymptotic directions at infinity are independent of the path of analytic continuation between these two directions. We shall show that the condition imposes the following constraints on \tilde l :

  • 2\tilde l+1 is a positive integer;
  • The allowed values of 2\tilde l+1 are those integers which cannot be written as \p s+\q t with s and t non-negative integers. In other words they form precisely the set of holes of the infinite sequence

    Equation (12)

    We shall call the integers (12) `representable' and denote the set of them by \Rth_{\p\q} .

  • As a consequence we shall see that, as \tilde l runs over its `allowed' values, the rational numbers

    Equation (13)

    precisely reproduce the set of conformal weights of the primary states lying in the Kac table of the minimal model {\cal M}_{\p \q} . Figures 1 and 2 illustrate the story for the Ising and Yang–Lee cases.

    Figure 1

    Figure 1. The holes (open circles) in the infinite sequence of integers defined in (12) for the critical Ising model {\cal M}_{3 4} . The holes are at 1, 2 and 5; the resulting conformal weights according to (13) are also shown, matching the primary field content of the Ising model.

    Figure 2

    Figure 2. Holes for the Lee–Yang model {\cal M}_{2 5} , at 1 and 3. Notation as in figure 1.

    To establish these claims we first note that the requirement that the general solution y(z) be projectively trivial means that χ1(z) and χ2(z) must have the same monodromy, which implies that the two roots of the indicial equation must differ by an integer

    Equation (14)

    This in turn restricts \tilde l to be an integer or half integer, so that, naïvely, the allowed solutions are even or odd under a 2π rotation around z  =  0. However it is well known that in such a circumstance, while χ1(z) keeps its power series expansion (9), χ2(z) generally acquires a logarithmic contribution:

    Equation (15)

    Unless D  =  0, this will spoil the projectively trivial monodromy of y(z). We now show that D  =  0 if and only if 2\tilde l+1 obeys the constraint (ii). In fact the logarithmic term is only absent when the recursion relations for dn's with D  =  0 admit a solution. These relations are

    Equation (16)

    with the initial conditions d0  =  1, dm<0  =  0.

    Consider first the situation when 2 \tilde l+1 \notin \Rth_{\p\q} . Then, starting from the given initial conditions, the recursion relation (16) generates a solution of the form

    Equation (17)

    where the only nonzero dn's are those for which the label n lies in the set \Rth_{\p\q} . Given that 2\tilde l+1\notin \Rth_{\p\q} , for these values of n the factor n(n-2\tilde l-1) on the lhs of (16) is never zero, and hence this procedure is well defined.

    If instead 2\tilde l+1\in\Rth_{\p\q} , then equation (16) taken at n=2\tilde l+1 yields the additional condition

    Equation (18)

    which is inconsistent for generic \skew3\tilde{E} , and so the logarithmic term is required1Note1.

    Given the characterization (12) of \Rth_{\p\q} , the set {\bb Z}^+ of non-negative integers can be written as a disjoint union

    Equation (19)

    where \NN_{\p\q} is the set of `nonrepresentable' integers. If the coprime integers \p and \q are larger than 1 then \NN_{\p\q} is non-empty; in fact |\NN_{\p\q}|=\frac{1}{2}(\p{-}1)(\q{-}1) , a result which goes back to Sylvester [18].

    To characterize \NN_{\p\q} more precisely, we start with the fact that given two coprime integers \p and \q , any integer n can be written as

    Equation (20)

    This is a classical result of number theory. An intuitive proof uses Euclid's algorithm for the greatest common divisor of two integers. Alternatively one can invoke the Euler totient function varphi2Note2 and the following theorem (see, for example, [19])

    Equation (21)

    which implies that \p^{\varphi(\q)} =1+h\q for some integer h. This immediately yields the solution s_0= n\p^{\varphi(\q)-1}, t_0=-n h . However, for any given n, this is not the only possibility. More precisely, we have \p s_0+\q t_0=\p s+\q t for some other pair of integers (s, t) if and only if s=s_0+\q k, t=t_0-\p k for some k\in{\bb Z} . This is easily proved: rearrange (20) as

    Equation (22)

    Since (\p,\q)=1, \q must be a factor of s0s and so s_0-s=-\q k for some k. Dividing (22) by \q then shows that t_0-t=\p k , as required. Hence the possible representatives for each integer n constitute the line of points (s,t)=(s_0+\q k,t_0-\p k), k\in{\bb Z} . For n to be a positive integer, \p s+\q t\gt0 , or s\gt-\q t/\p . If none of these points has both coordinates non-negative, then the corresponding n will be in \NN_{\p \q} . To keep t non-negative while making s as large as possible, we shift t by a multiple of \p so that 0\le t<\p . If s is still negative, then n will be in \NN_{\p \q} . The numbers we want are therefore represented by the points

    Equation (23)

    Negating t, the allowed (trivial monodromy) values of 2\tilde l+1 are therefore

    Equation (24)

    Figure 3 illustrates the argument.

    Figure 3

    Figure 3. A graphical representation of the nonrepresentable integers for \p=3, \q=5, \NN_{\p \q}=\{1,2,4,7\} . The elements of \NN_{\p \q} correspond to the four unshaded points. Three of the lines (s,t)=(s_0+\q k,t_0-\p k) have also been shown; each such line contains exactly one point in the region 0\le t<\p between the two dotted horizontal lines.

    Substituting back using (8) and (2), the allowed values for the conformal weights Δ precisely reproduce the Kac table for the minimal model {\cal M}_{\p \q} :

    Equation (25)

    It is striking that the full Kac table should emerge from such a simple consideration of the monodromy properties of the transformed differential equation (7). Finally, note that everything is symmetric in \p and \q so the same result can be obtained by starting from M+1=\p/\q instead.

    Another way to characterize \NN_{\p \q} is through the generating function

    Equation (26)

    It is straightforward to show that such a rational function is actually a polynomial, because all the zeros of the denominator are cancelled by zeros of the numerator; we would like to show that

    Equation (27)

    To see this, first note that Taylor expanding 1/(1-z^\p)(1-z^\q) yields

    Equation (28)

    In order to get rid of the unknown cn's, we combine the trivial identity

    Equation (29)

    with the similar one for 1/(1-z^\q) to write

    Equation (30)

    As a consequence the generating function of representable integers is

    Equation (31)

    Therefore the difference

    Equation (32)

    is the sought after formula (27).

    Before concluding this section we would like to mention that there is another ODE that can be associated with the same series of minimal models. This is the so-called A(2)2 description, related to phi12, phi21 and phi15 perturbations [20]. After a simple change of variable, the relevant third-order phi12-related ODE can be cast into the form

    Equation (33)

    Swapping \p and \q gives the phi21-related ODE, while replacing \p with 2\p and \q with \q/2 yields the phi15 equation. In all cases, the indicial equation is

    Equation (34)

    The zero relative monodromy condition among the three solutions to (33) requires (\tilde g+1) to be a (positive) integer, and to avoid logarithmic terms we should also simultaneously impose the following two conditions:

    Equation (35)

    with s, t  =  0, 1, 2,  .... For \p odd, it is easy to check that equations (35) lead to the same set of integers as the su(2)-related case discussed above, while only a subset is recovered for \p even. For the phi21-related case, the opposite situation occurs: for \p even the full Kac table is recovered, while for \p odd only a subset is found.

    3. Further generalizations and conclusions

    There are many possible generalizations of the above results. The existence of a simpler version, equation (7), of the basic ODE for minimal CFTs is not restricted to c < 1 Virasoro models, but generalizes to the higher su(2) coset CFTs of [21] and to the ABCD-related theories of [22]. The pseudo-differential equations listed in section 3 of [22] include the minimal models

    Equation (36)

    at fractional level L=K\p/(\q-\p)-h^\vee with \q-\p=K u , and u  =  1, 2,  ... (using the notation of appendix 18.B of [23]). We checked that after simple changes of variable, these equations reduce to equations similar in form to the originals, but for a change in the `potential', as follows:

    Equation (37)

    It is striking that when both L and K are integers the CFT is unitary and W(K,L)(z) simplifies further to

    Equation (38)

    Equation (38) motivates some simple comments and speculations. First we observe that the K leftrightarrow L invariance of (36) manifests itself in (38) as a shift in z. As a consequence of this symmetry, lateral quantization problems for the ground-state ODEs and the associated Stokes multipliers are, up to \skew3\tilde{E} \rightarrow -\skew3\tilde{E} , invariant under the exchange of L and K. However, this symmetry is explicitly broken in equations with extra Fuchsian singularities as in (7). A possible remedy is to treat the points z  =  0 and z^{(\q-\p)/K}= \tilde{E} more democratically. For instance, in unitary models the symmetry is globally restored after the addition of a second singularity at z=\skew3\tilde{E} , ensuring that the set of ODEs for the primary fields in a given CFT is mapped into itself by the transformation. We suspect that a similar modification may also resolve the problem of the missing states in the A(2)2 example discussed at the end of section 2.

    A further possibility is suggested by equation (38). Consider the following multi-parameter generalization of (38):

    Equation (39)

    where the constants ei (i  =  1, 2 ..., K  +  L – 1) are free parameters. To keep the discussion brief, we shall restrict attention to the ground-state equation for L  +  K  =  3 and Gg = su(2). Then if e1  =  0 and e_2=\skew3\tilde{E} the corresponding ODE is related to the tricritical Ising model {\cal M}_{4 5} , while for (e_1,e_2)=(\pm \sqrt{\skew3\tilde{E}},\mp \sqrt{\skew3\tilde{E}}) the equation corresponds to {\cal M}_{3 5} . Therefore this simple two-parameter model interpolates smoothly between equations associated with {\cal M}_{4 5} and {\cal M}_{3 5} . This phenomenon has a counterpart in the homogeneous sine-Gordon model corresponding to integrable perturbations of the \hat{su}(3)_2/U(1)^2 coset model. The thermodynamic Bethe ansatz (TBA) equations for this model [24, 25] have two independent scale parameters μ1 and μ2. If one of these parameters is set to zero, the TBA equations reduce to those for {\cal M}_{4 5}+\phi_{13} [26], while for μ1  =  μ2 the TBA equations map into a pair of identical equations for {\cal M}_{3 5}+\phi_{13} [27]. This, and other simple considerations, suggest that ODEs with multi-parameter potentials of the form (39) may have an interesting interpretation in terms of conformal field theory. Much more work will be needed in order to give this observation a more solid grounding, but we feel that it will be an interesting direction for future exploration.

    Acknowledgments

    We are very grateful to Junji Suzuki and Edward Frenkel for useful conversations. We would also like to take this opportunity to acknowledge the great debt over many years that our whole subject owes to Aliocha Zamolodchikov, of whose untimely passing we heard while working on this paper. PED and TCD thank Torino University and the INFN for hospitality at the beginning of this project, and PED thanks the Isaac Newton Institute for hospitality at the end. This project was also partially supported by grants from the Nuffield Foundation, grant number NAL/32601, and from the Leverhulme Trust.

    References
    [1] 
    Dorey P and Tateo R 1999 Anharmonic oscillators, the thermodynamic Bethe Ansatz, and nonlinear integral equations J. Phys. A: Math. Gen. 32 L419 (Preprint hep-th/9812211) 
    IOPsciencePreprint
    [2] 
    Bazhanov V V, Lukyanov S L and Zamolodchikov A B 1996 Integrable structure of conformal field theory, quantum KdV theory and thermodynamic Bethe ansatz Commun. Math. Phys. 177 381 (Preprint hep-th/9412229) 
    CrossRefPreprint
    [3] 
    Bazhanov V V, Lukyanov S L and Zamolodchikov A B 1997 Integrable structure of conformal field theory: II. Q-operator and DDV equation Commun. Math. Phys. 190 247 (Preprint hep-th/9604044) 
    CrossRefPreprint
    [4] 
    Sibuya Y 1975 Global Theory of a Second-order Linear Ordinary Differential Equation with Polynomial Coefficient (Amsterdam: North-Holland)  
    [5] 
    Voros A 1982 Semi-classical correspondence and exact results: the case of the spectra of homogeneous Schrödinger operators J. Phys. Lett. 43 L1 
    CrossRef
    Voros A 1983 The return of the quartic oscillator. The complex WKB method Ann. Inst. Henri Poincaré XXXIX 211 
    [6] 
    Bazhanov V V, Lukyanov S L and Zamolodchikov A B 2001 Spectral determinants for Schroedinger equation and Q-operators of conformal field theory J. Stat. Phys. 102 567 (Preprint hep-th/9812247) 
    CrossRefPreprint
    [7] 
    Gritsev V, Altman E, Demler E and Polkovnikov A 2006 Full quantum distribution of contrast in interference experiments between interacting one dimensional Bose liquids Nature Phys. 2 705 (Preprint cond-mat/0602475) 
    CrossRefPreprint
    [8] 
    Dorey P, Dunning C and Tateo R 2001 Spectral equivalences, Bethe Ansatz equations, and reality properties in \cal{P}\cal{T} -symmetric quantum mechanics J. Phys. A: Math. Gen. 34 5679 (Preprint hep-th/0103051) 
    IOPsciencePreprint
    [9] 
    Lukyanov S L, Vitchev E S and Zamolodchikov A B 2004 Integrable model of boundary interaction: the paperclip Nucl. Phys. B 683 423 (Preprint hep-th/0312168) 
    CrossRefPreprint
    [10] 
    Teschner J 2007 On the spectrum of the Sinh-Gordon model in finite volume Preprint hep-th/0702214  
    Preprint
    [11] 
    Feigin B and Frenkel E 2007 Quantization of soliton systems and Langlands duality Preprint 0705.2486  
    Preprint
    [12] 
    Dorey P, Dunning C and Tateo R 2007 The ODE/IM Correspondence J. Phys. A: Math. Theor. 40 R205 (Preprint hep-th/0703066) 
    IOPsciencePreprint
    [13] 
    Bazhanov V V, Lukyanov S L and Zamolodchikov A B 2004 Higher-level eigenvalues of Q-operators and Schroedinger equation Adv. Theor. Math. Phys. 7 711 (Preprint hep-th/0307108) 
    Preprint
    [14] 
    Dorey P and Tateo R 1999 On the relation between Stokes multipliers and the T-Q systems of conformal field theory Nucl. Phys. B 563 573 
    CrossRef
    Dorey P and Tateo R 2001 Nucl. Phys. B 603 581 (Preprint hep-th/9906219) (erratum) 
    CrossRefPreprint
    [15] 
    Klümper A and Pearce P A 1992 Conformal weights of RSOS lattice models and their fusion hierarchies Physica A 183 304 
    CrossRef
    [16] 
    Bazhanov V V, Lukyanov S L and Zamolodchikov A B 1997 Quantum field theories in finite volume: excited state energies Nucl. Phys. B 489 487 (Preprint hep-th/9607099) 
    CrossRefPreprint
    [17] 
    Dorey P, Millican-Slater A and Tateo R 2005 Beyond the WKB approximation in PT-symmetric quantum mechanics J. Phys. A: Math. Gen. 38 1305 (Preprint hep-th/0410013) 
    IOPsciencePreprint
    [18] 
    Sylvester J J 1884 Question 7382 Mathematical Questions from the Educational Times 41 21 
    [19] 
    Lejeune Dirichlet P G 1881  `Lezioni sulla teoria dei numeri' published and complemented by appendices by R Dedekind
    [20] 
    Dorey P and Tateo R 2000 Differential equations and integrable models: the SU(3) case Nucl. Phys. B 571 583 
    CrossRef
    Dorey P and Tateo R 2001 Nucl. Phys. B 603 582 (Preprint hep-th/9910102) (erratum) 
    CrossRefPreprint
    [21] 
    Lukyanov S L 2007 Notes on parafermionic QFT's with boundary interaction Nucl. Phys. B 784 151 (Preprint hep-th/0606155) 
    CrossRefPreprint
    [22] 
    Dorey P, Dunning C, Masoero D, Suzuki J and Tateo R 2007 Pseudo-differential equations, and the Bethe Ansatz for the classical Lie algebras Nucl. Phys. B 772 249 (Preprint hep-th/0612298) 
    CrossRefPreprint
    [23] 
    Di Francesco P, Mathieu P and Senechal D 1997 Conformal Field Theory (New York: Springer)  
    CrossRef
    [24] 
    Castro-Alvaredo O A, Fring A, Korff C and Miramontes J L 2000 Thermodynamic Bethe ansatz of the homogeneous sine-Gordon models Nucl. Phys. B 575 535 (Preprint hep-th/9912196) 
    CrossRefPreprint
    [25] 
    Dorey P and Miramontes J L 2004 Mass scales and crossover phenomena in the homogeneous sine-Gordon models Nucl. Phys. B 697 405 (Preprint hep-th/0405275) 
    CrossRefPreprint
    [26] 
    Zamolodchikov Al B 1991 Thermodynamic Bethe ansatz for RSOS scattering theories Nucl. Phys. B 358 497 
    CrossRef
    [27] 
    Ravanini F, Tateo R and Valleriani A 1993 Dynkin TBAs Int. J. Mod. Phys. A 8 1707 (Preprint hep-th/9207040) 
    CrossRefPreprint

    Notes

    Note1
     The reader might wonder whether the lhs of (18) could vanish identically as a result of an exceptional cancellation between the two terms. This can be ruled out by observing that the monodromy property in z is unaltered by the change of variable z \rightarrow \beta^{1/\q}z . This rescaling leads to a more general version of (16): n (n-2\tilde l-1) d_n=\alpha d_{n-\p}+\beta d_{n-\q} with \alpha=-\skew3\tilde{E} \beta^{\p/\q} . Then for the particular choice \alpha=(-1)^{\p+1}, \beta=(-1)^{\q+1} , one can see that \alpha d_{n-\p} and \beta d_{n-\q} must have the same sign for n \le 2\tilde l +1 ; hence, they cannot cancel identically.
    Note2  \varphi(\q) denotes the number of coprimes with \q in the set 1,2,\ldots,\q .



    Please login to access our web services, or create an account if you don't yet have one.

    You must have cookies enabled in your web browser to be able to login.

    Username
    Password

    Forgotten your password? Get a new one here.