Theorem (Periods integrals and type II solutions)
Consider the mean field equation4u+eu=ρδ0on E=C/Λ.
I If solutions exist forρ=8nπ, then there is a unique even solution within each type II scaling family. (`=2n, an+i=−ai.)
I The n-th equation is given byR
Lig∈√
−1R, which is equivalent to
∑
n i=1∇G(ai) =0.
I The n−1 algebraic equations:
Theorem (Green/polynomial system)
Forρ=8nπ, n∈N, the n equations for a={a1, . . . , an}are precisely
℘0(a1)℘r(a1) +· · · + ℘0(an)℘r(an) =0, where r=0, . . . , n−2, and∇G(a1) +· · · + ∇G(an) =0.
Theorem (Hyperelliptic geometry/Lam´e curve)
For xi := ℘(ai), yi:= ℘0(ai), the first n−1 algebraic equations
∑
yixri =0, r=0, . . . , n−2,defines an affine hyperelliptic curve under the 2 to 1 map a7→∑℘(ai): Xn:={(xi, yi)} ⊂SymnE−→ (x1+· · · +xn)∈P1.
I The proof relies on its relation to Lam´e equations:
I This is a long calculation via the polynomial system (omitted).
I Idea of proof on the hyperelliptic structure on Xn. I Consider y2=p(x) =4x3−g2x−g3, where
(x, y) = (℘(z),℘0(z)), and we set(xi, yi) = (℘(ai),℘0(ai)). Consider a basis of solutions to the Lam´e equation
w00= (n(n+1)℘(z) +B)w (for some B) given by wa(z)and w−a(z).
I Let X(z) =wa(z)w−a(z). By the addition theorem, X(z) = (−1)n
∏
n i=1σ(z+ai)σ(z−ai)
σ(z)2σ(ai)2 = (−1)n
∏
n i=1(℘(z)− ℘(ai)).
That is, X(x) = (−1)n∏ni=1(x−xi)is a polynomial in x.
I Key: X(z)satisfies the second symmetric power of the Lam´e equation:
d3X
dz3 −4(n(n+1)℘ +B)dX
dz −2n(n+1)℘0X=0.
I Hence X(x)is a polynomial solution, in variable x, to
p(x)X000+32p0(x)X00−4((n2+n−3)x+B)X0−2n(n+1)X=0.
I X is determined by B and certain initial conditions.
I Write X(x) = (−1)n(xn−s1xn−1+· · · + (−1)nsn), this translates to a linear recursive relation for µ=0,· · ·, n−1:
0=2(n−µ)(2µ+1)(n+µ+1)sn−µ
−4(µ+1)Bsn−µ−1
+12g2(µ+1)(µ+2)(2µ+3)sn−µ−2
−g3(µ+1)(µ+2)(µ+3)sn−µ−3. I We set s0=1.
I For µ=n−1 we get B= (2n−1)s1as expected.
I Thus all s2,· · ·, sn, X(z), are determined by s1, i.e. by B, alone.
I In fact, a slightly more work shows that the set a={ai}is also determined by B up to sign. Hence a7→Bais 2 to 1. QED
Theorem (Chai-Lin-W 2012)
I There is a natural projective compactification ¯Xn ⊂SymnE as a, possibly singular, hyperelliptic curve defined by
C2= `n(B, g2, g3) =4Bs2n+4g3sn−2sn−g2sn−1sn−g3s2n−1, in affine coordinates(B, C), where
sk=sk(B, g2, g3) =rkBk+· · · ∈Q[B, g2, g3]
is an universal polynomial of homogeneous degree k with deg g2=2, deg g3=3, and B= (2n−1)s1.
I Thus deg`n=2n+1 and ¯Xnhas arithmetic genus g=n.
I The curve ¯Xnis smooth except for a finite number ofτ, namely the discriminant loci of`n(B, g2, g3), so that`n(B)has multiple roots. In particular ¯Xnis smooth for rectangular tori.
(Continued.)
I The 2n+2 branch points a∈X¯n\Xnare characterized by−a=a.
I {−ai} ∩ {ai} 6=∅⇒ −a=a.
I Also 0∈ {ai} ⇒a= (0, 0,· · ·, 0).
I By setting(xi, yi) = (t3i, 2t2i), the limiting system at a=0n:
∑
ni=1t2r+1i =0, r=1, . . . , n−1,has a unique solution with ti6=0 and ti6= −tjinPn−1up to permutations.
I Meaning of C (I):Applying Cramer’s rule to the n−1 linear equations∑ni=1xkiyi=0 in yi’s, there is a constant C∈C×such that
yi= C
∏j6=i(xi−xj), i=1, . . . , n.
I Meaning of C (II):Let w1, w2be two ind. solutions of w00=Iw.
I From
I The constant terms give the hyperelliptic equation in(B, C). I In particular, C=0 if and only if wa=w−a, i.e. a=−a. These
are the branch points of ¯Xn.
I Definition:Denote by Yn=X¯n\{0n}the affine hyperelliptic curve defined by
C2= `n(B, g2, g3).
I Now we study the last equation on ¯Xn:
I Consider the rational function on En:
zn(a1, . . . , an):=ζ(a1+· · · +an)−
∑
n i=1ζ(ai). (It is periodic in each variable.)
I Let ai=tiω1+siω2, then
I It is thus crucial to study the branched covering map
σ : ¯Xn →E, a7→σ(a):=
∑
n i=1ai.
Theorem (Lin–W 2013, new pre-modular functions)
(1) The mapσ has degree equals 12n(n+1).
(2) There is a universal (weighted homogeneous) polynomial Wn(x)∈C[g2, g3,℘(σ),℘0(σ)][x]of degree 12n(n+1)such that
Wn(zn) =0.
In fact, zn∈K(X¯n)is a primitive generator for the field extension K(X¯n)over K(E).
(3) The function Zn(σ; τ):=Wn(Z)is pre-modular of weight
1
2n(n+1). That is, it is modular wrt.Γ(N)ifσ∈Eτ[N].
I Idea of proof for (1):Apply Theorem of the Cube: For any three morphisms f , g, h : Vn−→E and L∈Pic E,
(f+g+h)∗L∼= (f+g)∗L⊗ (g+h)∗L⊗ (h+f)∗L
⊗f∗L−1⊗g∗L−1⊗h∗L−1.
I Apply to the case Vn⊂Enwhich is the ordered n-tuples so that Vn/Sn=X¯n, and deg L=1. We prove inductively that the map
fk(a):=a1+· · · +ak
has degree 12k(k+1)n!. It is not hard to check for k=1, 2.
I From k to k+1, we let f =fk−1, g(a) =ak, and h(a) =ak+1. I Then fk+1has degree n! times
1
2k(k+1) +3+12k(k+1)−12(k−1)k−1−1= 12(k+1)(k+2).
I Idea of proof of (2):Major tool: tensor product of two Lam´e equations w00=I1w and w0=I2w, where I=n(n+1)℘(z), I1=I+Baand I2=I+Bb.
I For ¯Xn(τ)smooth, and a general point σ0∈E, we need to show that the 12n(n+1)points on the fiber of ¯Xn→E above σ0has distinct znvalues. It is enough to show that for σ(a) =σ(b) =σ0, the condition∑ ζ(ai) =∑ ζ(bi)implies Ba=Bb(and then a=b).
I If w001 =I1w1and w002 =I2w2, then the product q=w1w2satisfies q0000−2(I1+I2)q00−6I0q0+ ((Ba−Bb)2−2I00)q=0.
I If a6=b, by addition law we find that Q=waw−b+w−awbis an even elliptic function solution, namely a polynomial in x= ℘(z). This leads to strong constraints on the corresponding 4-th order ODE in variable x, and eventually leads to a contradiction for generic choices of σ0.
Indeed,
As an even elliptic function, Q takes the form
Q(x) =C The xn+2terms agree automatically. The xn+1degree gives
∑
℘(ci) =s1= 12Ba+Bb
2n−1 = 12(
∑
℘(ai) +∑
℘(bi)).I Inductively the xn+2−icoefficient in (4) gives recursive relations to solve siinterns of Ba+Bb,(Ba−Bb)2and g2, g3for i=1, . . . , n.
I Indeed
si=si(Ba+Bb,(Ba−Bb)2, g2, g3) =Ci(Ba+Bb)i+· · · is homogeneous of degree i if we assign deg Ba=deg Bb=1 and deg g2=2, deg g3=3.
I There are two remaining consistency equations F1=0, F0=0 coming from the x1and x0coefficients in (4).
I In fact(Ba−Bb)2is a factor of both equations and we may write F1(Ba, Bb) = (Ba−Bb)2d1G1(Ba, Bb)and
F0(Ba, Bb) = (Ba−Bb)2d0G0(Ba, Bb). I If Ba6=Bb(i.e∑℘(ai)6=∑℘(bi)), then
G1(Ba, Bb) =0, G0(Ba, Bb) =0,
which has only a finite number of solutions(Ba, Bb)’s, i.e. Eτ’s.
Example (of compatibility equations for n = 2)
For n=2 we have s1= 16(Ba+Bb)and
s2= 361(Ba+Bb)2+721(Ba−Bb)2−14g2.
The first compatibility equation from x1is (after substituting s1)
1
6(Ba−Bb)2(Ba+Bb) =0.
The second compatibility equation from x0is
(Ba−Bb)2(361(Ba+Bb)2+721(Ba−Bb)2−16g2) =0.
If Ba 6=Bbthen Bb=−Baand then we can solve Ba, Bb: B2a=3g2=⇒ ℘(a1) + ℘(a2) =±pg2/3.
Such a∈X¯2indeed lies at the branch loci of the Lam´e curve.
Example (of new pre-modular forms for n = 2)
For z2(a1, a2) =ζ(a1+a2)−ζ(a1)−ζ(a2), on X2: z32(a)−3℘(a1+a2)z2(a)− ℘0(a1+a2) =0.
On E2it has one more term−12(℘0(a1) + ℘0(a2)). Thus, Z2(σ; τ) =W2(Z) =Z3−3℘(σ)Z− ℘0(σ).
Example (n = 3)
For z=z3(a) =ζ(a1+a2+a3)−ζ(a1)−ζ(a2)−ζ(a3), on X3: z6−15℘z4−20℘0z3+ (274g2−45℘2)z2−12℘0℘z−54℘02 =0.
Thus, Z3(σ; τ) =W3(Z).
I Key point: Z1≡Z=−4π∇G is the Hecke modular function.
The critical point equation (⇐⇒type II solutions of MFE) is transformed into zero of pre-modular forms.
I For general n≥1, we have the equivalences:
• Solution u to MFE for ρ=8πn.
• Periods integral Z
Ljg∈√
−1R (=ωjcoordinates of∑ ai.)
• Green equation
∑
n i=1∇G(ai) =0 on Xn.
• zn(a) =Z(σ(a)).
• Non-trivial zero of Zn(σ; τ):=Wn(Z).
I Need to prove the last one. Notice that the branch point a∈Yn\Xn(a6= −a) satisfies the Green equation trivially.
I The second technique used in ρ=8π is to use the method of continuity to connect to the known case ρ=4π by establishing the non-degeneracy of linearized equations.
I For general ρ, such a non-degeneracy statement is out of reach.
However, since solutions uηalways exist for ρ=8πη, η 6∈N, it is natural to study the limiting behavior of uηas η→n. If the limit does not blow up, it converges to a solution u for ρ=8πn.
I For the blow-up case, we have the connection between the blow-up set and the hyperelliptic geometry of Yn→P1: I
Theorem
Suppose that S={p1,· · ·, pn}is the blow-up set of a sequence of solutions ukto withρk→8πn as k→∞, then S∈Yn. Moreover,
(1) Ifρk6=8πn then S is a branch point (a=−a) of Yn. (2) Ifρk=8πn for all k, then S is not a branch point of Yn.
I To go deeper, need to know the converse statement: for which p∈Yn\Xncan we construct a blow-up sequence with blow-up set p? The Morse type of p is fundamental.
Theorem
Suppose that the pair of non half-period critical points{±p}of G exists, the
±p are the minimal points of G.
I In fact our proof shows that any solution for ρ=8π must be a minimizer of the non-linear functional
J8π(u) =1 2 Z
E|∇u|2−8π log Z
Ee−8πG+u on u∈H1(E)∩ {u|R
Eu=0}.
Corollary
Forτ∈Ω5, all the three half periods are (non-degenerate) saddle points.
I If ukis a blow-up sequence with ρ=ρk→8π (as k→∞), ρk6=8π for large k, then the blow-up point q is a half period.
I Asymptotically
ρk−8π= (D(q) +o(1))e−λk (5) where λk=maxEτukand
D(q):= Z
Eτ
h(z)e8π( ˜G(z,q)−φ(q))−h(q)
|z−q|4 − Z
Ecτ
h(q)
|z−q|4. Here h(z) =e−8πG(z), ˜G(z, q)is the regular part of the Green function, and φ(q) =G˜(q, q).
I The sign of D(q)determines the direction where the bubbling may take place, namely ρk <8π or ρk>8π.
Theorem (Lin–W)
For any half period q∈Eτ,τ=a+bi, we have
D(q) =−4π2be−8πG(q)det D2G(q). (6)
I Hence D(q) >0 if q is a saddle point. In particular if τ∈Ω5
then D(q) >0 for all half-periods since they are all saddle.
I Since the extra critical point p (reps.−p) is a discrete minimal point, the index of∇G at p (reps.−p) is 1. By the Hopf–Poincar´e index theorem,
−1=χ(Eτ\{0}) =2+
∑
3 i=1ind1 2ωi∇G.
Since12ωiis non-degenerate,∇G has index±1 at it. Hence the index must be−1 for all i=1, 2, 3. This implies that12ωiis a saddle point for all i.
I Combining with a recent technique in analyzing uniqueness of blow-up solutions by Lin–Yan, we may classify all solutions to the mean field equation for ρ∈ (0, 8π+e0)for some e0>0:
Theorem (Lin–W)
(i) Ifτ∈Ω3then the MFE has only one solution forρ<8π, no solution forρ=8π, and two solutions for 8π<ρ<8π+e0. (ii) Ifτ∈Ω5then the MFE has only one solution forρ<8π,
infinitely many solutions forρ=8π, and four solutions for 8π<ρ<8π+e0.
I MFE with ρ=12π has exactly two solutions on Eτfor τ6=eπi/3. Hence when τ∈Ω5the bifurcation diagram is complicate for ρ∈ (8π, 12π). It is a natural question whether MFE has exactly two solutions for ρ∈ (8π, 16π)when τ∈Ω3.
I The Theorem also reflects the difficulty in the study the corresponding Lam´e equation for the case η6∈ 12N.
I The hyperelliptic curve Ynis parametrized by(B, C)with
I As in the case n=1, these two invariants are related to a geometric quantity D(p)derived from the blow-up analysis of solutions ukwith ρk→8πn. Let p= (p1,· · ·, pn)with
{p1,· · ·, pn}being the blow-up set of uk. Then ρk−8πn= (D(p) +o(1))e−λk, λk:=max
E uk. I The analytic expression of D(p)is rather complicate. However,
its geometric meaning is reflected in the following
Theorem
Moreover, cp>0 except for a finite set of tori.
(This has been verified for n=1, 2.)