Energy calibration of the NEXT-White detector with 1% resolution near Qββ of 136Xe

Excellent energy resolution is one of the primary advantages of electroluminescent high-pressure xenon TPCs. These detectors are promising tools in searching for rare physics events, such as neutrinoless double-beta decay (ββ0ν), which require precise energy measurements. Using the NEXT-White detector, developed by the NEXT (Neutrino Experiment with a Xenon TPC) collaboration, we show for the first time that an energy resolution of 1% FWHM can be achieved at 2.6 MeV, establishing the present technology as the one with the best energy resolution of all xenon detectors for ββ0ν searches.


Introduction
Searches for neutrinoless double beta decay (ββ0ν), the observation of which would imply total lepton number violation and would show that neutrinos are Majorana particles [1][2][3][4], require excellent energy resolution to eliminate background events that occur at energies similar to the Q-value of the decay (Q ββ ). The NEXT (Neutrino Experiment with a Xenon TPC) collaboration [5][6][7][8] intends to search for ββ0ν using of order 100 kg of xenon enriched to 90% in the candidate isotope 136 Xe (Q ββ = 2457.8 keV). In recent years, NEXT has developed and operated several gaseous xenon TPCs, including kg-scale detectors at Lawrence Berkeley National Lab (LBNL) and Instituto de Física Corpuscular (IFIC) [9,10] and more recently the 5 kg-scale NEXT-White 1 at the Canfranc Underground Laboratory (LSC) in the Spanish Pyrenees [11].
Previous analyses [12] of the NEXT-White energy resolution using gamma rays from 137 Cs and 232 Th sources showed an extrapolated 1% FWHM resolution at Q ββ . The relatively low pressure (7.2 bar) at which those data were taken meant that electron tracks of events with energy near Q ββ were not easily contained in the detector. Low statistics at the photopeak limited the highest energy at which a detailed analysis of energy resolution was performed to 1.6 MeV. More data has since been taken at a higher pressure (10.3 bar), and the results are reported in the present study. The experimental setup, similar to that of the previous study [12], is reviewed in section 2, and the analysis and obtained energy resolution are presented in section 3.

The NEXT-White electroluminescent TPC
The experimental setup is similar to that of the preceding study [12] and is summarized here. The detector NEXT-White is an electroluminescent (EL) time projection chamber  (TPC) filled with xenon gas and equipped with photosensors to detect the UV light emitted in interactions occurring within the active volume. Charged particles deposit energy within the drift region, producing a track of ionized and excited xenon atoms. The UV light emitted in the relaxation of the excited xenon atoms, called primary scintillation or S1, is detected immediately and the ionized electrons are drifted toward a readout plane consisting of a narrow region of high electric field, the EL gap. In passing through the EL gap, the electrons are accelerated to energies high enough to further excite, but not ionize, the atoms of the xenon gas, leading to the production of an amount of secondary scintillation photons (S2) proportional to the number of electrons traversing the gap. This amplification process, electroluminescence, allows for gains on the order of 1000 photons per electron with significantly lower fluctuations than avalanche gain. In addition, the time elapsed between the observation of S1 and the arrival of S2 can be used to determine the axial (z) coordinate at which the interaction took place. In NEXT-White (see figure 1 and also [11]), the primary (S1) and secondary (S2) scintillation are detected by an array of 12 Hamamatsu R11410-10 photomultiplier tubes (PMTs), called the "energy plane" placed 130 mm from a transparent wire mesh cathode held at negative high voltage. An electric field is established in the drift region defined by the cathode and another transparent mesh (the gate) located about 53 cm away. The EL region is defined by the mesh and a grounded quartz plate coated with indium tin oxide (ITO), placed 6 mm behind it. A grid (10 mm pitch) of 1792 SensL series-C silicon photomultipliers (SiPMs) is located behind the EL gap and measures the S2 scintillation, providing precise information on where the EL light was produced in (x, y). The active volume is shielded by an 60 mm thick ultra-pure inner copper shell, and the sensor planes are mounted on pure copper plates of thickness 120 mm. The sensor planes and active volume are enclosed in a pressure vessel constructed from the titanium-stabilized stainless steel alloy 316Ti. The vessel sits on top of a seismic table; and a lead shield that can be mechanically opened and closed surrounds the vessel. The vessel is connected to a gas system through which the xenon gas is continuously purified via the use of a hot getter. The entire experimental area, including gas system, electronics, pressure vessel, and seismic table, are stationed on an elevated tramex platform in the Laboratorio Subterráneo de Canfranc (LSC) in the Spanish Pyrenees.

Run configuration
As the goal of the present analysis was a detailed study of energy resolution, calibration sources were employed to yield energy peaks over a range of energies from several tens of keV up to and including Q ββ . 83m Kr was injected into the xenon gas, providing a uniform distribution of 41.5 keV point-like energy depositions used to map out the geometric variations in the sensor responses and electron lifetime of the detector [13]. 137 Cs and 228 Th calibration sources were also placed as shown in figure 1. The 137 Cs source provided 661.6 keV gamma rays, and 228 Th decays to 208 Tl which provides gamma rays of energy 2614.5 keV. In this study we focus on the energy peaks produced by interactions of these 137 Cs and 208 Tl gamma rays, and also the double-escape peak resulting from e + e − pair production interactions of the 208 Tl gamma ray in which the two 511 keV gamma rays escape. For the present analysis, the acquisition trigger was split into a lower-energy trigger seeking the 83m Kr events and a high-energy trigger aimed at capturing events with energy above approximately 150 keV. A summary of the datasets analyzed is given in table 1.
For each run, the low-energy triggers were used to compute the lifetime and geometric correction maps used to correct the events acquired with the high-energy trigger. The average electron lifetime determined over the course of the analyzed runs is also shown in Lifetime ( s) Figure 2. The average electron lifetime over the course of the analyzed runs determined using 83m Kr events. The errors on the measurements are smaller than the size of the points. Though consistently above 2 ms, the electron lifetime was unstable and therefore was monitored and the corresponding corrections determined for each run. outgassing of internal components, and for this reason corrections for each run needed to be determined using data taken over the same period of time.

Energy resolution
The signals from the SiPMs and PMTs were digitized in samples of width 1 µs and 25 ns respectively. Individual pulses in the energy plane waveform (summed over all PMTs, see figure 3) were selected and classified as either S1 or S2. Events with a single identified S1 were selected, and the S2 peaks were divided into "slices" of width 2 µs.
The pattern of light detected by the SiPMs of the tracking plane during the 2 µs interval of the slice was used to reconstruct the (x, y) location of the EL production, as done in [13], except multiple reconstructed positions sharing the energy E of a single slice were possible in order to allow for reconstruction of tracks that extend over greater distances (several cm or more) in the (x, y) plane within a single slice. The time elapsed since the S1 pulse was used to determine the z coordinate of each slice, and the energies E of the reconstructed depositions (x, y, z, E) were then multiplied by two correction factors: one accounting for the geometrical (x, y) dependence of the light collection over the EL plane, and another  accounting for losses due to the finite electron lifetime caused by attachment to impurities. This second factor depended on the drift length (z-coordinate) and the location in the EL plane (x, y), as the electron lifetime also varied in (x, y). Once fully reconstructed, fiducial cuts were made on each event as detailed in figure 4.
A final correction was applied for an empirically observed dependence on the track orientation. The origin of this dependence, whereby the measured energy of an event decreases with increasing axial (z) extent of the track, is still under study. However, we find it can be effectively corrected, as follows. The z-extent ∆z is defined as the difference between the maximum and minimum z-coordinates of all reconstructed slices in the event. The effect is shown in figure 5 along with the resolution obtained for each of the three peaks (662 keV, 1592 keV, and 2615 keV) after correcting for the effect using the average of the normalized slopes determined by a linear fit to each distribution, used to determine (m/b), and the remaining data, runs 6352, 6365, 6482, 6483, 6484, and 6485, were used to evaluate the energy resolution (see table 1).
Each peak was fitted to the sum of a Gaussian and a 2nd-order polynomial to account for the surrounding distribution of background events, and the resolution was computed using the width of the Gaussian. The obtained resolutions are: 1.20 ± 0.02% FWHM at 662 keV; 0.98 ± 0.03% FWHM at 1592 keV; and 0.91 ± 0.12% FWHM at 2615 keV. The total errors are estimated in each case based on the statistical errors of the fits (shown on the histograms in figure 5) and systematic effects including variations in the range of events included in the fit and (systematic) errors in the correction for the axial length effect. The energy conversion from detected photoelectrons to keV was determined (after application of all corrections) using a quadratic fit to the means of the three peaks of interest (662 keV, 1592 keV, 2615 keV) and the 29.7 keV K-α xenon X-ray peak. The X-rays had energies too low to be triggered on as individual events, but their energies were visible by examining the spectrum of isolated energy depositions within all events, which included small depositions due to xenon X-rays that traveled away from the main track before interacting.
The energy spectrum of high-energy triggers in the full active volume is shown in figure 6 after applying all corrections described in the present section. Unlike in the previous study [12], the 208 Tl photopeak at 2615 keV (near Q ββ ) is clearly resolved. The squared resolution is shown vs. the inverse energy in figure 7 for the three energy peaks studied and fit to a line, R 2 = a/E + c, where a = 548.52 ± 82.15 %FWHM 2 · keV and c = 0.62 ± 0.10 %FWHM 2 . The presence of a constant term in the resolution shows that detector-specific systematic effects have not been completely eliminated in our analysis, and there is still room for further improvement. Nevertheless, these results demonstrate  that excellent energy resolution is obtainable throughout the entire fiducial volume once correction for the axial length effect is made. Unlike corrections for electron lifetime, the correction for the axial length effect was relatively stable throughout the time (about 8 weeks) over which the data presented in this study was taken. This effect is thought to be a result of some internal nonlinearity in light production, possibly due to alterations of the EL field during track readout in a "charging-up" effect as electrons cross the EL gap, or due to electron loss in the EL gap from attachment to impurities produced by photoionizing the wavelength shifting material deposited on the surface behind the EL region. Other potential explanations such as PMT saturation, electron recombination with the ions of the original track, variations in electron lifetime throughout the detector, and emission of light by the SiPMs, have already been discredited. Further details are given in appendix A.

Summary
Energy resolution in the NEXT-White TPC has been further studied, and a resolution better than 1% FWHM is shown to be obtainable at 2615 keV, as predicted in the preceding study [12]. This resolution was obtained over nearly the entire active volume, demonstrating the effectiveness of the continuous 83m Kr-based calibration procedure implemented to correct for geometric and lifetime effects, and improved slightly with more restrictive fiducial cuts. Further study is required to understand the observed "axial length effect" in which the measured energy of extended tracks decreases with increasing track length in the axial (drift) direction. However, as high-pressure xenon TPCs provide detailed energy and topological information for each event, such effects can be remedied through careful calibration, and the outstanding resolution obtained highlights the strong potential of this detector technology to host a sensitive 0νββ search in which good energy resolution is essential.

JHEP10(2019)230
required to drift the electrons over the entire track length in z. For 208 Tl photopeak events (see figure 5, bottom), this would be about 120 µs, and integrating over this interval after the arrival of S1 for many such events, no evidence of the expected light was observed.
• Variations in electron lifetime: as the measured electron lifetime in NEXT-White is known to vary with location in the detector, there has been concern that small errors in the computation of the lifetime were giving rise to the observed effect when applied over long tracks. However, even after correcting Cs-photopeak events using a single average position (assuming pointlike tracks), the effect could still be observed by making a tight cut on average radius (effectively eliminating the error due to response variations in the xy-plane by considering only events that did not require significant xy correction).
• Light emitted from the SiPMs: the effect is also seen in the integrated charge of the SiPMs, and in fact is more dramatic (the normalized slopes m/b analogous to those shown in figure 5 are greater in magnitude). Therefore it was proposed that the SiPMs may be emitting additional light in a nonlinear manner during the production of EL. However, even after turning off the SiPM plane and using only information from the PMT plane for a less-precise xy reconstruction, the effect was still observed.
Several explanations for the effect have not yet been investigated in detail: • "Charge-up" effect at the EL plane: an electron crossing the EL gap may, at least locally, alter the electric field seen by the next electron crossing the gap for some amount of time. If this were to make the average gain somewhat dependent on track orientation -whether the electrons cross the gap more in "series" (more extended in z) or in "parallel" (more extended in xy) -this could give rise to the observed effect.
• Attachment to ionized impurities in the EL gap: the wavelength shifter tetraphenyl butadiene (TPB) is deposited on several components in NEXT-White including the quartz plate just behind the EL region, to shift the VUV scintillation produced by xenon to visible light that can be detected by the photosensors (the SiPMs are not VUV sensitive, and the PMTs are placed inside enclosures behind sapphire windows, which do not transmit VUV light, to shield them from the high pressure environment inside the detector). If the photons produced in electroluminescence are capable of photoionizing the TPB, the resulting ions would be drifted across the EL region, possibly capturing some of the electrons that arrived at later times before completely traversing the EL gap and thereby reducing the observed energy of the event.
The observed effect could also be a result of a nonlinearity in the light production process caused by some other internal component. Further investigation in future runs with NEXT-White, possibly involving alterations of the internal hardware and/or running systematically at different EL gains, will be necessary to understand this effect. Nevertheless, excellent resolution has been obtained due to the properties of high-pressure xenon TPCs, such as simultaneous energy and 3D position measurements, which allow for detailed calibration.

JHEP10(2019)230
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.