Pith. Sign in

REVIEW

Higgsino Above the Sea of Fog

High Energy Physics - Phenomenology · Submitted 1 Sep 2026

Pith's verdict: Well-argued theoretical speculation built on a 2.6-sigma event; honest about uncertainties, but the central assumption is fragile and the IceCube bound is a real tension.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 8, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.01504 (version 3)
Pith IDF7QKVG5B, Pith's permanent record of this paper and its review
CitationJiJi Fan, Matthew Reece. "Higgsino Above the Sea of Fog." arXiv:2609.01504 (2026). https://pith.science/paper/F7QKVG5B
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

The paper argues that if LZ's single 248 keV nuclear recoil is real dark matter, it is the inelastic scattering of a higgsino, with the consequence that electroweak gauginos and heavy Higgs scalars must sit near $10^6$–$10^7$ GeV.

The reading

The paper takes a single high-energy nuclear recoil event—$248 \pm 23$ keV, recorded by the LZ experiment in 2.84 tonne-years of data—and asks what it would mean if it is the first direct-detection sign of supersymmetric dark matter. Working under the paper's own stated assumption that the event is a genuine dark matter interaction rather than a background fluctuation, it argues that the natural explanation is inelastic scattering of higgsino dark matter: a roughly 1 TeV relic that scatters by jumping to a slightly heavier partner state, which requires a mass splitting of about 350 keV under the standard galactic halo model and closer to 500 keV if the Large Magellanic Cloud reshapes the local high-speed velocity tail. The central consequence is a concrete, testable spectrum: that small splitting forces the electroweak gauginos and the heavy Higgs scalars of supersymmetry up to $10^6$–$10^7$ GeV while keeping only the higgsino light. The paper then shows how such a spectrum can be built (suppressing the $\mu$ term without suppressing $b_\mu$), surveys the cosmological and experimental implications, and lays out the complementary probes—solar capture and IceCube neutrinos, gamma-ray telescopes, future colliders, and electron electric dipole moments—that would confirm or kill the picture.

What carries the argument

The load-bearing object is the split higgsino: a pair of nearly degenerate Majorana fermions, $\tilde\chi_1$ (the dark matter) and $\tilde\chi_2$ (its excited twin), split by a small mass $\delta$ and coupled to nuclei only through the off-diagonal $Z$-boson interaction, which makes every scattering an endothermic excitation—the inelastic dark matter mechanism. The identity that carries the argument is the tree-level splitting formula $\delta \approx m_Z^2(\sin^2\theta_W/M_1 + \cos^2\theta_W/M_2)$, which converts a measured recoil energy into a demand on the electroweak gaugino mass scale and gives $M_2 \approx 10^7$ GeV for $\delta = 600$ keV. Two further pieces do the quantitative work: the halo integral $\eta(v_{\min}) = \int_{v>v_{\min}} d^3v\, f(\vec v + \vec v_e)/v$, which turns the velocity requirement into an event rate and is where the astrophysical uncertainty enters, and the one-loop threshold correction to $\mu$, which ties the heavy Higgs mass $m_H$ to the gaugino scale and rules out a conventional split-SUSY spectrum.

Extended reading notes

Core claim

Assuming the LZ event at $E_R = 248 \pm 23\,(\mathrm{stat}) \pm 23\,(\mathrm{sys})$ keV is a dark matter nuclear recoil, the paper establishes that standard elastic scattering cannot account for it: a 250 keV event would be accompanied by thousands of lower-energy events, whereas the inelastic mechanism naturally produces a spectrum that starts only near the mass splitting. In the higgsino realization, the two nearly degenerate neutral mass eigenstates $\tilde\chi_1$ and $\tilde\chi_2$ are split by $\delta$, and $Z$-boson exchange couples them off-diagonally, so each scattering is an excitation. The minimum incoming speed, $v_{\min} \approx 350\sqrt{E_R/250\,\mathrm{keV}} + 425\sqrt{250\,\mathrm{keV}/E_R}\,(\delta/350\,\mathrm{keV})$ km/s, puts a $\delta \approx 350$ keV event on the extreme high-speed tail of the galactic halo, $v_{\min} \gtrsim 760$ km/s, where the halo integral is uncertain by more than an order of magnitude between the standard halo model and an LMC-modified distribution. A maximum-likelihood fit in the $(|\mu|,\delta)$ plane gives $\delta \approx 300$–$380$ keV for the standard halo model—$\delta \approx 350$ keV at the 1.1 TeV thermal benchmark—and $\delta \approx 480$–$500$ keV for the rough LMC model. The paper's decisive step is then the tree-level relation $\delta \approx m_Z^2(\sin^2\theta_W/M_1 + \cos^2\theta_W/M_2)$: a few-hundred-keV splitting demands $M_1, M_2 \sim 10^6$–$10^7$ GeV, and the one-loop correction to $\mu$ forces the heavy Higgs doublet to the same scale. The result is a heavy-but-not-split supersymmetric spectrum—wino, bino, and heavy Higgs near $10^6$–$10^7$ GeV, higgsino near a TeV—with mild tuning in either the $M_2/M_1$ ratio or the $m_H/M_2$ ratio.

Load-bearing premise

The whole chain rests on two stacked premises: that the single LZ recoil is a genuine dark matter interaction rather than a background or statistical fluctuation (the paper's own global significance is only 2.6 sigma), and that the extreme high-speed tail of the local dark matter velocity distribution, which simulations do not reliably resolve above roughly 950 km/s, determines the event rate.

Editorial extensions

If this is right

  • The supersymmetric spectrum is essentially fixed: wino, bino, and heavy Higgs doublet near $10^6$–$10^7$ GeV with only the higgsino around 1 TeV, and the measured 125 GeV Higgs mass then independently requires small $\tan\beta$ (roughly 1–2) at such scales, a condition the paper verifies with a dedicated benchmark calculation.
  • The LZ-favored splitting sits right against the strongest existing bound: solar capture followed by annihilation would produce high-energy neutrinos, and the paper's companion analysis [30] argues IceCube already excludes splittings below about 566 keV, so additional IceCube exposure can confirm or exclude the scenario.
  • The heavier higgsino state decays radiatively with a decay length of order a kilometer ($c\tau \approx 1600$ km at $\delta = 600$ keV, $\mu = 1.1$ TeV), so upscattering in the Earth followed by decay inside large neutrino detectors—'luminous dark matter'—plus heavy-element paleodetectors can probe the same splitting range with independent systematics.
  • Collider searches are immune to the halo-model uncertainty: a 100 TeV proton or 10 TeV muon collider could directly produce the higgsino and its charged partners up to the thermal benchmark, and the paper reviews that HL-LHC mono-X searches reach only about 200–300 GeV, far short of 1.1 TeV.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A free test the paper leaves implicit: inelastic scattering rates peak annually near June 2, and LZ logged its event on June 16, 2023; if more high-recoil events accumulate, the distribution of their event dates across the year would discriminate the excitation mechanism from background or elastic explanations.
  • If the event is real, it inverts a question in galactic astronomy: instead of asking what the high-velocity tail predicts for direct detection, the event demands that Milky Way–Large Magellanic Cloud simulations reliably resolve the local velocity distribution above roughly 950 km/s—a concrete deliverable from existing simulation programs that the paper identifies as missing.
  • Among the paper's UV completions, the single-modulus route to suppressing $\mu$ requires the Kähler potential exponent $c = 12$, yet the cited quantum-gravity bounds cap known examples at $c \le 7$; if the higgsino interpretation survives, the viable completion must therefore be multi-modulus with approximately integer F-term ratios, a restriction that could discriminate between string constructio
Original abstract

Higgsino dark matter is a classic, but still viable, weakly-interacting massive particle (WIMP) scenario that could point to supersymmetry, a new spacetime symmetry of nature. It is a simple realization of the "inelastic dark matter" mechanism, where its scattering off a nucleus occurs between two nearly degenerate states. Intriguingly, very recently the LZ collaboration reported a single nuclear recoil event at a high recoil energy. Assuming that this event is a dark matter candidate, it could be explained by higgsino dark matter, e.g., a 1.1 TeV thermal higgsino with a mass splitting of order hundreds of keV. The mass splitting required is very sensitive to the extreme high speed tail of the local dark matter velocity distribution, and thus susceptible to astrophysical uncertainties: the standard halo model prefers ~350 keV, while a model incorporating Large Magellanic Cloud effects prefers ~500 keV (a rough estimate based on partial information about a particular simulation). A much better understanding of the high speed tail is needed to evaluate whether the parameter space is consistent with an IceCube constraint on the annihilation of dark matter captured by the Sun. In any case, the mass splitting of several hundred keV implies that the electroweak gauginos and heavy Higgs scalars in supersymmetry have to be at a much higher scale ~10^{3-4} TeV. We discuss possible UV completions to generate a SUSY spectrum in which all superpartners are at a common scale, except for a much lighter higgsino mass $\mu$. Besides the thermal benchmark, we comment on possible non-thermal cosmological scenarios of higgsino dark matter with gravitino and modulus decays. The higgsino dark matter scenario could also be searched for at other complementary probes such as indirect detection and future colliders. One striking LZ event may herald an exciting future of further experimental developments.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 8,486 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.78 · human participants 0.01 · controlled comparison 0.05 · randomized assignment 0.01 · states quantitative result 0.90 · compares with prior results 0.45 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

Higgsino dark matter is a classic, but still viable, weakly-interacting massive particle (WIMP) scenario that could point to supersymmetry, a new spacetime symmetry of nature. It is a simple realization of the "inelastic dark matter" mechanism, where its scattering off a nucleus occurs between two nearly degenerate states. Intriguingly, very recently the LZ collaboration reported a single nuclear recoil event at a high recoil energy. Assuming that this event is a dark matter candidate, it could be explained by higgsino dark matter, e.g., a 1.1 TeV thermal higgsino with a mass splitting of order hundreds of keV. The mass splitting required is very sensitive to the extreme high speed tail of the local dark matter velocity distribution, and thus susceptible to astrophysical uncertainties: the standard halo model prefers ~350 keV, while a model incorporating Large Magellanic Cloud effects prefers ~500 keV (a rough estimate based on partial information about a particular simulation). A much better understanding of the high speed tail is needed to evaluate whether the parameter space is consistent with an IceCube constraint on the annihilation of dark matter captured by the Sun. In any case, the mass splitting of several hundred keV implies that the electroweak gauginos and heavy Higgs scalars in supersymmetry have to be at a much higher scale ~10^{3-4} TeV. We discuss possible UV completions to generate a SUSY spectrum in which all superpartners are at a common scale, except for a much lighter higgsino mass $\mu$. Besides the thermal benchmark, we comment on possible non-thermal cosmological scenarios of higgsino dark matter with gravitino and modulus decays. The higgsino dark matter scenario could also be searched for at other complementary probes such as indirect detection and future colliders. One striking LZ event may herald an exciting future of further experimental developments.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Higgsino dark matter is a classic, but still viable, weakly-interacting massive particle (WIMP) scenario that could point to supersymmetry, a new spacetime symmetry of nature. background 1.00 0.04 0.06 hedged possibility 0.82
s1 It is a simple realization of the "inelastic dark matter" mechanism, where its scattering off a nucleus occurs between two nearly degenerate states. background 1.00 0.07 0.04 asserted 1.00
s2 Intriguingly, very recently the LZ collaboration reported a single nuclear recoil event at a high recoil energy. background 0.90 0.04 0.03 asserted 0.99
s3 Assuming that this event is a dark matter candidate, it could be explained by higgsino dark matter, e.g., a 1.1 TeV thermal higgsino with a mass splitting of order hundreds of keV. interpretation 0.79 0.16 0.21 hedged possibility 1.00
s4 The mass splitting required is very sensitive to the extreme high speed tail of the local dark matter velocity distribution, and thus susceptible to astrophysical uncertainties: the standard halo model prefers ~350 keV, while a model incorporating Large Magellanic Cloud effects prefers ~500 keV (a rough estimate based on partial information about a particular simulation). limitation 0.49 0.18 0.06 asserted 0.78
s5 A much better understanding of the high speed tail is needed to evaluate whether the parameter space is consistent with an IceCube constraint on the annihilation of dark matter captured by the Sun. future work 0.41 0.09 0.05 asserted 0.71
s6 In any case, the mass splitting of several hundred keV implies that the electroweak gauginos and heavy Higgs scalars in supersymmetry have to be at a much higher scale ~10^{3-4} TeV. interpretation 0.90 0.35 0.34 asserted 0.92
s7 We discuss possible UV completions to generate a SUSY spectrum in which all superpartners are at a common scale, except for a much lighter higgsino mass $\mu$. contribution 0.98 0.59 0.14 hedged possibility 0.48
s8 Besides the thermal benchmark, we comment on possible non-thermal cosmological scenarios of higgsino dark matter with gravitino and modulus decays. contribution 0.95 0.42 0.04 hedged possibility 0.86
s9 The higgsino dark matter scenario could also be searched for at other complementary probes such as indirect detection and future colliders. future work 0.95 0.06 0.06 hedged possibility 1.00
s10 One striking LZ event may herald an exciting future of further experimental developments. interpretation 0.69 0.06 0.10 hedged possibility 0.88

What kind of paper the classifier reports this is: theoretical result 0.53 · evidence from argument 0.43.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. vmin at LZ as a function of δ for mχ = 500 GeV, 1.1 TeV and 10 TeV, fixing ER = 250 keV. Note that vmin as a function of ER at a given δ would be minimized at a recoil energy ER = µN mN δ , (2) which is approximately δ when mχ ≫ mN . This suggests that the typical recoil energy is essentially determined by δ, independently of the dark matter mass and the target nucleus [18]. Given the single event at 248±23(stas)±23(sys) keV, we expect that δ is also around this range. From
Figure 2. Averaged halo integral as a function of vmin for vmin > 600 km/s, for both SHM (black, solid) and LMC velocity distribution given in [25] (blue, dashed). equilibrium processes, could significantly modify the high-velocity tail. For example, it has been argued in [24–26] that the interaction between Large Magellanic Cloud (LMC) and Milky Way could have a major impact on the local dark matter velocity distribution. More specifically, simulations show that dark matter particles in the solar neighborhood from the LMC analogue could be the dominant source for the high speed tail of the…
Figure 3. Recoil spectra as a function of ER, assuming the SHM (left) and LMC (right) halo models. The blue curve corresponds to δ = 0 while the yellow one is for δ = 350 keV. We fix mχ = 1.1 TeV and σn = 7.4 × 10−39 cm2 . 2 Interestingly, using the traditional Helm form factor [29] will increase the needed mass splitting by about (5–10) keV, depending on the astrophysical model. It is a subdominant effect compared to the large astrophysi￾cal uncertainty, but may give a sense of how nuclear uncertainties could matter. 3As pointed out by [30], this value is corrected by a factor of 4 compared…
Figure 4. 1σ contour that could account for the LZ event in the (|µ|, δ) plane for higgsino dark matter with σn = 7.4 × 10−39 cm2 , assuming the SHM (left) and a rough LMC model (right). See the main text for more details of the LMC model and how to interpret it correctly. The red dashed line indicates the thermal higgsino benchmark at 1.1 TeV. We maximize the log of the extended likelihood and plot the 1σ contour in the (|µ|, δ) plane, as shown in
Figure 5. The amount of fine tuning in the ratio M2/M1 needed to produce a given splitting δ at a given value of M2. Without tuning, δ of a few hundred keV requires M2 ∼ 107 GeV. Allowing a factor of 10 tuning, M2 can be closer to the PeV scale. Notice that δ can be made arbitrarily small by tuning the ratio r21 = M2/M1 close to a critical value r ∗ 21 = − 1 tan2 θW ≈ −3.3. (14) Fitting the available data with relatively small M1,2 requires tuning r21 close to this ratio. A natural measure of this tuning is simply ∆r21 ≡
Figure 6. The magnitude of the 1-loop correction |∆µ| to the higgsino mass, for fixed choices of M2 (labeled curves), M1 chosen to reproduce the correct δ given M2, and varying ratio mH/M2. We have fixed tan β = 2 and δ = 600 keV. than the gauginos. In particular, the existence of a light Higgs boson and consistent electroweak symmetry breaking constrains how large bµ can be. The one-loop threshold correction is ∆µ ≈ − sin 2β 64π 2  3g 2M2 m2 H m2 H − M2 2 log m2 H M2 2 + g ′2M1 m2 H m2 H − M2 1 log m2 H M2 1  , (16) with mH the heavy higgs boson mass scale.5 The factor sin 2β becomes small…
Figure 7. Value of tan β needed for mh = 125 GeV, as a function of M2, with error band, computed with SUSYHD. We have fixed µ = 1.1 TeV, mQ3 = mU3 = M3 = 10M2, At = 0, and adjusted M1 as a function of M2 to achieve δ = 600 keV. to tune it against). This still requires tuning mH ≲ 0.1M1 to get small enough ∆µ, and so is no less tuned than the scenario where both the bino and wino play a role. The stability of higgsino dark matter implies the need for a symmetry to forbid decays like He0 1 → ν + Z. In particular, there are gauge invariant, lepton number violating terms of the form µLHuL that…
Figure 8. Running of the Standard Model gauge couplings at one loop in the model with µ = 1.1 TeV and other superpartners at the scale MSUSY = 3 × 106 GeV. be unsuppressed in our scenario. To see this, recall that the mass matrix for the Higgs doublets is
Figure 9. Predicted electron EDM de as a function of M2, with M1,2 taken to be real, |µ| = 1.1 TeV, arg(µ) = π/4, tan β = 2, and M1 chosen to achieve δ = 600 keV. (γh) diagram in the limit |M2| ≫ |µ| ≫ mh is [134] de e ≈ αg2me sin(2β) sin(ϕµ) 32π 3M2µ  log µ 2 m2 h + 2 ≈ 1.4 × 10−31 sin(ϕµ) 106 GeV M2 e cm. (51) A plot of de as a function of M2, fixing tan β = 2, ϕµ = π/4, and tuning M1 to achieve δ = 600 keV, is shown in
Figure 10. Normalized recoil spectra dR/dER for the three 4-fermi effective operators L1 = i( ¯χχ)(NN¯ ), L2 = ¯χχN γ¯ 5N, and L4 = ¯χγ5χN γ¯ 5N. The curves are normalized to have the same area. For a given operator coefficient, the rate of L4 is strongly suppressed relative to that of L2, which is in turn strongly suppressed relative to L1. they also strongly suppress the scattering rate. In

Figures reproduced from arXiv:2609.01504 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.