---
title: 'SUBMET: J-PARC Millicharged Particle Search'
url: https://www.emergentmind.com/topics/sub-millicharge-experiment-submet
type: topic
---

# SUBMET: J-PARC Millicharged Particle Search

The SUB-Millicharge ExperimenT (SUBMET) is a dedicated experimental program for the search for millicharged particles, $\chi$, at J-PARC, centered on the detection of extremely small ionization signals produced by particles with electric charge $Q=\epsilon e$ and sub-GeV to GeV-scale masses. In its J-PARC realization, SUBMET is a proton fixed-target experiment installed in the Neutrino Monitor building about $280$ m downstream of the graphite target, with detector concepts and projections developed for $N_{\rm POT}=10^{22}$ and an intended reach below $Q=10^{-3}e$ for $m_\chi<1.6~{\rm GeV}/c^2$ [2102.11493, 2507.19311].

## 1. Development of the program

SUBMET has been documented through a sequence of design stages rather than a single immutable hardware configuration. The 2020 Letter of Intent proposed a four-layer plastic-scintillator detector at the B2 location of the Neutrino Monitor building, about $280$ m from the target and about $30$ m underground, with projected $95\%$ CL sensitivity down to $\epsilon \simeq 3\times10^{-4}$ for $m_\chi \lesssim 0.2~{\rm GeV}/c^2$ and $\epsilon \simeq 1.5\times10^{-3}$ for $m_\chi \lesssim 1.6~{\rm GeV}/c^2$ at $N_{\rm POT}=10^{22}$ [2007.06329]. A 2021 feasibility study then reformulated SUBMET as a two-layer coincidence detector, explicitly using beam timing from the J-PARC bunch structure to suppress dark-count backgrounds, and quoted improved sensitivity to $\epsilon \simeq 5\times10^{-5}$ for $m_\chi<0.2~{\rm GeV}/c^2$ and $\epsilon \simeq 8\times10^{-4}$ for $m_\chi<1.6~{\rm GeV}/c^2$ [2102.11493]. A later detector-design paper describes the mechanically integrated apparatus, reports that data taking began in June 2024, and states that the experiment is optimized for particles with $m_\chi<1.6~{\rm GeV}/c^2$ and $Q<10^{-3}e$ [2507.19311].

| Stage | Detector concept | Quoted reach or status |
|---|---|---|
| 2020 LOI [2007.06329] | Four layers of BC-408 bars | $\epsilon \simeq 3\times10^{-4}$ below $0.2~{\rm GeV}/c^2$; $\epsilon \simeq 1.5\times10^{-3}$ below $1.6~{\rm GeV}/c^2$ |
| 2021 feasibility study [2102.11493] | Two-layer stacked scintillator bars with beam timing | $\epsilon \simeq 5\times10^{-5}$ below $0.2~{\rm GeV}/c^2$; $\epsilon \simeq 8\times10^{-4}$ below $1.6~{\rm GeV}/c^2$ |
| 2025 mechanical integration note [2507.19311] | Built two-layer EJ-200 detector with 160 modules | Data taking began June 2024 |

The differing layer counts, bar dimensions, and sensitivity numbers are therefore features of successive published design stages. A common misunderstanding is to treat all quoted SUBMET configurations as simultaneous descriptions of one finalized detector. The publications instead document optimization from the LOI to the feasibility study and then to an installed instrument.

## 2. Production channels and detection principle

In the J-PARC fixed-target program, SUBMET searches for millicharged particles produced by a $30$ GeV proton beam on graphite. The early sensitivity studies used neutral-meson and charmonium production as the baseline source term. The 2020 Letter of Intent treated $\pi^0$, $\eta$, and $J/\psi$ channels, with yields parameterized as
$$
N_\chi \propto c_{\mathfrak m}\,\epsilon^2 N_{\rm POT}\times f(m_\chi^2/m_{\mathfrak m}^2),
$$
where $f$ vanishes as $m_\chi \to m_{\mathfrak m}/2$ [2007.06329]. The 2021 feasibility study retained the meson-dominated production picture and used PYTHIA8 inputs $c_{\pi^0}\approx 1.9$, $c_\eta\approx 0.21$, and $c_{J/\psi}\approx 5\times10^{-9}$ per POT [2102.11493]. The later detector-design paper states the broader motivation in terms of proton-target production of $\pi^0$, $\eta$, $\rho$, $\omega$, $\phi$, and $J/\psi$, with the resulting $\chi$ flux forward-biased along the beam axis [2507.19311].

Detection relies on the fact that the ionization signal is quadratically suppressed by the fractional charge. The detector-design paper writes
$$
(dE/dx)_\chi \approx (Q/e)^2 \times (dE/dx)_e,
$$
which is the central reason SUBMET is optimized for single- to few-photoelectron response [2507.19311]. In the 2021 feasibility study, a GEANT4 optical simulation for a through-going unit-charge MIP gave
$$
N_{\rm PE}(\epsilon)=N_{\rm PE}(\epsilon=1)\,\epsilon^2 \approx 2.5\times10^5 \epsilon^2
$$
per layer, and the two-layer detection probability was modeled as
$$
P_{\rm det}(\epsilon)=[1-e^{-\mu}]^2,\qquad \mu=N_{\rm PE}(\epsilon).
$$
In the sub-millicharge regime $\mu\ll1$, this implies $P_{\rm det}\approx \mu^2$ and, because the accepted $\chi$ flux also scales as $\epsilon^2$, the signal yield scales as $N_{\rm sig}\propto\epsilon^6$ [2102.11493]. This strong scaling is why SUBMET optimization studies consistently identify photoelectron yield as the dominant experimental lever.

## 3. Installed detector, geometry, and mechanical integration

The mechanically integrated SUBMET detector described in 2025 is a two-layer system of long plastic scintillator bars coupled to photomultiplier tubes, positioned $280$ m downstream from the target and about $1^\circ$ below the beam axis. The bars are tilted upward by $4.5^\circ$ relative to the floor to align with the expected $\chi$ arrival direction, and the coincidence signature is a left-to-right traversal through both layers within a narrow time window [2507.19311].

Each layer contains $10\times8$ bars, for $80$ modules per layer and $160$ modules in total. The bars are EJ-200 plastic scintillators with dimensions $50\times50\times1500~{\rm mm}^3$, and the reported material properties are scintillation efficiency $10{,}000$ photons/MeV for a $1$ MeV electron, emission peak $425$ nm, attenuation length $3{,}800$ mm, rise time $0.9$ ns, and decay time $2.1$ ns [2507.19311]. Optical wrapping consists of two layers of $100~\mu$m Teflon, one layer of $18~\mu$m aluminum foil, and one layer of $160~\mu$m insulating tape, increasing the thickness by about $1$ mm. One Hamamatsu R7725 PMT is attached to one end of each bar using optical grease, and an LED circuit is housed at the opposite end for testing and calibration [2507.19311].

The mechanical structure is hierarchical. Modules are grouped into $2\times2$ supermodules, supermodules are mounted into a cage holding twenty such units, and two aluminum-extrusion tables support the two detector layers. The cage with supermodules has dimensions $1{,}850\times764\times620~{\rm mm}$ [2507.19311]. Finite-element analyses quoted in the design paper give a maximum von Mises stress of about $15$ MPa, a safety factor of about $15$, static displacement not exceeding $1$ mm for Table 1, and worst-case seismic displacements of about $1.14$ mm for Table 1 and $1.08$ mm for Table 2 [2507.19311]. These figures are presented as evidence that bar-to-bar registration and PMT coupling remain stable under static and seismic loads.

Magnetic shielding is implemented with two layers of $0.2$ mm mu-metal around the four PMTs in each supermodule. The reported attenuation of perpendicular external fields exceeds $95\%$, PMT gain stability was confirmed up to $2.5$ mT, and cosmic-muon pulse height was unaffected in that field range [2507.19311]. Before shipment, all $160$ modules were validated with $\mathcal O(100)$ reconstructed cosmic-muon tracks; the system was then disassembled, transported, and reassembled at the J-PARC Neutrino Monitor building [2507.19311].

## 4. Timing logic, backgrounds, and afterpulse modeling

The dominant background identified in the 2021 feasibility study is random coincidence between two layers due to PMT dark counts. The J-PARC beam structure is therefore integral to the SUBMET strategy. The beam has spills at $1.16$ s repetition, with eight bunches separated by $600$ ns, and site timing signals are used to gate the trigger and suppress random coincidences by a factor of order $10^{-6}$ [2102.11493]. In that study, a conservative total background budget of $N_{\rm bkg}=5$ over three years was adopted, and the main qualitative conclusion was that beam timing can render dark-count coincidences negligible at the analysis level [2102.11493].

The 2025 afterpulse study adds a more detailed treatment of PMT-intrinsic delayed backgrounds. It considers the installed detector as $160$ EJ-200 plus R7725 modules and reports that afterpulses follow large PMT pulses, especially primaries with peak height greater than $140$ mV. Because the J-PARC bunch spacing is about $600$ ns, afterpulses generated by one bunch can fall into later bunch windows and mimic single-photoelectron signals [2510.19724]. The waveform readout window is $5~\mu$s, and dedicated June 2024 data were taken with the trigger delay shifted so that the $\sim4~\mu$s region after the last bunch was recorded. The analysis used $50$k events, divided equally into parameter extraction and validation samples [2510.19724].

The modeling window is defined by
$$
t\in[t_0,t_1]=[720~{\rm ns},\,3475.2~{\rm ns}],
$$
excluding the first $720$ ns because of strong baseline fluctuations immediately after a large pulse. The instantaneous afterpulse rate is modeled as
$$
\frac{dn}{dt}=\frac{1}{\tau}n_0(A)e^{-t/\tau},
$$
where $A$ is the primary pulse area, $\tau$ is a per-module time constant extracted from the $\Delta t$ distribution, and $n_0(A)$ is fixed by the measured afterpulse yield in $[t_0,t_1]$ [2510.19724]. Two functional forms for the $A$-dependence were tested,
$$
n_{\rm pred}^{\rm lin}(A|t_0,t_1)=p_0^{\rm lin}+p_1^{\rm lin}A,
$$
and
$$
n_{\rm pred}^{\rm expo}(A|t_0,t_1)=\exp(p_0^{\rm expo}+p_1^{\rm expo}A).
$$
The resulting prediction for an arbitrary time interval $[t_i,t_f]$ with $t_i\ge720$ ns is obtained by integrating the exponential model over that window [2510.19724].

The paper reports that the model reproduces observed afterpulse rates with approximately $20\%$ precision and that both the linear and exponential-in-area parameterizations perform comparably, with the linear model adopted for simplicity [2510.19724]. Beam-off data showed dark current and radiation-induced rates of order $10^3$ Hz, negligible compared with the order-$10^6$ Hz afterpulse rates in the analysis windows. An operationally important result is that more than $70\%$ of collision events contain at least one large pulse; the afterpulse model therefore allows SUBMET to retain a substantial fraction of beam-on data that would otherwise be vetoed [2510.19724].

## 5. Sensitivity, optimization, and experimental context

The quoted SUBMET reach depends on which design stage is considered. The 2020 LOI projected $95\%$ CL sensitivity to $\epsilon \simeq 3\times10^{-4}$ for $m_\chi\lesssim0.2~{\rm GeV}/c^2$ and $\epsilon \simeq 1.5\times10^{-3}$ for $m_\chi\lesssim1.6~{\rm GeV}/c^2$ with $N_{\rm POT}=10^{22}$ [2007.06329]. The 2021 feasibility study, using a two-layer coincidence design and explicit beam timing, improved these numbers to $\epsilon \simeq 5\times10^{-5}$ for $m_\chi<0.2~{\rm GeV}/c^2$ and $\epsilon \simeq 8\times10^{-4}$ for $m_\chi<1.6~{\rm GeV}/c^2$ [2102.11493]. The 2025 detector-design paper repeats the same associated-JHEP sensitivity figures for a three-year run with $N_{\rm POT}=10^{22}$ and frames them as the design target of the installed apparatus [2507.19311].

Optimization studies in both early and later design stages converge on the same conclusion: photoelectron yield is the main driver. In the 2021 feasibility study, doubling $N_\chi$ in acceptance improved $\epsilon_{95}$ by less than $10\%$, while doubling $N_{\rm PE}$ improved $\epsilon_{95}$ by about $20\%$ [2102.11493]. The 2020 LOI reached a similar conclusion with a different baseline detector: doubling $N_\chi$ improved $\epsilon$ by about $10\%$, increasing the background by $100\times$ worsened $\epsilon$ by about $15\%$, and increasing $N_{\rm PE}$ by factors of $2$ and $3$ improved $\epsilon$ by about $25\%$ and $35\%$, respectively [2007.06329]. These comparisons are design-specific, but they all point to the same scaling logic in the sub-millicharge regime.

Within the broader experimental landscape, SUBMET occupies a proton fixed-target scintillation niche. The NA64 program at CERN is based on missing energy or missing momentum in $\sim100$ GeV electron and muon beams and projects $Q_\chi/e \approx \mathcal O(10^{-3})$ in muon mode and $\mathcal O(10^{-2})$ in electron mode around $m_\chi\approx200$ MeV [1810.06856]. The SENSEI search in the NuMI beam used a low-threshold silicon detector in the MINOS cavern and set world-leading constraints for masses between $30$ MeV and $380$ MeV, improving previous bounds by as much as a factor of $2$ in the $100$–$210$ MeV interval [2305.04964]. A later Jefferson Lab Hall A proposal with a $2\times14$ Skipper-CCD array reports that such a configuration can exceed the sensitivity of all existing searches below $1.5$ GeV and is competitive with SUBMET across the overlapping mass range [2412.09652]. SUBMET is therefore best understood as one member of a wider mCP search ecosystem that includes proton fixed-target scintillation detectors, missing-energy experiments, and ultralow-threshold silicon beam-dump searches.

## 6. Broader uses of the SUBMET label

Later literature uses the SUBMET label in a broader conceptual sense than the specific J-PARC scintillator apparatus. One strand is the terrestrial-accumulation program for Earth-bound millicharged particles. In that setting, underground electrostatic accelerators are proposed as devices that directly accelerate ambient thermalized mCPs by
$$
E_{\rm gain}=\epsilon e\,\Delta V,
$$
with MV-scale facilities such as LUNA, CASPAR, and JUNA cited as examples [2012.03957]. For a model beamline with $L=1$ m, radius $r=1$ mm, and depth of about $1$ km, the flux of accelerated particles above threshold is written
$$
\Phi[E>E_{\rm thr}] = 2\pi rL\left(1-\frac{E_{\rm thr}}{\epsilon e\Delta V}\right)\frac{f_Q\eta\rho_Q}{m_Q}v_{\rm th}\,
\min\!\left[1,\frac{r\,\epsilon e \Delta V}{L\sqrt{T\,E_{\rm thr}}}\right],
$$
and the study states that, for top-down accumulation with $m_Q\gtrsim1$ GeV, detectable rates can probe $f_Q$ down to about $10^{-20}$ with realistic backgrounds of order $10^{-3}$ Hz [2012.03957]. This is not the same configuration as the J-PARC detector, but the acronym is explicitly used there for an underground accelerator-based search strategy.

A second strand appears in the reinterpretation of Cavendish tests of Coulomb’s law as both accumulators and detectors for room-temperature mCP populations. That proposal specifies an outer accumulator shell with
$$
R_{\rm trap}=2~{\rm m},\qquad \phi_{\rm trap}=-1~{\rm MV},
$$
and an inner Cavendish shell with
$$
R_0=1~{\rm m},\qquad \phi_0=50~{\rm kV},\qquad \nu_0=250~{\rm Hz},
$$
together with a solid inner sphere of radius $0.5$ m and low vacuum $P\simeq10^{-6}$ atm [2510.25825]. The signal is an oscillating interior potential difference,
$$
\Delta\phi_\chi(t)=\frac{\rho_\chi(t)}{6}(r_2^2-r_1^2),
$$
and the paper argues that an accumulator-enhanced Cavendish configuration could reach ambient densities as low as $n_\chi\sim10^{-16}~{\rm cm}^{-3}$ and may detect the irreducible cosmic-ray mCP population for sub-GeV masses [2510.25825].

These later uses of the acronym show that “SUBMET” can denote either the specific J-PARC scintillator experiment or a broader family of terrestrial mCP search concepts. In the experimental literature, however, the unqualified name most commonly refers to the J-PARC program documented in the Letter of Intent, the feasibility study, the mechanical integration note, and the afterpulse-background analysis [2007.06329, 2102.11493, 2507.19311, 2510.19724].

Source: https://www.emergentmind.com/topics/sub-millicharge-experiment-submet