## Free energy from forward scattering in 1+1d

## Daniel Schubring

## Initiative for the Theoretical Sciences

The Graduate Center, CUNY

365 Fifth Ave, New York, NY 10016, USA

The free energy, or equivalently the ground state energy in finite volume, may be calculated from forward scattering amplitudes using a formula due to Dashen, Ma, and Bernstein [1]. However a naive treatment leads to singularities when considering the scattering of three or more particles. It is shown in detail how the approach can be applied to multi-particle scattering in various massive scalar theories in 1+1d, with or without integrability. The results for the sinh-Gordon, Lieb-Liniger, and O N ( ) non-linear sigma models are compared to exact results. It is shown how bound states can be considered in this approach by considering the attractive Lieb-Liniger model.

## Contents

| 1   | Introduction                          | 1                                                       |   2 |
|-----|---------------------------------------|---------------------------------------------------------|-----|
|     | 1.1                                   | Extending the L¨scher u formula via Beth-Uhlenbeck      |   3 |
|     | 1.2                                   | Forward scattering singularities . . . . . . . . . . .  |   4 |
|     | 1.3                                   | Outline . . . . . . . . . . . . . . . . . . . . . . . . |   6 |
| 2   | The DMB formula and integrable models | The DMB formula and integrable models                   |   6 |
|     | 2.1                                   | Derivation of the DMB formula . . . . . . . . . . .     |   6 |
|     | 2.2                                   | Integrable models: Sinh-Gordon and Lieb-Liniger .       |   8 |
|     | 2.3                                   | Applying DMB to integrable models . . . . . . . .       |   9 |
|     | 2.4                                   | Extending the argument to multiple particles . . .      |  11 |
| 3   | Cutting thermal field theory diagrams | Cutting thermal field theory diagrams                   |  12 |
|     | 3.1                                   | Single vertex corrections . . . . . . . . . . . . . . . |  12 |
|     | 3.2                                   | Cutting a general diagram . . . . . . . . . . . . . .   |  14 |
|     | 3.3                                   | The bubble and melon diagrams . . . . . . . . . .       |  16 |
|     | 3.3.1                                 | The bubble diagram . . . . . . . . . . . . .            |  16 |
|     | 3.3.2                                 | The melon diagram: two-body part . . . . .              |  17 |
|     | 3.3.3                                 | The melon diagram: three-body part . . . .              |  18 |

| 4 On-shell contributions and large N theories   | 4 On-shell contributions and large N theories   | 4 On-shell contributions and large N theories      |   19 |
|-------------------------------------------------|-------------------------------------------------|----------------------------------------------------|------|
|                                                 | 4.1                                             | Diagrams in the large N theory . . . . . . . .     |   20 |
|                                                 | 4.2                                             | On-shell contributions from the DMB formula        |   22 |
|                                                 | 4.2.1                                           | DMB calculation of the bubble diagram              |   23 |
|                                                 | 4.2.2                                           | Higher-order on-shell contributions . .            |   24 |
| 5 The TBA and bound states                      | 5 The TBA and bound states                      | 5 The TBA and bound states                         |   25 |
|                                                 | 5.1                                             | Perturbative expansion of the TBA . . . . . .      |   25 |
|                                                 |                                                 | 5.1.1 Perturbative expansion of the kernel .       |   26 |
|                                                 |                                                 | 5.1.2 The bubble and melon from the TBA            |   27 |
|                                                 | 5.2                                             | The attractive Lieb-Liniger model . . . . . .      |   28 |
|                                                 | 5.3                                             | Bound states in the DMB approach . . . . . .       |   30 |
| 6                                               | Discussion                                      | Discussion                                         |   31 |
| A                                               | More on the melon diagram                       | More on the melon diagram                          |   32 |
|                                                 | A.1                                             | Zero- and one-body parts . . . . . . . . . . .     |   33 |
|                                                 | A.2                                             | Cutting the melon diagram . . . . . . . . . .      |   33 |
|                                                 | A.3                                             | Direct calculation of the three-body part . . .    |   34 |
| B                                               | More on large N diagrams                        | More on large N diagrams                           |   35 |
|                                                 | B.1                                             | An expansion of x s . . . . . . . . . . . . . . .  |   35 |
|                                                 | B.2                                             | Calculation of f α 4 . . . . . . . . . . . . . . . |   36 |
| C                                               | The Lieb-Liniger model                          | The Lieb-Liniger model                             |   38 |

## 1 Introduction

One of the most basic physical quantities that may be calculated from a quantum field theory is the free energy as a function of temperature. It is well-understood how to calculate this perturbatively using compactified Euclidean time, an approach referred to here as thermal field theory. The free energy should in principle encode much information about the zero-temperature dynamics of the theory, but since the size β of the thermal circle is involved in each propagator of a diagram, the thermal field theory approach mixes statistics with dynamics in an essential way.

Much better in this respect is the thermodynamic Bethe ansatz (TBA) [2]. The zero-temperature dynamics are encoded in kernels K which involve derivatives of the logarithm of the two-body Smatrix (see e.g. equation (2.14)). These kernels then appear in an integral equation to determine the free energy. However, this approach only works for integrable systems.

There is a less widely known third approach to calculate the free energy introduced by Dashen, Ma, and Bernstein (DMB) [1] which extends a formula by Beth and Uhlenbeck [3]. This approach explicitly separates the zero temperature dynamics from the temperature by writing the partition

function in terms of the density of states ρ E ( ),

$$Z = \int d E e ^ { - \beta E } \rho ( E ).$$

The density of states may be calculated in terms of forward scattering amplitudes at zero-temperature, even for non-integrable systems. In fact the expression for ρ involves a derivative of a logarithm of an S-matrix operator, and it is suggestive of the kernel K in the TBA approach,

$$\rho ( E ) = \rho ^ { ( 0 ) } ( E ) + \frac { 1 } { 2 \pi i } \frac { d } { d E } \text{Tr} \log \hat { S } ( E ).$$

Here ρ (0) refers to the non-interacting density of states, and the operator ˆ will be defined carefully S in Sec 2 where the derivation of the DMB formula is reviewed.

The calculation of ρ in the DMB approach has certain subtleties, and the lack of explicit examples in the literature has perhaps hindered its adoption. In this paper, in order to clearly illustrate the DMB approach we will consider certain massive scalar field theories in two spacetime dimensions for which we can compare with calculations in the thermal field theory and TBA approaches.

## 1.1 Extending the L¨scher formula via Beth-Uhlenbeck u

A reconsideration of the DMB approach [1, 4, 5, 6], as well as related approaches to calculating the self-energy and higher correlation functions [7, 8, 9, 10], may ultimately be useful for extracting information on multi-particle scattering amplitudes from lattice simulations. Note that in 1+1d the free energy may alternatively be understood as a calculation of the ground state energy as a function of finite volume L . A better understanding of how to relate non-integrable scattering data to the energy levels as a function of L is relevant for the effective theory of Yang-Mills flux tubes [11, 12, 13]. There has been great success in using the TBA to relate the two-body S-matrix on the worldsheet of the flux tube to lattice calculations of the energy levels of the flux tube as a function of L [14, 15]. Even so, it is understood that the worldsheet theory for both 3D and 4D Yang-Mills is not actually integrable [16, 17].

There is an alternate approach due to L¨scher [18, 19] which is well-known in this context, and u which can be applied to the two-body energy levels in non-integrable systems. As will be explained below, the present investigation is an extension of the infinite volume limit of the L¨scher formula u to multi-particle scattering states (see also [20]). This is not yet directly applicable to the effective string theory program since we are dealing with massive particles, but it is clearly of value to better understand how multi-particle scattering affects the finite size dependence of the energy levels in a simple case.

The L¨ uscher formula in two dimensions is essentially the same thing as the two-particle asymptotic Bethe ansatz (ABA)

$$2 \pi i n = i k _ { n } L + 2 \log S ( k _ { n } ),$$

where n is an integer specifying the energy level 1 , k n is the relative momenta between the particles, and S k ( n ) is the two-body S-matrix. Strictly speaking this formula is valid for massive particles in the limit of large mL , and any additional terms which vanish in this limit are suppressed in this section. Since the energy levels are closely spaced in this limit, it makes sense to instead consider the density of two-particle states ρ (2) ,

$$\rho _ { ( 2 ) } ( E ) = \left ( \frac { d E } { d k } \right ) ^ { - 1 } \frac { 1 } { k _ { n + 2 } - k _ { n } } = \frac { L } { 4 \pi } \left ( \frac { d E } { d k } \right ) ^ { - 1 } + \frac { 1 } { 2 \pi i } \frac { d } { d E } \log S ( E ).$$

Here we are taking fixed center of mass momentum, so that the total energy E is a single-variable function of the relative momentum k , and the two body S-matrix is rewritten in terms of E . The first term is just the density of two-particle states in the non-interacting theory, and the remaining term δρ (2) is

$$\delta \rho _ { ( 2 ) } ( E ) = \frac { 1 } { 2 \pi i } \frac { d } { d E } \log S ( E ).$$

This formula was originally derived in ordinary quantum mechanics by Beth and Uhlenbeck [3].

Of course, the density of states also has contributions from states with more than two particles. For integrable systems we may continue to use the ABA in order to calculate the density of states [21]. But such an approach is not directly applicable to non-integrable systems.

Instead we can extend the Beth-Uhlenbeck formula following DMB [1]. The DMB formula (1.1) is superficially identical to the Beth-Uhlenbeck formula, but as will be discussed further in Sec 2 the operator ˆ ( S E ) is an S-matrix operator that may be applied to arbitrary multiparticle states, and the E argument need not be equal to the energy of the states involved in the trace.

## 1.2 Forward scattering singularities

To see roughly how (1.1) is to be interpreted, write the non-trivial part of the S-matrix operator as usual in terms of a T-matrix operator with an energy delta function factored out,

$$\hat { S } ( E ) = 1 - 2 \pi i \delta ( E - H _ { 0 } ) \hat { T } ( E ).$$

The factor log ˆ ( S E ) may be expanded to first order in powers of ˆ , T and after integrating by parts in E the partition function may be written as

$$\delta Z = - \beta \, \int d \alpha \, e ^ { - \beta E _ { \alpha } } T _ { \alpha \alpha } + \mathcal { O } ( T ^ { 2 } ) \,.$$

Here the trace was taken over eigenstates α of H 0 with energy E α , and dα schematically represents the appropriate measure for the trace. T αα is the forward scattering amplitude with both an incoming and outgoing state of α . As usual, we may consider instead the connected part of the T-matrix, T c , and the free energy density 2 f = -( βL ) -1 log Z ,

$$\delta f = L ^ { - 1 } \, \int d \alpha \, e ^ { - \beta E _ { \alpha } } T _ { \alpha \alpha } ^ { c } + \mathcal { O } ( T ^ { 2 } ) \,.$$

1 Since the equation is written only in terms of the relative momentum, n varies in steps of 2 for fixed center of mass momentum.

2 We use L for spatial volume, although the following formulas apply to arbitrary dimension.

Figure 1: A three-body scattering amplitude which may go on-shell in the forward scattering case.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-1c08c5dd-6324-4ab6-a2b4-c656d26012a7.jpg)

Here already we see a potential issue. A connected T-matrix amplitude T c βα may involve various delta functions, and in the forward scattering limit β = α the argument of these delta functions may go to zero. One delta function will just be due to momentum conservation and it will lead to an overall factor of L to cancel the L -1 in (1.5). But additional delta functions may arise from internal propagators in scattering amplitudes which can go on-shell.

For instance the lowest-order three-particle scattering amplitude in a relativistic scalar field theory in Fig 1 is just proportional to the internal propagator

$$\frac { 1 } { ( E _ { 1 } + E _ { 2 } - E _ { 1 } ^ { \prime } ) ^ { 2 } - ( k _ { 1 } + k _ { 2 } - k _ { 1 } ^ { \prime } ) ^ { 2 } - m ^ { 2 } + i \epsilon },$$

√

where the E are on-shell energies E = k 2 + m 2 associated to the spatial momentum k of a given external leg. This propagator goes on-shell when k ′ 1 equals k 1 or k 2 . The propagator can be broken into a real principal part, and an imaginary delta function. These delta functions can not be ignored and in fact for integrable models they are necessary to ensure that the set of final momenta are equal to the initial ones (see the discussion in [22]). So the forward scattering amplitude T c αα will certainly involve delta functions with a vanishing argument, and these forward scattering singularities need to be interpreted for the DMB formula to make sense.

As we will show perturbatively in Sec 2.4, this issue is resolved by including the higher order powers of T in the expansion in (1.5). The delta function parts of the internal propagators in diagrams like Fig 1 are canceled, and only the principal part remains. This introduces some new questions since for three-particle scattering in integrable theories the principal part is taken to vanish. Nevertheless we will show that there are two subtle types of contributions to the free energy due to the principal part: on-shell contributions and collinear contributions .

On-shell contributions are due to forward scattering diagrams that would have had an onshell divergence, such as Fig 1 with k ′ 1 taken to be k 1 or k 2 . The principal part of the internal propagator vanishes exactly on-shell, but it is large slightly off-shell, and the slightly off-shell region will contribute since the energy delta function involved in the S-matrix operator (1.4) actually has a width that depends on iϵ .

The calculation of the on-shell contributions has previously been discussed to a certain extent for the scattering of three or four particles [4, 5], as well as more generally using a different approach

- [6]. In section 4 we give an arguably simpler calculation of the on-shell contributions for an arbitrary number of particles in relativistic theories from both the thermal field theory and DMB perspectives.

Besides the on-shell contributions, there are also collinear contributions which appear at lowest order due to forward scattering diagrams such as Fig 1 with k ′ 1 = k 3 , so that the internal propagator need not be on-shell. If k , k 1 2 , k 3 are all distinct then the principal part vanishes as per the usual argument. But since we are integrating over all momentum in the DMB formula, we need to consider the cases where the particle momenta become collinear. It turns out that the principal part doesn't vanish in a small iϵ dependent region around the collinear point k 1 = k 2 = k 3 , and this is actually a significant contribution to the free energy.

In the ϕ 4 theory there are two free energy diagrams at second order in the thermal field theory approach, which may be called the bubble and melon diagrams (see Fig 4). We will show that the melon diagram in 1+1d is entirely the result of collinear contributions. Note that the massive [23, 24] and massless [25] ϕ 4 theory have been investigated before in 3+1d via the DMB approach and the related approach of Barton [26], but the essential features of the 1+1d case were not seen. The 1+1d melon diagram is calculated in Sec 3.3 and compared to a perturbative expansion of the TBA in Sec 5.1. In Sec 5.2 it is shown how these collinear contributions are necessary in perturbation theory since there are bound state contributions to the free energy upon flipping the sign of the coupling constant.

## 1.3 Outline

A brief outline of this paper is as follows. In Sec 2 we will review the general derivation of the DMB formula and apply it to integrable models involving a single massive scalar such as sinh-Gordon. In Sec 3 we consider sinh-Gordon and the non-integrable ϕ 4 theory at second order in perturbation theory, and show how free energy diagrams in the thermal field theory approach are connected to two- or three-body forward scattering diagrams in the DMB approach.

In Sec 4 diagrams involving only on-shell contributions are calculated at arbitrary order in perturbation theory, and it is shown how these give the leading contribution to the free energy for the large N extension of the ϕ 4 model. In Sec 5 we investigate the treatment of bound states in the DMB approach by considering the non-relativistic Lieb-Liniger model with attractive coupling. We show how the collinear contributions at repulsive coupling are related to the contributions due to bound states for attractive coupling. In Sec 6 we discuss future investigations into the DMB approach.

## 2 The DMB formula and integrable models

## 2.1 Derivation of the DMB formula

First let us review some material from scattering theory. The Lippmann-Schwinger equation expresses the in and out eigenstates Ψ ± of the Hamiltonian H = H 0 + V in terms of the eigenstates

Φ of H 0 with the same eigenvalue E ,

$$\Psi ^ { \pm } = \Phi + \frac { 1 } { E - H _ { 0 } \pm i \epsilon } V \Psi ^ { \pm }.$$

This can be written in terms of a quantity G 0 that is usually referred to as the propagator in this context,

$$G _ { 0 } ( E ) \equiv \frac { 1 } { E - H _ { 0 } + i \epsilon }.$$

An interacting propagator G may also be defined,

$$G ( E ) \equiv \frac { 1 } { E - H + i \epsilon } = \frac { 1 } { 1 - G _ { 0 } ( E ) V } G _ { 0 } ( E ).$$

Now we can define a T -matrix operator so that the Lippmann-Schwinger equation can be written as V Ψ + = ˆ ( T E )Φ,

$$\hat { T } ( E ) \equiv V \frac { 1 } { 1 - G _ { 0 } ( E ) V }.$$

An S -matrix operator may also be defined in terms of T ,

$$\hat { S } ( E ) \equiv 1 + \left ( G _ { 0 } ( E ) - G _ { 0 } ( E ) ^ { \dagger } \right ) \hat { T } ( E )$$

$$= \left ( 1 - G _ { 0 } ( E ) ^ { \dagger } V \right ) ( 1 - G _ { 0 } ( E ) V ) ^ { - 1 }$$

The first line is just the usual relation between S and T since G 0 -G † 0 → -2 πiδ E ( -H 0 ) in the ϵ → 0 limit.

The ˆ ( S E ) operator is defined even outside the usual limits in which scattering theory applies, but in order to interpret it as the usual S -matrix it must be acting on an H 0 eigenstate Φ with the same eigenvalue E as its argument, and the infinite volume and ϵ → 0 limits must be taken. Keeping these caveats in mind, we will suppress the E argument except when necessary.

Now the DMB formula follows from some simple manipulations using the formulas above and the fact that ∂ ∂E G 0 ( E ) = -G 0 ( E ) 2 .

$$\text{the fact that } \frac { \partial _ { \partial E } G _ { 0 } ( E ) = - G _ { 0 } ( E ) ^ { 2 }. \\ \frac { 1 } { 2 \pi i } \frac { d } { d E } \text{Tr} \log \hat { S } & = \frac { 1 } { 2 \pi i } \frac { d } { d E } \text{Tr} \left [ \log \left ( 1 - G _ { 0 } ^ { \dagger } V \right ) - \log \left ( 1 - G _ { 0 } V \right ) \right ] \\ & = \frac { 1 } { 2 \pi i } \text{Tr} \left [ \frac { \left ( G _ { 0 } ^ { \dagger } \right ) ^ { 2 } V } { 1 - G _ { 0 } ^ { \dagger } V } - \frac { G _ { 0 } ^ { 2 } V } { 1 - G _ { 0 } V } \right ] \\ & = \frac { 1 } { 2 \pi i } \text{Tr} \left [ \left ( \frac { 1 } { 1 - G _ { 0 } ^ { i } V } - 1 \right ) G _ { 0 } ^ { \dagger } - \left ( \frac { 1 } { 1 - G _ { 0 } V } - 1 \right ) G _ { 0 } \right ] \\ & = - \frac { 1 } { 2 \pi i } \text{Tr} \left [ \left ( G - G ^ { \dagger } \right ) - \left ( G _ { 0 } - G _ { 0 } ^ { \dagger } \right ) \right ] \\ & = \text{Tr} \left [ \delta ( E - H ) - \delta ( E - H _ { 0 } ) \right ],$$

where the last line holds in the ϵ → 0 limit. This is now in the form of the difference between the density of states ρ E ( ) -ρ (0) ( E ) as required in the DMB formula (1.1).

In situations where there is a conserved particle number N in both the free and interacting theories it further makes sense to introduce a chemical potential µ and corresponding fugacity y ≡ e βµ , in terms of which the grand canonical partition function is

$$\delta Z = \sum _ { N = 2 } ^ { \infty } y ^ { N } \int d E e ^ { - \beta E } \frac { 1 } { 2 \pi i } \frac { d } { d E } \text{Tr} _ { ( N ) } \log \hat { S }.$$

In what follows fugacity dependence will be used as a counting device to single out the trace over N -particle states Tr ( N ) , and we will set y = 1 at the end. In the theories of massive particles considered here the term in the partition function proportional to y N will have an asymptotic temperature dependence on the order of ∼ e -Nβm . So the fugacity expansion may equivalently be understood as a low temperature expansion in powers of e -βm .

We will use a notation where a subscript in parenthesis denotes the order in the fugacity expansion. On the other hand, a superscript in parentheses denotes the order in a coupling constant ( c or g 2 ) expansion, and a δ prefix denotes the difference of a quantity from the free theory. Later on we will introduce an n -body expansion that takes into account particle statistics, and this will be indexed by a subscript without parentheses.

## 2.2 Integrable models: Sinh-Gordon and Lieb-Liniger

It should be clear that the DMB formula above is quite general since it simply follows from the definitions of the operators ˆ and T ˆ , S but to illustrate it in practice we consider the sinh-Gordon model (SG) [28], with Lagrangian

$$\mathcal { L } _ { S G } = \frac { 1 } { 2 } \left ( \partial \phi \right ) ^ { 2 } + \frac { \left ( m ^ { 2 } + \delta m ^ { 2 } \right ) } { g ^ { 2 } } \left ( \cosh \left ( g \phi \right ) - 1 \right ).$$

Here δm 2 is a counterterm set to ensure that m is the physical mass.

In many respects the non-relativistic limit of SG may be understood to be the Lieb-Liniger model (LL) [27], which involves bosons in ordinary quantum mechanics interacting through a repulsive two-body delta function potential

$$V _ { L L } ( x ) = \frac { c } { m } \delta ( x ).$$

The parameter g in the SG Lagrangian is related to the parameter c by

$$c = m \sin \frac { \pi g ^ { 2 } } { 8 \pi + g ^ { 2 } } = \frac { m g ^ { 2 } } { 8 } - \frac { m g ^ { 4 } } { 6 4 \pi } + \mathcal { O } ( g ^ { 6 } ) \,.$$

Both SG and LL are integrable and the exact two-body S-matrix is given by

$$S ( p ) = \frac { p - i c } { p + i c }.$$

In LL, p is taken to be the momentum difference between the particles, whereas in SG, p = m sinh θ with θ being the rapidity difference.

Since both models involve a single particle type with no internal indices and no bound states, the TBA takes the simplest possible form (see [32] for a review). The free energy density f is calculated as

$$f = - T \int \frac { d p } { 2 \pi } \log \left ( 1 + e ^ { - \epsilon ( p ) / T } \right )$$

where the pseudo-energy ϵ is determined by solving the TBA equation

$$\epsilon = E - \mu - T \log \left ( 1 + e ^ { - \epsilon / T } \right ) * K.$$

Here E is the single particle dispersion relation, and the ⋆ denotes convolution 3

$$f \ast g ( \theta ) \equiv \int d \theta ^ { \prime } f ( \theta ^ { \prime } ) g ( \theta - \theta ^ { \prime } ).$$

K is a kernel that is derived from the exact S-matrix. The kernels for LL and SG are respectively,

$$K _ { L L } = \frac { 1 } { 2 \pi i } \frac { d } { d p } \log S ( p ) = \frac { 1 } { \pi } \frac { c } { p ^ { 2 } + c ^ { 2 } },$$

$$K _ { S G } = \frac { 1 } { 2 \pi i } \frac { d } { d \theta } \log S \left ( m \sinh \theta \right ) = m \cosh \theta K _ { L L }.$$

The formal similarity of the kernels to the DMB formula for the density of states is not an accident, since both approaches will reduce to the Beth-Uhlenbeck formula (1.3) in the two-particle limit, as shown in the section below.

## 2.3 Applying DMB to integrable models

Since SG and LL are integrable and only involve one species of particle, this implies that arbitrary multiparticle states α in the momentum basis are eigenstates of the S-matrix operator

$$\hat { S } ( E _ { \alpha } ) | \alpha \rangle = S ( \alpha ) | \alpha \rangle.$$

Here ˆ refers to the operator appearing in the DMB formula, and S S α ( ) is a complex eigenvalue which ends up being factorizable into two-body S-matrices.

The normalization of the eigenstates is fixed to be compatible with the measure in the sense that ⟨ α β | ⟩ acts like a delta function,

$$\int d \alpha \langle \alpha | \beta \rangle f ( \alpha ) = f ( \beta ).$$

Now let us use these eigenstates as the basis for the trace in the DMB formula,

$$\text{Now let us use these eigenstates as the basis for the trace in the DMB formula,} \\ \delta Z & = \frac { \beta } { 2 \pi i } \int d E e ^ { - \beta E } \text{Tr} \log \hat { S } ( E ) \\ & = - \frac { \beta } { 2 \pi i } \sum _ { k = 1 } \frac { 1 } { k } \int d E e ^ { - \beta E } \int d \alpha \langle \alpha | \left ( 2 \pi i \delta ( E - H _ { 0 } ) \hat { T } ( E ) \right ) ^ { k } | \alpha \rangle \\ & = - \beta \sum _ { k = 1 } \frac { 1 } { k } \int d \alpha \, e ^ { - \beta E _ { \alpha } } \langle \alpha | \hat { T } ( E _ { \alpha } ) \left ( 2 \pi i \delta ( E _ { \alpha } - H _ { 0 } ) \hat { T } ( E _ { \alpha } ) \right ) ^ { k - 1 } | \alpha \rangle. \\ \boxed { ^ { 3 } } \text{The convolution is taken over the rapidity variable and the momentum variable in SG and LL. respectively.}$$

3 The convolution is taken over the rapidity variable and the momentum variable in SG and LL, respectively.

Here the last step of using one delta function to evaluate the E integral is somewhat questionable since the delta function is really G 0 -G † 0 for finite iϵ , and there may be problems with the vanishing iϵ limit in the T-matrix amplitudes. But for now let us proceed. This may be evaluated further by inserting an integration over a complete set of states between each δ E ( -H T 0 ) , and using the idea that

$$\langle \beta | 2 \pi i \delta ( E _ { \alpha } - E _ { \beta } ) \ddot { T } ( E _ { \alpha } ) | \alpha \rangle = ( 1 - S ( \alpha ) ) \, \langle \beta | \alpha \rangle.$$

The result is

$$\intertext { u c u u } & = - \beta \sum _ { k = 1 } ^ { 1 } \frac { 1 } { k } \int d \alpha d \beta _ { 1 } \dots d \beta _ { k - 1 } \, e ^ { - \beta E _ { \alpha } } \langle \alpha | \hat { T } ( E _ { \alpha } ) | \beta _ { 1 } \rangle \dots \langle \beta _ { k - 1 } | 2 \pi i \delta ( E _ { \alpha } - E _ { k - 1 } ) \hat { T } ( E _ { \alpha } ) | \alpha \rangle \\ & = \beta \int d \alpha \, e ^ { - \beta E _ { \alpha } } \left ( \frac { T ( \alpha ) } { 1 - S ( \alpha ) } \right ) \log S ( \alpha ), \\ & \quad.$$

where T α ( ) ≡ ⟨ α T | ˆ ( E α ) | α ⟩ is the on-shell forward scattering T-matrix. Note that considering (2.17), T α ( ) appears to be proportional to 1 -S α ( ) and the factor in parenthesis in (2.18) appears to just be a kinematic factor.

This formula is able to reproduce the correct free energy for two-particle states. A quick way to see this is to use a normalization of two-particle states such that dα = dEdP with E and P being the total energy and momentum. If it is valid to cancel the energy delta function from both sides of (2.17) then we can calculate T ,

$$2 \pi i T ( \alpha ) = \left ( 1 - S ( \alpha ) \right ) \frac { L } { 2 \pi }.$$

Then after integrating by parts, and rewriting in terms of the two-particle free energy density f (2) , (2.18) becomes

$$\delta f _ { ( 2 ) } = - \frac { T } { 2 \pi } \int d P d E e ^ { - \beta E } \frac { 1 } { 2 \pi i } \frac { \partial } { \partial E } \, \log S ( \alpha ) | _ { P } \,.$$

This is just the Beth-Uhlenbeck formula (1.3) with the center of mass motion taken into account.

To write it more explicitly in SG, transform the integration measure to rapidity variables and use the kinematic relation E 2 = P 2 +4 m 2 cosh 2 θ to rewrite the E derivative as a θ derivative,

$$\delta f _ { ( 2 ) } = - \frac { T } { 2 \pi } \frac { 1 } { 2! } \int d \theta _ { 1 } d \theta _ { 2 } E e ^ { - \beta E } \frac { 1 } { 2 \pi i } \frac { \partial } { \partial \theta } \log S ( \theta ) = - T \int \frac { d p _ { 1 } } { 2 \pi } e ^ { - \beta E _ { 1 } } \int d \theta _ { 2 } e ^ { - \beta E _ { 2 } } K _ { S G } ( \theta ). \quad ( 2. 2 0 )$$

In the last equality this has been written in terms of a convolution integral over the SG kernel (2.15), and it is a short step to see agreement with the TBA.

According to the formula for the grand canonical partition function (2.8), the above expression should match the O ( y 2 ) term in a fugacity expansion of free energy density (2.13). This expansion may be carried out by expanding the pseudo-energy ϵ in powers of y and solving the TBA equation (2.13) at each order. The result is

$$\delta f _ { ( 2 ) } = \left ( f - f ^ { ( 0 ) } \right ) _ { ( 2 ) } = - T \int \frac { d p } { 2 \pi } e ^ { - \beta E } \left ( e ^ { - \beta E } \ast \delta K \right ).$$

Here we have subtracted off the the free energy density of a free boson gas,

$$f ^ { ( 0 ) } = + T \int \frac { d p } { 2 \pi$$

and δK is the interacting part of the kernel,

$$\delta K _ { L L } ( p ) \equiv K _ { L L } ($$

Clearly the TBA result (2.21) matches the corresponding DMB result for SG (2.20), and a similar agreement may be seen for the LL model. Agreement of the DMB formula with the TBA at O ( y 2 ) has also been noted for the O N ( ) model [29].

## 2.4 Extending the argument to multiple particles

Now we would like to apply the seemingly general formula (2.18) to the scattering of three or more particles, but it is difficult to see how this could succeed. If T/ (1 -S ) is still just a kinematic factor, then since S is factorizable into two-body S-matrices the log S factor just breaks up into a disconnected sum. In fact we were too cavalier in the treatment of the iϵ terms in the derivation of (2.18), and the contribution due to multi-particle scattering depends sensitively on this.

The contribution to the free energy due to three-particle scattering may be investigated by expanding the TBA to the third order in fugacity,

$$\dots _ { r } \dots \dots \dots \dots \dots$$

This expression has also been derived by Kato and Wadati using an approach based on the asymptotic Bethe ansatz (ABA) [21].

As will be discussed more at the end of Sec 3.1, the first line of (2.24) involving a single appearance of δK is better understood as a correction to the two-body contribution (2.21) due to particle exchange. The two-body contribution δf 2 to all orders in fugacity is

$$\delta f _ { 2 } = + T \int \frac { d p } { 2 \$$

The second line of (2.24) involving two factors of δK represents genuine three-particle scattering. Contributions due to connected scattering of even more particles may be derived by continuing the fugacity expansion of the TBA, or using the ABA approach of [21]. These expressions become increasingly complicated, and there is seemingly no way to derive them from the naive formula (2.18).

To better understand how multi-particle scattering is treated in the DMB formula, we consider the second order in perturbation theory, which is the lowest order for which three-particle scattering is relevant. The operator log ˆ may be expanded in powers of S V ,

$$\begin{array} {$$

Note that the internal propagator of ˆ T (2) is

$$G _ { 0 } = \frac { 1 } { E - H _ { 0 } + i \epsilon } = \mathcal { P } G _ { 0 } - i \pi \delta ( E - H _ { 0 } ).$$

For integrable theories, as long as all momenta of the particles are distinct (as is assumed in the ABA), the net contribution from the P G 0 terms vanishes when summed over all diagrams. Rather it is the delta function term which leads to an ˆ operator which is diagonal in the momentum basis, S which is assumed in the derivation of (2.18). However (2.26) asserts that only P G 0 contributes in the DMB formula, and not the delta function part.

One way in which P T (2) can be non-zero is if E is not strictly set to be equal to the energy of the external particles in the scattering amplitude. This is possible since the 'delta function' ( G 0 -ˆ G 0 ) has some width in terms of ϵ . These type of contributions will be referred to as on-shell contributions, since they involve diagrams for which the denominator of G 0 vanishes. On-shell contributions will be treated in DMB approach to arbitrary order in perturbation theory in Sec 4.

The other way in which P T (2) can be non-zero is if the particle momenta are collinear. This is not taken to be the case in the ABA, but it is encountered in the trace over free particle states in the DMB formula. We will treat the collinear contributions at second order in perturbation theory in Sec 3.3 and the associated Appendix A. As discussed in Sec 5.2, the collinear contributions are closely connected to bound state contributions, and they are expected to be important at higher order in perturbation theory as well.

## 3 Cutting thermal field theory diagrams

The DMB formula gives contributions to the free energy in terms of forward scattering amplitudes. In this section we will show how the forward scattering diagrams arise by cutting lines of vacuum bubble diagrams in thermal field theory. As a concrete example we will consider the second-order diagrams in sinh-Gordon (2.9) and its non-integrable truncation, the ϕ 4 theory,

$$\mathcal { L } = \frac { 1 } { 2 } \left ( \partial \phi \right ) ^ { 2 } + \frac { m ^ { 2 } + \delta m ^ { 2 } } { 2 } \phi ^ { 2 } + \frac { m ^ { 2 } g ^ { 2 } } { 4! } \phi ^ { 4 }.$$

## 3.1 Single vertex corrections

The simplest free energy corrections to consider involve a single vertex proportional to ϕ 2 n , which is contracted into a vacuum bubble diagram involving n tadpoles. See Fig 2 for the n = 2 and n = 3 cases, which are relevant in sinh-Gordon up to second-order. Besides the diagram involving the ϕ 2 n vertex, there will also be diagrams at the same order involving counterterm vertices proportional to δm 2 . For instance the total free energy density f (1) at first order is

$$f ^ { ( 1 ) } = \frac { m ^ { 2 } g ^ { 2 } } { 4! } 3 \, ( I ^ { 2 } - I _ { 0 } ^ { 2 } ) + \frac { \delta m ^ { 2 } } { 2 } ( I - I _ { 0 } ).$$

This is written in terms of the tadpole integral I ,

$$I \equiv \int _ { \beta } \frac { d ^ { 2 } p } { ( 2 \pi ) ^ { 2 } } \frac { 1 } { p ^ { 2 } + m ^ { 2 } } = \frac { 1 } { \beta } \sum _ { n } \int \frac { d k } { 2 \pi } \frac { 1 } { \omega _ { n } ^ { 2 } + k ^ { 2 } + m ^ { 2 } },$$

Figure 2: The free energy correction due to a ϕ 4 vertex and an associated two-point counterterm (top), and a ϕ 6 vertex and its associated counterterms (bottom).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-15a1477a-9824-477b-96f8-95a2a8510edf.jpg)

where the β subscript shorthand indicates that Euclidean time is compactified, and the integral may be written explicitly in terms of Matsubara frequencies ω n = 2 πn/β . This integral may be written as a sum of a zero-temperature part I 0 and a thermal part that involves the bosonic occupation number n k , which is defined as

$$n _ { k } ( \beta ) \equiv \frac { 1 } { e ^ { \beta E _ { k } } - 1 }, \quad E _ { k } \equiv \sqrt { k ^ { 2 } + m ^ { 2 } }.$$

To make this decomposition, the sum over ω n is transformed to an integral over ω by using the Poisson resummation formula

$$I & = \sum _ { l } \int \frac { d k } { 2 \pi } \frac { d \omega } { 2 \pi } \frac { e ^ { - i l \beta \omega } } { \omega ^ { 2 } + k ^ { 2 } + m ^ { 2 } } \\ & = I _ { 0 } + \int \frac { d k } { 2 \pi } \frac { n _ { k } ( \beta ) } { E _ { k } }. \\. & \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

̸

The l = 0 term is the same as the zero-temperature tadpole integral, and this is referred to as I 0 . In the l = 0 terms, the ω integral has been performed, and the geometric series summed.

Now noting that the mass counterterm at this order is δm 2 = -λ 2 I , 0 it is seen that all appearances of I 0 cancel in (3.2)

$$f ^ { ( 1 ) } & = \frac { m ^ { 2 } g ^ { 2 } } { 8 } \left ( \int \frac { d k } { 2 \pi } \frac { n _ { k } ( \beta ) } { E _ { k } } \right ) ^ { 2 } \\ & = \frac { 1 } { L } \, \frac { 1 } { 2! } \int \frac { d k _ { 1 } } { 2 \pi } \frac { 1 } { 2 E _ { 1 } } \frac { d k _ { 2 } } { 2 \pi } \frac { 1 } { 2 E _ { 2 } } \frac { e ^ { - \beta ( E _ { 1 } + E _ { 2 } ) } } { ( 1 - e ^ { - \beta E _ { 1 } } ) \, ( 1 - e ^ { - \beta E _ { 2 } } ) } m ^ { 2 } g ^ { 2 } L.$$

If we momentarily ignore the two factors of (1 -e -βE ) -1 , this is clearly in the form of the equation (1.5),

$$f ^ { ( 1 ) } = \frac { 1 } { L } \int d \alpha \, e ^ { - \beta E _ { \alpha } } T _ { \alpha \alpha } ^ { ( 1 ) },$$

where α represents a state of two identical relativistic particles with momenta k 1 and k 2 , and dα is the appropriate measure for a trace over these states. The first-order forward scattering T matrix element is just the vertex factor m g 2 2 times the volume L .

Figure 3: An example of a 'disconnected' forward scattering amplitude that nevertheless contributes to the free energy due to identical particle exchange. The labels at the bottom refer to incoming momenta, and the top outgoing momenta.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-743aacb4-07cc-4b52-bd7f-714e414dae1d.jpg)

Diagrams involving a single higher-order vertex, such as the second-order diagram involving the ϕ 6 vertex in Fig 2, follow similarly. The counterterm diagrams will again end up canceling the I 0 terms in the tadpole integrals, and the net contribution to the free energy density is

$$f _ { \phi ^ { 2 n } } = \frac { 1 } { L } \, \frac { 1 } { n! } \int \left ( \prod _ { i = 1 } ^ { n } \frac { d k _ { i } } { 2 \pi } \frac { 1 } { 2 E _ { i } } \frac { e ^ { - \beta E _ { i } } } { 1 - e ^ { - \beta E _ { i } } } \right ) m ^ { 2 } g ^ { 2 n - 2 } L.$$

This has the appropriate measure for n relativistic particles, and involves the n -body forward scattering amplitude T αα = m g 2 2 n -2 L , but once again there are additional factors of (1 -e -βE i ) -1 . These factors have to do with bosonic statistics, and they arise by considering auxiliary particles in the states α which do not participate in the two-body scattering. Due to identical particle exchange there will be parts of the forward scattering T matrix which are still connected in the sense that they involve a single factor of L and can not be written as a product of two or more free energy terms. An example is given in Fig 3. The auxiliary particle line leads to a delta function 2 E 3 2 πδ k ( 3 -k 1 ) that cancels the trace over the third particle momentum and leads to a correction which agrees with the first non-trivial term in the expansion of (1 -e -βE 1 ) -1 .

These (1 -e -βE i ) -1 factors involve arbitrarily high particle number N as far as the fugacity expansion (2.8) is concerned, but each term in the expansion involves only some finite number of n particles interacting through non-trivial scattering. It is useful to consider an expansion in the number of non-trivially interacting particles, which will be referred to as the n-body expansion. An ' n -body correction' like (3.7) will generically involve n factors of bosonic occupation numbers n i .

For later reference we will write the three-body correction due to the ϕ 6 vertex explicitly in terms of occupation numbers as a function of rapidity

$$f _ { \phi ^ { 6 } } = \frac { 1 } { 3 } \frac { m ^ { 2 } g ^ { 4 } } { 1 6 } \left ( \int \frac { d \theta } { 2 \pi } n ( \theta ) \right ) ^ { 3 },$$

where n θ ( ) = ( e βm cosh θ -1 ) -1 .

## 3.2 Cutting a general diagram

The previous example involved using the Poisson resummation formula to write finite temperature propagators in terms of an infinite sum over an index l . The l = 0 term which is present at zero-temperature was split from the rest. Now we will generalize this procedure.

A general r -loop free energy diagram δf may be written in terms of r independent frequencies, each with its own summation over exponentials due to the Poisson resummation formula,

$$\delta f = \int \prod _ { j = 1 } ^ { r } \frac { d k _ { j } } { 2 \pi } \frac { d \omega _ { j } } { 2 \pi } \frac { \sum _ { l _ { j } } e ^ { - i l _ { j } \beta \omega _ { j } } } { \omega _ { j } ^ { 2 } + E _ { j } ^ { 2 } } \mathcal { M } \left ( \omega, k \right ).$$

The factor M collects all of the combinatorial and vertex factors, as well as any additional propagators which are fixed by the delta functions at the vertices to be functions of the r independent ω , k j j .

̸

The summation over l j may be broken up into the l j = 0 term and the l j = 0 terms,

$$\sum _ { l _ { j } } e ^ { - i l _ { j } \beta \omega _ { j } } = 1 + \sum _ { \sigma _ { j } = \pm 1 } \sum _ { l _ { j } > 0 } e ^ { - i \sigma _ { j } l _ { j } \beta \omega _ { j } }.$$

̸

After decomposing the r propagators into two terms in this manner, there are a total of 2 r terms in δf . We refer to the propagators that keep the l = 0 terms as cut , and the those that keep the l = 0 term as uncut . We may index the 2 r terms in δf by a set π which contains the indices of the cut propagators. The other uncut propagators are integrated over and absorbed into the M factor, which is now denoted M π .

$$\delta f = \sum _ { \pi } \int \prod _ { j \in \pi } \frac { d k _ { j } } { 2 \pi } \frac { d \omega _ { j } } { 2 \pi } \sum _ { \sigma _ { j } = \pm 1 } \sum _ { l _ { j } > 0 } \frac { e ^ { - i \sigma _ { j } l _ { j } \beta \omega _ { j } } } { \omega _ { j } ^ { 2 } + E _ { j } ^ { 2 } } \mathcal { M } _ { \pi } \left ( \omega, k \right ).$$

The quantity M π is easy to understand diagramatically. It is just given by cutting the cut propagators of the original free energy diagram and taking them to be external legs. The remaining uncut propagators of the diagram are taken to be at zero-temperature. M π is an off-shell forward scattering amplitude depending on the frequency and spatial momenta of the cut lines.

If we momentarily assume that M π ( ω ) is analytic then only the pole at ω j = -iσ E j j due to the cut propagator will contribute to the ω j integral. Then after integration over ω and summation of the geometric series in l we get the formula,

$$\delta f = \sum _ { \pi } \int \prod _ { j \in \pi } \frac { d k _ { j } } { 2 \pi } \frac { n ( E _ { j } ) } { 2 E _ { j } } \sum _ { \sigma } \mathcal { M } _ { \pi } \left ( - i \sigma E, k \right ).$$

This is now roughly in the form of the DMB formula (1.5). M π has been put on-shell, and the signs σ determine whether the momentum of a cut line is considered to be incoming or outgoing. If there are n cut lines belonging to π then this is an n -body contribution to the free energy diagram δf , and the factors of n E ( ) automatically contain the higher-order corrections due to particle exchange.

Generically M π ( ω ) will not be analytic. There will be an additional contributions to the frequency integral from poles or branch cuts in M π . This will be easy enough to deal with in Sec 4 for the case of diagrams with only on-shell contributions. The result will agree with the full DMB formula that includes the higher order in ˆ corrections which are neglected in (1.5). T

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-805e304d-4cac-41ae-9e4c-2188e9ee94d9.jpg)

2

Figure 4: The 'bubble' and 'melon' corrections to the free energy of a ϕ 4 theory at second order, on the left and right of the top row, respectively. The two-body and three-body forward scattering amplitudes corresponding to these diagrams are shown in the second and third rows, with numbers to indicate distinct momenta.

## 3.3 The bubble and melon diagrams

For now let us consider the free energy diagrams at second-order in the sinh-Gordon and ϕ 4 theories. As discussed in Sec 2.4, perturbation theory at this order will clarify how three-particle scattering amplitudes are treated in the DMB approach. The diagram involving a single sextic vertex in sinh-Gordon has already been considered in (3.8). The remaining two diagrams -the bubble and melon - are shown in Fig 4.

The previous section has made clear which forward scattering amplitudes are associated to a given thermal field theory diagram. Lines of the thermal field theory diagram are cut to create incoming and outgoing external lines with the same momenta. Conversely, we may 'glue together' incoming and outgoing lines with the same momentum in a forward scattering amplitude (the lines labeled by the same number in Fig 4) to recover a free energy diagram.

## 3.3.1 The bubble diagram

The bubble diagram is shown on the left side of Fig 4. As in the single vertex case, free energy diagrams involving counterterms will cancel the uncut I 0 part of any tadpoles, so the two tadpoles on the sides of the diagram are always cut.

The two-body part of the diagram f bubble 2 , involves an uncut central bubble. Since in this case

the amplitude M π does not depend on the cut momenta, it can be readily calculated

$$f _ { \text{bubble}, 2 } = - \frac { m ^ { 4 } g ^ { 4 } } { 1 6 } \left ( \int \frac { d \theta } { 2 \pi } n ( \theta ) \right ) ^ { 2 } \frac { 1 } { 4 \pi m ^ { 2 } }.$$

The scattering amplitude that would correspond to the three-body part of the bubble diagram f bubble 3 , is given on the left of the bottom row of Fig 4. This diagram has a divergence due to the internal propagator being fixed by momentum-conservation to have the on-shell momentum labeled 3. This is an example of an 'on-shell contribution.'

Going back to the general cutting argument, this divergence just indicates that M π also has a pole as a function of ω 3 , and this needs to be taken into account in the integration over ω 3 . So dealing with the on-shell contribution is no major difficulty, it just involves integrating over a double pole

$$f _ { \text{bubble}, 3 } & = - \frac { m ^ { 4 } g ^ { 4 } } { 1 6 } \left ( \int \frac { d \theta } { 2 \pi } n ( \theta ) \right ) ^ { 2 } \int \frac { d k _ { 3 } } { 2 \pi } \sum _ { \sigma _ { 3 }, l _ { 3 } } \int \frac { d \omega _ { 3 } } { 2 \pi } \frac { e ^ { - i \sigma _ { 3 } l _ { 3 } \beta \omega _ { 3 } } } { \left ( \omega _ { 3 } ^ { 2 } + E _ { 3 } ^ { 2 } \right ) ^ { 2 } } \\ & = - \frac { m ^ { 2 } g ^ { 4 } } { 1 6 } \left ( \int \frac { d \theta } { 2 \pi } n ( \theta ) \right ) ^ { 2 } \int \frac { d k _ { 3 } } { 2 \pi } \frac { m ^ { 2 } } { 2 E _ { 3 } ^ { 2 } } \left ( - \frac { d } { d E _ { 3 } } + \frac { 1 } { E _ { 3 } } \right ) n ( E _ { 3 } ). \\ \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{matrix}$$

This expression will be derived from the DMB formula in Sec 4. By using similar manipulations to those in Appendix B.1 it can be brought to a simplified form

$$f _ { \text{bubble}, 3 } = \frac { 1 } { 3! } \frac { m ^ { 2 } g ^ { 4 } } { 1 6 } \beta \frac { \partial } { \partial \beta } \left ( \int \frac { d \theta } { 2 \pi } n ( \theta ) \right ) ^ { 3 }.$$

## 3.3.2 The melon diagram: two-body part

The forward scattering diagrams corresponding to cutting the melon diagram are shown on the right of Fig 4. In the following sections we will use these diagrams in the DMB approach to evaluate the melon diagram and show it involves only contributions with collinear momenta. Of course the melon diagram may also be evaluated directly from thermal field theory, and this is shown using methods similar to [23, 24] in Appendix A.2.

Note that the two-body amplitudes can be written in terms of the one-loop integral

$$J ( p ) \equiv \int \frac { d ^ { 2 } q } { ( 2 \pi ) ^ { 2 } } \frac { 1 } { q ^ { 2 } + m ^ { 2 } } \frac { 1 } { ( p - q ) ^ { 2 } + m ^ { 2 } }.$$

The bubble diagram was associated to the u-channel of two-body scattering, with no net momentum flowing through the loop, and (3.9) indeed involves the factor J (0) = (4 πm 2 ) -1 . The melon diagram instead involves the t-channel and s-channel. By the DMB approach, the two-body part of the melon should be

$$f _ { \text{melon}, 2 } = - \frac { m ^ { 4 } g ^ { 4 } } { 2 } \frac { 1 } { 2! } \int \left ( \prod _ { j = 1, 2 } \frac { d k _ { j } } { 2 \pi } \frac { n _ { j } } { 2 E _ { j } } \right ) \left [ J ( - i ( E _ { 1 } - E _ { 2 } ) ) + \mathcal { P } J ( - i ( E _ { 1 } + E _ { 2 } ) ) \right ].$$

The argument of the s-channel term J ( -i ( E 1 + E 2 )) sits on a branch cut, and some sort of iϵ prescription is needed to define it. The principal value prescription arising from the second-order DMB formula (2.26) is the correct one.

This can be calculated further, using the expression

$$J ( p ) = \frac { 1 } { 2 \pi p ^ { 2 } \xi } \log \left ( \frac { \xi + 1 } { \xi - 1 } \right ), \quad \xi \equiv \sqrt { 1 + 4 \frac { m ^ { 2 } } { p ^ { 2 } } }.$$

The sum of terms in (3.13) is proportional to a delta function,

$$\text{the sum of terms in (3.13) is proportional to a delta function, } J \left ( - i ( E _ { 1 } - E _ { 2 } ) \right ) + \frac { 1 } { 2 } \left [ J ( - i ( E _ { 1 } + E _ { 2 } ) + \epsilon ) + J ( - i ( E _ { 1 } + E _ { 2 } ) - \epsilon ) \right ] \\ & = \frac { \log e ^ { \theta } } { 4 \pi m ^ { 2 } \sinh \theta } + \frac { 1 } { 2 } \left [ - \frac { \log \left ( - e ^ { \theta } + i \epsilon \right ) } { 4 \pi m ^ { 2 } \sqrt { \sinh ^ { 2 } \theta - i \epsilon } } - \frac { \log \left ( - e ^ { \theta } - i \epsilon \right ) } { 4 \pi m ^ { 2 } \sqrt { \sinh ^ { 2 } \theta + i \epsilon } } \right ] \\ & = \frac { 2 \pi } { 1 6 m ^ { 2 } } \left [ \frac { i } { \pi \sqrt { \sinh ^ { 2 } \theta + i \epsilon } } - \frac { i } { \pi \sqrt { \sinh ^ { 2 } \theta - i \epsilon } } \right ] \\ & = \frac { 2 \pi } { 1 6 m ^ { 2 } } \delta ( \theta ). \\ \text{o in total } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

So in total

$$f _ { \text{melon}, \, 2 } = - \frac { 1 } { 1 6 } \frac { m ^ { 2 } g ^ { 4 } } { 1 6 } \int \frac { d \theta } { 2 \pi } n ( \theta ) ^ { 2 }.$$

## 3.3.3 The melon diagram: three-body part

The three-body part of the melon diagram is similarly calculated from the three-body amplitudes in Fig 4,

$$\text{in $freq 4$,} \\ f _ { \text{melon}, 3 } = - m ^ { 4 } g ^ { 4 } \frac { 1 } { 3! } \int \left ( \prod _ { j = 1, 2, 3 } \frac { d k _ { j } } { 2 \pi } \frac { n _ { j } } { 2 E _ { j } } \right ) \left [ \mathcal { P } \frac { 3 } { - ( E _ { 1 } + E _ { 2 } - E _ { 3 } ) ^ { 2 } + ( k _ { 1 } + k _ { 2 } - k _ { 3 } ) ^ { 2 } + m ^ { 2 } } \\ + \frac { 1 } { - ( E _ { 1 } + E _ { 2 } + E _ { 3 } ) ^ { 2 } + ( k _ { 1 } + k _ { 2 } + k _ { 3 } ) ^ { 2 } + m ^ { 2 } } \right ]. \\ \Lambda \, \text{on one $the commons are distinct $the non-zero one not each $and $the$rincimal $value$}$$

As long as the momenta are distinct the propagators are not on-shell, and the principal value regularization may be dropped.

The first propagator in (3.15) may be written in terms of rapidity as,

$$\frac { 1 } { - \left ( E _ { 1 } + E _ { 2 } - E _ { 3 } \right ) ^ { 2 } + \left ( k _ { 1 } + k _ { 2 } - k _ { 3 } \right ) ^ { 2 } + m ^ { 2 } } = \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } - e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 2 } } - e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) }.$$

There are two additional diagrams where index 3 is swapped with index 1 or 2 (hence the factor of three in (3.15)). Adding these terms with permuted indices leads to

$$\frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } - e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 2 } } - e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) } + \left ( 3 \leftrightarrow 1 \right ) + ( 2 \leftrightarrow 3 ) = \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 2 } } + e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) }.$$

But this simply cancels with the second propagator in (3.15),

$$\frac { 1 } { - \left ( E _ { 1 } + E _ { 2 } + E _ { 3 } \right ) ^ { 2 } + \left ( k _ { 1 } + k _ { 2 } + k _ { 3 } \right ) ^ { 2 } + m ^ { 2 } } = - \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 2 } } + e ^ { \theta _ { 3 } } \right ) \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) }.$$

The cancelation of the three-body forward scattering diagrams associated to the melon diagram is not surprising given the cancelations which must take place in the integrable sinh-Gordon model. 4 But actually the three-body part of the melon diagram does not vanish. The principal value regulation may not be ignored when the momentum become collinear, and a calculation in Appendix A.3 shows that the integrand of (3.15) is proportional to a delta function, much as in the two-body case. The final result is

$$f _ { \text{melon}, 3 } = - \frac { 1 } { 2 4 } \frac { m ^ { 2 } g ^ { 4 } } { 1 6 } \int \frac { d \theta } { 2 \pi } n ( \theta ) ^ { 3 }.$$

In Sec 5 these same results will be derived from the TBA. It will be discussed why similar collinear contributions are expected at higher order in perturbation theory, and how they agree with the appearance of bound states upon switching the sign of the coupling in the theory.

## 4 On-shell contributions and large N theories

The bubble diagram is part of a broader class of diagrams where the integrals over distinct loop momenta factorize. These are the same diagrams which would be singled out at leading order in a large N expansion. To study the net effect of these diagrams we consider extending the field in the ϕ 4 theory (3.1) to an N component field ϕ with Lagrangian

$$\mathcal { L } = \frac { 1 } { 2 } \left ( \partial \phi \cdot \partial \phi + \left ( m ^ { 2 } + \delta m ^ { 2 } \right ) \phi ^ { 2 } \right ) + \frac { \lambda } { 4 N } \left ( \phi ^ { 2 } \right ) ^ { 2 }.$$

The counterterm δm 2 = -λI 0 is chosen to ensure that m 2 is the physical mass at large N and zero temperature. The quartic term may be simplified by introducing a Hubbard-Stratonovich field α

$$\mathcal { L } = \frac { 1 } { 2 } \left ( \partial \phi \cdot \partial \phi + m ^ { 2 } \phi ^ { 2 } \right ) - \frac { N } { \lambda } \alpha ^ { 2 } + \alpha \left ( \phi ^ { 2 } - N I _ { 0 } \right ).$$

Here α has a quadratic term that vanishes in the strong coupling limit λ → ∞ . In this limit α can be understood as a Lagrange multiplier field, and the model reduces to the integrable O N ( ) non-linear sigma model.

Note that the counterterm was set to ensure that α has no expectation value at leading order for zero-temperature, but at finite temperature it does pick up an expectation value. We may absorb the expectation value into a temperature-dependent mass M β ( ), and a redefined field α β ≡ α -1 2 ( M 2 -m 2 ),

$$\mathcal { L } = \frac { 1 } { 2 } \left ( \partial \phi \cdot \partial \phi + M ( \beta ) ^ { 2 } \phi ^ { 2 } \right ) + \alpha _ { \beta } \left ( \phi ^ { 2 } - N I _ { 0 } \right ) - \frac { 1 } { 2 } \left ( M ( \beta ) ^ { 2 } - m ^ { 2 } \right ) N I _ { 0 }.$$

4 It can also be shown that the three-body diagrams associated with the bubble, displaced slightly from the forward direction, are canceled entirely by the three-body diagram given by the ϕ 6 vertex.

The leading-order free energy of this form of the O N ( ) non-linear sigma model was considered by Balog and Hegedus in [29]. They used an exact expression for the free energy density which follows directly from (4.2),

$$f = \frac { N m ^ { 2 } } { 2 \pi } \left [ - \frac { M ^ { 2 } } { m ^ { 2 } } \sum _ { j = 1 } ^ { \infty } K _ { 2 } \left ( j m \beta \frac { M } { m } \right ) + \frac { 1 } { 4 } \left ( \frac { M ^ { 2 } } { m ^ { 2 } } - 1 \right ) \right ].$$

This may be expanded in orders of fugacity by solving the gap equation [30] for M/m order-byorder,

$$\log \frac { M } { m } = 2 \sum _ { j = 1 } K _ { 0 } \left ( j \beta m \frac { M } { m } \right ).$$

Using this method, Balog and Hegedus were able to find the second-order in fugacity contribution f (2) and show it agrees with the DMB formula for two-particle scattering. In principle there is no obstacle to using this method to calculate arbitrary orders f ( j ) , but matching to the DMB even for f (3) encounters the problem of the forward scattering singularities. In the following sections we will give an alternate method to calculate the free energy density, and show how it may be matched to forward scattering amplitudes at arbitrary order.

However, note in passing that there is a direct connection between the thermal mass M β ( ) and forward scattering amplitudes along the same lines as the DMB relation. In [7, 8], L¨scher derives u a relation which in our present circumstances takes the form

$$\frac { M - m } { m } = - \frac { 1 } { 2 m ^ { 2 } } \int \frac { d p } { 2 \pi } \frac { 1 } { 2 E } e ^ { - \beta E } F + \mathcal { O } \left ( e ^ { - 2 m \beta } \right ).$$

Here F is the sum of the two-body forward scattering amplitudes with the other particles of the theory, which works out to be -8 πm 2 in the case of the large N limit of the O N ( ) non-linear sigma model. Expanding the gap equation (4.4) to first order in fugacity we recover this result. Higher order in fugacity corrections to M 2 can be derived from thermal self-energy diagrams calculated using the Feynman rules of the following section, and these correspond to multi-particle forward scattering amplitudes in the same way.

## 4.1 Diagrams in the large N theory

The equations (4.3) and (4.4) give a compact expression for the free energy, but actually solving them order-by-order quickly leads to a tedious calculation. We will organize the calculation of the free energy by an approach involving thermal field theory diagrams, which then can be quickly related to forward scattering amplitudes through the cutting rules discussed in Sec 3.2.

Let us return to the form of the Lagrangian involving the zero-temperature mass m 2 (4.1). We may keep finite λ dependence. There is a free boson contribution f 1 coming from the quadratic action of the ϕ field

$$f _ { 1 } = \frac { N } { 2 } \int _ { \beta } \frac { d ^ { 2 } p } { ( 2 \pi ) ^ { 2 } } \log \left ( p ^ { 2 } + m ^ { 2 } \right ) = - \frac { N m ^ { 2 } } { \pi } \sum _ { j = 1 } ^ { \infty } \frac { K _ { 1 } ( j \beta m ) } { j \beta m }.$$

There are further corrections coming from thermal field theory diagrams involving internal α propagators. Since each α propagator introduces a factor of 1 /N , at leading order it must be paired with a ϕ loop, which introduces a compensating factor of N . This implies that at leading order there will be no momentum flowing through the α propagators.

So to calculate the leading order diagrams, we must consider a class of integrals which may be derived from the basic tadpole integral I in (3.5),

$$I _ {, n } \equiv \int _ { \beta } \frac { d ^ { 2 } p } { ( 2 \pi ) ^ { 2 } } \frac { 1 } { ( p ^ { 2 } + m ^ { 2 } ) ^ { n + 1 } } = \frac { ( - 1 ) ^ { n } } { n! } \frac { d ^ { n } I } { d ( m ^ { 2 } ) ^ { n } }.$$

These integrals may be calculated by taking derivatives of the formula,

$$I = \frac { 1 } { 4 \pi } \log \frac { M _ { U V } ^ { 2 } } { m ^ { 2 } } + \frac { 1 } { \pi } \sum _ { j = 1 } ^ { \infty } K _ { 0 } ( j \beta m ),$$

where M UV is a UV cutoff. Using the Bessel function identity K α +1 ( x ) = ( α x -d dx ) K α ( x ), we can show

$$I _ {, n } = \frac { 1 } { 4 \pi n m ^ { 2 n } } \left [ 1 + \frac { x _ { n } } { 2 ^ { n - 2 } ( n - 1 )! } \right ],$$

$$x _ { n } \equiv \sum _ { j = 1 } ^ { \infty } ( j \beta m ) ^ { n } \, K _ { n } ( j \beta m ).$$

The integrals I ,n neatly decompose into an 'uncut' piece I 0 ,n = (4 πnm 2 n ) -1 , and a 'cut' piece involving the quantity x n . The uncut bubbles I 0 1 , will be used to form effective propagators for α .

An effective action for the zero mode of the α field may be found by integrating out the ϕ field in (4.1),

$$\mathcal { L } _ { \text{eff} } = - N \left ( \lambda ^ { - 1 } + I _ {, 1 } \right ) \alpha ^ { 2 } + N \sum _ { j = 3 } ^ { \infty } \frac { ( - 2 ) ^ { j - 1 } } { j } I _ {, j - 1 } \alpha ^ { j }.$$

The effective propagator may be written in terms of a geometric series of cut bubbles,

$$- \frac { \lambda } { 2 N } \frac { 1 } { 1 + \lambda I _ {, 1 } } = - \frac { 2 \pi m ^ { 2 } \Lambda } { N } \sum _ { j = 0 } ^ { \infty } \left ( - \frac { 2 \pi m ^ { 2 } \Lambda } { N } \frac { N x _ { 1 } } { \pi m ^ { 2 } } \right ) ^ { j }, \quad \Lambda \equiv \frac { \lambda } { \lambda + 4 \pi m ^ { 2 } }.$$

The renormalized coupling Λ is simply equal to 1 in the non-linear sigma model limit.

Now we may present the Feynman rules for large N free energy diagrams:

- · α propagators are drawn as dashed lines, each contributing -2 πm 2 Λ /N .
- · All diagrams are tree diagrams in terms of the α propagators.
- · ϕ loops are drawn as solid lines, and are always taken to be 'cut.'

Figure 5: Large N free energy diagrams up to four ϕ loops, corresponding to two-body (A), threebody (B), and four-body (C) scattering amplitudes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-cf5927e8-b3fd-47b9-8720-827ea25f12e4.jpg)

- · An uncut α j vertex is drawn as a solid point, and contributes a factor N ( -2) j -1 I 0 ,j -1 /j .
- · A cut α j vertex is drawn as a ϕ loop, and contributes an additional factor of 2 3 -j x j -1 / j ( -2)!.
- · A ϕ tadpole contributes Nx /π 0 .
- · ϕ bubbles may be inserted in an α propagator, contributing factors of Nx / πm 1 ( 2 ).

The free energy diagrams up to four ϕ loops are shown in Fig 5. Together with the 'onebody' contribution (4.6), they may be expanded up to fourth-order in fugacity to match with the calculation of f (4) via the gap equation (4.4).

## 4.2 On-shell contributions from the DMB formula

In general, a free energy diagram with n ϕ loops corresponds to the n -body forward scattering amplitude where each loop is cut. Diagrams that only involve ϕ tadpoles are straightforward to match to the DMB formula, as in Sec 3.1. However, a bubble insertion in an α propagator will correspond to a single internal propagator that goes on-shell, much like the three-body part of the bubble diagram in Sec 3.3.1. And a ϕ loop vertex with p attached α propagators will involve p -1 internal propagators going on-shell.

These are all what we have been referring to as on-shell contributions, and they are not specific to the large N theory. From the previous section it is clear that a forward scattering diagram with j -1 internal propagators going on-shell should involve a factor of x j -1 , which is a sum of Bessel functions K j -1 . So the problem of dealing with on-shell contributions in the DMB approach amounts to deriving this sum of Bessel functions from the DMB formula.

To do this it is helpful to use another expression for x s , which is derived in Appendix B.1,

$$x _ { s } = \int \frac { d k } { 2 E } \frac { m ^ { 2 s } } { E ^ { s } } \sum _ { l = 0 } ^ { s } \frac { ( s + l )! } { ( 2 E ) ^ { l } ( s - l )!! } \left ( - \frac { d } { d E } \right ) ^ { s - l } n ( E ).$$

Figure 6: Two distinct intermediate states contribute to the three-body part of the bubble diagram. The particles that cross the dotted line are acted upon by the propagator P G 0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-1152091f-695e-4fb4-8554-a791091d1e63.jpg)

In particular, the simplest case is

$$x _ { 1 } = \int d k \frac { m ^ { 2 } } { 2 E ^ { 2 } } \left ( - \frac { d } { d E } + \frac { 1 } { E } \right ) n ( E ),$$

which is clearly involved in the three-body part of the bubble diagram (3.10). We will first sketch how the bubble diagram is derived from the DMB formula, and then move on to the general case.

For simplicity we will ignore particle exchange corrections in the discussion that follows. This amounts to truncating the summation in x s to the first term and replacing n E ( ) by e -βE . For example (4.14) becomes

$$\beta m \, K _ { 1 } ( \beta m ) = \int d k \frac { m ^ { 2 } } { 2 E ^ { 2 } } \left ( \beta + \frac { 1 } { E } \right ) e ^ { - \beta E }.$$

## 4.2.1 DMB calculation of the bubble diagram

Given the second-order DMB formula (2.26), the bubble diagram will involve the amplitude ⟨ α | ( G 0 -¯ G 0 ) V P G V α 0 | ⟩ . The three-particle state is | α ⟩ = a † ( p 1 ) a † ( p 2 ) a † ( p 3 ) 0 | ⟩ and the interaction V = m g 2 2 4! ∫ d xϕ x 1 ( ) 4 . This involves the free field operator

$$\phi ( x ) = \int \frac { d k } { 2 \pi } \frac { 1 } { \sqrt { 2 E _ { k } } } \left ( a ^ { \dagger } ( k ) e ^ { i k x } + a ( k ) e ^ { - i k x } \right ).$$

This is essentially an old-fashioned perturbation theory calculation, and the rightmost copy of V can act on the state α in two different ways leading to an intermediate state with either three or five particles, see Fig 6. The diagram with a three-particle intermediate state is the source of the forward scattering divergence, and it involves a propagator with the same energy difference E -E α that appears in the would-be delta function G 0 -¯ G 0 .

As noted already in [1], this is equivalent to a derivative of the energy delta function

$$\left ( G _ { 0 } - \bar { G } _ { 0 } \right ) \mathcal { P } G _ { 0 } = \frac { - 2 i \epsilon } { ( E - E _ { \alpha } ) ^ { 2 } + \epsilon ^ { 2 } } \frac { E - E _ { \alpha } } { ( E - E _ { \alpha } ) ^ { 2 } + \epsilon ^ { 2 } } = - \frac { 1 } { 2 } \frac { \partial } { \partial E } \left ( G _ { 0 } - \bar { G } _ { 0 } \right ).$$

The E derivative will produce a factor of β by acting on e -βE after integration by parts, so effectively the propagator is replaced by

$$\mathcal { P } G _ { 0 } \rightarrow - \frac { 1 } { 2 } \beta.$$

Figure 7: The free energy diagram f α p (left). This corresponds to a forward scattering amplitude (right) involving p +1 particles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-ee2879e8-e5ef-4468-9dc5-a22d12542142.jpg)

On the other hand, in the diagram with a five-particle intermediate state, all of the combinatorics end up being the same but the propagator never goes on-shell so we can safely set E → E 1 + E 2 + E , 3

$$\mathcal { P } G _ { 0 } \rightarrow \mathcal { P } \frac { 1 } { E - ( E _ { 1 } + E _ { 2 } + 3 E _ { 3 } ) } = - \frac { 1 } { 2 } \frac { 1 } { E _ { 3 } }.$$

This is how the factor β +1 /E 3 in (4.15) arises from the DMB perspective.

## 4.2.2 Higher-order on-shell contributions

The basic large N free energy diagram involving a single cut α p vertex and p tadpoles is shown in Fig 7. It corresponds to a forward scattering amplitude among p +1 particles, where p -1 internal propagators all go on-shell with momentum k ′ . Following the Feynman rules in Sec 4.1, the free energy f α p associated to this diagram is

$$f _ { \alpha ^ { p } } & = - \frac { N m ^ { 2 } } { \pi } \frac { ( - 2 \Lambda ) ^ { p } } { p! } x _ { 0 } ^ { p } x _ { p - 1 }.$$

This may be calculated in old-fashioned perturbation theory as in the previous subsection. The α propagators lead to an effective interaction V = Λ πm 2 N ∫ d x ϕ x 1 ( ( ) 2 ) 2 which may be used to calculate ˆ in the DMB formula. For any given diagram with T p interaction vertices, the field operators may be contracted in p ! different orderings, which involve different energies of intermediate states.

The most singular contribution is due to the monotonic ordering such that each intermediate state has the same energy as the external particles E α . Each propagator G 0 evaluates to ( E -E α + iϵ ) -1 . This means that for the purposes of calculating this contribution we can treat G 0 like a c-number in the DMB formula,

$$\log \hat { S } = \log \left ( 1 + \frac { \left ( G _ { 0 } - \bar { G } _ { 0 } \right ) V } { 1 - G _ { 0 } V } \right ) = \log \left ( 1 - \bar { G } _ { 0 } V \right ) - \log \left ( 1 - G _ { 0 } V \right ).$$

Now the logarithms may be expanded up to p th order in V ,

$$\left ( \log \hat { S } \right ) ^ { ( p ) } = V ^ { p } \frac { 1 } { p } \left ( G _ { 0 } ^ { p } - \hat { G } _ { 0 } ^ { p } \right ) = V ^ { p } \frac { ( - 1 ) ^ { p - 1 } } { p! } \frac { d ^ { p - 1 } } { d E ^ { p - 1 } } \left ( G _ { 0 } - \bar { G } _ { 0 } \right ).$$

So the DMB formula for the most singular ordering leads to a ( p -1)th derivative of an energy delta function. This will correspond to the l = 0 term in the expansion (4.13) for x p -1 . The other orderings which are less singular will lead to lower-order derivatives of a delta function, and these will correspond to the l &gt; 0 terms in the expansion. This is shown in detail for p = 4 in Appendix B.2.

In the remainder of this section we will continue to focus on the most singular ordering, and show agreement with the expression (4.17) for the free energy f α p . After integrating by parts, the ( p -1)th derivative leads to a factor of β p -1 , and the DMB formula reduces to

$$f _ { \alpha ^ { p } } = \frac { 1 } { L } \frac { N ^ { p + 1 } } { ( p + 1 )! } \int \left ( \prod _ { \alpha } ^ { p + 1 } \frac { d k _ { \alpha } } { ( 2 \pi ) 2 E _ { \alpha } } \right ) e ^ { - \beta \sum _ { \alpha } E _ { \alpha } } \left \langle V ^ { p } \right \rangle _ { \alpha ^ { p } } \frac { ( - 1 ) ^ { p - 1 } } { p! } \beta ^ { p - 1 }.$$

⟨ V p ⟩ α p indicates the factors arising from the contraction of V p as in the forward scattering diagram in Fig 7. This factor is the same for all orderings,

$$\langle V ^ { p } \rangle _ { \alpha ^ { p } } = \left ( 8 \times \frac { \Lambda \pi m ^ { 2 } } { N } \right ) ^ { p } \frac { ( p + 1 )! } { ( 2 E ^ { \prime } ) ^ { p - 1 } } L.$$

There are p vertex factors from the effective interaction V , each introducing a combinatorial factor of 8. There are ( p +1)! permutations of the large N indices of the particles, a factor of L from the momentum conservation delta function, and factors of (2 E ′ ) -1 from contracting the fields in the p -1 internal propagators.

In total

$$f _ { \alpha ^ { p } } = - \frac { N m ^ { 2 } } { \pi } \frac { ( - 2 \Lambda ) ^ { p } } { p! } \left ( \int \frac { d k } { 2 E } e ^ { - \beta E } \right ) ^ { p } \int \frac { d k ^ { \prime } } { 2 E ^ { \prime } } \frac { m ^ { 2 ( p - 1 ) } } { ( E ^ { \prime } ) ^ { p - 1 } } \beta ^ { p - 1 } e ^ { - \beta E ^ { \prime } }$$

which indeed agrees with (4.17) and the l = 0 term of (4.13) at lowest order in fugacity.

## 5 The TBA and bound states

Although our motivation is coming from more general non-integrable theories like the ϕ 4 theory, it is useful to consider a perturbative expansion of the TBA in order to have an extra check on the validity of the DMB approach. In Sec 5.2 we will consider the TBA for the attractive Lieb-Liniger model and show that the perturbative expansion is consistent for both signs of the coupling. The collinear contributions to the free energy cancel delicately with the contributions coming from the bound states that appear for attractive coupling. In Sec 5.3 the DMB approach is used to calculate the density of states in the two-particle sector of the attractive LL model, and it is shown how the two-body bound state arises.

## 5.1 Perturbative expansion of the TBA

We have already compared the DMB formula and the TBA at second order in fugacity back in Sec 2.3. Now to consider the details of three-body forward scattering contributions it is useful to instead expand the TBA to second order in the coupling but to all orders in fugacity.

The two-body and three-body parts of the free energy density can be distinguished by the number of appearances of the kernel. The two-body part (2.25) has already been suggested from the fugacity expansion,

$$f _ { 2 } ^ { ( 2 ) } = + T \int \frac { d p } { 2 \pi } \frac { 1 } { e ^ { \beta E } - 1 } \left [ \log \left ( 1 - e ^ { - \beta E } \right ) * K ^ { ( 2 ) } \right ].$$

The second-order three-body part is

$$f _ { 3 } ^ { ( 2 ) } & = - T \int \frac { d p } { 2 \pi } \frac { 1 } { e ^ { \beta E } - 1 } \left ( \frac { 1 } { 2 } \frac { ( \log \left ( 1 - e ^ { - \beta E } \right ) * K ^ { ( 1 ) } \right ) ^ { 2 } } { 1 - e ^ { - \beta E } } - \left [ \frac { 1 } { e ^ { \beta E } - 1 } \left ( \log \left ( 1 - e ^ { - \beta E } \right ) * K ^ { ( 1 ) } \right ) \right ] * K ^ { ( 1 ) } \right ) \\ & = + \frac { T } { 2 } \frac { \partial } { \partial \beta } \int \frac { d \theta } { 2 \pi } \frac { 1 } { e ^ { \beta E } - 1 } \left ( \log \left ( 1 - e ^ { - \beta E } \right ) * K ^ { ( 1 ) } \right ) ^ { 2 }. \\ \text{----} \quad & \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \, - \,.$$

The second, simpler form of f 3 is valid for SG but not LL.

## 5.1.1 Perturbative expansion of the kernel

Some care is required in defining the perturbative expansion of the kernel. Recall that the kernel for both SG and LL is proportional to

$$K _ { L L } = \frac { 1 } { \pi } \frac { c } { p ^ { 2 } + c ^ { 2 } } \simeq \frac { c } { \pi p ^ { 2 } } \sum _ { j = 0 } ^ { \infty } \left ( - \frac { c ^ { 2 } } { p ^ { 2 } } \right ) ^ { j }.$$

̸

This perturbative expansion is valid for p = 0, but it fails at p = 0 and so can not be used for calculating the collinear contributions, such as the melon diagram. The divergence of the K (1) term at p = 0 must be regularized, and there will also be additional terms involving derivatives of delta functions for even powers of c .

A valid perturbative expansion of the kernel that takes the behavior near p = 0 into account may be found by simply Fourier transforming before expanding. This is easy to see for the LL kernel,

$$K _ { L L } ( p ) = \int \frac { d u } { 2 \pi } e ^ { - c | u | } e ^ { - i p u } = \sum _ { j = 0 } ^ { \infty } \frac { ( - 1 ) ^ { j } \, c ^ { j } } { j! } \int \frac { d u } { 2 \pi } | u | ^ { j } e ^ { - i p u }.$$

Even powers of j are just even derivatives of delta functions. Odd powers of j involve odd derivatives of the inverse Fourier transform of sgn u , which may be understood as the principal value of 1 /p ,

$$\mathcal { P } _ { \frac { 1 } { p } } = \frac { p } { p ^ { 2 } + \epsilon ^ { 2 } } = \pi i \int \frac { d u } { 2 \pi } \text{sgn} \, u \, e ^ { - \epsilon | u | } e ^ { - i p u }.$$

So this gives a definite regularization at p = 0 for the inverse powers of p 2 in the naive expansion of the kernel (5.3).

The situation is qualitatively the same for the SG kernel.

$$K _ { S G } ( \theta ) = \frac { 1 } { 2 \pi } \frac { 2 ( c / m ) \cosh \theta } { \sinh ^ { 2 } \theta + ( c / m ) ^ { 2 } } = \int \frac { d \nu } { 2 \pi } \frac { e ^ { - \gamma | \nu | } + e ^ { - ( \pi - \gamma ) | \nu | } } { 1 + e ^ { - \pi | \nu | } } e ^ { - i \nu \theta },$$

where γ ≡ arcsin( c/m ). This may be expanded to O ( c 2 ) ,

$$\frac { e ^ { - \gamma | \nu | } + e ^ { - ( \pi - \gamma ) | \nu | } } { 1 + e ^ { - \pi | \nu | } } = 1 - c \tanh \frac { \pi \nu } { 2 } \left ( \frac { \nu } { m } \right ) + \frac { c ^ { 2 } } { 2 } \left ( \frac { \nu } { m } \right ) ^ { 2 } + \mathcal { O } ( c ^ { 3 } ) \,.$$

Recall that c = mg 2 8 -mg 4 64 π + O ( g 6 ) , so the total second-order SG kernel is,

$$K _ { S G } ^ { ( 2 ) } = - \frac { g ^ { 2 } } { 8 \pi } K _ { S G } ^ { ( 1 ) } - \frac { g ^ { 4 } } { 1 2 8 } \delta ^ { \prime \prime } ( \theta ).$$

When this is substituted in (5.1), the first term in K (2) SG leads to the two-body part of the bubble diagram (3.9) and the second term leads to the two-body part of the melon diagram (3.14).

## 5.1.2 The bubble and melon from the TBA

The three-body contribution f (2) 3 involves the first-order kernel

$$K _ { S G } ^ { ( 1 ) } = - \frac { g ^ { 2 } } { 8 \pi } \frac { d } { d \theta } \mathcal { P } \frac { 1 } { \sinh \theta }.$$

The principal value regularization appears since the function πi tanh( πν/ 2) in (5.5) may be understood as the Fourier transform of P (sinh θ ) -1 5 .

First let us naively calculate f (2) 3 neglecting the details of the principal value regularization. We calculate (5.2), using the notation p 12 ≡ m sinh( θ 1 -θ 2 ),

$$f _ { 3 } ^ { ( 2 ) } & \simeq \frac { T } { 2 } \left ( \frac { m g ^ { 2 } } { 4 } \right ) ^ { 2 } \frac { \partial } { \partial \beta } \int \frac { d \theta _ { 1 } d \theta _ { 2 } d \theta _ { 3 } } { ( 2 \pi ) ^ { 3 } } \left ( e ^ { \beta E _ { 1 } } - 1 \right ) ^ { - 1 } \frac { 1 } { p _ { 1 2 } p _ { 1 3 } } \frac { 1 } { \partial \theta _ { 2 } \partial \theta _ { 3 } } \log \left ( 1 - e ^ { - \beta E _ { 2 } } \right ) \log \left ( 1 - e ^ { - \beta E _ { 3 } } \right ) \\ & \quad = \frac { T } { 2 } \left ( \frac { m g ^ { 2 } } { 4 } \right ) ^ { 2 } \frac { \partial } { \partial \beta } \int \frac { d \theta _ { 1 } d \theta _ { 2 } d \theta _ { 3 } } { ( 2 \pi ) ^ { 3 } } n _ { 1 } n _ { 2 } n _ { 3 } \frac { p _ { 2 } p _ { 3 } } { p _ { 1 2 } p _ { 1 3 } } \beta ^ { 2 }.$$

Now if the momentum factor is symmetrized it can be shown to reduce to 1,

$$\frac { p _ { 2 } p _ { 3 } } { p _ { 1 2 } p _ { 1 3 } } + \frac { p _ { 3 } p _ { 1 } } { p _ { 2 3 } p _ { 2 1 } } + \frac { p _ { 1 } p _ { 2 } } { p _ { 3 1 } p _ { 3 2 } } = 1,$$

$$f _ { 3 } ^ { ( 2 ) } \simeq \frac { T } { 3! } \left ( \frac { m g ^ { 2 } } { 4 } \right ) ^ { 2 } \frac { \partial } { \partial \beta } \int \frac { d \theta _ { 1 } d \theta _ { 2 } d \theta _ { 3 } } { ( 2 \pi ) ^ { 3 } } n _ { 1 } n _ { 2 } n _ { 3 } \beta ^ { 2 }.$$

The β derivative acting on β 2 produces the free energy due to the sextic vertex diagram (3.8), and the β derivative acting on one of the occupation numbers produces the free energy due to the bubble diagram (3.11).

The melon diagram is due to the details of the principal value prescription. Using identities like P 1 p 12 = 1 p 12 + iϵ + iπδ ( p 12 ), all momentum denominators may be adjusted so that they have the

5 The exact Fourier transform for finite ϵ is sinh θ sinh 2 θ + ϵ 2 = πi ∫ e -δν -e -( π -δ ν ) ( 1+ e -πν ) cosh δ e -iνθ with δ ≡ arcsin ϵ .

same sign of iϵ . The regularized version of (5.8) is

$$\begin{array} { c } \L \, & \, \cos \, & \, \sin \,, \\ p _ { 2 } p _ { 3 } \mathcal { P } \frac { 1 } { p _ { 1 2 } } \mathcal { P } \frac { 1 } { p _ { 1 3 } } + p _ { 3 } p _ { 1 } \mathcal { P } \frac { 1 } { p _ { 2 3 } } \mathcal { P } \frac { 1 } { p _ { 2 1 } } + p _ { 1 } p _ { 2 } \mathcal { P } \frac { 1 } { p _ { 3 1 } } \mathcal { P } \frac { 1 } { p _ { 3 2 } } \\ & = \frac { p _ { 1 2 } p _ { 1 3 } p _ { 2 3 } + \mathcal { O } ( \epsilon ) } { ( p _ { 1 2 } + i \epsilon ) ( p _ { 1 3 } + i \epsilon ) ( p _ { 2 3 } + i \epsilon ) } + 2 \pi i p _ { 1 } p _ { 2 } \delta ( p _ { 1 3 } ) \mathcal { P } \frac { 1 } { p _ { 1 2 } } + \frac { p _ { 1 } ^ { 2 } } { 4 } ( 2 \pi ) ^ { 2 } \delta ( p _ { 1 2 } ) \delta ( p _ { 1 3 } ). \\ \dots & \, \epsilon \, \cdots \, & \, \cdots \, & \, \epsilon \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \cdots \, & \, \end{array}$$

The first term reduces to 1, and leads to (5.9) as before. The second term is odd under the transformation p , p 1 2 → -p , 1 -p 2 and thus vanishes upon integration. It is the third term that leads to the three-body part of the melon diagram (3.19),

$$f _ { \text{melon}, \, 3 } = \frac { T } { 3! } \left ( \frac { m g ^ { 2 } } { 4 } \right ) ^ { 2 } \frac { \partial } { \partial \beta } \int \frac { d \theta } { 2 \pi } n ( \theta ) ^ { 3 } \frac { p ^ { 2 } } { 4 } \beta ^ { 2 } = - \frac { 1 } { 2 4 } \frac { m ^ { 2 } g ^ { 4 } } { 1 6 } \int \frac { d \theta } { 2 \pi } n ( \theta ) ^ { 3 }.$$

## 5.2 The attractive Lieb-Liniger model

The perturbative expansion of the TBA raises a minor puzzle. Upon flipping the sign of the coupling c , the kernel K LL flips sign, and thus the two-body correction f 2 in (2.25) flips sign. But in the thermal field theory approach, the even orders of perturbation theory do not depend on the sign of c , so how is this consistent?

The resolution of course is that the theory with negative c has bound states, and these must explicitly be taken into account in the formulation of the TBA. The contribution from the bound states should compensate the change in sign of the collinear contributions, so that the odd orders of perturbation theory flip sign and the even orders are invariant.

So far we have been focusing mostly on the sinh-Gordon model, and a change in the sign of c corresponds to considering the sine-Gordon model instead. The TBA for the sine-Gordon model [31] is complicated by the presence of the soliton sector. In the following sections we will focus on the simpler LL model instead (see Appendix C for the LL analogue of calculations in Sec 3.3)

The LL model with attractive coupling is not stable since there are bound states with arbitrarily large negative energy. However we may still make sense of the model by truncating at some finite order in the fugacity expansion. The TBA may be constructed using standard methods (see e.g. [32]). There is a distinct pseudoenergy ϵ j for each j -body bound state,

$$f = - T \sum _ { j } j \int \frac { d p } { 2 \pi } \log \left ( 1 + e ^ { - \beta \epsilon _ { j } } \right )$$

$$\epsilon _ { j } = E _ { j } - T \sum _ { k } \log \left ( 1 + e ^ { - \beta \epsilon _ { k } } \right ) \stackrel { \nu } { \ast } K _ { k j }, \quad K _ { k j } \equiv \frac { 1 } { 2 \pi i } \, \frac { d } { d u } \log S _ { k j } ( u ) \Big | _ { u = p _ { k } - p _ { j } } \,.$$

Here E j is the energy of the j -body bound state, and S kj is the two-body S-matrix between a k -body bound state and a j -body bound state. Only the lowest order quantities will be needed for our purposes

$$E _ { 1 } = \frac { p ^ { 2 } } { 2 m }, \quad E _ { 2 } = \frac { 2 p ^ { 2 } } { 2 m } - \frac { c ^ { 2 } } { 4 m }, \quad E _ { 3 } = \frac { 3 p ^ { 2 } } { 2 m } - \frac { c ^ { 2 } } { m },$$

$$S _ { 1 1 } ( u ) = \frac { u + i | c | } { u - i | c | }, \ \ S _ { 2 1 } ( u ) = S _ { 1 1 } ( u - i c / 2 ) S _ { 1 1 } ( u + i c / 2 ).$$

By expanding (5.12), the free energy at second-order in fugacity is

$$f _ { ( 2 ) } = - T \int \frac { d p } { 2 \pi } \left [ - \frac { 1 } { 2 } e ^ { - 2 \beta E _ { 1 } } + 2 e ^ { - \beta E _ { 2 } } + e ^ { - \beta E _ { 1 } } \left ( e ^ { - \beta E _ { 1 } } * K _ { 1 1 } \right ) \right ].$$

In the following section this will be derived in the DMB approach, but first let us consider the perturbative expansion in c .

At zeroth order in c , this ought to be given by the free boson gas (2.22),

$$f _ { ( 2 ) } ^ { ( 0 ) } = - T \int \frac { d p } { 2 \pi } \frac { 1 } { 2 } e ^ { - 2 \beta E _ { 1 } }.$$

Indeed E (0) 2 = 2 E 1 and K (0) 11 ( p ) = -δ p ( ), so this works out. Note that even though the theory is not interacting at zeroth order, the bound state contribution is needed to get the appropriate statistics of the free boson gas. The contributions to the free energy (5.11) due to ϵ 1 and ϵ 2 alone do not have the appropriate statistics.

At second order in c , the free energy is given by the non-relativistic analogue (C.4) of the two-body part of the melon diagram (3.14)

$$f _ { ( 2 ) } ^ { ( 2 ) } = - T \int \frac { d p } { 2 \pi } \frac { \beta c ^ { 2 } } { 4 m } e ^ { - 2 \beta E _ { 1 } }.$$

For positive c this is entirely due to the contribution of the kernel K (2) , but for negative c this flips sign

$$- T \int \frac { d p } { 2 \pi } e ^ { - \beta E _ { 1 } } \left ( e ^ { - \beta E _ { 1 } } * K _ { 1 1 } ^ { ( 2 ) } \right ) = + T \int \frac { d p } { 2 \pi } \frac { \beta c ^ { 2 } } { 4 m } e ^ { - 2 \beta E _ { 1 } }.$$

In order to get the appropriate free energy it must cancel with the binding energy of the two-body bound state ( e -βE 2 ) (2) = βc 2 4 m e -2 βE 1 . In this sense, the collinear contributions to the free energy have an intimate relation to the contributions due to bound states for negative c .

At higher orders in fugacity similar cancellations occur with the binding energy for higher order bound states. Introducing notation z j ≡ e -βE j , the next-order contribution is

$$\Phi _ { ( 3 ) } = - T \int \frac { d p } { 2 \pi } \left [ z _ { 1 } \left ( \frac { 1 } { 3 } z _ { 1 } ^ { 2 } + \left ( z _ { 1 } \left ( z _ { 1 } * K _ { 1 1 } \right ) \right ) * K _ { 1 1 } + \frac { 1 } { 2 } \left ( z _ { 1 } * K _ { 1 1 } \right ) ^ { 2 } \right ) \\ + \, z _ { 1 } \left ( - z _ { 1 } \left ( z _ { 1 } * K _ { 1 1 } \right ) - \frac { 1 } { 2 } \left ( z _ { 1 } ^ { 2 } * K _ { 1 1 } \right ) \right ) \\ + 3 z _ { 3 } + 2 z _ { 2 } \left ( z _ { 1 } * K _ { 2 1 } \right ) + z _ { 1 } \left ( z _ { 2 } * K _ { 2 1 } \right ) \right ]. \\ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ { \ } ^ { \ }. \ \Phi _ }.$$

The first line does not depend on the sign of c , and the second line flips sign. The third line involves the two-body and three-body bound states, and it is only present for negative c . Using the expansion of the kernel K 21 ,

$$K _ { 2 1 } ( p ) = - 2 \delta ( p ) + 2 K _ { 1 1 } ^ { ( 1 ) } ( p ) + \frac { 5 } { 4 } c ^ { 2 } \delta ^ { \prime \prime } ( p ) + \mathcal { O } ( c ^ { 3 } ) \,,$$

it can be shown that the perturbative expansion of f (3) indeed has the appropriate dependence on the sign of c at each order.

## 5.3 Bound states in the DMB approach

Now we will treat the two-particle sector of the attractive LL model from the DMB approach. This is arguably the simplest example involving bound states. The Hamiltonian is just

$$H = \frac { p ^ { 2 } } { m } + \frac { c } { m } \delta ( x ) + \frac { P ^ { 2 } } { 4 m },$$

where the the relative momentum p = ( p 1 -p 2 ) / 2 and position x = x 1 -x 2 form a canonically conjugate pair, and the total momentum P = p 1 + p 2 decouples. For the moment, we will ignore the decoupled P degree of freedom.

Using the notation | x 0 ⟩ for the state localized at x = 0, the delta function operator in the interaction may be written as ( c/m ) | x 0 ⟩⟨ x 0 | . The T-matrix operator (2.4) may be explicitly calculated

$$\hat { T } = \frac { c } { m - \langle x _ { 0 } | G _ { 0 } | x _ { 0 } \rangle } | x _ { 0 } \rangle \langle x _ { 0 } |,$$

and the trace-log of the S-matrix becomes

$$\text{Tr} \log \hat { S } = \log \left ( 1 + c \frac { \langle x _ { 0 } | G _ { 0 } - \bar { G } _ { 0 } | x _ { 0 } \rangle } { m - \langle x _ { 0 } | G _ { 0 } | x _ { 0 } \rangle } \right ).$$

For these calculations, note that ⟨ x 0 | p ⟩ = 1, and thus ∫ dp 2 π ⟨ p G | 0 | x 0 ⟩ = ⟨ x 0 | G x 0 | 0 ⟩ .

The expectation of the propagator is

$$\langle x _ { 0 } | G _ { 0 } | x _ { 0 } \rangle = \int \frac { d p } { 2 \pi } \frac { 1 } { E - \frac { v ^ { 2 } } { m } + i \epsilon } = - \frac { i m } { 2 \sqrt { m E ^ { + } } },$$

where √ mE + is defined to be the square root of mE + iϵ with positive imaginary part. Similarly, ⟨ x 0 | ¯ G x 0 | 0 ⟩ = im/ (2 √ mE -), where √ mE -is defined as the square root of mE -iϵ with negative imaginary part.

So (5.20) becomes,

$$\text{Tr} \log \hat { S } = \log \left ( 1 - \frac { i c / 2 } { \sqrt { m E ^ { + } } + i c / 2 } \left ( 1 + \frac { \sqrt { E ^ { + } } } { \sqrt { E ^ { - } } } \right ) \right ).$$

( √ )

In the limit of small ϵ , the factor 1 + E + √ E -evaluates to 2 for E &gt; 0, and vanishes for E &lt; 0. However if c is negative, then there is a pole at √ mE + = -ic/ 2, and this can compensate the vanishing of the factor. The argument of the logarithm winds around the complex plane in the vicinity of the pole, and this will lead to a delta function localized at the bound state energy after taking the E derivative in (1.1),

$$\frac { 1 } { 2 \pi i } \frac { d } { d E } \text{Tr} \log \hat { S } = \frac { d p } { d E } \frac { 1 } { 2 \pi } \frac { c } { p ^ { 2 } + c ^ { 2 } / 4 } \theta ( E ) + \delta \left ( E + \frac { c ^ { 2 } } { 4 m } \right ).$$

Let us now reinstate the total momentum P degree of freedom in the problem. All appearances of E become E -P / m 2 4 , and there is an overall factor of L and an additional P integral due to

the trace,

$$\delta Z = L \int d P \left [ e ^ { - \beta \left ( \frac { P ^ { 2 } } { 4 m } - \frac { c ^ { 2 } } { 4 m } \right ) } + \int d p \, e ^ { - \beta \left ( \frac { p ^ { 2 } } { m } + \frac { P ^ { 2 } } { 4 m } \right ) } \frac { 1 } { 2 \pi } \frac { 2 c } { 4 p ^ { 2 } + c ^ { 2 } } \right ].$$

After changing variables to ¯ p ≡ P/ 2 in the bound state contribution, to p , p 1 2 in the scattering contribution, and multiplying by -( βL ) -1 to obtain the free energy density, this takes the form of δf (2) calculated from the TBA in 5.15.

## 6 Discussion

In this paper we have investigated the DMB formula (1.1) for some massive scalar field theories in two spacetime dimensions. In particular we focused on integrable theories for which we may compare with the TBA or a large N expansion. For multi-particle scattering it was shown that there are two issues that complicate a naive treatment of the DMB formula. The first issue has to do with the presence of forward scattering divergences. This issue has been discussed from the outset of the DMB formula [1, 4, 5, 6], but we show here how by simply applying the DMB formula consistently the forward scattering divergences can be associated to derivatives of energy delta functions, and how the sum over all vertex time orderings corresponds to an easily calculable loop diagram in the thermal field theory approach.

The second issue is perhaps specific to two spacetime dimensions and has not to my knowledge been previously discussed in the literature. The issue has to do with the fact that the trace in the DMB formula has large contributions from states with collinear momenta. This issue arises already in a calculation of the melon diagram (A.3) at second order in the ϕ 4 theory. It is not just a difficulty of the DMB approach, and calculations of the melon diagram using thermal field theory or the TBA also involve subtleties at collinear momentum. The lesson is that regions of integration with collinear momenta should also be considered carefully for higher order diagrams, and for more complicated 1+1 dimensional theories.

There have not been many concrete examples of application of the DMB approach, so there are a number of possible directions to build on the initial investigations in this paper. With an eye towards the effective string theory program, clearly it would be illuminating to consider scattering in massless 1+1 dimensional theories. Also, one advantage of the TBA for that program is that higher energy levels in finite volume may be calculated [33], not just the ground state. It would be interesting to consider higher energy levels in the DMB approach as well. Something similar has already been done for the thermal self-energy [10], which gives information on the first excited state in finite volume. Indeed in section 4 the energy of first excited state M β ( ) was shown to be given at lowest order by the two-body forward scattering amplitude, as in the work of L¨scher [7, 8], and it u is straightforward to extend this to multi-particle scattering by calculating higher order corrections to the thermal self-energy using the rules of Sec 4.1.

We have mostly considered theories with one type of particle in the spectrum which corresponds neatly to the field in the Lagrangian, but of course the possibilities in quantum field theory are much richer. The particles of the free theory may form bound states or may decay into other stable

particles when an interaction is included. We have taken a preliminary look at bound states in the Lieb-Liniger model and seen that they may indeed be calculated in the DMB approach. The analysis of Sec 5 could be further extended to treating the breather bound states in the sine-Gordon model, although that model is more complicated due to the presence of solitons, and relatedly the presence of non-trivial semi-classical saddles in the thermal field theory approach.

Of course we may avoid these problems for sine-Gordon by considering the dual massive Thirring model instead. Treating fermions in the DMB approach is no major difficulty (see e.g. [6]), and although the calculation was not included here, indeed the lowest order correction to the free energy may be derived from the massive Thirring two-body forward scattering amplitudes. Even so, it may be enlightening to consider the bosonic description of the theory and investigate how the non-perturbative corrections due to the solitons can possibly arise in the DMB framework.

## Acknowledgements

I would like to thank Vladimir Rosenhaus for introducing reference [1], and collaboration on an early version of this project. I would also like to thank Bal´zs Pozsgay for discussions. a This work is supported in part by NSF grant PHY-2209116 and by the ITS through a Simons grant.

## A More on the melon diagram

This section deals more with the melon diagram integral M 4 ≡ ∫ dτdxD τ,x ( ) 4 , where D τ,x ( ) is the finite temperature propagator (with 0 ≤ τ &lt; β )

$$D ( \tau, x ) & = \frac { 1 } { \beta } \sum _ { n } \int \frac { d k } { 2 \pi } \frac { e ^ { - i ( k x + \omega _ { n } \tau ) } } { \omega _ { n } ^ { 2 } + k ^ { 2 } + m ^ { 2 } } \\ & = \int \frac { d k } { 2 \pi } e ^ { - i k x } \left ( \frac { e ^ { - E | \tau | } + e ^ { - E ( \beta - | \tau | ) } } { 2 E \left ( 1 - e ^ { - \beta E } \right ) } \right ).$$

For a theory with interaction term m g 2 2 4! ϕ 4 , the contribution to the free energy density is

$$f _ { \text{melon} } = - \frac { 1 } { 3 } \frac { m ^ { 4 } g ^ { 4 } } { 1 6 } M _ { 4 } - \text{counterterms},$$

where the counterterms subtract the zero-body and one-body parts.

The full result for M 4 is

$$M _ { 4 } = \frac { 7 } { 8 m ^ { 2 } } \frac { \zeta ( 3 ) } { ( 2 \pi ) ^ { 3 } } + \frac { 1 } { 1 6 m ^ { 2 } } \int \frac { d \theta } { 2 \pi } \left [ n ( \theta ) + 3 n ( \theta ) ^ { 2 } + 2 n ( \theta ) ^ { 3 } \right ].$$

The two-body part was calculated in (3.14). The three-body part (3.19) was calculated indirectly via the TBA in Sec 5.1, and will be calculated directly here.

## A.1 Zero- and one-body parts

For completeness we will also present the parts of the integral M 4 that are canceled by counterterms. At zero temperature ( β →∞ ), the propagator is just given by a Bessel function,

$$D _ { \infty } ( \tau, x ) = \frac { 1 } { 2 \pi } K _ { 0 } ( m \sqrt { \tau ^ { 2 } + x ^ { 2 } } ).$$

The 'zero-body' part is equivalent to the zero-temperature limit,

$$\mu ^ { - } _ { \mu } \mu ^ { - } _ { \mu } = \frac { 1 } { m ^ { 2 } ( 2 \pi ) ^ { 3 } } \int _ { 0 } ^ { \infty } d u \, u \, K _ { 0 } ( u ) ^ { 4 } = \frac { 7 } { 8 m ^ { 2 } } \frac { \zeta ( 3 ) } { ( 2 \pi ) ^ { 3 } }.$$

The one-body part involves the self-energy diagram given by cutting one line of the melon diagram. This self-energy diagram is often referred to as the 'sunset' diagram and it is given by

$$\lambda \, \stackrel { \sim } { \Pi } _ { \text{sunset} } ( p ) = \int d ^ { 2 } x e ^ { i p \cdot x } D _ { \infty } ( x ) ^ { 3 } = \frac { 1 } { ( 2 \pi ) ^ { 2 } } \int _ { 0 } ^ { \infty } d r \, r J _ { 0 } \left ( | p | r \right ) \left ( K _ { 0 } ( m r ) \right ) ^ { 3 }. \, \quad \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad$$

The one-body part of M 4 is then given by the on-shell ( p 2 = -m 2 ) value of the sunset diagram,

$$\leftarrow & \text{body part of } M _ { 4 } \text{ is then given by the on-shell } ( p ^ { \epsilon } = - m ^ { \epsilon } ) \text{ value of the sunset diagram,} \\ M _ { 4, \, 1 } & = \int \frac { d \theta } { 2 \pi } n ( \theta ) \, 4 \Pi _ { \text{sundset} } | _ { p ^ { 2 } = - m ^ { 2 } } \\ & = \int \frac { d \theta } { 2 \pi } n ( \theta ) \left ( \frac { 4 } { m ^ { 2 } \, ( 2 \pi ) ^ { 2 } } \int _ { 0 } ^ { \infty } d u \, u \, I _ { 0 } \, ( u ) \, ( K _ { 0 } ( u ) ) ^ { 3 } \right ) \\ & = \frac { 1 } { 1 6 m ^ { 2 } } \int \frac { d \theta } { 2 \pi } n ( \theta ). \\ \text{...} \quad \text{$-} \quad \text{$-} \quad \text{$-} \quad \text{$-}$$

## A.2 Cutting the melon diagram

For the sake of finding the two- and three-body parts, cutting the melon diagram according to the method of Sec 3.2 is not very convenient due to the singularities in M π . In Sec 3.3 we nevertheless considered the forward scattering amplitudes corresponding to cutting a line of the melon diagram, and regularized them by the principal value prescription suggested by the secondorder DMB formula (2.26).

In this subsection we will justify this directly from the integral M 4 . This involves integrating the vertices of the diagram over Euclidean time, rather than summing over Matsubara frequencies. There is a generalization of this method [6, 9] based on the work of Baym and Sessler [34], but M 4 is particularly simple, and it has been considered in [23, 24].

The propagator (A.1) may be written

$$D ( \tau, x ) = \int \frac { d k } { 2 \pi } e ^ { i k \cdot x } \int \frac { d \omega } { 2 E } \left ( \delta ( \omega - E ) - \delta ( \omega + E ) \right ) \left ( \frac { e ^ { - \omega | \tau | } } { 1 - e ^ { - \beta \omega } } \right ),$$

and a melon diagram M N with N propagators going between two vertices may be calculated,

$$M _ { N } \equiv \int d \tau d x D ( \tau, x ) ^ { N } = \int \frac { d k _ { 1 } } { 2 \pi } \dots \frac { d k _ { N } } { 2 \pi } 2 \pi \delta \left ( \sum _ { i } k _ { i } \right ) \tilde { M } _ { N },$$

$$\tilde { M } _ { N } = \int \prod _ { j = 1 } ^ { N } \left ( \frac { d \omega _ { j } } { 2 E _ { j } } \frac { \delta ( \omega _ { j } - E _ { j } ) - \delta ( \omega _ { j } + E _ { j } ) } { 1 - e ^ { - \beta \omega _ { j } } } \right ) \frac { \cdot \, \cdot \, \cdot } { \sum _ { i } \omega _ { i } }.$$

Let us focus on ˜ M 4 ,

$$\tilde { M } _ { 4 } = 2 \left ( \prod _ { j = 1 } ^ { 4 } \frac { n _ { j } } { 2 E _ { j } } \right ) \left [ \frac { e ^ { \beta \sum _ { i } E _ { i } } - 1 } { \sum _ { i } E _ { i } } + 4 \frac { e ^ { \beta ( E _ { 1 } + E _ { 2 } + E _ { 3 } ) } - e ^ { \beta E _ { 4 } } } { E _ { 1 } + E _ { 2 } + E _ { 3 } - E _ { 4 } } + 3 \frac { e ^ { \beta ( E _ { 1 } + E _ { 2 } ) } - e ^ { \beta ( E _ { 3 } + E _ { 4 } ) } } { E _ { 1 } + E _ { 2 } - E _ { 3 } - E _ { 4 } } \right ]. \quad ( \text{A.11} )$$

The final term of ˜ M 4 is actually analytic, since the numerator vanishes when E 1 + E 2 = E 3 + E 4 . Thus it is valid to adjust the denominator by an arbitrary infinitessimal ± iϵ , after which we may exchange integration variables 1 , 2 and 3 4, ,

$$3 \frac { e ^ { \beta ( E _ { 1 } + E _ { 2 } ) } - e ^ { \beta ( E _ { 3 } + E _ { 4 } ) } } { E _ { 1 } + E _ { 2 } - E _ { 3 } - E _ { 4 } \pm i \epsilon } \rightarrow \mathcal { P } \frac { 6 \, e ^ { \beta ( E _ { 1 } + E _ { 2 } ) } } { E _ { 1 } + E _ { 2 } - E _ { 3 } - E _ { 4 } }.$$

This is ultimately how the principal value prescription (2.26) arises in the thermal field theory approach.

The exponentials in the numerators may be written as e βE j = n -1 j + 1 and so some of the occupation numbers in the product ∏ j n j will be canceled. The r -body part of the diagram is just taken to be the part of the diagram with r overall occupation number factors. The three-body part is

$$\tilde { M } _ { 4, 3 } = \left ( \prod _ { j = 1 } ^ { 3 } \frac { n _ { j } } { E _ { j } } \right ) \left [ \frac { 1 } { - E _ { 4 } ^ { 2 } + ( E _ { 1 } + E _ { 2 } + E _ { 3 } ) ^ { 2 } } + \mathcal { P } \frac { 3 } { - E _ { 4 } ^ { 2 } + ( E _ { 1 } + E _ { 2 } - E _ { 3 } ) ^ { 2 } } \right ].$$

This obviously agrees with the DMB expression (3.15).

The two-body part is

$$\tilde { M } _ { 4, \, 2 } = \frac { 3 } { 4 } \frac { n _ { 1 } n _ { 2 } } { E _ { 1 } E _ { 2 } E _ { 3 } E _ { 4 } } \left [ \frac { 1 } { \sum _ { i } E _ { i } } + \frac { 2 } { E _ { 3 } + E _ { 4 } + E _ { 1 } - E _ { 2 } } + \mathcal { P } \frac { 1 } { E _ { 3 } + E _ { 4 } - E _ { 1 } - E _ { 2 } } \right ].$$

This also agrees with the DMB expression (3.13) after we write it in terms of the integral J p ( ), which may be written in terms of energy denominators

$$J ( p ) = \int \frac { d k _ { 3 } } { 2 \pi } \frac { d k _ { 4 } } { 2 \pi } \frac { 2 \pi \delta ( k _ { 3 } + k _ { 4 } - p ^ { 1 } ) } { 2 E _ { 3 } \, 2 E _ { 4 } } \left [ \frac { 1 } { i p ^ { 0 } + E _ { 3 } + E _ { 4 } } + \frac { 1 } { - i p ^ { 0 } + E _ { 3 } + E _ { 4 } } \right ].$$

## A.3 Direct calculation of the three-body part

This section builds on the discussion of Sec 3.3.3, but the iϵ terms are considered carefully.

For brevity introduce notation such as

$$e _ { 1 2 } \equiv e ^ { \theta _ { 1 } } - e ^ { \theta _ { 2 } }, \quad \delta _ { 1 2 } \equiv 2 \pi \delta ( e _ { 1 2 } ), \quad H _ { 1 2 } \equiv \frac { \text{sgn} ( e _ { 1 2 } ) + 1 } { 2 }.$$

As usual the Heaviside function H 12 is taken to equal 1 / 2 when θ 1 = θ 2 .

From (3.16), the principal value of the scattering amplitude corresponding to the middle diagram on the lower row of Fig 4 is

$$\mathcal { P } \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) e _ { 1 3 } e _ { 2 3 } } & = \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) } \text{Re} \frac { 1 } { \left ( e _ { 1 3 } + i \text{sgn} \left ( e _ { 2 3 } \right ) \epsilon \right ) \left ( e _ { 2 3 } + i \text{sgn} \left ( e _ { 1 3 } \right ) \epsilon \right ) } \\ & = \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) } \text{Re} \left [ \left ( \frac { 1 } { e _ { 1 3 } + i \epsilon } + i \text{H} _ { 1 3 } \delta _ { 1 3 } \right ) \left ( \frac { 1 } { e _ { 2 3 } + i \epsilon } + i \text{H} _ { 1 3 } \delta _ { 2 3 } \right ) \right ] \\ & = \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { m ^ { 2 } \left ( e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } \right ) } \left [ \left ( \frac { 1 } { e _ { 1 3 } + i \epsilon } \right ) \left ( \frac { 1 } { e _ { 2 3 } + i \epsilon } \right ) + \frac { 1 } { 4 } \delta _ { 1 3 } \delta _ { 2 3 } + \frac { i } { 2 } \left ( \delta _ { 1 3 } \mathcal { P } \frac { 1 } { e _ { 2 3 } } + \delta _ { 2 3 } \mathcal { P } \frac { 1 } { e _ { 1 3 } } \right ) \right ].$$

The last terms involving a single delta function will vanish upon integration since for instance,

$$\frac { e ^ { \theta _ { 1 } + \theta _ { 2 } + \theta _ { 3 } } } { e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } } \delta _ { 1 3 } \mathcal { P } \frac { 1 } { e _ { 2 3 } } = \frac { e ^ { \theta _ { 1 } + \theta _ { 2 } } } { e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } } 2 \pi \delta ( \theta _ { 1 } - \theta _ { 3 } ) \mathcal { P } \frac { e ^ { - \frac { \theta _ { 1 } + \theta _ { 2 } } { 2 } } } { 2 \sinh \frac { \theta _ { 1 } - \theta _ { 2 } } { 2 } } = 2 \pi \delta ( \theta _ { 1 } - \theta _ { 3 } ) \mathcal { P } \frac { 1 } { 2 \sinh \left ( \theta _ { 1 } - \theta _ { 2 } \right ) },$$

and this is an odd under the interchange of the signs of θ , θ 1 2 .

The first term in (A.16) will be symmetrized in order to cancel in the same way as in (3.16), but since the signs of the iϵ terms must be adjusted to be the same this will also produce some additional terms,

$$& \frac { 1 } { 3 } \left [ \frac { 1 } { e ^ { \theta _ { 1 } } + e ^ { \theta _ { 2 } } } \left ( \frac { 1 } { e _ { 1 3 } + i \epsilon } \right ) \left ( \frac { 1 } { e _ { 2 3 } + i \epsilon } \right ) + ( 3 \leftrightarrow 1 ) + ( 3 \leftrightarrow 2 ) \right ] \\ & \quad = \frac { 1 } { 3 } \left [ \dots - \frac { 1 } { e ^ { \theta _ { 3 } } + e ^ { \theta _ { 2 } } } i \tilde { \delta } _ { 1 3 } \left ( \frac { 1 } { e _ { 2 1 } + i \epsilon } \right ) + \frac { 1 } { e ^ { \theta _ { 1 } } + e ^ { \theta _ { 3 } } } \left ( i \tilde { \delta } _ { 1 2 } \left ( \frac { 1 } { e _ { 2 3 } + i \epsilon } \right ) + i \tilde { \delta } _ { 2 3 } \left ( \frac { 1 } { e _ { 2 1 } + i \epsilon } \right ) - \tilde { \delta } _ { 1 2 } \tilde { \delta } _ { 2 3 } \right ) \right ] \\ & \quad = \frac { 1 } { 3 } \left [ \dots - \frac { 1 } { 2 } \frac { 1 } { e ^ { \theta _ { 3 } } + e ^ { \theta _ { 2 } } } \tilde { \delta } _ { 1 3 } \tilde { \delta } _ { 1 2 } \right ], \\ & \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

where the ellipsis includes both the terms which will cancel as in (3.16) and the terms with a single delta function that were shown to vanish immediately above.

So in total the non-vanishing part of M 4 3 , in (A.12) involves a double delta function,

$$M _ { 4, 3 } = \int \frac { d \theta _ { 1 } d \theta _ { 2 } d \theta _ { 3 } } { ( 2 \pi ) ^ { 3 } } n _ { 1 } ^ { 3 } \frac { 3 e ^ { 3 \theta _ { 1 } } } { 2 m ^ { 2 } e ^ { \theta _ { 1 } } } \left ( \frac { 1 } { 4 } - \frac { 1 } { 6 } \right ) \delta _ { 1 3 } \delta _ { 2 3 } = \frac { 2 } { 1 6 m ^ { 2 } } \int \frac { d \theta } { 2 \pi } n ( \theta ) ^ { 3 }.$$

## B More on large N diagrams

## B.1 An expansion of x s

Recall that x s is defined as a sum of Bessel functions,

$$x _ { s } \equiv \sum _ { j = 1 } ^ { \infty } \left ( j \zeta \right ) ^ { s } K _ { s } ( j \zeta ), \quad \zeta \equiv \beta m.$$

There is a formula for x s in terms of energy derivatives of n which was used in Sec 4 to relate thermal field theory diagrams to the DMB formula,

$$x _ { s } = \int \frac { d k } { 2 E } \frac { m ^ { 2 s } } { E ^ { s } } \sum _ { l = 0 } ^ { s } \frac { ( s + l )! } { ( 2 E ) ^ { l } ( s - l )! l! } \left ( - \frac { d } { d E } \right ) ^ { s - l } n ( E ).$$

This formula may be proven inductively. The s = 0 case involves a well-known integral expression for Bessel functions

$$x _ { 0 } = \sum _ { j = 1 } ^ { \infty } K _ { 0 } ( j \zeta ) = \sum _ { j = 1 } ^ { \infty } \frac { 1 } { 2 } \int d \theta \, e ^ { - j \zeta \cosh \theta } = \int \frac { d k } { 2 E } \, n ( E ).$$

To prove (B.2) for x s +1 , first note that 6

$$- \zeta ^ { 2 s + 1 } \frac { d } { d \zeta } \zeta ^ { - 2 s } x _ { s } = \sum _ { j = 1 } ^ { \infty } ( j \zeta ) ^ { s + 1 } \left ( \frac { s } { j \zeta } - \frac { d } { d ( j \zeta ) } \right ) K _ { s } ( j \zeta ) = x _ { s + 1 }.$$

The same operator may be written as a derivative of m 2 at fixed β ,

$$- \zeta ^ { 2 s + 1 } \frac { d } { d \zeta } \zeta ^ { - 2 s } = - 2 m ^ { 2 s + 2 } \frac { d } { d m ^ { 2 } } m ^ { - 2 s }.$$

Now act with this operator on both sides of (B.2),

$$x _ { s + 1 } = - m ^ { 2 s + 2 } \int d k \frac { d } { d m ^ { 2 } } \left [ \frac { 1 } { E ^ { s + 1 } } \sum _ { l = 0 } ^ { s } \frac { ( s + l )! } { ( 2 E ) ^ { l } ( s - l )!! } \left ( - \frac { d } { d E } \right ) ^ { s - l } n ( E ) \right ].$$

Inside the integrand, the operator d dm 2 = 1 2 E dE d for fixed k . It is now a simple exercise to show that the acting with d/dE reproduces the terms of (B.2) for x s +1 , completing the proof.

Also note that the summation over j involved in x s is not essential to the proof. If n E ( ) is expanded like ∑ j =1 e -jβE , then the formula holds order by order in j . In particular it is true at the leading order in fugacity, j = 1.

## B.2 Calculation of f α 4

Recall from Sec 4 that the free energy f α p due to p +1 particles scattering in the manner of Fig 7 is given by (4.17),

$$f _ { \alpha ^ { p } } = - \frac { N m ^ { 2 } } { \pi } \frac { ( - 2 \Lambda ) ^ { p } } { p! } x _ { 0 } ^ { p } x _ { p - 1 }.$$

The leading order ( l = 0) term of this formula was matched to the DMB formula for arbitrary p in Sec 4.2. In this section we will fix p = 4, but match the full l expansion in (B.2) to the DMB formula. We will continue to work at lowest order in fugacity, corresponding to the leading order ( j = 1) terms in all appearances of x s .

6 Incidentally, the s = 0 case of this formula can be used to make the simplification (3.11).

f α 4 involves the factor x 3 , which is given by (B.2),

$$x _ { 3 } ^ { ( 1 ) } = m ^ { 6 } \int \frac { d k ^ { \prime } } { 2 E ^ { \prime } } \frac { e ^ { - \beta E ^ { \prime } } } { E ^ { \prime 3 } } \left ( \beta ^ { 3 } + \frac { 6 \beta ^ { 2 } } { E ^ { \prime } } + \frac { 1 5 \beta } { E ^ { \prime 2 } } + \frac { 1 5 } { E ^ { \prime 3 } } \right ).$$

Since the overall factors in f α 4 have already been established for the leading order β 3 term in Sec 4.2, it is sufficient here to reproduce the expansion in parentheses.

We can expand log ˆ in the DMB formula up to fourth order. Each term will take the schematic S form

$$\Delta G V O _ { 3 } V O _ { 2 } V O _ { 1 } V$$

where ∆ G ≡ G 0 -¯ G 0 , and the operator O i at position i represents either a propagator G 0 or a delta function ∆ G . When the trace is taken, the H 0 in leftmost copy of ∆ G will evaluate to the energy of the incoming or outgoing state E α , but the three operators O , O , O 1 2 3 may involve some other intermediate energy which will depend on the 'time ordering' in which the four copies of V are contracted (see e.g. Fig 6 for the p = 2 case).

There are 4! = 24 total time orderings of the vertices of the forward scattering diagram corresponding to f α 4 .

- · Of these, only the single ordering considered in Sec 4.2.2 will involve all three intermediate energies evaluating to E α . According to (4.18) this leads to a term proportional to

$$- \frac { 1 } { 4! } \frac { d ^ { 3 } } { d E ^ { 3 } } \delta \left ( E - E _ { \alpha } \right ).$$

- · There are 3 time orderings where two intermediate energies are on-shell, and the remaining energy is E α + 2 E ′ , where E ′ is the energy of the privileged particle in Fig 7. Any delta functions ∆ G vanish if they are evaluated at an 'off-shell' energy (i.e. an energy not equal to E α ), and the terms in the fourth-order expansion of log ˆ that survive are just proportional S to a third-order expansion of log ˆ with all intermediate energies set to S E α . So using (4.18) again, these 3 time orderings lead to a term

$$+ 3 \left ( \frac { 1 } { E - E _ { \alpha } - 2 E ^ { \prime } } \right ) \frac { 1 } { 3! } \frac { d ^ { 2 } } { d E ^ { 2 } } \delta \left ( E - E _ { \alpha } \right ).$$

- · Proceeding similarly, there are 7 time orderings with one energy equal to E α and two equal to E α +2 E ′ ,

$$- 7 \left ( \frac { 1 } { E - E _ { \alpha } - 2 E ^ { \prime } } \right ) ^ { 2 } \frac { 1 } { 2! } \frac { d } { d E } \delta \left ( E - E _ { \alpha } \right ).$$

- · There are 9 time orderings where all energies equal E α +2 E ′ , and 4 time orderings where two energies equal E α +2 E ′ , and one energy equals E α +4 E ′ ,

$$+ \left [ 9 \left ( \frac { 1 } { - 2 E ^ { \prime } } \right ) ^ { 3 } + 4 \left ( \frac { 1 } { - 2 E ^ { \prime } } \right ) ^ { 2 } \left ( \frac { 1 } { - 4 E ^ { \prime } } \right ) \right ] \delta \left ( E - E _ { \alpha } \right ).$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00729v1-figure-fb2dc90c-871c-4163-8371-df75f889db0d.jpg)

1

2

Figure 8: The 'bubble' and 'melon' corrections to the free energy for the Lieb-Liniger model.

All of these terms are added together in the density of states, where they multiply e -βE in the integrand of (1.1). After integrating the derivatives of delta functions by parts (being careful to also act on the appearance of E in the off-shell propagator), the sum of the above terms becomes

$$- \frac { 1 } { 4! } \left ( \beta ^ { 3 } + \frac { 6 \beta ^ { 2 } } { E ^ { \prime } } + \frac { 1 5 \beta } { \left ( E ^ { \prime } \right ) ^ { 2 } } + \frac { 1 5 } { \left ( E ^ { \prime } \right ) ^ { 3 } } \right ) \delta \left ( E - E _ { \alpha } \right ).$$

which gives the factor in parentheses in (B.6), which we set out to show.

## C The Lieb-Liniger model

In the text the discussion has mostly focused on diagrams in the relativistic ϕ 4 and sinh-Gordon models, but much the same story occurs for the non-relativistic Lieb-Liniger (LL) model, where the two-body interaction involves a repulsive delta function potential of strength c/m .

We may calculate with a coherent state path integral with Lagrangian

$$\mathcal { L } _ { L L } = \sum _ { k } a _ { k } ^ { \dagger } \left ( \partial _ { \tau } + E _ { k } \right ) a _ { k } + \frac { c } { 2 m } \sum _ { k + l = m + n } a _ { m } ^ { \dagger } a _ { n } ^ { \dagger } a _ { k } a _ { l },$$

where the indices k on the c-number field a k refer to spatial momentum. No counterterms are needed as long as the a † fields in the interaction vertex are understood to be at an infinitesimally later time than the a fields. The propagators ⟨ a τ ( ) a † (0) ⟩ (0) are directed, and at zero-temperature the propagator vanishes for negative τ . This greatly simplifies the perturbation theory since at zero-temperature all non-vanishing diagrams must involve arrows pointing in the same direction.

The diagrams correcting the free energy at O ( c 2 ) are shown in Fig 8. The strict requirements on the propagator at zero-temperature means that the bubble and melon diagrams correspond to far fewer scattering amplitudes than in Fig 4.

The free energy density due to the bubble diagram may be calculated straightforwardly in thermal field theory,

$$f _ { \text{bubble} } = - 2 \beta c ^ { 2 } \left ( \int \frac { d p } { 2 \pi } \frac { n ( p ) } { m } \right ) ^ { 2 } \int \frac { d q } { 2 \pi } \frac { e ^ { - \beta E _ { q } } } { \left ( 1 - e ^ { - \beta E _ { q } } \right ) ^ { 2 } }.$$

From the DMB perspective this involves only the on-shell divergent three-body scattering amplitude depicted on the lower left in Fig 8. The result is a straightforward non-relativistic limit of the relativistic result (3.11).

Note that the bubble diagram in LL does not have a two-body part. Recall that in the relativistic theory the two-body part of the melon was associated with expanding the first-order in c correction to the free energy density to the order g 4 , but in the present case there is no g 2 parameter at all.

As in the relativistic case, the bubble diagram may be found from the three-body part of the TBA (5.2) if the principal value signs in the first order expansion of the kernel are ignored. The free energy density calculated from the TBA is

$$f _ { 3 } ^ { ( 2 ) } = - 2 \beta \frac { c ^ { 2 } } { m ^ { 2 } } \int \frac { d p _ { 1 } d p _ { 2 } d p _ { 3 } } { ( 2 \pi ) ^ { 3 } } e ^ { \beta E _ { 1 } } n _ { 1 } ^ { 2 } n _ { 2 } n _ { 3 } \left [ p _ { 2 } p _ { 3 } \mathcal { P } \frac { 1 } { p _ { 1 2 } } \mathcal { P } \frac { 1 } { p _ { 1 3 } } + \dots _ { 3 } \right ].$$

In this case quantities such as p 12 should be understood as p 1 -p 2 , and once again the identity (5.8) holds.

The melon diagram involves a single two-body and a single three-body scattering amplitude depicted on the lower right of Fig 8. These can either be calculated directly, using the method of Appendix A, or indirectly from the perturbative expansion of the TBA using (5.1) and (C.3). The result is a straightforward non-relativistic limit of the relativistic case (3.14) and (3.19),

$$f _ { \text{melon} } = - \frac { c ^ { 2 } } { 1 2 m } \int \frac { d p } { 2 \pi } \left ( 3 n ( p ) ^ { 2 } + 2 n ( p ) ^ { 3 } \right ).$$

## References

- [1] R. Dashen, S. Ma, and H. J. Bernstein, 'S-matrix formulation of statistical mechanics,' Phys. Rev. 187 (1969) 345.
- [2] C. N. Yang and C. P. Yang, 'Thermodynamics of a one-dimensional system of bosons with repulsive delta-function interaction,' J. Math. Phys. 10 (1969) 1115.
- [3] E. Beth and G. E. Uhlenbeck, 'The quantum theory of the non-ideal gas. II. Behaviour at low temperatures,' Physica 4 (1937) 915.
- [4] R. Dashen and S. Ma, 'Singular three-body amplitudes in the theory of the third virial coefficient,' J. Math. Phys. 11 (1970) 1136.
- [5] R. Dashen and S. Ma, 'Singularities in forward multiparticle scattering amplitudes and the S-matrix interpretation of higher virial coefficients,' J. Math. Phys. , 11 (1971) 689.
- [6] R. E. Norton, 'Elementary particle scattering and statistical quasi-particles in quantum statistical mechanics,' Annals of Phys. 170 (1986) 18.
- [7] M. L¨ uscher, 'On a relation between finite size effects and elastic scattering processes,' in Progress in Gauge Field Theory , (1984) 451.

- [8] M. L¨ uscher, 'Volume dependence of the energy spectrum in massive quantum field theories. I. Stable particle states,' Commun. Math. Phys. 104 (1986) 177.
- [9] S. Jeon, 'Computing spectral densities in finite temperature field theory,' Phys. Rev. D 47 (1993) 4586, hep-ph/9210227.
- [10] S. Jeon and P. J. Ellis, 'Multiple scattering expansion of the self-energy at finite temperature,' Phys. Rev. D 58 (1998) 045013, hep-ph/9802246.
- [11] A. Athenodorou, B. Bringoltz, and M. Teper, 'Closed flux tubes and their string description in D = 3 + 1 SU N ( ) gauge theories,' JHEP (2011) 030, 1007.4720.
- [12] A. Athenodorou, B. Bringoltz, and M. Teper, 'Closed flux tubes and their string description in D = 2 + 1 SU N ( ) gauge theories,' JHEP (2011) 042, 1103.5854.
- [13] O. Aharony and Z. Komargodski, 'The effective theory of long strings,' JHEP (2013) 118, 1302.6257.
- [14] S. Dubovsky, R. Flauger, and V. Gorbenko, 'Evidence from lattice data for a new particle on the worldsheet of the QCD flux tube,' Phys. Rev. Lett. 111 (2013) 062006, 1301.2325.
- [15] S. Dubovsky, R. Flauger, and V. Gorbenko, 'Flux tube spectra from approximate integrability at low energies,' J. Exp. Theor. Phys. 120 (2015) 399, 1404.0037.
- [16] P. Cooper, S. Dubovsky, V. Gorbenko, A. Mohsen, and S. Storace, 'Looking for integrability on the worldsheet of confining strings,' JHEP (2015) 127, 1411.0703.
- [17] C. Chen, P. Conkey, S. Dubovsky, and G. Hernandez-Chifflet, 'Undressing confining flux tubes with TT ,' Phys. Rev. D 98 (2018) 114024, 1808.01339.
- [18] M. L¨ uscher, 'Volume dependence of the energy spectrum in massive quantum field theories: II. Scattering states,' Commun. in Math. Phys. 105 (1986) 153.
- [19] M. L¨ uscher and U. Wolff, 'How to calculate the elastic scattering matrix in two-dimensional quantum field theories by numerical simulation,' Nucl. Phys. B 339 (1990) 222.
- [20] M. T. Hansen and S. R. Sharpe, 'Relativistic, model-independent, three-particle quantization condition,' Phys. Rev. D 90 (2014) 116003, 1311.4848.
- [21] G. Kato and M. Wadati, 'Graphical representation of the partition function of a onedimensional δ -function Bose gas,' J. Math. Phys. 42 (2001) 4883, cond-mat/0212323.
- [22] P. Dorey, 'Exact S matrices,' in Eotvos Summer School in Physics: Conformal Field Theories and Integrable Models , (1996) 85, hep-th/9810026
- [23] A. Bugrij and V. Shadura, 'Three-loop contributions to the free energy of λϕ 4 QFT,' (1995) hep-th/9510232

- [24] J. O. Andersen, E. Braaten, and M. Strickland, 'Massive basketball diagram for a thermal scalar field theory,' Phys. Rev. D 62 (2000) 045004, hep-ph/0002048.
- [25] J. Frenkel, A. Saa, and J. Taylor, 'Pressure in thermal scalar field theory to three-loop order,' Phys. Rev. D 46 (1992) 3670.
- [26] G. Barton, 'On the finite-temperature quantum electrodynamics of free electrons and photons,' Annals Phys.s 200 (1990) 271.
- [27] E. H. Lieb and W. Liniger, 'Exact analysis of an interacting bose gas. I. The general solution and the ground state,' Phys. Rev. 130 (1963) 1605.
- [28] S. N. Vergeles and V. M. Gryanik, 'Two-dimensional quantum field theories having exact solutions,' Sov. J. Nucl. Phys. 23 (1976) 704.
- [29] J. Balog and A. Heged˝s, 'Virial expansion and TBA in ´ u O N ( ) sigma-models,' Phys. Lett. B 523 (2001) 211, hep-th/0108071
- [30] M. L¨ uscher, 'A new method to compute the spectrum of low-lying states in massless asymptotically free field theories,' Phys. Lett. B 118 (1982) 391.
- [31] B. C. Nagy, G. Tak´ acs, and M. Kormos, 'Thermodynamic Bethe Ansatz and generalised hydrodynamics in the sine-Gordon model,' SciPost Phys. 16 (2024) 145, 2312.03909.
- [32] S. J. van Tongeren, 'Introduction to the thermodynamic Bethe ansatz,' J. Phys. A: Math. Theor. 49 (2016) 323005, 1606.02951.
- [33] P. Dorey and R. Tateo, 'Excited states by analytic continuation of TBA equations,' Nucl. Phys. B 482 (1996) 639, hep-th/9607167.
- [34] G. Baym and A. M. Sessler, 'Perturbation-theory rules for computing the self-energy operator in quantum statistical mechanics,' Phys. Rev. 131 (1963) 2345.
Daniel Schubring
Free energy from forward scattering in 1+1d
The free energy, or equivalently the ground state energy in finite volume,
may be calculated from forward scattering amplitudes using a formula due to
Dashen, Ma, and Bernstein. However a naive treatment leads to singularities
when considering the scattering of three or more particles. It is shown in
detail how the approach can be applied to multi-particle scattering in various
massive scalar theories in 1+1d, with or without integrability. The results for
the sinh-Gordon, Lieb-Liniger, and $O(N)$ non-linear sigma models are compared
to exact results. It is shown how bound states can be considered in this
approach by considering the attractive Lieb-Liniger model.