16.3 Luescher analysis

The phase shift can be computed from the formula

cotδ=Z00(1,q2)π3/2γq where γ=E/ECM, q=kL/2π with k the scattering momentum k2=E2CM4m2=E2P24m2. The Energy in the center of mass is related to the one in a generic frame with total momentum P via E2CM=E2impP2. Eimp is the energy measured in the lattice E=EmeasuredEfreelatt+Efreecont Efreelatt=cosh1(cosh(m)+12(3i=14sin(p1i2)2))+cosh1(cosh(m)+12(3i=14sin(p2i2)2)).

For the Z function we use the rzeta package. The fit function for the phase shift is

kmcotδ=1a0m+r0m2k2m2

P[0]=am , P[1]=r0m

χ2/d.o.f.=0.976672P[0]=0.15212±(0.0014)P[1]=2.83566±(0.19) {C=(10.0009410.0009411)}

00.20.40.60.811.2−9−8.5−8−7.5−7−6.5−6−5.5
$k/m$$\frac{k}{m} \cot \delta $

16.3.1 ML in zeta func and M as normalization

χ2/d.o.f.=0.975176P[0]=0.152073±(0.0014)P[1]=2.83815±(0.19) {C=(10.0009380.0009381)}

00.20.40.60.811.2−9−8.5−8−7.5−7−6.5−6−5.5
$k/m$$\frac{k}{m} \cot \delta $

16.3.2 Using M

χ2/d.o.f.=2.31189P[0]=0.15966±(0.0036)P[1]=3.15533±(0.24) {C=(10.002830.002831)}

00.20.40.60.811.2−9−8.5−8−7.5−7−6.5−6−5.5−5
$k/m$$\frac{k}{m} \cot \delta $