Skip to content

Calculate the periods and effective coupling

The curve has supplied two contour integrals. We now turn them into functions that can be evaluated, differentiated, and continued between vacua. Two methods will meet: direct integration and an ordinary differential equation. Agreement between them tests more than either method alone.

Read the curve and its normalized differential first. We retain pure four-dimensional N=2 SU(2)SU(2) Yang–Mills, with eight real supercharges preserved by its Coulomb vacua and one generic low-energy Abelian vector multiplet. The answer concerns the exact two-derivative Wilsonian data. It does not determine the entire particle spectrum or justify integrating out a particle at a point where it becomes massless.

Write s=Λ2>0s=\Lambda^2>0 and C=2/(2π)C=\sqrt2/(2\pi). Our curve and differential are

y2=(x2−s2)(x−u),λ=C(x−u)dxy,Π=(aD,a)T.(SW11.1)y^2=(x^2-s^2)(x-u),\qquad \lambda=C(x-u)\frac{dx}{y},\qquad \Pi=(a_D,a)^T. \tag{SW11.1}

Initially u>su>s is real. The electric cycle goes from −s-s to ss where yy is negative real and returns on the other sheet. The magnetic cycle goes from ss to uu where yy is positive imaginary and returns on the other sheet. These are the preceding lesson’s cycles, with geometric intersection γe⋅γD=+1\gamma_e\cdot\gamma_D=+1; the original paper uses the opposite sign convention for its intersection symbol. No contour is reversed here.

Both legs contribute, so

a=2π∫−ssu−xs2−x2 dx,aD=i2π∫suu−xx2−s2 dx.(SW11.2)a=\frac{\sqrt2}{\pi}\int_{-s}^{s} \frac{\sqrt{u-x}}{\sqrt{s^2-x^2}}\,dx,\qquad a_D=\frac{i\sqrt2}{\pi}\int_s^u \sqrt{\frac{u-x}{x^2-s^2}}\,dx . \tag{SW11.2}

In particular a>0a>0 and aD/i>0a_D/i>0. These are the scale-restored original period integrals, (6.22)–(6.23). The units are [u]=[s]=2[u]=[s]=2 and [a]=[aD]=1[a]=[a_D]=1.

Define two complete elliptic integrals, which here mean the following ordinary definite integrals:

K(m)=∫0π/2dθ1−msin⁡2θ,E(m)=∫0π/21−msin⁡2θ dθ,m=2su+s.(SW11.3)K(m)=\int_0^{\pi/2}\frac{d\theta}{\sqrt{1-m\sin^2\theta}}, \qquad E(m)=\int_0^{\pi/2}\sqrt{1-m\sin^2\theta}\,d\theta, \qquad m=\frac{2s}{u+s}. \tag{SW11.3}

We use the parameter mm, not a modulus whose square is the argument. For u>su>s, 0<m<10<m<1, and all square roots in (SW11.3) are positive.

For the electric integral put x=−s+2ssin⁡2θx=-s+2s\sin^2\theta. Then dx/s2−x2=2 dθdx/\sqrt{s^2-x^2}=2\,d\theta and u−x=(u+s)(1−msin⁡2θ)u-x=(u+s)(1-m\sin^2\theta). For the magnetic integral put x=s+(u−s)sin⁡2θx=s+(u-s)\sin^2\theta, and then replace θ\theta by π/2−θ\pi/2-\theta. Its integral without i2/πi\sqrt2/\pi becomes 2(u−s)∫0π/2sin⁡2θ dθ/(u+s)[1−(1−m)sin⁡2θ]2(u-s)\int_0^{\pi/2}\sin^2\theta\,d\theta/\sqrt{(u+s)[1-(1-m)\sin^2\theta]}. The identity K(q)−E(q)=q∫sin⁡2θ dθ/1−qsin⁡2θK(q)-E(q)=q\int\sin^2\theta\,d\theta/\sqrt{1-q\sin^2\theta} therefore gives

a=22(u+s)πE(m),aD=2i2(u+s)π[K(1−m)−E(1−m)].(SW11.4)\boxed{ a=\frac{2\sqrt{2(u+s)}}{\pi}E(m),\qquad a_D=\frac{2i\sqrt{2(u+s)}}{\pi} \big[K(1-m)-E(1-m)\big].} \tag{SW11.4}

This evaluates the contours in terms of explicitly defined real functions; it has not replaced the cycles by a new normalization.

A prime below means d/dud/du at fixed Λ\Lambda. Differentiating the defining integrals gives dE/dm=(E−K)/(2m)dE/dm=(E-K)/(2m) and dK/dm=E/[2m(1−m)]−K/(2m)dK/dm=E/[2m(1-m)]-K/(2m). The second relation follows by adding the integral of d[sin⁡θcos⁡θ/1−msin⁡2θ]/dθd[\sin\theta\cos\theta/\sqrt{1-m\sin^2\theta}]/d\theta, whose endpoint contribution vanishes. Using m′=−m/(u+s)m'=-m/(u+s) in (SW11.4) yields

a′=2πu+sK(m),aD′=i2πu+sK(1−m),τeff=aD′a′=iK(1−m)K(m).(SW11.5)a'=\frac{\sqrt2}{\pi\sqrt{u+s}}K(m),\qquad a_D'=\frac{i\sqrt2}{\pi\sqrt{u+s}}K(1-m),\qquad \boxed{\tau_{\rm eff}=\frac{a_D'}{a'} =i\frac{K(1-m)}{K(m)}.} \tag{SW11.5}

Both real integrals on the right are positive. Thus Im⁡τeff>0\operatorname{Im}\tau_{\rm eff}>0 on this real ray, as required for a positive gauge kinetic term. On a general smooth fiber, the preceding lesson’s bilinear identity gives the same positivity statement in a valid local frame.

The derivative ratio follows from aD=∂F/∂aa_D=\partial\mathcal F/\partial a and τeff=∂2F/∂a2\tau_{\rm eff}=\partial^2\mathcal F/\partial a^2. It is not aD/aa_D/a. For a worked value choose u=3su=3s, so m=1/2m=1/2:

QuantityValue in the stated branch
a/Λa/\Lambda2.432013446852.43201344685
aD/Λa_D/\Lambda0.906493919846 i0.906493919846\,i
Λa′\Lambda a'0.4173134208370.417313420837
τeff\tau_{\rm eff}ii exactly

Here aD/a≃0.372734 ia_D/a\simeq0.372734\,i would give the wrong coupling. The equality of the two KK integrals fixes τ=i\tau=i without numerical approximation.

Tachikawa’s geometric argument, (4.3.5)–(4.3.10), also computes a derivative ratio. Its curve, electric units and periods differ: its weak electric period is u\sqrt u, and its displayed large-loop shear is 44, rather than our 2u\sqrt{2u} and 22. We use that discussion for the geometric method, not to import its coefficients.

A Picard–Fuchs equation is a differential equation for integrals of a differential over continuously transported cycles. We can derive the needed one directly, without assuming a course on such equations.

At fixed xx, the chosen representative obeys ∂uλ=−C dx/(2y)\partial_u\lambda=-C\,dx/(2y) and ∂u2λ=−C dx/[4(x−u)y]\partial_u^2\lambda=-C\,dx/[4(x-u)y]. Now calculate the following combination:

(u2−s2)∂u2λ+14λ=C4x2−2ux+s2(x−u)y dx=C2 dx ⁣(yx−u).(SW11.6)\begin{aligned} (u^2-s^2)\partial_u^2\lambda+\frac14\lambda &=\frac C4\frac{x^2-2ux+s^2}{(x-u)y}\,dx\\ &=\frac C2\,d_x\!\left(\frac{y}{x-u}\right). \end{aligned} \tag{SW11.6}

For the second equality, differentiate y2=(x2−s2)(x−u)y^2=(x^2-s^2)(x-u) with respect to xx. The function y/(x−u)y/(x-u) is single-valued and meromorphic on the surface. Its exact differential integrates to zero around a closed cycle, even though a projection of that cycle uses two sheets. We transport and, where needed, deform the cycle away from poles before differentiating.

Consequently each period p=ap=a or aDa_D satisfies

(u2−Λ4)p′′+14p=0.(SW11.7)\boxed{(u^2-\Lambda^4)p''+\frac14p=0.} \tag{SW11.7}

This is the equation obtained by the different reconstruction in Bilal, (6.4), (6.11)–(6.13), with the scale restored. All terms have dimension one. The equation is singular precisely at the two finite degenerations; its large-uu behavior also matters.

At an ordinary point u0≠±su_0\ne\pm s, a second-order linear equation needs two initial numbers, p(u0)p(u_0) and p′(u0)p'(u_0). It determines a unique local solution. The electric and magnetic contours fix two different sets of these data. An equation alone, without those data, would not identify the physical periods.

For example, put p=∑n≥0cn(u−u0)np=\sum_{n\ge0}c_n(u-u_0)^n. Substitution gives the practical recurrence

(u02−s2)(n+2)(n+1)cn+2+2u0n(n+1)cn+1+[n(n−1)+14]cn=0.(SW11.8)(u_0^2-s^2)(n+2)(n+1)c_{n+2} +2u_0n(n+1)c_{n+1} +\left[n(n-1)+\frac14\right]c_n=0. \tag{SW11.8}

Starting with c0=p(u0)c_0=p(u_0) and c1=p′(u0)c_1=p'(u_0) determines every higher coefficient. Evaluate within the distance to the nearest singularity, then use the resulting value and derivative as data for another disk. This continues the chosen solution, rather than resetting a square root or logarithm on each step.

There is also a conserved diagnostic. Define the Wronskian W=aaD′−aDa′W=a a_D'-a_D a'. Equation (SW11.7) gives W′=aaD′′−aDa′′=0W'=a a_D''-a_D a''=0 because it has no first-derivative term. The normalization below fixes

W=2iπ.(SW11.9)W=\frac{2i}{\pi}. \tag{SW11.9}

A drifting numerical Wronskian signals an integration problem. A constant but incorrect value can instead signal incorrect initial normalization.

Recover weak coupling with its fixed scale

Section titled “Recover weak coupling with its fixed scale”

For real u/s≫1u/s\gg1, expand the electric integrand in (SW11.2). Its odd powers of xx integrate to zero. The normalized even moments are ⟨x2⟩=s2/2\langle x^2\rangle=s^2/2 and ⟨x4⟩=3s4/8\langle x^4\rangle=3s^4/8, giving

a=2u[1−s216u2−15s41024u4+O ⁣(s6u6)].(SW11.10)a=\sqrt{2u}\left[ 1-\frac{s^2}{16u^2} -\frac{15s^4}{1024u^4} +O\!\left(\frac{s^6}{u^6}\right)\right]. \tag{SW11.10}

The logarithm in the magnetic period comes from an endpoint. To see its constant as well as its coefficient, write K(1−ϵ)=∫0∞dt/(1+t2)(1+ϵt2)K(1-\epsilon)=\int_0^\infty dt/\sqrt{(1+t^2)(1+\epsilon t^2)}. Split the integral at 1≪T≪ϵ−1/21\ll T\ll\epsilon^{-1/2}. The small-tt part tends to log⁡(2T)\log(2T) and the large-tt part to log⁡[2/(ϵT)]\log[2/(\sqrt\epsilon T)]. The sum is log⁡(4/ϵ)\log(4/\sqrt\epsilon). Expanding the regular remainder gives K(1−ϵ)=log⁡(4/ϵ)+O(ϵ∣log⁡ϵ∣)K(1-\epsilon)=\log(4/\sqrt\epsilon)+O(\epsilon|\log\epsilon|) and E(1−ϵ)→1E(1-\epsilon)\to1.

Apply this with ϵ=m\epsilon=m in (SW11.4). The leading result is aD=(i/π)2u[log⁡(8u/s)−2+o(1)]a_D=(i/\pi)\sqrt{2u}[\log(8u/s)-2+o(1)]. To obtain the next correction without a delicate endpoint expansion, insert aD=(i/π)2u[L−2+(cL+d)s2/u2+⋯ ]a_D=(i/\pi)\sqrt{2u}[L-2+(cL+d)s^2/u^2+\cdots] into (SW11.7), where L=log⁡(8u/s)L=\log(8u/s). The coefficients of LL and of the constant give 4c+1/4=04c+1/4=0 and −4c+4d−1/2=0-4c+4d-1/2=0. Hence

aD=i2uπ[L−2−s216u2(L−1)+O ⁣(s4u4L)],τeff=iπ[L−5s216u2+O ⁣(s4u4)].(SW11.11)\begin{aligned} a_D&=\frac{i\sqrt{2u}}{\pi} \left[L-2-\frac{s^2}{16u^2}(L-1) +O\!\left(\frac{s^4}{u^4}L\right)\right],\\ \tau_{\rm eff}&=\frac{i}{\pi} \left[L-\frac{5s^2}{16u^2} +O\!\left(\frac{s^4}{u^4}\right)\right]. \end{aligned} \tag{SW11.11}

Differentiating both period expansions verifies the second line; the logarithmic correction cancels in the ratio. Their leading terms give W=2i/πW=2i/\pi. The constants 88 and −2-2 belong to our already fixed curve scale and period frame. An earlier running-coupling formula that leaves a finite scale matching unspecified cannot fix these constants by itself.

The expansion parameter is Λ4/u2\Lambda^4/u^2, not merely the inverse of a dimensionful uu. For complex large uu, these expansions use the square root and logarithm continued from the positive ray in a fixed sector avoiding their cuts. The displayed error estimates are asymptotic there; they are not estimates near u=±su=\pm s.

Approach the monopole point in the same branch

Section titled “Approach the monopole point in the same branch”

Let δ=u−s>0\delta=u-s>0 with δ/s≪1\delta/s\ll1. In the magnetic integral use x=s+δtx=s+\delta t:

aD=i2 δπ∫01(1−t)/t2s+δt dt=iδ2Λ[1−δ16s+3δ2256s2+O ⁣(δ3s3)].(SW11.12)a_D=\frac{i\sqrt2\,\delta}{\pi} \int_0^1\frac{\sqrt{(1-t)/t}}{\sqrt{2s+\delta t}}\,dt =\frac{i\delta}{2\Lambda} \left[1-\frac{\delta}{16s} +\frac{3\delta^2}{256s^2} +O\!\left(\frac{\delta^3}{s^3}\right)\right]. \tag{SW11.12}

For example ∫01(1−t)/t dt=π/2\int_0^1\sqrt{(1-t)/t}\,dt=\pi/2, and its weighted averages of t,t2t,t^2 are 1/4,1/81/4,1/8. These give every displayed coefficient. In particular aD′(s)=i/(2Λ)≠0a_D'(s)=i/(2\Lambda)\ne0: the vanishing period is a good local coordinate.

The electric period stays finite. At δ=0\delta=0 its integral gives a(s)=4Λ/πa(s)=4\Lambda/\pi. In (SW11.5), 1−m=δ/(2s+δ)1-m=\delta/(2s+\delta), so the same endpoint logarithm now gives a′=[log⁡(32s/δ)+O((δ/s)∣log⁡(δ/s)∣)]/(2πΛ)a'=[\log(32s/\delta)+O((\delta/s)|\log(\delta/s)|)]/(2\pi\Lambda). Integrating from δ=0\delta=0 fixes the constant term:

a=4Λπ+δ2πΛ[log⁡32sδ+1]+O ⁣(δ2Λ3∣log⁡(δ/s)∣),τeff∼iπlog⁡(32s/δ),τD=−1τeff∼iπlog⁡32sδ.(SW11.13)\begin{aligned} a&=\frac{4\Lambda}{\pi} +\frac{\delta}{2\pi\Lambda} \left[\log\frac{32s}{\delta}+1\right] +O\!\left(\frac{\delta^2}{\Lambda^3}|\log(\delta/s)|\right),\\ \tau_{\rm eff}&\sim\frac{i\pi}{\log(32s/\delta)}, \qquad \tau_D=-\frac1{\tau_{\rm eff}} \sim\frac{i}{\pi}\log\frac{32s}{\delta}. \end{aligned} \tag{SW11.13}

The limits agree with the original local calculation, (6.27)–(6.30). The dual electric frame becomes weakly coupled. At the point itself the extra light monopole hypermultiplet must be retained, as explained in the strong-coupling lesson. A small electric Im⁡τ\operatorname{Im}\tau is not a negative kinetic term and is not permission to keep a photon-only action through its singularity.

Numerical evaluation should compare independently formulated quantities. Removing square-root endpoints makes the contour quadrature well behaved: a=(2/π)∫0πu−scos⁡θ dθa=(\sqrt2/\pi)\int_0^\pi\sqrt{u-s\cos\theta}\,d\theta and aD=(2i2/π)(u−s)∫0π/2cos⁡2θ dθ/2s+(u−s)sin⁡2θa_D=(2i\sqrt2/\pi)(u-s)\int_0^{\pi/2}\cos^2\theta\,d\theta/\sqrt{2s+(u-s)\sin^2\theta}. These formulas and their differentiated integrals provide initial data at, for example, u=3su=3s.

Then continue those data with (SW11.8), independently of elliptic-integral evaluation. Compare at real points and at complex points reached within Re⁡u>s\operatorname{Re}u>s, where these interval integrals continue from their positive-real roots without ambiguity. Also compare WW, Im⁡τ\operatorname{Im}\tau, and the errors of the two asymptotic expansions as their small parameters decrease. Increasing arithmetic precision does not by itself control a truncated series: decrease the disk step or increase the retained order too.

Neither a finite collection of numerical points nor the ODE’s existence proves a stable-particle spectrum. The checks establish agreement of the normalized period calculations in their stated domain. Closed paths around punctures require the explicit continuation and loop conventions of the monodromy lesson.

Derive the s4/u4s^4/u^4 coefficient in (SW11.10) directly from the contour integral. Then invert the result through order s4/a6s^4/a^6 to find u(a)u(a) on the positive weak branch. Why does u=a2/2u=a^2/2 fail beyond the classical approximation?

Hint 1

The fourth term of the binomial expansion of 1−z\sqrt{1-z} is −5z4/128-5z^4/128. Average powers of x=scos⁡θx=s\cos\theta.

Hint 2

Square the period first and use an ansatz u=a2/2+As2/a2+Bs4/a6u=a^2/2+A s^2/a^2+B s^4/a^6. Every term must have dimension two.

Solution

The fourth moment is 3s4/83s^4/8, giving (−5/128)(3/8)=−15/1024(-5/128)(3/8)=-15/1024. Squaring the period gives a2=2u−s2/(4u)−13s4/(256u3)+O(s6/u5)a^2=2u-s^2/(4u)-13s^4/(256u^3)+O(s^6/u^5). Substituting the proposed inverse, the s2/a2s^2/a^2 coefficient requires 2A−1/2=02A-1/2=0, and the s4/a6s^4/a^6 coefficient requires 2B+A−13/32=02B+A-13/32=0. Thus

u(a)=a22+Λ44a2+5Λ864a6+O ⁣(Λ12a10).u(a)=\frac{a^2}{2}+\frac{\Lambda^4}{4a^2} +\frac{5\Lambda^8}{64a^6} +O\!\left(\frac{\Lambda^{12}}{a^{10}}\right).

All terms have dimension two. The coordinate uu is a gauge-invariant vacuum expectation value, whereas aa is a quantum period. Their classical relation sets a normalization at infinity; it does not remove quantum corrections. This inversion is valid for ∣a∣≫Λ|a|\gg\Lambda on the declared branch.

2. Test the local coordinate and the scalar metric

Section titled “2. Test the local coordinate and the scalar metric”

Using (SW11.12), find uu as a series in aDa_D through quadratic order. Then estimate G(u)=Im⁡τeff∣a′(u)∣2G(u)=\operatorname{Im}\tau_{\rm eff}|a'(u)|^2 on the ray δ>0\delta>0 near the monopole point. Here GG omits any overall action normalization. Does its divergence place the singular point at infinite distance in this two-derivative metric?

Hint 1

Set w=−2iΛaDw=-2i\Lambda a_D. Then w=δ−δ2/(16s)+⋯w=\delta-\delta^2/(16s)+\cdots.

Hint 2

Use the logarithmic formulas for both τ\tau and a′a'. A radial distance involves G dδ\sqrt{G}\,d\delta, not G dδG\,d\delta.

Solution

Inverting gives δ=w+w2/(16s)+O(w3/s2)\delta=w+w^2/(16s)+O(w^3/s^2), hence u=s−2iΛaD−aD2/4+O(aD3/Λ)u=s-2i\Lambda a_D-a_D^2/4+O(a_D^3/\Lambda). In particular du/daD∣0=−2iΛdu/da_D|_0=-2i\Lambda has the required dimension one and is nonzero.

With Lδ=log⁡(32s/δ)L_\delta=\log(32s/\delta), G∼Lδ/(4πΛ2)>0G\sim L_\delta/(4\pi\Lambda^2)>0. The distance integral is proportional to ∫0ϵlog⁡(32s/δ) dδ\int_0^\epsilon\sqrt{\log(32s/\delta)}\,d\delta. Set δ=32se−t\delta=32s e^{-t}: its tail is an exponentially convergent integral ∫∞e−tt dt\int^\infty e^{-t}\sqrt t\,dt. The point is at finite distance although the photon-only metric coefficient diverges. This is a statement about the limiting effective metric, not a reason to omit the light hypermultiplet.

3. Diagnose a reassuring but incomplete numerical test

Section titled “3. Diagnose a reassuring but incomplete numerical test”

A program returns a^=a\widehat a=a and a^D=aD+κa\widehat a_D=a_D+\kappa a, where κ\kappa is a real constant. It reports that both functions solve (SW11.7), their Wronskian is 2i/π2i/\pi, and their derivative ratio has positive imaginary part. Are these tests enough to establish the specific contour answer? Find a discriminating boundary test. Contrast this error with a^D=−aD\widehat a_D=-a_D.

Hint 1

The equation is linear. Compute the Wronskian and derivative ratio after adding a multiple of the other solution.

Hint 2

The assigned magnetic cycle shrinks at u=su=s, while the electric period approaches 4Λ/π4\Lambda/\pi.

Solution

Linearity preserves the equation. The added terms cancel in aa^D′−a^Da′a\widehat a_D'-\widehat a_Da', and τ^=τ+κ\widehat\tau=\tau+\kappa has the same imaginary part. But a^D(s)=4κΛ/π\widehat a_D(s)=4\kappa\Lambda/\pi, so the assigned vanishing-period condition forces κ=0\kappa=0. An integral κ\kappa could instead describe a different integral cycle basis, accompanied by the corresponding charge transformation. It is still not the contour specification of this problem. Arbitrary real κ\kappa need not even preserve the integral cycle lattice.

Changing only the magnetic sign gives W^=−2i/π\widehat W=-2i/\pi and Im⁡τ^<0\operatorname{Im}\widehat\tau<0 on the real ray. These diagnostics detect that orientation error immediately. Thus the ODE, positivity, Wronskian, and boundary normalization test different parts of the calculation; no one test replaces the rest.

We can calculate aDa_D, aa and the effective coupling at a smooth vacuum, and approach a degeneration in a controlled local coordinate. BPS masses and the stability of particles asks how Zγ=nmaD+neaZ_\gamma=n_ma_D+n_ea constrains particle masses and decays. A computable central charge for an integer row is necessary information; deciding whether that row labels an existing stable BPS particle requires more.

For a different test of the same low-energy data, recover a prepotential coefficient by instanton counting. Evaluate the first equivariant contribution, translate its electric and magnetic normalization, and compare it with the weak-coupling periods derived here.