Skip to content

What does a fourth derivative of a sphere partition function measure? It contains connected four-point information, but it also remembers how the action depends nonlinearly on the source. Ignoring that dependence can produce a finite, reproducible wrong answer. We will expose the discrepancy in the smallest ABJM theory, then change the sphere shape while keeping the mass source fixed.

Prepare with the reproduction method, ABJM masses and squashing, and current normalizations. This is a reproduction and controlled training extension. It does not determine an unknown interacting four-point function or establish a new duality.

Use compact U(1)1×U(1)−1U(1)_1\times U(1)_{-1} ABJM: N=1N=1, levels exactly +1,−1+1,-1, both gauge centers retained, and no extra quotient. Its local undeformed conformal theory is the free maximally supersymmetric multiplet, containing four complex scalars and their fermions. There are sixteen real Poincaré supercharges, with sixteen conformal partners at the conformal point. The sphere calculation uses the manifest three-dimensional N=2\mathcal N=2 construction; the generic toric background supplies one rigid label direction of each R charge. This is a statement about the constructed subalgebra, not a count of every symmetry at a special undeformed point.

Keep the prerequisite’s toric metric, with ℓ1=Lb\ell_1=Lb, ℓ2=L/b\ell_2=L/b and f2=ℓ12sin⁡2χ+ℓ22cos⁡2χf^2=\ell_1^2\sin^2\chi+\ell_2^2\cos^2\chi. Its complete rigid sources are

H=if,V=C=0,A(R)=12(1−ℓ1f)dφ+12(1−ℓ2f)dϑ.(D13.1)H=\frac{i}{f},\qquad V=C=0,\qquad A^{(R)}=\frac12\left(1-\frac{\ell_1}{f}\right)d\varphi +\frac12\left(1-\frac{\ell_2}{f}\right)d\vartheta . \tag{D13.1}

The space is the unquotiented sphere with its usual spin structure. In particular, the round limit uses the positive auxiliary branch H=i/LH=i/L. Take two dimensionless real masses x,yx,y, meaning physical sources x/L,y/Lx/L,y/L. In the ordered chiral list (A1,A2,B2,B1)(A_1,A_2,B_2,B_1) choose

T1=12diag⁡(1,1,−1,−1),T2=12diag⁡(1,−1,1,−1),Tr⁡4(TaTb)=δab.(D13.2)T_1=\frac12\operatorname{diag}(1,1,-1,-1),\qquad T_2=\frac12\operatorname{diag}(1,-1,1,-1),\qquad \operatorname{Tr}_4(T_aT_b)=\delta_{ab}. \tag{D13.2}

These are the first two Cartan generators of Binder–Chester–Pufu (3.7), after the displayed field ordering. They are flavor generators, not gauge generators. Each superpotential word has total charge zero. Freeze each background multiplet at aa=0a_a=0, σa=(x,y)a/L\sigma_a=(x,y)_a/L and Da=−HσaD_a=-H\sigma_a, with vanishing gauginos. Thus changing a mass also changes its auxiliary source; it is not a change of only the fermion mass.

The charge eigenvalues of xT1+yT2xT_1+yT_2 are p,q,−q,−pp,q,-q,-p, where p=(x+y)/2p=(x+y)/2 and q=(x−y)/2q=(x-y)/2. In the integer-charge cover used in the prerequisite, the one-direction source is M=x/2M=x/2. This fixes the factor of two when comparing derivatives.

Use the inherited double-sine prescription and Abel limit of the Cartan integral, with no added integer source counterterms or new dimensionful subtraction scale. Set Fb=−log⁡ZbF_b=-\log Z_b on the branch continuous from zero real mass. The rank-one common-center integration already established in the prerequisite now gives

Zb(x,y)=Db(p)Db(q).(D13.3)\boxed{Z_b(x,y)=D_b(p)D_b(q).} \tag{D13.3}

Indeed the two matter factors before that integration are Db(s−t+p)D_b(s-t+p) and Db(t−s−q)D_b(t-s-q); the normalized delta sets s−t=0s-t=0, and DbD_b is even. This is Kubo–Nosaka–Pang (4.11) with their ABJM masses (m1,m2,m3)=(x,y,0)(m_1,m_2,m_3)=(x,y,0) and dual variables (ζ,m)=(p,−q)(\zeta,m)=(p,-q). Their dual theory is a three-dimensional N=4\mathcal N=4 gauge theory, not four-dimensional SYM. No large-rank approximation has entered (D13.3).

At the round point D1(z)=1/(2cosh⁡πz)D_1(z)=1/(2\cosh\pi z), so

F1(x,y)=log⁡4+log⁡cosh⁡(πp)+log⁡cosh⁡(πq)=log⁡4+π24(x2+y2)−π496(x4+6x2y2+y4)+O((∣x∣+∣y∣)6).(D13.4)\begin{aligned} F_1(x,y)&=\log4+\log\cosh(\pi p)+\log\cosh(\pi q)\\ &=\log4+\frac{\pi^2}{4}(x^2+y^2) -\frac{\pi^4}{96}(x^4+6x^2y^2+y^4)+O((|x|+|y|)^6). \end{aligned} \tag{D13.4}

The expansion follows from log⁡cosh⁡z=z2/2−z4/12+O(z6)\log\cosh z=z^2/2-z^4/12+O(z^6); the factors 1/21/2 in the generators must remain inside p,qp,q. At x=y=0x=y=0,

Fxx=π22,Fxy=0,Fxxxx=Fxxyy=−π44.(D13.5)F_{xx}=\frac{\pi^2}{2},\qquad F_{xy}=0,\qquad \boxed{F_{xxxx}=F_{xxyy}=-\frac{\pi^4}{4}.} \tag{D13.5}

The boxed result reproduces BCP (3.43). A fourth derivative with respect to the corresponding dimensional mass has an extra factor L4L^4 at fixed radius. A negative fourth derivative does not contradict positivity of a current norm; that positivity constrains the relevant two-point coefficient.

For a normalized generator (D13.2), define τF\tau_F by ⟨jμ(x)jν(0)⟩=τF(δμν∂2−∂μ∂ν)∣x∣−2/(16π2)\langle j_\mu(x)j_\nu(0)\rangle=\tau_F(\delta_{\mu\nu}\partial^2-\partial_\mu\partial_\nu)|x|^{-2}/(16\pi^2) at separated points. Here the four free chirals give τF=∑i(Ta)ii2=1\tau_F=\sum_i(T_a)_{ii}^2=1. In the enhanced theory the sourced Ward normalization is τF=cT/16\tau_F=c_T/16, with cT=1c_T=1 for one real scalar. Equivalently, in this two-mass slice cT=32Re⁡Fxx/π2=16c_T=32\operatorname{Re}F_{xx}/\pi^2=16; Chester–Kalloor–Sharon (2.8) uses m±=y±xm_\pm=y\pm x, accounting for its factor 64. Comparing its stress-tensor tensor structure with the prerequisite gives

τRR=cT16=1,CT=3cT32π2=32π2.(D13.6)\tau_{RR}=\frac{c_T}{16}=1,\qquad C_T=\frac{3c_T}{32\pi^2}=\frac{3}{2\pi^2}. \tag{D13.6}

The equality of these flavor and R coefficients uses this normalized enhanced multiplet. It is not a rule for an arbitrary flavor generator or a generic N=2\mathcal N=2 theory.

Write a one-parameter action as S(x)=S0+xV+x2W/2S(x)=S_0+x\mathcal V+x^2\mathcal W/2. The symbols V=S,x(0)\mathcal V=S_{,x}(0) and W=S,xx(0)\mathcal W=S_{,xx}(0) are integrated operators; the latter is sometimes called a seagull insertion. Expanding −log⁡⟨e−xV−x2W/2⟩0-\log\langle e^{-x\mathcal V-x^2\mathcal W/2}\rangle_0 in connected cumulants gives

Fxxxx=−⟨V4⟩c+6⟨V2W⟩c−3⟨W2⟩c.(D13.7)F_{xxxx}=-\langle\mathcal V^4\rangle_c +6\langle\mathcal V^2\mathcal W\rangle_c -3\langle\mathcal W^2\rangle_c. \tag{D13.7}

For example, the cubic cumulant has three placements of W\mathcal W, each supplying −x2/2-x^2/2; multiplication by 4!4! produces the coefficient six after taking −log⁡Z-\log Z. The formula assumes no cubic or quartic term in this chosen action. A more general source family requires those derivatives too.

Here the conformal scalar operator at b=1b=1, in units L=1L=1, is −∇2+3/4+w2−iw-\nabla^2+3/4+w^2-iw. The term w2w^2 produces W\mathcal W; deleting it changes the theory being differentiated. BCP’s opposite sign for the imaginary linear term is translated by w↦−ww\mapsto-w. Our paired masses make the resulting even derivatives identical.

An independent mode calculation makes the size of the error visible. Put an=n2−1/4a_n=n^2-1/4 for n=1,2,…n=1,2,\ldots. Scalar modes have degeneracy n2n^2; the Dirac eigenvalues are ±i(n+1/2)\pm i(n+1/2) with degeneracy n(n+1)n(n+1) for each sign. A pair of chirals of masses ±w\pm w contributes the following fourth-order coefficient, which is convergent even though the undifferentiated determinant needs its prescription:

[w4]Fpair=∑n=1∞[n2(1an2−(2an+1)22an4)+n(n+1)(n+1/2)4].(D13.8)[w^4]F_{\rm pair}= \sum_{n=1}^\infty\left[ n^2\left(\frac1{a_n^2}-\frac{(2a_n+1)^2}{2a_n^4}\right) +\frac{n(n+1)}{(n+1/2)^4}\right]. \tag{D13.8}

The first bracket comes from log⁡[(an+w2)2+w2]\log[(a_n+w^2)^2+w^2]; the second from −2n(n+1)log⁡[(n+1/2)2+w2]-2n(n+1)\log[(n+1/2)^2+w^2]. After multiplication by 4!4!, partial fractions telescope to −192∑n≥1(2n−1)−4=−2π4-192\sum_{n\geq1}(2n-1)^{-4}=-2\pi^4. Four chirals of masses ±x/2\pm x/2 supply 2/24=1/82/2^4=1/8 of this pair response, recovering (D13.5).

If the scalar w2w^2 term is omitted, its bracket becomes −n2/(2an4)-n^2/(2a_n^4). The difference between the complete and linear-source sums is

248∑n=1∞(−n2an2−2n2an3)=−9π24.(D13.9)\frac{24}{8}\sum_{n=1}^\infty \left(-\frac{n^2}{a_n^2}-\frac{2n^2}{a_n^3}\right) =-\frac{9\pi^2}{4}. \tag{D13.9}

Thus the linear-insertion result is −π4/4+9π2/4-\pi^4/4+9\pi^2/4, not the physical derivative. This reproduces the correction in BCP (3.40)–(3.43). These are source-action insertions, not freely adjustable background counterterms. Nor should the free-theory singlet contribution be carried unchanged into the strongly coupled formula (3.45).

Now keep the same charges, x,yx,y, scale LL and prescription, but vary bb. The prerequisite’s finite determinant difference uses

Kb(t)=sinh⁡[(b+b−1)t/2]sinh⁡(bt)sinh⁡(t/b),log⁡Db(w)Db(0)=∫0∞dttKb(t)[cos⁡(2wt)−1].(D13.10)K_b(t)=\frac{\sinh[(b+b^{-1})t/2]}{\sinh(bt)\sinh(t/b)},\qquad \log\frac{D_b(w)}{D_b(0)}= \int_0^\infty\frac{dt}{t}K_b(t)[\cos(2wt)-1]. \tag{D13.10}

Consequently Fxx(b;0,0)=2∫0∞tKb(t) dtF_{xx}(b;0,0)=2\int_0^\infty tK_b(t)\,dt. Exchange of the two circle scales makes its first bb derivative vanish at one. Differentiating both denominator factors gives

∂b2Kb(t)∣1=−tcoth⁡t+2t2csch⁡2tsinh⁡t,Fbbxx(1;0,0)=2π2−π44.(D13.11)\left.\partial_b^2K_b(t)\right|_1 =\frac{-t\coth t+2t^2\operatorname{csch}^2t}{\sinh t}, \qquad \boxed{F_{bbxx}(1;0,0)=2\pi^2-\frac{\pi^4}{4}.} \tag{D13.11}

To evaluate it, use ∫0∞t/sinh⁡t dt=π2/4\int_0^\infty t/\sinh t\,dt=\pi^2/4 and ∫0∞t3/sinh⁡t dt=π4/8\int_0^\infty t^3/\sinh t\,dt=\pi^4/8, obtained by summing positive exponentials. The identity (csch⁡t)′′=csch⁡t+2csch⁡3t(\operatorname{csch}t)''=\operatorname{csch}t+2\operatorname{csch}^3t and integration by parts then give ∫t3csch⁡3t dt=3π2/4−π4/16\int t^3\operatorname{csch}^3t\,dt=3\pi^2/4-\pi^4/16. All boundary terms vanish. On a compact neighborhood of b=1b=1 the differentiated kernels have integrable small-tt limits and uniform exponential tails, justifying differentiation under the integral.

Equation (D13.11) is also the rank-one specialization of CKS (3.7): Fm+m+=π2/4F_{m_+m_+}=\pi^2/4, Fm+m+m+m+=−π4/8F_{m_+m_+m_+m_+}=-\pi^4/8, and the mixed m+,m−m_+,m_- derivative vanishes. It is a definite integrated response to this complete background family, not an arbitrary metric variation or a separated four-point function by itself.

The executed mode sums and real-kernel integrals agree with the formulas at 70-digit working precision; finite mode cutoffs and finite differentiation steps have separate errors. For example, the uncorrected pair-mode cutoff n≤100n\leq100 still misses the fourth response by about 0.0012040.001204. More working digits do not remove that cutoff error.

A dated reading route and a bounded next calculation

Section titled “A dated reading route and a bounded next calculation”

The literature check for this dossier is dated 20 October 2026. BCP explains how mass derivatives constrain integrated stress-multiplet correlators, including source terms. CKS derives relations between squashing and mass derivatives using a specified complex mass path. KNP studies levels ±1\pm1, with the dual adjoint mass tuned to mb=i(b2−3)/(4b)m_b=i(b^2-3)/(4b), and obtains the perturbative large-rank Airy expansion on its stated locus. Its derivation for odd b2≥3b^2\geq3 and numerical continuation to other b>1b>1 have different evidential status; neither is a generic exact answer at arbitrary level and real mass. The tuned b=1b=1 endpoint instead has mb=−i/2m_b=-i/2 and a divergent determinant, as its footnote 4 explains. Its dual variables and its three-dimensional matter theory must remain in the comparison.

A concrete next assignment is to extract the same fourth derivative at a chosen higher rank from a convergent regulated integral, comparing changes of step size, cutoff and contact prescription separately. The rank-one results above are calibration values, not predictions for that rank. An apparent discrepancy must survive those checks and the source dictionary before it can be interpreted as new CFT information. This selected reading route makes no claim that such a variation is novel or that the literature ends with these papers.

Set x=tcos⁡αx=t\cos\alpha, y=tsin⁡αy=t\sin\alpha at b=1b=1. Show that the generator has unit squared trace. Find the second and fourth tt derivatives at zero, and explain why equal two-point normalization does not force equal fourth response.

Hint 1

Use the orthogonality of the two generators, then substitute into the polynomial in (D13.4).

Hint 2

The fourth invariant contains 6cos⁡2αsin⁡2α6\cos^2\alpha\sin^2\alpha, whereas the square of the quadratic invariant contains coefficient two.

Solution

For T=cos⁡α T1+sin⁡α T2T=\cos\alpha\,T_1+\sin\alpha\,T_2, orthogonality gives Tr⁡T2=1\operatorname{Tr}T^2=1. Thus

Ftt(0)=π22,Ftttt(0)=−π44[1+sin⁡2(2α)].F_{tt}(0)=\frac{\pi^2}{2},\qquad F_{tttt}(0)=-\frac{\pi^4}{4} \bigl[1+\sin^2(2\alpha)\bigr].

Indeed Tr⁡T4=[1+sin⁡2(2α)]/4\operatorname{Tr}T^4=[1+\sin^2(2\alpha)]/4 is not fixed by Tr⁡T2\operatorname{Tr}T^2. At α=0\alpha=0 the fourth response is −π4/4-\pi^4/4; at α=π/4\alpha=\pi/4 it is −π4/2-\pi^4/2. The latter direction has eigenvalues 1/2,0,0,−1/21/\sqrt2,0,0,-1/\sqrt2. These different eigenvalue sets are not conjugate flavor generators. The comparison is a local source-direction calculation; a generic irrational rotation is not being declared a new compact unit-period subgroup.

For f(x)=F1(x,0)f(x)=F_1(x,0), define

Dh=f(2h)−4f(h)+6f(0)−4f(−h)+f(−2h)h4.\mathscr D_h= \frac{f(2h)-4f(h)+6f(0)-4f(-h)+f(-2h)}{h^4}.

Find its leading bias, construct a fourth-order improved estimate, and bound amplification of an absolute input uncertainty δ\delta in each value of ff. Explain whether taking hh as small as possible is always sensible.

Hint 1

Expand through x6x^6. The stencil annihilates the constant and quadratic terms, but not the sixth-order term.

Hint 2

Combine the estimates at hh and h/2h/2 to cancel a term proportional to h2h^2. For the input bound, add absolute stencil coefficients.

Solution

Taylor expansion gives Dh=f(4)(0)+h2f(6)(0)/6+O(h4)\mathscr D_h=f^{(4)}(0)+h^2f^{(6)}(0)/6+O(h^4). Since f=log⁡4+2log⁡cosh⁡(πx/2)f=\log4+2\log\cosh(\pi x/2), f(6)(0)=π6/2f^{(6)}(0)=\pi^6/2, and

Dh=−π44+π612h2+O(h4),Rh=4Dh/2−Dh3=f(4)(0)+O(h4).\mathscr D_h=-\frac{\pi^4}{4}+\frac{\pi^6}{12}h^2+O(h^4), \qquad \mathscr R_h=\frac{4\mathscr D_{h/2}-\mathscr D_h}{3} =f^{(4)}(0)+O(h^4).

For h=0.08,0.04,0.02h=0.08,0.04,0.02, the direct biases are approximately 0.502627,0.127543,0.03200600.502627,0.127543,0.0320060; the improved biases are 0.00251550,0.000160266,0.00001006530.00251550,0.000160266,0.0000100653. They exhibit the predicted orders, not arbitrary precision. Input uncertainty contributes at most 16δ/h416\delta/h^4 to the direct estimate. Decreasing hh suppresses truncation bias while amplifying numerical or integral error. One must vary both the step and the input accuracy, retaining the source prescription throughout.

Use KNP’s rank-one dual variables, for which F(b;ζ,m)=−log⁡[Db(ζ)Db(m)]\mathcal F(b;\zeta,m)=-\log[D_b(\zeta)D_b(m)]. Set ζ=0\zeta=0 and m=mb=i(b2−3)/(4b)m=m_b=i(b^2-3)/(4b). At b0=3b_0=\sqrt3 the mass vanishes. Given D3(m)=[1+2cosh⁡(2πm/3)]−1D_{\sqrt3}(m)=[1+2\cosh(2\pi m/\sqrt3)]^{-1}, compare d2F(b;0,mb)/db2d^2\mathcal F(b;0,m_b)/db^2 with the fixed-mass partial derivative Fbb\mathcal F_{bb} at that point.

Hint 1

Differentiate the composite function twice before setting m=0m=0. Evenness in mm removes two terms, but not the term containing (mb′)2(m_b')^2.

Hint 2

At b0b_0, mb′=i/2m_b'=i/2. Compute Fmm\mathcal F_{mm} from the supplied elementary determinant, retaining i2=−1i^2=-1.

Solution

The chain rule gives

d2Fdb2=Fbb+2mb′Fbm+(mb′)2Fmm+mb′′Fm.\frac{d^2\mathcal F}{db^2} =\mathcal F_{bb}+2m_b'\mathcal F_{bm} +(m_b')^2\mathcal F_{mm}+m_b''\mathcal F_m.

At zero mass, Fm=Fbm=0\mathcal F_m=\mathcal F_{bm}=0 by evenness. The elementary formula gives Fmm=8π2/9\mathcal F_{mm}=8\pi^2/9, so

d2F(b;0,mb)db2∣b0=Fbb∣b0,m=0−2π29.\left.\frac{d^2\mathcal F(b;0,m_b)}{db^2}\right|_{b_0} =\left.\mathcal F_{bb}\right|_{b_0,m=0}-\frac{2\pi^2}{9}.

The negative addition comes from the imaginary tangent to the mass path. Near this point the determinant has no zero or pole, so the local analytic continuation is well defined; this does not authorize an arbitrary global contour continuation. The two derivatives probe different source families even though they pass through the same background. Neither should be identified with the round-point response (D13.11).

Throughout, allowed quadratic flavor Chern–Simons contact changes shift the complex mass Hessian but not these fourth mass derivatives. For example, the integer-charge cover gives ΔF=iπκ(x/2)2\Delta F=i\pi\kappa(x/2)^2 with κ∈Z\kappa\in\mathbb Z. Its sign is fixed by the positive-HH background, and its fourth derivative vanishes. It is independent of bb, so it also leaves FbbxxF_{bbxx} unchanged. This statement presumes the declared conformal prescription; introducing arbitrary new dimensionful source counterterms would change the observable under comparison.