derivation of integral representations of Jacobi ϑ functions


By rearranging the Fourier series of cos(ux), one obtains the series

πcos(ux)2usin(πu)=12u2+n=1(-1)ncos(nx)u2-n2

This equation which is valid for all real values of x such that -πxπ and all non-integral complex values of u. By comparison with the convergent series n=01/n2, it follows that this series is absolutely convergent. Note that this series may be viewed as a Mittag-Leffler partial fraction expansion.

Let y be a positive real number. Multiply both by 2ue-yu2 and integrate.

i-i+πcos(ux)e-yu2sin(πu)𝑑v=2i-i+e-yu2[12u2+n=0(-1)ncos(nx)u2-n2]u𝑑u

Because of the exponential, the integrand decays rapidly as ui± provided that u>0, and hence the integral converges absolutely. Make a change of variables v=u2

=Pe-yv[12v+n=1(-1)ncos(nx)v-n2]𝑑v

The contour of integration P is a parabolaPlanetmathPlanetmath in the complex v-plane, symmetric about the real axis with vertex at v=-1, which encloses the real axis. Its equation is v+1=2(v)2

Let Sm (m is an integer) be the straight line segment joining the points v=(i+m+1/2)2 and v=(i-m-1/2)2. Along this line segmentMathworldPlanetmath, we may bound the integrand in absolute valueMathworldPlanetmathPlanetmathPlanetmath as follows:

|n=1(-1)ncos(nx)v-n2|n=1(-1)n|v-n2|n=1(-1)n|vm-n2|

where vm=m2+m-3/4 is the point of intersectionMathworldPlanetmath of Sm with the real axis. To proceed further, we break up the last summation into two parts.

Since the squares closest in absolute value to vm are m2 and (m+1)2=m2+2m+1, it follows that |vm-n2||m-3/4| for all m,n. Hence, we have

i=12m1|vm-n2|2mm-3/48

When n>2m, we have n2(2m+1)2=4m2+4m+1>4m2+4m-3=4vm. Hence, |n2-vm|>3n2/4 and

n=2m+11|vm-n2|<43n=2m+11n2<43n=11n2=2π9

Finally 1/(2vm)<1/2 since vm>1 when m1. Also, |e-yv|=e-yv-e-yvm<e-ym2. From these observations, we conclude that

|Sme-yv[12v+n=1(-1)ncos(nx)v-n2]𝑑v|<e-ym2(1+8+2π9)Sm𝑑v=(4m+2)(9+2π9)e-ym2

Note that this quantity approaches 0 in the limit m.

Let Pm be the arc of the parabola P boundedPlanetmathPlanetmathPlanetmathPlanetmath by the endpointsMathworldPlanetmath of Sm. Together, Sm and Pm form a closed contour which encloses poles of the integrand. Hence, by the residue theoremMathworldPlanetmath , we have

Pme-yv[12v+n=1(-1)ncos(nx)v-n2]𝑑v+Sme-yv[12v+n=1(-1)ncos(nx)v-n2]𝑑v=
2πin=1m(-1)ncos(nx)e-n2y

Taking the limit m we obtain

Pe-yv[12v+n=1(-1)ncos(nx)v-n2]𝑑v=2πi(12+n=1(-1)ncos(nx)e-n2y)

Going back to the beginning of the proof, where the integral on the left hand side was expressed as an integral with respect to u, we obtain

i-i+πcos(ux)e-yu2sin(πu)𝑑v=2πi(12+n=1(-1)ncos(nx)e-n2y)

Making a change of variables x=2z,y=-iπτ and tidying up some, we obtain

i-i+cos(2uz)eiπτu2sin(πu)dv=i(1+2n=1(-1)neiπn2τcos(2nz))=iϑ4(z|τ)

Because of the initial assumptionPlanetmathPlanetmath about the Fourier series, we only know that this formulaMathworldPlanetmathPlanetmath is valid when τ is purely imaginary with strictly positive imaginary part and z is real and π/2<z<π/2. However, we can use analytic continuation to extend the domain of its validity. On the one hand, the theta function on the right-hand side is analyticPlanetmathPlanetmath for all z and all τ such that τ>0.

On the other hand, I claim that the integral on the left hand side is also an analytic function of z and τ whenever τ>0. To validate this claim, we need to examine the behaviour of the integrand as ui±. The contribution of the denominator is bounded;

|1sinπu|<c

for some constant c whenever u=1. The absolute value of the cosine in the numerator is easy to bound:

|cos(2uz)|e2|u||z|

To bound the remaining term, let us examine the argument of the exponential carefully:

(τu2)=2τu+τ(u)2-τ=τ((u+ττ)2-1-(ττ)2)

Therefore, if |u|>1+3|τ|/(τ), it will be the case that (τu2)τ(u)2/9, and so

|eiπτu2|=e-π(τu2)e-πτ(u)2/9

Taken together, the estimates of the last paragraph imply that

|i+Ri+cos(2uz)eiπτu2sin(πu)|<ci+Ri+e2|u||z|-πτ(u)2/9

when R>1+3|τ|/(τ). If we impose the further conditions

R>180|z|πτ  R2>180|z|πτ  ,

it will be the case that

2|u||z|-πτ(u)2/9<2u|z|+2|z|-πτ(u)2/9<
(2u|z|-πτ(u)2/180)+(2|z|-πτ(u)2/180)-πτ(u)2/10<
-πτ(u)2/10  ,

and hence

|i+Ri+cos(2uz)eiπτu2sin(πu)𝑑u|<ci+Ri+e-πτ(u)2/10𝑑u<5cπτRe-πτR2/10  .

Likewise, under the same restrictionPlanetmathPlanetmathPlanetmath on R,

|i-i-Rcos(2uz)eiπτu2sin(πu)𝑑u|<ci+Ri+e-πτ(u)2/10𝑑u<5cπτRe-πτR2/10  .

Since the contour of integration is compactPlanetmathPlanetmath and the integrand is analytic in a neighborhoodMathworldPlanetmathPlanetmath of the contour,

i-Ri+Rcos(2uz)eiπτu2sin(πu)𝑑u

will be an analytic function of z and τ. Suppose that z and τ are restricted to bounded regions of the complex plane and that, furthermore, Imτ is positive and bounded away from zero. Then the inequalities of the last paragraph imply that the integral converges uniformly as R, and hence

i-i+cos(2uz)eiπτu2sin(πu)𝑑u

is an analytic function of u and z in the domain τ>0.

Thus, by the fundamental theorem of analytic continuation, we may conclude that

i-i+cos(2uz)eiπτu2sin(πu)dv=i(1+2n=1(-1)neiπn2τcos(2nz))=iϑ4(z|τ)

throughout this domain.

Finally, integral representations of the remaining three theta functions may be easily obtained from this one by adding the appropriate half-quasiperiods to z.

Title derivation of integral representations of Jacobi ϑ functions
Canonical name DerivationOfIntegralRepresentationsOfJacobivarthetaFunctions
Date of creation 2013-03-22 14:39:54
Last modified on 2013-03-22 14:39:54
Owner rspuzio (6075)
Last modified by rspuzio (6075)
Numerical id 26
Author rspuzio (6075)
Entry type Derivation
Classification msc 33E05