removed thermal mass from ManySinglets, Yukawa and SM

Author: jorindevandevis

Models/ManySinglets/manySinglets.py

@@ -123,15 +123,11 @@ def topMsqDerivative(fields: Fields) -> Fields:
                 ]
             )
 
-        def topMsqThermal(T: float) -> float:
-            return self.modelParameters["g3"] ** 2 * T**2 / 6.0
-
         topQuark = Particle(
             "top",
             index=0,
             msqVacuum=topMsqVacuum,
             msqDerivative=topMsqDerivative,
-            msqThermal=topMsqThermal,
             statistics="Fermion",
             totalDOFs=12,
         )

Models/StandardModel/standardModel.py

@@ -119,15 +119,11 @@ def topMsqVacuum(fields: Fields) -> Fields:
         def topMsqDerivative(fields: Fields) -> Fields:
             return self.modelParameters["yt"] ** 2 * fields.getField(0)
 
-        def topMsqThermal(T: float) -> float:
-            return self.modelParameters["g3"] ** 2 * T**2 / 6.0
-
         topQuarkL = Particle(
             name="TopL",
             index=0,
             msqVacuum=topMsqVacuum,
             msqDerivative=topMsqDerivative,
-            msqThermal=topMsqThermal,
             statistics="Fermion",
             totalDOFs=6,
         )
@@ -138,7 +134,6 @@ def topMsqThermal(T: float) -> float:
             index=1,
             msqVacuum=topMsqVacuum,
             msqDerivative=topMsqDerivative,
-            msqThermal=topMsqThermal,
             statistics="Fermion",
             totalDOFs=6,
         )
@@ -151,15 +146,11 @@ def WMsqVacuum(fields: Fields) -> Fields:  # pylint: disable=invalid-name
         def WMsqDerivative(fields: Fields) -> Fields:  # pylint: disable=invalid-name
             return self.modelParameters["g2"] ** 2 * fields.getField(0) / 2
 
-        def WMsqThermal(T: float) -> float:  # pylint: disable=invalid-name
-            return self.modelParameters["g2"] ** 2 * T**2 * 3.0 / 5.0
-
         wBoson = Particle(
             name="W",
             index=4,
             msqVacuum=WMsqVacuum,
             msqDerivative=WMsqDerivative,
-            msqThermal=WMsqThermal,
             statistics="Boson",
             totalDOFs=9,
         )

Models/Yukawa/yukawa.py

@@ -64,16 +64,11 @@ def psiMsqDerivative(fields: Fields) -> Fields:
                 )
             )
 
-        # Asymptotic thermal mass (twice the static thermal mass)
-        def psiMsqThermal(T: float) -> float:
-            return 1 / 8 * self.modelParameters["y"] ** 2 * T**2
-
         psiL = Particle(
             "psiL",
             index=1,
             msqVacuum=psiMsqVacuum,
             msqDerivative=psiMsqDerivative,
-            msqThermal=psiMsqThermal,
             statistics="Fermion",
             totalDOFs=2,
         )
@@ -82,7 +77,6 @@ def psiMsqThermal(T: float) -> float:
             index=2,
             msqVacuum=psiMsqVacuum,
             msqDerivative=psiMsqDerivative,
-            msqThermal=psiMsqThermal,
             statistics="Fermion",
             totalDOFs=2,
         )