Skip to content

Just-in-Time Mathematical Toolkit

This toolkit develops only the topology and spectral theory used repeatedly in the book. Its purpose is operational: after reading it, a loop should mean a monodromy matrix, a lifted contour should mean a period or Voros symbol, and a boundary condition should mean part of an operator domain. The later chapters return to each construction where more precision is needed.

Fix a connected Riemann surface Σ\Sigma, a finite singular set DD, and the punctured surface

X=ΣD.X=\Sigma\setminus D.

The same space supports two related descriptions of a linear ODE:

  • topology records how solutions continue along paths in XX;
  • analysis records which continued solutions satisfy the boundary or growth conditions of a spectral problem.

Keeping these roles separate prevents two common category errors: treating monodromy as a basis-independent matrix, and treating every zero of a connection coefficient as an eigenvalue of a self-adjoint operator.

A path from aa to bb is a continuous map γ:[0,1]X\gamma:[0,1]\to X with γ(0)=a\gamma(0)=a and γ(1)=b\gamma(1)=b. Two such paths are equivalent when one can be continuously deformed into the other while both endpoints remain fixed. Analytic continuation on XX depends on this fixed-endpoint homotopy class, not on the particular parametrization.

Choose a base point zXz_*\in X. Homotopy classes of loops beginning and ending at zz_* form the fundamental group π1(X,z)\pi_1(X,z_*). The identity is the constant loop, and the inverse loop traverses the same curve backwards. We use the path-product convention fixed on the notation page: γ1γ2\gamma_1\gamma_2 traverses γ2\gamma_2 first and then γ1\gamma_1.

For the sphere punctured at nn points a1,,ana_1,\ldots,a_n, one may choose positively oriented based loops with a single relation

γ1γ2γn=1.\gamma_1\gamma_2\cdots\gamma_n=1.

The ordering is part of the choice of cuts and generators. For X=C{0,1}X=\mathbb C\setminus\{0,1\}, the loops about 00 and 11 freely generate π1(X,z)\pi_1(X,z_*); the loop about infinity is then fixed by the sphere relation.

Paths between several endpoints form the fundamental groupoid rather than a single group. This is the natural language for connection problems: monodromy uses loops at one base point, whereas a connection matrix compares frames normalized in two different regions.

A rank-rr complex local system L\mathcal L on XX is, informally, a collection of rr-dimensional solution spaces that can be canonically identified along sufficiently short paths. A nonsingular linear system

 ⁣dY ⁣dz=A(z)Y\frac{\dd Y}{\dd z}=A(z)Y

defines such a local system: its local horizontal sections are its solutions. Continuation around a based loop gives a representation

ρ:π1(X,z)GL(r,C),ρ(γ)=Mγ.\rho:\pi_1(X,z_*)\longrightarrow GL(r,\mathbb C), \qquad \rho(\gamma)=M_\gamma.

With the book’s right-action convention, Φγ=ΦMγ\Phi^\gamma=\Phi M_\gamma, the chosen path product gives

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

A fundamental matrix Φ(z)\Phi(z_*) is a framing of the fiber at the base point. Replacing it by Φ(z)G\Phi(z_*)G changes every monodromy matrix by

MγG1MγG.M_\gamma\longmapsto G^{-1}M_\gamma G.

Consequently, an unframed local system determines a representation only up to simultaneous conjugation. Traces, characteristic polynomials, and conjugacy classes are invariant; individual matrix entries are not.

At an irregular singularity, the local system on the punctured neighborhood remembers ordinary monodromy but not the full asymptotic classification. Formal exponential factors, formal monodromy, sectorial normalizations, and Stokes matrices are additional data developed in Chapters 1 and 2.

If Φα\Phi_\alpha and Φβ\Phi_\beta are normalized frames in two regions, then

Φβ=ΦαCαβ.\Phi_\beta=\Phi_\alpha C_{\alpha\beta}.

The matrix CαβC_{\alpha\beta} belongs to groupoid data: it depends on the continuation path and both endpoint frames. Monodromy data and connection matrices are related, but neither is a basis-free synonym for the other.

For a Fuchsian equation on P1{0,1,}\mathbb P^1\setminus\{0,1,\infty\}, choose generators satisfying

γ0γ1γ=1.\gamma_0\gamma_1\gamma_\infty=1.

Then

M0M1M=I.M_0M_1M_\infty=I.

This identity checks matrix order, but it does not determine the connection matrix between a basis at 00 and a basis at 11. That comparison also needs the local normalizations and a continuation path.

Riemann surfaces and multivalued functions

Section titled “Riemann surfaces and multivalued functions”

A Riemann surface is a one-complex-dimensional manifold with holomorphic coordinate changes. It turns a multivalued expression on the zz-plane into a single-valued object on a suitable covering surface.

For example, y2=zy^2=z defines a two-sheeted cover of the punctured plane. Continuing once around z=0z=0 exchanges the sheets; continuing twice returns to the original value. A branch cut is a planar bookkeeping device that selects one sheet over a cut domain. It is not an intrinsic boundary of the Riemann surface.

In WKB analysis, let

φ=Q(z)( ⁣dz)2\varphi=Q(z)(\dd z)^2

be the quadratic differential in a chosen normal-form coordinate. The coordinate-invariant spectral cover over U=ΣPoles(φ)U=\Sigma\setminus\operatorname{Poles}(\varphi) is the subspace of TUT^*U defined by

Σ^:λ2=φ.\widehat\Sigma: \qquad \lambda^2=\varphi.

Writing λ=y ⁣dz\lambda=y\,\dd z gives the familiar local equation

y2=Q(z),y^2=Q(z),

One then normalizes this cover and, when useful, compactifies it over the poles. A zero or pole of φ\varphi of odd order is a branch point of the normalized double cover. Even-order zeros and poles are not branch points. Zeros of any positive order remain turning points or higher WKB critical points, while an even-order pole generally lifts to poles of the WKB one-form λ\lambda.

If a degree-two cover π:Σ^Σ\pi:\widehat\Sigma\to\Sigma has rr simple branch points, the Riemann–Hurwitz formula gives

2g(Σ^)2=2(2g(Σ)2)+r.2g(\widehat\Sigma)-2 =2\bigl(2g(\Sigma)-2\bigr)+r.

For a cover of the sphere, this reduces to g(Σ^)=(r2)/2g(\widehat\Sigma)=(r-2)/2. Thus a generic quartic curve

y2=j=14(zaj)y^2=\prod_{j=1}^{4}(z-a_j)

has genus one, while y2=z21y^2=z^2-1 has genus zero after compactification.

Let SS be the page-specific set of lifted poles and other points removed from the normalized cover, and set

X^=Σ^S.\widehat X=\widehat\Sigma\setminus S.

An oriented closed curve on X^\widehat X determines a class in the absolute homology group

H1(X^,Z).H_1(\widehat X,\mathbb Z).

Two closed curves represent the same class when their difference is the boundary of an oriented two-chain. If a meromorphic one-form λ\lambda is holomorphic throughout the swept region, homologous cycles have equal periods:

γ1λ=γ2λ.\oint_{\gamma_1}\lambda =\oint_{\gamma_2}\lambda.

If the deformation crosses a pole, the periods differ by the appropriate 2πi2\pi\ii times residue. Homology on the compact surface alone does not erase this analytic information.

Open paths whose endpoints lie in a marked set PP define relative classes

H1(X^,P;Z).H_1(\widehat X,P;\mathbb Z).

The boundary map records their endpoints:

:H1(X^,P;Z)H0(P;Z).\partial: H_1(\widehat X,P;\mathbb Z) \longrightarrow H_0(P;\mathbb Z).

For connected X^\widehat X, the relevant part of the long exact sequence is

H1(X^;Z)H1(X^,P;Z)  H0(P;Z)H0(X^;Z).\begin{aligned} H_1(\widehat X;\mathbb Z) &\longrightarrow H_1(\widehat X,P;\mathbb Z)\\ &\xrightarrow{\ \partial\ }H_0(P;\mathbb Z) \longrightarrow H_0(\widehat X;\mathbb Z). \end{aligned}

Thus a relative path has a signed endpoint sum of zero. Closed cycles lie in the kernel of \partial and come from absolute homology, up to the preceding map. In exact WKB, a closed cycle typically supports a quantum period, while a path between turning points, poles, or asymptotic ends is relative and may require endpoint subtraction.

For a path-connected space,

H1(X^;Z)π1(X^,z)ab.H_1(\widehat X;\mathbb Z) \cong \pi_1(\widehat X,z_*)^{\mathrm{ab}}.

Homology therefore remembers winding and intersection data but forgets the order of noncommuting loops. It cannot replace the full fundamental group in a general monodromy representation.

There is no universal deleted set for every WKB problem. One common exact-WKB choice uses absolute cycles in H1(Σ^P^;Z)H_1(\widehat\Sigma\setminus\widehat P;\mathbb Z) and paths in a group of the form H1(Σ^P^0,P^;Z)H_1(\widehat\Sigma\setminus\widehat P_0,\widehat P_\infty;\mathbb Z), where the lifted turning points, simple poles, higher poles, and asymptotic ends are assigned roles dictated by the differential. Singular endpoints may instead require a real blow-up or a regularized path space. Each later page declares the group it actually uses.

A contour drawn only in the zz-plane is incomplete. One must also state:

  1. its lift to the chosen sheet;
  2. how it crosses the branch cuts;
  3. its orientation;
  4. whether its endpoints or punctures are fixed;
  5. which poles lie inside a deformation.

For a two-sheeted cover with involution ι:yy\iota:y\mapsto-y, the WKB differential changes sign:

ι(y ⁣dz)=y ⁣dz.\iota^*(y\,\dd z)=-y\,\dd z.

The same projected loop on opposite sheets therefore has the opposite classical action.

Two transverse oriented cycles on an oriented surface have an algebraic intersection number. At each crossing, assign +1+1 when the ordered tangent vectors of the first and second cycles agree with the surface orientation, and 1-1 otherwise; then add the local signs. The result is bilinear, homology-invariant, and antisymmetric:

γ1,γ2=γ2,γ1.\langle\gamma_1,\gamma_2\rangle =-\langle\gamma_2,\gamma_1\rangle.

For a compact surface of genus gg, choose a symplectic basis A1,,Ag,B1,,BgA_1,\ldots,A_g,B_1,\ldots,B_g with

Ai,Aj=0,Bi,Bj=0,Ai,Bj=δij,Bi,Aj=δij.\begin{aligned} \langle A_i,A_j\rangle&=0, & \langle B_i,B_j\rangle&=0,\\ \langle A_i,B_j\rangle&=\delta_{ij}, & \langle B_i,A_j\rangle&=-\delta_{ij}. \end{aligned}

Periods in a different symplectic basis are related by an integral symplectic transformation. A formula written in terms of “the AA period” therefore has no invariant meaning until the cycle basis and orientation are declared.

Intersection numbers govern later discontinuity and wall-crossing formulae. Those chapters state the precise sign convention again, because reversing either a cycle or the order in the pairing reverses the exponent or jump.

There is also an absolute–relative pairing

H1(X^;Z)×H1(X^,P;Z)Z,H_1(\widehat X;\mathbb Z) \times H_1(\widehat X,P;\mathbb Z) \longrightarrow\mathbb Z,

defined by perturbing representatives so their interior crossings are transverse and away from PP. It is the pairing most directly used when a closed WKB cycle acts on an open Voros path. On a punctured surface the absolute pairing can be degenerate—for example, small loops around punctures may lie in its radical—so a “symplectic basis” is asserted only for the nondegenerate quotient or for compact absolute homology.

A differential expression becomes a spectral problem only after its function space and domain have been specified. For a real Sturm–Liouville expression, write

τy=1w(x)[(p(x)y(x))q(x)y(x)],p>0,w>0.\tau y =-\frac1{w(x)} \left[\bigl(p(x)y'(x)\bigr)'-q(x)y(x)\right], \qquad p>0,\quad w>0.

It acts naturally in L2((a,b),w(x) ⁣dx)L^2((a,b),w(x)\,\dd x), but the expression τ\tau alone is not yet an operator. An operator LL also needs a domain D(L)\mathcal D(L) encoding regularity and boundary conditions.

For sufficiently regular ff and gg and interior points a<c<d<ba<c<d<b, the Lagrange identity is

cd[(τf)gfτg]w ⁣dx=[p(fgfg)]cd.\begin{aligned} &\int_c^d \left[ (\tau f)\,\overline g -f\,\overline{\tau g} \right]w\,\dd x\\ &\qquad =\left[ p\left( f\,\overline{g'} -f'\,\overline g \right) \right]_c^d. \end{aligned}

Self-adjoint boundary conditions make the boundary form on the right vanish in the appropriate endpoint limits for every pair in the domain and make the domain equal to that of the adjoint. At a singular endpoint, the individual limits need not exist until the maximal domain and boundary form have been defined. Merely having real coefficients is not enough.

At a singular endpoint, the limit-point/limit-circle alternative decides whether one must impose a separated boundary condition. The endpoint is limit-circle when every solution of τy=λy\tau y=\lambda y for one—and hence every—nonreal λ\lambda is locally square-integrable there; otherwise it is limit-point.

  • in the limit-point case, no boundary condition is imposed there for a separated self-adjoint realization;
  • in the limit-circle case, one boundary condition is required there.

The precise test depends on the differential expression and measure. “Choose the regular solution” is not a substitute when both local solutions lie in the Hilbert space or when a physical boundary condition is not self-adjoint. More general self-adjoint domains can couple the two endpoints.

For a closed densely defined operator LL, the resolvent set consists of λ\lambda for which

(Lλ)1(L-\lambda)^{-1}

exists as a bounded operator on the whole Hilbert space. Its complement is the spectrum. An eigenvalue is a spectral point with a nonzero solution of

(Lλ)y=0,yD(L).(L-\lambda)y=0, \qquad y\in\mathcal D(L).

A self-adjoint operator has real spectrum, but it need not have only eigenvalues. A purely discrete spectrum follows under additional compactness conditions, such as compact resolvent. For nonselfadjoint operators, the spectrum may be complex, eigenvectors need not be orthogonal or complete, and geometric and algebraic multiplicities can differ.

An analytic operator pencil L(λ)L(\lambda) is more general than LλL-\lambda. Its spectral points are values at which L(λ)L(\lambda) fails to be invertible. Before applying the analytic Fredholm theorem, one needs a holomorphic Fredholm family on a fixed domain, normally of index zero, that is invertible for at least one parameter value.

Suppose yL(z,λ)y_L(z,\lambda) and yR(z,λ)y_R(z,\lambda) are solutions selected by the left and right boundary conditions, analytically normalized in λ\lambda. At a fixed ordinary match point zmz_m, their boundary function is

D(λ)=Wr[yL,yR]z=zm.D(\lambda) =\Wr[y_L,y_R]\big|_{z=z_m}.

Abel’s identity determines how this Wronskian varies with zz; in normal form it is independent of zz. In Sturm–Liouville form the naturally constant quantity is p(z)Wr[yL,yR]p(z)\Wr[y_L,y_R]. Suppose each boundary condition selects a nonzero one-dimensional solution subspace. Then

D(λ)=0D(\lambda)=0

means that the two one-dimensional solution spaces coincide, so one solution satisfies both boundary conditions.

Rescaling the boundary solutions by nonvanishing analytic factors gives

yLa(λ)yL,yRb(λ)yR,Da(λ)b(λ)D.y_L\mapsto a(\lambda)y_L, \qquad y_R\mapsto b(\lambda)y_R, \qquad D\mapsto a(\lambda)b(\lambda)D.

The zero set is unchanged, but the function and its derivatives are not. This is why a Jost or Evans function is usually defined only up to a nonvanishing analytic factor unless a normalization is fixed.

For scattering resonances and black-hole quasinormal modes, the selected solutions are typically outgoing, ingoing, or decaying in different asymptotic regions. The resulting problem is generally nonselfadjoint. Resonances are often defined as poles of a meromorphically continued resolvent or scattering matrix; a Wronskian zero represents the same resonance only after that equivalence and the continuation sheet have been established.

Later objectToolkit interpretationExtra data still required
Monodromy matrixParallel transport around a based loopA framing and loop convention
Connection matrixTransport between two normalized framesEndpoint bases and continuation path
Spectral curveBranched cover carrying the globally defined one-form λ\lambdaNormalization, compactification, and punctures
Closed classical or quantum periodIntegral along an absolute cycleSheet, orientation, differential, and regularization
Voros path or cycle symbolExponential of a regularized WKB integral along a relative path or absolute cycleEndpoint treatment, subtraction, and lateral sum
Boundary Wronskian or Jost–Evans functionAnalytic Wronskian of selected solutionsBoundary normalizations and operator interpretation
EigenvalueSpectral point with a vector in the operator domainHilbert space, closed operator, and boundary conditions
Resonance or QNMPole or equivalent boundary zero on a continued sheetContinuation, radiation conditions, and multiplicity convention

Confusing a loop with its drawing. Two curves that look similar in a chosen cut plane can lift to different paths or sheets. Record the based homotopy class and the lift, not only the sketch.

Calling a monodromy matrix invariant. A matrix requires a framed solution basis. The unframed data are a conjugacy class of representations, together with any additional local conjugacy constraints.

Integrating on the base plane. The square root in a WKB differential is single-valued on the spectral cover. A projected contour without sheet and cut-crossing data does not determine the sign of its period.

Replacing monodromy by homology. First homology is the abelianization of the fundamental group. It is ideal for additive period data but discards the ordering information carried by noncommuting monodromy matrices.

Equating an ODE with an operator. The differential expression does not specify a Hilbert space, domain, or boundary conditions. Spectral language is used only after those choices are stated.

1. Monodromy at three punctures. On P1{0,1,}\mathbb P^1\setminus\{0,1,\infty\}, suppose the chosen generators satisfy γ0γ1γ=1\gamma_0\gamma_1\gamma_\infty=1. Express MM_\infty in terms of M0M_0 and M1M_1, and show that a change of framing preserves the relation.

Solution

The representation property gives

M0M1M=I,M=(M0M1)1.M_0M_1M_\infty=I, \qquad M_\infty=(M_0M_1)^{-1}.

Under ΦΦG\Phi\mapsto\Phi G, every matrix changes to M~j=G1MjG\widetilde M_j=G^{-1}M_jG. Hence

M~0M~1M~=G1(M0M1M)G=I.\widetilde M_0\widetilde M_1\widetilde M_\infty =G^{-1}(M_0M_1M_\infty)G =I.

2. Genus of a polynomial spectral curve. Let QQ be a polynomial of degree d1d\geq1 with simple zeros. Compactify and normalize y2=Q(z)y^2=Q(z). Show that

g=d12.g=\left\lfloor\frac{d-1}{2}\right\rfloor.
Solution

Every finite zero is a simple branch point. Infinity is a branch point exactly when dd is odd, so the number of branch points is

r={d,d even,d+1,d odd.r= \begin{cases} d, & d\ \text{even},\\ d+1, & d\ \text{odd}. \end{cases}

Riemann–Hurwitz for a double cover of the sphere gives g=(r2)/2g=(r-2)/2, which is (d1)/2\lfloor(d-1)/2\rfloor in both cases.

3. Absolute or relative? Compactify y2=z21y^2=z^2-1, and then remove the two points +\infty_+ and \infty_- above infinity. Classify a lifted path joining the two branch points and the closed lift obtained by going along the cut on one sheet and returning on the other. Explain what changes when the orientation is reversed.

Solution

The path joining the branch points has nonzero boundary and is a relative class with endpoints in P={1,1}P=\{-1,1\}. Traversing from one branch point to the other on one sheet and returning on the other produces a closed absolute cycle. It separates +\infty_+ from \infty_- and generates H1(Σ^{+,};Z)ZH_1(\widehat\Sigma\setminus\{\infty_+,\infty_-\};\mathbb Z)\cong\mathbb Z. On the compact genus-zero cover it becomes null homologous. Reversing either path negates its class in the relevant absolute or relative group and therefore negates the integral of any fixed one-form along it.

4. Boundary-function normalization. If D(λ)=Wr[yL,yR]D(\lambda)=\Wr[y_L,y_R], show that rescaling yLy_L and yRy_R by nonvanishing analytic functions preserves the zeros of DD. Why can such a rescaling still matter in a residue or norming-constant formula?

Solution

Bilinearity gives

D~(λ)=a(λ)b(λ)D(λ).\widetilde D(\lambda) =a(\lambda)b(\lambda)D(\lambda).

Because aa and bb never vanish, DD and D~\widetilde D have the same zeros with the same multiplicities. Their derivatives at a zero λn\lambda_n satisfy

D~(λn)=a(λn)b(λn)D(λn),\widetilde D'(\lambda_n) =a(\lambda_n)b(\lambda_n)D'(\lambda_n),

so any residue or norming formula involving D(λn)D'(\lambda_n) must use the same normalization as its numerator.

5. One expression, two operators. On [0,π][0,\pi], compare L= ⁣d2/ ⁣dx2L=-\dd^2/\dd x^2 with Dirichlet conditions y(0)=y(π)=0y(0)=y(\pi)=0 and with Neumann conditions y(0)=y(π)=0y'(0)=y'(\pi)=0. Derive boundary functions and their spectra.

Solution

Write k=λk=\sqrt{\lambda}. For the Dirichlet problem, the solution normalized by y(0)=0y(0)=0 and y(0)=1y'(0)=1 is

yD(x,λ)=sin(kx)k,y_D(x,\lambda)=\frac{\sin(kx)}{k},

with its removable value at k=0k=0. A boundary function is

DD(λ)=sin(πλ)λ.D_D(\lambda) =\frac{\sin(\pi\sqrt\lambda)}{\sqrt\lambda}.

Its zeros are λ=n2\lambda=n^2 for n=1,2,n=1,2,\ldots. For the Neumann problem, use yN(x,λ)=cos(kx)y_N(x,\lambda)=\cos(kx) and impose the right derivative condition:

DN(λ)=yN(π,λ)=λsin(πλ).D_N(\lambda) =y_N'(\pi,\lambda) =-\sqrt\lambda\sin(\pi\sqrt\lambda).

Its zeros are λ=n2\lambda=n^2 for n=0,1,2,n=0,1,2,\ldots. The differential expression is the same, but the domains—and therefore the operators and spectra—are different.