Chapter 13: Beyond the Standard Model

Despite its extraordinary success, the Standard Model leaves fundamental questions unanswered. The hierarchy problem asks why the Higgs mass is so much lighter than the Planck scale. The strong CP problem demands an explanation for the vanishing neutron electric dipole moment. Neutrino masses, dark matter, and the failure of gauge coupling unification all point toward new physics beyond the Standard Model.

Derivation 1: The Hierarchy Problem

The Higgs boson mass receives quadratically divergent radiative corrections. The dominant contribution comes from the top quark loop:

$\delta m_H^2 = -\frac{3y_t^2}{8\pi^2}\Lambda^2 + \ldots$

where $\Lambda$ is the UV cutoff (the scale of new physics) and $y_t \approx 1$is the top Yukawa coupling. If the Standard Model is valid up to the Planck scale$M_\text{Pl} \approx 10^{19}$ GeV, then:

$\delta m_H^2 \sim -(10^{17} \text{ GeV})^2$

Yet the physical Higgs mass is $m_H = 125$ GeV, requiring a cancellation between the bare mass and radiative corrections to 34 decimal places. This extreme fine-tuning is the hierarchy problem.

Naturalness Criterion

't Hooft's naturalness criterion states that a small parameter is natural only if setting it to zero increases the symmetry of the theory. For the Higgs mass, there is no such symmetry in the Standard Model (unlike fermion masses, which are protected by chiral symmetry). Proposed solutions include:

Supersymmetry: Fermion and boson loops cancel due to a new fermion-boson symmetry, stabilizing $m_H$

Composite Higgs: The Higgs is a pseudo-Goldstone boson of a new strong interaction at the TeV scale

Extra dimensions: The fundamental Planck scale is actually at the TeV scale; gravity appears weak because it propagates in extra dimensions

Status: The LHC has not found evidence for any of these solutions up to scales of several TeV, creating a tension known as the "little hierarchy problem." Whether naturalness is the right guiding principle remains actively debated.

Derivation 2: The Strong CP Problem

QCD permits a CP-violating term in the Lagrangian:

$\mathcal{L}_\theta = \frac{\theta}{32\pi^2} G_{\mu\nu}^a \tilde{G}^{a\mu\nu}$

where $\tilde{G}^{\mu\nu} = \frac{1}{2}\epsilon^{\mu\nu\rho\sigma}G_{\rho\sigma}$ is the dual field strength. This term is a total derivative but has physical consequences due to the non-trivial topology of the QCD vacuum (instantons). The effective parameter is $\bar{\theta} = \theta + \arg\det(M_q)$ where $M_q$ is the quark mass matrix.

Neutron EDM Constraint

The $\theta$ term induces a neutron electric dipole moment:

$d_n \approx 3.6 \times 10^{-16}\,\bar{\theta} \text{ e}\cdot\text{cm}$

The experimental bound $|d_n| < 1.8 \times 10^{-26}$ e$\cdot$cm requires$|\bar{\theta}| < 10^{-10}$. Why this dimensionless parameter is so extraordinarily small is the strong CP problem. There is no symmetry reason within the SM for $\bar{\theta} = 0$.

Peccei-Quinn Solution

Peccei and Quinn proposed a new global $U(1)_\text{PQ}$ symmetry that is spontaneously broken at a high scale $f_a$. The associated pseudo-Goldstone boson — the axion — dynamically relaxes $\bar{\theta} \to 0$:

$\bar{\theta}_\text{eff} = \bar{\theta} + \frac{a}{f_a} \to 0 \text{ (minimized dynamically)}$

Axion as dark matter: For $f_a \sim 10^{10}\text{-}10^{12}$ GeV, the axion is a viable cold dark matter candidate with mass $m_a \sim 10^{-5}\text{-}10^{-3}$ eV. Experiments like ADMX are actively searching for axion dark matter through its coupling to photons in a strong magnetic field.

Derivation 3: Neutrino Masses and the Seesaw Mechanism

Neutrino oscillation experiments have established that neutrinos have nonzero masses, but the Standard Model (with only left-handed neutrinos) predicts them to be exactly massless. The simplest extension adds right-handed neutrinos $N_R$ with both Dirac and Majorana mass terms:

$\mathcal{L}_\text{mass} = -m_D \bar{\nu}_L N_R - \frac{1}{2}M_R \overline{N_R^c} N_R + \text{h.c.}$

where $m_D = y_\nu v/\sqrt{2}$ is the Dirac mass and $M_R$ is the Majorana mass. The mass matrix in the $(\nu_L, N_R)$ basis is:

$M = \begin{pmatrix} 0 & m_D \\ m_D & M_R \end{pmatrix}$

For $M_R \gg m_D$, diagonalization gives the light neutrino mass:

$m_\nu \approx \frac{m_D^2}{M_R} = \frac{y_\nu^2 v^2}{2M_R}$

This is the Type-I seesaw: the light neutrino mass is suppressed by the large scale$M_R$. For $m_\nu \sim 0.05$ eV and $y_\nu \sim 1$:$M_R \sim 10^{14}$ GeV, tantalizingly close to the GUT scale.

Leptogenesis: The heavy right-handed neutrinos can generate a lepton asymmetry through their CP-violating decays. This lepton asymmetry is then converted to a baryon asymmetry by sphaleron processes, potentially explaining why the universe contains more matter than antimatter.

Derivation 4: Dark Matter and the WIMP Miracle

Cosmological and astrophysical observations require that approximately 27% of the energy density of the universe is non-baryonic dark matter. In the thermal relic scenario, dark matter particles were in equilibrium in the early universe and froze out when the annihilation rate fell below the Hubble expansion rate.

Relic Abundance

The relic density of a thermal relic is approximately:

$\Omega_\chi h^2 \approx \frac{3 \times 10^{-27} \text{ cm}^3/\text{s}}{\langle\sigma v\rangle}$

For a weakly interacting massive particle (WIMP) with mass $m_\chi \sim 100$ GeV and weak-scale coupling:

$\langle\sigma v\rangle \sim \frac{\alpha_W^2}{m_\chi^2} \sim \frac{(1/30)^2}{(100 \text{ GeV})^2} \sim 3 \times 10^{-26} \text{ cm}^3/\text{s}$

This gives $\Omega h^2 \sim 0.1$, remarkably close to the observed value$\Omega_\text{DM} h^2 = 0.120 \pm 0.001$. This numerical coincidence — that a particle with weak-scale mass and coupling naturally produces the observed dark matter abundance — is the WIMP miracle.

Detection strategies: Direct detection (nuclear recoils in underground detectors), indirect detection (annihilation products in cosmic rays), and collider production (missing energy signatures at the LHC). Despite decades of searches, no conclusive WIMP signal has been found, pushing the parameter space toward heavier masses or weaker couplings.

Derivation 5: Gauge Coupling Unification

The three Standard Model gauge couplings run with energy according to the one-loop renormalization group equations:

$\alpha_i^{-1}(Q) = \alpha_i^{-1}(M_Z) - \frac{b_i}{2\pi}\ln\frac{Q}{M_Z}$

where the one-loop beta function coefficients are:

SM: $b_1 = 41/10, \; b_2 = -19/6, \; b_3 = -7$

MSSM: $b_1 = 33/5, \; b_2 = 1, \; b_3 = -3$

In the Standard Model, the three couplings approach each other at high energies but do not meet at a single point. In the MSSM (with superpartners at ~1 TeV), the couplings unify precisely at $M_\text{GUT} \approx 2 \times 10^{16}$ GeV with $\alpha_\text{GUT}^{-1} \approx 24$.

Grand Unified Theories

The simplest GUT is $SU(5)$ (Georgi-Glashow, 1974), which embeds the SM gauge group:

$SU(3)_c \times SU(2)_L \times U(1)_Y \subset SU(5)$

GUTs predict proton decay ($p \to e^+\pi^0$) with lifetime$\tau_p \sim M_\text{GUT}^4 / (m_p^5 \alpha_\text{GUT}^2)$. The Super-Kamiokande bound $\tau_p > 2.4 \times 10^{34}$ years rules out minimal $SU(5)$ but is consistent with SUSY GUTs (where the dominant mode is $p \to K^+ \bar{\nu}$).

Proton Decay

In GUTs, quarks and leptons are unified in the same multiplet, and the heavy GUT gauge bosons $X, Y$ mediate baryon number violation. The dominant proton decay mode in minimal SU(5) is $p \to e^+\pi^0$ with rate:

$\tau_p \sim \frac{M_X^4}{\alpha_\text{GUT}^2 m_p^5} \sim 10^{31\text{-}36} \text{ years}$

Super-Kamiokande has set the bound $\tau(p \to e^+\pi^0) > 2.4 \times 10^{34}$ years, ruling out minimal SU(5) but consistent with SUSY GUTs where dimension-5 operators from colored Higgsino exchange dominate, giving $p \to K^+\bar{\nu}$ as the primary mode with a longer predicted lifetime.

Charge quantization: GUTs explain why$|Q_e| = |Q_p|$ exactly — quarks and leptons live in the same multiplet of the unified group, and the electric charge generator is a traceless matrix with eigenvalues that automatically satisfy charge quantization.

Computational Analysis

This simulation demonstrates gauge coupling running in the SM versus MSSM (showing unification only in the latter), quantifies the hierarchy problem by computing radiative Higgs mass corrections as a function of the cutoff scale, estimates WIMP relic abundances, and maps neutrino masses from the seesaw mechanism.

Beyond the SM: Gauge Unification, Hierarchy Problem & Dark Matter

Python

script.py238 lines

#!/usr/bin/env python3
"""
beyond_standard_model.py - BSM Physics: Hierarchy Problem & Gauge Unification
Computes gauge coupling running, hierarchy problem scales, and dark matter constraints
"""
import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import base64
from io import BytesIO

print("=" * 72)
print("  BEYOND THE STANDARD MODEL - NUMERICAL ANALYSIS")
print("=" * 72)
print()

# ============================================================
# Part 1: Gauge Coupling Unification
# alpha_i^{-1}(Q) = alpha_i^{-1}(M_Z) - b_i/(2*pi) * ln(Q/M_Z)
# ============================================================

M_Z = 91.188
alpha1_mz = 1.0 / 59.0  # U(1)_Y (GUT normalized: 5/3 * alpha')
alpha2_mz = 1.0 / 29.6  # SU(2)_L
alpha3_mz = 0.1181       # SU(3)_c

# SM one-loop beta coefficients
b1_SM = 41.0/10.0
b2_SM = -19.0/6.0
b3_SM = -7.0

# MSSM beta coefficients (with SUSY at 1 TeV)
b1_MSSM = 33.0/5.0
b2_MSSM = 1.0
b3_MSSM = -3.0

print("Gauge Coupling Unification:")
print("-" * 50)
print(f"  alpha_1^-1(M_Z) = {1/alpha1_mz:.1f}")
print(f"  alpha_2^-1(M_Z) = {1/alpha2_mz:.1f}")
print(f"  alpha_3^-1(M_Z) = {1/alpha3_mz:.1f}")
print()
print(f"  SM beta coefficients:   b1={b1_SM:.1f}, b2={b2_SM:.2f}, b3={b3_SM:.1f}")
print(f"  MSSM beta coefficients: b1={b1_MSSM:.1f}, b2={b2_MSSM:.1f}, b3={b3_MSSM:.1f}")

# Find intersection points (SM)
logQ = np.linspace(np.log10(M_Z), 18, 2000)
Q = 10**logQ

def alpha_inv(alpha_mz, b, Q, M_Z):
    return 1.0/alpha_mz - b/(2*np.pi) * np.log(Q/M_Z)

a1_inv_SM = alpha_inv(alpha1_mz, b1_SM, Q, M_Z)
a2_inv_SM = alpha_inv(alpha2_mz, b2_SM, Q, M_Z)
a3_inv_SM = alpha_inv(alpha3_mz, b3_SM, Q, M_Z)

print()
print("SM: Couplings at 10^16 GeV:")
Q_GUT = 1e16
print(f"  alpha_1^-1 = {alpha_inv(alpha1_mz, b1_SM, Q_GUT, M_Z):.1f}")
print(f"  alpha_2^-1 = {alpha_inv(alpha2_mz, b2_SM, Q_GUT, M_Z):.1f}")
print(f"  alpha_3^-1 = {alpha_inv(alpha3_mz, b3_SM, Q_GUT, M_Z):.1f}")
print("  -> DO NOT unify precisely in the SM!")
print()

# ============================================================
# Part 2: Hierarchy Problem
# ============================================================

print("Hierarchy Problem:")
print("-" * 50)
v = 246.0   # Higgs VEV (GeV)
M_Pl = 1.22e19  # Planck mass (GeV)
print(f"  Electroweak scale: v = {v} GeV")
print(f"  Planck scale: M_Pl = {M_Pl:.2e} GeV")
print(f"  Hierarchy: M_Pl/v = {M_Pl/v:.2e}")
print(f"  Fine-tuning: (v/M_Pl)^2 = {(v/M_Pl)**2:.2e}")
print()

# ============================================================
# Part 3: Dark Matter Relic Density
# ============================================================

print("Dark Matter - WIMP Miracle:")
print("-" * 50)
# Omega h^2 ~ 0.1 pb / <sigma v>
# For weak-scale particles: <sigma v> ~ alpha^2 / m^2 ~ pb for m ~ 100 GeV
m_wimp = np.array([10, 50, 100, 500, 1000])
alpha_weak = 1.0/30.0
sigma_v = alpha_weak**2 / (m_wimp**2) * 0.389e6  # convert GeV^-2 to pb (approx)
omega_h2 = 0.1 / (sigma_v / 1.0)  # rough scaling

for m, sv, oh2 in zip(m_wimp, sigma_v, omega_h2):
    print(f"  m_WIMP = {m:5d} GeV: <sigma v> ~ {sv:.2e} pb, Omega h^2 ~ {oh2:.3f}")

print(f"  Observed: Omega_DM h^2 = 0.120 +/- 0.001 (Planck)")
print()

fig, axes = plt.subplots(2, 2, figsize=(14, 11))
fig.patch.set_facecolor('#0f0a1a')

# --- Plot 1: Gauge coupling running (SM vs MSSM) ---
ax1 = axes[0, 0]
ax1.set_facecolor('#1a1025')
for spine in ax1.spines.values():
    spine.set_color('#4c1d95')
ax1.tick_params(colors='#c4b5fd')
ax1.grid(True, color='#4c1d95', alpha=0.3)

# SM running
ax1.plot(logQ, a1_inv_SM, linewidth=2, color='#a78bfa', label='alpha_1^-1 (SM)')
ax1.plot(logQ, a2_inv_SM, linewidth=2, color='#f97316', label='alpha_2^-1 (SM)')
ax1.plot(logQ, a3_inv_SM, linewidth=2, color='#22d3ee', label='alpha_3^-1 (SM)')

# MSSM running (assume SUSY threshold at 1 TeV)
M_SUSY = 1000
logQ_mssm = np.linspace(np.log10(M_SUSY), 18, 1000)
Q_mssm = 10**logQ_mssm

a1_msusy = alpha_inv(alpha1_mz, b1_SM, M_SUSY, M_Z)
a2_msusy = alpha_inv(alpha2_mz, b2_SM, M_SUSY, M_Z)
a3_msusy = alpha_inv(alpha3_mz, b3_SM, M_SUSY, M_Z)

a1_inv_MSSM = 1.0/(1.0/a1_msusy + b1_MSSM/(2*np.pi)*np.log(M_SUSY/Q_mssm))
a2_inv_MSSM = 1.0/(1.0/a2_msusy + b2_MSSM/(2*np.pi)*np.log(M_SUSY/Q_mssm))
a3_inv_MSSM = 1.0/(1.0/a3_msusy + b3_MSSM/(2*np.pi)*np.log(M_SUSY/Q_mssm))
# Fix: use proper running
a1_inv_MSSM = a1_msusy - b1_MSSM/(2*np.pi)*np.log(Q_mssm/M_SUSY)
a2_inv_MSSM = a2_msusy - b2_MSSM/(2*np.pi)*np.log(Q_mssm/M_SUSY)
a3_inv_MSSM = a3_msusy - b3_MSSM/(2*np.pi)*np.log(Q_mssm/M_SUSY)

ax1.plot(logQ_mssm, a1_inv_MSSM, linewidth=2, color='#a78bfa', linestyle='--', alpha=0.7)
ax1.plot(logQ_mssm, a2_inv_MSSM, linewidth=2, color='#f97316', linestyle='--', alpha=0.7)
ax1.plot(logQ_mssm, a3_inv_MSSM, linewidth=2, color='#22d3ee', linestyle='--', alpha=0.7)
ax1.axvline(x=np.log10(M_SUSY), color='white', linewidth=1, linestyle=':', alpha=0.3)
ax1.text(3.2, 55, 'MSSM (dashed)', color='white', fontsize=10, alpha=0.7)

ax1.set_xlabel('log10(Q / GeV)', color='#a78bfa', fontsize=11)
ax1.set_ylabel('alpha_i^{-1}(Q)', color='#a78bfa', fontsize=11)
ax1.set_title('Gauge Coupling Unification: SM vs MSSM', color='#e9d5ff', fontsize=13, fontweight='bold')
ax1.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=9)
ax1.set_ylim(0, 70)

# --- Plot 2: Higgs mass corrections (hierarchy problem) ---
ax2 = axes[0, 1]
ax2.set_facecolor('#1a1025')
for spine in ax2.spines.values():
    spine.set_color('#4c1d95')
ax2.tick_params(colors='#c4b5fd')
ax2.grid(True, color='#4c1d95', alpha=0.3)

Lambda_vals = np.logspace(3, 19, 200)  # cutoff from TeV to Planck
# delta m_H^2 ~ (lambda_t^2 / 16 pi^2) Lambda^2 (top quark loop)
yt = 1.0  # top Yukawa
delta_mH2 = (3 * yt**2 / (8 * np.pi**2)) * Lambda_vals**2

ax2.plot(np.log10(Lambda_vals), np.log10(np.sqrt(delta_mH2)), linewidth=2.5, color='#a78bfa',
         label='delta m_H from top loop')
ax2.axhline(y=np.log10(125), color='#22d3ee', linewidth=2, linestyle='--', label='Physical m_H = 125 GeV')
ax2.fill_between(np.log10(Lambda_vals), np.log10(125), np.log10(np.sqrt(delta_mH2)),
                 alpha=0.1, color='#ef4444')
ax2.text(12, 4, 'Fine-tuning
required', color='#ef4444', fontsize=12, fontweight='bold', ha='center')

ax2.set_xlabel('log10(Lambda / GeV)', color='#a78bfa', fontsize=11)
ax2.set_ylabel('log10(delta m_H / GeV)', color='#a78bfa', fontsize=11)
ax2.set_title('Hierarchy Problem: Radiative Higgs Mass', color='#e9d5ff', fontsize=13, fontweight='bold')
ax2.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=10)

# --- Plot 3: WIMP relic density ---
ax3 = axes[1, 0]
ax3.set_facecolor('#1a1025')
for spine in ax3.spines.values():
    spine.set_color('#4c1d95')
ax3.tick_params(colors='#c4b5fd')
ax3.grid(True, color='#4c1d95', alpha=0.3)

m_dm = np.logspace(0, 4, 500)
# Approximate relic density for thermal WIMP
# Omega h^2 ~ 3e-27 cm^3/s / <sigma v>
# <sigma v> ~ g^4 / (16 pi m^2) for s-wave
g_weak = 0.65
sigma_v_wimp = g_weak**4 / (16 * np.pi * m_dm**2) * (0.197e-13)**2 * 3e10  # -> cm^3/s
omega_wimp = 3e-27 / sigma_v_wimp

ax3.plot(m_dm, omega_wimp, linewidth=2.5, color='#a78bfa', label='Thermal WIMP')
ax3.axhline(y=0.120, color='#22d3ee', linewidth=2, linestyle='--', label='Planck: Omega h^2 = 0.120')
ax3.fill_between(m_dm, 0.119, 0.121, alpha=0.2, color='#22d3ee')

ax3.set_xscale('log')
ax3.set_yscale('log')
ax3.set_xlabel('WIMP mass (GeV)', color='#a78bfa', fontsize=11)
ax3.set_ylabel('Omega h^2', color='#a78bfa', fontsize=11)
ax3.set_title('WIMP Dark Matter Relic Density', color='#e9d5ff', fontsize=13, fontweight='bold')
ax3.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=10)
ax3.set_ylim(1e-4, 1e4)

# --- Plot 4: Neutrino mass scales ---
ax4 = axes[1, 1]
ax4.set_facecolor('#1a1025')
for spine in ax4.spines.values():
    spine.set_color('#4c1d95')
ax4.tick_params(colors='#c4b5fd')
ax4.grid(True, color='#4c1d95', alpha=0.3)

# Seesaw mechanism: m_nu ~ y^2 v^2 / M_R
M_R_range = np.logspace(9, 16, 200)
y_nu = [1e-1, 1e-2, 1e-3]
colors_nu = ['#a78bfa', '#f97316', '#22d3ee']

for y, col in zip(y_nu, colors_nu):
    m_nu = y**2 * v**2 / M_R_range  # in GeV
    m_nu_eV = m_nu * 1e9  # convert to eV
    ax4.plot(np.log10(M_R_range), np.log10(m_nu_eV), linewidth=2.5, color=col, label=f'y_nu = {y}')

# Observed neutrino masses
ax4.axhline(y=np.log10(0.05), color='#ef4444', linewidth=2, linestyle='--', label='m_nu ~ 0.05 eV (atmospheric)')
ax4.axhline(y=np.log10(0.009), color='#ef4444', linewidth=1.5, linestyle=':', alpha=0.7, label='m_nu ~ 0.009 eV (solar)')

ax4.set_xlabel('log10(M_R / GeV)', color='#a78bfa', fontsize=11)
ax4.set_ylabel('log10(m_nu / eV)', color='#a78bfa', fontsize=11)
ax4.set_title('Seesaw Mechanism: Neutrino Masses', color='#e9d5ff', fontsize=13, fontweight='bold')
ax4.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=9)

plt.tight_layout(pad=2.0)

buf = BytesIO()
plt.savefig(buf, format='png', dpi=150, bbox_inches='tight', facecolor='#0f0a1a')
buf.seek(0)
img_b64 = base64.b64encode(buf.read()).decode()
print(f'<img src="data:image/png;base64,{img_b64}" alt="Beyond Standard Model Calculations" />')
buf.close()

print()
print("Key insight: The SM has several shortcomings - the hierarchy problem, lack of")
print("gauge unification, neutrino masses, and dark matter - all pointing to new physics.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Additional BSM Scenarios

Extra Dimensions

The ADD model (Arkani-Hamed, Dimopoulos, Dvali, 1998) proposes $n$ large extra dimensions in which gravity propagates, while SM particles are confined to a 3+1 dimensional brane. The fundamental Planck scale $M_*$ is related to the observed Planck scale by:

$M_\text{Pl}^2 = M_*^{n+2} R^n$

For $M_* \sim 1$ TeV (eliminating the hierarchy problem), $n = 2$ extra dimensions of radius $R \sim 0.1$ mm are needed. The Randall-Sundrum model (1999) uses a single warped extra dimension to generate the hierarchy through an exponential warp factor $e^{-kR\pi}$, requiring only a modest $kR \approx 12$.

Composite Higgs

In composite Higgs models, the Higgs boson is a pseudo-Nambu-Goldstone boson (pNGB) of a new strong sector that confines at the TeV scale. The coset structure$G/H$ determines the Higgs quantum numbers — for example,$SO(5)/SO(4)$ delivers exactly one Higgs doublet. The Higgs mass is naturally light (protected by the shift symmetry) while the compositeness scale$f \sim 1$ TeV sets the cutoff for new resonances.

The key parameter is $\xi = v^2/f^2$, which controls deviations of Higgs couplings from their SM values:

$g_{hVV} = g_{hVV}^\text{SM}\sqrt{1-\xi}, \quad g_{hff} = g_{hff}^\text{SM}\frac{1-2\xi}{\sqrt{1-\xi}}$

LHC measurements constrain $\xi \lesssim 0.1$, implying $f \gtrsim 800$ GeV.

Baryogenesis

The observed baryon asymmetry $\eta_B = (n_B - n_{\bar{B}})/n_\gamma \approx 6 \times 10^{-10}$requires (Sakharov conditions): baryon number violation, C and CP violation, and departure from thermal equilibrium. The SM satisfies all three in principle (through sphalerons, CKM phase, and a first-order phase transition) but the effects are far too small. BSM sources of CP violation and/or a stronger electroweak phase transition are needed.

Experimental landscape: The search for BSM physics proceeds on multiple fronts: direct production at the LHC (new particles up to ~3 TeV), precision measurements (Higgs couplings, electroweak observables, flavor physics), rare processes (proton decay, lepton flavor violation, EDMs), and cosmological probes (dark matter detection, gravitational waves, CMB).

The Deepest Open Questions

Quantum Gravity

General relativity is non-renormalizable as a quantum field theory. New degrees of freedom must appear at or below the Planck scale $M_\text{Pl} \sim 10^{19}$ GeV. String theory and loop quantum gravity are the leading candidates, but neither has made testable predictions at accessible energies.

Cosmological Constant Problem

The observed dark energy density $\rho_\Lambda \sim (10^{-3} \text{ eV})^4$ is 120 orders of magnitude smaller than the naive QFT estimate $\rho_\text{QFT} \sim M_\text{Pl}^4$. This is arguably the worst fine-tuning problem in physics, and no convincing solution exists.

Flavor Puzzle

Why are there three generations? Why do the fermion masses span six orders of magnitude ($m_e/m_t \sim 3 \times 10^{-6}$)? Why are the CKM and PMNS mixing patterns so different? No principle within the Standard Model explains the Yukawa coupling structure.

Matter-Antimatter Asymmetry

The observed universe is overwhelmingly composed of matter. Standard Model CP violation is insufficient by many orders of magnitude. New sources of CP violation are required, perhaps connected to leptogenesis or electroweak baryogenesis.

Experimental Probes of New Physics

Energy Frontier (LHC and Beyond)

The LHC at $\sqrt{s} = 13.6$ TeV directly searches for new particles up to ~3 TeV. The HL-LHC (2029+) will increase the dataset by a factor of 10, and proposed future colliders (FCC-hh at 100 TeV, muon colliders at 10+ TeV) would dramatically extend the reach. A 100 TeV collider could probe new physics scales up to ~30 TeV and measure the Higgs self-coupling to 5%.

Precision Frontier

The muon $g-2$ anomaly ($\Delta a_\mu \sim 5\sigma$ tension with some SM calculations) could signal new physics at the TeV scale. Electric dipole moment searches probe CP violation from BSM sources up to scales of $\sim 10^3$ TeV. Lepton flavor violation searches ($\mu \to e\gamma$, $\mu \to eee$, $\mu$-to-$e$conversion) probe BSM scales up to $\sim 10^4$ TeV.

Cosmological Probes

Dark matter direct detection (LZ, XENONnT), indirect detection (CTA, Fermi-LAT), and gravitational wave observatories (LIGO/Virgo/KAGRA, LISA) probe BSM physics from vastly different angles. The CMB and large-scale structure constrain the sum of neutrino masses to $\sum m_\nu < 0.12$ eV.

Rare Processes

Proton decay searches at Hyper-Kamiokande will reach lifetimes of $\sim 10^{35}$years, probing SUSY GUT predictions. Neutrinoless double-beta decay experiments test whether neutrinos are Majorana particles. Flavor-changing neutral currents in $B$ and $K$ meson decays provide indirect sensitivity to new particles at scales up to $\sim 10^3$ TeV.

Summary: Open Questions Beyond the Standard Model

Hierarchy Problem

$\delta m_H^2 \sim \Lambda^2$ requires fine-tuning of $\sim (v/M_\text{Pl})^2 \sim 10^{-34}$. Solutions: SUSY, composite Higgs, extra dimensions.

Strong CP Problem

$|\bar{\theta}| < 10^{-10}$ with no SM explanation. Peccei-Quinn symmetry and the axion provide the most elegant solution.

Neutrino Masses

Seesaw mechanism: $m_\nu \sim m_D^2/M_R$. Right-handed neutrino mass$M_R \sim 10^{14}$ GeV points toward the GUT scale.

Dark Matter & Unification

WIMP miracle: weak-scale particles naturally give $\Omega h^2 \sim 0.1$. Gauge coupling unification works in the MSSM at $M_\text{GUT} \sim 2 \times 10^{16}$ GeV.

Extra Dimensions & Composite Higgs

ADD and Randall-Sundrum models address the hierarchy problem through extra spatial dimensions. Composite Higgs models treat the Higgs as a pseudo-Goldstone boson of a new strong sector, with $\xi = v^2/f^2$ controlling deviations from the SM.

Baryogenesis

The observed baryon asymmetry $\eta_B \approx 6 \times 10^{-10}$ requires BSM sources of CP violation and/or a stronger electroweak phase transition. Leptogenesis via heavy right-handed neutrino decay is a leading candidate mechanism.

Experimental Probes

The search proceeds on multiple fronts: energy frontier (LHC, future colliders), precision frontier (muon $g-2$, EDMs, rare decays), and cosmological probes (dark matter detection, gravitational waves, CMB). Each approach is sensitive to different BSM scenarios and energy scales.

← The Higgs Mechanism Next: Supersymmetry →