shellRadius_ratio

formal source

  29
  30noncomputable def shellRadius (k : ℕ) : ℝ := phi ^ k
  31
  32theorem shellRadius_ratio (k : ℕ) : shellRadius (k + 1) / shellRadius k = phi := by
  33  unfold shellRadius
  34  have hpos : (0 : ℝ) < phi ^ k := pow_pos phi_pos k
  35  rw [div_eq_iff hpos.ne', pow_succ]
  36  ring
  37
  38theorem shellRadius_pos (k : ℕ) : 0 < shellRadius k := pow_pos phi_pos k
  39
  40structure SolvationShellCert where
  41  five_shells : Fintype.card SolvationShell = 5
  42  phi_ratio : ∀ k, shellRadius (k + 1) / shellRadius k = phi
  43  radius_always_pos : ∀ k, 0 < shellRadius k
  44
  45noncomputable def solvationShellCert : SolvationShellCert where
  46  five_shells := solvationShell_count
  47  phi_ratio := shellRadius_ratio
  48  radius_always_pos := shellRadius_pos
  49
  50end IndisputableMonolith.Chemistry.SolvationShellsFromConfigDim