SM + CxSM + other fixes

docs/index.rst

@@ -54,8 +54,10 @@ BSMPT - Beyond the Standard Model Phase Transitions documentation
 
    models/c2hdm
    models/cp_in_the_dark
-   models/r2hdm
+   models/cxsm
    models/n2hdm
+   models/sm
+   models/r2hdm
 
 .. toctree::
    :maxdepth: 1

docs/models/c2hdm.rst

@@ -22,7 +22,7 @@ The potential is given by
 .. math::
    \begin{aligned}V_{\text {C2HDM}}= & m_{11}^2 \Phi_1^{\dagger} \Phi_1+m_{22}^2 \Phi_2^{\dagger} \Phi_2-\left[m_{12}^2 \Phi_1^{\dagger} \Phi_2+\text { h.c. }\right]+\frac{1}{2} \lambda_1\left(\Phi_1^{\dagger} \Phi_1\right)^2+\frac{1}{2} \lambda_2\left(\Phi_2^{\dagger} \Phi_2\right)^2 \\& +\lambda_3\left(\Phi_1^{\dagger} \Phi_1\right)\left(\Phi_2^{\dagger} \Phi_2\right)+\lambda_4\left(\Phi_1^{\dagger} \Phi_2\right)\left(\Phi_2^{\dagger} \Phi_1\right)+\left[\frac{1}{2} \lambda_5\left(\Phi_1^{\dagger} \Phi_2\right)^2+\text { h.c. }\right] .\end{aligned}
 
-where all coupling constants are real except :math:`m_{12}` and :math:`\lambda_5` which, if non-real, explicitly break CP in the model at :math:`T = 0`.
+where all parameters are real except :math:`m_{12}` and :math:`\lambda_5` which, if non-real, explicitly break CP in the model at :math:`T = 0`.
 
 The :math:`T=0` vacuum is chosen to be
 

docs/models/cp_in_the_dark.rst

@@ -5,7 +5,7 @@ CP in the dark
 
 This implementation was based on [`1807.10322 <https://arxiv.org/abs/1807.10322>`_] and [`2204.13425 <https://arxiv.org/abs/2204.13425>`_]. 
 
-On the scalar sector there are two doublets 
+On the scalar sector there are two doublets and a real scalar field
 
 .. math::
    \Phi_1=\frac{1}{\sqrt{2}}\binom{\rho_1+\mathrm{i} \eta_1}{\zeta_1+\omega_1+\mathrm{i} \psi_1}, \quad \Phi_2=\frac{1}{\sqrt{2}}\binom{\rho_2+\omega_{\mathrm{CB}}+\mathrm{i} \eta_2}{\zeta_2+\omega_2+\mathrm{i}\left(\psi_2+\omega_{\mathrm{CP}}\right)}, \quad \Phi_S=\zeta_3+\omega_S
@@ -24,7 +24,7 @@ The potential is given by
 .. math::
     \begin{aligned}V_\text{CP in the Dark} &= m_{11}^2\left|\Phi_1\right|^2+m_{22}^2\left|\Phi_2\right|^2+\frac{1}{2} m_S^2 \Phi_S^2+\left(A \Phi_1^{\dagger} \Phi_2 \Phi_S+\text { h.c. }\right) \\& +\frac{1}{2} \lambda_1\left|\Phi_1\right|^4+\frac{1}{2} \lambda_2\left|\Phi_2\right|^4+\lambda_3\left|\Phi_1\right|^2\left|\Phi_2\right|^2+\lambda_4\left|\Phi_1^{\dagger} \Phi_2\right|^2+\frac{1}{2} \lambda_5\left[\left(\Phi_1^{\dagger} \Phi_2\right)^2+\text { h.c. }\right] \\& +\frac{1}{4} \lambda_6 \Phi_S^4+\frac{1}{2} \lambda_7\left|\Phi_1\right|^2 \Phi_S^2+\frac{1}{2} \lambda_8\left|\Phi_2\right|^2 \Phi_S^2,\end{aligned}
 
-where all coupling constants are real except :math:`A`, whose complex phase is directly linked with CP violation.
+where all parameters are real except :math:`A`, whose complex phase is directly linked with CP violation.
 
 The :math:`T=0` vacuum is chosen to be
 

docs/models/cxsm.rst

@@ -0,0 +1,24 @@
+.. _cxsm:
+
+CxSM
+==============
+
+This implementation was based on [`1512.05355 <https://arxiv.org/abs/1512.05355>`_]. 
+
+On the scalar sector there are a doublet and a complex scalar field
+
+.. math::
+   H=\frac{1}{\sqrt{2}}\binom{G^{+}}{v+h+i G^0} \quad \text { and } \quad \mathbb{S}=\frac{1}{\sqrt{2}}\left[v_S+s+i\left(v_A+a\right)\right]
+
+The potential is given by
+
+.. math::
+   V_{\mathrm{CxSM}}=\frac{m^2}{2} H^{\dagger} H+\frac{\lambda}{4}\left(H^{\dagger} H\right)^2+\frac{\delta_2}{2} H^{\dagger} H|\mathbb{S}|^2+\frac{b_2}{2}|\mathbb{S}|^2+\frac{d_2}{4}|\mathbb{S}|^4+\left(\frac{b_1}{4} \mathbb{S}^2+a_1 \mathbb{S}+c . c .\right)
+
+where all parameters are real except :math:`a_1` and :math:`b_1` (one of the phases can be absorbed by a phase rotation).
+
+The :math:`T=0` vacuum is chosen to be
+
+.. math::
+   \left.\left\langle\Phi_1\right\rangle\right|_{T=0}=\frac{1}{\sqrt{2}}\binom{0}{v_1},\left.\quad\left\langle\mathbb{S}\right\rangle\right|_{T=0}=\frac{1}{\sqrt{2}}(v_s+i v_a)\,.
+

docs/models/n2hdm.rst

@@ -5,7 +5,7 @@ N2HDM
 
 This implementation was based on [`1612.01309 <https://arxiv.org/abs/1612.01309>`_] and [`1912.10477 <https://arxiv.org/abs/1912.10477>`_]. 
 
-On the scalar sector there are two doublets 
+On the scalar sector there are two doublets and a real scalar field
 
 .. math::
    \Phi_1=\frac{1}{\sqrt{2}}\binom{\rho_1+\mathrm{i} \eta_1}{\zeta_1+\omega_1+\mathrm{i} \psi_1}, \quad \Phi_2=\frac{1}{\sqrt{2}}\binom{\rho_2+\omega_{\mathrm{CB}}+\mathrm{i} \eta_2}{\zeta_2+\omega_2+\mathrm{i}\left(\psi_2+\omega_{\mathrm{CP}}\right)}, \quad \Phi_S=\zeta_3+\omega_S
@@ -33,7 +33,7 @@ The potential is given by
 .. math::
    \begin{aligned}V_{\mathrm{N} 2 \mathrm{HDM}}= & m_{11}^2 \Phi_1^{\dagger} \Phi_1+m_{22}^2 \Phi_2^{\dagger} \Phi_2-m_{12}^2\left(\Phi_1^{\dagger} \Phi_2+\text { h.c. }\right)+\frac{\lambda_1}{2}\left(\Phi_1^{\dagger} \Phi_1\right)^2+\frac{\lambda_2}{2}\left(\Phi_2^{\dagger} \Phi_2\right)^2 \\& +\lambda_3 \Phi_1^{\dagger} \Phi_1 \Phi_2^{\dagger} \Phi_2+\lambda_4 \Phi_1^{\dagger} \Phi_2 \Phi_2^{\dagger} \Phi_1+\frac{\lambda_5}{2}\left(\left(\Phi_1^{\dagger} \Phi_2\right)^2+\text { h.c. }\right) \\& +\frac{1}{2} m_S^2 \Phi_S^2+\frac{\lambda_6}{8} \Phi_S^4+\lambda_7\left(\Phi_1^{\dagger} \Phi_1\right) \Phi_S^2+\lambda_8\left(\Phi_2^{\dagger} \Phi_2\right) \Phi_S^2,\end{aligned}
 
-where all coupling constants are real to preserve CP in the model at :math:`T = 0`.
+where all parameters are real to preserve CP in the model at :math:`T = 0`.
 
 The :math:`T=0` vacuum is chosen to be
 

docs/models/r2hdm.rst

@@ -23,7 +23,7 @@ The potential is given by
 .. math::
    \begin{aligned}V_{\text {R2HDM}}= & m_{11}^2 \Phi_1^{\dagger} \Phi_1+m_{22}^2 \Phi_2^{\dagger} \Phi_2-\left[m_{12}^2 \Phi_1^{\dagger} \Phi_2+\text { h.c. }\right]+\frac{1}{2} \lambda_1\left(\Phi_1^{\dagger} \Phi_1\right)^2+\frac{1}{2} \lambda_2\left(\Phi_2^{\dagger} \Phi_2\right)^2 \\& +\lambda_3\left(\Phi_1^{\dagger} \Phi_1\right)\left(\Phi_2^{\dagger} \Phi_2\right)+\lambda_4\left(\Phi_1^{\dagger} \Phi_2\right)\left(\Phi_2^{\dagger} \Phi_1\right)+\left[\frac{1}{2} \lambda_5\left(\Phi_1^{\dagger} \Phi_2\right)^2+\text { h.c. }\right] .\end{aligned}
 
-where all coupling constants are real. :math:`m_{12}` and :math:`\lambda_5` could, in principle, be complex but are set to real to that there is no explicit CP violation in the model at :math:`T = 0`.
+where all parameters are real. :math:`m_{12}` and :math:`\lambda_5` could, in principle, be complex but are set to real to that there is no explicit CP violation in the model at :math:`T = 0`.
 
 The :math:`T=0` vacuum is chosen to be
 

docs/models/sm.rst

@@ -0,0 +1,24 @@
+.. _SM:
+
+SM
+==============
+
+This implementation was based on [`1512.05355 <https://arxiv.org/abs/1512.05355>`_]. 
+
+On the scalar sector there is a doublet 
+
+.. math::
+   H=\frac{1}{\sqrt{2}}\binom{G^{+}}{v+h+i G^0}
+
+The potential is given by
+
+.. math::
+   V_{\mathrm{CxSM}}=\frac{m^2}{2} H^{\dagger} H+\frac{\lambda}{4}\left(H^{\dagger} H\right)^2
+
+where all parameters are real.
+
+The :math:`T=0` vacuum is chosen to be
+
+.. math::
+   \left.\left\langle\Phi_1\right\rangle\right|_{T=0}=\frac{1}{\sqrt{2}}\binom{0}{v_1}\,.
+