1008 Page 4 of 18 Eur. Phys. J. C (2019) 79 :1008
Let us now present separately the parameterizations we
have used to describe the two partial waves of interest for
this work.
3.1 S0-wave parameterization
As explained above, our parameterizations will be consistent
with the dispersive data analysis of [28], which extends up
to 1.43 GeV. Above that we will only provide three phe-
nomenological fits to three sets of incompatible data, care-
fully matched to our parameterizations below. Let us discuss
both regions separately.
3.1.1 S0-wave parameterization below 1.4 GeV
The σ/f
0
(500) and f
0
(980) resonances dominate the behav-
ior of the S0 partial wave in this region. The somewhat con-
troversial f
0
(1370) couples very weakly to two pions and its
effect in this region can be treated as background. For our
purposes it is important to remark that the σ/f
0
(500) has an
associated pole very deep in the complex plane that produces
a wide structure increasing monotonously from threshold up
to roughly 900 MeV, reaching a phase-shift of 90
◦
around 800
MeV, as seen in Fig. 1. It is known [40] that the σ/f
0
(500)
pole can be generated in the S0 partial wave by a simple trun-
cated conformal expansion, that we will call t
0
0,conf
(s).How-
ever, above 900 MeV the f
0
(980) pole adds a further sharp
increase that makes the phase larger than 200
◦
right below
the K
¯
K threshold. The phase then keeps growing slower but
monotonously up to 2 GeV.
It is worth noticing that the interplay between the
σ/f
0
(500) and f
0
(980) poles produces a sharp dip in the
modulus of the amplitude and the elasticity right above the
K
¯
K threshold. In order to describe the f
0
(980) effects accu-
rately and consistently with the dispersive results, we will
factorize in the S matrix the f
0
(980) shape separately from
the conformal expansion that contains the σ/f
0
(500) pole.
In other words, we will use S
0
0
= S
0
f
0
S
0
0,conf
, where
S
0
0
= 1 + 2iσ(s)t
0
0
,
S
0
0,conf
= 1 + 2iσ(s)t
0
0,conf
,
S
0
f
0
= 1 + 2iσ(s)t
0
f
0
. (3)
This factorization ensures elastic unitarity for the S0 wave,
i.e., |S
0
0
|=1, when both the conformal and the f
0
(980)
contributions fulfill elastic unitarity independently, i.e.,
|S
0
0,conf
|=|S
0
f 0
|=1. This will be the case in the elastic
region below the
¯
KK threshold, s < 4m
2
K
.
For our purposes, we are interested in the amplitude
partial-wave
t
0
0
(s) = t
0
0,conf
(s) + t
0
f
0
(s) + 2i σ(s)t
0
0,conf
(s)t
0
f
0
(s). (4)
Now, the conformal factor of the partial wave is con-
structed by analogy to the elastic formulation in Eq. (1)
t
0
0,conf
(s) =
1
σ(s)
1
0
0
(s) − i
, s < 1.4GeV (5)
where, building on [40]
0
0
(s) =
√
s
2q(s)
m
2
π
s − z
2
0
/2
z
2
0
m
π
√
s
+
N
n=0
B
n
ω(s)
n
. (6)
Let us remark that
0
0
(s) = cot δ
0
0
(s) in the elastic region
s ≤ 4m
2
K
, implying |S
0
0,conf
|=1. The s −z
2
0
/2 denominator
provides the so-called Adler zero at s
Ad ler
= z
2
0
/2 required
by chiral symmetry [59]. For z
0
= m
π
one would recover
the Current-Algebra result, namely the leading order ChPT
value. However, for us it will be a free parameter, fundamen-
tal to describe the subthreshold region. As we will see, it
comes out from the fits consistent with the dispersive evalua-
tion, which in turn is consistent with higher order ChPT eval-
uations (see [28,60]). Note also that if we did not include the
∼ 1/
√
s term added to the conformal series then, for the val-
ues of B
n
obtained from the fit, we would find in the unphys-
ical region that
0
0
(s
Ad ler
) →+i ∞ and
0
0
(0) → i 0
+
,so
that necessarily
0
0
(s) = i somewhere between those two
subthreshold points. This would yield a spurious pole, i.e., a
ghost, in the partial wave. As explained in [40] these ghosts
are mostly harmless and have little relevance in the fit quality
and the pole positions, but as a matter of principle it is bet-
ter to remove them. Actually, as explained there, adding the
1/
√
s term to the conformal series simply shifts the value
0
0
(0) so that the spurious pole disappears. The 1/
√
s cut
does not introduce any additional singularity in the partial
wave since its square-root cut falls right on the left cut. At
the same time this term is suppressed in the physical region
and thus barely affects the fit and resonance pole positions.
As shown below, for this wave it will be enough to set N = 5
to obtain a good overall χ
2
/d.o. f. in the elastic region.
The conformal variable is defined as
ω(s) =
√
s − α
√
s
0
− s
√
s + α
√
s
0
− s
, (7)
where s
0
corresponds to the highest value of s where the
expansion is real and then α sets the center of the conformal
expansion. We have found in practice that the S0 wave is
more conveniently described if the conformal expansion, by
becoming imaginary, introduces some inelasticity above the
K
¯
K threshold [50]). Thus, we choose s
0
= 4m
2
K
with α =
1 for simplicity, so that the expansion center lies near 0.7
GeV. Hence, between K
¯
K threshold and 1.4 GeV, the
0
0
(s)
function will be complex, which effectively introduces an
inelasticity. We tried otherwise, with a higher s
0
, so that no
123