Fact-checked by Grok 2 weeks ago

Neutralino

In supersymmetric extensions of the Standard Model, such as the Minimal Supersymmetric Standard Model (MSSM), the neutralino is a hypothetical elementary particle that arises as a mass eigenstate from the mixing of the superpartners of the neutral electroweak gauge bosons (the bino and neutral wino) and the neutral components of the Higgsinos.^[1] It is a spin-1/2 Majorana fermion, meaning it is its own antiparticle and carries no electric charge, with its properties determined by model parameters including the gaugino masses M_1 and M_2, the Higgsino mass parameter \mu, and the ratio of Higgs vacuum expectation values \tan\beta.^[1] In R-parity conserving supersymmetric models, the neutralino is often the lightest supersymmetric particle (LSP), rendering it stable and unable to decay into Standard Model particles, which allows it to escape detectors while interacting weakly.^[1] This stability positions the lightest neutralino as a leading weakly interacting massive particle (WIMP) candidate for cold dark matter, potentially accounting for the observed cosmological dark matter density if its relic abundance matches \Omega h^2 \approx 0.12.^[1] The idea of neutralino dark matter was first proposed in 1983, highlighting its thermal relic production in the early universe via annihilation into Standard Model particles. Neutralinos come in four flavors in the MSSM, ordered by increasing mass, and their composition can be bino-like (dominated by the U(1)_Y gaugino superpartner), wino-like (SU(2)_L gaugino), higgsino-like, or mixed, influencing their couplings and production cross-sections at colliders.^[1] Experimental searches at the Large Hadron Collider (LHC) by ATLAS and CMS collaborations have set lower mass limits on neutralinos exceeding 300–1000 GeV in simplified models, depending on the decay channels and assumptions about other supersymmetric particles, with no direct evidence observed as of 2025.^[1] Indirect detection efforts, such as gamma-ray observations from Fermi-LAT or antimatter searches by AMS-02 on the International Space Station, probe neutralino annihilation signals in galactic dark matter halos, while direct detection experiments like XENONnT and LZ constrain spin-independent scattering cross-sections below $10^{-47} cm² for neutralino masses around 30–100 GeV.^[1] The neutralino's viability as a dark matter constituent remains a cornerstone of supersymmetry phenomenology, though tensions with collider constraints and the lack of supersymmetric particle discoveries have motivated extensions like the pMSSM (phenomenological MSSM) to accommodate lighter higgsino-like neutralinos or non-universal gaugino masses.^[1] Future prospects include high-luminosity LHC runs, which could probe neutralino masses up to 1 TeV in multi-jet plus missing energy signatures, and next-generation direct detection experiments aiming for sensitivities to $10^{-49} cm².^[1]

Definition and Basics

Definition

In supersymmetric theories, which extend the Standard Model by introducing superpartners to each known particle—pairing bosons with fermions and vice versa to achieve symmetry between matter and force carriers—the neutralino emerges as a key hypothetical particle. These superpartners address issues like the hierarchy problem and provide candidates for dark matter, with the neutralino specifically arising in the neutral sector of the theory.^[2] The neutralino is defined as a neutral, massive, spin-1/2 fermion that is a linear mixture of the superpartners of the U(1)_Y gauge boson (the bino, denoted \tilde{B}), the neutral SU(2)_L gauge boson (the neutral wino, denoted \tilde{W}^3), and the two neutral Higgsinos (denoted \tilde{H}_d^0 and \tilde{H}_u^0) from the Higgs doublets required for electroweak symmetry breaking. These four neutral fermionic states mix via a 4×4 mass matrix to form the physical neutralino mass eigenstates, labeled \tilde{\chi}_i^0 for i=1,2,3,4, ordered by increasing mass. As a Majorana fermion—meaning it is its own antiparticle—the neutralino carries no electric charge, color charge, or lepton/baryon number, making it uncolored and weakly interacting.^[2] In many supersymmetric models, particularly the Minimal Supersymmetric Standard Model (MSSM), the lightest neutralino \tilde{\chi}_1^0 is the lightest supersymmetric particle (LSP). If R-parity—a discrete symmetry conserving baryon and lepton numbers modulo 2—is conserved, the LSP is stable and cannot decay to Standard Model particles, rendering the neutralino a viable cold dark matter candidate. The term "neutralino" was coined to reflect its electric neutrality combined with the "-ino" suffix conventionally used for fermionic superpartners, analogous to names like wino, gluino, or selectron (though the latter denotes a scalar).^[2]^[3]

Historical Context

Supersymmetry emerged in the early 1970s as a theoretical framework linking bosons and fermions, initially proposed by Pierre Ramond in 1971 within the context of dual resonance models for string theory. This concept was extended to four-dimensional quantum field theories by Julius Wess and Bruno Zumino, who constructed the first supersymmetric Lagrangian in 1974, demonstrating its consistency with non-Abelian gauge interactions. The neutralino, a Majorana fermion arising as a mixture of superpartners to the gauge and Higgs bosons, was conceptualized within the Minimal Supersymmetric Standard Model (MSSM) developed in the early 1980s. This model extended the Standard Model by introducing supersymmetric partners to all particles, with the neutralino sector formalized in a comprehensive review by Howard E. Haber and Gordon L. Kane in 1985, which outlined the mass matrix and mixing for these states.^[2] Initial motivations for supersymmetry included addressing the gauge hierarchy problem, where quantum corrections would otherwise drive the Higgs mass to the Planck scale unless fine-tuned, a issue alleviated by the cancellation between bosonic and fermionic loops in SUSY models.^[4] Additionally, SUSY facilitated grand unification of the strong, weak, and electromagnetic gauge couplings at high energies, a feature absent in the non-supersymmetric Standard Model.^[5] By the late 1980s, advances in cosmology positioned the lightest supersymmetric particle—often the neutralino—as a natural candidate for non-baryonic dark matter, stable due to R-parity conservation and capable of relic densities matching observational constraints.^[6] The 1990s saw heightened interest in neutralino dark matter following the 1992 Cosmic Background Explorer (COBE) satellite detection of cosmic microwave background anisotropies, which supported the cold dark matter paradigm and required a weakly interacting massive particle to form large-scale structure. This spurred detailed calculations of neutralino relic abundances and annihilation cross-sections, with seminal work by Griest and Seckel in 1992 establishing the neutralino's viability as the primary cold dark matter component within the MSSM. Through the 2000s, these studies evolved alongside precision cosmology, reinforcing the neutralino's role in supersymmetric extensions amid ongoing searches for indirect detection signals.

Theoretical Framework

Supersymmetry Fundamentals

Supersymmetry (SUSY) is a theoretical framework in particle physics that extends the Standard Model by introducing a symmetry relating bosons and fermions, the two fundamental classes of particles distinguished by their integer and half-integer spin values, respectively.^[7] This symmetry predicts the existence of superpartners, or sparticles, for each Standard Model particle: fermionic superpartners (sfermions) for bosons and bosonic superpartners (gauginos and Higgsinos) for fermions, ensuring that the theory pairs particles of different statistics in representations of the supersymmetry algebra.^[7] In the context of models like the Minimal Supersymmetric Standard Model (MSSM), the lightest sparticle, often the neutralino, can serve as a stable dark matter candidate if certain conditions are met.^[7] The foundational structure of supersymmetry is captured by its algebra, which extends the Poincaré algebra of spacetime symmetries. For N=1 supersymmetry in four dimensions, the key anticommutation relation is

\{ Q_\alpha, \bar{Q}_{\dot{\beta}} \} = 2 (\sigma^\mu)_{\alpha \dot{\beta}} P_\mu,

where Q_\alpha and \bar{Q}_{\dot{\beta}} are the supercharges (spin-1/2 generators transforming fermions into bosons and vice versa), \sigma^\mu are the Pauli matrices extended to four dimensions, and P_\mu is the momentum operator corresponding to translations. This algebra closes under two supersymmetry transformations, yielding a spacetime translation, and implies that supersymmetric theories are invariant under these operations unless explicitly broken. The Lagrangian of a supersymmetric theory is constructed from a superpotential W(\Phi), a holomorphic function of chiral superfields \Phi, which generates fermion-boson interactions, supplemented by soft SUSY-breaking terms that introduce explicit breaking at low energies without reintroducing quadratic divergences in scalar masses.^[8]^[9] Supersymmetry addresses several shortcomings of the Standard Model, providing motivations rooted in theoretical consistency and unification. One primary motivation is the resolution of the hierarchy problem, where radiative corrections to the Higgs boson mass would otherwise receive large quadratic divergences from loops involving top quarks and gauge bosons; in SUSY, these divergences cancel precisely between contributions from bosons and their fermionic superpartners, stabilizing the electroweak scale around 100 GeV without fine-tuning.^[7] Additionally, SUSY facilitates grand unification by predicting that the three gauge couplings of the Standard Model converge at a high energy scale near 10^{16} GeV in minimal models, a feature absent in the non-supersymmetric Standard Model due to logarithmic running discrepancies.^[5] Finally, supersymmetry is essential for the consistency of string theory, where it prevents tachyonic instabilities and enables supersymmetric vacua that unify gravity with other forces at the Planck scale. A crucial feature of many supersymmetric models is the conservation of R-parity, a discrete Z_2 symmetry under which all Standard Model particles have R-parity +1 and all sparticles have R-parity -1, prohibiting sparticle decay into solely Standard Model particles and ensuring the stability of the lightest supersymmetric particle (LSP).^[7] This conservation arises naturally in certain grand unified extensions but can be imposed ad hoc in minimal models to avoid rapid proton decay.^[10] To reconcile supersymmetry with the observed absence of sparticles at low energies and the lack of exact degeneracy between particles and superpartners, soft SUSY-breaking terms—such as gaugino masses, scalar masses, and trilinear couplings—are introduced in the Lagrangian, parameterized by a few dozen parameters in the MSSM and arising from higher-scale dynamics like supergravity or string theory.^[9] These terms break supersymmetry spontaneously or explicitly while preserving the theory's ultraviolet finiteness and predictive power.^[11]

Neutralinos in the MSSM

The Minimal Supersymmetric Standard Model (MSSM) represents the simplest supersymmetric extension of the Standard Model, incorporating two Higgs doublets to generate masses for all fermions while preserving supersymmetry and ensuring gauge anomaly cancellation.^[3] Unlike the single Higgs doublet of the Standard Model, the MSSM requires separate doublets H_d (with hypercharge Y = -1/2) for down-type fermion masses and H_u (with Y = +1/2) for up-type fermion masses, leading to four neutral fermionic components in the superpartner sector: the bino \tilde{B} (the fermionic partner of the U(1)_Y gauge boson), the neutral wino \tilde{W}^3 (the third component of the SU(2)_L gaugino triplet), and the neutral Higgsinos \tilde{H}_d^0 and \tilde{H}_u^0 (the fermionic partners of the neutral components of H_d and H_u).^[3] These fields mix through electroweak symmetry breaking, parameterized by the vacuum expectation values v_d and v_u of the Higgs doublets, to form four neutralinos, the Majorana mass eigenstates of the theory.^[12] The mixing among these neutral components is governed by a symmetric 4×4 mass matrix M_N in the gauge-eigenstate basis (\tilde{B}, \tilde{W}^3, \tilde{H}_d^0, \tilde{H}_u^0), with diagonal elements set by the soft supersymmetry-breaking gaugino masses M_1 (for the bino) and M_2 (for the wino), and the off-diagonal elements originating from the Higgsinos' interactions via the supersymmetric Higgsino mass parameter \mu and the ratio \tan\beta = v_u / v_d.^[3] The electroweak contributions to the off-diagonal terms involve the Z-boson mass m_Z and the weak mixing angle \theta_W, reflecting the breaking of SU(2)_L × U(1)_Y symmetry. This matrix encapsulates the full neutralino sector dynamics within the MSSM at tree level.^[12] The neutralino mass matrix takes the explicit form

M_N = \begin{pmatrix} M_1 & 0 & -m_Z s_W c_\beta & m_Z s_W s_\beta \\ 0 & M_2 & m_Z c_W c_\beta & -m_Z c_W s_\beta \\ -m_Z s_W c_\beta & m_Z c_W c_\beta & 0 & -\mu \\ m_Z s_W s_\beta & -m_Z c_W s_\beta & -\mu & 0 \end{pmatrix},

where s_W = \sin\theta_W, c_W = \cos\theta_W, s_\beta = \sin\beta, and c_\beta = \cos\beta.^[3] The eigenvalues and eigenvectors of M_N determine the neutralino masses and compositions, respectively. To obtain the physical states, M_N is diagonalized via a unitary matrix N satisfying N^* M_N N^\dagger = \operatorname{diag}(m_{\tilde{\chi}_1^0}, m_{\tilde{\chi}_2^0}, m_{\tilde{\chi}_3^0}, m_{\tilde{\chi}_4^0}), where the masses m_{\tilde{\chi}_i^0} (for i = 1 to $4) are ordered increasingly, and \tilde{\chi}_1^0 denotes the lightest neutralino, often the lightest supersymmetric particle (LSP) in models with conserved R-parity.^[12] The mixing matrix N parameterizes the extent to which each neutralino is bino-like, wino-like, or Higgsino-like, depending on the relative sizes of M_1, M_2, \mu, and \tan\beta.

Physical Properties

Mass and Mixing

In the Minimal Supersymmetric Standard Model (MSSM), the masses of the four neutralinos are obtained by diagonalizing the neutralino mass matrix, which depends primarily on the soft supersymmetry-breaking gaugino mass parameters M_1 (for the U(1)_Y bino) and M_2 (for the SU(2)_L wino), the Higgsino mass parameter \mu, and the ratio \tan \beta = v_u / v_d of the vacuum expectation values of the up-type and down-type Higgs doublets.^[13] These parameters typically range from hundreds of GeV to a few TeV, with M_1 and M_2 often related by grand unification assumptions such as M_1 \approx 0.5 M_2 at the electroweak scale, while \tan \beta spans values from about 2 to 60 to accommodate electroweak symmetry breaking and Higgs phenomenology.^[13] The lightest neutralino mass m_{\tilde{\chi}_1^0}, frequently considered the lightest supersymmetric particle (LSP), has a theoretical range from roughly 1 GeV up to several TeV, but collider experiments impose stringent lower bounds.^[13] Experimental constraints from the Large Hadron Collider (LHC), analyzed by ATLAS and CMS collaborations, exclude m_{\tilde{\chi}_1^0} \lesssim 100 GeV for stable bino-like neutralinos in simplified models where sleptons or other superpartners mediate production, based on searches for events with missing transverse energy, jets, and leptons using up to 140 fb^{-1} of 13 TeV data from Run 2, with ongoing Run 3 analyses as of 2025.^[14] In more general phenomenological MSSM scans, the lower limit on m_{\tilde{\chi}_1^0} can reach 200–300 GeV or higher depending on the superpartner spectrum and decay assumptions, with no signals observed as of 2025.^[14] Earlier limits from LEP experiments set a model-independent bound of m_{\tilde{\chi}_1^0} > 46 GeV for stable neutralinos, but LHC results have significantly tightened constraints in viable SUSY scenarios.^[14] The neutralino mass eigenstates \tilde{\chi}_i^0 (with i = 1 to $4, ordered by increasing mass) are mixtures of the gauge eigenstates, expressed as

\tilde{\chi}_i^0 = N_{i1} \tilde{B} + N_{i2} \tilde{W}^3 + N_{i3} \tilde{H}_d^0 + N_{i4} \tilde{H}_u^0,

where N is the $4 \times 4 unitary mixing matrix that diagonalizes the neutralino mass matrix, and the coefficients N_{ij} determine the composition and thus the interaction strengths of each eigenstate.^[13] The mixing is governed by the relative magnitudes of M_1, M_2, and |\mu| compared to the electroweak scale; for instance, if |M_1| \ll |\mu|, M_2, the LSP is predominantly bino-like (|N_{11}| \approx 1), leading to a relatively light state with suppressed couplings to gauge bosons and fermions due to the bino's hypercharge nature.^[13] Common neutralino scenarios include the bino-like LSP, which favors lighter masses (often 100–500 GeV in viable parameter space) and weaker electroweak interactions, making it a motivated dark matter candidate but challenging to produce directly at colliders.^[13] In contrast, a higgsino-like LSP (when |\mu| \ll M_1, M_2) results in heavier masses (typically above 300–1000 GeV to evade bounds) with enhanced couplings to W and Z bosons owing to the Higgsino components (|N_{i3}|, |N_{i4}| \approx 1/\sqrt{2}), leading to nearly degenerate multiplets.^[13] Mixed compositions arise when M_1, M_2 \sim |\mu|, yielding intermediate properties where the mixing elements N_{ij} balance gaugino and Higgsino contributions, influencing both mass splittings and phenomenology across the spectrum.^[13]

Composition and Interactions

Neutralinos are spin-1/2 Majorana fermions, meaning they are self-conjugate particles that obey Fermi-Dirac statistics and possess both particle and antiparticle properties within the same state.^[3] In the Minimal Supersymmetric Standard Model (MSSM), each neutralino is a linear superposition of the neutral gaugino states—the bino \tilde{B} (superpartner of the U(1)_Y gauge boson) and the neutral wino \tilde{W}^3 (superpartner of the SU(2)_L gauge boson)—and the two neutral higgsino states \tilde{H}_d^0 and \tilde{H}_u^0 (superpartners of the Higgs doublets).^[13] This mixing arises from the neutralino mass matrix, which is diagonalized to yield the physical mass eigenstates \tilde{\chi}_i^0 (for i=1,2,3,4), with the composition determined by the unitary mixing matrix N, such that \tilde{\chi}_i^0 = N_{i1} \tilde{B} + N_{i2} \tilde{W}^3 + N_{i3} \tilde{H}_d^0 + N_{i4} \tilde{H}_u^0.^[3] The composition of a neutralino varies based on the relative scales of the supersymmetry-breaking parameters, particularly the gaugino masses M_1 and M_2, and the higgsino mass parameter [\mu](/page/MU). For low-mass neutralinos, particularly the lightest one (\tilde{\chi}_1^0), a bino-dominated composition is common when |M_1| \ll |[\mu](/page/MU)|, |M_2|, reflecting the weaker weak hypercharge coupling g' associated with the U(1)_Y gauge group.^[3] In contrast, wino-dominated neutralinos emerge when |M_2| \ll |M_1|, |[\mu](/page/MU)|, leading to stronger electroweak interactions via the SU(2)_L coupling g.^[13] Higgsino-dominated neutralinos, prevalent when |[\mu](/page/MU)| \ll |M_1|, |M_2|, exhibit Yukawa-like couplings similar to those of the Higgs sector, enhancing interactions with Higgs bosons and fermions.^[3] Neutralinos participate in tree-level interactions with the Z-boson (axial-vector type), Higgs bosons (scalar type), and sfermions (via gaugino-fermion-sfermion vertices), but they have no tree-level coupling to the photon due to their electric neutrality.^[3] The Z-boson coupling strength for neutralinos \tilde{\chi}_i^0 and \tilde{\chi}_j^0 is proportional to the higgsino mixing terms, specifically g_{\tilde{\chi}_i^0 \tilde{\chi}_j^0 Z} \propto (N_{i3} N_{j3} - N_{i4} N_{j4}), reflecting the difference between the down-type and up-type higgsino asymmetries in the weak current.^[13] This interaction is captured in the effective Lagrangian term \mathcal{L} \supset g_{\tilde{\chi}_i^0 \tilde{\chi}_j^0 Z} Z_\mu \bar{\tilde{\chi}}_i^0 \gamma^\mu \gamma^5 \tilde{\chi}_j^0, where the axial-vector structure arises from the Majorana nature of the neutralinos.^[3] Higgs and sfermion couplings further depend on the respective bino, wino, or higgsino fractions, with strengths scaled by the gauge couplings g' or g for gauginos and by Yukawa couplings for higgsinos.^[13]

Phenomenological Aspects

Production at Colliders

Neutralinos, as the lightest supersymmetric partners in many models, are primarily produced at high-energy colliders through indirect cascade decays of heavier supersymmetric particles or via direct electroweak pair production processes.^[1] In the Minimal Supersymmetric Standard Model (MSSM), associated production often involves the creation of squarks or sleptons, which subsequently decay into a neutralino and a standard model fermion; for example, a squark \tilde{q} decays via \tilde{q} \to q + \tilde{\chi}^0, where q is a quark and \tilde{\chi}^0 denotes a neutralino.^[15] Similarly, gluino pair production, a dominant strong-interaction process at the Large Hadron Collider (LHC), leads to neutralinos through cascade decays such as \tilde{g} \to q \bar{q} \tilde{\chi}^0, where \tilde{g} is the gluino, contributing significantly to event topologies with multiple jets and missing transverse energy.^[16] These cascades are simulated using event generators like MadGraph for matrix elements and Pythia for parton showers, providing theoretical cross sections at next-to-leading order (NLO) plus next-to-leading logarithmic (NLL) accuracy.^[17] Direct pair production of neutralinos, such as pp \to \tilde{\chi}^0_i \tilde{\chi}^0_j + X (where i, j = 1, 2, \ldots), proceeds via electroweak s-channel processes mediated by W^\pm, Z, or Higgs bosons, with cross sections typically suppressed compared to strong production modes.^[1] For instance, the cross section for \tilde{\chi}_1^0 \tilde{\chi}_2^0 production approximates \sigma(pp \to \tilde{\chi}_1^0 \tilde{\chi}_2^0) \sim \alpha^2 s / M^2 in the high-energy limit for s-channel dominance, where \alpha is the fine-structure constant, s is the center-of-mass energy squared, and M represents the relevant mediator or neutralino mass scale; numerical evaluations using tools like Resummino yield values on the order of picobarns for TeV-scale masses at LHC energies.^[18] These electroweak processes are particularly relevant for heavier neutralino states and are also modeled with MadGraph/Pythia frameworks to incorporate higher-order corrections.^[17] Kinematic considerations impose thresholds for neutralino production, requiring the collider center-of-mass energy \sqrt{s} to exceed $2 m_{\tilde{\chi}} for on-shell pair production, though cascades from heavier particles relax this for the lightest neutralino.^[1] At the LHC, Run 3 operations (2022–2025) at \sqrt{s} = 13.6 TeV enable probing of neutralino masses up to approximately 1 TeV in certain cascade scenarios involving colored superpartners, with sensitivity diminishing for compressed mass spectra where decay products carry low momentum.^[1]

Decay Modes

In the Minimal Supersymmetric Standard Model (MSSM), neutralinos heavier than the lightest supersymmetric particle (LSP), denoted as \tilde{\chi}_i^0 for i \geq 2, primarily decay via two-body channels into the LSP \tilde{\chi}_1^0 plus a Standard Model boson, provided the kinematics allow it (i.e., m_{\tilde{\chi}_i^0} > m_{\tilde{\chi}_1^0} + m_V where V is the boson mass).^[1] The dominant modes include \tilde{\chi}_i^0 \to \tilde{\chi}_1^0 + Z, \tilde{\chi}_i^0 \to \tilde{\chi}_1^0 + h (lightest Higgs), or \tilde{\chi}_i^0 \to \tilde{\chi}_1^0 + \gamma (radiative), with branching ratios strongly influenced by the neutralino mixing compositions, such as bino-, wino-, or higgsino-like dominance.^[19] If two-body decays are kinematically suppressed, three-body channels like \tilde{\chi}_i^0 \to \tilde{\chi}_1^0 + f \bar{f} (via virtual Z or sfermions) become relevant, though their rates are typically smaller.^[1] The decay widths for these processes depend on the mass splitting \Delta m = m_{\tilde{\chi}_i^0} - m_{\tilde{\chi}_1^0} and the effective couplings from the neutralino-Higgsino-gaugino mixing matrix. For higgsino-like neutralinos, the partial width for \tilde{\chi}_2^0 \to \tilde{\chi}_1^0 + h is approximately given by

\Gamma(\tilde{\chi}_2^0 \to \tilde{\chi}_1^0 h) \approx \frac{g^2}{16\pi} \frac{(m_{\tilde{\chi}_2^0}^2 - m_{\tilde{\chi}_1^0}^2)^2}{m_{\tilde{\chi}_2^0}^3},

valid in the limit where the Higgs mass is negligible compared to the neutralino masses; this yields widths on the order of 0.1–0.5 GeV for \Delta m \sim 100 GeV and m_{\tilde{\chi}_2^0} \sim 500 GeV, with one-loop corrections modifying the tree-level result by up to 10–30% in complex MSSM scenarios.^[19] Overall lifetimes for these non-LSP neutralinos are extremely short, typically \tau \sim 10^{-24}–$10^{-20} s for typical mass splittings, leading to prompt decays within collider detectors.^[1] In scenarios with R-parity violation (RPV), the LSP neutralino \tilde{\chi}_1^0 itself can decay, enabling invisible channels such as \tilde{\chi}^0 \to \nu \bar{\nu} (via bilinear RPV) or \tilde{\chi}^0 \to e^+ e^- \nu (via slepton exchange), though these are suppressed in minimal RPV extensions of the MSSM due to small coupling constants \lambda, \lambda'.^[20] The resulting widths are typically \Gamma \lesssim 10^{-15} GeV, far below standard MSSM rates, and branching ratios remain small unless RPV parameters are tuned large.^[1] Long-lived or metastable neutralinos arise in parameter regions with small mass splittings (\Delta m \sim 1–10 GeV) or suppressed couplings, such as in gauge-mediated supersymmetry breaking (GMSB) models where the next-to-LSP is nearly degenerate with the LSP. In these cases, decay lengths can reach c\tau \sim 0.1–20 m for certain mixings and masses below 300 GeV, producing displaced vertices observable at colliders; for example, ATLAS and CMS exclude lifetimes of 10–2000 cm for m_{\tilde{\chi}_1^0} = 300 GeV in such scenarios.^[1]

Dark Matter Role

Relic Density

In the standard cosmological model, the lightest neutralino (χ) is a viable cold dark matter candidate due to its stability from R-parity conservation and weak-scale interactions, allowing it to decouple via thermal freeze-out in the early universe. Initially in thermal equilibrium with the plasma, neutralinos remain in chemical equilibrium as long as their annihilation rate exceeds the Hubble expansion rate H. Freeze-out occurs when the temperature drops to T_f ≈ m_χ / 20 (for m_χ the neutralino mass), at which point the number density n_χ becomes Boltzmann-suppressed, leaving a relic abundance that matches the observed dark matter density Ω h^2 ≈ 0.120 ± 0.001 if the thermally averaged annihilation cross-section times relative velocity ⟨σ v⟩ ≈ 3 × 10^{-9} GeV^{-2}.^[21] The required ⟨σ v⟩ is achieved through neutralino pair annihilation into Standard Model particles, with dominant channels including χχ → W^+ W^-, ZZ, hh (where h is the Higgs boson), and ff̄ (fermion-antifermion pairs) for light quarks and leptons. These processes proceed via t-channel sfermion exchange, s-channel gauge or Higgs boson exchange, and gauge interactions, with the cross-section depending sensitively on the neutralino's bino-wino-higgsino composition. In scenarios where the neutralino is nearly degenerate with other supersymmetric particles, coannihilation effects—such as with light staus (τ̃) or stops (t̃)—enhance the effective annihilation rate by including processes like χ τ̃ → W τ or χ t̃ → q χ̃^±, allowing viable relic densities even for lighter neutralinos.^[22] The relic density is computed by solving the Boltzmann equation for the neutralino number density evolution:

\frac{d n_\chi}{d t} = -3 H n_\chi - \langle \sigma v \rangle (n_\chi^2 - (n_\chi^\mathrm{eq})^2),

where H is the Hubble parameter and n_χ^eq is the equilibrium density; this integro-differential equation is typically solved numerically over cosmic time using dedicated codes like micrOMEGAs, which incorporate full annihilation and coannihilation networks in the Minimal Supersymmetric Standard Model (MSSM).^[23] Achieving the observed relic density requires tuning supersymmetric parameters to ensure efficient annihilation without excessive fine-tuning. For instance, a nearly pure higgsino-like neutralino with mass m_χ ≈ 1.1 TeV yields Ω h^2 ≈ 0.12 through coannihilation with the nearly degenerate chargino and next-to-lightest neutralino, while pure wino-like neutralinos, which would require m_χ ≈ 3 TeV to account for Sommerfeld-enhanced electroweak annihilations, are disfavored by recent indirect detection constraints from gamma-ray observations.^[24] The neutralino mass and mixing, which dictate the coupling strengths and available channels, thus play a central role in matching cosmological observations.

Detection Prospects

Direct detection of neutralinos as dark matter candidates primarily involves observing elastic scattering processes, denoted as \chi N \to \chi N, where \chi represents the neutralino and N a nucleus in the target material. These interactions can proceed via spin-independent channels, dominated by Higgs boson exchange in the t-channel, or spin-dependent channels mediated by Z-boson exchange.^[25] For bino-like neutralinos, the spin-independent cross-section is typically on the order of \sigma_{SI} \sim 10^{-47} cm², though it varies with model parameters such as the neutralino mass and mixing angles. Experiments employing xenon or other noble gases, such as XENONnT and LZ, aim to probe these low cross-sections by achieving sensitivities down to approximately $10^{-48} cm² with projected exposures, while next-generation detectors like DARWIN or XLZD are expected to reach $10^{-49} cm² in the 2030s.^[26]^[27] Indirect detection strategies target annihilation products from neutralino pairs in dense astrophysical environments, such as galactic halos or the Sun.^[28] Gamma-ray signals from processes like \chi \chi \to \gamma \gamma produce monochromatic lines at energy E = m_\chi, observable by the Fermi-LAT telescope, while continuum emissions arise from quark and lepton final states.^[29] Positron excesses potentially attributable to neutralino annihilations into lepton pairs are probed by the AMS-02 experiment on the International Space Station.^[30] Neutrino telescopes like IceCube search for high-energy neutrinos from neutralino annihilations captured in the Sun or Earth, focusing on muon neutrino signatures.^[31] Neutralinos consistent with the observed cosmic relic density provide a benchmark for these searches, as their annihilation rates must align with thermal freeze-out requirements. Future collider experiments offer complementary precision probes of neutralino properties through indirect signatures, such as measurements of the Higgs boson's invisible decay width, which can constrain light neutralino scenarios. Facilities like the International Linear Collider (ILC) and the electron-positron stage of the Future Circular Collider (FCC-ee) are projected to achieve percent-level precision on electroweak parameters sensitive to supersymmetric extensions, enabling tests of neutralino mixing and masses. Additionally, astrophysical effects such as Sommerfeld enhancement—arising from long-range forces between neutralinos at low velocities—can boost annihilation signals in regions like the galactic center by factors of up to $10^2 or more, improving detection prospects in indirect searches.

Experimental Status

Collider Constraints

Collider searches provide some of the strongest experimental constraints on neutralinos, primarily through the absence of supersymmetric signals in missing transverse energy (MET) signatures from cascade decays ending in the stable lightest neutralino \tilde{\chi}_1^0. These limits are derived from data at e^+e^- colliders like LEP and hadron colliders including the LHC, focusing on production modes such as pair production of colored superpartners (gluinos, squarks) or electroweakinos decaying to \tilde{\chi}_1^0 plus visible particles like jets or leptons. At LEP, operating at center-of-mass energies up to \sqrt{s} = 209 GeV, the ALEPH, DELPHI, L3, and OPAL experiments set a model-independent lower mass limit of m_{\tilde{\chi}_1^0} > 46.3 GeV at 95% confidence level from direct searches for e^+e^- \to \tilde{\chi}_1^0 \tilde{\chi}_1^0 Z and associated productions, assuming R-parity conservation and a stable \tilde{\chi}_1^0. This bound, close to the kinematic threshold m_Z/2 \approx 45.3 GeV, is robust in the minimal supersymmetric Standard Model (MSSM) with gaugino mass unification, reaching up to 94 GeV for scenarios with M_2 < 1 TeV and |\mu| \leq 2 TeV without third-generation mixing. Pure bino-like neutralinos face even stronger exclusions, with m_{\tilde{\chi}_1^0} > 100 GeV in certain parameter spaces, derived from the lack of acoplanar lepton or jet events with significant missing energy. The LHC's ATLAS and CMS experiments have extended these constraints using proton-proton collisions at \sqrt{s} = 13-13.6 TeV and integrated luminosities up to \sim 140 fb^{-1} from Runs 1 and 2, with Run 3 data collected through 2025 further refining the bounds without observing any supersymmetric signals. In simplified models, where gluinos or squarks decay promptly to \tilde{\chi}_1^0 plus quarks or gluons, no excesses in MET + jets channels exclude gluino masses up to 2.4 TeV and squark masses up to 1.9 TeV for m_{\tilde{\chi}_1^0} \lesssim 500 GeV, translating to effective lower limits on m_{\tilde{\chi}_1^0} > 200-500 GeV depending on the mass hierarchy and branching ratios. For electroweak production, such as \tilde{\chi}_1^\pm \tilde{\chi}_2^0 pairs decaying to \tilde{\chi}_1^0 + leptons + MET, ATLAS and CMS exclude next-to-lightest neutralino masses up to 600 GeV assuming a massless \tilde{\chi}_1^0, with compressed spectra (mass splittings \Delta m \sim 10 GeV) challenging detection but closing the slepton gap from LEP by excluding sleptons up to 250 GeV. In the phenomenological MSSM (pMSSM), these results exclude bino-like \tilde{\chi}_1^0 masses below ~100 GeV across broad parameter scans. By 2025, Particle Data Group summaries confirm no viable light bino-like neutralinos below 100 GeV in key channels, with ongoing Run 3 analyses at higher luminosities expected to probe deeper into compressed regions.

Direct and Indirect Searches

Direct detection experiments targeting neutralino dark matter as weakly interacting massive particles (WIMPs) have yielded null results, imposing stringent limits on spin-independent (SI) scattering cross-sections. The XENONnT experiment, utilizing over 1 tonne of liquid xenon, reported no excess events in its 2025 analysis of dark matter-electron and nuclear interactions, excluding SI cross-sections above approximately $1.7 \times 10^{-47} cm² for a 30 GeV neutralino, with limits around $10^{-47} cm² across 30-100 GeV, approaching the neutrino fog background.^[32] Similarly, the LUX-ZEPLIN (LZ) collaboration's 2025 results from 4.2 tonne-years of exposure set world-leading constraints, surpassing prior exclusions by a factor of four for WIMP masses above 9 GeV/c² and ruling out σ_SI > 2.2 × 10^{-48} cm² at 40 GeV in the 30-100 GeV range for neutralino-like candidates.^[33] The PandaX-4T experiment's ~1.5 tonne-year dataset from 2025 further corroborates these null findings, with comparable exclusions around $10^{-47} cm² for low-mass neutralinos through searches for both nuclear recoils and light dark matter interactions.^[34] Indirect detection efforts probe neutralino annihilation products, primarily gamma rays and neutrinos, from astrophysical sources. The Fermi Large Area Telescope (Fermi-LAT) analysis of dwarf spheroidal galaxies in 2025, combining over 16 years of data, established upper limits on the velocity-averaged annihilation cross-section ⟨σv⟩ below 10^{-25} cm³/s at 95% confidence level for neutralino masses around 10-100 GeV, depending on annihilation channels.^[35] No significant gamma-ray excess attributable to neutralino annihilation has been confirmed in the galactic center, with recent morphological studies challenging dark matter interpretations of the observed GeV excess due to inconsistencies with expected neutralino signals. Complementarily, IceCube's 2025 search for neutrinos from WIMP annihilation in the Sun and Earth's core, using ten years of data, provides bounds on ⟨σv⟩ for neutralino masses above 100 GeV, with spin-independent scattering limits competitive with direct experiments in the multi-TeV regime.^[36] As of 2025, viable parameter space remains for neutralino dark matter, particularly Higgsino-like candidates around 100 GeV and wino-like at the TeV scale, which evade current bounds while matching relic density requirements through co-annihilation or non-standard cosmology. However, tensions arise if the galactic center gamma-ray excess is ascribed to neutralino annihilation, as the required ⟨σv⟩ exceeds limits from dwarf galaxies by up to an order of magnitude. Multi-messenger approaches, including LIGO's gravitational wave constraints on dark matter distributions around black holes, offer marginal additional limits on neutralino-primordial black hole mixtures but do not significantly restrict the core parameter space.

References

[1]
[PDF] 89. Supersymmetry, Part II (Experiment) - Particle Data Group
May 31, 2024 · A NLSP neutralino will decay to a gravitino and a SM particle whose nature is determined by the neutralino composition. Final states with ...
[2]
[PDF] arXiv:1307.3561v2 [hep-ph] 4 Nov 2013
Nov 4, 2013 · In the mini- mal supersymmetric standard model (MSSM), the light- est neutralino is a Majorana fermion, and is a mixture of the ...
[3]
[PDF] PROBING PHYSICS BEYOND THE STANDARD MODEL
behaves like a gluino or a neutralino points strongly toward a supersymmetric ... Haber and G.L. Kane, Thesearch for supersymmetry: Probing physics beyond the ...
[4]
[1302.6587] Naturalness and the Status of Supersymmetry - arXiv
Feb 26, 2013 · In this review, we begin by recalling the theoretical motivations for weak-scale supersymmetry, including the gauge hierarchy problem, grand unification, and ...
[5]
[hep-ph/9308233] Supersymmetry with Grand Unification - arXiv
Aug 5, 1993 · SUSY stabilizes scalar mass corrections (the hierarchy problem), greatly reduces the number of free parameters, facilitates gauge coupling unification.Missing: motivations | Show results with:motivations
[6]
[hep-ph/9604232] Weak-Scale Supersymmetry: Theory and Practice
Apr 3, 1996 · These lectures contain an introduction to the theory and practice of weak-scale supersymmetry. They begin with a discussion of the hierarchy problem.
[7]
[PDF] 88. Supersymmetry, Part I (Theory) - Particle Data Group
May 31, 2024 · 88.1 Introduction. Supersymmetry (SUSY) is a generalization of the space-time symmetries of quantum field the- ory that transforms fermions ...Missing: paper | Show results with:paper
[8]
[PDF] Introduction to Supersymmetry - CERN Indico
The single-particle states of the theory fall into groups called supermultiplets, which are turned into each other by Q and Q†. Fermion and boson members of a ...
[9]
Soft supersymmetry breaking terms from supergravity ... - Inspire HEP
We review the origin of soft supersymmetry-breaking terms in N=1 supergravity models of particle physics. We first consider general formulae for those terms ...
[10]
[hep-ph/9312232] Natural Conservation of R Parity in Supersymmetry
Dec 6, 1993 · For a more satisfactory realization of supersymmetry, we propose here a new model which conserves R automatically. It is unifiable under SO(14) ...
[11]
[2005.13577] A string landscape guide to soft SUSY breaking terms
May 27, 2020 · We examine several issues pertaining to statistical predictivity of the string theory landscape for weak scale supersymmetry (SUSY).
[12]
https://pdg.lbl.gov/2023/reviews/rpp2023-rev-susy-1-theory.pdf
[13]
[PDF] 88. Supersymmetry, Part I (Theory) - Particle Data Group
Dec 1, 2023 · Consequently, there are five neutralino mass eigenstates that are obtained by a Takagi- diagonalization of the 5 × 5 neutralino mass matrix.
[14]
None
Summary of each segment:
[15]
MSSM neutralino properties at the LHC from gluino/squark cascade ...
Dec 5, 2005 · Here the heavy states are neutralinos (i = 2,3,4) - note that is excluded - which are produced in gluino/squark ( /\tilde{q}) cascade decay ...
[16]
[hep-ph/0501157] Pair-produced heavy particle topologies:MSSM ...
Pair-produced heavy particle topologies:MSSM neutralino properties at the LHC from gluino/squark cascade decays. Authors:M. Bisset, N. Kersting, J. Li (Tsinghua ...
[17]
SUSYCrossSections < LHCPhysics < TWiki
### Summary of Neutralino Production Cross Sections at LHC
[18]
Pair production of neutralinos and charginos at the LHC
Dec 14, 2011 · Separate cross sections for neutralino-pair production σ ( χ ˜ 1 0 χ ˜ 2 0 ) ( p b ) at the LHC with s = 14 TeV as a function of tan ⁡ β (left) ...
[19]
[1208.4106] Neutralino Decays in the Complex MSSM at One-Loop
Aug 20, 2012 · Abstract:We evaluate two-body decay modes of neutralinos in the Minimal Supersymmetric Standard Model with complex parameters (cMSSM).
[20]
Neutralino Decay Rates with Explicit R-parity Violation - hep-ph - arXiv
Sep 24, 1997 · We compute the neutralino decay rate in the minimal supersymmetric standard model with the addition of explicit R-parity violation. We include ...
[21]
[1807.06209] Planck 2018 results. VI. Cosmological parameters - arXiv
Jul 17, 2018 · A combined analysis gives dark matter density \Omega_c h^2 = 0.120\pm 0.001, baryon density \Omega_b h^2 = 0.0224\pm 0.0001, scalar spectral ...
[22]
[hep-ph/9704361] Neutralino Relic Density including Coannihilations
Apr 20, 1997 · We evaluate the relic density of the lightest neutralino, the lightest supersymmetric particle, in the Minimal Supersymmetric extension of the Standard Model ( ...
[23]
A program for calculating the relic density in the MSSM - hep-ph - arXiv
Dec 20, 2001 · We present a code that calculates the relic density of the LSP in the MSSM. All tree-level processes for the annihilation of the LSP are included.
[24]
[astro-ph/0402033] Wimp/Neutralino Direct Detection - arXiv
Feb 2, 2004 · This review describes the most competitive and promising experiments with different detection techniques.
[25]
Projected WIMP Sensitivity of the XENONnT Dark Matter Experiment
Jul 17, 2020 · With the exposure goal of 20 t\,y, the expected sensitivity to spin-independent WIMP-nucleon interactions reaches a cross-section of 1.4\times10 ...Missing: LZ | Show results with:LZ
[26]
Indirect detection of dark matter with γ rays - PNAS
Indirect detection searches for the products of WIMP annihilation or decay. This is generally done through observations of γ-ray photons or cosmic rays.
[27]
In Wino Veritas? Indirect Searches Shed Light on Neutralino Dark ...
Jul 16, 2013 · By combining results from Fermi-LAT and HESS, we show that: light nonthermal wino dark matter is strongly excluded; thermal wino dark matter is ...
[28]
Dark matter for excess of AMS-02 positrons and antiprotons - arXiv
Apr 29, 2015 · We propose a dark matter explanation to simultaneously account for the excess of antiproton-to-proton and positron power spectra observed in the AMS-02 ...Missing: neutralino | Show results with:neutralino
[29]
[1312.4445] IceCube, DeepCore, PINGU and the indirect search for ...
Dec 16, 2013 · In this study, we re-evaluate and give an update on the prospects for detecting supersymmetric neutralinos with neutrino telescopes, focussing ...
[30]
[2410.17036] Dark Matter Search Results from 4.2 Tonne-Years of ...
Oct 22, 2024 · The LUX-ZEPLIN experiment found no evidence for an excess over expected backgrounds in its dark matter search, with world-leading constraints ...Missing: limits neutralino
[31]
Combined dark matter search towards dwarf spheroidal galaxies ...
Aug 27, 2025 · Depending on the DM content modeling, the 95% confidence level observed limits reach 1.5\times10^{-24} cm^3s^{-1} and 3.2\times10^{-25} cm^3s^{- ...
[32]
[2412.12972] Search for dark matter from the center of the Earth with ...
Dec 17, 2024 · ... IceCube ... Our upper limits on the spin-independent WIMP scattering are world-leading among neutrino telescopes for WIMP masses m_{\chi}>100~GeV.Missing: 2025 | Show results with:2025