MCF-stack · Numerical experiment supplied by Ao Sun

From a perturbed Angenent torus to another genus-one shrinker

Mean curvature flow from the torus perturbed in its first angular mode (ε = 0.2), from the Chu–Sun numerical note. The surface is shown rescaled about the singular point, so a flow that approaches a self-shrinker appears to settle down. Drag to rotate, scroll or pinch to zoom, and move the time slider or press ▶.

Readouts at this time

Distance to candidateGaussian RMS, |Y| < 3–
Shrinker residualGaussian RMS of H·ν + ⟨Y,ν⟩/2, |Y| < 3–
Gaussian density Fabout the singular point (c, T)–

Display

How close is the flow to a self-shrinker?

Log scale. Blue: distance to the independent stationary candidate. Orange: shrinker-equation residual. Click or drag to set the time.

Gaussian density about the singular point

Falls from near the torus value (1.8512) toward the candidate’s 1.7143. Computed on the part of each stage’s box away from its edge, so small wiggles reflect truncation.

What am I looking at?

The flow starts at t = 0 from the rotationally symmetric shrinking torus perturbed by a normal graph ε[−cos(πs) φ₀ + sin(πs) φ₁]ν with ε = 0.2 and s ≈ 0.50041. Here φ₀ is the positive principal unstable Jacobi eigenfunction and φ₁ ∝ e(r²+z²)/4 cos θ / r is the explicit first-angular-mode eigenfunction chosen in the note. The one tuned parameter s sits between “inward tube collapse” (smaller s) and “outward hole closure” (larger s), the interpolation mechanism of Chu–Sun’s Theorem 1.1. The flow develops a genus-one singularity on the x-axis at time T ≈ 0.90768.

The picture shows the rescaled surface Y = (X − c)/√(T − t) with c ≈ (−0.3345, 0, 0). The slider runs in τ = −ln(T − t), which stretches the final instants so that all eight spatial zoom stages of the computation are visible. The torus axis is vertical (z); the surface is mirror-symmetric in y and z, which is why it is drawn by reflecting one quarter.

Where the pictures and numbers come from

I replayed the note’s exact ε = 0.2 critical trajectory (eight adaptive spatial zooms) with the archive’s own level-set solver and saved the field at 129 times equally spaced in τ. The replayed stage fields match the archived ones to within 2×10⁻⁴ of a grid cell. Surfaces are extracted from the unsmoothed scalar field (marching cubes).

The residual and distance use the note’s own estimators (quartic Savitzky–Golay derivatives with window 7; Gaussian-weighted over the ball |Y| < 3; face-centroid distances to the stationary mesh). At the last stage my common-frame residual is 0.00844, matching the 0.008439 reported in the note for the same fitted centre and time. All frames use that one common (c, T); per-frame fits shift the numbers slightly.

Caveats

Data: Chu–Sun numerical note and companion archive (commit aa66772). Rendering with three.js; frames decoded in the browser.