feat: Add Basic

This commit is contained in:
jstoobysmith 2024-04-17 15:12:20 -04:00
parent e9600cc68b
commit 28cc0ff712
2 changed files with 553 additions and 0 deletions

View file

@ -0,0 +1,552 @@
/-
Copyright (c) 2024 Joseph Tooby-Smith. All rights reserved.
Released under Apache 2.0 license.
Authors: Joseph Tooby-Smith
-/
import Mathlib.Tactic.FinCases
import Mathlib.Algebra.Module.Basic
import Mathlib.Tactic.Ring
import Mathlib.Algebra.GroupWithZero.Units.Lemmas
import HepLean.AnomalyCancellation.Basic
import Mathlib.Algebra.BigOperators.Fin
import Mathlib.Logic.Equiv.Fin
/-!
# The MSSM with 3 families and RHNs
We define the system of ACCs for the MSSM with 3 families and RHNs.
We define the system of charges for 1-species. We prove some basic lemmas about them.
-/
universe v u
open Nat
open BigOperators
/-- The vector space of charges corresponding to the MSSM fermions. -/
@[simps!]
def MSSMCharges : ACCSystemCharges := ACCSystemChargesMk 20
/-- THe vector spaces of charges of one species of fermions in the MSSM. -/
@[simps!]
def MSSMSpecies : ACCSystemCharges := ACCSystemChargesMk 3
namespace MSSMCharges
/-- An equivalence between `MSSMCharges.charges` and the space of maps
`(Fin 18 ⊕ Fin 2 → )`. The first 18 factors corresponds to the SM fermions, whils the last two
are the higgsions. -/
@[simps!]
def toSMPlusH : MSSMCharges.charges ≃ (Fin 18 ⊕ Fin 2 → ) :=
((@finSumFinEquiv 18 2).arrowCongr (Equiv.refl )).symm
/-- An equivalence between `Fin 18 ⊕ Fin 2 → ` and `(Fin 18 → ) × (Fin 2 → )`. -/
@[simps!]
def splitSMPlusH : (Fin 18 ⊕ Fin 2 → ) ≃ (Fin 18 → ) × (Fin 2 → ) where
toFun f := (f ∘ Sum.inl , f ∘ Sum.inr)
invFun f := Sum.elim f.1 f.2
left_inv f := by
aesop
right_inv f := by
aesop
/-- An equivalence between `MSSMCharges.charges` and `(Fin 18 → ) × (Fin 2 → )`. This
splits the charges up into the SM and the additional ones for the MSSM. -/
@[simps!]
def toSplitSMPlusH : MSSMCharges.charges ≃ (Fin 18 → ) × (Fin 2 → ) :=
toSMPlusH.trans splitSMPlusH
/-- An equivalence between `(Fin 18 → )` and `(Fin 6 → Fin 3 → )`. -/
@[simps!]
def toSpeciesMaps' : (Fin 18 → ) ≃ (Fin 6 → Fin 3 → ) :=
((Equiv.curry _ _ _).symm.trans
((@finProdFinEquiv 6 3).arrowCongr (Equiv.refl ))).symm
/-- An equivalence between `MSSMCharges.charges` and `(Fin 6 → Fin 3 → ) × (Fin 2 → ))`.
This split charges up into the SM and additional fermions, and further splits the SM into
species. -/
@[simps!]
def toSpecies : MSSMCharges.charges ≃ (Fin 6 → Fin 3 → ) × (Fin 2 → ) :=
toSplitSMPlusH.trans (Equiv.prodCongr toSpeciesMaps' (Equiv.refl _))
/-- For a given `i ∈ Fin 6` the projection of `MSSMCharges.charges` down to the
corresponding SM species of charges. -/
@[simps!]
def toSMSpecies (i : Fin 6) : MSSMCharges.charges →ₗ[] MSSMSpecies.charges where
toFun S := (Prod.fst ∘ toSpecies) S i
map_add' _ _ := by aesop
map_smul' _ _ := by aesop
lemma toSMSpecies_toSpecies_inv (i : Fin 6) (f : (Fin 6 → Fin 3 → ) × (Fin 2 → )) :
(toSMSpecies i) (toSpecies.symm f) = f.1 i := by
change (Prod.fst ∘ toSpecies ∘ toSpecies.symm ) _ i= f.1 i
simp
/-- The `Q` charges as a map `Fin 3 → `. -/
abbrev Q := toSMSpecies 0
/-- The `U` charges as a map `Fin 3 → `. -/
abbrev U := toSMSpecies 1
/-- The `D` charges as a map `Fin 3 → `. -/
abbrev D := toSMSpecies 2
/-- The `L` charges as a map `Fin 3 → `. -/
abbrev L := toSMSpecies 3
/-- The `E` charges as a map `Fin 3 → `. -/
abbrev E := toSMSpecies 4
/-- The `N` charges as a map `Fin 3 → `. -/
abbrev N := toSMSpecies 5
/-- The charge `Hd`. -/
@[simps!]
def Hd : MSSMCharges.charges →ₗ[] where
toFun S := S ⟨18, by simp⟩
map_add' _ _ := by aesop
map_smul' _ _ := by aesop
/-- The charge `Hu`. -/
@[simps!]
def Hu : MSSMCharges.charges →ₗ[] where
toFun S := S ⟨19, by simp⟩
map_add' _ _ := by aesop
map_smul' _ _ := by aesop
lemma charges_eq_toSpecies_eq (S T : MSSMCharges.charges) :
S = T ↔ (∀ i, toSMSpecies i S = toSMSpecies i T) ∧ Hd S = Hd T ∧ Hu S = Hu T := by
apply Iff.intro
intro h
rw [h]
simp only [forall_const, Hd_apply, Fin.reduceFinMk, Fin.isValue, Hu_apply, and_self]
intro h
apply toSpecies.injective
apply Prod.ext
funext i
exact h.1 i
funext i
match i with
| 0 => exact h.2.1
| 1 => exact h.2.2
lemma Hd_toSpecies_inv (f : (Fin 6 → Fin 3 → ) × (Fin 2 → )) :
Hd (toSpecies.symm f) = f.2 0 := by
rfl
lemma Hu_toSpecies_inv (f : (Fin 6 → Fin 3 → ) × (Fin 2 → )) :
Hu (toSpecies.symm f) = f.2 1 := by
rfl
end MSSMCharges
namespace MSSMACCs
open MSSMCharges
/-- The gravitational anomaly equation. -/
@[simp]
def accGrav : MSSMCharges.charges →ₗ[] where
toFun S := ∑ i, (6 * Q S i + 3 * U S i + 3 * D S i
+ 2 * L S i + E S i + N S i) + 2 * (Hd S + Hu S)
map_add' S T := by
simp only
repeat rw [map_add]
simp [mul_add]
repeat erw [Finset.sum_add_distrib]
ring
map_smul' a S := by
simp only
repeat rw [(toSMSpecies _).map_smul]
erw [Hd.map_smul, Hu.map_smul]
simp [HSMul.hSMul, SMul.smul]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
--rw [show Rat.cast a = a from rfl]
ring
/-- Extensionality lemma for `accGrav`. -/
lemma accGrav_ext {S T : MSSMCharges.charges}
(hj : ∀ (j : Fin 6), ∑ i, (toSMSpecies j) S i = ∑ i, (toSMSpecies j) T i)
(hd : Hd S = Hd T) (hu : Hu S = Hu T) :
accGrav S = accGrav T := by
simp only [accGrav, MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue,
Fin.reduceFinMk, LinearMap.coe_mk, AddHom.coe_mk]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
repeat erw [hj]
rw [hd, hu]
rfl
/-- The anomaly cancelation condition for SU(2) anomaly. -/
@[simp]
def accSU2 : MSSMCharges.charges →ₗ[] where
toFun S := ∑ i, (3 * Q S i + L S i) + Hd S + Hu S
map_add' S T := by
simp only
repeat rw [map_add]
simp [mul_add]
repeat erw [Finset.sum_add_distrib]
ring
map_smul' a S := by
simp only
repeat rw [(toSMSpecies _).map_smul]
erw [Hd.map_smul, Hu.map_smul]
simp [HSMul.hSMul, SMul.smul]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
--rw [show Rat.cast a = a from rfl]
ring
/-- Extensionality lemma for `accSU2`. -/
lemma accSU2_ext {S T : MSSMCharges.charges}
(hj : ∀ (j : Fin 6), ∑ i, (toSMSpecies j) S i = ∑ i, (toSMSpecies j) T i)
(hd : Hd S = Hd T) (hu : Hu S = Hu T) :
accSU2 S = accSU2 T := by
simp only [accSU2, MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue,
Fin.reduceFinMk, LinearMap.coe_mk, AddHom.coe_mk]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
repeat erw [hj]
rw [hd, hu]
rfl
/-- The anomaly cancelation condition for SU(3) anomaly. -/
@[simp]
def accSU3 : MSSMCharges.charges →ₗ[] where
toFun S := ∑ i, (2 * (Q S i) + (U S i) + (D S i))
map_add' S T := by
simp only
repeat rw [map_add]
simp [mul_add]
repeat erw [Finset.sum_add_distrib]
ring
map_smul' a S := by
simp only
repeat rw [(toSMSpecies _).map_smul]
simp [HSMul.hSMul, SMul.smul]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
--rw [show Rat.cast a = a from rfl]
ring
/-- Extensionality lemma for `accSU3`. -/
lemma accSU3_ext {S T : MSSMCharges.charges}
(hj : ∀ (j : Fin 6), ∑ i, (toSMSpecies j) S i = ∑ i, (toSMSpecies j) T i) :
accSU3 S = accSU3 T := by
simp only [accSU3, MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue, LinearMap.coe_mk,
AddHom.coe_mk]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
repeat erw [hj]
rfl
/-- The acc for `Y²`. -/
@[simp]
def accYY : MSSMCharges.charges →ₗ[] where
toFun S := ∑ i, ((Q S) i + 8 * (U S) i + 2 * (D S) i + 3 * (L S) i
+ 6 * (E S) i) + 3 * (Hd S + Hu S)
map_add' S T := by
simp only
repeat rw [map_add]
simp [mul_add]
repeat erw [Finset.sum_add_distrib]
ring
map_smul' a S := by
simp only
repeat rw [(toSMSpecies _).map_smul]
erw [Hd.map_smul, Hu.map_smul]
simp [HSMul.hSMul, SMul.smul]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
-- rw [show Rat.cast a = a from rfl]
ring
/-- Extensionality lemma for `accGrav`. -/
lemma accYY_ext {S T : MSSMCharges.charges}
(hj : ∀ (j : Fin 6), ∑ i, (toSMSpecies j) S i = ∑ i, (toSMSpecies j) T i)
(hd : Hd S = Hd T) (hu : Hu S = Hu T) :
accYY S = accYY T := by
simp only [accYY, MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue,
Fin.reduceFinMk, LinearMap.coe_mk, AddHom.coe_mk]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
repeat erw [hj]
rw [hd, hu]
rfl
/-- The symmetric bilinear form used to define the quadratic ACC. -/
@[simps!]
def quadBiLin : BiLinearSymm MSSMCharges.charges where
toFun S := ∑ i, (Q S.1 i * Q S.2 i + (- 2) * (U S.1 i * U S.2 i) +
D S.1 i * D S.2 i + (- 1) * (L S.1 i * L S.2 i) + E S.1 i * E S.2 i) +
(- Hd S.1 * Hd S.2 + Hu S.1 * Hu S.2)
map_smul₁' a S T := by
simp only
rw [mul_add]
congr 1
rw [Finset.mul_sum]
apply Fintype.sum_congr
intro i
repeat erw [map_smul]
simp only [HSMul.hSMul, SMul.smul, toSMSpecies_apply, Fin.isValue, neg_mul, one_mul]
ring
simp only [map_smul, Hd_apply, Fin.reduceFinMk, Fin.isValue, smul_eq_mul, neg_mul, Hu_apply]
ring
map_add₁' S T R := by
simp only
rw [add_assoc, ← add_assoc (-Hd S * Hd R + Hu S * Hu R) _ _]
rw [add_comm (-Hd S * Hd R + Hu S * Hu R) _]
rw [add_assoc]
rw [← add_assoc _ _ (-Hd S * Hd R + Hu S * Hu R + (-Hd T * Hd R + Hu T * Hu R))]
congr 1
rw [← Finset.sum_add_distrib]
apply Fintype.sum_congr
intro i
repeat erw [map_add]
simp only [ACCSystemCharges.chargesAddCommMonoid_add, toSMSpecies_apply, Fin.isValue, neg_mul,
one_mul]
ring
rw [Hd.map_add, Hu.map_add]
ring
swap' S L := by
simp only [MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue, neg_mul, one_mul,
Hd_apply, Fin.reduceFinMk, Hu_apply]
congr 1
rw [Fin.sum_univ_three, Fin.sum_univ_three]
simp only [Fin.isValue]
ring
ring
/-- The quadratic ACC. -/
@[simp]
def accQuad : HomogeneousQuadratic MSSMCharges.charges := quadBiLin.toHomogeneousQuad
/-- Extensionality lemma for `accQuad`. -/
lemma accQuad_ext {S T : (MSSMCharges).charges}
(h : ∀ j, ∑ i, ((fun a => a^2) ∘ toSMSpecies j S) i =
∑ i, ((fun a => a^2) ∘ toSMSpecies j T) i)
(hd : Hd S = Hd T) (hu : Hu S = Hu T) :
accQuad S = accQuad T := by
simp only [accQuad, BiLinearSymm.toHomogeneousQuad_toFun]
erw [← quadBiLin.toFun_eq_coe]
rw [quadBiLin]
simp only
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
ring_nf
have h1 : ∀ j, ∑ i, (toSMSpecies j S i)^2 = ∑ i, (toSMSpecies j T i)^2 := by
intro j
erw [h]
rfl
repeat rw [h1]
rw [hd, hu]
/-- The function underlying the symmetric trilinear form used to define the cubic ACC. -/
@[simp]
def cubeTriLinToFun
(S : MSSMCharges.charges × MSSMCharges.charges × MSSMCharges.charges) : :=
∑ i, (6 * (Q S.1 i * Q S.2.1 i * Q S.2.2 i)
+ 3 * (U S.1 i * U S.2.1 i * U S.2.2 i)
+ 3 * (D S.1 i * D S.2.1 i * D S.2.2 i)
+ 2 * (L S.1 i * L S.2.1 i * L S.2.2 i)
+ E S.1 i * E S.2.1 i * E S.2.2 i
+ N S.1 i * N S.2.1 i * N S.2.2 i)
+ (2 * Hd S.1 * Hd S.2.1 * Hd S.2.2
+ 2 * Hu S.1 * Hu S.2.1 * Hu S.2.2)
lemma cubeTriLinToFun_map_smul₁ (a : ) (S T R : MSSMCharges.charges) :
cubeTriLinToFun (a • S, T, R) = a * cubeTriLinToFun (S, T, R) := by
simp only [cubeTriLinToFun]
rw [mul_add]
congr 1
rw [Finset.mul_sum]
apply Fintype.sum_congr
intro i
repeat erw [map_smul]
simp only [HSMul.hSMul, SMul.smul, toSMSpecies_apply, Fin.isValue]
ring
simp only [map_smul, Hd_apply, Fin.reduceFinMk, Fin.isValue, smul_eq_mul, Hu_apply]
ring
lemma cubeTriLinToFun_map_add₁ (S T R L : MSSMCharges.charges) :
cubeTriLinToFun (S + T, R, L) = cubeTriLinToFun (S, R, L) + cubeTriLinToFun (T, R, L) := by
simp only [cubeTriLinToFun]
rw [add_assoc, ← add_assoc (2 * Hd S * Hd R * Hd L + 2 * Hu S * Hu R * Hu L) _ _]
rw [add_comm (2 * Hd S * Hd R * Hd L + 2 * Hu S * Hu R * Hu L) _]
rw [add_assoc]
rw [← add_assoc _ _ (2 * Hd S * Hd R * Hd L + 2 * Hu S * Hu R * Hu L +
(2 * Hd T * Hd R * Hd L + 2 * Hu T * Hu R * Hu L))]
congr 1
rw [← Finset.sum_add_distrib]
apply Fintype.sum_congr
intro i
repeat erw [map_add]
simp only [ACCSystemCharges.chargesAddCommMonoid_add, toSMSpecies_apply, Fin.isValue]
ring
rw [Hd.map_add, Hu.map_add]
ring
lemma cubeTriLinToFun_swap1 (S T R : MSSMCharges.charges) :
cubeTriLinToFun (S, T, R) = cubeTriLinToFun (T, S, R) := by
simp only [cubeTriLinToFun, MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue, Hd_apply,
Fin.reduceFinMk, Hu_apply]
congr 1
rw [Fin.sum_univ_three, Fin.sum_univ_three]
simp only [Fin.isValue]
ring
ring
lemma cubeTriLinToFun_swap2 (S T R : MSSMCharges.charges) :
cubeTriLinToFun (S, T, R) = cubeTriLinToFun (S, R, T) := by
simp only [cubeTriLinToFun, MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue, Hd_apply,
Fin.reduceFinMk, Hu_apply]
congr 1
rw [Fin.sum_univ_three, Fin.sum_univ_three]
simp only [Fin.isValue]
ring
ring
/-- The symmetric trilinear form used to define the cubic ACC. -/
def cubeTriLin : TriLinearSymm MSSMCharges.charges where
toFun S := cubeTriLinToFun S
map_smul₁' := cubeTriLinToFun_map_smul₁
map_add₁' := cubeTriLinToFun_map_add₁
swap₁' := cubeTriLinToFun_swap1
swap₂' := cubeTriLinToFun_swap2
/-- The cubic ACC. -/
@[simp]
def accCube : HomogeneousCubic MSSMCharges.charges := cubeTriLin.toCubic
/-- Extensionality lemma for `accCube`. -/
lemma accCube_ext {S T : MSSMCharges.charges}
(h : ∀ j, ∑ i, ((fun a => a^3) ∘ toSMSpecies j S) i =
∑ i, ((fun a => a^3) ∘ toSMSpecies j T) i)
(hd : Hd S = Hd T) (hu : Hu S = Hu T) :
accCube S = accCube T := by
simp [cubeTriLin, cubeTriLinToFun]
erw [← cubeTriLin.toFun_eq_coe]
rw [cubeTriLin]
simp only [cubeTriLinToFun]
repeat erw [Finset.sum_add_distrib]
repeat erw [← Finset.mul_sum]
ring_nf
have h1 : ∀ j, ∑ i, (toSMSpecies j S i)^3 = ∑ i, (toSMSpecies j T i)^3 := by
intro j
erw [h]
rfl
repeat rw [h1]
rw [hd, hu]
end MSSMACCs
open MSSMACCs
/-- The ACCSystem for the MSSM without RHN. -/
@[simps!]
def MSSMACC : ACCSystem where
numberLinear := 4
linearACCs := fun i =>
match i with
| 0 => accGrav
| 1 => accSU2
| 2 => accSU3
| 3 => accYY
numberQuadratic := 1
quadraticACCs := fun i =>
match i with
| 0 => accQuad
cubicACC := accCube
namespace MSSMACC
open MSSMCharges
lemma quadSol (S : MSSMACC.QuadSols) : accQuad S.val = 0 := by
have hS := S.quadSol
simp [MSSMACCs.accQuad, HomogeneousQuadratic.toFun] at hS
exact hS 0
/-- A solution from a charge satisfying the ACCs. -/
@[simp]
def AnomalyFreeMk (S : MSSMACC.charges) (hg : accGrav S = 0)
(hsu2 : accSU2 S = 0) (hsu3 : accSU3 S = 0) (hyy : accYY S = 0)
(hquad : accQuad S = 0) (hcube : accCube S = 0) : MSSMACC.Sols :=
⟨⟨⟨S, by
intro i
simp at i
match i with
| 0 => exact hg
| 1 => exact hsu2
| 2 => exact hsu3
| 3 => exact hyy⟩, by
intro i
simp at i
match i with
| 0 => exact hquad
⟩ , by exact hcube ⟩
lemma AnomalyFreeMk_val (S : MSSMACC.charges) (hg : accGrav S = 0)
(hsu2 : accSU2 S = 0) (hsu3 : accSU3 S = 0) (hyy : accYY S = 0)
(hquad : accQuad S = 0) (hcube : accCube S = 0) :
(AnomalyFreeMk S hg hsu2 hsu3 hyy hquad hcube).val = S := by
rfl
/-- A `QuadSol` from a `LinSol` satisfying the quadratic ACC. -/
@[simp]
def AnomalyFreeQuadMk' (S : MSSMACC.LinSols) (hquad : accQuad S.val = 0) :
MSSMACC.QuadSols :=
⟨S, by
intro i
simp at i
match i with
| 0 => exact hquad
/-- A `Sol` from a `LinSol` satisfying the quadratic and cubic ACCs. -/
@[simp]
def AnomalyFreeMk' (S : MSSMACC.LinSols) (hquad : accQuad S.val = 0)
(hcube : accCube S.val = 0) : MSSMACC.Sols :=
⟨⟨S, by
intro i
simp at i
match i with
| 0 => exact hquad
⟩ , by exact hcube ⟩
/-- A `Sol` from a `QuadSol` satisfying the cubic ACCs. -/
@[simp]
def AnomalyFreeMk'' (S : MSSMACC.QuadSols) (hcube : accCube S.val = 0) : MSSMACC.Sols :=
⟨S , by exact hcube ⟩
lemma AnomalyFreeMk''_val (S : MSSMACC.QuadSols)
(hcube : accCube S.val = 0) :
(AnomalyFreeMk'' S hcube).val = S.val := by
rfl
/-- The dot product on the vector space of charges. -/
@[simps!]
def dot : BiLinearSymm MSSMCharges.charges where
toFun S := ∑ i, (Q S.1 i * Q S.2 i + U S.1 i * U S.2 i +
D S.1 i * D S.2 i + L S.1 i * L S.2 i + E S.1 i * E S.2 i
+ N S.1 i * N S.2 i) + Hd S.1 * Hd S.2 + Hu S.1 * Hu S.2
map_smul₁' a S T := by
simp only [MSSMSpecies_numberCharges]
repeat rw [(toSMSpecies _).map_smul]
rw [Hd.map_smul, Hu.map_smul]
rw [Fin.sum_univ_three, Fin.sum_univ_three]
simp only [HSMul.hSMul, SMul.smul, Fin.isValue, toSMSpecies_apply, Hd_apply, Fin.reduceFinMk,
Hu_apply]
ring
map_add₁' S T R := by
simp only [MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue,
ACCSystemCharges.chargesAddCommMonoid_add, map_add, Hd_apply, Fin.reduceFinMk, Hu_apply]
repeat erw [AddHom.map_add]
rw [Fin.sum_univ_three, Fin.sum_univ_three, Fin.sum_univ_three]
simp only [Fin.isValue]
ring
swap' S L := by
simp only [MSSMSpecies_numberCharges, toSMSpecies_apply, Fin.isValue, Hd_apply, Fin.reduceFinMk,
Hu_apply]
rw [Fin.sum_univ_three, Fin.sum_univ_three]
simp only [Fin.isValue]
ring
end MSSMACC