Normal ordering normal modes

In a soliton sector of a quantum field theory, it is often convenient to expand the quantum fields in terms of normal modes. Normal mode creation and annihilation operators can be normal ordered, and their normal ordered products have vanishing expectation values in the one-loop soliton ground state. The Hamiltonian of the theory, however, is usually normal ordered in the basis of operators which create plane waves. In this paper we find the Wick map between the two normal orderings. For concreteness, we restrict our attention to Schrodinger picture scalar fields in 1+1 dimensions, although we expect that our results readily generalize beyond this case. We find that plane wave ordered n-point functions of fields are sums of terms which factorize into j-point functions of zero modes, breather and continuum normal modes. We find a recursion formula in j and, for products of fields at the same point, we solve the recursion formula at all j.


Introduction
In perturbation theory about a translation-invariant vacuum, it is customary to decompose the quantum fields into operators a † p and a p which create and annihilate plane wave excitations. The free vacuum is annihilated by a p and is the initial state in the perturbative expansion. This perturbation theory is simplest when the Hamiltonian is normal ordered, so that all a † p appear to the left of all a p . At the same leading order, 1 the ground state of a quantum soliton is given by a coherent state formed by shifting the fields by the functions corresponding to their classical solutions [3][4][5]. The normal modes of the quantum soliton are, at linear order, described by quantum harmonic oscillators. The one-loop ground state of the soliton sector consists of the ten- 1 Even at weak coupling quantum corrections can affect the existence itself of the solution [1,2]. a e-mail: jarah@impcas.ac.cn (corresponding author) sor product of the ground states of these oscillators [6]. If the fields are decomposed into the normal modes of the soliton, with operators b † k and b k corresponding to the raising and lowering operators in the corresponding quantum harmonic oscillators, then the one-loop soliton ground state is, after the shift operator noted above, the state annihilated by all of the b k .
Again a sensible perturbation theory exists which describes the spectrum of the one soliton sector. It is similar to that of the vacuum sector, except that treating zero modes requires special care [7][8][9]. In particular the calculation is simplest if the Hamiltonian is normal ordered by placing all b † k to the left of b k .
However, the Hamiltonian is usually given normal ordered in terms of plane waves, as is convenient for vacuum sector perturbation theory. Therefore the first step in soliton perturbation theory is to convert plane wave normal ordering to normal mode normal ordering. The goal of the present note is to describe how this can be done in the case of a scalar field theory in 1+1 dimensions. The fact that the theory is scalar and only in 1+1 dimensions does not appear to play a central role in our analysis, and so we expect that the approach in this paper can be trivially generalized to more complicated theories in more dimensions.
We find that the problem of converting plane wave normal ordering into normal mode normal ordering can be achieved in two steps. First, as will be described in Sect. 3, we show that plane wave normal ordered products of the form : φ n (x) : a can be decomposed into sums of products of factors of the form : φ is a decomposition into different kinds of normal modes, such as even and odd breather modes. Next, in Sect. 4, we find that these factors can each be converted according to the Wick formula where the contraction, except for the case of zero modes, is schematically with ω k and ω p the energy of a normal mode and plane wave respectively. Combining the decomposition (1.1) with the Wick formula (1.2) for each component, one arrives at our Wick's theorem where I(x) can be found using the identity [9] (1.5) in terms of continuum normal modes g k (x) and breathers g i (x) together with the condition that it vanish at as |x| tends to infinity. We begin in Sect. 2 with a review of our formalism.

The setup
In this section we will review the one loop description of kinks developed in Refs. [10][11][12] using the formalism developed in Refs. [6,13,14], which has the advantage that it resolves the ambiguity noted in Ref. [15]. The key elements of our notation are summarized in Table 1.
For concreteness, we consider a theory of a real scalar field φ(x) and its canonical momentum π(x) in 1+1 dimensions, described by a Hamiltonian where M has dimensions of mass and g has dimensions of action −1/2 . The perturbative expansion will be an expansion in g 2h and we will seth = 1. The plane-wave normalordering :: a will be defined momentarily.
We assume that the potential V has degenerate minima so that the classical equations of motion admit a timeindependent kink solution φ(x, t) = f (x). (2.2) In the Schrodinger picture of the quantum theory, the translation operator satisfies the identity [13] : for any functional F and maps the vacuum sector to the kink sector. For example, the kink ground state may be written where | is the free scalar vacuum state and O may be calculated in perturbation theory. As |K is a Hamiltonian eigenstate, O| is an eigenstate of its similarity transform where Q 0 is the classical kink mass and H I consists of higher order terms in the g expansion. Note that g f (x) is dimensionless and so contains no powers ofh and so no powers of g. We assume that the field has been shifted so that φ = 0 is a global minimum of the potential V . Defining one easily evaluates It is a mess. The perturbative ground state, which is annihilated by a p , is not an eigenstate. As our shift function f (x) has explicit x dependence, this expression does not even have manifest translation invariance. This motivates us to search for a new basis of the operator algebra in terms of which This is a wave equation for a particle in a potential and its solutions are the normal modes of the field theory in the kink background. It generally has bound state and continuum solutions. We will refer to even and odd bound state solutions as g B E,i (x) and g B O,i (x) respectively, where the index i runs over distinct solutions if there is more than one. The corresponding frequencies, defined using (2.9), will be denoted ω B E,i and ω B O,i . There will always be an even bound state solution corresponding to the translation symmetry, which we call As it corresponds to a symmetry, it is a zero mode ω B = 0. The other bound state solutions correspond to breather modes. Let n e and n o be the number of even and odd breather modes. We will name the continuum states g k (x) where k is defined by ω 2 k = k 2 +m 2 and the sign of k is fixed by demanding that asymptotically it becomes the corresponding plane wave. All of these solutions are clearly mutually orthogonal and we normalize them such that (2.11) We also impose (2.14) where I runs over all n e + n o + 1 bound state field labels {B, {B E, i}, {B O, i}} and the symbol (−1) I is equal to −1 for odd breathers and otherwise is 1. This entire paper will be in the Schrodinger picture, and so the field φ(x) and π(x) will be independent of time and of any choice of Hamiltonian, and of any decomposition of the Hamiltonian into a free and imaginary part. They satisfy canonical commutation relations exactly, in the full interacting theory. 2 We define the linear combinations From the canonical algebra satisfied by φ(x) and π(x) one easily finds the algebra satisfied by their components with all other commutators within each decomposition vanishing.
Using the completeness relations (2.14) one may invert (2.15) to arrive at two distinct expansions of the same fields. The first is in terms of plane waves Here the the mass m is defined by The other decomposition is in terms of normal modes 3 For each g here there is a corresponding ω, defined by Eq. (2.9). In the case of the continuum modes in φ C , these are equal to It may appear surprising that m, which was defined in (2.18) using the other expansion, appears here. However (2.20) is merely our choice of parametrization of |k|, while its sign is fixed by (2.12). The only parametrization-independent content of (2.20) is that the minimum ω k for a continuum mode is ω 0 = m which arises because this delocalized mode has only measure zero support on the kink, and so its mass agrees with that of the lowest mass vacuum sector mode, m. Normal ordering may be defined with respect to either decomposition. Plane wave normal ordering :: a places all a † to the left of each a. Normal mode normal ordering :: b places all b † and φ 0 on the left of all b and π 0 .
Finally one may simplify H 2 where Q 1 is the one-loop correction to the kink energy. One recognizes this system as the sum of a free quantum mechanical particle with position φ 0 and momentum π 0 plus an infinite set of quantum harmonic oscillators. The one-loop vacuum therefore is annihilated by π 0 and also by all operators b We can now describe the entire Fock space of H 2 . It consists of the vacuum, plus an arbitrary center of mass momentum proportional to π 0 , where relativistic corrections appear only at higher orders in pertubation theory [10]. Also there is an infinite tower of oscillator excitations, created by b † k with mass ω k . Equation (2.22) implies that normal mode normal ordered operators : A : b , with vanishing c-number component, have vanishing expectation values at one loop This is one motivation for considering normal mode normal ordering. Another is that it allows an efficient computation of states and energies beyond one loop [9].

Factorization
The Hamiltonian H is plane wave normal ordered and a similarity transform by D f preserves the normal ordering [13]. Therefore the Hamiltonian H is also plane wave normal ordered. However for several applications, normal mode normal ordering is most efficient. In this paper we will study how to relate the two. The Hamiltonian H , at nth order, for n > 2 is where V (n) is the nth functional derivative of the potential V with respect to its argument. To calculate the soliton spectrum and energy corrections in perturbation theory, beginning with O 1 | , it is easiest to normal mode normal order the Hamiltonian. The plane wave normal ordering is defined in terms of a † and a and so, to evaluate these terms, we must use the plane wave expansion (2.17) To rewrite this in terms of normal mode operators, one need only insert (2.19) into the inverse of (2.17) to obtain the Bogoliubov transformations The key simplification comes from the fact that the modes from distinct oscillators commute with each other and they all commute with the zero modes. Thus, after inserting (3.3) into (3.2), one can separate the modes of each oscillator and the zero modes Note that in a local Hamiltonian, normal-ordered products appear in the combination (3.2) where this product is integrated over a kernel which is symmetric with respect to permutations of the p i . Thus only the symmetric part of the product contributes to the Hamiltonian. This depends on the subsets J M only via their cardinalities j M = |J M | which sum to n S : where S symmetrizes all values of p i . Where the letter M appears in the limits of the sum, it is understood that we have numbered the n o + n e + 2 values of M from 1 to n o + n e + 2.
The ordering chosen does not matter. For example, (3.5) becomes where the three terms correspond to To avoid clutter, below the operator S will not be written explicitly, but we will write in the text when we symmetrize. If we decompose φ(x) similarly to the plane wave operators then we can use (3.6) to decompose Note that the symmetrization is automatic here because of the symmetric kernel of the p integration in (3.2). The normal ordering on the right hand side of (3.4) is defined to be whatever one obtains from (3.2) when all of the different oscillators are separated. This is well-defined. But is it a normal ordering?

The problem
To simplify this question, let us restrict our attention momentarily to the case n e = n o = 0, as in the Sine-Gordon theory.
The generalization to other values is trivial. Clearly, whatever :: a on the a M means, it is linear since the factorization above can be performed separately for each summand. So consider one summand in (3.5) The simplest guess would be that :: a places the a † B on the left, and so the answer could be a Similarly, in the case of continuum modes (3.12) Thus the action of :: a on a † M, p and a M,− p is more complicated than simply putting all a † M on the left, since their order matters. This was not the case with the undecomposed plane wave oscillator modes because where we used the completeness relations (2.14) and the product of zero and a delta function vanishes at p 1 = p 2 because this is the commutator of an operator with itself.

Conclusion:
One may freely interchange the undecomposed plane wave mode operators a † and also a inside of :: a , for example in Eq. (3.2). However, this shuffling fixes the order of the components a † M and also a M in (3.4). In particular, the ordering of the components must be the same for all M, as this ordering is that chosen for the undecomposed operators.
Note that the symmetrized commutators vanish, and so this problem does not arise in the symmetrized products relevant to the computations of products of fields at the same point, as appear for example in the Hamiltonian.

A practical convention
In the previous section we learned that we need to make a choice. We need to choose the ordering of the a † p i and also of the a − p i in (3.2). This choice does not affect our answer but it fixes the orderings of each component in (3.4). In this section we will choose an ordering which will facilitate the computations in the next section.
Let us define the shorthand (3.14) We choose the ordering defined by We remind the reader that the value of N k does not depend on this choice of ordering, as all a † commute with each other as do all a. However it does affect the definition of the normal ordering of the components. Our strategy will be the following. First we will guess a formula for the normal ordering of the components Then we will show that, using the factorization formula (3.4) the guess yields the correct value of N k . Recall that our definition of :: a on components is that it satisfies (3.4) and so once we have shown this, we will have verified that our guess indeed satisfies the definition and so corresponds to a valid convention.
Our guess is The factorization formula (3.4) in the case n e = n o = 0 is Here we have adopted the shorthand p S for the ordered set of all p j with j ∈ S. The ordering is just the ascending order, since that appeared on the left hand side of the equation. We need to show that our guess (3.17) satisfies (3.18).
Our proof will be by induction. The base case, k = 0 is trivial as the only term in the sum is J = ∅ and so (3.18) becomes 1 = 1. Next assume that (3.18) is satisfied for some value of k and definê (3.19) where the right hand side is defined using (3.17). We need to prove thatN = N to complete the induction. Each J either does or does not contain the element {k + 1} and so we may respectively divide the sum in two parts, redefining the dummy set J in the first sum by removing {k + 1} completing the induction. In summary, we have shown that if we adopt the definition (3.17) for the plane wave normal ordering of component fields a † M and a M , then the factorization formula (3.18) is satisfied and so these components N M k can be assembled to determine the plane wave normal ordered product N k of the undecomposed operators. Although our proof was for the case with no breather modes n e = n o = 0, the equation  where N B n was defined in (3.16). Then using (3.17) we can find the next product Dividing through by the product on the left one finds Finally matching terms with the same power of φ 0 we arrive at the recursion relation which, together with the initial condition α 0m = δ m,0 fixes all of the coefficients α. The recursion relation (4.5) has a simple interpretation in terms of a Wick's theorem. m is the number of contractions. The (n + 1)st operator may either not contract, leading to the first term on the right hand side, or else it may contract. If it does contract, since there are m contractions in all, the first n operators have m − 1 contractions. Therefore the n + 1st operator may contract with any one of the n − 2(m − 1) uncontracted operators, yielding the factor of n − 2m + 2 in the second term. Each contraction yields a factor of −1/(2ω p n+1 ). Note that this contraction factor is not symmetric with respect to a permutation of the p i , since it depends only on the p i with the highest value of i among the two contracted operators, which is p n+1 .

Solving the recursion formula
Recall that to compute the Hamiltonian we only need the symmetrized N n . In this case the choice of ω p i is irrelevant, it is only important that no N have two ω p i with the same i. Said differently, adding an antisymmetric piece to N will not change the symmetrized N and so will not change H . We can thus shift N to be of the form where a is a pure number which simply counts the number of ways to make m contractions. a satisfies the recursion relation a n+1,m = a nm + (n − 2m + 2)a n,m−1 . (4.7) As the contractions are interchangeable, a nm contains a factor of 1/m!. This is multiplied by the number of choices for the jth contraction, which is (n − 2 j + 2)(n − 2 j + 1)/2, for each j from 1 to m. In all one finds In the case with no breathers, the decomposition of the fields (3.9) becomes Assembling the results above, we have evaluated the first factor in (4.9) where we have introduced the contraction factor (4.11)

Odd breathers
Similarly to the plane wave ordered products N n ( p) we will define the normal mode ordered products Our goal in this section is to learn how to expand N n ( p) in terms of B B O n (k). Using the identity one readily derives the anticommutator (4.14) and the commutator which will be useful momentarily.
Proceeding as for the zero mode, we define coefficients α nm by (4.16) Then using (3.17) Matching coefficients we obtain the recursion relation So far we have not used symmetrization, and so our recursion relation may be applied to computing any n-point function. Again, for calculating n-point functions at the same point, as in our interaction terms, we may shift N B O by an operator which vanishes when symmetrized (4.20) Note that the product on the right can be rewritten so that it resembles the contraction terms in (4.6). Proceeding as above, the recursion formula satisfied by the a nm is again (4.7) and so the a nm are given by (4.8).
where we have defined the contraction factor The contraction factor is similar to I B (x) except that the contraction contains two terms 1/(2ω B O ) and 1/(2ω p ) with a relative sign. These are respectively the contraction arising from the normal mode normal ordering and the plane wave normal ordering. In the case of I B (x) the normal mode normal ordering was fundamentally different, as it was a rule for the placement of the canonical variables φ 0 and π 0 and not for the oscillator modes. The occurrence of a difference of contractions in I B O is reminiscent of the general contraction defined in Ref. [18].

Even breathers
The normal ordering of even breathers is identical to that of odd breathers except for a few sign differences. Defining and using the identity Then defining (4.27) The same computation as in the odd case yields the recursion relation (4.28) Comparing (4.19) and (4.28) one sees that the contractions of even and odd breathers differ by an overall sign.
In the symmetric case one may shift N B E to where a bm again satisfies (4.7) and so we conclude that where we have defined the contraction factor The relative sign in the recursion relation has indeed translated into a relative sign in the contraction factor with respect to I B O . As g B O (x) is imaginary and g B E (x) is real due to our convention (2.12), the relative sign may be absorbed by taking the complex conjugate of g(x) in the definition of I(x).
We will now see in the continuum case that this definition arises quite naturally.

Continuum modes
Define Using the identity and similarly the anticommutator We will need the integrals of these identities, where the integral over k is performed using the Dirac delta function for arbitrary matrices α nm . We will define the matrices α nm by (4.38) Then (3.17) implies (4.39) Summarizing, we find wherek j indicates that k j is omitted. Matching yields the recursion relation (4.41) Symmetrizing we may write where again a nm satisfies (4.7) and so is given by (4.8). We therefore conclude While the algebra leading up to our result seemed more complicated than in the case of the bound states, our final result is essentially the same. The only difference is that I C is integrated over normal modes k. However, even in the case of breathers, there will be a sum over breather modes i, and so this distinction is superficial.

Wick's Theorem
Finally we may use our decomposition (3.9) to reassemble : φ j (x) : a by inserting (4.10), (4.22), (4.30) and (4.43) The contraction factor I(x) can be efficiently evaluated using the completeness relations (2.14) It is finite. The first term is x-independent, and integrates to infinity. Now we take the derivative with respect to x, and so the first term vanishes Notice that the relative sign between the even and odd breathers has disappeared, because g i (x) 2 is positive in the even case and negative in the odd case. This, together with the boundary condition that I(x) vanishes as x goes to infinity, completely determines I(x).

Normal ordering the Hamiltonian
In the Sine-Gordon theory the interaction Hamiltonian density in H is [19] The contribution arising from bound states is Using (4.10) the plane wave normal ordering may be evaluated explicitly To simplify this sum, we will include the terms at n = 0 and n = 1, which are present in the Hamiltonian although they are not the only terms at their orders. These terms only affect the noninteracting part of the Hamiltonian, which is known to be the Poschl-Teller Hamiltonian. So we definẽ The n = 0 and n = 1 terms that we have added arẽ Similarly, including the n = 0 term, Substituting this back into (5.2) and letting the tilde remind the reader that the noninteracting (quadratic in φ and below) terms are incorrect, we find As the terms in (5.9) are at most O(λ 0 ), and H I anyway does not include all terms at this order, we will not be interested in these terms. The terms at these leading orders are given in this example in Ref. [19], and more generally are of the form (2.21). In particular they include no interactions and the tadpole vanishes in the full expression. Equation (5.8) has a straightforward interpretation. The combination φ B (x) + f (x) is just the D f translated field, brutally truncated to the zero mode part. The prefactor and the cosine term are thus just the original Sine-Gordon action, translated and truncated. However we see that the plane wave normal ordering is now gone, indeed it was our goal to eliminate it, and instead there is an exponential of a contraction term. Thus plane wave normal ordering is equivalent to multiplication by the exponent of the bound state contraction. Of course only the bound state contraction appeared because we have truncated our Hamiltonian by only considering the bound component of the field. Our result is trivially normal mode normal ordered as it only involves the operator φ 0 .
More generally we may expect the exponential to include the sum of the contractions of the various normal modes The appearance of contractions in an exponential in (5.10) is also similar to the generalized Wick's theorem postulated in Ref. [18]. It would be useful to understand this connection more precisely, as the generalized Wick's theorem may provide a simple extension of our results to more complicated and interesting models.

Application: time independent perturbation theory
We have already noted that Eq. (5.10) is only correct beyond O(λ 0 ), as we did not include all terms at lower orders. Recalling that f is of order λ −1/2 , the cos term may be expanded as in (5.1) In particular the O λ 1/2 terms in H I are (5.12) This was reported already in Ref. [9], but now using Wick's theorem a calculation that originally required several days has been reduced to less than an hour. Using (5.11) we can easily go to higher orders. For example at O(λ) we find that H I contains To apply this formula, first expand the soliton ground state |0 n (5.14) where |0 0 = O 1 | is the one-loop state described above and each successive term comes with another power of λ 1/2 . Then, similarly expanding the Hamiltonian the Schrodinger equation can be expanded order by order, beginning with The first equation was solved in Ref. [16] and the second in [9]. However Q 2 and |0 2 are still unknown. We can see from the third equation that Q 2 receives three contributions, one from each term in the decomposition of |0 0 . From (5.13) we see that the |0 0 contribution to the two-loop correction to the soliton mass Q 2 is The |0 1 contribution Q (1) 2 can be found by acting (5.12) on the expression for |0 1 in [9]. To compute the |0 2 contribution Q (2) 2 one needs to use the translation invariance of the ground state as described in [9]. The full calculation appears in Ref. [20].

Another application: correlation functions
The unsymmetrized Wick's theorem can be used to calculate Schrodinger picture correlation functions of plane wave normal ordered operators, should one be interested in such objects and φ C (x) commute, this product can be factorized into a sum over all subsets S of [1, n] A complete set of states in the one-soliton sector is generated by functions ψ(φ 0 ) times products of b † k acting on |0 0 . Inserting the identity written in terms of such a complete set of states in between φ C (x n ) and φ B (x 1 ), one sees that all terms with b † operators vanish because the right matrix element contains no other b † operators. Similarly, if one considers an orthogonal basis of functions of φ 0 , then only the trivial function 1 will contribute as for all other functions the left matrix element will vanish by orthogonality, as none of the operators on the left contain φ 0 or π 0 and they all commute with φ 0 . Thus we have argued where for any set S = i 1 , .., i n we have adopted the shorthand that x S is x i 1 . . . x i n . These products of fields can be rewritten in terms of our symbols N as in (4.10) but with no symmetrization As all matrix elements of normal mode normal ordered operators vanish, only the completely contracted operators (n = 2m) contribute. These depend on the coefficients α n,n/2 . Immediately we see that the n-point functions vanish unless n is even. Let us begin with the bound state components. Solving the recursion relation (4.5) in a few cases we find α n,0 = 1, The completely contracted terms can be written immediately in the first few cases It is remarkably asymmetrical in the x i . This is to be expected, as we made an arbitrary choice in our convention for the bound state normal ordering. The correlator F B depends on that choice. Only the product F is independent. Now let us turn our attention to the continuum modes φ C . Solving the recursion relation (4.41) in the first few cases we find α k 1 ...k n n,0 = n j=1g k j ( p j ), (5.24) The first completely contracted term and so, using the fact that for the empty set S, F M (x S ) = 1 we find This is looks like the completeness relation (2.14) but is different because of the hats. As we are working in the Schrodinger picture, all operators are time-independent and so these are equal time correlators. When x 1 = x 2 the points are spacelike separated. Nonetheless this trivial calculation shows that the correlator does not vanish for spacelike separated x 1 and x 2 . Also it is not translation invariant. Neither of these properties should be a surprise. This is not a correlation function in the vacuum, but in the 1-loop kink state |0 1 which is not translation invariant. It is, to one-loop, a Hamiltonian eigenstate and this condition implies that it has nonlocal correlations, just as the phase in a superconductor varies slowly over macroscopic distances or as parton wave functions are correlated in a hadron.
What about causality? Causality demands that for spacelike separated x 1 and x 2 . This expression is not normal ordered and one needs to be careful when normal ordering a commutator as commutation and normal ordering do not commute If we start by using the equal time canonical commutation relation [φ(x 1 ), φ(x 2 )] = 0 then (5.28) is trivially satisfied and we have not succeeded in testing our correlation function. So instead rewrite (5.28) as Now we can normal order both sides. So long as we do not to use the commutation relations, equality before normal ordering will imply equality after normal ordering and thus we will obtain a weaker but still necessary condition for causality Is this satisfied by (5.27)?
We used the completeness relation to arrive at the second line.
The last line has two terms. The first is symmetric under x 1 ↔ x 2 because ω p is an even function of p and so the odd piece of the exponential does not contribute to the integral. To see that the second term is also symmetric, it suffices to relabel the dummy variable k → −k . Therefore the causality condition (5.31) is satisfied. The continuum contribution to the four-point function is The total four-point function is then given in terms of Eqs. (5.23), (5.26) and (5.33).

Remarks
We have found that plane wave normal ordering can be converted into normal mode normal ordering by following a simple rule, playing the role of Wick's theorem. After decomposing a product of n fields into products of j field components, where each component corresponds to a set of normal modes, the components can be decomposed by summing over all possible contractions. For each contraction one replaces the pair of field components with the difference between the inverse plane wave energy ω p and inverse normal mode energy, suitably normalized over the spectrum. Intuitively the first term arises from eliminating the plane wave normal ordering and the second from imposing the normal mode normal ordering.
Of course with no normal ordering at all, one expects divergences. However the difference between these two energies is, when suitably averaged, quite small and thus all expressions are finite given either normal ordering scheme. Once we go beyond scalar theories and 1+1 dimensions there will be other divergences which must be regularized and renormalized. In Ref. [9] the conversion between normal orderings was the most complicated part of the perturbation theory treatment of the one soliton sector. Now that we have treated this problem at all orders, and in a much more general class of theories, we expect that it will be easier to extend that calculation to two loops or beyond. In Ref. [20], we have used the Wick's theorem here to find the two-loop ground state and mass of a scalar kink in a theory with an arbitrary potential. Previously the scalar kink mass was only found in the Sine-Gordon model [21][22][23], and our results agree in that case. However it is still not obvious that the solution to the zero mode problem in Ref. [9] also solves the problem at higher loops. If it does not, then it may be necessary to use other formalisms such as that of [7] and [8].
To go beyond perturbation theory, we will eventually need supersymmetry. In this context, coherent states have been constructed in Refs. [24,25]. This will require a fermionic generalization of the Wick's theorem found here. Perhaps the generalized Wick's theorem of Ref. [18] can provide an efficient derivation.
Another application of our Wick's theorem is the calculation of plane wave normal ordered correlation functions in soliton states. So far our calculations have been in the Schrodinger picture. As a result we can only evaluate equal time correlators. However we suspect that a generalization to the interaction picture would be straightforward, as interaction picture operators evolve according to the free Hamiltonian. This would allow us to compute correlators of operators at arbitrary time. We note that there are two interaction pictures, corresponding to the free Hamiltonians H and H which describe perturbations in the zero-soliton and onesoliton sectors. For the calculation of one-soliton state correlation functions, one would use the latter as the one-loop one-soliton ground state is an eigenstate of the noninteracting part. As the one-soliton ground state is an eigenstate of H , it is even plausible that an LSZ reduction formula exists in the one-soliton sector. This would be the starting point for a treatment of scattering in this sector.
The recent discovery of spectral walls [26] caused by transitions between breather and continuum states has rekindled interest in kink scattering [27,28]. The treatment of this phenomenon has so far been largely classical. While the current methodology is most straightforwardly applied to the one kink sector, it could nonetheless allow an understanding of the role played by breathers in fully quantum scattering. In particular the scattering of a kink with a plane wave or wave packet could be treated in the one kink sector.