Stability Counter Example not potential Is Stable
TwoHiggsDoublet.stabilityCounterExample_not_potentialIsStable
Plain-language statement
The potential stabilityCounterExample is not stable.
Exact Lean statement
lemma stabilityCounterExample_not_potentialIsStable :
¬ PotentialIsStable .stabilityCounterExampleFormal artifact
Lean source
lemma stabilityCounterExample_not_potentialIsStable : ¬ PotentialIsStable .stabilityCounterExample := by simp [PotentialIsStable] intro c /- The angle t and properties thereof. -/ let t := Real.arctan (2 * Real.sqrt (|c| + 1))⁻¹ have t_pos : 0 < t := arctan_pos.mpr (by positivity) have t_le_pi_div_2 : t ≤ Real.pi / 2 := by simpa [t] using le_of_lt <| arctan_lt_pi_div_two ((√(|c| + 1))⁻¹ * 2⁻¹) have t_ne_zero : t ≠ 0 := t_pos.ne' have sin_t_pos : 0 < sin t := Real.sin_pos_of_pos_of_lt_pi t_pos (by nlinarith [Real.pi_pos]) have cos_t_pos : 0 < cos t := cos_arctan_pos (2 * Real.sqrt (|c| + 1))⁻¹ have t_mul_sin_t_nonneg : 0 ≤ 2 * t * sin t - t ^ 2 := by have hs := Real.mul_le_sin t_pos.le t_le_pi_div_2 rw [div_mul_eq_mul_div, div_le_iff₀ Real.pi_pos] at hs nlinarith [mul_le_mul_of_nonneg_left hs t_pos.le, Real.pi_le_four] /- The Two Higgs doublet violating stability. The two Higgs doublet is constructed so that for the gram vector `v` we have: - `v₀ = cos t/(2 * t * (sin t)^2)` - `v₁/v₀ = (1 - t * sin t)` - `v₂/v₀ = - t * cos t` - `v₃ = 0` -/ let H : TwoHiggsDoublet := { Φ1 := !₂[√(cos t/(4 * t * (sin t)^2)), 0] Φ2 := √(cos t/(4 * t * (sin t)^2)) • !₂[1 - t * sin t - Complex.I * t * cos t, √(2 * t * sin t - t ^ 2)] } have Φ1_norm_sq : ‖H.Φ1‖ ^ 2 = cos t/(4 * t * (sin t)^2) := by simp [H, PiLp.norm_sq_eq_of_L2] rw [sq_sqrt (by positivity)] have Φ2_norm_sq : ‖H.Φ2‖ ^ 2 = cos t/(4 * t * (sin t)^2) := by simp [H, norm_smul, mul_pow] rw [sq_sqrt (by positivity)] simp [PiLp.norm_sq_eq_of_L2] rw [sq_sqrt (by positivity)] have h0 : ‖1 - ↑t * Complex.sin ↑t - Complex.I * ↑t * Complex.cos ↑t‖ ^ 2 = 1 + t ^ 2 - 2 * t * sin t := by rw [← Complex.normSq_eq_norm_sq, Complex.normSq_apply] simp only [Complex.sub_re, Complex.sub_im, Complex.mul_re, Complex.mul_im, Complex.I_re, Complex.I_im, Complex.ofReal_re, Complex.ofReal_im, Complex.one_re, Complex.one_im, Complex.sin_ofReal_re, Complex.sin_ofReal_im, Complex.cos_ofReal_re, Complex.cos_ofReal_im] nlinarith [Real.sin_sq_add_cos_sq t] rw [h0] field_simp ring have Φ1_inner_Φ2 : ⟪H.Φ1, H.Φ2⟫_ℂ = Complex.ofReal (cos t/(4 * t * (sin t)^2) * (1 - t * sin t)) + Complex.I * Complex.ofReal (cos t/(4 * t * (sin t)^2) * (- t * cos t)) := by simp [H, PiLp.inner_apply] trans Complex.ofReal ((√(cos t / (4 * t * sin t ^ 2))) ^ 2) * (1 - ↑t * Complex.sin ↑t - Complex.I * ↑t * Complex.cos ↑t) · simp ring rw [sq_sqrt (by positivity)] simp only [Complex.ofReal_div, Complex.ofReal_cos, Complex.ofReal_mul, Complex.ofReal_ofNat, Complex.ofReal_pow, Complex.ofReal_sin] ring have Φ1_inner_Φ2_re : (⟪H.Φ1, H.Φ2⟫_ℂ).re = cos t/(4 * t * (sin t)^2) * (1 - t * sin t) := by rw [Φ1_inner_Φ2, Complex.add_re, Complex.ofReal_re, Complex.re_mul_ofReal] simp have Φ1_inner_Φ2_im : (⟪H.Φ1, H.Φ2⟫_ℂ).im = cos t/(4 * t * (sin t)^2) * (- t * cos t) := by rw [Φ1_inner_Φ2, Complex.add_im, Complex.im_mul_ofReal, Complex.ofReal_im] simp have potential_H_cos_sin : potential .stabilityCounterExample H = - (cos t) ^ 2/ (4 * (sin t)^2) := by rw [potential, massTerm_stabilityCounterExample, quarticTerm_stabilityCounterExample] rw [Φ1_norm_sq, Φ2_norm_sq, Φ1_inner_Φ2_re, Φ1_inner_Φ2_im] field have potential_eq_c : potential .stabilityCounterExample H = - (|c| + 1) := by have htan : sin t / cos t = (2 * √(|c| + 1))⁻¹ := by rw [← tan_eq_sin_div_cos] exact tan_arctan _ rw [potential_H_cos_sin, show -cos t ^ 2 / (4 * sin t ^ 2) = -1 / (4 * (sin t / cos t) ^ 2) by field, htan] field_simp rw [sq_sqrt (by positivity)] ring /- Proving potential is unbounded. -/ use H rw [potential_eq_c] grind- Project
- Physlib
- License
- Apache-2.0
- Commit
- dd43e9e65791
- Source
- Physlib/Particles/BeyondTheStandardModel/TwoHDM/Potential.lean:446-527
Reuse this declaration
Bring the exact result into your workflow
The import identifies the source module. Your project still needs the pinned package dependency shown on this page.
What this badge means
This completion status comes from the project or community source. It has not yet been represented here as an independent rebuild and axiom audit.
Continue in this project
Related declarations
Adiabatic relation log
adiabatic_relation_log
Plain-language statement
Adiabatic relation in logarithmic form: If S(Ua,Va,N) = S(Ub,Vb,N) with N fixed, then c * log (Ua/Ub) + log (Va/Vb) = 0.
Source project: Physlib
Person-level attribution pending.
Adiabatic relation Ua Ub Va Vb
adiabatic_relation_UaUbVaVb
Plain-language statement
Adiabatic relation in product form: If S(Ua,Va,N) = S(Ub,Vb,N) with N fixed, then (Ua/Ub)^c * (Va/Vb) = 1.
Source project: Physlib
Person-level attribution pending.
Deriv Within mean Energy Beta eq neg variance
CanonicalEnsemble.derivWithin_meanEnergy_Beta_eq_neg_variance
Plain-language statement
(∂U/∂β) = -Var(E) for finite systems.
Source project: Physlib
Person-level attribution pending.