text:
# 1 Introduction

## 1 Introduction

The sine-Gordon (SG) and supersymmetric sine-Gordon (SSG)  models belong to an interesting set of integrable quantum field theories that are said to have fractional supersymmetry , thus attracting the name fractional supersymmetric sine-Gordon models (FSSG). The field content  consists of a single boson interacting with a $`Z_L`$ generating parafermion, the choices $`L=1`$ and $`2`$ yielding the SG and SSG models respectively. The models correspond to the fractional conformal field theory obtained from the Wess-Zumino-Witten model $`SU(2)_L`$ by compactifying the free boson on a circle of radius $`2R=\beta /\sqrt{4\pi }`$, perturbed by

$$\mathrm{\Phi }(z,\overline{z})=\frac{g\sqrt{4\pi }}{\beta }\mathrm{\Psi }_1(z)\overline{\mathrm{\Psi }}_1(z)e^{i\frac{\beta }{\sqrt{4}\pi }\phi (z,\overline{z})}+\mathrm{c}.\mathrm{c}..$$
(1)
The mass scale is set by the dimensionful parameter $`g`$, and the real dimensionless parameter $`\beta `$ determines the particle content. Since the boson $`\phi (z,\overline{z})`$ and generating parafermion $`\mathrm{\Psi }_1(z,\overline{z})`$ have conformal dimensions $`1`$ and $`11/L`$ respectively, the perturbing term has conformal dimension $`\beta ^2/8\pi +11/L`$. The ultraviolet effective central charge is $`1+2(L1)/(L+2)=3L/L+2`$.

In the limit $`\beta ^2/8\pi 1/L`$ the perturbed model becomes the current-current perturbation of the WZW model $`SU(2)_L`$ . Instead if we tune the parameter

$$\xi =\frac{L\beta ^2/8\pi }{1/L\beta ^2/8\pi }$$
(2)
to $`K+2`$ we obtain the unitary coset $`SU(2)_L\times SU(2)_K/SU(2)_{L+K}`$ perturbed by the operator $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ . The indices indicate which representation of $`SU(2)_L,SU(2)_K`$ and $`SU(2)_{L+K}`$ is associated to the corresponding weights. This choice of $`\xi `$ when $`L=1`$ corresponds to the well-known quantum group restriction  of the sine-Gordon model to the $`\varphi _{13}`$ perturbed minimal models. In fact, the FSSG models were first proposed as the models obtained by ‘unrestricting’ the S-matrices describing the $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ perturbation of the cosets $`_{L,K}SU(2)_L\times SU(2)_K/SU(2)_{L+K}`$.

The existence of higher spin conserved charges in 1+1 dimensional integrable quantum field theories strongly constrain the allowed scattering processes. There is no particle production, the individual particle momenta are conserved, and all $`n`$-particle scattering processes decompose into products of two-particle scattering amplitudes . The S-matrix of a supersymmetric model should also commute with the supersymmetry generators. This typically leads to a further factorisation of the two-particle S-matrix into a direct product in which one factor encodes the supersymmetry dependence . The S-matrices of the FSSG models exactly fit this pattern :

$$S(\beta )=S_LS_{\mathrm{SG}}(\beta ).$$
(3)
The supersymmetry-related factor $`S_L`$ describes the scattering of the massive $`\varphi _{13}`$ perturbation of the minimal model $`_{1,L}`$ , and the bosonic factor $`S_{\mathrm{SG}}(\beta )`$ is the S-matrix of the sine-Gordon model at coupling $`\beta `$  . The particle content of the FSSG consists of a soliton of mass $`M`$ and its antisoliton together with a number of breathers of mass

$$M_j=2M\mathrm{sin}\left(\frac{\pi j\xi }{2}\right),j=1\mathrm{}<\frac{1}{\xi }.$$
(4)
There are no breathers in the repulsive regime $`\xi <1`$. Once the spectrum and associated S-matrices have been conjectured the thermodynamic Bethe ansatz  is usually employed to check the proposed scattering theory. However when the scattering is nondiagonal, as above, deriving the TBA equations from the S-matrices becomes rather difficult.

Fendley and Intriligator noted that one can immediately write down a TBA system for a model with a tensor product S-matrix of the form $`S_GS_H`$ if the TBA equations are already known for the models described by S-matrices $`S_G`$ and $`S_H`$ . If the individual TBA equations are encoded on Dynkin-like diagrams of type $`G`$ and $`H`$, each having a single massive node, then the TBA equations for the model with S-matrix $`S_GS_H`$ are found by gluing the individual TBA equations at the massive node . In the compact notation of  this corresponds to $`(GH)_L`$, where $`L`$ indicates the massive node. It is easy to check that the TBA equations found for the FSSG model at the special values of the coupling

$$\frac{\beta ^2}{8\pi }=\frac{1}{L}\frac{2}{m},m^+,L=1,2,\mathrm{}\mathrm{int}[m/2]1,$$
(5)
and those for the $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ perturbation of the $`_{L,K}`$ models  all fit this pattern. TBA equations for the models arising as the $`\mathrm{\Phi }_{13}`$ perturbation of the supersymmetric models $`𝒮_{2,4n+4}`$ are known , but currently lack an interpretation in terms of Dynkin-like diagrams.

A $`Y`$-system, from which the TBA equations can be obtained via Fourier transform , for the FSSG model at any coupling has been proposed using the gluing idea and the sine-Gordon $`Y`$-system . The SG $`Y`$-system depends on writing the smaller of $`\xi `$ or $`1/\xi `$ as a simple continued fraction. The number of TBA equations is found to be equal to the sum of the leading terms of the denominators in the continued fraction plus one further equation. Though there is an neat depiction in terms of Dynkin diagrams  it is hard to deal with a large number of equations. Moreover those $`\xi `$ for which the continued fraction does not terminate are associated to a TBA system with an infinite number of equations.

There is an alternative method of calculating the finite size effects of the sine-Gordon model, known as the nonlinear integral equation (NLIE) technique. The result is a single equation, first proposed for the groundstate of the sine-Gordon model by Destri and De Vega in 1992 . A complete derivation to include the excited states of even topological charge followed in , and an equation for the excited states of odd topological charge was proposed in . The conformal limit of the groundstate equation had been proposed in 1991 in an integrable lattice model context by Klümper, Batchelor and Pearce . Further motivation for this work came from the the lattice model connection: the $`XXX`$ spin chain at spin $`S`$ is a model with underlying field theory $`SU(2)_{2S}`$. The spin-$`\frac{1}{2}`$ case corresponds to the lattice model mentioned above; the spin-$`1`$ problem is discussed in . Recently Suzuki has generalised the technique of  to obtain a a finite set of nonlinear integral equations that describe the thermodynamics of the spin-$`S`$ chain at finite temperature . Given the relations between the spin chains and $`SU(2)_{2S}`$, and the FSSG models and perturbed $`SU(2)_L`$, we should expect the conformal limit of any integral equation describing the FSSG at $`\beta ^2=1/L`$ to be closely related to the zero temperature version of Suzuki’s NLIE.

Motivated by Suzuki’s equations and the gluing idea we propose a set of nonlinear integral equations to describe the groundstate energy of the fractional supersymmetric sine-Gordon model at any value of the coupling constant $`\beta ^2`$. In §2 a short description of the component equations leads to the proposed NLIEs, and then a number of tests are carried out in §3.

The $`N=2`$ supersymmetric sine-Gordon model is not a member of the FSSG series. However the groundstate energy of this model can be studied by means of a further set of NLIEs, obtained by a simple generalisation of the FSSG equations. The details are in §4 along with a short discussion of related models: the sausage models and the parafermionic perturbation of the $`Z_M`$ models. Directions for future work can be found in §5.

## 2 Nonlinear integral equations

There are two integrable quantum field theories associated with the $`\varphi _{13}`$ perturbation of the minimal model $`_{1,L}`$. If the coupling $`\lambda `$ is positive the action

$$𝒜=𝒜_{\mathrm{CFT}}+\lambda d^2x\varphi _{13}$$
(6)
describes a ‘massless flow’, or renormalisation group trajectory, from $`_{1,L}`$ in the ultraviolet to $`_{1,L1}`$ in the infrared . On the other hand, if the coupling is negative (6) describes a massive integrable quantum field theory with S-matrix $`S_L`$ . The TBA equations for the groundstate energy of the massive scattering theory were found by Al.B. Zamolodchikov . As we explain below one small change allows the same equations to unexpectedly describe the groundstate energy of the massless flows . The equations are written in terms of the incidence matrix $`I_{ab}`$ of the $`A_L`$ Dynkin diagram . The TBA system is***Note that the TBA equations are usually written in terms of pseudoenergy $`\epsilon =f`$.

$`f_a(\theta )`$ $`=`$ $`g_a(\theta )+{\displaystyle \underset{b=1}{\overset{L}{}}}{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑\theta ^{}I_{ab}\chi (\theta \theta ^{})\mathrm{ln}(1+e^{f_b(\theta ^{})})`$
$`c_{\mathrm{eff}}(r)`$ $`=`$ $`{\displaystyle \frac{6}{\pi }}{\displaystyle \underset{a=1}{\overset{L}{}}}{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑\theta g_a(\theta )\mathrm{ln}(1+e^{f_a(\theta )}),`$ (7)
with kernel $`\chi (\theta )=1/(2\pi \mathrm{cosh}\theta ).`$ Here $`r=MR`$ where $`M`$ is the mass of the elementary kink, or in the massless theory the crossover scale, and $`R`$ is the circumference of the cylinder on which the model is considered. The massive theory is found by setting the driving term to

$$g_a(\theta )=\delta _{a,1}r\mathrm{cosh}\theta $$
(8)
while the massless theory appears if

$$g_a(\theta )=\frac{r}{2}(\delta _{a,L}e^\theta +\delta _{a,1}e^\theta ).$$
(9)
Both theories have the same ultraviolet effective central charge: $`c_{\mathrm{eff}}(0)=16/(L+2)(L+3)`$, while the infrared value is zero or $`16/(L+1)(L+2)`$ for the massive/massless models respectively.

To avoid dealing with a different and possibly infinite set of equations for each choice of $`\beta ^2`$ we consider the second type of nonlinear integral equation mentioned in the introduction. The equation involves a complex function $`a(\theta )`$ and its conjugate $`\overline{a}(\theta )`$ coupled via

$`\mathrm{ln}a(\theta )`$ $`=`$ $`ir\mathrm{sinh}\theta +{\displaystyle _{𝒞_1}}𝑑\theta ^{}\phi (\theta \theta ^{})\mathrm{ln}(1+a(\theta ^{})){\displaystyle _{𝒞_2}}𝑑\theta ^{}\phi (\theta \theta ^{})\mathrm{ln}(1+a^1(\theta ^{})).`$ (10)
The integration contours $`𝒞_1`$ and $`𝒞_2`$ run from $`\mathrm{}`$ to $`+\mathrm{}`$, just below and just above the real $`\theta `$-axis respectively. The kernel $`\phi (\theta )`$ is proportional to the logarithmic derivative of the soliton-soliton scattering amplitude of the sine-Gordon model

$$\phi (\theta )=\frac{dk}{2\pi }e^{ik\theta }\frac{\mathrm{sinh}(\xi 1)\frac{\pi k}{2}}{2\mathrm{sinh}\frac{\pi \xi k}{2}\mathrm{cosh}\frac{\pi k}{2}}.$$
(11)
The effective central charge found using

$$c_{\mathrm{eff}}(r)=\frac{3ir}{\pi ^2}\left(_{𝒞_1}𝑑\theta \mathrm{sinh}\theta \mathrm{ln}(1+a(\theta ))_{𝒞_2}𝑑\theta \mathrm{sinh}\theta \mathrm{ln}(1+a^1(\theta ))\right)$$
(12)
has ultraviolet value $`1`$. With constant $`i\pi \alpha `$  added to the RHS of (10) the effective central charge becomes $`c_{\mathrm{eff}}(0)=16\alpha ^2\xi /(\xi +1)`$. The choice $`\xi =p/(qp)`$ and $`\alpha =\pm 1/p`$  implements the restriction  of the SG model to the $`\varphi _{13}`$ perturbed minimal models.

The explanation for Fendley and Intriligator’s gluing proposal lies in the derivation of the TBA equations . The key is the factorisability of the S-matrix: it implies each factor can be diagonalised separately and each generates an individual set of pseudoparticles. These are seen to only couple through the massive particle, with the result being the ‘factorised’ or glued TBA system described in the introduction. There is no equivalent of the pseudoparticles in the derivation of the sine-Gordon NLIE. Moreover, in further contrast to the TBA, the S-matrix plays no role in the NLIE derivation. Instead one notes a posteriori that the kernel function of the resulting NLIE is related as described above to the sine-Gordon S-matrix. It therefore seemed unlikely that the gluing idea would be of any use in finding a NLIE for a model with an S-matrix of the form (3). Fortunately this turns out not to be the case.

We obtain the NLIEs for the FSSG model by gluing the minimal model TBA equations to the sine-Gordon NLIE at the massive node following the Fendley-Intriligator prescription . However the resulting equations are not quite right and it is necessary to make a few simple modifications. Due to singularities in $`\mathrm{ln}(1+y_1(\theta ))`$ for real $`\theta `$ a TBA-like convolution $`𝑑\theta ^{}\chi (\theta \theta ^{})\mathrm{ln}(1+y_1(\theta ^{}))`$ along the real axis cannot be used. The solution is to replace this term with a convolution along the contours $`𝒞_1`$ and $`𝒞_2`$ encircling the real axis as in the SG NLIE. To ensure for real $`\theta `$ the reality of the $`L1`$ TBA-like functions $`y_2,\mathrm{},y_L`$, and the complex nature of the NLIE-like function $`y_1`$ we find we should couple $`y_1`$ and $`y_2`$ to each other with an extra shift of $`i\pi /2`$ in the kernel $`\chi `$. The conjectured equations are

$`\mathrm{ln}y_1(\theta )`$ $`=`$ $`\nu _1(\theta )+i\pi \alpha `$
$`+{\displaystyle _{𝒞_1}}𝑑\theta ^{}\phi (\theta \theta ^{})\mathrm{ln}(1+y_1(\theta ^{})){\displaystyle _{𝒞_2}}𝑑\theta ^{}\phi (\theta \theta ^{})\mathrm{ln}(1+y_1^1(\theta ^{}))`$
$`+{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑\theta ^{}\chi (\theta \theta ^{}+\frac{i\pi }{2})\mathrm{ln}(1+y_2(\theta ^{}));`$
$`\mathrm{ln}y_2(\theta )`$ $`=`$ $`\nu _2(\theta )+{\displaystyle _{𝒞_1}}𝑑\theta ^{}\chi (\theta \theta ^{}\frac{i\pi }{2})\mathrm{ln}(1+y_1(\theta ^{})){\displaystyle _{𝒞_2}}𝑑\theta ^{}\chi (\theta \theta ^{}\frac{i\pi }{2})\mathrm{ln}(1+y_1^1(\theta ^{}))`$
$`+{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑\theta ^{}\chi (\theta \theta ^{})\mathrm{ln}(1+y_3(\theta ^{}));`$
$`\mathrm{ln}y_a(\theta )`$ $`=`$ $`\nu _a(\theta )+{\displaystyle \underset{b=1}{\overset{L}{}}}I_{ab}{\displaystyle _{\mathrm{}}^{\mathrm{}}}d\theta ^{}\chi (\theta \theta ^{})\mathrm{ln}(1+y_b(\theta ^{})),a=3\mathrm{}L,`$ (13)
where $`\chi (\theta )`$ remains as above. The driving terms are

$$\nu _a(\theta )=\{\begin{array}{cc}ir\mathrm{sinh}(\theta )& a=1\hfill \\ 0& \mathrm{otherwise}\hfill \end{array}.$$
(14)
Borrowing the TBA notation of  our shorthand for equations of this type will be $`(A_1A_L)_a^{[i]}`$. The superscript $`[i]`$ indicates the NLIE-like function and the subscript $`a`$ the function on which the mass is placed. The equations above are therefore represented by $`(A_1A_L)_1^{[1]}`$. The effective central charge defined in terms of the solution of (13) can be calculated using

$$c_{\mathrm{eff}}(r)=\frac{3}{\pi ^2}\underset{a=1}{\overset{L}{}}_{𝒞_1}𝑑\theta \nu _a(\theta )\mathrm{ln}(1+y_a(\theta ))_{𝒞_2}𝑑\theta \nu _a(\theta )\mathrm{ln}(1+y_a^1(\theta )).$$
(15)
At $`L=1`$ the system (13,14,15) trivially reduces to the NLIE describing the sine-Gordon equation. As expected these equations become equivalent to Suzuki’s  for the spin-$`S`$ $`XXX`$ spin chain in the limit $`\beta ^21/L`$ ($`\xi \mathrm{}`$) . The NLIE-like functions are approximately related via $`y_1(\theta )\mathrm{ln}𝔞(\frac{2}{\pi }\theta +i)`$, provided the spin chain driving term is replaced with the one relevant for the field theory.

Based on previous experience with NLIEs we expect to find the effective central charge for the FSSG model at level $`L`$ and coupling $`\beta ^2`$ by tuning $`\xi `$ according to (2) and setting $`\alpha =0`$. The $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ perturbation of the coset $`_{L,K}`$ should be found at

$$\xi =K+2,\alpha =\pm 1/(K+2).$$
(16)
In fact we hope the equations go a little further: even though the quantum group restriction of the S-matrices may not work it is thought that the FSSG models also describe the nonunitary cosets perturbed by $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$. If we set

$$\xi =\frac{Lp}{qp},\alpha =\pm \frac{1}{p},$$
(17)
then the NLIEs should yield the effective central charge of the $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ perturbation of $`_{L,Lp/(qp)2}`$. This is a coset at rational, and possibly negative, level. These models make sense provided the integers $`p,q`$ are such that $`p`$ and $`(qp)/L`$ are coprime integers. At $`L=1`$ this is the $`\varphi _{13}`$ perturbation of the minimal models $`(p,q)`$, and at $`L=2`$ the supersymmetry preserving perturbation $`\mathrm{\Phi }_{13}`$ of the superminimal models $`𝒮(p,q)`$.

In the same way that TBA equations for the massive perturbation of the minimal models can be converted into those describing the massless perturbations, we can describe some massless flows using the proposed NLIEs with new driving terms

$$\nu _a(\theta )=\{\begin{array}{cc}\frac{i}{2}re^\theta \hfill & a=1\hfill \\ \frac{1}{2}re^\theta \delta _{a,LK+1}\hfill & \mathrm{otherwise}\hfill \end{array}.$$
(18)
If $`(\xi ,\alpha )`$ are fixed to the values prescribed in (16) the equations reproduce the flow $`_{L,K}+\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}_{L,KL}`$. Given the different nature of $`y_1`$ and $`y_{LK+1}`$ it is rather surprising that the resulting NLIEs do in fact describe massless flows. The massless flows within the sine-Gordon model ($`L=1`$) can alternatively be studied via two coupled NLIE functions .

Using the symmetry between $`L`$ and $`K`$ we can minimise the number of NLIEs for the coset models $`_{L,K}`$ by choosing $`L<K`$. However the massless flows only make sense if $`L>K`$. Both possibilities lead to a smaller number of equations than the relevant TBA equations. Therefore the new NLIEs are a computationally more efficient way of calculating the finite size effects of the FSSG models. In the next section we take a closer at the equations and obtain some results confirming the above conjectures.

## 3 Testing the NLIEs

We begin by examining the behaviour of the individual functions at $`r=0`$ and $`r=\mathrm{}`$. By differentiating the first equation of (13) with respect to $`\theta `$ we see the function $`i\mathrm{ln}(y_1)(\theta )`$ becomes approximately constant in the region $`|\theta |\mathrm{ln}(2/r)`$ as $`r0`$. For $`\theta \mathrm{ln}(2/r)`$ it behaves as $`e^{r\theta }`$, and for $`\theta \mathrm{ln}(2/r)`$ as $`e^{r\theta }`$. Therefore $`y_1`$ is a complex constant, which we denote $`x_1^0`$, inside the the central region $`\mathrm{ln}(2/r)\theta \mathrm{ln}(2/r)`$. The real and imaginary parts of $`y_1`$ oscillate outside this region but because $`\eta `$ is nonzero $`y_1(\theta i\eta )`$ and $`y_1^1(\theta +i\eta )`$ ultimately tend to zero. The functions $`y_2,y_3,\mathrm{},y_L`$ behave exactly as if they were TBA Y-functionsin standard notation $`Y_a(\theta )=\mathrm{exp}[\epsilon _a(\theta )]`$: they are symmetric in $`\theta `$ and real for real $`\theta `$. From the behaviour of $`y_1`$ it follows that in the ultraviolet limit they become constants $`x_a^{},x_a^0,x_a^+`$ in the regions $`\theta \mathrm{ln}(2/r),\mathrm{ln}(2/r)\theta \mathrm{ln}(2/r)`$ and $`\theta \mathrm{ln}(2/r)`$ respectively, while interpolating smoothly in between. A similar conclusion can be reached for the stationary solutions of the massless equations, apart from a small change in the behaviour of $`y_1`$ as it now becomes a nonzero constant in the region $`\theta \mathrm{ln}(2/r)`$. It is simple to check this picture by plotting some numerical solutions to the NLIEs. The exact values of the stationary solutions, derived from the NLIEs in the usual way , are summarised for both massive and massless cases in table 1.

Using the standard methods  we find the ultraviolet effective central charge takes the form

$`c_{\mathrm{eff}}(0)`$ $`=`$ $`{\displaystyle \frac{3}{\pi ^2}}[2\left({\displaystyle \frac{1}{1+\overline{x}_1^0}}\right)\left({\displaystyle \frac{1}{1+\overline{x}_1^+}}\right)\left({\displaystyle \frac{1}{1+\overline{x}_1^{}}}\right)`$ (19)
$`+{\displaystyle \underset{a=1}{\overset{L}{}}}2\left({\displaystyle \frac{1}{1+x_a^0}}\right)\left({\displaystyle \frac{1}{1+x_a^+}}\right)\left({\displaystyle \frac{1}{1+x_a^{}}}\right)`$
$`{\displaystyle \frac{i\pi \alpha }{2}}(\mathrm{ln}{\displaystyle \frac{(1+x_1^0)^2}{(1+x_1^{})(1+x_1^+)}}\mathrm{ln}{\displaystyle \frac{(1+\overline{x}_1^0)^2}{(1+\overline{x}_1^{})(1+\overline{x}_1^+)}})],`$
where $``$ is Rogers dilogarithm function

$$(x)=\frac{1}{2}_0^x𝑑t\frac{\mathrm{ln}1t}{t}+\frac{\mathrm{ln}t}{1t}.$$
(20)
To evaluate (19) we use the standard properties and known sum rules of the dilogarithm (see for example )

$$\frac{6}{\pi ^2}\underset{a=2}{\overset{L}{}}\left(\frac{\mathrm{sin}^2[\pi /(L+2)]}{\mathrm{sin}^2[\pi a/(L+2)]}\right)=\frac{2(L1)}{L+2},\frac{6}{\pi ^2}(1)=1,$$
(21)
and one further sum rule which involves both real and complex arguments

$$\frac{6}{\pi ^2}\left(\left(\frac{1}{1+\overline{x}_1^0}\right)+\underset{a=1}{\overset{L}{}}\left(\frac{1}{1+x_a^0}\right)\right)=1.$$
(22)
Numerical calculations confirm that (22) holds for arbitrary real numbers $`\xi `$ and $`\alpha `$, $`L`$ is by definition an integer. The ultraviolet effective central charge predicted from the NLIEs takes the form

$$c_{\mathrm{eff}}(0)=\frac{3L}{L+2}\frac{6\alpha ^2\xi L}{L+\xi }.$$
(23)
Does this reproduce the expected results? Yes, it does: with $`\xi `$ related to $`\beta ^2`$ as in (2) and $`\alpha =0`$, the formula predicts $`c_{\mathrm{eff}}(0)=3L/(L+2)`$, whereas the choice (16) leads to the effective central charge of the unitary coset models $`_{L,K}`$:

$$c_{\mathrm{eff}}(0)=\frac{3L}{L+2}\frac{6L}{(L+K+2)(K+2)},$$
(24)
and the possibility (17) yields the effective central charge for the nonunitary coset models

$$c_{\mathrm{eff}}(r)=\frac{3L}{L+2}\frac{6L}{pq}.$$
(25)
If we use the massless driving terms, the individual parts of (19) combine in a different way, but the final result is again (23). Since we have fixed the parameters $`\xi =K+2`$ and $`\alpha =\pm 1/(K+2)`$ we actually obtain (24). We need to check the infrared limit to be sure the NLIEs really do describe these massless flows. As $`R\mathrm{}`$ the contribution to $`c_{\mathrm{eff}}(\mathrm{})`$ from $`y_{L+\xi +3},\mathrm{},y_L`$ becomes vanishingly small. The remaining functions have stationary solutions given by the massless ones in table 1 but with $`L`$ replaced by $`LK`$. Plugging these values into (19) gives

$$c_{\mathrm{eff}}(\mathrm{})=\frac{3L}{L+2}\frac{6L}{(K+2)(KL+2)},$$
(26)
which is exactly the ultraviolet effective central charge of $`_{L,KL}`$.

Before obtaining the conformal dimension of the perturbing operator we comment on a surprising relation between the NLIEs at the coset points and the associated TBA equations. The massive $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ perturbation of $`_{L,K}`$ has a TBA system based on $`(A_1A_{L+K1})_L`$. The symmetry between $`L`$ and $`K`$ means that we can either take a NLIE with $`L`$ functions (below denoted $`y_a`$) and $`\xi =K+2`$, or one with $`K`$ functions (denoted $`\stackrel{~}{y}_a`$) and $`\xi =L+2`$, in either case $`c_{\mathrm{eff}}(r)`$ is exactly the same. Curiously the magnonic parts of both NLIEs are closely related to the TBA functions. Numerical evidence indicates

$$\mathrm{exp}(f_a)=\{\begin{array}{cc}y_{La+1}& a=1,\mathrm{},L1\hfill \\ \stackrel{~}{y}_{aK+1}& a=L+1,\mathrm{},L+K1\hfill \end{array}.$$
(27)
The functions with nonzero driving term, $`\mathrm{exp}f_L,y_1,\stackrel{~}{y}_1`$, are not equal, but the integral equations imply they must be related. One such relation is

$`{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑\theta ^{}\chi (\theta \theta ^{})\mathrm{ln}(1+e^{f_{L2}(\theta ^{})})`$ $`=`$ (28)
$`{\displaystyle _{𝒞_1}}𝑑\theta ^{}\chi (\theta \theta ^{}\frac{i\pi }{2})\mathrm{ln}(1+y_1(\theta ^{})){\displaystyle _{𝒞_2}}𝑑\theta ^{}\chi (\theta \theta ^{}\frac{i\pi }{2})\mathrm{ln}(1+y_1^1(\theta ^{})).`$
The question of the conformal dimension of the perturbing operator, predicted as usual from the periodicity of the functions entering the NLIEs , can be partially answered via the relation described above with the coset TBA systems. For $`\xi =K+2`$ and $`\alpha =\pm 1/(K+2)`$ the periodicity of the TBA-like functions $`y_2,\mathrm{},y_L`$ inferred from that of the TBA systems for the $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ perturbation of the coset $`_{L,K}`$ is $`i\pi (\xi +L)`$ . It is not expected to depend on $`\alpha `$ and so one only needs to check what happens at other values of $`\xi `$. Doing so numerically leads us to believe all of the entire functions $`y_a(\theta )`$ have periodicity $`i\pi (\xi +L)`$, and therefore a Laurent expansion in powers of $`\mathrm{exp}(2/(\xi +L)\theta )`$. The quantities $`\mathrm{ln}(1+y_a(\theta ))`$ appearing in the kernels of the NLIEs have essentially the same shape as the quantities $`\mathrm{ln}(1+\mathrm{exp}(\epsilon _a(\theta ))`$ that appear in a pure TBA system, that is a central plateau which quickly goes to zero for $`\theta `$ outside the region $`|\theta |\mathrm{ln}(2/r)`$. Therefore we may employ Zamolodchikov’s argument  to predict that $`c_{\mathrm{eff}}(r)`$ expands in powers of $`r^{4/(\xi +L)}`$. Note that for some $`\alpha `$, $`x_1^0=0`$ and Zamolodchikov’s argument acquires a small modification, as explained in, with the result that $`c_{\mathrm{eff}}(r)`$ instead expands in powers of $`r^{2/(\xi +L)}`$. However for the choices of $`\alpha `$ needed here $`x_1^0`$ is never identically zero. If the perturbing operator is odd we would expect an expansion in $`r^{2(1\mathrm{\Delta })}`$ leading to the prediction $`\mathrm{\Delta }=12/(\xi +L)`$. Inserting the appropriate value of $`\xi `$ (17) we find $`\mathrm{\Delta }=12/(K+L+2),`$ which exactly matches the conformal dimension of $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ for both unitary and nonunitary (with $`K=Lp/(qp)`$) cosets. On the other hand if the perturbing operator is even we would expect an expansion in $`r^{4(1\mathrm{\Delta })}`$ which implies $`\mathrm{\Delta }=11/(\xi +L)`$. Writing $`\xi `$ in terms of $`\beta ^2`$ using (2) we find the conformal dimension of the FSSG models: $`\mathrm{\Delta }=\beta ^2/8\pi +11/L`$.

In the infrared limit the groundstate energy of the massive model expands in powers of $`e^{MR}`$ . Adapting the arguments of  we obtain from the NLIEs

$$E(R)4\mathrm{cos}\frac{\pi }{L+2}\mathrm{cos}\pi \alpha \frac{d\theta }{2\pi }\mathrm{cosh}\theta e^{MR\mathrm{cosh}\theta }.$$
(29)
With the choices of $`L`$ and $`\alpha `$ as given in the text this agrees perfectly with TBA results obtained for the perturbed cosets $`_{L,K}+\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ , the $`\mathrm{\Phi }_{13}`$ perturbation of the $`N=2`$ supersymmetric minimal models  and the $`SU(2)`$ Gross-Neveu model .

We now test the NLIEs numerically, comparing against the effective central charge found using various TBA equations . These TBA equations correspond either to $`_{L,K}+\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ for integer $`K`$, or to the FSSG model at the specific coupling given in (5). We also compare against the magnonic TBA system based on $`(A_1T_n)_L`$, identified in  with $`_{L,n3/2}+\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$. However this is only true at $`L=1`$ and for general $`L`$ these TBA equations correspond to $`_{L,nL1/2}+\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$. As no data for the chosen cases has previously been reported in the literature we solved both the TBA and NLIE equations numerically, using an iterative method. We checked stability of the solution by varying the number of iterations and the integration step size. The results are shown in table 2. The tables include a diagrammatic version of the TBA system with a solid circle representing the massive node. We also reproduce one massless flow whose TBA diagram has crosses to indicate the nodes on which the massless driving terms are placed. In all cases the agreement is very good.

We can also enquire to which models the NLIEs based on $`(A_1A_L)_j^{[1]}`$ correspond. If we set the $`j^{\mathrm{th}}`$ driving term to $`r\mathrm{cosh}\theta `$ and all others to zero, the ultraviolet effective central charge becomes

$$c_{\mathrm{eff}}(0)=\frac{3(Lj+1)}{Lj+3}\frac{6\xi \alpha ^2(Lj+1)}{(\xi +L)(\xi +j1)}.$$
(30)
For $`\xi =K+3j`$ and $`\alpha =\pm 1/\xi `$ we rediscover the massive $`\varphi _{\mathrm{id},\mathrm{id},\mathrm{adj}}`$ perturbation of a familiar model: $`_{Lj+1,K}`$.

## 4 $`N=2`$ supersymmetric sine-Gordon model

Recent interest in the $`N=2`$ supersymmetric sine-Gordon model  has arisen in two different contexts, both related to string theory. Though little is currently known about the excited state spectrum of the $`N=2`$ sine-Gordon model it is predicted to provide nontrivial information on the spectrum of states of closed superstrings in certain backgrounds . Also, boundary versions of the $`N=2`$ SSG model (see for example ) are relevant for open string theory.

The S-matrix of the $`N=2`$ SSG model has the usual factorised form

$$S_{SG}^{N=2}(\stackrel{~}{\beta })=S^{N=2}S_{\mathrm{SG}}(\beta ).$$
(31)
The supersymmetric piece $`S^{N=2}`$ is the sine-Gordon S-matrix at its $`N=2`$ supersymmetric point $`\beta ^2/8\pi =2/3`$. The bosonic piece is the sine-Gordon S-matrix at coupling $`\beta `$. The $`N=2`$ coupling constant $`\stackrel{~}{\beta }`$ is related to the sine-Gordon coupling via

$$\stackrel{~}{\beta }^2=\beta ^2/(1\beta ^2/8\pi )$$
(32)
TBA equations derived from the $`S^{N=2}`$ S-matrix are encoded on $`(A_1A_3)_2`$ . As noted in , the TBA equations for the $`N=2`$ SSG model can, in theory, be obtained by gluing the $`(A_1A_3)_2`$ TBA system to that of the sine-Gordon model at coupling $`\beta `$. In practice this has been done only for $`\stackrel{~}{\beta }^2/8\pi =k+1`$, since at these points the SG TBA reduces to a simple Dynkin diagram: $`D_{k+2}`$. A second set of simple TBA equations exist for some related models: the minimal $`N=2`$ conformal field theories perturbed by the superfield $`\mathrm{\Phi }_{13}`$ .

To describe the groundstate energy of the $`N=2`$ SSG model at any coupling we replace the complicated sine-Gordon TBA equations with a single nonlinear integral equation, as we did above for the FSSG models. We propose equations similar to (13,14,15) but of type $`(A_1D_3)_1^{[1]}`$ with $`\alpha =0`$, rather than $`(A_1A_L)_1^{[1]}`$. The nodes of $`D_n`$ are labelled such that $`1`$ is the first node of the tail, and $`n,n1`$ label the fork nodes.The stationary values $`x_a^0`$ are all infinite, which suggests that numerical solutions for very small $`r`$ will converge extremely slowly, if at all (this is equally true for any TBA equations with infinite stationary values). However combined with $`x_2^\pm =x_3^\pm =1`$ and $`x_1^\pm =0`$ we find the correct ultraviolet central charge of $`3`$. The leading logarithmic corrections to this value have been proposed for some of the models in . It is simple to check the proposed NLIEs have the expected infrared limit.

If the mass is moved onto either of the magnonic nodes so that the NLIEs are of type $`(A_1D_3)_2^{[1]}`$, we find the effective central charge of the sausage models $`SST_\lambda ^+`$ . Again the possibly infinite set of TBA equations are in principle known. The simplest, occurring at $`\lambda =1/M`$, are represented by an extended $`D_M`$ Dynkin diagram . A match is found with the NLIEs if we set $`\xi =M2`$ and $`\alpha =0`$. Alternatively if the twist $`\alpha `$ is set to $`\pm 1/(M2)`$ we find the $`\psi _1\overline{\psi }_1+\psi _1^{}\overline{\psi }_1^{}`$ perturbation of the $`Z_{M2}`$ models, usually denoted $`H_{M2}^{(0)}`$ . The models denoted $`H_{M2}^{(\pi )}`$ refer to the massless perturbation of the $`Z_{M2}`$ models. The flow of the effective central charge from these models to the minimal models $`(M1,M)`$ can be reproduced within the NLIEs provided the magnonic nodes are given driving terms of the form $`\frac{1}{2}re^\theta `$ and $`\frac{1}{2}re^\theta `$. Upon resetting $`\alpha `$ to zero we find the effective central charge of the ‘massless’ sausage models $`SST_{1/M}^{}`$.

Table 3 shows a numerical comparison of results found using this second set of NLIEs against some of the simplest TBA systems.

## 5 Summary

The SS model  is a two-parameter family of integrable models whose name follows from the S-matrix, which is a direct product of two sine-Gordon S-matrices. All of the models discussed above can be found as an SS model at specific values of the parameters . We leave for future investigation the possibility of gluing two sine-Gordon NLIEs together to describe the groundstate energy of this family of models.

Motivated by the existing results for the sine-Gordon model both sets of new NLIE should be generalised to encode the full finite size spectrum of the associated models. Furthermore, it is curious that the TBA driving term trick to find equations for the massless flows also works for the NLIEs; though only for the unitary coset models. A rather different method was used in  for the sine-Gordon model, also , to obtain NLIEs describing massless flows at any value of $`\beta ^2`$. It remains to be seen whether a similar idea can be used on the new NLIEs. Finally it would be extremely interesting to develop the finite size effects of the fractional sine-Gordon models defined on a half-line.

Acknowledgments – I’m grateful to Zoltan Bajnok, Andreas Klümper, Paul Pearce, Francesco Ravanini and Gerard Watts for useful discussions, and to Patrick Dorey and Roberto Tateo for useful suggestions and discussions. I thank the EPSRC for a Postdoctoral Fellowship and the Department of Mathematics and Statistics, University of Melbourne for hospitality.