Skip to content

Hopf fibers and transverse geometry

A circle drawn on a sphere is not automatically a fiber, and a closed curve on the base need not lift to a closed horizontal curve. This laboratory makes both distinctions calculable. You will draw two fibers, check the geometry without choosing local coordinates, and transport a point around a base loop that passes through the missing pole of a familiar section.

Prepare with circle bundles and the Hopf fibration and almost contact metric structures. Those lessons supply the analytic construction and its geometric meaning. Here we turn their formulas into a portable calculation, including tests that fail when a component is omitted.

Write a point of the unit three-sphere as y=(y0,y1,y2,y3)y=(y_0,y_1,y_2,y_3), with yTy=1y^Ty=1, and assemble two complex coordinates z1=y0+iy1z_1=y_0+iy_1, z2=y2+iy3z_2=y_2+iy_3. The column zz is not a supersymmetry parameter. The physical round metric is g=ℓ2dyTdyg=\ell^2dy^Tdy on tangent vectors, with radius ℓ>0\ell>0.

The Hopf map and the circle action are

n=(2Re⁡zˉ1z2,  2Im⁡zˉ1z2,∣z1∣2−∣z2∣2),z⟼eiχz,χ∼χ+2π.(LAB05.1)n=(2\operatorname{Re}\bar z_1z_2,\;2\operatorname{Im}\bar z_1z_2, |z_1|^2-|z_2|^2),\qquad z\longmapsto e^{i\chi}z,\qquad \chi\sim\chi+2\pi. \tag{LAB05.1}

Every orbit lies over one point n∈S2n\in S^2. Increasing χ\chi fixes its orientation. We orient S3S^3 as the outward boundary of the unit ball in the ordered real coordinates (y0,y1,y2,y3)(y_0,y_1,y_2,y_3), and S2S^2 outward in its three nn coordinates.

The normalized principal connection is

Θ=−iz†dz=y0dy1−y1dy0+y2dy3−y3dy2,Θ(ξ0)=1,ξ0=(−y1,y0,−y3,y2).(LAB05.2)\Theta=-iz^\dagger dz =y_0dy_1-y_1dy_0+y_2dy_3-y_3dy_2, \qquad \Theta(\xi_0)=1, \quad \xi_0=(-y_1,y_0,-y_3,y_2). \tag{LAB05.2}

The metric-unit objects are instead η=ℓΘ\eta=\ell\Theta and ξ=ξ0/ℓ\xi=\xi_0/\ell. A fiber has circumference 2πℓ2\pi\ell. In P (5.11), the round limit has r=ℓr=\ell, h=1h=1 and the paper’s angular coordinate has period 4π4\pi: locally ψ=2χN\psi=2\chi_N, and the paper’s connection bracket is 2Θ2\Theta. Missing this factor changes both the metric and the apparent bundle number.

This is a calculation in three-dimensional positive metric geometry. We choose no supergravity sources or dynamical gauge theory, so no preserved supercharge count or gauge rank is assigned by the laboratory.

Use northern and southern sections

sN=(cos⁡(θ/2)eiϕsin⁡(θ/2)),sS=(e−iϕcos⁡(θ/2)sin⁡(θ/2)).(LAB05.3)s_N=\begin{pmatrix}\cos(\theta/2)\\e^{i\phi}\sin(\theta/2)\end{pmatrix}, \qquad s_S=\begin{pmatrix}e^{-i\phi}\cos(\theta/2)\\\sin(\theta/2)\end{pmatrix}. \tag{LAB05.3}

They obey sN=eiϕsSs_N=e^{i\phi}s_S. For the same total-space point, therefore, χS=χN+ϕ\chi_S=\chi_N+\phi modulo 2π2\pi. Pulling back (LAB05.2) gives

Θ=dχN+1−cos⁡θ2dϕ=dχS−1+cos⁡θ2dϕ,F=dΘ=sin⁡θ2dθ∧dϕ.(LAB05.4)\Theta=d\chi_N+\frac{1-\cos\theta}{2}d\phi =d\chi_S-\frac{1+\cos\theta}{2}d\phi, \qquad \mathcal F=d\Theta=\frac{\sin\theta}{2}d\theta\wedge d\phi. \tag{LAB05.4}

Here the curvature on the right is understood as a form on the base and then pulled back. Its integral is exactly ∫S2F=2π\int_{S^2}\mathcal F=2\pi, so the first Chern number is +1+1. The transition function has the same unit winding. The program integrates this formula symbolically; it does not infer an integer by rounding a sampled flux.

The two sections are useful, but neither covers the entire base. The coordinate column zz in the ambient C2\mathbb C^2 description has no such local-section singularity.

Let J0J_0 denote multiplication by ii on both complex coordinates: J0y=ξ0J_0y=\xi_0. Define the tangent and horizontal projectors

Π=I4−yyT,Πh=Π−ξ0ξ0T,Φ=−J0Πh.(LAB05.5)\Pi=I_4-yy^T,\qquad \Pi_h=\Pi-\xi_0\xi_0^T, \qquad \Phi=-J_0\Pi_h. \tag{LAB05.5}

“Horizontal” means annihilated by Θ\Theta. At y=(1,0,0,0)y=(1,0,0,0) the oriented tangent basis is (E1,E2,E3)(E_1,E_2,E_3), with EjE_j in the yjy_j direction. Thus ΦE2=−E3\Phi E_2=-E_3 and ΦE3=E2\Phi E_3=E_2, precisely the convention of the prerequisite lesson.

Because J02=−I4J_0^2=-I_4 and J0J_0 preserves the horizontal plane, multiplication gives

Φ2=−Πh,Φξ0=0,ξ0TΦ=0,ΦTΦ=Πh.(LAB05.6)\Phi^2=-\Pi_h,\quad \Phi\xi_0=0,\quad \xi_0^T\Phi=0,\quad \Phi^T\Phi=\Pi_h. \tag{LAB05.6}

Restricted to the tangent space, these are P (B.1)–(B.2), including metric compatibility. The radial direction is an extra zero direction of our ambient four-by-four matrix, not an extra tangent direction on S3S^3.

The script also differentiates the Hopf map itself. For a horizontal tangent vector XX it finds ∣dn(X)∣2=4∣X∣2|dn(X)|^2=4|X|^2, confirming the base radius ℓ/2\ell/2. Finally, the circle flow Rχ=eχJ0R_\chi=e^{\chi J_0} obeys

Φ(Rχy)=RχΦ(y)Rχ−1.(LAB05.7)\Phi(R_\chi y)=R_\chi\Phi(y)R_\chi^{-1}. \tag{LAB05.7}

This says that the transverse rotation is preserved along fibers. Its infinitesimal version is the full tensor Lie derivative Lξ0Φ=0\mathcal L_{\xi_0}\Phi=0, including the index-rotation terms. It is stronger than the projected condition in P (B.8). Meanwhile Θ∧dΘ\Theta\wedge d\Theta never vanishes: the horizontal planes are not tangent to surfaces. A transverse complex structure on the base and nonintegrable horizontal planes are compatible statements.

Choose the two base points n=(±1,0,0)n=(\pm1,0,0) and their lifts z=(1,±1)/2z=(1,\pm1)/\sqrt2. Their fibers are the literal curves eiχ(1,±1)/2e^{i\chi}(1,\pm1)/\sqrt2. The script projects them to R3\mathbb R^3 by

r=(−y1,y2,y3)1−y0.(LAB05.8)\boldsymbol r=\frac{(-y_1,y_2,y_3)}{1-y_0}. \tag{LAB05.8}

The omitted point is (1,0,0,0)(1,0,0,0), which belongs to neither displayed fiber. The minus sign makes this chart preserve the declared outward orientation. The resulting circles have centers (0,±1,0)(0,\pm1,0), radius 2\sqrt2, and planes r3=∓r1r_3=\mp r_1.

Blue and orange Hopf fibers over opposite equatorial base points form two linked circles in a stereographic view. Arrows mark increasing common phase; gaps distinguish depth at planar crossings.

The view is a fixed orthographic projection of those three-dimensional circles, at ℓ=1\ell=1. Its apparent crossings are not intersections of fibers. Distances in the picture are not round-sphere distances. Both circles close at χ=2π\chi=2\pi and have linking number +1+1 with the declared orientations. One can establish the sign without a plot: the second circle pierces the oriented planar disk of the first once positively. A separate Gauss linking integral checks the numerical parametrization.

Lift a loop that crosses a section’s missing point

Section titled “Lift a loop that crosses a section’s missing point”

Let c∈(−1,1)c\in(-1,1) be a latitude cosine and set s=1−c2s=\sqrt{1-c^2}. Use the oriented orthonormal triple

a=(s,0,−c),b=(c,0,s),d=a×b=(0,−1,0),n(t)=ca+s(bcos⁡t+dsin⁡t),0≤t≤2π.(LAB05.9)a=(s,0,-c),\quad b=(c,0,s),\quad d=a\times b=(0,-1,0), \qquad n(t)=ca+s(b\cos t+d\sin t),\quad 0\le t\le2\pi. \tag{LAB05.9}

This tilted latitude passes through the south pole at t=πt=\pi. Its cap about the axis aa has positive area 2π(1−c)2\pi(1-c) on the unit base sphere, hence curvature flux π(1−c)\pi(1-c).

For transport, use the ordinary Pauli matrices σ=(σ1,σ2,σ3)\sigma=(\sigma_1,\sigma_2,\sigma_3) with σ3=diag⁡(1,−1)\sigma_3=\operatorname{diag}(1,-1). They encode the coordinate identity z†σz=nz^\dagger\sigma z=n; they are not a change to the site’s spacetime gamma-matrix convention. The rank-one projector onto the coordinate line is Pn=(I2+n⋅σ)/2P_n=(I_2+n\cdot\sigma)/2. A globally regular horizontal-lift equation is

z˙=P˙nz=12n˙⋅σ z,z(0)=(sc).(LAB05.10)\dot z=\dot P_n z=\frac12\dot n\cdot\sigma\,z, \qquad z(0)=\begin{pmatrix}s\\c\end{pmatrix}. \tag{LAB05.10}

To see why this works, differentiate Pn2=PnP_n^2=P_n to get PnP˙nPn=0P_n\dot P_nP_n=0. On Pnz=zP_nz=z, (LAB05.10) preserves that constraint and gives z†z˙=0z^\dagger\dot z=0, preserving the norm and imposing Θ(z˙)=0\Theta(\dot z)=0. All coefficients stay smooth at t=πt=\pi.

There is also a check that does not integrate this equation numerically. Rotating about aa gives U(t)=cos⁡(t/2)I2−isin⁡(t/2)a⋅σU(t)=\cos(t/2)I_2-i\sin(t/2)a\cdot\sigma. Its lift has z†z˙=−ic/2z^\dagger\dot z=-ic/2, so multiplication by a compensating phase gives

zH(t)=eict/2U(t)z(0),zH(2π)=e−iπ(1−c)z(0).(LAB05.11)z_H(t)=e^{ict/2}U(t)z(0),\qquad z_H(2\pi)=e^{-i\pi(1-c)}z(0). \tag{LAB05.11}

The return phase agrees with e−i∫Fe^{-i\int\mathcal F}. It is the phase of a horizontally transported point in the bundle. A charge-mm field component with convention D=d−imaD=d-ima has the opposite-sign parallel-transport phase e+im∮ae^{+im\oint a}; do not interchange these objects.

For the default c=1/2c=1/2, the return phase is −i-i. The base closes once while its horizontal lift does not. This is different from the fibers drawn above, whose base point never moves.

Download the complete lab (ZIP) with the program, pinned dependencies, reference output and run instructions.

Download the Python script, pinned requirements, expected output, and the reproducible SVG into one directory. The tested environment is Python 3.13.7, SymPy 1.13.3 and mpmath 1.3.0. No plotting package, network access, or site checkout is needed after installation.

python -m pip install -r requirements.txt
python hopf_fibration.py
python hopf_fibration.py --figure hopf-fibers.svg
python hopf_fibration.py --cosine 1/3
Mathematical objectVariable or functionWhat it checks
J0,Π,Πh,ΦJ_0,\Pi,\Pi_h,\PhiJ, P, H, Phi in exact_checksExact identities modulo yTy=1y^Ty=1; H here is a projector, not a supergravity source
n,n˙,zHn,\dot n,z_Hbase, velocity, exact in dataSmooth loop and an analytic reference curve
(LAB05.10)rhs inside integrateFour-stage Runge–Kutta transport
(LAB05.8)point inside figureNative SVG generated from the actual fibers

The default run passes 23 grouped exact assertions. The transport uses 50 decimal digits and does not renormalize after steps, so norm drift remains observable. The maximum vector error over all step endpoints is

Steps around the loopMaximum vector errorMaximum norm-squared drift
642.61152×10−72.61152\times10^{-7}4.54939×10−74.54939\times10^{-7}
1281.61382×10−81.61382\times10^{-8}2.80134×10−82.80134\times10^{-8}
2561.00294×10−91.00294\times10^{-9}1.73771×10−91.73771\times10^{-9}

The error ratios are 16.1816.18 and 16.0916.09, approaching the factor 242^4 expected for fourth-order global error. More decimal digits alone do not remove the finite-step error. The tiny approximately 10−5210^{-52} real part in the printed exact −i-i phase is rounding; it is not a geometric correction. The displayed error is measured at step endpoints, not a certified continuous-time bound.

The exact tensor tests, analytic cap flux and rotating-frame solution address different possible errors. A further check replaces the loop by short geodesic edges, transports by successive normalized projectors Pnj+1P_{n_{j+1}}, and compares with the polygon’s oriented solid angle. That method converges to the same smooth-loop phase at second order, rather than fourth: it approximates a different curve at finite resolution.

1. A transverse test that misses a vertical error

Section titled “1. A transverse test that misses a vertical error”

At y=(1,0,0,0)y=(1,0,0,0), let u=E2u=E_2 and change the tensor to Φε=Φ+εξ0uT\Phi_\varepsilon=\Phi+\varepsilon\xi_0u^T. Does testing only ΠhΦεΠh\Pi_h\Phi_\varepsilon\Pi_h detect the error? Evaluate the full vertical and norm tests for ε=1/3\varepsilon=1/3. What would you add to a test suite?

Hint 1

The added output points along ξ0\xi_0, which the left horizontal projector removes.

Hint 2

Apply the modified tensor to uu. Compare both its Θ\Theta component and its squared length with the required values.

Solution

Since Πhξ0=0\Pi_h\xi_0=0, the projected test gives exactly ΠhΦεΠh=Φ\Pi_h\Phi_\varepsilon\Pi_h=\Phi for every ε\varepsilon. But

Φεu=−E3+εE1,Θ(Φεu)=ε,∣Φεu∣2=1+ε2.\Phi_\varepsilon u=-E_3+\varepsilon E_1,\qquad \Theta(\Phi_\varepsilon u)=\varepsilon,\qquad |\Phi_\varepsilon u|^2=1+\varepsilon^2.

At ε=1/3\varepsilon=1/3 the last two results are 1/31/3 and 10/910/9, instead of zero and one. The modified tensor fails to map horizontal vectors into horizontal vectors and fails metric compatibility. Test the full left kernel ΘΦ=0\Theta\Phi=0 and ΦTΦ=Πh\Phi^T\Phi=\Pi_h, in addition to projected identities. The released script includes this negative control: its successful projected check is deliberately insufficient.

For the default c=1/2c=1/2, show that the north section has different limiting values as the loop approaches t=πt=\pi from either side, although zH(t)z_H(t) is smooth. Give the limiting azimuths and show how the transition to sSs_S removes the discrepancy. Is this a change of the connection or of its description?

Hint 1

Near t=πt=\pi, use n1=3(1+cos⁡t)/4n_1=\sqrt3(1+\cos t)/4 and n2=−3sin⁡t/2n_2=-\sqrt3\sin t/2.

Hint 2

At the south pole the north section tends to (0,eiϕ)(0,e^{i\phi}), whereas the south section tends to (0,1)(0,1). Remember χS=χN+ϕ\chi_S=\chi_N+\phi.

Solution

Writing t=π+δt=\pi+\delta, we have n1=3δ2/8+O(δ4)n_1=\sqrt3\delta^2/8+O(\delta^4) and n2=3δ/2+O(δ3)n_2=\sqrt3\delta/2+O(\delta^3). Thus the azimuth tends to −π/2-\pi/2 on approach and +π/2+\pi/2 on departure. The north section tends to (0,−i)(0,-i) and (0,+i)(0,+i): the two local representatives differ by a sign.

The exact lift is finite and smooth. Direct substitution in (LAB05.11) gives zH(π)=(0,e−iπ/4)z_H(\pi)=(0,e^{-i\pi/4}). Its north fiber coordinate compensates for the changing section. In the south description, sS→(0,1)s_S\to(0,1) and eiχS→e−iπ/4e^{i\chi_S}\to e^{-i\pi/4} on both sides, because χS=χN+ϕ\chi_S=\chi_N+\phi. The local potentials differ by aN−aS=dϕa_N-a_S=d\phi, while dχN+aN=dχS+aSd\chi_N+a_N=d\chi_S+a_S. The global connection is unchanged. An algorithm that treats the north section’s jump as a jump of the transported point has confused a coordinate representative with the point itself.

Change the latitude cosine to c=1/3c=1/3 and run the final command above. Predict the return phase and the least number of complete base circuits needed to close the horizontal lift. Compare with the default. What happens for irrational cc? Determine the length of the closed default lift and compare it with a single Hopf fiber.

Hint 1

After qq repetitions, the accumulated phase is e−iqπ(1−c)e^{-iq\pi(1-c)}. Ask when this is one.

Hint 2

The base radius is ℓ/2\ell/2. A horizontal lift has the same length as its base curve in that metric; a single latitude circuit has length πℓ1−c2\pi\ell\sqrt{1-c^2}.

Solution

For c=1/3c=1/3 the phase is e−2πi/3=−1/2−i3/2e^{-2\pi i/3}=-1/2-i\sqrt3/2, so three circuits close the lift. At c=1/2c=1/2 the phase is −i-i and four circuits are required. In general the condition is q(1−c)∈2Zq(1-c)\in2\mathbb Z; irrational cc admits no positive integer qq. This is an exact arithmetic condition, not a numerical decision based on a phase close to one.

The default one-circuit length is π3ℓ/2\pi\sqrt3\ell/2, so its four-circuit closed lift has length 2π3ℓ2\pi\sqrt3\ell. A Hopf fiber has length 2πℓ2\pi\ell and projects to a point. The longer closed horizontal curve projects to a moving latitude four times. Closing a curve does not make it a fiber, and changing this loop does not change the bundle number.

The useful outcome is a distinction among global geometry, local descriptions and numerical representations. The smooth fibers and horizontal lift are global; the local section omits a pole; the finite-step trajectory approximates a differentially defined curve. These are separate facts, and the checks should keep them separate.