Twisting formula: a push-off winding times around a knot in links it more times
ProvedHryniewiczCriterion.gaussLinkingIntegral_twisted_pushOffLet be a , -periodic loop, injective modulo , and let be , -periodic vector fields along with
so that span a positively oriented complement of in (oriented as the boundary of the ball). Let and be with and . For put
Then there are a pole and such that for the push-off is well defined and disjoint from , misses , , , and the Gauss linking integrals from satisfy
Proof idea: take off and so small that all push-offs stay in a thin tube around away from . The homotopy , stays in the tube minus , so by homotopy invariance of the Gauss integral we may take , . Pulled back to the torus , the Biot–Savart form of is closed, so the curve of class has integral equal to that of the class plus times the meridian integral. The meridian integral is (Ampère's law), as the rescaling limit shows: the knot becomes its tangent line and the meridian a circle around it, whose Gauss integral is . The sign is fixed by the orientation hypothesis.
import Definitions.Def_HryniewiczCriterion_SelfLinking open HryniewiczCriterion open scoped ContDiff
theorem HryniewiczCriterion.gaussLinkingIntegral_twisted_pushOff (γ e₁ e₂ : ℝ → R4)
(hγ : ContDiff ℝ 2 γ) (he₁ : ContDiff ℝ 2 e₁) (he₂ : ContDiff ℝ 2 e₂)
(hγper : ∀ s, γ (s + 1) = γ s) (he₁per : ∀ s, e₁ (s + 1) = e₁ s)
(he₂per : ∀ s, e₂ (s + 1) = e₂ s) (hγunit : ∀ s, euclidNorm (γ s) = 1)
(hinj : ∀ s t, γ s = γ t → ∃ n : ℤ, t = s + n)
(hdet : ∀ s, 0 < Matrix.det (Matrix.of ![γ s, deriv γ s, e₁ s, e₂ s]))
(r θ : ℝ → ℝ) (hr : ContDiff ℝ 2 r) (hθ : ContDiff ℝ 2 θ) (hr0 : ∀ s, 0 < r s)
(hrper : ∀ s, r (s + 1) = r s) (k : ℤ) (hθper : ∀ s, θ (s + 1) = θ s + 2 * Real.pi * k) :
∃ N : R4, euclidNorm N = 1 ∧ ∃ ε₀ : ℝ, 0 < ε₀ ∧ ∀ ε : ℝ, 0 < ε → ε < ε₀ →
(∀ s, γ s + ε • e₁ s ≠ 0) ∧
(∀ s t, γ s ≠ radialNormalize (γ t + ε • e₁ t)) ∧
(∀ s, γ s ≠ N ∧ radialNormalize (γ s + ε • e₁ s) ≠ N ∧
radialNormalize (γ s + ε • (r s • (Real.cos (θ s) • e₁ s + Real.sin (θ s) • e₂ s))) ≠ N) ∧
gaussLinkingIntegral N γ
(fun s => radialNormalize (γ s + ε • (r s • (Real.cos (θ s) • e₁ s + Real.sin (θ s) • e₂ s)))) =
gaussLinkingIntegral N γ (fun s => radialNormalize (γ s + ε • e₁ s)) + k := by sorry