Skip to content

From Blocks to Full Connection Coefficients

A classical conformal block can determine the accessory parameter of an ODE without determining a single normalized connection coefficient. The missing information is not cosmetic: it specifies which two solutions are being compared, how their leading terms are normalized, which continuation path reaches the target, and which finite fusion or sectorial kernel acts on the degenerate probe.

This page turns that list into one matrix identity. It then explains why derivatives of the classical block survive the b0b\to0 limit, why an order-one background term usually cancels from the leading unit-normalized ODE matrix, and when that term must be restored. The discussion stays in the classical cVirc_{\mathrm{Vir}}\to\infty regime; the cVir=1c_{\mathrm{Vir}}=1 Fourier sum of page 2 is not an extra factor in the matrix below.

Four outputs that should not share one name

Section titled “Four outputs that should not share one name”

The word “coefficient” is used for several mathematically different objects. The following ladder keeps their data requirements separate.

OutputWhat determines itWhat it still does not determine
Oper accessoryA modulus derivative such as tf\partial_t f or LLFL\partial_LF on a chosen block branchEndpoint frames or a path
Block connection kernelDegenerate fusion, braiding, collision, or Whittaker identitiesUnit-leading ODE normalization
ODE connection matrixKernel, accessory inversion, two endpoint conversions, and a path or sectorA physical pairing or flux convention
Physical coefficientODE matrix plus vertex, reflection, pairing, boundary, or flux normalizationNothing beyond the declared physical problem

The accessory specifies the differential equation. A connection matrix specifies the same two-dimensional solution space together with two ordered bases. A physical amplitude contracts that matrix with additional boundary data. These operations are related, but none is a synonym for another.

The frame equation is the complete assembly rule

Section titled “The frame equation is the complete assembly rule”

Work first away from resonance. Let

Bi=(Bi,1,Bi,2),Bj=(Bj,1,Bj,2)\boldsymbol{\mathcal B}_i = \left( \mathcal B_{i,1}, \mathcal B_{i,2} \right), \qquad \boldsymbol{\mathcal B}_j = \left( \mathcal B_{j,1}, \mathcal B_{j,2} \right)

be two row frames of semiclassical block representatives. Suppose the block continuation is indexed source first:

Bi,s=rKsr[p]Bj,r,\mathcal B_{i,s} = \sum_r \mathcal K_{sr}^{[\mathfrak p]} \mathcal B_{j,r},

where [p][\mathfrak p] records the continuation path and, for an irregular endpoint, the lateral sector. The target-row matrix is

(Kji[p])rs:=Ksr[p].\left( K_{ji}^{[\mathfrak p]} \right)_{rs} := \mathcal K_{sr}^{[\mathfrak p]}.

Thus

Bi=BjKji[p].\boldsymbol{\mathcal B}_i = \boldsymbol{\mathcal B}_j K_{ji}^{[\mathfrak p]}.

This is the transpose that is easily lost when a paper lists coefficients with the source sign first.

Let P(z)P(z) be the common scalar gauge from the normal-form BPZ solution to the chosen standard ODE. Declare the ODE frames by

Hi=P(z)BiNi,Hj=P(z)BjNj.\begin{aligned} \boldsymbol H_i &= P(z) \boldsymbol{\mathcal B}_i \mathscr N_i, \\ \boldsymbol H_j &= P(z) \boldsymbol{\mathcal B}_j \mathscr N_j. \end{aligned}

The matrices Ni\mathscr N_i and Nj\mathscr N_j are independent of zz. For generic Frobenius frames they are diagonal. They may contain:

  • the leading coefficient of the scalar gauge at a finite endpoint;
  • a power of a collision or modulus parameter;
  • a finite derivative of the classical block;
  • a phase from the chosen local coordinate;
  • a relative sign fixing the declared frame orientation.

At an irregular endpoint, Nj\mathscr N_j is defined by an asymptotic comparison inside a sector. There is no meaningful instruction to “evaluate P()P(\infty).”

Substitution gives

Hi=PBjKji[p]Ni=HjNj1Kji[p]Ni.\begin{aligned} \boldsymbol H_i &= P\boldsymbol{\mathcal B}_j K_{ji}^{[\mathfrak p]} \mathscr N_i \\ &= \boldsymbol H_j \mathscr N_j^{-1} K_{ji}^{[\mathfrak p]} \mathscr N_i. \end{aligned}

Therefore

Cji[p]=Nj1Kji[p]Ni.\boxed{ C_{ji}^{[\mathfrak p]} = \mathscr N_j^{-1} K_{ji}^{[\mathfrak p]} \mathscr N_i. }

The common scalar function P(z)P(z) cancels. Its endpoint constants do not cancel after they have been assigned to different unit-leading frames; those constants live in Ni\mathscr N_i and Nj\mathscr N_j.

Taking determinants gives

detCji[p]=detNidetNjdetKji[p].\det C_{ji}^{[\mathfrak p]} = \frac{ \det\mathscr N_i }{ \det\mathscr N_j } \det K_{ji}^{[\mathfrak p]}.

For any ordinary matching point zz_*,

Cji[p]=Φj(z)1Φi(z),Φ=(H,1H,2H,1H,2).C_{ji}^{[\mathfrak p]} = \Phi_j(z_*)^{-1}\Phi_i(z_*), \qquad \Phi_\ell = \begin{pmatrix} H_{\ell,1}&H_{\ell,2}\\ H_{\ell,1}'&H_{\ell,2}' \end{pmatrix}.

Hence the same determinant must equal

detCji[p]=Wr[Hi,1,Hi,2]Wr[Hj,1,Hj,2].\det C_{ji}^{[\mathfrak p]} = \frac{ \operatorname{Wr} \left[ H_{i,1},H_{i,2} \right] }{ \operatorname{Wr} \left[ H_{j,1},H_{j,2} \right] }.

This check is independent of conformal blocks. A wrong transpose can occasionally preserve the determinant, but a missing endpoint power, a single-column sign error, or a wrong gamma numerator usually does not.

If the declared endpoint frames are changed by constant invertible matrices,

H~i=HiMi,H~j=HjMj,\widetilde{\boldsymbol H}_i = \boldsymbol H_iM_i, \qquad \widetilde{\boldsymbol H}_j = \boldsymbol H_jM_j,

then

C~ji=Mj1CjiMi.\widetilde C_{ji} = M_j^{-1}C_{ji}M_i.

Thus an individual entry is not basis invariant. Zero conditions become meaningful only after the source and target boundary vectors have been transformed with the frames. This is why “the connection coefficient” is incomplete unless its two bases are named.

Every factor has one provenance and one audit

Section titled “Every factor has one provenance and one audit”

The assembly can be implemented as a data contract.

FactorOriginIndependent audit
Internal lift a(q)a(q) or m(q)m'(q)Invert the accessory relation on a declared branchSubstitute back; monitor the inverse Jacobian
KjiK_{ji}Exact degenerate fusion/braiding or sectorial kernel, followed by its classical limitRigid special-function reduction or kernel identities
External derivativesVanishing degenerate shifts of heavy external momentaFinite-difference test before taking b0b\to0
Internal derivativeShifted intermediate channel in a composed continuationSymmetry under the two intermediate signs
Coordinate powersUnit-leading Frobenius or asymptotic conventionDirect local series
Gauge constantsLeading scalar gauge at each finite endpointSubstitute the local gauge asymptotic
Path phaseDeclared logarithms and continuation pathRepeat with the opposite lateral path
Sector sign or Stokes actionOrientation of the sectorial target frameWronskian and adjacent-sector law
Whole matrixAll preceding factorsTwo-point matching and Abel’s determinant

The inverse step is part of the answer. Derivatives of FF are first taken at fixed internal lift and fixed remaining parameters; only then is the selected a=a(q)a=a(q) substituted. Differentiating the composite F(a(q),)F(a(q),\ldots) would add a chain-rule term that does not arise from the degenerate momentum shift.

The accessory branch, block continuation, and endpoint ledgers feed a target-row connection-matrix assembly, which is checked by Abel's identity, direct matching, and frame covariance.

The accessory inverse selects the internal branch used by the exact kernel, while the two endpoint ledgers supply powers, gauges, derivative dressings, phases, and frame signs. Abel’s determinant is a necessary audit; matching all four entries at ordinary points completes the check.

A vanishing shift leaves a finite derivative

Section titled “A vanishing shift leaves a finite derivative”

The derivative dressings are sometimes described as “subleading,” but they are order one in the normalized probe. Let λ=bαL\lambda=b\alpha_{\rm L} be any scaled heavy momentum. Fusion with a light (2,1)(2,1) degenerate field shifts it by

λθ=λθb22,θ=±1.\lambda_\theta = \lambda-\frac{\theta b^2}{2}, \qquad \theta=\pm1.

Assume a differentiable even-power expansion of the relevant background block,

\LogFbg(λ;b)=fcl(λ)b2+gbg(λ)+b2hbg(λ)+O(b4).\begin{aligned} \Log\mathfrak F_{\mathrm{bg}}(\lambda;b) ={}& \frac{ f_{\mathrm{cl}}(\lambda) }{b^2} + g_{\mathrm{bg}}(\lambda) \\ &+ b^2h_{\mathrm{bg}}(\lambda) + O(b^4). \end{aligned}

Taylor expansion before taking the limit gives

\LogFbg(λθ;b)Fbg(λ;b)=θ2λfcl+b2[18λ2fclθ2λgbg]+O(b4).\begin{aligned} \Log \frac{ \mathfrak F_{\mathrm{bg}}(\lambda_\theta;b) }{ \mathfrak F_{\mathrm{bg}}(\lambda;b) } ={}& -\frac{\theta}{2} \partial_\lambda f_{\mathrm{cl}} \\ &+ b^2 \left[ \frac18 \partial_\lambda^2f_{\mathrm{cl}} - \frac{\theta}{2} \partial_\lambda g_{\mathrm{bg}} \right] \\ &+ O(b^4). \end{aligned}

Consequently,

limb0Fbg(λθ;b)Fbg(λ;b)=exp[θ2λfcl].\lim_{b\to0} \frac{ \mathfrak F_{\mathrm{bg}}(\lambda_\theta;b) }{ \mathfrak F_{\mathrm{bg}}(\lambda;b) } = \exp \left[ -\frac{\theta}{2} \partial_\lambda f_{\mathrm{cl}} \right].

The shift vanishes, but it acts on an exponent of order b2b^{-2}. This is the origin of factors such as

esF0/2,esF1/2,eσFa/2,erFm/2.e^{-sF_0/2}, \qquad e^{-sF_1/2}, \qquad e^{-\sigma F_a/2}, \qquad e^{-rF_m/2}.

They are finite normalization data, not changes in the leading local exponents.

Three derivatives have different jobs:

DerivativeRole
tf\partial_t f or LLFL\partial_LFFixes the oper accessory
aiF\partial_{a_i}FNormalizes a degenerate branch at endpoint ii
aF\partial_aF, mF\partial_mF, or mFD\partial_{m'}F_DNormalizes an intermediate or irregular channel

The modulus derivative and the momentum derivatives cannot be exchanged. They differentiate different data while holding different variables fixed.

Put

Rθ(b):=Fbg(λθ;b)Fbg(λ;b).R_\theta(b) := \frac{ \mathfrak F_{\mathrm{bg}}(\lambda_\theta;b) }{ \mathfrak F_{\mathrm{bg}}(\lambda;b) }.

The two sign combinations obey

\LogR+\LogR=λfclb2λgbg+O(b4),\LogR++\LogR=b24λ2fcl+O(b4).\begin{aligned} \Log R_+ - \Log R_- ={}& -\partial_\lambda f_{\mathrm{cl}} - b^2\partial_\lambda g_{\mathrm{bg}} + O(b^4), \\ \Log R_+ + \Log R_- ={}& \frac{b^2}{4} \partial_\lambda^2f_{\mathrm{cl}} + O(b^4). \end{aligned}

These formulas give a clean finite-bb diagnostic. After subtracting the leading derivative, the residual should scale by a factor of four under bb/2b\mapsto b/2. The antisymmetric residual measures λgbg\partial_\lambda g_{\mathrm{bg}}; the symmetric residual measures the Hessian of the classical action.

For a generic, nonzero entry of a unit-leading classical ODE matrix, gbgg_{\mathrm{bg}} does not contribute at order b0b^0 to the shifted ratio above. Its difference begins at order b2b^2. The leading matrix therefore needs:

  • the complete classical action in the chosen normalization;
  • its relevant modulus and momentum derivatives;
  • the leading classical limit of the exact degenerate kernel;
  • exact local-coordinate, gauge, path, and sector data.

To compute the first b2b^2 correction, one must add all contributions at that order, not only one convenient term:

first correction=Hessian of fcl+gbg+finite-b kernel correction+corrected probe and gauge normalization.\begin{aligned} \text{first correction} = {}& \text{Hessian of }f_{\mathrm{cl}} \\ &+ \partial g_{\mathrm{bg}} + \text{finite-}b\text{ kernel correction} \\ &+ \text{corrected probe and gauge normalization}. \end{aligned}

Absolute chiral or physical coefficients can contain order-one three-point, reflection, pairing, or flux factors even when gbgg_{\mathrm{bg}} cancels from the unit-leading ODE ratio. A saddle evaluation of an integral fusion transform can also contribute a Hessian determinant. Such quantities belong to the declared chiral or physical observable, not automatically to the unit-leading scalar ODE matrix.

If a leading entry vanishes, two exponents coalesce, or a generic Verma-module inverse develops a pole, the power counting can change and subleading terms may become decisive. Page 7 treats those singular limits rather than substituting them into the generic formulas here.

A regular target is one diagonal conversion on each side

Section titled “A regular target is one diagonal conversion on each side”

The general-Heun endpoint assembly is the regular-singular prototype. In its notation,

H0=HtCt0,Ct0=Nt1MfrN0.\boldsymbol H_0 = \boldsymbol H_tC_{t0}, \qquad C_{t0} = \mathscr N_t^{-1} M_{\mathrm{fr}} \mathscr N_0.

Here (Mfr)rs=Msr(M_{\mathrm{fr}})_{rs}=\mathcal M_{sr} is the transpose of the source-first degenerate fusion array. The matrices N0,Nt\mathscr N_0,\mathscr N_t contain local powers, scalar-gauge constants, and the external derivatives of the classical block. Expanding this product gives all four gamma-function entries in the general-Heun fusion core.

This compact form also says what the gamma kernel cannot know. It does not know which Frobenius coordinate was assigned leading coefficient one, which logarithm defines a fractional power, or which tt-independent normalization was added to the classical block. Those choices enter only through the two N\mathscr N matrices.

The sectorial case makes the normalization budget more visible. Use the standard confluent-Heun equation

y+(γz+δz1+L)y+αzqz(z1)y=0y'' + \left( \frac{\gamma}{z} + \frac{\delta}{z-1} + L \right)y' + \frac{\alpha z-q}{z(z-1)}y =0

and the common normal-form gauge

P(z)=eLz/2zγ/2(1z)δ/2.P(z) = e^{-Lz/2} z^{-\gamma/2} (1-z)^{-\delta/2}.

Put

a0=1γ2,a1=1δ2,a_0=\frac{1-\gamma}{2}, \qquad a_1=\frac{1-\delta}{2},

and keep the full accessory dictionary, logarithms, and frame definitions of the confluent-Heun gauge ledger. The ordered regular frames are

H0=(H0,,H0,+),H1=(H1,,H1,+),\boldsymbol H_0 = \left( H_{0,-},H_{0,+} \right), \qquad \boldsymbol H_1 = \left( H_{1,-},H_{1,+} \right),

while the sectorial infinity frame is

Y(η)=(Y+(η),Y(η)).\boldsymbol Y_\infty^{(\eta)} = \left( Y_+^{(\eta)},Y_-^{(\eta)} \right).

For a lateral ray on which

arg(1z)=ηπ,η=±1,\arg(1-z)=-\eta\pi, \qquad \eta=\pm1,

write

χη:=eηiπδ/2.\chi_\eta := e^{\eta\ii\pi\delta/2}.

The square χη2\chi_\eta^2 is the phase that enters the Wronskian ratio.

The sign and chart conventions used in the remaining formulas are:

SymbolMeaning and order
H0,H1\boldsymbol H_0,\boldsymbol H_1Regular frames ordered (,+)(-,+)
Y\boldsymbol Y_\inftyInfinity frame ordered (+,)(+,-)
ss in a source-first arraySource sign
rr in a source-first arrayTarget sign
KjiK_{ji}Target-row matrix: target index first, source index second
η\etaChoice of the two displayed lateral determinations
kkGeneral sector index under repeated continuation
mmFixed irregular/background momentum
mm'Internal coordinate of the dual large-LL chart
μ,ρ\mu,\rhoAbbreviations mmm'-m and 2mm2m'-m

Every block derivative is taken at fixed remaining arguments and fixed logarithms before the accessory inverse is substituted. With

A(z):=eLzzγ(1z)δ,\mathcal A(z) := e^{-Lz} z^{-\gamma} (1-z)^{-\delta},

the three declared unit-leading frames satisfy

Wr[H0,,H0,+]=2a0A(z),Wr[H1,,H1,+]=2a1eLA(z),Wr[Y+(η),Y(η)]=Lχη2A(z).\begin{aligned} \operatorname{Wr} \left[ H_{0,-},H_{0,+} \right] &= 2a_0\mathcal A(z), \\ \operatorname{Wr} \left[ H_{1,-},H_{1,+} \right] &= -2a_1e^L\mathcal A(z), \\ \operatorname{Wr} \left[ Y_+^{(\eta)},Y_-^{(\eta)} \right] &= -L\chi_\eta^{-2}\mathcal A(z). \end{aligned}

These exact Abel normalizations are the target that both asymptotic block charts must reproduce.

The block entries in this subsection are all-orders formal small-LL identities, conditional on the derivative-compatible semiclassical block limit stated on the preceding page’s scope statement. The unit-leading endpoint conversions and the Abel Wronskian laws are exact for the declared ODE frames.

Let

Fν:=νF,ν{a0,a1,a,m},F_\nu:=\partial_\nu F, \qquad \nu\in \left\{ a_0,a_1,a,m \right\},

with the abbreviations F0,F1,Fa,FmF_0,F_1,F_a,F_m used below. The block-to-ODE conversions forced by the unit-leading local terms are

H0=PB0N0,H1=PB1N1,Y(η)=PD(η)N(η),\begin{aligned} \boldsymbol H_0 &= P\boldsymbol B_0 \mathscr N_0, \\ \boldsymbol H_1 &= P\boldsymbol B_1 \mathscr N_1, \\ \boldsymbol Y_\infty^{(\eta)} &= P\boldsymbol D_\infty^{(\eta)} \mathscr N_\infty^{(\eta)}, \end{aligned}

where

N0=(eF0/200eF0/2),\mathscr N_0 = \begin{pmatrix} e^{-F_0/2}&0\\ 0&e^{F_0/2} \end{pmatrix}, N1=eL/2(eF1/200eF1/2),\mathscr N_1 = e^{L/2} \begin{pmatrix} e^{-F_1/2}&0\\ 0&e^{F_1/2} \end{pmatrix},

and, in the target order (+,)(+,-),

N(η)=χη1×(eFm/2L1/2+m00eFm/2L1/2m).\begin{aligned} \mathscr N_\infty^{(\eta)} ={}& \chi_\eta^{-1} \\ &\times \begin{pmatrix} e^{F_m/2}L^{1/2+m}&0\\ 0&e^{-F_m/2}L^{1/2-m} \end{pmatrix}. \end{aligned}

All powers use the same selected \LogL\Log L. The common factor eL/2e^{L/2} in N1\mathscr N_1 is the inverse of the gauge value eL/2e^{-L/2} at z=1z=1. It is not generated by a conformal block and cannot be inferred from the fusion kernel.

Let

(Mfr)rs=Msr(a0,a1;a),s,r=,+,\left( M_{\mathrm{fr}} \right)_{rs} = \mathcal M_{sr}(a_0,a_1;a), \qquad s,r=-,+,

with M\mathcal M defined in the first-kind block connection. Within the formal small-LL chart, the complete matrix is

C10small=N11MfrN0.\boxed{ C_{10}^{\mathrm{small}} = \mathscr N_1^{-1} M_{\mathrm{fr}} \mathscr N_0. }

Equivalently, entrywise,

(C10small)rs=eL/2exp[sF0rF12]Msr(a0,a1;a).\left( C_{10}^{\mathrm{small}} \right)_{rs} = e^{-L/2} \exp \left[ \frac{ sF_0-rF_1 }{2} \right] \mathcal M_{sr}(a_0,a_1;a).

The signs s,rs,r are read as the numbers 1,+1-1,+1. The accessory equation is inverted for a=a(q)a=a(q) before the entries are evaluated, but after the partial derivatives have been taken.

Gamma reflection gives

detMfr=a0a1.\det M_{\mathrm{fr}} = -\frac{a_0}{a_1}.

Since

detN0=1,detN1=eL,\det\mathscr N_0=1, \qquad \det\mathscr N_1=e^L,

the exact Abel determinant of the declared ODE frames is

detC10small=a0a1eL.\det C_{10}^{\mathrm{small}} = -\frac{a_0}{a_1}e^{-L}.

Abel’s identity gives the same result from the unit-leading zero- and one-frames. In particular, the missing factor eL/2e^{-L/2} in each entry would produce a determinant error by eLe^L.

To keep the target-row convention unambiguous, denote the source-first collision array in the small-LL overlap factorization by Krawsf,(η)\mathcal K_{\mathrm{raw}}^{\mathrm{sf},(\eta)}. It is

(Krawsf,(η))sr=σ=±Msσ(a1,a;a0)Nσ,r(η)(a,m)×LσaeσFa/2.\begin{aligned} \left( \mathcal K_{\mathrm{raw}}^{\mathrm{sf},(\eta)} \right)_{sr} ={}& \sum_{\sigma=\pm} \mathcal M_{s\sigma}(a_1,a;a_0) \mathcal N_{-\sigma,r}^{(\eta)}(a,m) \\ &\times L^{\sigma a} e^{-\sigma F_a/2}. \end{aligned}

Its identification with the declared second-kind frame requires

Ksf,(η)=Krawsf,(η)J,J=diag(1,1).\mathcal K^{\mathrm{sf},(\eta)} = \mathcal K_{\mathrm{raw}}^{\mathrm{sf},(\eta)}J_\infty, \qquad J_\infty = \operatorname{diag}(1,-1).

Equivalently, set

j+=1,j=1,j_+=1, \qquad j_-=-1,

so that

Ksrsf,(η)=jr(Krawsf,(η))sr.\mathcal K_{sr}^{\mathrm{sf},(\eta)} = j_r \left( \mathcal K_{\mathrm{raw}}^{\mathrm{sf},(\eta)} \right)_{sr}.

The corresponding target-row matrix is

K1(η):=(Ksf,(η))T.K_{\infty1}^{(\eta)} := \left( \mathcal K^{\mathrm{sf},(\eta)} \right)^{\mathsf T}.

With

H1=Y(η)C1(η),\boldsymbol H_1 = \boldsymbol Y_\infty^{(\eta)} C_{\infty1}^{(\eta)},

the normalized matrix is

C1(η)=(N(η))1K1(η)N1.\boxed{ C_{\infty1}^{(\eta)} = \left( \mathscr N_\infty^{(\eta)} \right)^{-1} K_{\infty1}^{(\eta)} \mathscr N_1. }

Its individual target-row entries are

(C1(η))rs=χηeL/2jrL1/2rm×exp[sF1rFm2]×σ=±Msσ(a1,a;a0)Nσ,r(η)(a,m)×LσaeσFa/2.\begin{aligned} \left( C_{\infty1}^{(\eta)} \right)_{rs} ={}& \chi_\eta e^{L/2} j_r L^{-1/2-rm} \\ &\times \exp \left[ \frac{ sF_1-rF_m }{2} \right] \\ &\times \sum_{\sigma=\pm} \mathcal M_{s\sigma}(a_1,a;a_0) \mathcal N_{-\sigma,r}^{(\eta)}(a,m) \\ &\hspace{6em}\times L^{\sigma a} e^{-\sigma F_a/2}. \end{aligned}

The formula has four logically distinct layers: the MN\mathcal M \mathcal N kernel, the internal-channel dressing, the source normalization, and the sectorial target normalization. The sign jrj_r belongs to the last identification; it does not alter the exact Whittaker identity defining N\mathcal N.

The source-first and target-row matrices have the same determinant:

detKsf,(η)=detK1(η)=2a1.\det\mathcal K^{\mathrm{sf},(\eta)} = \det K_{\infty1}^{(\eta)} = 2a_1.

Moreover,

detN1=eL,detN(η)=χη2L.\det\mathscr N_1=e^L, \qquad \det\mathscr N_\infty^{(\eta)} = \chi_\eta^{-2}L.

Therefore the formal block matrix reproduces the exact Abel determinant

detC1(η)=2a1LeLχη2.\det C_{\infty1}^{(\eta)} = \frac{ 2a_1 }{L} e^L \chi_\eta^2.

The entrywise block representation is an all-orders formal small-LL identity. The determinant statement is exact for the declared ODE frames and is reproduced directly by their Abel Wronskians. Without JJ_\infty, the sign would be wrong; without eL/2e^{L/2}, the exponential factor would be wrong.

Composition supplies a third audit:

C0(η)=C1(η)C10small,detC0(η)=2a0Lχη2.C_{\infty0}^{(\eta)} = C_{\infty1}^{(\eta)} C_{10}^{\mathrm{small}}, \qquad \det C_{\infty0}^{(\eta)} = -\frac{2a_0}{L}\chi_\eta^2.

The factors e±Le^{\pm L} cancel because the source has moved from z=1z=1 to z=0z=0.

The block entries in this subsection are formal large-LL asymptotic identities under the corresponding dual-block exponentiation hypothesis. The endpoint conversions follow exactly from the stated leading terms, and the Abel determinant laws of the actual ODE frames remain exact. No equality of resummed small- and large-LL charts is assumed.

The large-LL chart uses a different internal coordinate mm' and the dual block FDF_D. Write

FD,0=a0FD,FD,1=a1FD,FD,m=mFD,FD,m=mFD.F_{D,0} = \partial_{a_0}F_D, \quad F_{D,1} = \partial_{a_1}F_D, \quad F_{D,m} = \partial_mF_D, \quad F_{D,m'} = \partial_{m'}F_D.

For compact formulas also set

μ:=mm,ρ:=2mm.\mu:=m'-m, \qquad \rho:=2m'-m.

The unit-leading conversions now contain powers absent from the small-LL chart:

N0D=(La0eFD,0/200La0eFD,0/2),\mathscr N_0^{D} = \begin{pmatrix} L^{a_0}e^{-F_{D,0}/2}&0\\ 0&L^{-a_0}e^{F_{D,0}/2} \end{pmatrix}, N1D=eL/2(La1eFD,1/200La1eFD,1/2),\mathscr N_1^{D} = e^{L/2} \begin{pmatrix} L^{a_1}e^{-F_{D,1}/2}&0\\ 0&L^{-a_1}e^{F_{D,1}/2} \end{pmatrix},

and

ND,(η)=χη1×(eL/2eFD,m/2L1/2(mm)00eL/2eFD,m/2L1/2+(mm)).\begin{aligned} \mathscr N_\infty^{D,(\eta)} ={}& \chi_\eta^{-1} \\ &\times \begin{pmatrix} e^{L/2} e^{F_{D,m}/2} L^{1/2-(m'-m)} &0 \\ 0& e^{-L/2} e^{-F_{D,m}/2} L^{1/2+(m'-m)} \end{pmatrix}. \end{aligned}

These matrices follow from the displayed one-end and infinity asymptotics in the large-LL block chart, together with the zero-end analogue obtained by endpoint exchange. They also show why replacing FF by FDF_D in a small-LL component formula is insufficient: the powers of LL and the endpoint exponentials change at the same time.

Define the source-first block kernel

(QD,sf)sr=σ=±Γ(1+2sa0)Γ(2ra1)Γ(12+sa0σm)Γ(12ra1σμ)×eσL/2LσρeσFD,m/2×exp[iπ2(1σ)(sa0ra1ρ)].\begin{aligned} \left( \mathcal Q^{D,\mathrm{sf}} \right)_{sr} ={}& \sum_{\sigma=\pm} \frac{ \Gamma(1+2s a_0) \Gamma(-2r a_1) }{ \Gamma\left( \frac12+s a_0-\sigma m' \right) \Gamma\left( \frac12-r a_1-\sigma\mu \right) } \\ &\times e^{\sigma L/2} L^{-\sigma\rho} e^{-\sigma F_{D,m'}/2} \\ &\times \exp \left[ \frac{\ii\pi}{2} (1-\sigma) \left( s a_0-r a_1-\rho \right) \right]. \end{aligned}

This is the expanded product of a Whittaker kernel and an inverse Whittaker kernel. Writing it explicitly removes any ambiguity about the row order of that inverse. The displayed phase uses the same clockwise Whittaker branch and continuation homotopy, which we denote [p+][\mathfrak p_+]; C10DC_{10}^{D} is shorthand for C10D,[p+]C_{10}^{D,[\mathfrak p_+]}. Other homotopy classes require the corresponding monodromy transport. Define the target-row transpose by

K10D:=(QD,sf)T.K_{10}^{D} := \left( \mathcal Q^{D,\mathrm{sf}} \right)^{\mathsf T}.

The normalized matrix is

C10D=(N1D)1K10DN0D,C_{10}^{D} = \left( \mathscr N_1^{D} \right)^{-1} K_{10}^{D} \mathscr N_0^{D},

or, entrywise,

(C10D)rs=eL/2Lsa0+ra1exp[sFD,0rFD,12](QD,sf)sr.\boxed{ \left( C_{10}^{D} \right)_{rs} = e^{-L/2} L^{-s a_0+r a_1} \exp \left[ \frac{ sF_{D,0}-rF_{D,1} }{2} \right] \left( \mathcal Q^{D,\mathrm{sf}} \right)_{sr}. }

The exact kernel identity gives

detQD,sf=detK10D=a0a1.\det\mathcal Q^{D,\mathrm{sf}} = \det K_{10}^{D} = -\frac{a_0}{a_1}.

Therefore the large-LL formal matrix obeys the same ODE determinant law as the small-LL matrix:

detC10D=a0a1eL.\det C_{10}^{D} = -\frac{a_0}{a_1}e^{-L}.

Agreement of the determinant does not identify the two asymptotic series entry by entry. It only confirms that both charts use the same unit-leading ODE frames.

Regular to infinity in one reference sector

Section titled “Regular to infinity in one reference sector”

The printed large-LL block relation fixes the clockwise Whittaker determination

\Log ⁣(eiπx)=\Logxiπ.\Log\!\left(e^{-\ii\pi}x\right) = \Log x-\ii\pi.

In the convention used here this is the reference lateral ray arg(1z)=π\arg(1-z)=-\pi, hence η=+1\eta=+1. Its rigid source-first kernel is

(K1D,sf,(+))sr=Ns,r(+)(a1,μ).\left( \mathcal K_{\infty1}^{D,\mathrm{sf},(+)} \right)_{sr} = \mathcal N_{s,-r}^{(+)} \left( a_1,\mu \right).

Explicitly,

Ns,r(+)(a1,μ)=Γ(1+2sa1)Γ(12+sa1+rμ)×exp[iπ2(1+r)(12μ+sa1)].\begin{aligned} \mathcal N_{s,-r}^{(+)}(a_1,\mu) ={}& \frac{ \Gamma(1+2s a_1) }{ \Gamma\left( \frac12+s a_1+r\mu \right) } \\ &\times \exp \left[ \frac{\ii\pi}{2} (1+r) \left( \frac12-\mu+s a_1 \right) \right]. \end{aligned}

Define the target-row matrix

K1D,(+):=(K1D,sf,(+))T.K_{\infty1}^{D,(+)} := \left( \mathcal K_{\infty1}^{D,\mathrm{sf},(+)} \right)^{\mathsf T}.

Consequently,

C1D,(+)=(ND,(+))1K1D,(+)N1D,C_{\infty1}^{D,(+)} = \left( \mathscr N_\infty^{D,(+)} \right)^{-1} K_{\infty1}^{D,(+)} \mathscr N_1^D,

with entries

(C1D,(+))rs=χ+e(1r)L/2L1/2sa1+rμ×exp[sFD,1rFD,m2]×Ns,r(+)(a1,μ).\boxed{ \begin{aligned} \left( C_{\infty1}^{D,(+)} \right)_{rs} ={}& \chi_+ e^{(1-r)L/2} L^{-1/2-s a_1+r\mu} \\ &\times \exp \left[ \frac{ sF_{D,1}-rF_{D,m} }{2} \right] \\ &\times \mathcal N_{s,-r}^{(+)} \left( a_1,\mu \right). \end{aligned} }

Here no JJ_\infty is inserted:

det[Ns,r(+)]=2a1\det \left[ \mathcal N_{s,-r}^{(+)} \right] = 2a_1

in the stated row and column orders. Adding the small-LL collision-frame sign to this independently normalized large-LL kernel would violate the Wronskian determinant. In the reference sector the result is

detC1D,(+)=2a1LeLχ+2.\det C_{\infty1}^{D,(+)} = \frac{ 2a_1 }{L} e^L \chi_+^2.

The opposite lateral matrix is not obtained by replacing χ+\chi_+ with χ\chi_- in the displayed entries. The scalar phase χη\chi_\eta records only the gauge branch; the sectorial jump belongs to a different kernel N(η)\mathcal N^{(\eta)}, equivalently to Stokes transport. Once that kernel is continued consistently, the pathwise exact Abel law is

detC1D,(η)=2a1LeLχη2.\det C_{\infty1}^{D,(\eta)} = \frac{2a_1}{L} e^L\chi_\eta^2.

The next section gives the matrix transport law that constructs the adjacent sector.

A comparison ledger for the printed translations

Section titled “A comparison ledger for the printed translations”

Readers checking the formulas against the source can use the following diagnostic ledger. Each row is an apparent transcription or normalization inconsistency, not an official erratum.

Printed featureUnit-leading reconstructionAudit
Small-LL one-end conversion lacks a common endpoint factorRestore eL/2e^{L/2}detC10\det C_{10}
Infinity series is displayed with argument zz in a connection formulaUse its defining inverse-power argument z1z^{-1}Large-zz series
Small-LL collision frame has the opposite relative orientationApply JJ_\infty only in that collision identificationdetC1\det C_{\infty1}
A small-LL infinity component has an inconsistent FmF_m sign or lateral phaseDerive rFm/2-rF_m/2 and the phase from N1KT\mathscr N_\infty^{-1}K^{\mathsf T}Exact Whittaker slice
Large-LL finite-end conversions omit powersRestore Lsa0L^{-s a_0} and Lsa1L^{-s a_1}Local leading coefficients
Large-LL exponential infinity branch repeats the ++ sign labelUse the r=r=- branchFormal exponential
A component expression places a factor eL/2e^{\mp L/2} inside a power of LLKeep exponentials and powers separateDimensions and Abel determinant
A large-LL exponential coefficient uses the opposite sign of mm' in one gamma denominatorUse the sign fixed by Ns,r(+)(a1,μ)\mathcal N_{s,-r}^{(+)}(a_1,\mu)Whittaker-kernel definition

The safe procedure is structural: start with a block relation whose sign indices are explicit, derive both N\mathscr N matrices from the local leading terms, and multiply them in the frame equation. Copying four expanded components independently discards precisely the redundancy that makes the result auditable.

Suppose the sectorial infinity frames satisfy the book’s right-action law

Y(k+1)=Y(k)Sk.\boldsymbol Y_\infty^{(k+1)} = \boldsymbol Y_\infty^{(k)}S_k.

For one fixed regular source frame,

Hi=Y(k)Ci(k)=Y(k+1)Ci(k+1).\boldsymbol H_i = \boldsymbol Y_\infty^{(k)} C_{\infty i}^{(k)} = \boldsymbol Y_\infty^{(k+1)} C_{\infty i}^{(k+1)}.

Therefore

Ci(k+1)=Sk1Ci(k).C_{\infty i}^{(k+1)} = S_k^{-1} C_{\infty i}^{(k)}.

The block-to-one-sector kernel and the Stokes matrix have different jobs. The former constructs Ci(k)C_{\infty i}^{(k)}; the latter changes its target frame. Folding SkS_k into an unlabeled “irregular connection coefficient” destroys the information needed to continue it. The phase χη\chi_\eta tracks only the scalar-gauge branch, so changing χη\chi_\eta alone cannot implement this matrix jump. In particular, when δ=0\delta=0 the two phases coincide even though SkS_k can be nontrivial. On a fixed analytic branch detSk=1\det S_k=1, as required by the unchanged Abel determinant.

A factor ablation catches what a determinant cannot

Section titled “A factor ablation catches what a determinant cannot”

The genuinely four-singular general-Heun slice provides an entrywise test of the assembly theorem. A direct Frobenius match at z=t/2z=t/2 is compared with the level-two classical-block matrix: qHq_{\mathrm H} and both endpoint derivatives are truncated consistently through t2t^2. The same calculation is then repeated after setting both endpoint derivatives to zero while retaining the fusion kernel and all coordinate powers.

ttFully dressed level-two matrixEndpoint derivatives omitted
0.080.089.12×1059.12\times10^{-5}1.40×1021.40\times10^{-2}
0.040.041.09×1051.09\times10^{-5}6.83×1036.83\times10^{-3}
0.020.021.33×1061.33\times10^{-6}3.37×1033.37\times10^{-3}
0.010.011.64×1071.64\times10^{-7}1.68×1031.68\times10^{-3}

The fully dressed error decreases by approximately eight when tt is halved, as expected from the omitted O(t3)O(t^3) terms. The ablated error decreases only by approximately two because the missing derivative normalization starts at O(t)O(t).

At t=0.04t=0.04, all three computed determinants agree to the 18 digits reported here:

detCdirectdetCblockdetCablated1.03319670758600654.\det C_{\mathrm{direct}} \approx \det C_{\mathrm{block}} \approx \det C_{\mathrm{ablated}} \approx -1.03319670758600654.

The block and ablated determinants obey the Abel value algebraically; the direct-series determinant is numerical. The derivative factors have reciprocal determinant and cancel from the Abel law. This example makes the hierarchy of tests explicit:

  1. a determinant failure proves the matrix is wrong;
  2. a determinant success does not prove its entries are right;
  3. matching at two ordinary points tests the entire matrix.

The general-Heun reproducibility script generates the table and the ablation. For one row, run:

Terminal window
python3 public/code/advanced-ode/general-heun-connection-check.py \
--t 0.04 --order 2 --terms 120 --dps 50 \
--no-table --show-ablation --check-second-point

The reported run used Python 3.10.16, mpmath 1.3.0, 120 Frobenius terms, 50 decimal digits, and positive-real branches with 0<t<2/30<t<2/3; its measured wall time was 0.060.06 seconds on the audit machine. At t=0.04t=0.04, changing the series length from 80 through 180 terms leaves the displayed digits stable. Repeating the direct match at z=t/3z=t/3 changes the matrix by only 2.28×10222.28\times10^{-22} relative to the z=t/2z=t/2 result. The script rejects under-resolved precision or series-length inputs and warns when the direct-series determinant signals poor convergence near a disk boundary.

The exact Whittaker benchmark on the preceding page’s exact benchmark independently checks a sectorial kernel, its derivatives, Wronskian determinant, and adjacent-sector matrix law.

For a new conformal-block connection problem:

  1. Fix the scalar equation. State its gauge, parameters, and the exact accessory coordinate used by the block.
  2. Declare both ODE frames. Give their order, local coordinate, unit-leading behavior, logarithms, and continuation path.
  3. Choose one block chart. Do not mix FF with FDF_D, or a cVir=1c_{\mathrm{Vir}}=1 Fourier block with a classical block.
  4. Invert the accessory relation. Select a(q)a(q) or m(q)m'(q) and verify the relevant inverse Jacobian is nonzero.
  5. Differentiate before substitution. Hold the other block data fixed while computing the external and internal derivatives.
  6. Build the target-row kernel. Transpose a source-first coefficient array exactly once.
  7. Derive the endpoint ledgers. Read every factor in Ni,Nj\mathscr N_i,\mathscr N_j from a declared local or asymptotic leading term.
  8. Multiply the three factors. Use Cji=Nj1KjiNiC_{ji}=\mathscr N_j^{-1}K_{ji}\mathscr N_i before expanding components.
  9. Audit independently. Check the Abel determinant, match at two ordinary points, and verify any sector-change law.
  10. Add physical normalization last. Only then contract with boundary vectors or include reflection, pairing, flux, or vertex factors.

Record the claimed accuracy together with the result. “Leading classical,” “through LNL^N,” “through LNL^{-N},” and “including the first b2b^2 correction” are different deliverables.

Calling the fusion kernel the connection matrix. The kernel acts on block representatives. The two endpoint conversions are required before its entries refer to unit-leading ODE functions.

Using the accessory derivative as an endpoint derivative. tf\partial_t f selects the equation, while aif\partial_{a_i}f fixes a probe normalization. They are partial derivatives with different held data.

Differentiating after solving for the internal lift. The shifted block holds the internal channel fixed. Applying the chain rule to F(a(q),)F(a(q),\ldots) inserts a term not present in the fusion limit.

Evaluating the gauge at infinity. An irregular target is normalized by an asymptotic comparison in a sector. There is no finite endpoint value P()P(\infty).

Carrying a collision-frame sign into every chart. The small-LL kernel needs JJ_\infty for the declared frame; the large-LL kernel already has the correct column orientation. The Wronskian decides.

Assuming the order-one block term always contributes. gbgg_{\mathrm{bg}} cancels from a generic leading shifted ratio, but it enters the first b2b^2 correction and can survive in absolute chiral or physical normalizations.

Calling an ODE entry a physical amplitude. A unit-leading connection entry changes under endpoint rescaling. A physical observable must include the associated boundary or flux conventions.

A continuation formula is printed source first:

Bi,s=r=12AsrBj,r,A=(abcd).\mathcal B_{i,s} = \sum_{r=1}^2 \mathcal A_{sr}\mathcal B_{j,r}, \qquad \mathcal A = \begin{pmatrix} a&b\\ c&d \end{pmatrix}.

Let

Ni=diag(p,q),Nj=diag(u,v).\mathscr N_i=\operatorname{diag}(p,q), \qquad \mathscr N_j=\operatorname{diag}(u,v).

Construct the target-row connection matrix CjiC_{ji}. Which coefficient multiplies Hj,2H_{j,2} in the expansion of Hi,1H_{i,1}? Show why a determinant test alone cannot detect the error made by omitting the transpose.

Solution

The first index of Asr\mathcal A_{sr} labels a source solution, whereas the row of a frame matrix labels a target solution. Thus

Kji=AT=(acbd).K_{ji} = \mathcal A^{\mathsf T} = \begin{pmatrix} a&c\\ b&d \end{pmatrix}.

The two-sided conversion gives

Cji=Nj1KjiNi.C_{ji} = \mathscr N_j^{-1} K_{ji}\mathscr N_i.

Hence

Cji=(u100v1)(acbd)(p00q)=(ap/ucq/ubp/vdq/v).\begin{aligned} C_{ji} &= \begin{pmatrix} u^{-1}&0\\ 0&v^{-1} \end{pmatrix} \begin{pmatrix} a&c\\ b&d \end{pmatrix} \begin{pmatrix} p&0\\ 0&q \end{pmatrix} \\ &= \begin{pmatrix} ap/u&cq/u\\ bp/v&dq/v \end{pmatrix}. \end{aligned}

The first column expands Hi,1H_{i,1}, so the requested coefficient is bp/vbp/v. Omitting the transpose would instead give cp/vcp/v. Yet

det ⁣(Nj1ATNi)=det ⁣(Nj1ANi),\det\!\left( \mathscr N_j^{-1}\mathcal A^{\mathsf T}\mathscr N_i \right) = \det\!\left( \mathscr N_j^{-1}\mathcal A\mathscr N_i \right),

because detAT=detA\det\mathcal A^{\mathsf T}=\det\mathcal A. Direct entrywise matching, not the Abel determinant alone, catches this transpose error.

2. Recover the first correction to a shifted block

Section titled “2. Recover the first correction to a shifted block”

Let

λθ=λθb22,\lambda_\theta = \lambda-\frac{\theta b^2}{2},

and use the expansion of \LogFbg\Log\mathfrak F_{\mathrm{bg}} on this page. Find \LogRθ\Log R_\theta through order b2b^2.

Solution

Taylor expansion gives

f(λθ)f(λ)=θb22f+b48f+O(b6),g(λθ)g(λ)=θb22g+O(b4).\begin{aligned} f(\lambda_\theta)-f(\lambda) ={}& -\frac{\theta b^2}{2}f' + \frac{b^4}{8}f'' + O(b^6), \\ g(\lambda_\theta)-g(\lambda) ={}& -\frac{\theta b^2}{2}g' + O(b^4). \end{aligned}

After dividing the first line by b2b^2,

\LogRθ=θ2f+b2(18fθ2g)+O(b4).\Log R_\theta = -\frac{\theta}{2}f' + b^2 \left( \frac18f'' - \frac{\theta}{2}g' \right) + O(b^4).

The b2hb^2h term changes only at order b4b^4.

3. Recover both confluent-Heun determinant laws

Section titled “3. Recover both confluent-Heun determinant laws”

Use the determinants of the three small-LL endpoint matrices and the kernel identities

detMfr=a0a1,detK=2a1.\det M_{\mathrm{fr}} = -\frac{a_0}{a_1}, \qquad \det K=2a_1.

Find detC10\det C_{10} and detC1(η)\det C_{\infty1}^{(\eta)}.

Solution

For the two finite endpoints,

detN0=1,detN1=eL.\det\mathscr N_0=1, \qquad \det\mathscr N_1=e^L.

Therefore

detC10=1eL(a0a1)=a0a1eL.\det C_{10} = \frac1{e^L} \left( -\frac{a_0}{a_1} \right) = -\frac{a_0}{a_1}e^{-L}.

At infinity,

detN(η)=χη2L.\det\mathscr N_\infty^{(\eta)} = \chi_\eta^{-2}L.

Hence

detC1(η)=eLχη2L(2a1)=2a1LeLχη2.\det C_{\infty1}^{(\eta)} = \frac{ e^L }{ \chi_\eta^{-2}L } (2a_1) = \frac{2a_1}{L} e^L\chi_\eta^2.

Given

Y(k+1)=Y(k)Sk,\boldsymbol Y_\infty^{(k+1)} = \boldsymbol Y_\infty^{(k)}S_k,

derive the relation between Ci(k+1)C_{\infty i}^{(k+1)} and Ci(k)C_{\infty i}^{(k)}.

Solution

Write the same source frame in the two target frames:

Hi=Y(k)Ci(k)=Y(k)SkCi(k+1).\boldsymbol H_i = \boldsymbol Y_\infty^{(k)} C_{\infty i}^{(k)} = \boldsymbol Y_\infty^{(k)} S_kC_{\infty i}^{(k+1)}.

Linear independence gives

Ci(k)=SkCi(k+1),C_{\infty i}^{(k)} = S_kC_{\infty i}^{(k+1)},

and therefore

Ci(k+1)=Sk1Ci(k).C_{\infty i}^{(k+1)} = S_k^{-1}C_{\infty i}^{(k)}.

5. Show that a boundary scalar is frame independent

Section titled “5. Show that a boundary scalar is frame independent”

Let a scalar boundary function be

Q=jTCjivi.\mathcal Q = \ell_j^{\mathsf T} C_{ji}v_i.

If the endpoint frames change by Mi,MjM_i,M_j, find the transformations of viv_i and j\ell_j that leave Q\mathcal Q unchanged.

Solution

The connection matrix changes to

C~ji=Mj1CjiMi.\widetilde C_{ji} = M_j^{-1}C_{ji}M_i.

Choose

v~i=Mi1vi,~jT=jTMj.\widetilde v_i = M_i^{-1}v_i, \qquad \widetilde\ell_j^{\mathsf T} = \ell_j^{\mathsf T}M_j.

Then

Q~=jTMjMj1CjiMiMi1vi=jTCjivi=Q.\begin{aligned} \widetilde{\mathcal Q} &= \ell_j^{\mathsf T}M_j M_j^{-1}C_{ji}M_i M_i^{-1}v_i \\ &= \ell_j^{\mathsf T}C_{ji}v_i = \mathcal Q. \end{aligned}

The scalar can be invariant even though every matrix entry changes.

Why does deleting the endpoint derivative factors in the general-Heun test leave the determinant unchanged while spoiling the entries?

Solution

The derivative dressings are diagonal with reciprocal entries, so each has determinant one. Deleting them therefore leaves

detNidetNjdetKji\frac{ \det\mathscr N_i }{ \det\mathscr N_j } \det K_{ji}

unchanged. Their individual diagonal entries are not one, however, so they rescale different rows and columns and change all four connection coefficients. Abel’s determinant is necessary but not sufficient; direct matching supplies the missing entrywise test.

The next page separates Frobenius resonance from Kac or Zamolodchikov poles. That distinction is essential precisely when the generic endpoint matrices or shifted-block expansion used here cease to be uniform.