Chemistry Crystal Structure Hcp Ratio Near Phi

The HCP ratio and the golden mean

Hexagonal close-packed (HCP) is one of the three common ways atoms pack into a crystal. In the ideal version of this packing, each atom touches twelve neighbors, and the structure has a characteristic shape: a hexagonal base with a height that is not arbitrary. The ratio of that height to the side of the hexagon is called the c/a ratio. For a perfect HCP crystal, geometry fixes this ratio at the square root of 8/3, which is approximately 1.633. This value is not an empirical accident; it follows from the condition that the spheres are as close as possible.

The number 1.633 is close to the golden ratio, often written as φ and equal to (1 + √5)/2 ≈ 1.618. The two numbers differ by about 0.015, which is less than one percent of either value. This proximity has drawn attention because the golden ratio appears in many natural patterns, from plant growth to spiral shells. In a crystal, however, the c/a ratio is not a matter of aesthetics; it determines how the layers of atoms sit on top of one another. Real HCP metals such as magnesium and zinc have c/a ratios that differ from the ideal value, and those deviations reflect real bonding forces between atoms.

The machine-checked library of formal theorems in the Recognition Science framework includes a statement about this ratio. The declaration hcp_ratio_near_phi establishes a precise, limited fact: the absolute difference between the ideal HCP ratio and the golden ratio is less than 0.03. That is a theorem in the framework's library, proved by computation from the definitions of the two numbers. The framework also defines the ideal HCP ratio as the square root of 8/3 and proves that this value lies between 1.63 and 1.64. These are formal statements about the real numbers, not measurements of any physical crystal.

What the declaration does not claim is just as important. It does not say that the golden ratio causes HCP crystals to form, nor that real metals must have c/a near φ. The closeness is a numerical coincidence within a tolerance, not a physical law. The framework's library also contains definitions for packing efficiency, coordination numbers, and energy scales, but those are separate declarations with their own proofs. The hcp_ratio_near_phi theorem alone says nothing about why atoms choose HCP over other structures, and it does not predict the c/a ratio of any specific element. It is a statement about two real numbers, and nothing more.

For a reader meeting HCP for the first time, the useful fact is the classical one: the ideal c/a ratio is √(8/3) ≈ 1.633, a fixed geometric value. The golden ratio connection is a curiosity that the framework's formal library records as a true inequality. The distinction matters because a theorem about numbers is not a theory of crystals. The framework's own documentation treats the closeness as a numerical observation, not as a mechanism for why HCP exists.

/-- The ideal HCP ratio is close to φ ≈ 1.618.
    √(8/3) ≈ 1.633, φ ≈ 1.618, difference ≈ 0.015.
    Using available bounds: 1.63 < √(8/3) < 1.64, 1.61 < φ < 1.62.
    This gives |√(8/3) - φ| < 1.64 - 1.61 = 0.03. -/
theorem hcp_ratio_near_phi : |idealHCPRatio - phi| < 0.03 := by
  simp only [idealHCPRatio]
  -- First establish that √(8/3) > φ, so |√(8/3) - φ| = √(8/3) - φ
  have h_phi_lt : phi < 1.62 := phi_lt_onePointSixTwo
  have h_163_lt_sqrt : (1.63 : ℝ) < Real.sqrt (8/3) := by
    rw [Real.lt_sqrt (by norm_num : (0 : ℝ) ≤ 1.63)]
    norm_num
  have h_sqrt_gt_phi : Real.sqrt (8/3) > phi := by linarith
  rw [abs_of_pos (by linarith : Real.sqrt (8/3) - phi > 0)]
  -- Now show √(8/3) - φ < 0.03
  -- √(8/3) < 1.64 and φ > 1.61, so √(8/3) - φ < 1.64 - 1.61 = 0.03
  have h_sqrt_lt : Real.sqrt (8/3) < 1.64 := by
    rw [Real.sqrt_lt' (by norm_num : (0 : ℝ) < 1.64)]
    norm_num
  have h_phi_gt : phi > 1.61 := phi_gt_onePointSixOne
  linarith

/-- Ideal c/a ratio for HCP: √(8/3) ≈ 1.633. -/
def idealHCPRatio : ℝ := Real.sqrt (8/3)

/-- The ideal HCP c/a ratio is approximately 1.633. -/
theorem ideal_hcp_ratio_value : 1.63 < idealHCPRatio ∧ idealHCPRatio < 1.64 := by
  -- √(8/3) ≈ 1.6329931..., so 1.63 < √(8/3) < 1.64
  -- We verify: 1.63² = 2.6569 < 8/3 ≈ 2.6667 < 2.6896 = 1.64²
  simp only [idealHCPRatio]
  have h_sq_lo : (1.63 : ℝ)^2 < 8/3 := by norm_num
  have h_sq_hi : (8 : ℝ)/3 < (1.64 : ℝ)^2 := by norm_num
  constructor
  · -- 1.63 < √(8/3) ⟺ 1.63² < 8/3 (for positive values)
    rw [Real.lt_sqrt (by norm_num : (0 : ℝ) ≤ 1.63)]
    exact h_sq_lo
  · -- √(8/3) < 1.64 ⟺ 8/3 < 1.64² (for positive values)
    rw [Real.sqrt_lt' (by norm_num : (0 : ℝ) < 1.64)]
    exact h_sq_hi

What this page does not claim

The golden ratio causes HCP crystal formation. Real HCP metals have c/a ratios near the golden ratio. The hcp_ratio_near_phi theorem predicts the crystal structure of any element.

Verify this page

Every tagged claim above names its theorem. To check one yourself rather than trust this page, elaborate the source module with Lean 4 and audit its axiom basis:

$ lake env lean IndisputableMonolith/Chemistry/CrystalStructure.lean
expected axiom basis: [propext, Classical.choice, Quot.sound] (the Lean kernel's standard three; no RS-specific axioms)

A page whose claims cannot be reproduced this way does not ship. In production, every anchor links to the exact declaration in the public source release, and this block carries the build receipt for the page itself.

Derived articles

This page is generated by a question-recursion engine: the questions its answers raise become the next pages. The current agenda, with open targets marked red: