Skip to content

Global-Data Problems: Monodromy, Basis Changes, and Resonance

The chapter has spent seven pages separating objects that are often denoted by the same letter: monodromy matrices, Stokes factors, connection coefficients, boundary Wronskians, and several determinant constructions. This capstone asks whether those distinctions survive an actual calculation. The problems progress from convention checks to resonant limits, operator domains, and a reproducible hypergeometric computation.

Complete solutions are provided for the closed problems. The computational lab gives pseudocode and acceptance criteria rather than software-specific output, and the final research problem gives milestones rather than pretending that every model has a universal closed-form answer.

Use the chapter conventions throughout:

Φγ=ΦMγ,Φβ=ΦαCαβ,C~αβ=Hα1CαβHβ,Cαγ=CαβCβγ.\begin{aligned} \Phi^\gamma &= \Phi M_\gamma,\\ \Phi_\beta &= \Phi_\alpha C_{\alpha\beta},\\ \widetilde C_{\alpha\beta} &= H_\alpha^{-1}C_{\alpha\beta}H_\beta,\\ C_{\alpha\gamma} &= C_{\alpha\beta}C_{\beta\gamma}. \end{aligned}

The path product γ1γ2\gamma_1\gamma_2 traverses γ2\gamma_2 first and then γ1\gamma_1, so

Mγ1γ2=Mγ1Mγ2.M_{\gamma_1\gamma_2} = M_{\gamma_1}M_{\gamma_2}.

For a scalar equation in monic form,

y+p(z)y+q(z)y=0,y''+p(z)y'+q(z)y=0,

the book’s Wronskian and Abel identity are

Wr[f,g]=fgfg,Wr=pWr.\Wr[f,g]=fg'-f'g, \qquad \Wr'=-p\Wr.

Every answer involving a connection coefficient should name its ordered bases, direction, branches, and continuation path. Every spectral answer should additionally name the Hilbert space or resonance sheet and the operator domain or asymptotic boundary condition.

ProblemLevelPrincipal audit
1CoreRequested datum and missing hypotheses
2CoreBased-loop order and conjugacy
3IntermediateFour different transformations
4IntermediateConnection cocycle and boundary entry
5IntermediateA commuting Riemann–Hilbert reconstruction
6AdvancedWild factorization and residual gauge
7AdvancedCoalescing bases and Jordan monodromy
8AdvancedHurwitz moves and marked trace data
9AdvancedRobin boundary function and Fredholm ratio
10AdvancedResonant Legendre quantization
11Computational labExact and near-resonant numerical audits
12ResearchA complete global-data dossier

For each request below, identify the primary mathematical object and at least one datum or hypothesis that must be added before the request is well posed.

  1. Continue a normalized germ once around a puncture.
  2. Compare unit-leading Frobenius bases at two singular points.
  3. Measure the jump between two lateral asymptotic sums at an irregular singularity.
  4. Find the values of a parameter for which a left-selected solution also satisfies a right endpoint condition.
  5. Evaluate det(IK(λ))\det(I-K(\lambda)).
  6. Evaluate detζT\det_\zeta T.
Solution

The objects and their minimum ledgers are:

RequestPrimary objectData or hypothesis still required
Continued germA based monodromy matrix, or its representation classBase point, based homotopy class, ordered frame, loop orientation, and path-composition convention
Two local basesA path-labelled connection matrixExact basis normalizations, matrix direction, branches, and continuation path
Lateral asymptotic jumpA Stokes factor inside wild monodromy dataFormal normal form, sector labels, lateral summation directions, and orientation of the crossing
Two-end conditionA selected connection entry or boundary WronskianHilbert space and operator domain for eigenvalues, or time convention, asymptotic condition, and sheet for resonances
det(IK)\det(I-K)A Fredholm determinant, or possibly a regularized determinantA trace-class hypothesis for the ordinary determinant; if only KSpK\in\mathcal S_p, the appropriate detp\det_p must be declared
detζT\det_\zeta TA zeta-regularized operator determinantA closed operator with suitable discrete spectral behavior, a spectral cut when needed, a summability half-plane, and regular continuation to the evaluation point

The first four objects arise directly from solution spaces and analytic continuation. The last two are operator constructions with separate existence theorems. A boundary Wronskian is not promoted to a Fredholm or zeta determinant merely because its zeros are spectral.

2. A based-loop product in an explicit frame

Section titled “2. A based-loop product in an explicit frame”

Let

M0=(1a01),M1=(10b1),a,b0.M_0 = \begin{pmatrix} 1&a\\ 0&1 \end{pmatrix}, \qquad M_1 = \begin{pmatrix} 1&0\\ b&1 \end{pmatrix}, \qquad a,b\neq0.

The distinguished loops obey γ0γ1γ=1\gamma_0\gamma_1\gamma_\infty=1 in the book’s convention.

  1. Compute MM_\infty.
  2. Verify the full product and its determinant.
  3. Change the common base frame by H=diag(h,h1)H=\operatorname{diag}(h,h^{-1}). Compute the new M0M_0 and M1M_1, and identify a quantity built from a,ba,b that survives.
Solution

Because the representation respects the declared path product,

M=(M0M1)1.M_\infty = (M_0M_1)^{-1}.

Now

M0M1=(1+abab1),M_0M_1 = \begin{pmatrix} 1+ab&a\\ b&1 \end{pmatrix},

whose determinant is one. Therefore

M=(1ab1+ab).M_\infty = \begin{pmatrix} 1&-a\\ -b&1+ab \end{pmatrix}.

Direct multiplication gives M0M1M=IM_0M_1M_\infty=I. Each matrix has determinant one, so the determinant relation is also satisfied.

The frame change acts by simultaneous conjugation:

M~j=H1MjH.\widetilde M_j = H^{-1}M_jH.

Hence

M~0=(1h2a01),M~1=(10h2b1).\widetilde M_0 = \begin{pmatrix} 1&h^{-2}a\\ 0&1 \end{pmatrix}, \qquad \widetilde M_1 = \begin{pmatrix} 1&0\\ h^2b&1 \end{pmatrix}.

The separate off-diagonal entries depend on the frame, but their product abab does not. Equivalently,

trM=2+ab\operatorname{tr}M_\infty=2+ab

is a simultaneous-conjugacy invariant. The loop relation itself is also preserved:

M~0M~1M~=H1(M0M1M)H=I.\widetilde M_0\widetilde M_1\widetilde M_\infty = H^{-1} \bigl( M_0M_1M_\infty \bigr) H =I.

3. Four transformations that should not be conflated

Section titled “3. Four transformations that should not be conflated”

For a flat frame Φ\Phi, a connection matrix CαβC_{\alpha\beta}, and a boundary function E(λ)E(\lambda) built from selected solution lines, compare the following operations:

  1. a constant change of the common global frame;
  2. independent changes of normalized local bases;
  3. a single-valued invertible left gauge G(z)G(z);
  4. a multivalued scalar gauge g(z)g(z) with gγ=χ(γ)gg^\gamma=\chi(\gamma)g.

State how monodromy and connection data transform, and which spectral information in EE is invariant under a nowhere-zero analytic renormalization.

Solution

A constant global frame change Φ~=ΦH\widetilde\Phi=\Phi H gives

M~γ=H1MγH.\widetilde M_\gamma = H^{-1}M_\gamma H.

Thus traces, determinants, Jordan form, and all simultaneous-conjugacy invariants survive, while individual entries generally do not.

Independent local changes

Φ~α=ΦαHα,Φ~β=ΦβHβ\widetilde\Phi_\alpha=\Phi_\alpha H_\alpha, \qquad \widetilde\Phi_\beta=\Phi_\beta H_\beta

give the two-sided law

C~αβ=Hα1CαβHβ.\widetilde C_{\alpha\beta} = H_\alpha^{-1}C_{\alpha\beta}H_\beta.

This is not usually a conjugation because the two endpoint normalizations are independent. If HαH_\alpha and HβH_\beta preserve the selected boundary lines and rescale their generators by nowhere-zero analytic factors, then the associated scalar boundary function changes by a nowhere-zero analytic factor. Its zero set and zero multiplicities are unchanged.

For a single-valued invertible left gauge,

Φ~=GΦ,\widetilde\Phi=G\Phi,

analytic continuation gives

Φ~γ=GΦγ=Φ~Mγ.\widetilde\Phi^\gamma = G\Phi^\gamma = \widetilde\Phi M_\gamma.

The numerical right monodromy matrices are unchanged when the same solution frame is retained. Consistently gauged local frames also retain the same connection matrices.

A multivalued scalar gauge instead yields

M~γ=χ(γ)Mγ.\widetilde M_\gamma = \chi(\gamma)M_\gamma.

It twists the GL(n)GL(n) representation by a character. Projective monodromy is unchanged, but traces and determinants of a chosen lift generally are not. These four operations coincide only in special circumstances and should never be inferred from the visual similarity of their formulas.

4. A cocycle, a path change, and one boundary zero

Section titled “4. A cocycle, a path change, and one boundary zero”

Suppose three normalized frames satisfy

Φ0=Φ1A,Φ1=ΦB.\Phi_0=\Phi_1A, \qquad \Phi_1=\Phi_\infty B.
  1. Find C0C_{\infty0} in Φ0=ΦC0\Phi_0=\Phi_\infty C_{\infty0}.

  2. Let Φ~j=ΦjHj\widetilde\Phi_j=\Phi_jH_j. Verify the transformed cocycle.

  3. If a new comparison path replaces the endpoint frames by ΦD\Phi_\infty D_\infty and Φ0D0\Phi_0D_0, find the new C0C_{\infty0}.

  4. Write

    C0=(pqrs).C_{\infty0} = \begin{pmatrix} p&q\\ r&s \end{pmatrix}.

    If the first solution at zero is required to lie in the first selected solution line at infinity, identify the boundary equation. Explain why diagonal, nowhere-zero analytic renormalizations preserve it.

Solution

Composition gives

Φ0=ΦBA,\Phi_0 = \Phi_\infty BA,

so

C0=BA.C_{\infty0}=BA.

The two factors transform as

A~=H11AH0,B~=H1BH1.\widetilde A = H_1^{-1}AH_0, \qquad \widetilde B = H_\infty^{-1}BH_1.

Their product is

B~A~=H1BAH0=C~0,\widetilde B\widetilde A = H_\infty^{-1}BAH_0 = \widetilde C_{\infty0},

and the middle normalization H1H_1 cancels exactly as a cocycle requires.

For the changed endpoint branches, solve

Φ0D0=ΦDC0new.\Phi_0D_0 = \Phi_\infty D_\infty C_{\infty0}^{\mathrm{new}}.

This yields

C0new=D1C0D0.C_{\infty0}^{\mathrm{new}} = D_\infty^{-1}C_{\infty0}D_0.

Finally, the first column says

f0=ph1+rh2.f_0=p\,h_1+r\,h_2.

Requiring f0f_0 to lie in the line spanned by h1h_1 is therefore the boundary equation

r(λ)=0.r(\lambda)=0.

Diagonal rescalings multiply rr by nonzero analytic factors from the source and target normalizations, so its zeros and their orders survive. A general nondiagonal basis change can mix the selected line with its complement and describes a different boundary condition.

Reconstruction, wild data, and marked representations

Section titled “Reconstruction, wild data, and marked representations”

5. A commuting Riemann–Hilbert reconstruction

Section titled “5. A commuting Riemann–Hilbert reconstruction”

Let

A0=(α00α),A1=(β00β),A_0= \begin{pmatrix} \alpha&0\\ 0&-\alpha \end{pmatrix}, \qquad A_1= \begin{pmatrix} \beta&0\\ 0&-\beta \end{pmatrix},

and consider the prescribed local monodromies

M0=e2πiA0,M1=e2πiA1.M_0=\ee^{2\pi\ii A_0}, \qquad M_1=\ee^{2\pi\ii A_1}.
  1. Construct a Fuchsian system on P1{0,1,}\mathbb P^1\setminus\{0,1,\infty\} with these monodromies.
  2. Determine its residue and local monodromy at infinity.
  3. Show that the same monodromy representation does not determine unique residue logarithms.
  4. Explain geometrically what the ambiguity changes.
Solution

Because A0A_0 and A1A_1 commute, chosen branches of the matrix powers give

Y(z)=zA0(z1)A1.Y(z)=z^{A_0}(z-1)^{A_1}.

It is a fundamental matrix of

 ⁣dY ⁣dz=(A0z+A1z1)Y.\frac{\dd Y}{\dd z} = \left( \frac{A_0}{z} + \frac{A_1}{z-1} \right)Y.

A positive loop about zero adds 2πi2\pi\ii to \Logz\Log z and gives M0=e2πiA0M_0=\ee^{2\pi\ii A_0}; the analogous loop about one gives M1=e2πiA1M_1=\ee^{2\pi\ii A_1}.

In the local coordinate w=1/zw=1/z, the residue at infinity is

A=(A0+A1).A_\infty=-(A_0+A_1).

Therefore

M=e2πiA=(M0M1)1,\begin{aligned} M_\infty &= \ee^{2\pi\ii A_\infty}\\ &= \left( M_0M_1 \right)^{-1}, \end{aligned}

and the based product M0M1M=IM_0M_1M_\infty=I holds. Here commutativity removes the ordering difficulty; it does not remove the need to state the distinguished loop system.

Now choose traceless integral diagonal matrices

Kj=(mj00mj),mjZ.K_j= \begin{pmatrix} m_j&0\\ 0&-m_j \end{pmatrix}, \qquad m_j\in\mathbb Z.

Then

e2πi(Aj+Kj)=e2πiAj.\ee^{2\pi\ii(A_j+K_j)} = \ee^{2\pi\ii A_j}.

The single-valued meromorphic gauge

YzK0(z1)K1YY\longmapsto z^{K_0}(z-1)^{K_1}Y

shifts the finite residues by K0K_0 and K1K_1 while leaving the local system on the punctured sphere unchanged. Thus the representation determines the exponents only modulo integers. Choosing particular logarithms amounts to choosing how the flat bundle is extended across the punctures—a logarithmic lattice, such as a Deligne extension after a preferred strip for the real parts has been fixed.

This example is deliberately abelian. In a noncommuting reconstruction, existence of a logarithmic connection with a prescribed pole structure is an additional Riemann–Hilbert problem; one cannot obtain it by choosing three unrelated matrix logarithms.

6. Wild factors hidden by ordinary monodromy

Section titled “6. Wild factors hidden by ordinary monodromy”

At a rank-two irregular singularity, assume that the formal exponential type has two distinct, labelled one-dimensional blocks. Block swaps are not allowed: residual formal gauges must preserve this exponential grading. Take

Mf=(q00q1),S0=I+sE12,S1=I+tE21,M_{\mathrm f} = \begin{pmatrix} q&0\\ 0&q^{-1} \end{pmatrix}, \qquad S_0=I+sE_{12}, \qquad S_1=I+tE_{21},

with q0q\neq0. In the orientation and fixed-frame convention of the wild-monodromy page, the actual local monodromy is

Mloc=MfS11S01.M_{\mathrm{loc}} = M_{\mathrm f}S_1^{-1}S_0^{-1}.
  1. Compute MlocM_{\mathrm{loc}}, its determinant, and its trace.
  2. Determine how ss and tt transform under the residual formal-frame change H=diag(h,h1)H=\operatorname{diag}(h,h^{-1}).
  3. Which combination of the two Stokes multipliers can be recovered from trMloc\operatorname{tr}M_{\mathrm{loc}} when qq is known?
  4. Explain precisely what ordinary monodromy forgets.
Solution

The inverse Stokes factors give

S11S01=(1st1+st),S_1^{-1}S_0^{-1} = \begin{pmatrix} 1&-s\\ -t&1+st \end{pmatrix},

so

Mloc=(qqsq1tq1(1+st)).M_{\mathrm{loc}} = \begin{pmatrix} q&-qs\\ -q^{-1}t&q^{-1}(1+st) \end{pmatrix}.

Consequently,

detMloc=1,trMloc=q+q1(1+st).\det M_{\mathrm{loc}}=1, \qquad \operatorname{tr}M_{\mathrm{loc}} = q+q^{-1}(1+st).

Because HH preserves the exponential grading—and commutes with MfM_{\mathrm f}—it is a residual normalization of the formal eigenlines. Conjugation gives

H1S0H=I+h2sE12,H1S1H=I+h2tE21.\begin{aligned} H^{-1}S_0H &= I+h^{-2}sE_{12},\\ H^{-1}S_1H &= I+h^2tE_{21}. \end{aligned}

Thus ss and tt are framed quantities, whereas stst is invariant under this centralizer action. If qq and the trace are known, then

st=q(trMlocqq1).st = q\left( \operatorname{tr}M_{\mathrm{loc}} -q-q^{-1} \right).

Even this invariant does not reconstruct the direction-labelled factorization. For example, if st=0st=0 and q21q^2\neq1, the cases

(s,t)=(0,0),(s,0),(0,t)(s,t)=(0,0), \qquad (s,0), \qquad (0,t)

all have ordinary monodromy conjugate to diag(q,q1)\operatorname{diag}(q,q^{-1}), provided the displayed nonzero multiplier is allowed. The three wild data sets differ: one has no jump, and the other two place a jump in different Stokes groups. Ordinary conjugacy retains only the product matrix; wild monodromy retains the formal monodromy, ordered singular directions, sector normalizations, and Stokes factors.

7. Coalescing Euler modes and Jordan monodromy

Section titled “7. Coalescing Euler modes and Jordan monodromy”

For ε0\varepsilon\neq0, consider

x2y+(1ε)xy=0x^2y''+(1-\varepsilon)xy'=0

on the slit plane with a chosen \Logx\Log x. A nonresonant basis is

fε(x)=1,gε(x)=xε.f_\varepsilon(x)=1, \qquad g_\varepsilon(x)=x^\varepsilon.
  1. Explain why this basis becomes defective as ε0\varepsilon\to0.
  2. Construct a basis with a finite, nonzero Wronskian limit.
  3. Transform the monodromy into that basis and take the limit.
  4. Check the answer directly from the limiting differential equation.
Solution

The Wronskian of the raw pair is

Wr[fε,gε]=εxε1,\Wr[f_\varepsilon,g_\varepsilon] = \varepsilon x^{\varepsilon-1},

which tends to zero because both columns tend to 11. The singularity is in the chosen basis, not in the two-dimensional solution space.

Set

uε=fε,vε=gεfεε.u_\varepsilon=f_\varepsilon, \qquad v_\varepsilon = \frac{ g_\varepsilon-f_\varepsilon }{\varepsilon}.

Then

Wr[uε,vε]=xε1x1,\Wr[u_\varepsilon,v_\varepsilon] = x^{\varepsilon-1} \longrightarrow x^{-1},

and vε\Logxv_\varepsilon\to\Log x. If

Φε=(fε,gε),Ψε=(uε,vε),\Phi_\varepsilon = \left( f_\varepsilon,g_\varepsilon \right), \qquad \Psi_\varepsilon = \left( u_\varepsilon,v_\varepsilon \right),

then

Φε=ΨεSε,Sε=(110ε).\Phi_\varepsilon = \Psi_\varepsilon S_\varepsilon, \qquad S_\varepsilon = \begin{pmatrix} 1&1\\ 0&\varepsilon \end{pmatrix}.

In the raw basis, positive continuation around zero is

Dε=diag(1,e2πiε).D_\varepsilon = \operatorname{diag} \left( 1,\ee^{2\pi\ii\varepsilon} \right).

The right monodromy in the convergent basis is therefore

M^ε=SεDεSε1=(1e2πiε1ε0e2πiε).\begin{aligned} \widehat M_\varepsilon &= S_\varepsilon D_\varepsilon S_\varepsilon^{-1}\\ &= \begin{pmatrix} 1& \dfrac{ \ee^{2\pi\ii\varepsilon}-1 }{\varepsilon}\\ 0&\ee^{2\pi\ii\varepsilon} \end{pmatrix}. \end{aligned}

Taking the limit gives the nontrivial Jordan matrix

M^0=(12πi01).\widehat M_0 = \begin{pmatrix} 1&2\pi\ii\\ 0&1 \end{pmatrix}.

At ε=0\varepsilon=0, the equation is x2y+xy=0x^2y''+xy'=0 and the limiting basis is (1,\Logx)(1,\Log x). Since a positive loop sends \Logx\Log x to \Logx+2πi\Log x+2\pi\ii, the direct continuation gives the same matrix. Taking the limit of the diagonal eigenvalues alone would have produced the identity and lost the logarithmic extension data.

8. A Hurwitz move on marked character data

Section titled “8. A Hurwitz move on marked character data”

Let a four-puncture monodromy tuple satisfy ABCD=IABCD=I, with

A=(1a01),B=(10b1),C=(r00r1),A= \begin{pmatrix} 1&a\\ 0&1 \end{pmatrix}, \quad B= \begin{pmatrix} 1&0\\ b&1 \end{pmatrix}, \quad C= \begin{pmatrix} r&0\\ 0&r^{-1} \end{pmatrix},

where abr(rr1)0abr(r-r^{-1})\neq0, and set D=(ABC)1D=(ABC)^{-1}. Perform the Hurwitz move

(A,B,C,D)(B,B1AB,C,D).(A,B,C,D) \longmapsto \left( B,B^{-1}AB,C,D \right).
  1. Verify that the product constraint survives and that the first two local conjugacy classes are exchanged.
  2. Compute the three pair traces tr(AB)\operatorname{tr}(AB), tr(BC)\operatorname{tr}(BC), and tr(AC)\operatorname{tr}(AC) before the move.
  3. Compute the corresponding traces after the move.
  4. Why is this a change of marking rather than an ordinary simultaneous conjugation?
Solution

Write

A=B,B=B1AB.A'=B, \qquad B'=B^{-1}AB.

Then

AB=AB,A'B'=AB,

so ABCD=ABCD=IA'B'CD=ABCD=I. Moreover, AA' is the old BB, while BB' is conjugate to the old AA; hence their local conjugacy classes are exchanged.

Before the move,

tr(AB)=2+ab,tr(BC)=r+r1,tr(AC)=r+r1.\begin{aligned} \operatorname{tr}(AB) &= 2+ab,\\ \operatorname{tr}(BC) &= r+r^{-1},\\ \operatorname{tr}(AC) &= r+r^{-1}. \end{aligned}

Direct calculation gives

B1AB=(1+abaab21ab).B^{-1}AB = \begin{pmatrix} 1+ab&a\\ -ab^2&1-ab \end{pmatrix}.

Therefore the new pair traces are

tr(AB)=2+ab,tr(AC)=r+r1,tr(BC)=r+r1+ab(rr1).\begin{aligned} \operatorname{tr}(A'B') &= 2+ab,\\ \operatorname{tr}(A'C) &= r+r^{-1},\\ \operatorname{tr}(B'C) &= r+r^{-1} + ab(r-r^{-1}). \end{aligned}

The last value differs from both old traces under the stated assumptions. A simultaneous conjugation preserves every trace of every labelled word, so it cannot produce this change. The Hurwitz move instead changes the distinguished generators of the punctured-sphere fundamental group. It is a braid or mapping-class action on the marked character variety: local conjugacy classes and the total product remain compatible, but the trace coordinates attached to labelled loop words transform.

9. A Robin boundary function and a Fredholm ratio

Section titled “9. A Robin boundary function and a Fredholm ratio”

Fix L>0L>0 and hRh\in\mathbb R. Let

Th= ⁣d2 ⁣dx2T_h=-\frac{\dd^2}{\dd x^2}

in L2(0,L)L^2(0,L), with

D(Th)={yH2(0,L):y(0)=0, y(L)+hy(L)=0}.\mathcal D(T_h) = \left\{ y\in H^2(0,L): y(0)=0,\ y'(L)+hy(L)=0 \right\}.

Let k2=λk^2=\lambda and use the entire-in-λ\lambda left-normalized solution

u(x,λ)=sin(kx)k,u(0,λ)=0,u(0,λ)=1.u(x,\lambda) = \frac{\sin(kx)}{k}, \qquad u(0,\lambda)=0, \qquad u'(0,\lambda)=1.
  1. Construct an entire boundary function and identify its zeros.
  2. Prove directly that its eigenvalue zeros are simple.
  3. When is its normalized ratio an ordinary Fredholm determinant?
  4. What fails at the exceptional value h=1/Lh=-1/L?
Solution

Applying the right boundary functional gives

Fh(λ)=u(L,λ)+hu(L,λ)=cos(kL)+hsin(kL)k.\begin{aligned} F_h(\lambda) &= u'(L,\lambda)+hu(L,\lambda)\\ &= \cos(kL) + h\frac{\sin(kL)}{k}. \end{aligned}

Both apparent functions of λ\sqrt\lambda have power series in λ\lambda, so FhF_h is entire. Its zeros are exactly the eigenvalues of the specified self-adjoint domain: the left-normalized solution then also satisfies the Robin condition.

Differentiate

u=λu-u''=\lambda u

with respect to λ\lambda, writing u˙=λu\dot u=\partial_\lambda u. The Wronskian identity is

 ⁣d ⁣dxWr[u,u˙]=u2.\frac{\dd}{\dd x} \Wr[u,\dot u] = -u^2.

At an eigenvalue λn\lambda_n, the left normalization gives zero Wronskian at x=0x=0, while the Robin condition gives

Wr[u,u˙](L)=u(L,λn)Fh(λn).\Wr[u,\dot u](L) = u(L,\lambda_n)F_h'(\lambda_n).

Hence

u(L,λn)Fh(λn)=0Lu(x,λn)2 ⁣dx.u(L,\lambda_n)F_h'(\lambda_n) = -\int_0^L u(x,\lambda_n)^2\,\dd x.

The endpoint value u(L,λn)u(L,\lambda_n) cannot vanish, since the Robin condition would then force both Cauchy data at LL to vanish. The integral is positive, so Fh(λn)0F_h'(\lambda_n)\neq0.

At zero,

Fh(0)=1+hL.F_h(0)=1+hL.

If h1/Lh\neq-1/L, then ThT_h is invertible. Its inverse is trace class, because the eigenvalues of ThT_h grow quadratically. Since FhF_h has order 1/21/2 as an entire function of λ\lambda, its normalized genus-zero factorization is

Fh(λ)Fh(0)=n(1λλn).\frac{F_h(\lambda)}{F_h(0)} = \prod_n \left( 1-\frac{\lambda}{\lambda_n} \right).

The same convergent product defines the Fredholm determinant, so

Fh(λ)Fh(0)=detF(IλTh1).\frac{F_h(\lambda)}{F_h(0)} = \det\nolimits_{\mathrm F} \left( I-\lambda T_h^{-1} \right).

This equality uses the declared domain, trace-class inverse, growth, and normalization; it does not follow from the zero set alone.

For h=1/Lh=-1/L, the function u(x,0)=xu(x,0)=x satisfies both boundary conditions. Thus zero is an eigenvalue, Fh(0)=0F_h(0)=0, and neither the displayed quotient nor Th1T_h^{-1} exists. The unnormalized boundary function remains valid. One may factor out the zero mode or normalize at a shifted spectral parameter, but that is a different determinant statement. An absolute zeta determinant would require its own continuation, normalization, and, for a non-self-adjoint problem, spectral-cut data.

Let

τy= ⁣d ⁣dx[(1x2)y]\tau y = -\frac{\dd}{\dd x} \left[ (1-x^2)y' \right]

and begin with the maximal domain

Dmax={y:  y, (1x2)yACloc(1,1),τyL2(1,1)}.\begin{aligned} \mathcal D_{\max} = \bigl\{ y:\;& y,\ (1-x^2)y'\in AC_{\mathrm{loc}}(-1,1),\\ & \tau y\in L^2(-1,1) \bigr\}. \end{aligned}

Consider the endpoint-regular self-adjoint restriction of τ\tau defined by

 ⁣d ⁣dx[(1x2)y]=λy-\frac{\dd}{\dd x} \left[ (1-x^2)y' \right] = \lambda y

in L2(1,1)L^2(-1,1), with the separated conditions

limx1+(1x2)y(x)=0,limx1(1x2)y(x)=0.\lim_{x\to-1^+} (1-x^2)y'(x)=0, \qquad \lim_{x\to1^-} (1-x^2)y'(x)=0.

Write λ=ν(ν+1)\lambda=\nu(\nu+1) and choose the solution regular at x=1x=1,

Pν(x)=2F1(ν, ν+11;1x2).P_\nu(x) = {}_2F_1 \left( \begin{matrix} -\nu,\ \nu+1\\ 1 \end{matrix} ; \frac{1-x}{2} \right).
  1. Show that the hypergeometric equation is zero-balanced at x=1x=-1.
  2. Compute the coefficient of the logarithmic branch there.
  3. Derive the spectrum selected by regularity at both endpoints.
  4. Explain why “the equation is resonant” is not itself a quantization condition.
Solution

Set

z=1x2,t=1z=1+x2.z=\frac{1-x}{2}, \qquad t=1-z=\frac{1+x}{2}.

The hypergeometric parameters are

a=ν,b=ν+1,c=1=a+b.a=-\nu, \qquad b=\nu+1, \qquad c=1=a+b.

Thus the exponent difference at z=1z=1, corresponding to x=1x=-1, is zero for every ν\nu. For νZ\nu\notin\mathbb Z, the zero-balanced continuation formula has leading term

2F1(a,b;a+b;z)=\Log(1z)+O(1)Γ(a)Γ(b){}_2F_1(a,b;a+b;z) = \frac{ -\Log(1-z)+O(1) }{ \Gamma(a)\Gamma(b) }

because Γ(a+b)=Γ(1)=1\Gamma(a+b)=\Gamma(1)=1. Euler’s reflection formula gives

Γ(ν)Γ(ν+1)=πsin(πν).\Gamma(-\nu)\Gamma(\nu+1) = -\frac{\pi}{\sin(\pi\nu)}.

Therefore, as x1+x\to-1^+,

Pν(x)=sin(πν)π\Log(1+x2)+O(1).P_\nu(x) = \frac{\sin(\pi\nu)}{\pi} \Log\left( \frac{1+x}{2} \right) + O(1).

The logarithmic coefficient sin(πν)/π\sin(\pi\nu)/\pi extends to integer ν\nu by the parameter limit, even though the two gamma factors in the intermediate formula then have poles. The endpoint-regular domain removes the logarithmic branch, so

sin(πν)=0.\sin(\pi\nu)=0.

Using the symmetry νν1\nu\leftrightarrow-\nu-1, choose ν=nZ0\nu=n\in\mathbb Z_{\ge0}. The eigenvalues and eigenfunctions are

λn=n(n+1),yn(x)=Pn(x).\lambda_n=n(n+1), \qquad y_n(x)=P_n(x).

The singular point was resonant for every value of ν\nu; most of those values produce a logarithm. Quantization occurs only when the parameter-dependent logarithmic connection coefficient vanishes. At integer ν\nu, the hypergeometric series terminates and the limiting solution is a Legendre polynomial. This is the distinction between a local exponent collision and a global two-end boundary condition.

11. Reproduce a hypergeometric global-data audit

Section titled “11. Reproduce a hypergeometric global-data audit”

Build one high-precision computation with

a=13,b=25,c=76,δ=1330.a=\frac13, \qquad b=\frac25, \qquad c=\frac76, \qquad \delta=\frac{13}{30}.

Use the unit-leading bases and upper/lower phase conventions from the hypergeometric benchmark. Your computation must do all of the following:

  1. form C10C_{10} and the four path-labelled matrices C0(s)C_{\infty0}^{(s)}, C1(s)C_{\infty1}^{(s)} for s=±1s=\pm1 from gamma quotients;
  2. verify the connection cocycle and the three exact determinant targets;
  3. assemble the lower- and upper-stem based monodromy products;
  4. repeat with c=a+b+εc=a+b+\varepsilon and compare the raw z=1z=1 basis with its logarithmic recombination as ε0\varepsilon\to0.
Solution and acceptance criteria

The exact determinant targets are

detC10=513,detC0(+)=52e17πi/30,detC1(+)=132e13πi/30.\begin{aligned} \det C_{10} &= \frac5{13},\\ \det C_{\infty0}^{(+)} &= \frac52\ee^{17\pi\ii/30},\\ \det C_{\infty1}^{(+)} &= -\frac{13}{2}\ee^{-13\pi\ii/30}. \end{aligned}

The upper matrices must obey

C0(+)=C1(+)C10.C_{\infty0}^{(+)} = C_{\infty1}^{(+)}C_{10}.

With the lower stem, transport the three diagonal local matrices into the zero frame and check

M0(0)M1(0)M(0)=I.M_0^{(0)}M_1^{(0)}M_\infty^{(0)}=I.

With the upper stem and its corresponding infinity matrix, check instead

M1(0)M0(0)M(+)=I.M_1^{(0)}M_0^{(0)}M_\infty^{(+)}=I.

Testing the wrong word against the right matrices is expected to fail.

A language-neutral implementation can follow this outline:

set working precision P and guard digits G
define gamma quotients A_f, A_g, B_f, B_g
construct C10 and Cinf0(s), Cinf1(s) for s = -1, +1
record every logarithm branch and every lateral phase
for each claimed identity:
compute a scale-aware matrix residual
repeat at P, 2P, and 4P digits
compare determinants separately with their exact targets
construct D0, D1, Dinf
transport them with the declared connection matrices
test the lower word and the upper word independently
for epsilon in 10^(-2), 10^(-4), 10^(-6), ...:
set c = a + b + epsilon
form the raw C10(epsilon)
set S = [[1, 1], [0, epsilon]]
form Chat(epsilon) = S * C10(epsilon)
record condition numbers of C10 and Chat
compare Chat with its gamma–digamma limit

The numerical report passes only if:

  • relative determinant errors and scaled matrix residuals decrease as the working precision increases;
  • the upper cocycle and both correctly ordered monodromy words converge to zero residual;
  • the raw near-resonant condition number grows on the scale 1/ε1/|\varepsilon|, while the transformed matrix remains bounded and converges when a+b1a+b\neq1;
  • an independent check—Wronskian evaluation, direct ODE integration along the declared paths, or overlapping local series—agrees with the matrix calculation;
  • the report records software and version, precision, branch conventions, path geometry, matrix norm, and all parameter values.

A table of residual versus precision is evidence; a single string of machine digits is not. Near resonance, evaluate the recombined expressions directly or with sufficient guard digits, since subtracting two O(1/ε)O(1/\varepsilon) terms at ordinary precision recreates the instability that the basis change was meant to remove.

Choose one nontrivial model: a Bessel, Airy, Whittaker, Heun, or compactly-supported Schrödinger equation is suitable. Formulate one global question—connection, scattering, bound-state, resonance, or quasinormal-mode—and answer it as completely as the model permits.

Your dossier must separate the following layers:

  1. differential equation, parameter domain, and singularity types;
  2. normalized local or asymptotic bases, including branches and sectors;
  3. base point, path system, and ordinary or wild monodromy data;
  4. the connection entry or boundary Wronskian that answers the question;
  5. Hilbert space and operator domain, or time convention, outgoing condition, and resonance sheet;
  6. the status of any Fredholm, zeta, or canonical-product determinant;
  7. at least two independent analytic or numerical checks;
  8. a ledger marking each statement as proved, numerically supported, or conjectural.
Guidance and known milestones

A successful dossier is reproducible before it is ambitious. Begin with a one-page convention ledger. Normalize every basis by a leading coefficient or Cauchy datum; draw or describe every cut and path; and reserve different symbols for local monodromy, Stokes factors, connection matrices, and spectral functions.

Then meet these milestones:

  1. Derive local exponents or formal exponential parts directly from the equation.
  2. Compute at least one connection quantity in two independent ways, such as a Wronskian ratio and numerical continuation.
  3. Transport every local monodromy matrix to one base frame before testing the global product.
  4. If an irregular point is present, give the ordered Stokes factorization and its residual centralizer action.
  5. Derive the spectral condition from the stated domain or asymptotics, rather than importing a familiar spectrum from another realization.
  6. State separately whether the resulting scalar function is merely a boundary function, a canonical product, a Fredholm determinant, or a zeta determinant. Cite the comparison theorem for every asserted equality.
  7. Include machine-readable parameters and a precision-convergence table. One check should be structurally different from the formula being tested.

For a Bessel model, the integer-order logarithmic limit is a natural resonance test. Airy supplies a clean Stokes audit. Whittaker combines regular and irregular singularities. Heun exposes accessory-parameter dependence, while a compactly-supported Schrödinger model makes the distinction between physical and nonphysical sheets unavoidable. A negative result—such as proving that the available data do not determine a claimed determinant—is a valid research conclusion when the missing hypothesis is identified precisely.