Letter of Intent: Search for sub-millicharged particles at J-PARC

We propose a new experiment sensitive to the detection of millicharged particles produced at the $30$ GeV proton fixed-target collisions at J-PARC. The potential site for the experiment is B2 of the Neutrino Monitor building, $280$ m away from the target. With $\textrm{N}_\textrm{POT}=10^{22}$, the experiment can provide sensitivity to particles with electric charge $3\times10^{-4}\,e$ for mass less than $0.2$ $\textrm{GeV}/\textrm{c}^2$ and $1.5\times10^{-3}\,e$ for mass less than $1.6$ $\textrm{GeV}/\textrm{c}^2$. This brings a substantial extension to the current constraints on the charge and the mass of such particles.

1 Motivation for the experiment Electric charge quantization is a long-standing question in particle physics.While fractionally charged particles (millicharged particles hereafter) have typically been thought to preclude the possibility for Grand Unified Theories (GUTs), well-motivated dark-sector models have been proposed to predict the existence of millicharged particles while preserving the possibility for unification.Such models can contain a rich internal structure, providing candidate particles for dark matter.One well motivated mechanism that leads to millicharged particles is to introduce a new U (1) in dark sector with a massless dark-photon (A ) and a massive dark-fermion (χ) [1], where the last term shows that A and B kinetically mix with the mixing parameter κ.After replacing / A with / A +κB µ , the coupling between the dark fermion and B becomes κe ( κe χ / Bχ).
The charge of χ varies by the size of mixing, so χ can be a millicharged particle.Hereafter, χ is used to denote millicharged particles.
A number of experiments have searched for millicharged particles, including in an electron fixed-target experiment [2], proton-proton colliders [3][4][5], and neutrino experiments [6].A comprehensive review is in Reference [7].In the parameter space of the charge (Q) and the mass (m χ ), the region of m χ > 0.1 GeV/c 2 and Q < 10 −2 e is largely unexplored.
Proton fixed-target experiments provide a solid testing ground for χs.The particle flux is much larger than the collider experiments and they can reach higher energy regime than electron fixed-target experiments.The sensitivity of such experiments to χs can reach beyond Q χ ∼ 10 −3 e for a wide mass range from a few MeV/c 2 to a couple of GeV/c 2 .We propose an experiment, SUBMET (SUB-Millicharge ExperimenT), which utilizes the 30 GeV J-PARC proton beam to search for χs in this unexplored region. 2 Production of millicharged particles at J-PARC We consider χs with charge Q produced from π 0 , η and J/ψ neutral meson decays.The Υ production is irrelevant because the center-of-mass energy √ s is 7.5 GeV for the collisions between the 30 GeV proton beam and the fixed target.The lighter mesons (m = π 0 , η) decay through photons (π 0 , η → γχ χ), while J/ψ mostly decay to χs directly (J/ψ → χ χ).In both cases, m χ up to m m /2 is kinematically allowed.The number of produced χs N χ can be calculated by the equation in [8], where c m is the number of mesons produced per proton-on-target (POT), N POT is the total number of POT, = Q/e, and f is a phase space related integral that goes to 0 as m χ approaches m m /2. Figure 1 shows the expected number of χs that reach the detector with face area of 0.5 m 2 , located 280 m from the target.Assuming N POT = 10 22 that corresponds to running the experiment for ∼ 3 years, the expected number of χs is in the order of 10 16 at = 1 and it is in the order of 10 9 at = 10 −4 .
3 Site and detector concept  If χs are produced, because of their low interaction probability, they will penetrate the beam dump facility, the muon monitor, and the sand without a significant energy loss.Therefore, they can be detected at the NM building if a detector sensitive to identifying such particles is installed.We consider the area behind the V-INGRID on B2 (∼ 30 m underground) as a potential detector site.The distance from the axis of the neutrino beam is ∼ 5 m.As described below, the signal acceptance is only slightly smaller than the on-axis location.
The detector concept proposed here is based on a similar proposal made in [9], namely, using a segmented detector with large scintillator bars as the core detector technique.In order to be sensitive to charges down to 10 −3 e, a thick sensitive volume is needed.It is advantageous to segment the large volume because it reduces backgrounds due to dark currents and shower particles from cosmogenic muons to a negligible level.It also allows for utilizing the directionality of the incident χs to further suppress non-pointing particles.The detector, as shown in Figure 3, is composed of 4 layers of stacked 5 × 5 × 80 cm 3 plastic scintillator bars.They are aligned such that the produced χs pass through all layers in a narrow time window.In each layer there are 14 × 14 scintillator bars, so the area of the detector face is about 0.5 m 2 .A prototype of a detector with a similar design has been installed at the LHC, and shown robustness and sensitivity to χs [5].
At the end of each scintillator bar a photodetector is attached to convert the photons to an electronic signal.There are several options, but we consider Photomultipliers (PMTs), because of their large area coverage, low cost, and low dark current.Since the size of a bar is 5 × 5 cm 2 , a fast commercial 50 mm diameter PMT is suitable.Including PMTs, the total volume of the detector is approximately 0.7 × 0.7 × 4 m 3 .The signal acceptance rate, the fraction of χs that go into the detector area of 0.5 m 2 at 280 m from the target, is calculated.Since the photons are massless, the direction of the photons and that of χs should be aligned.Therefore, the direction of photons from π 0 's is used as a proxy for the direction of χs. Figure 4 shows the acceptance rate in 8 × 8 m 2 plane where the axis of the neutrino beam is at the center.The acceptance rate does not depend on the location within this range strongly, i.e., near the corners the rate is only ∼ 10% lower than the central region because of the large distance from the target.This provides some flexibility in selecting the location of the detector.
The operation of SUBMET does not interfere with the operation of the neutrino program and dedicated beam time will not be required.The detector is located behind the V-INGRID detector, so we don't expect any interference with the measurements of the neutrino beam properties and the operation of the INGRID detector.The detector will be triggered by its own trigger system and the experiment will have an independent data stream.Therefore, the operation of this new experiment does not have any conflict with the existing programs at J-PARC.

Background sources
The χs that reach the detector will go through all layers within a ∼ 20 ns time window (coincidence signal).In this section, the background sources that can mimic this coincidence signal are discussed.They can be divided into three categories, random coincidence, beam-induced, and cosmic-induced backgrounds.
In PMTs, spurious current pulses can be produced by thermal electrons liberated from the photocathode.If the rate of such pulses (dark count rate, DCR) is high, the rate of random coincidence can be substantially large even if the time window for coincidence signal is 20 ns.
The random coincidence rate can be calculated by R = nN n τ n−1 where n is the number of layers, N is the DCR, and τ is the coincidence time window.Using a typical DCR of 500 Hz, n = 4, and τ = 20 ns, the random coincidence rate is 6 × 10 −5 per year.Since there can be 14 × 14 = 196 such coincidence signals, the total coincidence rate is ∼ 0.01 per year.
Muons are produced from the pion decays in the decay volume together with neutrinos.The density of quartz, which typically takes up the largest fraction of sand, is 2.65 g/cm 3 and dE/dx = 1.699MeV/cm 2 g [10], so the energy loss of a MIP in > 100 m of sand is much larger than 30 GeV.Therefore, such beam-induced muons can't reach the detector.Although the muons from the pion decays can't reach the detector, neutrinos can and may interact with the scintillator material to produce small signals.The number of neutrino interaction events in INGRID is ∼ 1.5 × 10 8 for N POT = 10 22 [11].Since a large fraction of INGRID material is iron, we can use the rate of INGRID as a upper bound for SUBMET.One layer of SUBMET is approximately 30 times smaller, so the rate is ∼ 5×10 6 for N POT = 10 22 in one layer of SUBMET.
Requiring coincidence in 4 layers, the expected number of this background is negligible.The interaction of the neutrinos and the material of the wall of the NM building in front of the detector can produce muons that go through the detector.These muons can be identified and rejected by installing scintillator plates between the wall and the detector or by using the very large scintillation yield of a muon that can be separated from millicharge signal.
Cosmic muons that penetrate the cavern or the materials above the detector can produce a shower of particles that is large enough to hit all layers simultaneously.In such events, the hits in multiple layers can be within the coincidence time window and will look like a signal event.The particles in the shower generate more photons than χs, so the signals from cosmic muon showers can be rejected by vetoing large pulses.As done to tag the muons produced in the wall of the NM building in front of the detector, scintillator plates can be installed covering the whole detector to tag any ordinary-charge particles or photons incident from top and sides of the detector.These auxiliary components were proven to be effective in rejecting events with such particles [5].In addition, the cosmic shower penetrates the detector sideways, leaving hits in multiple bars in the same layer, while χs will cause a smaller number of hits.A cosmic shower and signals form radioactive decays overlapping with dark current can be another source of background.Since the rate of this background depends on the environment strongly, a precise measurement can be performed in situ only.
To estimate the sensitivity of the experiment, we assume that the total background over three years of running is O(1) event.

Sensitivity
The probability of detecting a χ in a n-layer detector is given by Poisson distribution P = (1 − e −N PE ) n where N PE is the number of photoelectrons.N PE is proportional to the quantum efficiency (QE) of PMT, 2 , and the number of photons that reach the end of the scintillator (N γ ).The 2 term comes from the fact that the energy loss of a charged particle in matter is proportional to 2 .In order to get N γ a GEANT4 [12] simulation is performed.Using a 5 × 5 × 80 cm 3 Saint-Gobain BC-408 plastic scintillator with surface reflectivity of 98%, we get N γ = 6.3 × 10 5 .Using QE = 30% and taking the area of the PMT into account, we get N PE = 1.3 × 10 5 2 .We can now calculate the total number of signal events s = N χ P measured by the detector.Figure 5 shows the number of detected signal events s in the plane of and m χ .Figure 6 shows the 95% CL exclusion curve for N POT = 10 22 .SUBMET provides exclusion down to = 3 × 10 −4 in m χ < 0.2 GeV/c 2 and = 1.5 × 10 −3 in m χ < 1.6 GeV/c 2 .We do not consider a systematic uncertainty on b because it does not have a significant impact on the exclusion limit (Table 1).Sensitivity of SUBMET is comparable with that of FerMINI in m χ < 0.2 GeV/c 2 .Above that mass range the performance of SUBMET is worse because the production rate of J/ψ is much lower with the 30 GeV proton beam.
As shown in Fig 5, the number of signal events recorded by the detector drops rapidly in < 10 −3 due to small N PE .Therefore, one can expect that increasing N PE or N χ would not have a large impact on the sensitivity.This will be discussed in a quantitative way in the next section.

Alternative detector design
In this section, we discuss how sensitivity to sub-millicharged particles depends on the configuration of the detector.If further optimization of the detector is needed, we can use this quantitative study as a guidance.
Silicon Photomultipliers (SiPMs) can be considered instead of PMTs because of its excellent sensitivity to single photons and small size.A challenge is that they have a typical DCR of 5 MHz/cm 2 at room temperature, about three orders of magnitude larger than that of PMTs.This leads to a significantly larger random coincidence rate.Since the dark current is a thermal effect, it has a strong dependence on temperature.Typical reduction rate is a factor of ∼ 2 Table 1: Various detector configurations and their sensitivity.N χ (relative) is the yields of signal events within the acceptance relative to the baseline, N PE (relative) is the number of photon electrons relative to the baseline N PE = 1.3 × 10 5 , and b is the background yields.
is chosen to be 1 × 1 cm due to difficulty in covering the large area of the scintillator.In this case the light attenuation effect increases due to more reflections and the photon yields decreases by a factor of 2.
Searches in the sub-millicharge regime rely on the Poissonian fluctuation of small N PE .
Approximating P N PE to the first order for small , we arrive at the following relation where the subscript = 1 refers to the values at = 1 and s 95% is the number of signal events that provides 95% exclusion limit.Reordering in terms of , we recognize the exclusion limit at Not surprisingly, the -power term induces a sharp cutoff in signal at ∼ O(10 −4 ) as shown in Figure 5.This limits the sensitivity to the same regime regardless of the detector configuration.
Table 1 shows different detector configurations and the corresponding exclusion limit for m χ = 10 MeV/c 2 .The default configuration is in the first row, N χ (relative)=1, N PE (relative)=1, b = 1, and n = 4 where "relative" means relative to the baseline.
The improvements we can achieve by extending the duration of data collection or making the detector area larger (adding more bars to each layer) is very modest.Extending the duration or the detector area by a factor of 2 increases N χ by the same amount.Making an assumption that the background remains same, the sensitivity is only improved by 10% (2 nd row in Table 1).
The impact of b is also small.If b is increased by a factor of 100, the limit is degraded by 15% (3 rd row in Table 1).
As the 4 th and 5 th rows in Table 1 show, the most effective component to enhance sensitivity is N PE .Increasing N PE by a factor of 2(3) improves the sensitivity limit by 25(35)%.This can be achieved by using an inorganic crystal with large Z such as BGO or by making the scintillator bars longer.GEANT4 simulation shows that a 5 × 5 × 20 cm3 BGO bar provides N γ = 1.2 × 10 6 , about a factor of 2 larger than the 5 × 5 × 80 cm 3 BC-408 scintillator.Further increasing the length of the bar to 40 cm, one can achieve N γ = 1.7 × 10 6 .However, cost should be considered as well because BGO is less cost-effective than plastic scintillators.
Though N PE plays the main role in enhancing sensitivity to sub-millicharged particles, the exclusion limit is still in the range of = 1.5 − 3.0 × 10 −4 with the variations considered.This indicates that the sensitivity does not strongly depend on the configuration of the detector.
With 3 layers the number of detected signal events is larger because P = (1 − e −N PE ) 3 instead of (1 − e −N PE ) 4 .However, the background increases as well; the random coincidence rate is approximately 3000 with 3 layers.Using n = 3 and b = 3000, the exclusion limit is slightly improved for m χ < 0.2 GeV/c 2 as shown in Fig 6 .However, the experience from the milliQan demonstrator [5] gives a lesson that the random coincidence and the correlated backgrounds (induced by either cosmic showers or radiation within the scintillator bars) overlapping with dark current are potentially the major background sources.Adding one more layer (4-layer configuration) helps understanding and rejecting these backgrounds so that the experiment is carried out in a low background environment.This will make the observation of a potential signal more robust.

Summary
We studied the feasibility of an experiment searching for millicharged particles using 30 GeV proton fixed-target collisions at J-PARC.With N POT = 10 22 , the experiment provides a sensitivity to χs with charge down to in m χ < 1.6 GeV/c 2 ; this is the regime largely uncovered by the previous experiments.We also explored a few detector designs to achieve an optimal sensitivity to χs.The photoelectron yields is the main driver, but the sensitivity does not have a strong dependence on the detector configurations in the sub-millicharge regime.

Figure 1 :
Figure 1: Expected number of χs that reach a 0.5 m 2 detector located at 280 m from the target.N POT = 10 22 is assumed.
In J-PARC, a 30 GeV proton beam is incident on a graphite target to produce hadrons that subsequently decay to a muon and muon neutrino in the decay volume.The remaining hadrons are then dumped in the beam dump facility.Since they are MIPs, muons penetrate the beam dump and are identified by the muon monitor located behind the beam dump facility.The on-axis near detector (INGRID) is inside the Neutrino Monitor (NM) building located at 280 m from the target.The space between the muon monitor and INGRID is filled with sand.

Figure 2 :
Figure 2: Illustration of the experimental site.χs are produced near the target and reach SUBMET after penetrating the beam dump, the muon monitor and the sand.

Figure 3 :
Figure 3: Demonstration of the SUBMET detector.There are four layers of scintillator bars (blue).At one end of each bar a PMT (black) is attached A χ will penetrate 4 layers in a narrow time window.

Figure 4 :
Figure 4: Signal acceptance rate at 280 m from the target.One box corresponds to 0.5 × 0.5 m 2 .

Figure 5 :
Figure 5: Expected number of detected signal events s in N POT = 10 22 .There is a sharp drop in < 10 −3 due to small N PE .

Figure 6 :
Figure 6: Exclusion at 95% CL for N POT = 10 22 .For comparison, the constraints from previous experiments are shown as shaded area and the expected sensitivity of the proposed experiments are drawn in gray dotted line.Note that FerMINI curve is from 3-layer detector design.

per 10 •
C. So, to suppress the DCR to the similar level of PMT's, the temperature should be lowered by of order of 100 • C. In the scenario where SiPMs are used, the size of scintillator bars N χ (relative) N PE (relative) b Exclusion limit on for m χ = 10 MeV/c 2