Gamma-induced background in the KATRIN main spectrometer

TL;DR

This study investigates whether gamma radiation from external sources significantly contributes to background noise in the KATRIN experiment's main spectrometer, finding that shielding effectively reduces such background events.

Contribution

The paper provides the first detailed simulation and experimental analysis quantifying gamma-induced background in the KATRIN main spectrometer, establishing an upper limit on this background source.

Findings

Gamma radiation contributes less than 0.006 counts/sec to background.

Electrostatic and magnetic shielding effectively suppress gamma-induced secondary electrons.

Gamma-induced background is not a significant limiting factor for KATRIN's sensitivity.

Abstract

The KATRIN experiment aims to measure the effective electron antineutrino mass $m_{\overline{\nu}_e}$ with a sensitivity of 0.2 eV/c$^2$ using a gaseous tritium source combined with the MAC-E filter technique. A low background rate is crucial to achieving the proposed sensitivity, and dedicated measurements have been performed to study possible sources of background electrons. In this work, we test the hypothesis that gamma radiation from external radioactive sources significantly increases the rate of background events created in the main spectrometer (MS) and observed in the focal-plane detector. Using detailed simulations of the gamma flux in the experimental hall, combined with a series of experimental tests that artificially increased or decreased the local gamma flux to the MS, we set an upper limit of 0.006 count/s (90% C.L.) from this mechanism. Our results indicate the effectiveness of the electrostatic and magnetic shielding used to block secondary electrons emitted from the inner surface of the MS.

Tables

Table 1. Table 1 : Specific activities (in units of Bq / kg Bq kg \mathrm{Bq}\mathrm{/}\mathrm{k}\mathrm{g} ) determined from radioassay measurements for materials in the KATRIN spectrometer hall. The activity of the wall concrete was only measured for the basement walls; for the upper walls, the same distribution of activities was assumed. The steel activities were measured at the Oroville Low Background Facility  Smith2015 .

Isotope	Wall concrete	Floor concrete	MS steel
40K	$409\pm 22$	$61\pm 3$
238U	$23\pm 5$	$6.0\pm 0.7$	$\left(5.0\pm 2.5\right)\text{×}{10}^{-2}$
232Th	$18\pm 1$	$3.3\pm 0.3$
235U	$1.7\pm 0.5$	$0.6\pm 0.1$
60Co			$\left(1.8\pm 0.3\right)\text{×}{10}^{-3}$

Table 2. Table 2 : Simulated total gamma fluxes integrated at the inner surface of the MS. For each component, the flux was calculated based on the activities listed in Table 1 , after summing the individual contributions from the different radioisotopes. The errors are purely systematic and indicate the uncertainty of the component activities.

Component	MS flux (gammas/s)	Percentage
Concrete walls	$\left(2.06\pm 0.08\right)\text{×}{10}^6$	80.2 %
Concrete floors	$\left(4.1\pm 0.2\right)\text{×}{10}^5$	15.8 %
Air	$\left(9.0\pm 2.7\right)\text{×}{10}^4$	3.5 %
Steel hull	$\left(1.1\pm 0.5\right)\text{×}{10}^4$	0.4 %
Total	$\left(2.56\pm 0.09\right)\text{×}{10}^6$

Table 3. Table 3 : The voltage settings used during the gamma-induced background measurements.

Component	Voltage	Value/range ($V$)
MS hull	$U_{\text{0}}$	$-\mathrm{18\hspace{0.17em}500}\text{ to }-\mathrm{18\hspace{0.17em}600}$
MS inner electrode	$\Delta U_{\text{IE}}$	$0\text{ to }-100$
Post-acceleration electrode	$U_{\text{PAE}}$	$\mathrm{10\hspace{0.17em}000}$
Bias ring of detector wafer	$U_{\text{bias}}$	$120$

Table 4. Table 4 : Results from the FPD measurements with the 60 Co source, where the rates are measured in units of counts per second (cps). Rate is the raw result from the detected counts and measurement duration Δ ​ t Δ 𝑡 \Updelta t , while Rate ∗ is corrected for the broken preamplifier module. The errors on the rates are statistical. The magnetic field settings used are those shown in Fig.  4 . For the Asym. M measurements, all rates exclude the inner 16 detector pixels.

B-field	$U_{\text{0}}$ (kV)	$\Delta U_{\text{IE}}$ (V)	Source	Counts	$\Delta t$ (s)	Rate ($\mathrm{cps}$)	$\Delta$Rate ($\mathrm{cps}$)	Rate∗ ($\mathrm{cps}$)	$\Delta$Rate∗ ($\mathrm{cps}$)
Asym. M	$-18.5$	$0$	open	1500629	1790	$838.3\pm 0.7$	$222.4\pm 0.9$	$871.9\pm 0.7$	$231.7\pm 0.9$
			closed	1102527	1790	$615.9\pm 0.6$		$640.2\pm 0.6$
Asym. M	$-18.5$	$-100$	open	164532	1790	$91.9\pm 0.2$	$25.8\pm 0.3$	$95.6\pm 0.2$	$26.8\pm 0.3$
			closed	118376	1790	$66.1\pm 0.2$		$68.8\pm 0.2$
Sym.	$-18.5$	$-100$	open	18254	33500	$0.545\pm 0.004$	$0.005\pm 0.006$	$0.566\pm 0.004$	$0.005\pm 0.006$
			closed	14761	27320	$0.540\pm 0.004$		$0.561\pm 0.005$

Table 5. Table 5 : Total gamma fluxes through the interior of the MS as determined from simulations. Φ 0 subscript Φ 0 \Phi_{0} is either the flux with the closed 60 Co source (second column) or with no water shielding (third column). Similarly, Φ 1 subscript Φ 1 \Phi_{1} is either the flux with the open 60 Co source or with water shielding. The flux includes all crossings (ingoing and outgoing) of the inner surface. The errors include both statistical and systematic uncertainties, although the latter only includes the uncertainty on the activity of the gamma sources (see Table  1 ). Other systematic effects, such as the accuracy of the included Geant4 physics processes and the correctness of the simulation geometry, were not calculated.

Flux (${10}^6\mathrm{gammas}/s$)	60Co	Water
${\Phi }_0$	$2.8\pm 0.1$	$2.564\pm 0.090$
${\Phi }_1$	$20.9\pm 0.9$	$2.411\pm 0.088$
${\Phi }_1-{\Phi }_0$	$18.2\pm 0.9$	$-0.153\pm 0.002$

Table 6. Table 6 : Results from the FPD measurements with(out) the water shielding. “Asym. U” and “Asym. D” indicate that the asymmetric magnetic field images the upstream ( z = 𝑧 absent z= − 3.9   m  to  − 1.0   m range times -3.9 meter times -1.0 meter -3.9\text{\,}\mathrm{m}-1.0\text{\,}\mathrm{m} ) and downstream ( z = 𝑧 absent z= 1.3   m  to  3.7   m range times 1.3 meter times 3.7 meter 1.3\text{\,}\mathrm{m}3.7\text{\,}\mathrm{m} ) region of the MS surface, respectively (see Fig. 4 ). Rate is the raw result from the detected counts and measurement duration Δ ​ t Δ 𝑡 \Updelta t , while Rate ∗ is corrected for the missing detector pixels. For the Asym. U (Asym. D) measurements, all rates exclude the inner 28 (4) pixels. The errors on the rates are statistical.

B-field	$U_{\text{0}}$ (kV)	$\Delta U_{\text{IE}}$ (V)	Shielding	Counts	$\Delta t$ (s)	Rate ($\mathrm{cps}$)	$\Delta$Rate ($\mathrm{cps}$)	Rate∗ ($\mathrm{cps}$)	$\Delta$Rate∗ ($\mathrm{cps}$)
Asym. U	$-18.6$	$0$	no water	21319281	35900	$593.9\pm 0.1$	$-2.2\pm 0.2$	$620.1\pm 0.1$	$-2.3\pm 0.2$
			water	21241447	35900	$591.7\pm 0.1$		$617.8\pm 0.1$
Asym. D	$-18.6$	$0$	no water	13767972	21990	$626.1\pm 0.2$	$-2.6\pm 0.2$	$653.7\pm 0.2$	$-2.8\pm 0.2$
			water	24620183	39490	$623.5\pm 0.1$		$650.9\pm 0.1$
Sym.	$-18.5$	$-100$	no water	19411	43080	$0.451\pm 0.003$	$0.0008\pm 0.0043$	$0.469\pm 0.003$	$0.0008\pm 0.0045$
	$-18.6$	$-100$	water	25927	57440	$0.451\pm 0.003$		$0.469\pm 0.003$

Table 7. Table 7 : Values used in the calculation of the secondary electron yield ( Y 𝑌 Y ) and fraction of secondaries induced by environmental gamma radiation ( f env subscript 𝑓 env f_{\text{env}} ). See the text for details about each listed parameter. Φ env subscript Φ env \Phi_{\text{env}} and Δ ​ Φ Δ Φ \Delta\Phi differ for each asymmetric field setting (since different regions of the spectrometer surface are measured for each setting), and these values also differ from the values given in Table 5 (since the latter considered the entire spectrometer surface).

	60Co source	60Co source	Water shielding	Water shielding
Magnetic field	Asym. M	Asym. M	Asym. U	Asym. D
z-range (m)	$-2.44\text{ to }1.58$	$-2.44\text{ to }1.58$	$-3.88\text{ to }-0.99$	$1.32\text{ to }3.70$
$\Delta U_{\text{IE}}$ (V)	$0$	$-100$	$0$	$0$
$R$ (cps)	$640.2\pm 0.6$	$68.8\pm 0.2$	$620.1\pm 0.1$	$653.7\pm 0.2$
$\Delta R$ (cps)	$231.7\pm 0.9$	$26.8\pm 0.3$	$-2.3\pm 0.2$	$-2.8\pm 0.2$
${\Phi }_{\text{env}}$ (${10}^5$ gammas/s)	$4.5\pm 0.2$	$4.5\pm 0.2$	$3.3\pm 0.1$	$2.7\pm 0.1$
$\Delta \Phi$ (${10}^5$ gammas/s)	$86\pm 4$	$86\pm 4$	$-0.257\pm 0.004$	$-0.209\pm 0.003$
$P_{\text{arrival}}$ (%)	$8.8\pm 0.1$	–	$11.2\pm 0.1$	$16.6\pm 0.1$
$Y$ (${10}^{-4}$ e-/gamma)	$3.1\pm 0.2$	–	$7.9\pm 0.7$	$8.0\pm 0.6$
$f_{\text{env}}\left({10}^{-2}\right)$	$1.9\pm 0.1$	$2.1\pm 0.1$	$4.6\pm 0.4$	$5.5\pm 0.5$

Table 8. Table 8 : Background rate R env subscript 𝑅 env R_{\text{env}} induced by environmental gamma radiation under standard conditions (symmetric magnetic field and Δ ​ U IE = − 100   V Δ subscript 𝑈 IE times -100 volt \Updelta U_{\textrm{IE}}=$-100\text{\,}\mathrm{V}$ ). The relevant values used to calculate this rate are also listed; the values of Δ ​ Φ Δ Φ \Delta\Phi and Φ env subscript Φ env \Phi_{\text{env}} come from Table 5 . A corrective factor of 2.6 was applied to the upper limit on R env subscript 𝑅 env R_{\text{env}} for the 60 Co source.

	60Co source	Water shielding
$\Delta \Phi$ (${10}^5\gamma /s$)	$182\pm 9$	$-1.53\pm 0.02$
${\Phi }_{\text{env}}$ (${10}^5\gamma /s$)	$25.6\pm 0.9$	$25.6\pm 0.9$
$\Delta R$ (mcps)	$5\pm 6$	$0.8\pm 4.5$
$R_{\text{env}}$ (mcps)	$0.7\pm 0.9$	$-13\pm 75$
Upper limit ($90\%C.L.$) on $R_{\text{env}}$ (mcps)	$5.6$	$110$

Equations

$$E_{\text{expected}}=E_{\text{0}}+e[U_{\text{PAE}}+U_{\text{bias}}-U_{\text{0}}-\text{Updelta}{U}_{\text{IE}}],$$

$$R$$

$$S$$

$$\Phi$$

$$\Delta R_{\text{cobalt}}=R_{\text{open}}-R_{\text{closed}},$$

$$\Delta R_{\text{shielding}}=R_{\text{water}}-R_{\text{no water}}.$$

$$Y=\frac{\Delta S}{\Delta \Phi }.$$

$$\Delta S=\frac{\Delta R}{ϵ\cdot P_{\text{arrival}}},$$

$$F(E)\propto \frac{E}{(E+W{)}^4},$$

$$f_{\text{env}}=\frac{S_{\text{env}}}{S}\approx \frac{R_{\text{env}}}{R}.$$

$$R_{\text{env}}\approx \frac{\Delta R}{\Delta \Phi }{\Phi }_{\text{env}}.$$

$$f_{\text{env}}\approx \frac{{\Phi }_{\text{env}}}{\Delta \Phi }\frac{\Delta R}{R}.$$

$$R_{\text{env}}\leq$5.6\text{\,}\mathrm{mcps}$~{}($90\text{\,}\mathrm{\char 37\relax}\text{\,}\mathrm{C.L.}$)$$

Full text

††journal: Journal

1 Introduction

The Karlsruhe Tritium Neutrino (KATRIN) experiment has been designed to reach a sensitivity of $0.2\text{\,}\mathrm{eV}\mathrm{/}\mathrm{c}^{2}$ (90% C.L.) on the effective electron antineutrino mass $m_{\overline{\nu}_{e}}$ KATRIN2005 . As the successor to the Mainz Kraus2005 and Troitsk Aseev2011 experiments, KATRIN will precisely measure the energy of electrons produced from tritium $\upbeta$-decay and determine $m_{\overline{\nu}_{e}}$ by fitting the shape of the $\upbeta$-spectrum near the endpoint energy ($18.6\text{\,}\mathrm{keV}$). Reaching the sensitivity goal requires a thorough understanding and mitigation of background sources along the entire beamline of the experiment. Background electrons produced in the main spectrometer (MS), the high-resolution MAC-E filter Beamson1980 ; Lobashev1985 ; Picard1992 that analyzes the energy of the $\upbeta$-particles, are of particular importance. Several sources of background in the MS have already been analyzed in detail, including cosmic-ray muons Altenmueller2019 and radon decays RadonPaper .

X-rays and gammas can also contribute to the background and have been previously studied with other MAC-E filter spectrometers. The Mainz spectrometer PhDFlatt2004 and the KATRIN pre-spectrometer DthLammers2009 were both irradiated using an X-ray tube with a peak energy of $E_{\text{X-ray}}=$70\text{,}\mathrm{keV}$$. These tests showed the effectiveness of the electrostatic and magnetic shielding measures in place, which suppressed the background contribution caused by the X-ray tube by up to a factor of 50. However, the total background was still strongly elevated while the irradiation took place, indicating that gamma radiation can potentially have a large contribution to the background rate for MAC-E filter spectrometers. In this paper, the effect of environmental gamma radiation on the MS background rate is investigated through a combination of simulation and measurement data, and the effectiveness of the shielding in the MS is demonstrated.

The experimental setup of the KATRIN experiment is described in Section 2, focusing on the MS and the focal-plane detector (FPD) system, which measures the electron rate at the exit side of the MS. Section 3 discusses the background generation due to environmental gammas and gives details about their simulation in the spectrometer hall. In Section 4, the simulation results are compared to background measurements under conditions of gamma enhancement and suppression. The contribution of environmental gamma radiation to the KATRIN background rate is derived in Section 5, and some final remarks on the MS background are provided in Section 6.

2 Experimental apparatus

The KATRIN experiment is located at the Karlsruhe Institute of Technology (KIT), Campus North, near Karlsruhe, Germany. The approximately $70$-$\mathrm{m}$ long beamline is shown in Fig. 1. The windowless gaseous tritium source has been designed to produce more than ${10}^{11}$ $\upbeta$-particles per second Babutzka2012 . Superconducting magnets adiabatically guide emitted electrons through the differential PhDKosmider2012 and cryogenic Gil2010 ; PhDJansen2015 pumping sections, where the tritium flow must be reduced by 14 orders of magnitude KATRIN2005 .

Two MAC-E filters are available to analyze the electron energy. The smaller pre-spectrometer can be used as a pre-filter to reflect low-energy $\upbeta$-particles back toward the source and only allow the highest energy electrons, close to the endpoint energy, to enter the MS. The precise energy filtering is performed in the MS, which operates with an energy resolution of about $1\text{\,}\mathrm{eV}$ near the tritium endpoint KATRIN2005 . Electrons that pass through the MS are measured with the FPD system. Due to their importance in studying the gamma-induced background, the MS and FPD system will be discussed in detail in the following subsections.

2.1 Main spectrometer

The largest component of the KATRIN beamline is the MS (see Fig. 1). The stainless steel vessel has an inner diameter of $9.8\text{\,}\mathrm{m}$ at its center and a total length of $23.3\text{\,}\mathrm{m}$ KATRIN2005 . The walls of the MS vary in thickness between $32\text{\,}\mathrm{mm}$ for the central cylindrical region and $25\text{\,}\mathrm{mm}$ for the conical sections KATRIN2005 . During standard KATRIN operation, the retarding voltage $U_{\textrm{0}}$ applied to the MS hull will be varied systematically around $-18.6\text{\,}\mathrm{kV}$ in order to analyze electrons with energies close to the tritium endpoint energy. To prevent scattering with residual gas, the MS is designed to operate under ultra-high vacuum conditions, with a pressure close to ${10}^{-11}\text{\,}\mathrm{mbar}$ Arenz2016 .

Superconducting solenoids placed at the entrance and exit of the MS generate the magnetic flux tube that is responsible for adiabatically guiding electrons through the vessel Arenz2018a . Air-coils surrounding the MS are used to fine-tune the magnetic field inside the vessel and compensate for the Earth’s magnetic field Glueck2013 ; Erhard2017 . By manipulating the magnitude and polarity of the air-coil currents, it is possible to create non-standard magnetic field settings inside the vessel (see Section 4).

The magnetic field inside the MS acts as a passive shield against charged particles emitted from the vessel surface; the particles are deflected by the Lorentz force back toward the walls of the vessel, or they follow magnetic field lines that do not reach the detector. Additional shielding is provided by a wire electrode system installed at the inner surface of the MS Valerius2010 . The inner electrodes (IE) can be placed at an offset potential $\Delta U_{\textrm{IE}}$ relative to the vessel (i.e. $U_{\textrm{IE}}=U_{0}+\Updelta U_{\textrm{IE}}$). If $\Updelta U_{\textrm{IE}}$ is set to a negative value, most low-energy electrons emitted from the vessel walls will be repelled by the wires back toward the surface.

2.2 Focal-plane detector

The FPD Amsbaugh2015 is a monolithic, segmented silicon PIN diode that detects particles emerging from the MS. These particles pass through an ultra-high vacuum system mated to the exit of the MS. Charged particles can be further accelerated using a post-acceleration electrode (PAE), with voltage $U_{\textrm{PAE}}$, immediately preceding the detector wafer; this acceleration helps distinguish signal electrons from background originating within the FPD system. A superconducting solenoid focuses charged particles onto the wafer, which has a sensitive area with a diameter of $90\text{\,}\mathrm{mm}$ Amsbaugh2015 . The dartboard-segmentation pattern of the wafer into 148 equal-area pixels provides sensitivity to the transverse spatial distribution of the particles. The wafer is divided into 12 concentric rings, each with 12 pixels, and a central bullseye with four pixels. The PIN diode is biased with a voltage $U_{\textrm{bias}}$; signals pass through preamplifiers (located in vacuum) before readout.

Energy and timing for each event are analyzed online in the data-acquisition system using a cascaded pair of trapezoidal filters. Energy calibration is provided by gammas from an 241Am calibration source. The electron response is characterized using a photo-electron source with adjustable energy Amsbaugh2015 . Each source can be inserted into the line of sight of the FPD for a periodic, dedicated calibration run. An energy resolution of $1.52\text{\,}\mathrm{keV}$ (full width at half-maximum, FWHM) and timing resolution of $246\text{\,}\mathrm{ns}$ (FWHM) have been achieved with the system using $18.6\text{\,}\mathrm{keV}$ electrons and a $6.4\text{\,}\mathrm{\SIUnitSymbolMicro s}$ shaping time Amsbaugh2015 .

3 Environmental gamma radiation in the spectrometer hall

During standard KATRIN operation, $\upbeta$-particles that are detected by the FPD must overcome the retarding potential applied to the MS. This potential reaches its largest value near the middle of the vessel, approximately equidistant from both ends of the MS. Because the value of the retarding potential will be scanned close to the endpoint energy, $\upbeta$-particles traveling through the middle region of the MS will have low kinetic energies, below about $30\text{\,}\mathrm{eV}$ KATRIN2005 . Because of the relatively poor energy resolution of the FPD compared to the MS, any low-energy secondary electrons produced in this region cannot be energetically distinguished from the signal $\upbeta$-particles.

When passing through the MS steel, environmental gamma radiation (i.e., primary radiation) can cause the emission of secondary particles, including secondary electrons. Of particular concern is the generation of “true-secondary” electrons, which are defined to have energies below $50\text{\,}\mathrm{eV}$ Furman2002 . The inner surface of the MS provides a large area for electron emission: $690\text{\,}\mathrm{m}^{2}$ for the steel hull and $532\text{\,}\mathrm{m}^{2}$ for the wire electrode system (although the effective surface area for emission from the latter is reduced due to the two-layer structure of the IE) Arenz2016 . Due to imperfections in the magnetic and electrostatic shielding, secondary electrons emitted from these surfaces may have a small probability to enter the sensitive magnetic flux tube that connects to the FPD. The validity of this background-generating mechanism has been confirmed with other MAC-E-filter spectrometers PhDFlatt2004 ; DthLammers2009 . Because of the increased size of the spectrometer, there is the potential for a significant background contribution from gamma-induced surface electrons in the MS.

3.1 Radioactivity in the spectrometer hall

The MS was constructed using low-radioactivity materials in order to limit background production from environmental gammas. The central cylindrical portion of the MS vessel was built from $32\text{\,}\mathrm{mm}$-thick sheet metal, composed of type 316LN stainless steel. A sample of this metal, from the same batch used for the MS vessel, was measured for radioisotopes; the results of this study are given in Table 1. However, the primary source of gammas in the KATRIN spectrometer hall is not the vessel itself but rather the intrinsic radioactivity of the concrete used to construct the walls and floors.

The walls of the spectrometer hall are made of standard concrete. To reduce the effect of ambient gamma radiation, the floors were built using a low-activity concrete due to their close proximity to the MS. During construction of the hall, samples from every load of low-activity concrete were monitored with a NaI detector to ensure the activity fell within specifications. Additionally, the gamma spectra of several samples of the concrete were measured at KIT using a shielded HPGe detector. The concrete activities shown in Table 1 were derived from these measurements. 40K was found to be the largest contributor to the activity in the concrete, followed by the 238U and 232Th decay chains.

3.2 Simulation of the gamma flux

To better understand the background due to environmental gamma radiation, simulations of the gamma flux in the MS were performed with the Geant4 simulation toolkit geant:2003 ; geant:2006 ; geant:2016 , version 10.4.p02. A simplified reproduction of the KATRIN spectrometer hall was implemented; the geometry included the steel MS vessel and the concrete walls and floors (Fig. 2). Radioactive isotopes were uniformly created in the walls, floors, and steel with the relative production rates determined from the radioassay measurements given in Table 1. The effect of 222Rn in the air was also included, assuming an activity of $49\pm 15\text{\,}\mathrm{Bq}\mathrm{/}\mathrm{m}^{3}$, which is the average indoor radon level in Germany WHO2007 . The decays of the isotopes and any subsequent daughters were handled by the Geant4 Radioactive Decay Module, and secular equilibrium was assumed. This module has recently been updated by the Geant4 collaboration to include newer versions of the Evaluated Nuclear Structure Data File (ENSDF) datasets ENSDF and to better ensure energy conservation for decays geant:2016 . Available physics processes were set by the Shielding_EMZ physics list, which uses the most accurate electromagnetic physics models and is well-suited for shielding simulations geant:PhysicsList .

The gamma fluxes determined from simulation are listed in Table 2 for each geometrical component. As expected, the concrete walls are the primary contributor to the gamma flux inside the MS, followed by the concrete floors. The gammas originating from the steel vessel and the air, combined together, make up less than $4\text{\,}\%$ of the flux inside the MS.

Gamma spectra were measured at various locations in the spectrometer hall using a HPGe detector. One of these spectra, collected for the detector facing the western wall in the basement of the hall, is shown in Fig. 3. The dominant contributor to the gamma spectrum is the $1461\text{\,}\mathrm{keV}$ line from the decay of 40K Chen2017 . To compare with the measurement, a very simple germanium detector was implemented in the Geant4 simulation. For reasons of computational expense, only gammas from a small (0.6\text{\,}\mathrm{m}$\times$0.6\text{\,}\mathrm{m}) region of the concrete wall were simulated, and detector efficiencies were not accounted for in the simulation. A scaling factor was therefore applied to match the measured rate. Nonetheless, the simulated spectrum (see Fig. 3) qualitatively replicates the important features of the measured gamma spectrum. Overall, the simulation is able to adequately approximate the spectrum of environmental gammas in the spectrometer hall.

4 Background measurements

During the summer of 2015, KATRIN proceeded with a measurement phase with only the MS and the FPD system, lasting several months. One of the primary goals of this commissioning campaign was to study background events originating from the MS.

As described in the following subsections, the electron rate was measured by the FPD for several configurations that modified the flux of environmental gammas at the MS. The voltage settings used for these measurements are given in Table 3. The electron rates were measured for the energy region of interest (ROI), which was defined as the range between $3\text{\,}\mathrm{keV}$ below and $2\text{\,}\mathrm{keV}$ above the expected electron energy at the detector. An asymmetric energy window is applied since the shape of the electron peak has a low-energy tail. For an electron with initial energy $E_{\textrm{0}}\approx$0\text{,}\mathrm{eV}$$ emitted at the MS surface, the expected energy is set by the change in the electric potential between the detector wafer and the MS inner electrode:

[TABLE]

From the values given in Table 3, the expected electron energy is $\sim 28.7\text{\,}\mathrm{keV}$. Due to a broken preamplifier module in the detector readout, six FPD pixels (each located in a different pixel ring) were not functional. The rate for each detector pixel ring with a missing pixel is linearly scaled by a corrective factor of $\frac{12}{11}$.

Several magnetic field configurations were utilized when studying the gamma-induced background. A symmetric magnetic field, similar to the planned field for standard KATRIN operation, is shown in the upper panel of Fig. 4. In this setting, the vast majority of electrons emitted from the MS surface should not reach the detector due to the inherent magnetic shielding of the flux tube PhDWandkowsky2013 . With this symmetric field setting, a total background rate of $0.561\pm 0.005$ counts per second (cps) was measured.

In order to directly study electrons emitted from the surface, an asymmetric field configuration can be used by adjusting the currents applied to the air-coils surrounding the MS; an example field is shown in the lower panel of Fig. 4. Because the magnetic field lines directly connect the vessel walls to the FPD, this setting allows a significantly larger fraction of surface electrons to be detected. To ensure the measured electrons are emitted from a well-defined surface region of the MS, certain detector pixels are excluded from the rate analysis. For the Asym. M configuration, the inner 16 detector pixels are excluded, while the inner 28 (4) pixels are excluded for the Asym. U (Asym. D) configuration.

4.1 Enhancement of gamma flux

To increase the gamma-induced background, a 60Co source with a total activity of $53.3\pm 2.7\text{\,}\mathrm{MBq}$ was positioned in the vicinity of the MS (Fig. 5). The source was originally used for geological surveys along underground piping and is therefore equipped with a shielded scintillator detector. It is housed in a lead-shielded transport container that allows convenient handling and transport when not in active use. Measurements were completed with the source partially outside its lead shielding (“open” configuration) and completely inside its lead shielding (“closed” configuration). The decay of 60Co primarily results in the cascade emission of two gammas at $1173\text{\,}\mathrm{keV}$ and $1332\text{\,}\mathrm{keV}$ Browne2013 .

Measurements of the electron rate with the FPD system were performed with the source located at different locations near the MS. However, extended measurements with the “closed” configuration were only performed at one position: the 60Co source located under the west side of the MS, approximately equidistant from the ends of the vessel (see Fig. 2). Therefore, only results from this position are presented here.

The FPD rate was measured for several magnetic-field and electrostatic shielding configurations; the rates can be found in Table 4. The effect of the 60Co source can clearly be seen in the bottom-right portion of the detector wafer in the left and center panels of Fig. 6 (asymmetric magnetic field), but is absent in the right panel (symmetric magnetic field).

The measurements with the 60Co source demonstrate the effectiveness of the shielding inside the MS against gamma-induced electron backgrounds. The background rate due to the 60Co source under an asymmetric magnetic configuration dropped from $231.7\pm 0.9\text{\,}\mathrm{cps}$ to $26.8\pm 0.3\text{\,}\mathrm{cps}$ with the addition of electrostatic shielding (changing $\Updelta U_{\textrm{IE}}$ from $0\text{\,}\mathrm{V}$ to $-100\text{\,}\mathrm{V}$). A further reduction in the rate by at least three orders of magnitude occurred when switching to the symmetric magnetic configuration ($0.005\pm 0.006\text{\,}\mathrm{cps}$, corresponding to less than $1\text{\,}\%$ of the total background rate). Though a large number of 60Co-induced secondary electrons were emitted from the MS surface, no significant rate effect was observed with the nominal magnetic field setting.

A simulation of the gamma flux from the 60Co source was performed using the geometry described in Section 3.2. Table 5 shows the simulated fluxes for gammas traversing the inner surface of the MS vessel. The presence of the 60Co source increases the gamma flux through the entire MS surface by about a factor of 8. The change in flux due to the open 60Co source is plotted as a function of axial position along the MS in Fig. 7. The distributions for the measured electron rate (asymmetric field setting) and the simulated gamma flux exhibit similar shapes and peak at the same axial position.

4.2 Suppression of gamma flux

In an effort to reduce the gamma flux originating from the bottom floor of the spectrometer hall, water shielding was temporarily added below the MS. The basin beneath the MS ($24.1\text{\,}\mathrm{m}$ long by $5.6\text{\,}\mathrm{m}$ wide) was filled with water to a depth of $20\text{\,}\mathrm{c}\mathrm{m}$. Additionally, a total of four flexible water tanks (each approximately $6.5\text{\,}\mathrm{m}$ long by $3.2\text{\,}\mathrm{m}$ wide by $0.6\text{\,}\mathrm{m}$ high) were installed next to the basin to increase the shielded area. The water tanks can be seen in Fig. 8.

The background rates for measurements with and without water shielding are shown in Table 6. Two asymmetric field settings, with field lines intersecting different regions of the MS surface, were implemented; similar reductions in the electron rate due to the water shielding were found for both ($\sim 0.4\text{\,}\mathrm{\char 37\relax}$). For the symmetric magnetic field setting, the shielding had no significant effect on the electron rate.

To investigate the effect of the shielding on the gamma flux, water was added to the simulation geometry using the dimensions cited above (see Fig. 2). The simulated fluxes are shown in Table 5. The addition of water shielding reduced the gamma flux inside the MS by about $7\text{\,}\mathrm{\char 37\relax}$. The change in flux due to the water shielding is plotted as a function of axial position along the MS in Fig. 9. The simulated gamma flux is mostly flat across the measured range and roughly matches the measured electron rate distribution (which is statistics-limited).

5 Gamma-induced background contribution

By combining the asymmetric magnetic field measurements and simulation results, it is possible to compute two quantities of interest: the secondary electron yield for gamma radiation traversing the inner surface of the MS, and the fraction of secondary electrons caused by environmental gamma radiation. The symmetric field measurements provide a way to determine the gamma-induced background contribution under standard operating conditions in future $m_{\overline{\nu}_{e}}$ measurements.

The relevant values used to calculate these quantities are listed in Table 7. Three different rates are to be distinguished: the measured FPD electron rate ($R$), the calculated gamma-induced MS electron emission rate ($S$), and the calculated gamma flux through the MS ($\Phi$). The values of $R$, $S$, and $\Phi$ can be defined as follows:

[TABLE]

where “env” indicates the contribution from environmental gamma emitters (e.g. concrete) and “other” indicates the contribution from other backgrounds (e.g. cosmic-ray muons). Eq. (missing) 2 requires that the 60Co source is closed and no water shielding is present.

The change in rate due to the 60Co source is defined to be

[TABLE]

while for water shielding the change in rate is

[TABLE]

Identical formulations hold for $\Delta S$ and $\Delta\Phi$.

5.1 Secondary electron yield

The yield $Y$ is the gamma-induced electron rate divided by the gamma flux through the same surface. This can be computed from the effect of the 60Co source or water shielding according to the following equation:

[TABLE]

$\Delta S$ can be computed from $\Delta R$ after including electron transport and detection efficiencies:

[TABLE]

where $\epsilon=$0.950\pm 0.028$$ is the FPD detection efficiency Amsbaugh2015 (ignoring the effect of backscattering from the detector surface PhDRenschler2011 ) and $P_{\text{arrival}}$ is the average arrival probability for electrons.

Because of the magnetic mirror effect, electrons emitted from the MS surface have a small probability to reach the FPD, which depends on their initial energy and emission angle relative to the magnetic field direction. $P_{\text{arrival}}$ was calculated using Kassiopeia Furse2017 , the particle-tracking simulation software developed by the KATRIN collaboration. For each magnetic field configuration, $1.6\text{\times}{10}^{5}$ electrons were started on the MS surface with emission angles sampled from a cosine angular distribution Henke1977 ; Seiler1983 ; Furman2002 . The electron energy spectrum was assumed to have the following form Chung1974 ; Seiler1983 ; Joy2004 :

[TABLE]

where $E$ is the electron energy and $W$ = $3.5\text{\,}\mathrm{eV}$ PhDBehrens2016 is the work function of the MS surface. The energy spectrum for true-secondary electrons is treated independently of the primary particle energy; the shape of the spectrum is known not to vary significantly with respect to the incident photon energy Henke1977 ; Seiler1983 . The validity of these assumptions on the electrons’ initial properties was shown in the context of the background analysis of cosmic-ray muons Altenmueller2019 .

Using the 60Co measurement, one finds $Y=$3\text{\times}{10}^{-4}\text{,}\mathrm{e}^{-}\mathrm{/}\upgamma$$. However, the yields derived from the water-shielding measurements give consistent values which are a factor of 2.6 larger than the 60Co result (see Table 7). One can attempt to explain the discrepancy between the two measurements by recognizing that although the shape of the secondary electron spectrum is independent of the gamma energy, the scaling factor for the electron spectrum is energy dependent Grudskii1984 . Because the gammas emitted from the 60Co source and the gammas blocked by the water shielding have ostensibly different spectral shapes, an accompanying disparity in electron yields would be expected. However, an analysis of the simulation results does not substantiate this expectation; inside the spectrometer, the gammas from the 60Co source and those blocked by the water shielding have very similar average energies ($676\text{\,}\mathrm{k}\mathrm{e}\mathrm{V}$ and $674\text{\,}\mathrm{k}\mathrm{e}\mathrm{V}$, respectively).

The difference in measured yields is likely the result of an incorrect value for $\Delta\Phi$ obtained from simulations; the tension can be alleviated by decreasing the simulated effect of the 60Co source or by increasing the simulated effect of the water shielding. The latter possibility would entail that the gamma emission from the concrete floor is underestimated in the simulations.

5.2 Fraction of secondary electrons induced by gammas

The fraction of secondary electrons emitted from the MS surface which are caused by environmental gamma radiation can be computed in the following manner:

[TABLE]

Assuming that $R_{\text{env}}$ is proportional to the flux of gammas in the MS, the following relation applies:

[TABLE]

The gamma-induced fraction can thus be obtained by combining Eqns. 8 and 9:

[TABLE]

Table 7 shows the values of $f_{\text{env}}$ calculated from the 60Co measurements under the two electrostatic shielding conditions, as well as from the gamma suppression measurements with water shielding. The results indicate that less than $6\text{\,}\mathrm{\char 37\relax}$ of secondary electrons emitted from the MS surface are induced by environmental gammas. However, the values from the 60Co and water shielding measurements differ by a factor of 2.5. The scale of this discrepancy is equivalent to the difference in the electron yields between the two types of measurements, as discussed in the previous section.

5.3 Gamma-induced background rate under standard conditions

Similar to the asymmetric field measurements, it is possible to use Eqn. 9 to determine the effects of environmental gamma radiation under symmetric field conditions. Applying the measured and simulated rates listed in Table 8, one finds that $R_{\text{env}}=$0.7\pm 0.9\text{,}\mathrm{mcps} (millicount per second), which is consistent with zero. Assuming the rate is Gaussian, one can follow the unified approach [Feldman1997 ](#bib.bib40) and set an upper limit on the gamma-induced background rate, obtaining $R_{\text{env}}\leq$2.2\text{\,}\mathrm{mcps} ($90\text{\,}\mathrm{\char 37\relax}\text{\,}\mathrm{C.L.}$).

However, one must account for the discrepancy in the results between the 60Co source and water shielding measurements, as mentioned in the previous sections. A conservative approach is to allow for the possibility that the simulation overestimates the flux of gammas through the MS from the 60Co source by a factor of $2.6$. In this case, one finds

[TABLE]

Given a nominal rate of $561\text{\,}\mathrm{mcps}$, this result indicates that less than $\sim 1\text{\,}\mathrm{\char 37\relax}$ of the MS background rate can be attributed to environmental gamma radiation.

A similar procedure was followed with respect to the water shielding data, giving a limit of $R_{\text{env}}\leq$110\text{,}\mathrm{mcps}$~{}($90\text{,}\mathrm{\char 37\relax}\text{,}\mathrm{C.L.}$)$111 Here, a corrective factor of $(2.6)^{-1}$ (which accounts for the possibility that the simulation underestimates the effect of the water shielding) was not applied, in order to obtain a conservative result., which is a significantly weaker limit than that obtained from the 60Co measurements. The weaker limit obtained with the water shielding configuration is a result of the small simulated value of $\Delta\Phi$ and the large measured uncertainty on $\Delta R$. Therefore, the measurements with the 60Co source are more sensitive to the effect of environmental gammas than the water shielding measurements.

6 Discussion and conclusions

Low-energy background electrons produced inside the MS are indistinguishable from signal $\upbeta$-particles. Thus, it is necessary to understand and limit the various sources of electrons in the MS. Measurements using the asymmetric magnetic field, combined with simulation, indicate that less than $6\text{\,}\mathrm{\char 37\relax}$ of secondary electrons emitted from the MS surface are induced by environmental gammas.

The measurements with the 60Co source show that the electrostatic and magnetic shielding are highly effective in mitigating the effect from secondary electrons. Changes in the flux of environmental gammas have little if any effect on the MS background rate under standard operating conditions. When combined with simulation, the results indicate that less than $5.6\text{\,}\mathrm{mcps}$ ($90\text{\,}\mathrm{\char 37\relax}\text{\,}\mathrm{C.L.}$) of the MS background rate is gamma-induced, corresponding to less than $\sim 1\text{\,}\mathrm{\char 37\relax}$ of the total rate.

The remaining (and much more dominant) background rate seems to be caused by the ionization of Rydberg atoms PhDTrost2018 ; Fraenkle2017 . These atoms are highly excited neutral atoms which can penetrate the electromagnetic shielding implemented in KATRIN. Furthermore, thermal radiation at room temperature is energetic enough to ionize these atoms and create background electrons. The production of these Rydberg atoms is thought to primarily originate from the decay of 210Pb within the MS walls PhDHarms2015 . Plans to mitigate this background are currently being studied.

Acknowledgements.

We acknowledge the support of Helmholtz Association (HGF), Ministry for Education and Research BMBF (5A17PDA, 05A17PM3, 05A17PX3, 05A17VK2, and 05A17WO3), Helmholtz Alliance for Astroparticle Physics (HAP), and Helmholtz Young Investigator Group (VH-NG-1055) in Germany; Ministry of Education, Youth and Sport (CANAM-LM2011019), cooperation with the JINR Dubna (3+3 grants) 2017–2019 in the Czech Republic; and the Department of Energy through grants DE-FG02-97ER41020, DE-FG02-94ER40818, DE-SC0004036, DE-FG02-97ER41033, DE-FG02-97ER41041, DE-AC02-05CH11231, DE-SC0011091, and DE-SC0019304 in the United States.