Relative-energy dependence of the Bethe–Salpeter amplitude in two-dimensional massless QED

The exact qq¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q\bar{q}$$\end{document} Bethe–Salpeter bound state amplitude is investigated in the space of relative energy E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} for a fixed value of the relative position. By means of an approximate analysis it is shown to possess singularities in E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E$$\end{document} whenever one of the quarks reaches the energy threshold for the creation of a certain number of Schwinger bosons. These results are confirmed by plots of numerical results.


Introduction
One of the essential and unsolved problems in quantum field theory (QFT) is the correct description of bound states in a truly relativistic way, with retardation effects taken into account. On one hand it belongs to the less understood aspects of QFT, and on the other it is one of the most important issues, necessary for the explanation of the whole spectrum of hadrons. Over 60 years ago the equation for the bound state amplitude-the so-called Bethe-Salpeter (B-S) equationwas proposed [1,2] but its solutions are still lacking even in model studies.
Solving the B-S equation is not an easy task: it is an integral equation in several dimensions but, what is worse, it requires from the very beginning the knowledge of the explicit form of the nonperturbative propagators for bound particles, as well as their interaction kernel. Because none of these quantities is known in realistic field theory or even expected to be known, one is inevitably led to strong simplifications as for instance that consisting in replacing full propagators with free ones and using the so-called ladder approximation. But even then the problem is still very complicated and has found its solution only in few model calculations [3][4][5]. In hadronic physics a step forward has been done in some numerical studies exploiting the truncated (i.e. again approximated) QCD [6]. a e-mail: t.radozycki@uksw.edu.pl A novel and troublesome feature of the B-S equation is the appearance of the additional variable, absent in nonrelativistic studies: the so-called relative time or equivalently the relative energy. There have been undertaken many attempts to simplify the equation by neglecting the relative time dependence of the amplitude [7][8][9][10][11][12][13][14]. Such approaches, however, lead to oversimplification: by imposing the instantaneity condition they throw aside the most important feature, which differentiates the relativistic description from the nonrelativistic one. In such an approach the most interesting part of the information is then lost because of the approximation itself.
The analysis of that nontrivial aspect of a bound system has been undertaken in a couple of works [15][16][17][18][19][20][21] both in euclidean [22][23][24][25][26][27][28] and in Minkowski space [28][29][30][31][32][33][34][35] (where it is more difficult), but a complete description of any B-S amplitude in its full complexity is missing, with one exception, to be mentioned below. That has been mainly due to the lack of an appropriate model, which on one hand would constitute a nontrivial QFT and on the other would allow for analytic and exact solutions. However, a couple of years ago the B-S amplitude in the full analytic form was found in position space in massless quantum electrodynamics in two dimensions (in the so-called Schwinger model (SM) [36]) [37,38], allowing ipso facto for a more detailed investigation of the role played by relative variables in a nontrivial field theory.
The model in question is defined through the Lagrangian density in two-dimensional (2D): where g is the coupling constant and ξ the gauge fixing parameter (later set to be equal to infinity, which corresponds to the Landau gauge). The Dirac gamma matrices may be chosen as 2D objects, for instance in the following form: and for the metric tensor we have g 00 = −g 11 = 1. In the last 50 years this model has become a perfect testing laboratory for various nonperturbative properties of QFT. Particular attention deserve, in this context, confinement, topological sectors, instantons, and condensates. Among other interesting features one should mention the existence of an anomaly and nonzero mass generation of the gauge boson without the necessity of introducing any auxiliary Higgs field.
In this work we concentrate on another remarkable property of this model, spoken of above. Thanks to the analytic determination of the nonperturbative four-point Green's function [39][40][41] it allows for the investigation of the formation of bound states. It has been found that the SM model possesses a state, which might be called a 'meson', constituting a bound system of a quark and an antiquark (the fundamental fermion of the theory is often called a quark, due to the similarity between SM and QCD in some respects, as for instance confinement) and which is also known as the Schwinger boson of mass μ = g/ √ π . In our former papers we concentrated on this issue. Firstly, we analyzed the structure of the full four-point (two-fermion) Green's function and found the exact and explicit form of the B-S bound state amplitude in all contributing instanton sectors [37,38] (i.e. for k = 0, ±1). This result seems to be of a certain significance, since, to the best of our knowledge, it is a unique exact B-S function, which is known exactly in any nontrivial field theory. This result was obtained not by solving the B-S equation, which we are even unable to write down, but by analyzing the residue in the pole corresponding to the Schwinger boson (cf. [42,43]).
Secondly, a subsequent paper [44] was devoted to the investigation of the properties of the previously found B-S function in the space of relative variables. The present work continues this research concentrating on the relative-energy variable.
The B-S function, which was found from the t-channel of the two-fermion Green's in the residue of the pole corresponding to P 2 = μ 2 (where P is the center-of-mass twomomentum), i.e. to the Schwinger boson, has the form where x = [t, r ] denotes the relative quark-antiquark coordinates and the indices 0 and 1 refer to the k = 0 and k = ±1 instanton sectors. Higher sectors do not contribute. The Green's function in question is a vacuum (strictly speak-ing θ -vacuum) expectation value of four fermion fields and the residue factorizes into two bilinear functions. It is well known, however, that to change the topological index of the vacuum by 1 an operator which is quadratic in the fermion fields is needed [45,46]. This is connected with the masslessness of fermions in the theory and the absence of the chirally noninvariant term m in the Hamiltonian, which results in the suppression of tunneling between topological vacua. Consequently only chirality changing operators can have nonvanishing matrix elements between various topological vacua, with the 2:1 correspondence between chirality of the operator and instantonic number.
The following expressions for P have been found in the above mentioned works: (1) where θ is the vacuum parameter. S 0 (x) here denotes the free fermion propagator: (6) and the function β, appearing in many quantities of the Schwinger model, is defined as The symbol γ E is here the Euler constant and H (1) 0 and K 0 are Hankel function of the first kind and the Basset function, respectively.
The present paper is organized as follows. In the next section we consider (0) P , the exact B-S amplitude in the instanton sector, k = 0. We investigate its relative-energy dependence for fixed values of relative position in the centerof-mass frame. Using some approximations we identify the location and character of singularities in relative energy. This analysis has a qualitative rather than a quantitative character. The obtained results are confirmed by the appropriate plots. In Sect. 3 we concentrate on the B-S amplitude in the sector k = 1. We carry out the same program as in Sect. 2 and find similar singularities for the function (1) P . The last section contains a summary of the results obtained and some conclusions.

The B-S function in the noninstantonic sector
The relative-energy B-S function in the instanton sector k = 0 is defined through the Fourier transform over relative time: Fourier-transformed quantities will be, throughout this paper, identified only by their arguments, without the change of the appropriate symbols. We hope it will not be confusing. The one half in the exponent comes from the definition of the center-of-mass variables and relative ones, we use for the two-body system and hence With this convention the Fourier exponent in question becomes where P = [P 0 , P 1 ] = p 1 + p 2 is the total (centerof-mass) two-momentum satisfying P 2 = μ 2 , and Q = [E, q 1 ] = p 1 − p 2 is the relative one. The inverse relations are The quantities C 0 P (E, r ) and C 1 P (E, r ) may be given the following form: where we introduced the notation The coefficient functions A and B after some rearrangements are expressed through the Cauchy principal value integrals as follows: From the asymptotics of the Hankel function of the first kind it may easily be seen that for large t the factor e −ig 2 β(t,r ) behaves as and, by virtue of the Dirichlet rule, the integral for A owes its convergence at infinity to the interplay between the oscillatory functions sin ωt and H (1) 0 in (7). This is an important observation for our conclusions. The convergence of B, which has one power of t more in the denominator, is absolute, but it may be improved by an oscillatory factor too. However, due to the complicated form of the function β(t, r ) these integrals cannot be found explicitly and one is doomed to use approximations or perform a numerical analysis. It is a common feature of the Schwinger model that analytic solutions can be calculated only in coordinate space. It refers to various physical quantities, an exception being purely bosonic Green's functions.
It is clear that the function A(ω, r ) is odd in the argument ω and even in relative position r , and B(ω, r ) conversely. Consequently C 0,1 P (E, r ) have the following symmetry properties: In the center-of-mass frame, where P 1 = 0 and P 0 = μ, these quantities satisfy Let us now try to identify the possible singularities of A(ω, r ) in ω for fixed r . According to what was said above, they arise when the Dirichlet rule fails (or rather may not be applied due to the cancelation among oscillating factors) and therefore it is sufficient to consider the integral from a certain point λ|r | to infinity (with λ > 1), which will be called A(ω, r ). The last term in (16) is inessential for this analysis, since it does not contribute to any singularities. Using (7), we havẽ Since for large values of t the Hankel function behaves as [47] H (1) each following term of the above expansion as well as that of the exponent function in the integrand of (22) has better convergence properties, contributing thereby less to the singularities in ω. Therefore the strongest ones, which can eventually be clearly identified on the plots, come only from the first few terms: Below we analyze them starting from the most singular one: For large values of λ one can neglect r 2 in the denominator, and the value of the integral can be found to bẽ where S(x) is the Fresnel integral: We are interested in the behavior at ω ≈ 0 (it is the only possible singular point) with large but fixed λ, which corresponds to a small value of σ . In such a case we can use the expansion and the first two terms of (25) becomẽ The quantityÃ a is then seen to possess a strong singularity at ω = 0, which should be represented with an infinite peak on the plot of (0) The subsequent term of (22) after some rearrangements has the form Again r 2 may be omitted for large λ and one can obtain for A b the approximated result: where Ei is the exponential integral function. It has the following expansion for an imaginary argument: which shows thatÃ b (ω, r ) contains the logarithmic singularities at ω = ±μ. They should again manifest themselves as infinite peaks, although less pronounced than those of (28). The next term may be given the form One can show in the same way as above that it has branch points of the kind √ ω ± 2μ, but the integral is convergent even for ω = ±2μ.
The quantity B is smoother, but it also contributes to singularities. We defineB in an analogous way toÃ and find that the most singular term is This integral is convergent in an absolute way and it may be shown to contain a singularity of the kind √ |ω|, which might eventually appear in the form of slight cusps, but they will be hidden by the peaks ofÃ a in the same place. The following terms inB(ω, r ) need not be considered. Because (0) P is a combination of A(ω + , r ) and A(ω − , r ) via the Eqs. (8), (13), and (14), we expect the manifestations of the following well-defined singularities in the form of peaks in the B-S amplitude: The singularities coming from (32) occurring at ω ± = 2μ and ω ± = −2μ, which correspond to E = −5μ, −3μ, 3μ, 5μ are expected to be less visible. What is important, the positions of all the singularities are not affected by the value of the relative coordinate r . These analytical although approximate results are confirmed below by the appropriate numerical graphs.
Since (0) P is a 2 × 2 complex matrix, it cannot be directly represented on the plot. To show the relative-energy dependence of the B-S amplitude, we therefore define a scalar, similarly to what was done in [44]: It is a gauge invariant, dimensionless quantity and can serve as a measure for the strength of the B-S amplitude. It may be simply expressed in terms of A and B, as where bars over symbols refer to complex conjugation.
To investigate the relative-energy dependence of the B-S function, the most natural thing to do is to perform the plot of | (0) P (E, r )| in the CM frame, thereby putting P 1 = 0. In this frame it takes the form (36) and becomes an even function in both relative energy and relative position. In a boosted frame the plotted function is no longer even, but the positions of the singularities in relative energy remain invariant. Figure 1 shows the dependence of (36) on the relative energy for certain fixed relative positions. The main noticeable feature is the presence of strong peaks for the relative energy constituting an odd multiple of the Schwinger mass μ. This is what we expected from our approximations.
Some explanations require the behavior at E = ±3μ, where, according to our former analysis, the function should display a logarithmic nature. This ought to lead to infinite peaks, and this really appears to be the case, when investigated carefully. It may be, however, estimated that they are extremely narrow. For instance at the value of 1 the width becomes almost 10 times narrower than the thickness of the plotting line, and therefore it cannot be properly represented on the drawing. The singularities at E = ±5μ can also be identified, but it requires very high resolution plots.
The appearance of the singularities at relative energies corresponding to an odd number of Schwinger boson masses μ may be understood from the simple picture. Let us consider a quark-antiquark pair of energies E 1 and E 2 which form a bound state, i.e. a meson of the CM energy equal to μ. If, due to some fluctuation, one quark acquires the energy sufficient to create a meson, i.e. E 1 = μ, becoming of course a deeply off-shell particle (if one at all can talk about on-shell particles in the bound state) then we must have E 2 = 0 and the relative Now assume that a quark gained an energy sufficient to create n bosons. Then E 1 = nμ and E 2 = (1 − n)μ in order that E 1 + E 2 = μ (in this rough image). In such a case for the relative energy we obtain E = E 1 − E 2 = (2n + 1)μ. Such fluctuations are, however, much less probable and hence higher singularities are softer.
The other observation, which can be done by analyzing the E and r dependence of | (0) P | is that for small relative energies the amplitude decreases with increasing r . However, for E exceeding the first threshold, i.e. for E > μ, the described behavior becomes opposite: the amplitude grows with relative distance. This observation has yet some limitations: we were not able to verify it for very big values of r since then the oscillatory integrals (16) and (17) become extremely slowly convergent. As we will see in the following section, this behavior will not be the same for the k = ±1 instanton sector.

Instantonic contribution to the B-S function
The relative-energy B-S function in the instanton sector k = ±1 has the form The tensor structure of (1) P is simpler than that of (0) P . It contains only one (scalar) coefficient D P (E, r ), which may be written in terms of a new function F(ω, r ), defined as in the form For large t the exponent under the integral has the following behavior: which again guarantees the convergence (via the Dirichlet rule) of the improper integral (38), except for certain particular values of ω corresponding to divergent singularities. The function F(ω, r ) is, in an obvious way, an even function in both arguments, which gives for D P : Consequently in the CM frame D P (E, r ) has the following symmetry properties: Let us now concentrate on the possible singularities of F(ω, r ) as the function of ω with fixed value of relative position r . Similarly to the quantity A(ω, r ) we focus on the behavior of the integrand function at infinity. We therefore definẽ assuming that the parameter λ is large. Using again the expansion (23) we find the less convergent term in the form For λ 1 one can omit r 2 as before, and the approximate value of the above integral is where C(x) is the second Fresnel integral, defined as For small values of σ it may be approximated, according to the well-known formula: and we find inF a (ω, r ) a similar singularity at ω = 0 as was identified inÃ a (ω, r ): For the next term of the expansion of (43) we obtaiñ Using (31), we see that (50) is logarithmically divergent as ω → ±μ. It is a result identical to the case ofÃ b (ω, r ). The last term we are interested in is Due to the sufficiently high power of t in the denominator, this integral is convergent even for ω = ±2μ. It may be shown to contain a slight singularity of the kind √ ω ± 2μ.
One can then summarize that the instanton sector k = ±1 reveals identical singularities at the same values of ω as for k = 0. It is also confirmed by the appropriate plots.
Similarly to (34) we define now the quantity, which represents the 'value' of the B-S function in this sector: It is independent both on the choice of gauge and on the vacuum θ parameter and may be given the form which in the CM frame reduces to the simple expression Similarly to | (0) P (E, r )| it becomes an even function of both relative position and relative energy. In Fig. 2 the plot of (54) is presented as a function of E for various fixed values of r .
The general relative-energy dependence reflects the principal features of the trivial topological sector. The peaks at E = (2n + 1)μ for small values of n are clearly visible. However, contrary to Fig. 1, there is no apparent enhancement with increasing r for relative energy exceeding the first threshold.
In Fig. 3 we plot for completeness the dependence on E of the full B-S function | P |, defined as

Summary
In summary, one should emphasize that the Schwinger model turns out to be an exceptional tool to investigate the properties of the bound states in QFT (apart from another nonperturbative aspects). The fact that the exact form of the B-S amplitude is explicitly known, constitutes a unique opportunity to analyze its dependence on relative variables and particularly the role played by the relative energy. The results obtained in the present work show that the 'strength' of the B-S function defined by Eqs. (34) and (52) exhibits singularities for odd number of Schwinger boson mass μ. The appearance of this threshold structure remains in some correspondence with the results of the approximated scalar-scalar model [15]. This effect finds its justification in a simple picture, in which we attribute energies E 1 and E 2 to both quarks and keep the total energy fixed, which in the center-of-mass frame is equal to μ. The location of singularities is identical in all instanton sectors. The other effect, present only in the instanton sector, k = 0, is the enhancement of the 'value' of the B-S function for medium values of r (in units of μ −1 ) for a relative energy exceeding the first threshold for production of Schwinger bosons, and opposite behavior for small r . The region of large relative positions was not investigated due to obstacles of numerical origin.
Our formulas show that, in a boosted frame, the distribution of singularities remains unaltered, although plots lose in this case their symmetric character with respect to the replacement E → −E.
It should be noted that the results of Sects. 2 and 3 stay in agreement with those of [44], where it was found that the singularities in the two-momentum space are located on hyperbolas [see Eq. (17)]. In the CM frame and in the notation of the present work, they have the form 1 4 (E + μ) 2 − (q 1 ) 2 = n 2 μ 2 , and may be rewritten as Now, in the Fourier integral over q 1 the singularities in E arise, according to the Landau procedure [42,48,49], when two singularities of the integrand functions merge. This happens for E = (2n − 1)μ or E = −(2n + 1)μ in the first case, and for E = −(2n−1)μ or E = (2n+1)μ in the second one. They all correspond to an odd number of Schwinger boson masses, as found in the present work. This is only a rough argumentation, since the explicit form on the B-S amplitude in two-momentum space is unknown.
The results obtained here may be of some interest for studies of bound states in more realistic models, both in particle and nuclear physics. It is believed that the retardation effects may play an essential role in true QFT bound states [50]. One should also mention that processes in which bound states appear in internal lines require Feynman integrations over relative coordinates or momenta and are strongly affected by the presence of singularities. The B-S function appears also in the matrix elements and scattering processes involving bound states [43,51] and its role is significant for properly describing the meson decay properties. Therefore any insight into its structure seems to deserve particular attention .
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3 / License Version CC BY 4.0.