Neutrino PhysicsView this Special Issue
Review Article | Open Access
G. Drexlin, V. Hannen, S. Mertens, C. Weinheimer, "Current Direct Neutrino Mass Experiments", Advances in High Energy Physics, vol. 2013, Article ID 293986, 39 pages, 2013. https://doi.org/10.1155/2013/293986
Current Direct Neutrino Mass Experiments
In this contribution, we review the status and perspectives of direct neutrino mass experiments, which investigate the kinematics of -decays of specific isotopes (3H, 187Re, 163Ho) to derive model-independent information on the averaged electron (anti)neutrino mass. After discussing the kinematics of -decay and the determination of the neutrino mass, we give a brief overview of past neutrino mass measurements (SN1987a-ToF studies, Mainz and Troitsk experiments for 3H, cryobolometers for 187Re). We then describe the Karlsruhe Tritium Neutrino (KATRIN) experiment currently under construction at Karlsruhe Institute of Technology, which will use the MAC-E-Filter principle to push the sensitivity down to a value of 200 meV (90% C.L.). To do so, many technological challenges have to be solved related to source intensity and stability, as well as precision energy analysis and low background rate close to the kinematic endpoint of tritium -decay at 18.6 keV. We then review new approaches such as the MARE, ECHO, and Project8 experiments, which offer the promise to perform an independent measurement of the neutrino mass in the sub-eV region. Altogether, the novel methods developed in direct neutrino mass experiments will provide vital information on the absolute mass scale of neutrinos.
The various experiments with atmospheric, solar, accelerator, and reactor neutrinos [1–5] provide compelling evidence that neutrino flavor states are nontrivial superpositions of neutrino mass eigenstates and that neutrinos oscillate from one flavor state into another during flight. By these neutrino oscillation experiments, we can determine the neutrino mixing angles and the differences between the squares of neutrino masses. In the case of the so-called solar or reactor mass splitting , we not only know the modulus of this difference but also its sign. Clearly these findings prove that neutrinos have nonzero masses, but neutrino oscillation experiments being a kind of interference experiment cannot determine absolute masses. We may parameterize our missing knowledge by a free parameter , the mass of the smallest neutrino mass eigenstate (see Figure 1).
We should note that throughout this paper, we will not distinguish between the mass of a neutrino and the mass of an antineutrino, which should be the same if the CPT theorem holds. Therefore, we will use the term neutrino when we speak of neutrinos and of antineutrinos. But we will explain for each measurement whether the result is obtained for neutrinos or antineutrinos.
The absolute value of the neutrino masses is very important for astrophysics and cosmology because of the role of neutrinos in structure formation due to the huge abundance of relic neutrinos left over in the universe from the big bang (336/cm3) . In addition, the key role of neutrino masses in understanding, which of the possible extensions or new theories beyond the Standard Model of particle physics is the right one [7, 8], makes the quest for the absolute value of the neutrino mass among of the most urgent questions of nuclear and particle physics.
Three different approaches can lead to the absolute neutrino mass scale as follows.
(i) Cosmology. Today’s visible structure of the universe has been formed out of fluctuations of the very early universe. Due to the large abundance of relic neutrinos and their low masses they acted as hot dark matter: neutrinos have smeared out fluctuations at small scales. How small or large these scales are is described by the free streaming length of the neutrinos which depends on their mass. By determining the early fluctuations imprinted on the cosmic microwave background with the WMAP satellite  and mapping out today’s structure of the universe by large galaxy surveys like SDSS  conclusions on the sum of the neutrino masses can be drawn. Up to now, only upper limits on the sum of the neutrino masses have been obtained around eV , which are to some extent model and analysis dependent .
(ii) Neutrinoless Double -Decay (0νββ). A neutrinoless double -decay (two -decays in the same nucleus at the same time with emission of two -electrons (positrons) while the antineutrino (neutrino) emitted at one vertex is absorbed at the other vertex as a neutrino (antineutrino)) is forbidden in the Standard Model of particle physics. It could exist, if the neutrino is its own antiparticle (“Majorana-neutrino” in contrast to “Dirac-neutrino”) . Furthermore, a finite neutrino mass is the most natural explanation to produce in the chirality-selective interaction a neutrino with a small component of opposite handedness on which this neutrino exchange subsists. Then the decay rate will scale with the absolute square of the so-called effective Majorana neutrino mass, which takes into account the neutrino mixing matrix Here represents the sum of the neutrino masses contribution coherently to the -decay. Hence this coherent sum carries their relative phases (the usual CP-violating phase of an unitary mixing matrix plus two so-called Majorana-phases). A significant additional uncertainty which enters the relation of and the decay rate is the nuclear matrix element of the neutrinoless double -decay . There is one claim for evidence at eV by part of the Heidelberg-Moscow collaboration , which is being challenged by limits from different experiments in the same range, for example, very recently by the EXO-200 experiment .
(iii) Direct Neutrino Mass Determination. The direct neutrino mass determination is based purely on kinematics without further assumptions. Essentially, the neutrino mass is determined by using the relativistic energy-momentum relationship . Therefore it is sensitive to the neutrino mass squared . In principle there are two methods: time-of-flight measurements and precision investigations of weak decays. The former requires very long baselines and therefore very strong sources, which only cataclysmic astrophysical events like a core-collapse supernova could provide. The supernova explosion SN1987a in the Large Magellanic Cloud gave limits of 5.7 eV (95% C.L.)  or of 5.8 eV (95% C.L.)  on the neutrino mass depending somewhat on the underlying supernova model. Unfortunately nearby supernova explosions are too rare and seem to be not well enough understood to allow to compete with the laboratory direct neutrino mass experiments.
Therefore, aiming for this sensitivity, the investigation of the kinematics of weak decays and more explicitly the investigation of the endpoint region of a -decay spectrum (or an electron capture) is still the most sensitive model-independent and direct method to determine the neutrino mass. Here the neutrino is not observed but the charged decay products are precisely measured. Using energy and momentum conservation, the neutrino mass can be obtained. In the case of the investigation of a -spectrum usually the “average electron neutrino mass” is determined (see (20) in the next subsection) In contrast to in neutrinoless double -decay (see (1)), this sum averages over all neutrino mass states contributing to the electron neutrino and no phases of the neutrino mixing matrix enter. The decay into the different neutrino mass eigenstates add incoherently, which we will discuss in more detail for the neutrino mixture to sterile neutrinos in Section 2.1.
Figure 2 demonstrates that the different methods are complementary to each other and compares them.
This paper is structured as follows. In Section 2, the neutrino mass determination from the kinematics of -decay is described. Section 3 presents the analysis of the spectrum of neutrinos from supernova SN1987a and the recent -decay experiments in search for the neutrino mass scale. In Section 4, an overview of the present KATRIN experiment is given. New approaches to directly determine the neutrino mass are presented in Section 5. This paper ends with a conclusion in Section 6. For a more detailed and complete overview on this subject, we would like to refer to the reviews [18–23].
2. -Decay and -Mass
According to Fermi’s Golden Rule, the decay rate for a -decay is given by the square of the transition matrix element summed and integrated over all possible discrete and continuous final states (we use the convention for simplicity). Let us first calculate the density of the final states. The number of different final states of outgoing particles inside a normalization volume into the solid angle with momenta between and , or, respectively, with energies in the corresponding interval around the total energy , is This gives a state density per energy interval and solid angle of
Since the mass of the nucleus is much larger than the energies of the two emitted leptons, we can use for the next steps the following simplification: the nucleus takes nearly no energy but balances all momenta (we will consider the recoil energy of the nucleus later.). Therefore we need to count the state density of the electron and the neutrino only
The transition matrix element can be divided into a leptonic part, , and a nuclear one, . Usually the coupling is written separately and expressed in terms of Fermi’s coupling constant and the Cabibbo angle
2.1. Allowed and Superallowed Transitions
We first discuss the case of allowed or superallowed decays like that of tritium. Here, none of the leptons has to carry away angular momentum. Hence, the leptonic part essentially results in the probability of the two leptons to be found at the nucleus, which is for the neutrino and for the electron, yielding The Fermi function takes into account the final electromagnetic interaction of the emitted -electron with the daughter nucleus of nuclear charge . The Fermi function is approximately given by  with the Sommerfeld parameter .
For an allowed or superallowed transition the nuclear matrix element is independent of the kinetic energy of the electron. The coupling of the lepton spins to the nuclear spin is usually contracted into the nuclear matrix element. This nuclear matrix element of an allowed or superallowed transition can be divided into a vector current or Fermi part () and into an axial current or Gamov-Teller part ( but no ). In the former case, the spins of electron and neutrino couple to , in the latter case to . What remains is an angular correlation of the two outgoing leptons. Since charge current weak interactions like -decay maximally violate parity they prefer—depending on velocity—negative helicities for particles and positive helicities for antiparticles. Thus the momenta or directions of the leptons are correlated with respect to their spins and therefore to each other. This results are an () angular correlation factor with the electron velocity and the neutrino velocity . The angular correlation coefficient amounts to for pure Fermi transitions and to for pure Gamov-Teller transitions within the Standard Model .
The phase space density (6) is distributed over a surface in the two-particle phase space which is defined by a -function conserving the decay energy. With this prescription, we can integrate (3) over the continuum states and get the partial decay rate into a single channel; for instance, the ground state of the daughter system with probability
In direct neutrino mass measurements usually the formulas are given in terms of the kinetic energy of the electron The maximal kinetic energy of the electron for the case of zero neutrino mass zero is called endpoint energy which is defined by a vanishing neutrino energy
A correct integration over the unobserved neutrino variables in (11) has to respect the () angular correlation factor (10), which also has to be considered in calculating the exact recoil energy of the nucleus . If we consider that the -electrons of interest have a certain minimal kinetic energy then we can calculate the range of recoil energies of the daughter nucleus of mass : the recoil energy is bound upwards by the case, in which the outgoing electron takes the maximum kinetic energy and downwards by the case, in which the electron of kinetic energy is emitted opposite to the direction of the neutrino, which has in this case a momentum (neglecting for a moment the nonzero value of the neutrino mass) Due to the largeness of even for sizeable electron energy ranges below the endpoint , according to (14) the recoil energy does not change much (numbers are given for the example of tritium in the next section). Therefore, for the region of interest below the endpoint , we can apply a constant recoil energy correction . and (13) becomes (We should note here that we cannot use (15) to derive the endpoint energy from a measured nuclear -value with the precision required for the direct measurement of the neutrino mass due to the uncertainty of , which is eV at best. Therefore has to be fitted from the -spectrum together with the neutrino mass squared. For the case of tritium there is currently a large experimental effort to improve significantly the precision on the -value of tritium -decay by ultrahigh precision ion cyclotron resonance mass spectroscopy in a multi-Penning trap setup measuring the 3He-3H mass difference  with the final goal to use the measured -value in the neutrino mass fit.)
Further integration over the angles yields through (10) an averaged nuclear matrix element, as mentioned above. Besides integrating over the ()-continuum, we have to sum over all other final states. For a -decaying atom or molecule it is a double sum: one summation runs over all neutrino mass eigenstates with probabilities which are kinematically accessible (). The second summation has to be done over all of the electronic final states of the daughter system with probabilities and excitation energies . These comprise excitations of the electron shell but also—in the case of -decaying molecules—rotational and vibrational excitations. These excitations are caused by the sudden change of the nuclear charge from to which requests a rearrangement of the electronic orbitals of the daughter atom or molecule and the interatomic distances in case of a molecule. They give rise to shifted endpoint energies. Introducing the definition the total neutrino energy now amounts for this excitation to . The -electrons are leaving the nucleus on a time scale much shorter than the typical Bohr velocities of the shell electrons of the mother isotope. Therefore, the excitation probabilities of electronic states—and of vibrational-rotational excitations in the case of molecules—can be calculated in the so-called sudden approximation from the overlap of the primary electron wave function with the wave functions of the daughter ion
Rather than in the total decay rate, we are interested in its energy spectrum , which we can read directly from (11) without performing the second integration over the -energy. Using and summing up over the final states it reads
The -function confines the spectral components to the physical sector . This causes a technical difficulty in fitting mass values smaller than the sensitivity limit of the data, as statistical fluctuations of the measured spectrum might occur which can no longer be fitted within the allowed physical parameter space. Therefore, one has to define a reasonable mathematical continuation of the spectrum into the region which leads to -parabolas around (see, e.g., ).
Assuming unitarity of the kinematic accessible neutrino mass states (), we can expand the second line of (18) for This average over the squared masses of the neutrino mass states in (20) defines what we called the electron neutrino mass in (2). This simplification always applies, if we cannot resolve the different neutrino mass states.
Figure 3 shows the -spectrum at the endpoint according to (18). The influence of the neutrino mass on the -spectrum shows only at the upper end below , where the neutrino is not fully relativistic and can exhibit its massive character. The relative influence decreases in proportion to (see Figure 3), which leads far below the endpoint—according to (20)—to a small constant offset proportional to .
Concerning the various neutrino mass states, we can also assume that there is one heavy neutrino mass state (this heavy state might comprise more than one heavy state, which are experimentally not distinguishable) and one light one (again this could be the sum of more than one light neutrino). Such a situation could arise, if 3 light active neutrinos and one heavy sterile neutrino are mixed. With and , we can rewrite the last line of (18) for this case into
Figure 3 defines the requirements of a direct neutrino mass experiment which investigates a -spectrum: the task is to resolve the tiny change of the spectral shape due to the neutrino mass in the region just below the endpoint , where the count rate is going to vanish. Therefore a high sensitivity experiment requires high energy resolution, large -decay source strength and acceptance, and low background rate.
Now we should firstly discuss, what is the best -emitter for such a task. Figure 4 shows the total count rate of a superallowed -emitter as function of the endpoint energy. Of course, the total count rate rises strongly with , while the relative fraction in the last 10 eV below decreases. Interestingly, the total count rate in the last 10 eV below , which we can take as a measure of our energy region of interest for determining the neutrino mass, is increases with regard to . This increase is caused by the larger phase space for th -electron. From Figure 4 one might argue that the endpoint energy does not play a significant role in selecting the right -isotope, but we have to consider the fact that we need a certain energy resolution to determine the neutrino mass. Experimentally it makes a huge difference, whether we have to achieve a certain at a low energy or at a higher one. Secondly, the -electrons of no interest with regard to the neutrino mass could cause experimental problems (e.g., as background or pileup) and again this argument favors a low .
2.2. Tritium -Decay
The heaviest of the hydrogen isotopes tritium undergoes -decay with a half-life of 12.3 y. Tritium and Helium-3 are mirror nuclei of the same isospin doublet; therefore, the decay is superallowed. Thus the nuclear matrix element for tritium is close to that of the -decay of the free neutron and amounts to 
With an endpoint energy of 18.6 keV, it has one of the lowest endpoints of all emitters together with a reasonable long half-life. Its superallowed shape of the -spectrum and its simple electronic structure allow the tritium -spectrum to be measured with small systematic uncertainties.
The recoil correction for tritium is not an issue. Up to now all tritium -decay experiments used molecular tritium, which give a maximal recoil energy to the daughter molecular ion of eV. Even for the most sensitive tritium -decay experiment, the upcoming KATRIN experiment (see Section 4.1), the maximum variation of over the energy interval of investigation (the last 30 eV below the endpoint) amounts to meV only. It was checked  that this variation can be neglected and the recoil energy can be replaced by a constant value of eV, yielding a fixed endpoint according to (15).
Furthermore, one may apply radiative corrections to the spectrum [28, 29]. However, they are quite small and would influence the result on even for the KATRIN experiment only by a few percent of its present systematic uncertainty. One may also raise the point of whether possible contributions from right-handed currents might lead to measurable spectral anomalies [30, 31]. It has been checked that the present limits on the corresponding right-handed boson mass  rule out a sizeable contribution within present experimental uncertainties. Even the forthcoming KATRIN experiment will hardly be sensitive to this problem [32, 33].
Concerning the calculation of the electronic final states according to (17), we have to consider molecular tritium since all tritium -decay experiments so far have been using molecular tritium sources, containing the molecule . The wave functions of the tritium molecule are much more complicated, since in addition to two identical electrons they comprise also the description of rotational and vibrational states, which may be excited during the -decay as well. Figure 5 shows a recent numerical calculation of the final states of the molecule. The transition to the electronic ground state of the daughter ion as well as the transition to higher excited electronic states are not sharp in energy, but broadened due to rotational-vibrational excitations. More recent calculations agree with these results [35, 36].
The first group of excited electronic states starts at around eV . Therefore excited states play almost no role for the energy interval considered for the upcoming KATRIN experiment: only the decay to the ground state of the daughter molecule, which is populated with about 57% probability, has to be taken into account. Due to the nuclear recoil, however, a large number of rotational-vibrational states with a mean excitation energy of 1.7 eV and a standard deviation of 0.4 eV are populated. These values hold for a pure T2 source without contamination by other hydrogen isotopes. But a contamination of the T2 molecules by DT or HT molecules does not matter in first order: the shift of the mean rotational-vibrational excitation of HT with respect to T2 is compensated by a corresponding change of the nuclear recoil energy of HT with respect to the 1.5 times heavier T2 molecule  (see Figure 5(b)).
Summarizing the properties of tritium for direct neutrino mass measurements: it is the standard isotope for this kind of study due to its low endpoint of 18.6 keV, its rather short half-life of 12.3 y, its superallowed shape of the -spectrum, and its simple electronic structure. Tritium -decay experiments using a tritium source and a separated electron spectrometer have been performed in search for the neutrino mass for more than 60 years. But when the electron spectrometer is identical to the -source, the situation is different and a -isotope with an even lower endpoint energy is preferred, even if it does not have an allowed decay.
2.3. Forbidden Transitions Like 187Re
The isotope 187Re exhibits the lowest endpoint energy with keV of all known -emitters decaying to the ground-state of the daughter nucleus (there is even a decay of the into an excited nuclear -state of the daughter nucleus with a much lower -value of eV [37, 38]. But there are two reasons, why such a decay cannot be used for a direct neutrino mass measurement: this partial decay has an ultra-long half-life of y  and the signature of the neutrino mass is hidden in the -spectrum of the decay into the ground state of ). The ground state of the mother isotope 187Re has spin and parity . The -decay goes to the ground state of the daughter 187Os with spin and parity . Therefore, the decay is a first unique forbidden transition. The lepton pair, electron and antineutrino, has to carry away two units of angular momentum and has to change parity. The two leptons will couple to spin in this case, one unit of orbital momentum has to be carried away by either the electron or the antineutrino. Therefore the half life of y is huge and about as long as the age of the universe. The advantage of the 7 times lower endpoint energy = 2.47 keV of 187Re with respect to tritium does not compensate the fact that one needs a large number of 187Re atoms to obtain enough count rate near the endpoint to measure the neutrino mass. Therefore, a classical source ≠ spectrometer arrangement like for the tritium experiments is not feasible for 187Re, because the -electrons will undergo too many inelastic scattering processes within the 187Re source. Secondly, the isotope 187Re has a complicated electron shell and the electronic final states might not be calculable precisely enough.
Therefore, the -decay of 187Re can only be exploited if the -spectrometer measures the entire released energy, except that of the neutrino but including the energy loss by inelastic scattering processes and electronic excitations. This situation can be realized by using a cryogenic bolometer as the -spectrometer, which at the same time contains the -emitter 187Re  (see Section 3.4).
We discuss now the consequences for the -spectrum. Since either the electron or the antineutrino has to be emitted with orbital angular momentum we cannot, expand the plain wave of this lepton to zeroth order anymore as we did in (8), but we have to go to first order In contrast to an allowed decay, the matrix element will become dependent on energy. According to (25), an additional factor proportional to or to will occur, depending on whether the electron or the antineutrino carries away the unit of orbital angular momentum. In the decay rate the square of these momenta will appear. For both cases, a Fermi function or needs to be considered which describes the Coulomb interaction of the outgoing electron in a or state with the remaining osmium ion  where is the nuclear radius. The first term of (26) proportional to is by 4 orders of magnitude larger than the second term proportional to . The nuclear matrix element is more complex than that of an allowed -decay.
In contrast to the case of tritium, we do not account in (26) for excited electronic final states, since we assume that all losses by electromagnetic excitations will be added to the energy of the -electron as well as to the recoil energy in the source = detector arrangement by the signal integration of the cryobolometer. Thus the -spectrum looks simpler (we will discuss later that some excited states may live longer than the signal integration time of the cryobolometer and that the corresponding excitation energy may not be measured, which would cause systematic uncertainties). But this is only true up to first order. In second order, the electronic final states with excitation energy and probability have to be considered since the modification of the phase space of the outgoing electron and the squared matrix element have to be taken into account.
We will discuss the influence of electronic final states for the case of cryobolometers in some detail: when an electronic final state takes the excitation energy , in the calculation for an allowed decay (see (18)), we had just shifted the effective endpoint energy by this amount () and multiplied this fraction of the -spectrum with its probability . Thus the whole -spectrum including the phase space of the outgoing leptons was calculated correctly up to possible electron energies . When we measure with a cryobolometer the sum energy in the case of an electronic excitation , the true kinetic energy of the electron amounts only to . The residual energy release detected in the cryobolometer stems from the deexcitation of the electronic excitation. Therefore the -decay probability, or the corresponding phase space factor and the squared matrix element have to be calculated for the true electron kinetic energy and the reduced endpoint energy as We can expand the relevant parameters for the phase space and squared matrix element calculation to first order assuming and Equation (26) then becomes For a typical electronic excitation eV and a typical kinetic energy of the electron keV the correction factor in (32) might have enough influence on the shape of the -spectrum that it needs to be considered in future high precision cryobolometer experiments.
Secondly for a 187Re atom within a crystalline environment the so-called beta-environmental fine structure (BEFS) has to be taken into account (see Section 3.5), which leads to a modulation of the -spectrum. Similarly to the extended X-ray absorption fine structure (EXAFS) this fine structure is caused by an interference of the outgoing wave of the -electron with waves scattered on the neighboring atoms.
3. Past Direct Neutrino Mass Experiments
3.1. Neutrino Mass Limit from Supernova SN1987a
On February 23, 1987, neutrinos from the supernova SN1987a in the Large Magellanic Cloud (LMC) have been observed by the three neutrino detectors Kamiokande II , IMB , and at Baksan  (see Figure 6). This core-collapse supernova emitted a total energy of about erg = J, 99% of this energy was released by neutrinos . These neutrinos traveled over a distance of about kpc (165,000 lyr). Although core collapse supernovae emit neutrinos and antineutrinos of all flavors at different phases of the collapse only electron antineutrinos from SN1987a have been detected via the famous inverse -decay reaction on the proton . Other neutrino reactions were suppressed by too high thresholds or too low cross-sections.
To calculate the time-of-flight, we first express the velocity of a relativistic particle by its mass and total energy as
The delay of the arrival of a supernova neutrino at earth with regard to a particle at speed of light can be expressed as function of the neutrino mass and its total energy as
Therefore we would expect the neutrino energy to follow hyperbolas for each neutrino mass state as function of the square root of the neutrino arrival time. Of course this only holds, if the neutrino emission is sharp in time with respect to the spread of arrival times on earth. The energy versus time spectra of the three experiments which detected neutrinos from supernova SN1987a do not exhibit a dependence following equation (34) (see Figure 6). Therefore only upper limits of 5.7 eV (95% C.L.)  or of 5.8 eV (95% C.L.)  on the neutrino mass can be deduced which depend somewhat on the underlying supernova model.
Nowadays more and larger neutrino detectors are online and capable to measure neutrinos from a galactic or nearby core-collapse supernova. They are interconnected by the SuperNova Early Warning System SNEWS . No core-collapse supernova from our galaxy or from a satellite galaxy like the LMC has been observed since then. The expected rate for galactic core-collapse supernovae is 2-3 per century. Even if a new supernova will yield many more neutrino events than supernova SN1987a it would be difficult getting a much more precise limit on the neutrino mass. The limiting factors are the time and energy spectra of the neutrino emission of a core-collapse supernova, which are not known well enough. One possibility to bypass this problem exists for core-collapse supernovae, if the core-collapse supernova forms a black-hole. Then the neutrino emission will stop abruptly  and this stamp on the latest neutrino start time could be used in the analysis. A possible limit of the present neutrino detectors would be in the eV range. Whether new detectors with larger masses and possibly lower energy thresholds would allow a higher sensitivity on the neutrino mass still needs to be studied.
3.2. Laboratory Direct Neutrino Mass Limits
The majority of the published direct laboratory results on originate from the investigation of -decays, which are sensitive to the average of the antineutrino mass states contributing to the electron antineutrino. The mass of the neutrino could in principle be accessed by the investigation of decays, but measuring electron capture decays is much more sensitive . By the investigation of the electron capture of 163Ho two groups obtained upper limits on the average mass of the electron neutrino of eV at 95% C.L.  and of eV at 68% C.L. . These experiments will be discussed in some detail in Section 5.3.
An exotic way of a direct neutrino mass measurement was using bound-state -decay, where the outgoing -electron is captured into a bound electronic state: totally ionized ions were circulating in a storage ring and undergoing bound-state -decay although neutral 163Dy atoms are stable. The measurement of this bound state -decay resulted in a limit on the neutrino mass of 410 eV at 68% C.L. .
Except 187Re -decay, on which the investigations have been started only within the last decade, all direct neutrino mass experiments using -decays were done with the isotope tritium. In the long history of these tritium -decay experiments, about a dozen results have been reported starting with the experiment of Curran in the late forties yielding keV .
In the beginning of the eighties a group from the Institute of Theoretical and Experimental Physics (ITEP) at Moscow [51, 52] claimed the discovery of a nonzero neutrino mass of around 30 eV. The ITEP group used as -source a thin film of tritiated valine combined with a new type of magnetic spectrometer, the Tretyakov spectrometer (a Tretyakov spectrometer uses a toroidal magnet field with a dependent strength like it is used for magnetic horns of neutrino beam facilities at accelerators to focus the secondary pions behind the proton target). The first tests of the ITEP claim came from the experiments at the University of Zürich  and the Los Alamos National Laboratory (LANL) . Both experiments applied similar Tretyakov-type spectrometers, but more advanced tritium sources with respect to the ITEP group. The Zürich experiment used a solid source of tritium implanted into carbon and later a self-assembling film of tritiated hydrocarbon chains. The LANL group was the first to develop a gaseous molecular tritium source avoiding solid state corrections. Both experiments disproved the ITEP result. The reason for the mass signal at ITEP was twofold: the energy loss correction was probably overestimated, and a 3He–T mass difference measurement  confirming the endpoint energy of the ITEP result, turned out only later to be significantly wrong [56, 57].
Also in the nineties tritium -decay experiments yielded controversially discussed results: Figure 7 shows the final results of the experiments at LANL and Zürich together with the results from other more recent measurements with magnetic spectrometers at University of Tokyo, Lawrence Livermore National Laboratory, and Beijing. The sensitivity on the neutrino mass had improved a lot but the values for the observable populated the nonphysical negative region. In 1991 and 1994 two new experiments started data taking at Mainz and at Troitsk, which used a new type of electrostatic spectrometer, so-called MAC-E-Filters, which were superior in energy resolution and luminosity with respect to the previous magnetic spectrometers. However, even their early data were confirming the large negative values of the LANL and Livermore experiments when being analyzed over the last 500 eV of the -spectrum below the endpoint . But the large negative values of disappeared when analyzing only small intervals below the endpoint . This effect, which could only be investigated by the high luminosity MAC-E-Filters, pointed towards an underestimated or missing energy loss process, seemingly to be present in all experiments. The only common feature of the various experiments seemed to be the calculations of the electronic excitation energies and excitation probabilities of the daughter ions. Different theory groups checked these calculations in detail. The expansion was calculated to one order further and new interesting insight into this problem was obtained, but no significant changes were found [34, 36].
Then the Mainz group found the origin of the missing energy loss process for its experiment. The Mainz experiment used as tritium source a film of molecular tritium quench-condensed onto aluminum or graphite substrates. Although the film was prepared as a homogenous thin film with flat surface, detailed studies showed  that the film was not stable: it underwent a temperature-activated roughening transition into an inhomogeneous film by forming microcrystals. Thus, unexpected large inelastic scattering probabilities were obtained, which were not taken into account in previous analyses. This extra energy losses were only significant when analyzing larger energy intervals below the endpoint.
The Troitsk experiment on the other hand used a windowless gaseous molecular tritium source, similar to the LANL apparatus. Here, the influence of large-angle scattering of electrons magnetically trapped in the tritium source was not considered in the first analysis. After correcting for this effect the negative values for disappeared.
The fact that more experimental results of the early nineties populate the region of negative values (see Figure 7) can be understood by the following consideration . For , (19) can be expanded into On the other hand the convolution of a -spectrum (18) at with a Gaussian of width leads to Therefore, in the presence of a missed experimental broadening with Gaussian width one expects a shift of the result on of which gives rise to a negative value of .
The significant improvement in the neutrino mass sensitivity by the Troitsk and the Mainz experiments are due to MAC-E-Filters (Magnetic Adiabatic Collimation with an Electrostatic Filter). This new type of spectrometer—based on early work by Kruit and Read —was developed for the application to the tritium -decay at Mainz and Troitsk independently [70, 71]. The MAC-E-Filter combines high luminosity at low background and a high energy resolution, which are essential features to measure the neutrino mass from the endpoint region of a -decay spectrum.
The main features of the MAC-E-Filter are illustrated in Figure 8: two superconducting solenoids are producing a magnetic guiding field. The -electrons, starting from the tritium source in the left solenoid into the forward hemisphere, are guided magnetically on a cyclotron motion along the magnetic field lines into the spectrometer yielding an accepted solid angle of nearly . On the way of an electron into the center of the spectrometer the magnetic field decreases smoothly by several orders of magnitude keeping the magnetic orbital moment of the electron invariant (equation given in nonrelativistic approximation) Therefore nearly all cyclotron energy is transformed into longitudinal motion (see Figure 8 bottom) giving rise to a broad beam of electrons flying almost parallel to the magnetic field lines. This is just the opposite of the so-called magnetic mirror or magnetic bottle effect.
This parallel beam of electrons is energetically analyzed by applying an electrostatic barrier created by a system of one or more cylindrical electrodes. The relative sharpness of this energy high-pass filter is only given by the ratio of the minimum magnetic field reached at the electrostatic barrier in the so-called analyzing plane and the maximum magnetic field between -electron source and spectrometer
It is beneficial to place the electron source in a magnetic field somewhat lower than the maximum magnetic field . Thus the magnetic mirror effect based on the adiabatic invariant (38) hinders electrons with large starting angles at the source and long paths inside the source to enter the MAC-E-Filter. Only electrons are able to pass the pinch field which exhibit starting angles at of In principle, the pinch magnet could also be installed between the MAC-E-Filter and the detector, which counts the electrons transmitted by the MAC-E-Filter, as long as the electron transport is always adiabatically. Such an arrangement has been realized at the KATRIN experiment.
The exact shape of the transmission function can be calculated analytically. For an isotropically emitting monoenergetic source of particles with kinetic energy and charge it reads as function of the retarding potential as
The -electrons are spiralling around the guiding magnetic field lines in zeroth approximation. Additionally, in non-homogeneous electrical and magnetic fields they feel a small drift , which reads in first order :
The clear advantages of the MAC-E-Filter of large angular acceptance and high energy resolution come together with the danger to store charged particles in Penning, magnetic mirror and combined traps . This problem and countermeasures will be discussed later at the example of the KATRIN experiment (see Sections 4.4 and 4.5).
A very interesting application is using the MAC-E-Filter in time-of-flight mode. This mode has the advantage to be nonintegrating but requests to measure to or restrict the start time of the electron under investigation. The analysis can be done on cutting on the time-of-flight  or fully making use of the individual measured time-of-flights .
The two tritium -decay experiments at Mainz and at Troitsk used similar MAC-E-Filters with an energy resolution of 4.8 eV (3.5 eV) at Mainz (Troitsk). The spectrometers differed slightly in size: the diameter and length of the Mainz (Troitsk) spectrometer are 1 m (1.5 m) and 4 m (7 m). The major differences between the two setups are the tritium sources: The Mainz Neutrino Mass Experiment used a novel condensed solid tritium source, whereas the experiment at Troitsk applied a windowless gaseous molecular tritium source similar to the ones of the experiments at Los Alamos and at Livermore before.
3.2.2. The Mainz Neutrino Mass Experiment
The tritium source was a thin film of molecular tritium, which was quench-condensed on a cold graphite substrate. By laser ellipsometry the film thickness was determined. Typically film thicknesses of 20 to 40 monolayers were applied before and of 150 monolayers after the upgrade of the experiment in 1995–1997. The retarding potential of the MAC-E-Filter was created by a complex system of cylindrical electrodes. The upgrade of the Mainz setup in 1995–1997 (Mainz phase II, see Figure 10) includes the installation of a new tilted pair of superconducting solenoids between the tritium source and the spectrometer and the use of a new cryostat to keep the temperature of the tritium film below 2 K. The first measure eliminated source-correlated background and allowed the source strength to be increased significantly. From this upgrade on, the Mainz experiment ran with a similar signal and similar background rate as the Troitsk experiment. The second measure avoided the roughening transition of the homogeneously condensed tritium films with time , which before had given rise to negative values of when the data analysis used large intervals of the -spectrum below the endpoint . The upgrade was completed by the application of HF pulses on one of the electrodes in between measurements every 20 s, and a full automation of the apparatus and remote control. The former improvement lowered and stabilized the background, the latter one allowed long-term measurements.
Figure 11 shows the endpoint region of the Mainz phase II data (from 1998, 1999, and 2001) in comparison with the former Mainz 1994 data. An improvement of the signal-to-background ratio by a factor 10 by the upgrade as well as a significant enhancement of the statistical quality of the data by long-term measurements are clearly visible. The main systematic uncertainties of the Mainz experiment are the inelastic scattering of -electrons within the tritium film, the excitation of neighbour molecules due to sudden change of the nuclear charge during -decay, and the self-charging of the tritium film as a consequence of its radioactivity. As a result of detailed investigations in Mainz [64, 76–78]—mostly by dedicated experiments—the systematic corrections became much better understood and their uncertainties were reduced significantly. The high-statistics Mainz phase II data from 1998–2001 allowed the first determination of the probability of the neighbour excitation to occur in of all -decays  in good agreement with the theoretical expectation .
The analysis of the last 70 eV below the endpoint of the phase II data resulted in  which—using the Feldman-Cousin method —corresponds to an upper limit of An analysis of the Mainz phase II data with respect to setting a limit on the contribution of a light sterile neutrino is underway .
3.3. The Troitsk Neutrino Mass Experiment
The windowless gaseous tritium source of the Troitsk experiment  is essentially a tube of 5 cm diameter filled with resulting in a column density of molecules/cm2. The source is connected to the ultrahigh vacuum of the MAC-E-Filter type spectrometer by a series a differential pumping stations (see Figure 12).
From their first measurement in 1994 the Troitsk group had reported for more than a decade the observation of a small, but significant anomaly in its experimental spectra starting a few eV below the -endpoint . This anomaly appeared as a sharp step of the count rate . Because of the integrating property of the MAC-E-Filter, this step should correspond to a narrow line in the primary spectrum with a relative intensity of about of the total decay rate. In 1998, the Troitsk group claimed that the position of this line was oscillating with a frequency of 0.5 years between 5 eV and 15 eV below . By 2000 the anomaly no longer followed the 0.5 year periodicity, but still existed in most data sets. The reason for such an anomaly with these features was not clear but gave rise to many speculations. In presence of this problem, the Troitsk experiment corrected for this anomaly by fitting an additional line to the -spectrum run-by-run.
In 2011 the Troitsk group repeated the analysis of their data . Special care was taken for calculating the time-dependent column density of the windowless tritium source and applying these values to the analysis. Secondly the data were very carefully selected with regard to data quality and stability of the experiment. Thus a time and intensity varying anomaly was not any more needed to describe the Troitsk -spectra. Combining all the selected Troitsk data from 1994 to 2004 gave  from which using the Feldman-Cousin method  an upper limit can be deduced
Due to the complicated electronic structure of 187Re and its -decay (compare to Section 2) the advantage of the 7 times lower endpoint energy of 187Re with respect to tritium can only be exploited if the -spectrometer measures the entire released energy, except that of the neutrino. This situation requires a -spectrometer, which at the same time contains the -electron emitting 187Re. The advantage of this approach is that many systematic uncertainties—like the electronic final state spectrum, energy losses by inelastic scattering, and so forth—drop out in first order, since all released energy except that of the neutrino is measured in the same way and summed up automatically (we assume that all deexcitation processes are faster than the integration time of a detector signal).
To reach the required energy resolution such a source = detector arrangement can be ideally realized by a cryobolometer (see Figure 13) . The energy release by the -decay results in a temperature rise of the crystal. This temperature increase can be measured by the thermometer, if is large enough. The temperature rise depends on the energy release and on the heat capacity In order to obtain a large temperature increase in the presence of a very small energy release (it should be noted, that the sensitivity on the neutrino mass requires a very low endpoint energy ), the heat capacity has to be extremely small. The first measure is using tiny cryobolometers of typical masses of mg. Secondly, the temperature of the cryobolometer has to be as low as possible. The Debye model states that the phonon part of the heat capacitance of a crystal consisting of atoms scales with the third power on temperature in units of the Debye temperature
There is a second reason why the crystals should be small: The cryobolometer is not an integral spectrometer like the MAC-E-Filter but measures always the entire spectrum. Therefore pileup of two random events may pollute the endpoint region of a -decay on which the neutrino mass is imprinted. The pileup rate of a detector of random rate requiring a minimal time interval to distinguish two events of amounts to . With a half-life y and a natural abundance of 62.6% the specific activity of pure rhenium amounts to about 1 Bq mg−1. The rise time typically scales with the mass of the cryobolometers. Rise times of cryobolometers with temperature read-out of mg mass are typically in the μs range (new cryobolometer approaches like Metallic Magnetic Calorimeters (MMC) as discussed in Section 5.4 could be much faster). Therefore cryobolometer detectors with mg masses are required to suppress pileup by 4 or more orders of magnitude. This has to be compared to the fraction of the -spectrum which contains the information on the neutrino mass. Even for the lowest known endpoint energy of 187Re with keV, the relative fraction of 187Re -decay events in the last eV below is of order only (compare to Figure 3). Therefore arrays of many bolometers are required to reach a high sensitivity on the neutrino mass.
Another technical challenge is the energy resolution of the cryobolometers. Although cryogenic bolometers with an energy resolution of a few eV have been produced with other absorbers, this resolution has not yet been achieved with rhenium.
3.5. 187Re -Decay Experiments
Two groups have started the field of 187Re -decay experiments at Milan (MiBeta) and Genoa (MANU): The MANU experiment used a single metallic rhenium crystal of about 1.6 mg as absorber read out by a neutron transmutation doped (NTD) germanium thermistor . Sensor and absorber were connected using epoxy glue and were suspended by aluminum wires used for readout of the device. These lines also provided the thermal link to a heat bath cooled down to temperatures below 100 mK using a 3He-4He dilution refrigerator. At a detector threshold of 350 eV, an event rate of about 1.1 counts per second was observed with this absorber. The rise time of this early detector was of the order of 1 ms with a decay time of the signals of several tens of ms. This clearly limits the amount of activity allowed per absorber crystal, in order to avoid pileup problems. The energy resolution of the detector was determined to be 96 eV FWHM. While the main objective of the initial measurement was a determination of the endpoint energy and the half life of 187Re -decay, an interesting side effect was the first observation of the so-called -environmental fine structure (BEFS) in pure rhenium . This fine structure is caused by an interference of the outgoing direct wave of the beta-electron with incoming waves reflected by the neighboring atoms and leads to a modulation of the shape of the -spectrum that is most pronounced at low electron energies. (see Figure 14(a)). This effect is of concern to future calorimetric neutrino mass experiments, as it produces difficulties to calculate distortions of the 187Re -spectrum which will become important for experiments aiming at a sensitivities in the sub-eV region . The effect was later also observed in the MiBeta experiment , albeit with a different characteristic due to the use of a different absorber material. In the subsequent analysis the data taken by MANU set an upper limit on the neutrino mass of eV/c2 at 95% CL .
The second pioneering 187Re -decay experiment was set up by the MiBeta collaboration who were working with AgReO4 absorbers with a mass of about 0.25–0.30 mg each, read out by silicon implanted thermistors. The group was the first to work with an array of detectors to circumvent the problem of the low maximum activity allowed for the individual crystals. Their setup contained 10 detectors with an average energy resolution of 28.5 eV FWHM . The rise times of the detector signals were of the order 0.5 ms. Figure 14(b) shows a Kurie plot of the accumulated spectrum obtained during one year of data taking. The analysis of the spectrum near the endpoint resulted in an upper limit on the electron neutrino mass of eV at 90% CL.
4. The KATRIN Experiment
The Karlsruhe Tritium Neutrino (KATRIN) experiment is a next-generation direct neutrino mass experiment which is currently under construction by a large international collaboration with groups from Lawrence Berkeley National Laboratory (LBNL), Bonn University, Fulda University of Applied Sciences, Max-Planck-Institute for Nuclear Physics Heidelberg (MPIK), Karlsruhe Institute of Technology (KIT), Mainz University, Massachusetts Institute of Technology (MIT), Münster University, University of North Carolina at Chapel Hill, Academy of Sciences of the Czech Republic, University of Santa Barbara, University of Washington (UW) Seattle, Swansea University and Institute for Nuclear Research (INR) Troitsk, with associated groups from Aarhus University, University College London (UCL), and Federal University of Paraná. The experiment is housed at Tritium Laboratory Karlsruhe (TLK) at KIT's Campus North site. KATRIN has been designed to substantially increase the sensitivity of the Mainz and Troitsk forerunners, while employing the same general spectroscopic principles. It will push the MAC-E-filter technology and tritium process technology, as well as many other methods, to their technological limits to press forward into the sub-eV sensitivity regime of the averaged electron (anti)neutrino mass .
KATRIN has been designed to improve the sensitivity by a factor of 10 from the present value of 2 eV to 200 meV at 90% C.L. This increase in sensitivity by one order of magnitude will allow covering almost the entire region of quasi-degenerate neutrino masses (nuclear and particle physics motivation) and to directly probe the neutrino hot dark matter fraction in the universe down to a contribution of the total matter-energy budget of the universe (cosmological motivation) . As the experimental observable in -decay is the square of the averaged electron (anti)neutrino mass , this requires an improvement of the experimental precision in -spectroscopy by two orders of magnitude as compared to Mainz and Troitsk. This in turn requires significant, major improvements to key experimental parameters such as source activity, energy resolution, and background rate. In addition, it requires a much better control of systematic effects (by one order of magnitude). The systematic effects mainly arise from parameters related to the tritium source such as fluctuations of the source temperature or the injection pressure or source composition (such as different hydrogen isotopologues or ions), but other systematics such as fluctuations of the retarding potential have to be limited to very small values as well.
Extensive design work by the KATRIN Collaboration reported in  has revealed that the key statistical parameters to reach a neutrino mass sensitivity of 200 meV are an increase of the source strength by a factor of 100 and of the measurement time by a factor of 10. The required energy resolution is less critical, as the information on the neutrino mass is derived from an energy interval which typically extends to about 30 eV below the endpoint, so that an improvement by about a factor of four to a value of eV at KATRIN is sufficient. On the other hand, the background rate close to the endpoint has to be limited to a very low value of <10−2 cps (counts per second) to achieve the design sensitivity. This is identical to the background value obtained at the Mainz and Troitsk experiments, so that novel background reduction mechanisms have to be developed for the much larger KATRIN spectrometer section.
The 70 m long KATRIN setup has been designed along these general criteria (see Figure 15). The experimental configuration can be grouped into the following major components: (a) the Windowless Gaseous Tritium Source (WGTS), where electrons are produced per second by the -decay of molecular high-purity tritium gas at a temperature of 30 K, (b) an electron transport and tritium elimination section, comprising an active differential pumping (DPS) followed by a passive cryo-pumping section (CPS), where the tritium flow is reduced by more than 14 orders of magnitude (c) the electrostatic prespectrometer of MAC-E-Filter type, which offers the option to prefilter the low-energy part of the tritium -decay spectrum, (d) the large electrostatic main spectrometer of MAC-E-Filter type which represents the precision energy filter for electrons and which is, with its dimensions of 10 m diameter and 24 m length, the largest UHV recipient in the world (not shown here are its inner electrode system and its outer air coil system) and (e) a segmented Si-PIN diode array to count the transmitted electrons. The experiment is completed by a rear section which allows to control and monitor key source-related parameters as well as the monitor spectrometer for redundant monitoring of the retarding potential. Finally, there are extensive infrastructure facilities.
The unique properties of the gaseous molecular tritium source (high activity and stability) and of the large main spectrometer (ppm precision of the retarding high voltage) allow KATRIN to extent its physics reach from its main goal of measuring the neutrino mass in the sub-eV range to look for contributions by possible sterile neutrinos from the sub-eV up to the multi-keV range. The sensitivity of KATRIN for (sub-)eV sterile neutrinos, as suggested by [87–90] is expected to be very high [91–93]. Interestingly, a high precision search of Warm Dark Matter in form of keV sterile neutrinos looks very attractive, as a heavy sterile neutrino would manifest itself as a tiny kink and subsequent spectral distortion deep in the -spectrum further away from the endpoint. Finally, KATRIN will allow performing stringent tests of other physics issues beyond the Standard Model such as large extra dimensions , right handed currents  and Lorentz violation .
4.2. KASSIOPEIA: A Full Monte Carlo Simulation Software for KATRIN
The main principles of the MAC-E filter and its application in the KATRIN experiment can be understood analytically via the adiabatic approximation. However, in order to understand the behavior of electrons not only in perfect vacuum and beyond the assumption of adiabatic motion, a precise and fast computational tool is required. Over the past years the KATRIN Collaboration has developed the universal code package KASSIOPEIA , which is based on the ground-laying work described in [97, 98]. The software package includes various particle generators for signal electrons as well as background processes.
KASSIOPEIA comprises a number of modules for the creation of particles, the calculation of their subsequent trajectories in electromagnetic fields, and the detection of particles in Si-based detectors. To do so, KASSIOPEIA provides a detailed WGTS model, a number of electric and magnetic field calculation methods and different methods for the calculation of particle trajectories. The particle detection module includes backscattering of electrons on the detector surface and a comprehensive number of physical phenomena of low-energy electrons in silicon .
At KATRIN very large dimensions of the order of tens of meters (main spectrometer) concur with very small dimensions of the order of μm (e.g., wire electrode). This fact constitutes the biggest challenge for the electromagnetic field calculation, that cannot be handled adequately by commercial programs. A fast and precise field calculation is achieved by a variety of field calculation methods, ranging from very fast axisymmetric field calculations [97, 98], to fully three-dimensional field calculation methods  based on the boundary element methods which is most suitable for large dimensional differences.
The particle trajectory calculations are based on explicit Runge-Kutta methods described in [101–103] and include physical processes like synchrotron radiation and elastic, electronic excitation and ionization collisions of electrons with molecular hydrogen [104–106].
The detailed tritium source model allows to simulate the actual neutrino mass measurement with high precision. The source model includes, among other things, the final state distribution of tritium [35, 36], the Doppler broadening, and scattering in the source [107, 108].
KASSIOPEIA was successfully used to study background processes  and transmission properties of the spectrometer , to optimize the electromagnetic design of the spectrometer section and to study systematic effects related to the source section of KATRIN .
4.3. Source and Transport Section
The tritium source-related parts of KATRIN comprise the Windowless Gaseous Tritium Source (WGTS), the Differential (DPS) and Cryogenic (CPS) Pumping Sections, as well as the Calibration and Monitoring System (CMS) at the rear end. The main tasks of these elements are to (a) control and monitor the tritium column density in the source tube to a precision of better than and (b) reduce the tritium partial pressure from the source inlet to the entrance into the spectrometers by more than 14 orders of magnitude, while at the same time transporting the -decay electrons adiabatically towards the spectrometers.
The extensive tritium infrastructure required for the continuous operation of the KATRIN tritium source is provided by Tritium Laboratory Karlsruhe (TLK), which has developed key technologies for experiments involving large amounts of tritium since the early 1990s. As of today, TLK has an amount 24 g of tritium on-site and holds a license for up to 40 g. Tritium process technology at TLK is based on closed tritium cycles with their central elements of tritium storage, isotope recovery and separation, and tritium retention, among others. The yearly throughput of 10 kg of high purity tritium in the KATRIN source will be equivalent to the ITER operation in the late 2020s . The KATRIN tritium cycle is subdivided into two closed loops, with an Inner Loop being optimized for a highly stabilized injection pressure into the WGTS, and an Outer Loop designed to maintain a high tritium purity by withdrawing a small tritium fraction (<1%) for clean-up and isotope separation (the amount withdrawn is reinjected to the Inner Loop after purification).
4.3.1. Gas Dynamics of the WGTS
The geometry of the tritium source of KATRIN is defined by a cylindrical stainless steel beam tube of 90 mm diameter and 10 m length, which is housed inside the WGTS cryostat. At the center of the beam tube, high purity molecular tritium gas is injected via capillaries with an inlet pressure of about mbar. At both ends of the beam tube the injected flow is pumped out by large turbomolecular pumps (TMP). The isotopic composition of the injected hydrogen isotopologues (, , , HT, DT, HD), with dominating by a large factor, is constantly being monitored by Laser Raman spectroscopy. In this windowless geometry the systematic effects for electrons close to the 18.6 keV end point of tritium -decay are minimized, allowing to adiabatically transport electrons out of the source by a homogeneous magnetic field of T. This field is provided by a system of three large superconducting solenoids surrounding the beam tube. The operating temperature of the beam tube is in the 27–30 K temperature range to minimize the tritium throughput through the beam tube as well as the contribution of thermal Doppler broadening of electron energies due to the molecular motion. Finally, this temperature regime also mitigates systematic effects due to clustering and condensation of hydrogen isotopologues.
Over the past years, a detailed 2D/3D model of the gas dynamical characteristics of the tritium source has been developed, which is much more detailed than the initial 1D model described in , where perfect temperature homogeneity and isotropy of the source was assumed. The new gas dynamical model of the source has been used to derive precise values for the integral WGTS column density and its radial variation and the pressure distribution along the beam tube . The stability of the column density, which is of the order of -molecules/cm2, plays a central role in the sensitivity of KATRIN, as it impacts not only the statistical accuracy of the measurement, but also defines the energy losses of the -decay electrons inside the source due to inelastic scattering . An important aspect in the detailed modelling of the source characteristics is the fact that the gas rarefaction parameter along the tube varies from the strong hydrodynamic regime at the injection point to free molecular flow at the beginning of the differential pumping chambers. Accordingly, a significant part of the beam tube exhibits values in the intermediate regime, where the phenomenological intermediate conductance formula of Knudsen is applied .
4.3.2. The WGTS Cryostat and Associated Control and Monitoring Systems
Initial estimates in  have revealed that a stabilisation of the column density on the 0.1% level is required to obtain a neutrino mass sensitivity of 200 meV. This translates into a stabilisation of the inlet pressure and the beam tube temperature on the same level, which certainly is a major technological challenge, given the macroscopic WGTS beam tube dimensions. A stable gas feed into the WGTS is achieved by injecting high-purity tritium gas from a pressure controlled buffer vessel over a 5 m long tritium capillary with constant conductance, while pumping out the gas at both ends of the beam tube with large TMPs. To fulfil the requirements on temperature stability, a novel beam tube cooling system based on a 2-phase boiling fluid has been designed and tested successfully (a two phase system, when used as a cooling method, has the well-known distinct advantage that it can absorb heat without changing its temperature). The choice of neon as cooling fluid has the advantage that it coexists in two phases at temperatures around 30 K for moderate pressures in the range from 1–3 bar.
In combination with the requirements of adiabatic electron transport and safe tritium throughput in a closed loop system, the required stability of the tritium column density on the level has resulted in a technical design of a complex cryostat (see Figure 16) with a length of 16.1 m and a weight of 27 tons  and a cryogenic system comprising 13 fluid circuits where 6 cryogenic fluids are processed. The cryostat also houses an extensive array of more than 500 sensors (magnetic field, pressure, gas flow, liquid levels, voltage taps) as part of a dedicated measurement and control system. The precision temperature readings of the sensor array will later on be used as input to model minute fluctuations of the actual source characteristics over typical time scales of about 1 min.
In order to test the advanced cryogenic technologies required for the cooling of the beam tube, a large-scale test unit, the WGTS demonstrator, was built. This 10 m long (1 : 1 scale) unit makes use of many original components which will be used later on in the final WGTS cryostat. In 2011 extensive measurements with the WGTS demonstrator were performed at TLK targeted at a proof-of-principle and optimization of the beam tube cooling system. To obtain a precise measurement of the beam tube temperature, an array of 24 metallic resistance temperature sensors (Pt500) were distributed along the beam tube. The sensors feature a sub-mK resolution and, via calibration with sensors measuring the saturation pressure of neon , an uncertainty for the absolute temperature of 4 mK only. During the measurements, a peak to peak variation of mK over a time period of 4 h at K was obtained, which is one order of magnitude better than required [107, 116]. This result offers the potential to substantially reduce the systematic errors from the tritium source during the long-term measurements. At present, the WGTS demonstrator is being reassembled to the final WGTS cryostat. These extensive works will be completed by the end of 2014.
A dedicated Laser Raman (LARA) system, based on inelastic Raman scattering  of photons from gas molecules, is connected to the inner tritium loop and allows for inline and near time monitoring of the gas composition. The different tritiated gas species (T2, DT, HT) have to be distinguished not only because of their different properties in the gas dynamical calculations, but also due to their slightly different endpoint energies as well as rotational and vibrational final states. Test measurements with a Raman cell suitable for operation with tritium [118, 119] show that a precision of for monitoring can be achieved in measurement intervals of <250 s under KATRIN operating conditions [120, 121]. With an optimized setup, the acquisition time recently has been reduced to 60 s only.
4.3.3. Calibration and Monitoring System Upstream of the WGTS
To monitor the stability of the tritium activity and of the column density of the windowless gaseous tritium source, the rear side of the WGTS is connected to a calibration and monitoring system. An angular-selective electron gun is used to measure the energy loss function of the electrons due to inelastic scattering in the WGTS and to monitor the integral column density at different source radii. This novel UV-laser-based e-gun is based on the photoelectric effect producing monoenergetic electrons with stable intensity and with a well-defined starting angle with respect to the magnetic -field [122–125]. Secondly this electron source will be used to determine the transmission function of the main spectrometer and the response function of the whole setup with high precision.
The calibration and monitoring system has another very important task related to the source operation. The -electrons will be guided from the source to the spectrometer and mapped onto the detector within a magnetic flux tube of 191 T cm2. The guiding by superconducting solenoids implies a tight transverse confinement by the Lorentz force to all charged particles. This includes the daughter ions per second, which emerge from -decay in the source tube, as well as the electron-ion pairs per second produced therein by the -electron-flux through ionization of T2 molecules. The strong magnetic field of 3.6 T within the source is confining this plasma strictly in the transverse direction so that charged particles cannot diffuse to the conducting wall of the source tube where they would get neutralized . The charges within the windowless gaseous tritium source will thus be neutralized by a gold-plated crystalline rear wall which is part of the calibration and monitoring system, and on which all magnetic field lines from the windowless gaseous tritium source end. Accordingly, the high longitudinal conductance of the plasma will define and also stabilise the electric potential within the tritium source.
Finally, the source activity will be monitored by -induced X-ray spectroscopy (BIXS), where the potential-defining rear wall acts as an X-ray converter. An encapsulated X-ray detector positioned behind the rear wall will monitor the source activity by recording the X-ray intensity, while the fluorescence lines from the surrounding structures would allow an inline calibration of the setup.
4.3.4. Tritium Retention Techniques
Since the spectrometer section must be essentially tritium-free for background reasons (see next chapter), the tritium flow along the beam line must be reduced from its initial value at the injection point of 1.8 mbar·l/s to a level of mbar·l/s at the entry to the spectrometers. This unprecedented large suppression factor will be achieved by a combination of differential and cyrogenic pumping [127–130]. While reducing the tritium flux by 14 orders of magnitude, these sections also have to maintain adiabatic guidance of -decay electrons from the WGTS to the spectrometers over a distance of more than 15 m.
4.3.5. Differential Pumping Fection (DPS2-F)
Differential pumping by large TMPs is the first tritium retention technique in use at the KATRIN beam line, with the initial stages being performed inside the WGTS cryostat by pump ports at the rear (DPS1-R) and front (DPS1-F) section. Thus, the DPS2-F cryostat is the second pumping unit in the forward direction. Its scientific objectives are the following: (a) active pumping of tritium molecules with TMPs, (b) reduction of the tritium flow by a factor larger than , (c) maintaining a stable tritium circulation in the Outer Loop, and (d) diagnostics of the composition and suppression of ion species from the WGTS.
The DPS2-F has been designed as a single large cryostat of 6.96 m length (see Figure 17). It was the first KATRIN source-related component on site at KIT and was successfully commissioned in 2010/11 . The unit features a beam line with an -shaped chicane to avoid the molecular beaming effect, in which case gas particles would not be pumped out due to the alignment between gas particle momenta and beam tube axis. The cascaded differential pumping system consists of four large TMPs mounted on pump ducts, which minimize the effect of thermal radiation from the TMPs operated at room temperature to the beam tube kept at 77 K, while also limiting TMP rotor heating due to induced eddy-currents . During the commissioning measurements a retention factor for tritium (based on measurements for other hydrogen-like and noble gas species) of ≈2.5 × 104 was obtained, demonstrating the excellent tritium retention characteristics of differential pumping already at room temperature .
The DPS2-F beam tube instrumentation is used for diagnostics and suppression of ions which result from ionization processes of primary -decay-electrons with the neutral gas molecules of the source. As ions are guided by the field lines (and thus are not being pumped out) they remain trapped between a ring electrode at the end of the DPS2-F, where they are reflected electrostatically, and the WGTS, where they are reflected back by the gas pressure. The number density and composition of ions will be measured on a near-time basis by two Fourier Transform-Ion Cyclotron Resonance (FT-ICR) devices , which discriminate Penning-trapped ions by detecting differences in Ion Cyclotron frequencies in the 5.6 T -field at the center of the DPS2-F solenoids. The ions stored between the WGTS and the end of the DPS2-F are eliminated by a system of three electrostatic dipoles by an induced -drift, irrespective of their mass, charge sign, or motion with respect to the electric field, within a time period of 20 ms or less, which is much shorter than the allowed upper time limit of 10 s for ion storage.
The magnetic flux tube of 191 T cm2 of -decay electrons which is used for the neutrino mass analysis is guided through the beam tube chicane with a minimum margin of 10 mm from the walls of the inner vacuum system. The system of 15 superconducting solenoids is cooled by liquid helium supplied from an external Linde TCF 50 refrigerator, with a measured DPS2-F cryogenic stand-alone time of about eight hours . Unfortunately, due to a malfunction of one of the protective diodes of the magnet system during the commissioning measurements, further tests with the DPS2-F had to be stopped in mid-2011. At present different options to ensure a future fail-safe operation of the differential pumping system during long-term tritium running are being studied.
4.3.6. Cryogenic Pumping Section (CPS)
The Cryogenic Pumping Section (CPS) is the final tritium retention unit in the KATRIN beam line, featuring a large-area cryosorption pump based on a thin layer of argon snow frosted onto gold-plated steel as adsorbent (see Figure 18) ). It is designed for a tritium flow reduction factor of ≥3 × 107.
The 6.73 m long CPS cryostat features an even more complex beam tube design than DPS2-F with a beam tube segmented into 7 sections, as shown in Figure 18. To prevent residual tritium absorption onto the walls, the inner surfaces of the beam tube elements are gold-plated. The beam tubes of the cryopump section are tilted by an angle of 15° with respect to each other and are covered by a thin argon frost layer to enlarge the inner surface for sorption. The CPS design allows the insertion of a Condensed Source (CKrS)  for calibration purposes, as well as a small Si-diode for source intensity monitoring. The CPS is currently being manufactured by an industrial partner in Genoa (Italy) and will be shipped to KIT at the end of 2013. After commissioning and initial bake-out, the CPS will regularly cycle through several modes of operation. In a first step the beam tube elements will be cooled to their respective operational temperature, which is followed by the preparation of an argon frost layer. This requires a temperature in the range of 6 K to favour the growth of small to medium sized crystallites. After argon frost preparation, the cryopump temperature will be lowered again to 3 K. After a data collection cycle of about 60 days, the CPS has to be regenerated to remove the adsorbed tritium activity of about 2 Ci by purging 100 K gaseous helium through the cryopumping segments.
By combining differential pumping by TMPs in the DPS with cryosorption on a large-area argon frost pump in the CPS, the number of migrating HT-molecules to the pre- and main spectrometer is kept at a minimum. When including all relevant tritium retention techniques, that is, not only DPS and CPS, but also tritium interactions with the spectrometer walls and tritium pumping by the large nonevaporable getter (NEG) strips in the spectrometers, a tritium partial pressure in the main spectrometer of about mbar can be expected . This ultralow level should result in a tritium-related background rate of only cps, which is two orders of magnitude below the design specification. It is only by this huge tritium reduction factor that the KATRIN experiment can combine an ultra-luminous windowless gaseous tritium source of -decay electrons with a high-resolution electrostatic spectrometer operated with an extremely low background level of 0.01 cps .
4.4. Spectrometer and Detector Section
In the following the non-tritium-related parts of KATRIN will be described. First the tandem spectrometer system performing the energy analysis of the -electrons is addressed and then the detector and data acquisition system are described.
4.4.1. Overview: Prefiltering and Precision Scanning
The spectrometer section consists of two electrostatic retarding filters of the MAC-E-Filter type: the prespectrometer and the main spectrometer. The prespectrometer offers the option to act as a prefilter, reflecting all electrons with energies up to 300 eV below the endpoint, while transmitting the interesting part of the spectrum undisturbed. All electrons transmitted through this first stage are guided to the main spectrometer for precise energy analysis close to the endpoint. A third spectrometer, the monitor spectrometer, is used to monitor precisely the high voltage of the main spectrometer.
The prespectrometer has a length of 3.4 m and a diameter of 1.7 m. At both ends a superconducting magnet is installed providing a magnetic field of 4.5 T at the center of the magnet and 15.6 mT in the center of the spectrometer. As a novel electromagnetic design feature, if compared to the Mainz and Troitsk set-ups, the tank itself is set on high voltage. An inner electrode system consisting of two conical full electrodes and a central cylindrical wire electrode can be set to a different potential as the vessel. Figure 19 shows a schematic view of the prespectrometer.
One of the prespectrometer’s major tasks has been to serve as a prototype for advanced technologies and experimental methods later applied to the much larger main spectrometer. Many basic concepts, such as the design of the ultra high vacuum (UHV) system and high voltage stabilization were successfully tested at the prespectrometer. Especially important was the detailed investigation of background processes in MAC-E-filters. Two new classes of background were identified and studied at the prespectrometer: background due to small Penning traps and radon induced background . The knowledge gained during the extensive measurement period at the prespectrometer has proven to be extremely valuable for the design and operation of the main spectrometer.
4.4.3. Main Spectrometer
The exceedingly large dimensions of the main spectrometer are essential to operate it as an extremely precise high-energy filter with an energy resolution of eV (0%–100% transmission). This sharp resolution requires that the stability of the retarding high voltage is in the ppm range. The scanning voltage is varied in steps of 0.5–1 V around a narrow region close to the endpoint at keV. The highest electrostatic potential is located in the central plane of the spectrometer, perpendicular to the beam axis, the so-called “analyzing plane.”
To achieve a very high energy resolution, the cyclotron motion of the electrons, being isotropically created in the WGTS, needs to be almost fully transformed into longitudinal motion parallel to the magnetic field lines, since only the latter is analyzed by the retarding potential (see Section 3.2). To obtain adequate parallelization, the magnetic field drops by four orders of magnitude from the entrance to the center of the main spectrometer. Since the magnetic flux is conserved, the cross-section of the flux tube in the center plane has to be four orders of magnitude larger than at the entrance. This scaling explains the huge size of the main spectrometer (length m, diameter m, surface area m2, and volume m3) (see Figure 20).
To compensate for the distorting earth magnetic field and to fine-tune the magnetic field inside the main spectrometer, the vessel is surrounded by a large external air coil system (see Figure 21). It consists of 10 horizontal current loops and 16 vertical ones to compensate the earth magnetic field (EMCS). Additionally, a system of 15 Helmholtz-like coils with individually adjustable currents allows for precise adjusting of the gradients and the overall strength of the magnetic field (LFCS). Two arrays of magnetic field sensors, one at fixed positions, and the other a mobile one attached to robots , allow to measure and monitor precisely the magnetic field.
The main background source of the Mainz experiment was secondary electrons from the electrodes at high potential. These electrons are created by cosmic muons or by environmental radioactivity. The shielding effect of the magnetic field will prohibit most of these electrons to enter the magnetic flux tube and thus to reach the detector. However, measurements at Mainz  and Troitsk yielded a transmission rate of through the magnetic shielding of a MAC-E-filter. A large shielding factor is important considering the large surface of 650 m2 of the KATRIN main spectrometer. To suppress this kind of background, a new idea of additional background suppression has been developed and successfully tested at the Mainz spectrometer [142, 143] and optimized at the prespectrometer. This technique will also be applied to the KATRIN main spectrometer: The vessel walls at high potential will be covered by a system of nearly massless wire electrodes, which are put to a slightly more negative potential (see Figure 22). Secondary electrons ejected by cosmic rays or environmental radioactivity from the vessel wall will thus be repelled and prohibited from entering the magnetic flux tube. To achieve a high suppression factor of 10, as demonstrated in the Mainz experiment, the KATRIN main spectrometer has been instrumented by a two-layer wire electrode system [144–146] (see Figure 23).
This system consists of 248 modules of a typical size of m2. The wires of the inner (outer) layer have a diameter of 200 μm (300 μm) and a distance of 22 cm (15 cm) from the surface of the spectrometer vessel. Essentially the inner (outer) wire layer will be on a potential of −200 V (−100 V) with respect to the main spectrometer vessel. The total number of wires amounts to 23440. Stringent demands are required for this electrode system: The two wire layers in the central part have to be mounted with a precision of 200 μm to avoid distortions of the electrical potential. The wires have to withstand baking at C. The outgassing of the wire electrode system has to be extremely low to allow a final residual gas pressure of mbar. The mounting of the wire electrode system within the huge KATRIN main spectrometer under clean-room conditions was a big challenge on its own (see Figure 23). The installation of the wire electrode system was successfully completed in January 2012.
In addition to the background reduction, the wire electrode serves three other purposes: fine-shaping of the electric retarding potential is achieved by applying slightly different voltages on the different axial rings of the wire modules to guarantee full adiabatic motion of the electrons. Reducing fluctuations of the electric retarding potential is done by decoupling the low-noise high-precision retarding voltage of the inner-most wire layer from the high-voltage of the vessel, on which electrical devices such as TMPs are mounted. The multilayer system acts as a Faraday cage. Finally, by splitting the wire module arrangement into an eastern and a western half an electric dipole voltage can be applied to eject low-energy stored particles by an induced drift (see, (42)).
4.4.4. High Voltage and Monitoring Layout
The energy analysis of -decay electrons based on the MAC-E filter principle [70, 71] relies primarily on the stability of the electrostatic filter potential [147, 148]. The latter has been identified as one of the five main contributions to the KATRIN uncertainty budget . Its contribution to the systematic error of the energy filter potential has to be restricted to eV2 in order to meet the desired level of sensitivity.
We showed earlier in (37) that any unknown fluctuation of the retarding voltage with variance leads to a shift of the squared neutrino rest mass (the electron charge takes care of the conversion from voltage to energy.). This relation and the systematic error limit given above have been taken into account in order to define the stability requirements for the KATRIN high voltage monitoring system, which is mV . At the tritium endpoint energy keV, this corresponds to a relative voltage stability of , which has to be maintained for the whole KATRIN measurement time.
Therefore KATRIN will apply two redundant monitoring methods, one being based on two high voltage dividers with highest possible precision, that is, relative precision and long-term stability for up to 35 kV , and the other one being based on a monitoring beam-line which observes a monoenergetic electron source . To do so in both cases, the low-noise retarding high voltage applied to the inner layer of the central wire electrode modules is distributed to the precision high voltage dividers and to the monitor spectrometer.
The monitor spectrometer consists of the refurbished and upgraded MAC-E filter spectrometer of the Mainz experiment (see Section 3.2) and observes an implanted source [150, 151]. For the calibration source is produced either at the ISOLDE facility at CERN or at the cyclotrons at Nuclear Research Institute (Rez) in the Czech Republic or at Bonn University in Germany [152, 153]. The decay provides monoenergetic conversion electrons of the K shell with an energy close to the tritium endpoint at keV . While applying the MAC-E filter principle as well as the main spectrometer retarding potential, the source has to be elevated by about 800 eV in order to monitor the conversion electrons. This method relies on atomic and nuclear physics standards and provides monitoring with a stability at the 50 meV level as required by the systematic uncertainty limit.
Based on a self-compensating principle , two precision high voltage dividers have been developed in cooperation with PTB Braunschweig. These HV dividers rank among the most stable ones in the world. For voltages of up to 35 kV, they provide a long-term stability of <10−6 per month as well as negligible temperature and voltage dependencies . The two dividers have served recently to calibrate the ISOLDE facility at CERN . The permanent voltage monitoring with the precision voltage divider and frequent calibration measurements at the monitor spectrometer will provide reliable and stable calibration and monitoring for the retarding potential of the main spectrometer over long periods of time.
4.4.5. Focal Plane Detector
The focal plane detector is a semi-conductor-based silicon PIN diode . Its main goal is to detect transmitted electrons with a detection efficiency of >90%.
The detector is subdivided into 148 pixels to achieve good spatial resolution (see Figure 24). This is important, as electrons passing the analyzing plane at different radii will experience slightly different retarding potentials, due to the radial inhomogeneity of the electric potential. To account for this effect the detector consists of 12 concentric rings subdivided azimuthally into 12 pixels each. The innermost part is composed of 4 quarter circle pixels. This pixel arrangement allows for a precise mapping of inhomogeneities of the retarding potential as well as of the magnetic field. Each detector pixel measures an independent tritium -spectrum which has to be corrected for the actual retarding potential and magnetic field value at the analyzing plane.
The signal of each detector pixel is connected by pogo-pins to a charge-sensitive pre-amplifier. The amplified signals are digitized by FADC- and FPGA-boards  developed for the Pierre Auger Observatory. The Data Acquisition System is based on the software package ORCA (Object-oriented Real-time Control and Acquisition), which was originally developed at the University of Washington as a general purpose, highly modular, object-oriented, acquisition and control system. Since 2008 further development continues at the University of North Carolina at Chapel Hill .
The detector is cooled by . To suppress external background, it is surrounded by low radioactivity lead shield. Additionally an active muon veto system consisting of a cylindrical plastic scintillator read out by wavelength shifting fibers is installed. There is an option to use postacceleration to increase the kinetic energy of the -electrons to about 30 keV or above before they hit the detector. This could help in discriminating signal from fluorescence background. Additionally, the impact angle of the electrons relative to the detector surface would be increased, which decreases the probability of backscattering from the detector surface .
The detector is situated in a superconducting magnet providing a magnetic field of 3–6 T. The detector magnet is adjacent to the so-called pinch magnet which provides the maximal magnetic field of T of the entire KATRIN setup. All electrons that started in the source with an angle larger than will be reflected by the pinch magnet before they reach the detector (see (40)). This is advantageous, as electrons emitted under a large angle perform a lot of cyclotron motion, which in turn increases their total path length and therefore their scattering probability and also their synchrotron losses. To exclude those electrons, the magnetic field of the source is smaller than the maximal field.
4.5. Signal and Sensitivity
By counting -electrons at different retarding potentials, KATRIN measures the integral -spectrum (18) which is a convolution of the response function , incorporating the energy loss in the WGTS and the spectrometer transmission function, with the differential energy spectrum , describing the number of decays per second, nucleus and energy bin. In contrast to (18), here we have explicitly stated the dependence of on the endpoint energy and on the squared neutrino mass . The parameters and denote the total number of tritium nuclei and the measurement time at a certain retarding potential, respectively. Assuming a constant background rate the fit function to the measured spectrum is given by where and are the relative fraction of signal and background. In the fit , , , and are free parameters.
4.5.1. Sources of Systematic Errors
Sources of systematic uncertainties arise from(i)unconsidered corrections to the -spectrum ,(ii)unaccounted variations of experimental parameters, for example retarding potential and number of decays , (iii)an imprecise knowledge of the response function , (iv)and a non-constant background in time and energy.
Most of the systematic effects are related to the tritium source. They include energy losses of -electrons, primarily due to inelastic scattering in the source, fluctuations of the column density, potential charging effects due to remaining ions, and the accuracy of the quantum-chemical calculations of the final state distribution of molecular tritium (see Section 4.3). The main systematic effect related to the spectrometer section is the high voltage stability (see Section 4.4).
New physics effects such as right handed currents , sterile neutrinos [91–93, 160], extra dimensions , and so forth. also affect the shape of the tritium spectrum (see (21)). However, the influence of these phenomena typically becomes important only further away from the endpoint. Another new physics contribution that might affect the spectrum also close to the endpoint is violation of Lorentz invariance .
4.5.2. Background Sources
Only a small fraction of of the tritium decays per second in the WGTS will produce -electrons in the interesting energy region close to the endpoint (1 eV below the endpoint). This leads to a generic rather low count rate of only cps in this energy region. To achieve a sensitivity of 200 meV (90% C.L.), the background must be of the same order of magnitude, or smaller. Of major concern in this context are backgrounds which are not constant in time or energy.
The background originates mainly from the spectrometer section and partly from the detector. With a detector energy resolution of the order of 1 keV, all nonsignal electrons detected in the region-of-interest (i.e., from approximately 15 keV to 20 keV) will contribute to the background level. Possible detector background sources are electrons produced by cosmic muons (and subsequent neutrons and gammas), high energetic gammas of environmental radioactivity (mainly from the thorium and uranium decay chains in the surrounding area) and decays of radionuclei in the detector material.
To reduce the backgrounds, the materials used in the construction of the detector system were radioassayed very thoroughly. In addition, the detector is surrounded by a muon veto system. Finally, there is the option to use postacceleration to increase the energy of signal electrons up to 30 keV, which would allow for better signal-background discrimination. The total expected detector background is about cps, which has been verified in the on-going commissioning measurements of the focal plane detector system.
The spectrometer-related background is especially challenging since all low-energy electrons created in the spectrometer are accelerated to the retarding potential before hitting the detector. Consequently, they cannot by distinguished from signal electrons. In Section 4.4, we have already discussed that secondary electrons produced at the spectrometer walls are prevented from entering the sensitive magnetic flux tube by the dominant inherent magnetic shielding of the MAC-E filter, and by subdominant electrostatic shielding of the wire electrode system. However, secondary electrons can also be generated inside the sensitive flux tube as a result of ionisation processes of residual gas. Of particular concern are electrons that are stored, either in Penning traps, or by the inherent magnetic bottle property of the MAC-E-configuration. These processes have been studied in the prespectrometer set-up in an extended measurement program.
In [137, 138], it was shown that small Penning traps of the size of cm3 can lead to background rates easily exceeding cps. Electrons being stored in Penning traps produce background via messenger particles, such as positive ions or photons, which can leave the trap and subsequently can create secondaries in the sensitive volume. This background has been tackled by a very precise and careful electromagnetic design [161, 162] that avoids the creation of Penning traps ab initio.
The simultaneous operation of two electrostatic retarding spectrometers (prespectrometer and main spectrometer) next to each other in the beam line, and with ground potential between them, would create a large Penning trap for electrons . Signal electrons that experience no energy loss on their way through the spectrometers are not affected. However, all electrons that are created between the spectrometers with less kinetic energy than the prespectrometer retarding energy remain trapped. A stored electron will fill the trap in an avalanche effect by continually ionizing residual gas molecules. Simulations have shown that a primary stored electron together with all its secondaries may create up to positive ions, which are not trapped and can freely propagate into the main spectrometer. There, they can ionize residual hydrogen molecules thereby producing further electrons that can reach the detector and contribute to the background level.
A possible countermeasure would be the installation of a wire scanner between the spectrometers which regularly wipes through the trap, thereby removing the stored electrons and suppressing the background production [163, 164]. Another most promising option would be to operate the prespectrometer at reduced filter energy [99, 110], even down to zero potential. This mode would avoid the creation of the Penning trap. Test measurements and simulations have proven the effectiveness of both approaches.
Another major background source results from magnetically stored electrons [165, 166]. At the main spectrometer, this will cause all electrons with (transversal) energies larger than 1 eV (i.e., larger than the width of the transmission function) to be stored. Owing to the excellent UHV conditions in the KATRIN main spectrometer ( mbar) the storage times of electrons in the multi-keV-range can reach several hours. During its cooling time, a stored electron produces several hundred secondary electrons, which eventually leave the spectrometer and hit the detector. The main source of primary electrons originate from ,  and tritium decays inside the spectrometer volume. This background source features large time fluctuations and thereby a strong non-Poissonian characteristic. This fact has a large impact on the neutrino mass sensitivity .
As a counter measure liquid nitrogen cooled baffles will passively shield the spectrometer volume from , which is mainly emanating from the large NEG pumps . The tritium decay rate inside the main spectrometer volume can be reduced by the option to use the prespectrometer as an additional tritium pump. Finally, a number of active background reduction techniques, which are targeted at removing the stored electrons, have been successfully investigated experimentally and via extensive MC simulations. A method based on stochastic heating with an RF field by the well known technique of Electron Cyclotron Resonance (ECR)  and a method based on nulling the magnetic field by magnetic pulsing  seem most promising. Figure 25 shows the strong background-reducing effect of a very short RF-pulse fed into the main spectrometer on magnetically stored electrons as calculated by KASSIOPEIA (even for very moderate RF-amplitudes). These short RF-pulses would empty all stored particles with negligible duty-cycle.
4.5.3. Neutrino Mass Sensitivity of KATRIN
For the reference setup of the KATRIN experiment, the quadratic sum of all known systematic uncertainties is expected to be eV2. After three “full beam” years of measurement time the statistical error is about as small as the systematic error. This yields a total error of eV2. Figure 26 shows the discovery potential of KATRIN as a function of beam time for different neutrino masses. Accordingly, a neutrino mass of meV would be seen with significance. The figure also shows the 90% C.L. upper limit as a function of measurement time in case that no neutrino mass signal is seen. In this case, an upper limit of meV can be stated at 90% C.L after three full beam years of measurements.
4.6. Status and Outlook
The commissioning of the KATRIN beam line elements is progressing well. On one side, the multifaceted technical hurdles are challenging, causing the commissioning of the overall setup to take longer than originally planned. On the other side, the excellent results of the large-scale test units have opened the possibility to substantially reduce the systematic effects during the long-term measurements with the final configuration, thus improving the neutrino mass sensitivity of the experiment.
Beginning at the source-related components, an important breakthrough was the verification of the novel beam tube cooling system. In a dedicated setup, consisting largely of original components, the WGTS demonstrator, a temperature stabilization of the 10 m long beam tube of of ≈10−4 was achieved by using two-phase neon fluid as cooling agent. This is one order of magnitude better than specified. Together with the achieved pressure stabilisation of the inner loop mock-up of of ≈10−4, this opens up the possibility of reduced systematic errors from column density fluctuations (this is one of the largest overall systematic errors). At present the WGTS demonstrator is being reassembled to the final WGTS cryostat. The final WGTS assembly at KIT is expected to be completed by the end of 2014. Further progress has been made with regard to tritium analytics (LARA setup), as well as the design of the rear section which will include extensive control and monitoring units.
Major progress has also been achieved in the field of large-scale tritium retention. After the successful commissioning of the DPS2-F cryostat, first tritium retention measurements with a beam tube at room temperature have yielded experimental flow suppression factors which are in very good agreement with corresponding MC simulations. Due to the malfunction of a protective diode of the superconducting magnet system of DPS2-F, a new magnet safety concept for all s.c. solenoids has been designed. This concept is currently being implemented for WGTS and CPS, as well as a fail-safe differential pumping section. The manufacture of the cryopump CPS is well under way with assembly works expected to be finished by the end of 2013.
In the spectrometer section, the extensive measurement program with the prespectrometer facility has given important insights into background reduction techniques, precision electromagnetic layout, vacuum technologies, and high voltage stability. At present the prespectrometer is ready for beam line integration.
The main spectrometer together with its external air coil system and its inner electrode system, which was completed at the beginning of 2012, is currently being prepared for test measurements. These measurements will be focused first on extensive background studies, with the objective to remove any remaining small-scale Penning traps, to quantify the contribution of cosmic muon induced background and to study its signature by making use of external muon detectors. An important aspect of the background studies will be the identification of background due to stored electrons following nuclear decays, and the optimisation of active and passive background reduction techniques to limit the spectrometer background to a level of <10−2 cps. Another important task will be to map the transmission properties of the spectrometer with an angular-selective electron gun. In all these investigations the recently commissioned focal plane detector system with its excellent properties will be of vital importance. Finally, the extensive software developments for simulation and analysis tools are in an advanced state and the software packages are continually being refined and extended.
After integration of all source-related and spectrometer-related components, the first runs in the final KATRIN configuration are expected in the second half of 2015.
5. New Approaches
While spectrometer experiments based on the MAC-E filter principle  currently provide the highest sensitivities in direct neutrino mass experiments, there are alternative approaches that aim for comparable performance and better scalability in the study of weak decays.
A very recent development is promoted by the Project 8 team (see Section 5.1) where tritium technology from the KATRIN experiment is used in conjunction with microwave antennas to detect coherent cyclotron radiation emitted by individual decay electrons in a magnetic field. The aim is to extract a -decay spectrum without the need for a large electrostatic spectrometer.
Most of the work on alternative experimental methods is, however, focused on using microcalorimeters to study rhenium -decays (see Section 5.2) or holmium electron capture decays (see Sections 5.3 and 5.4). The main advantage of using microcalorimeters lies in the source = detector principle that allows to measure the complete decay energy (excluding the energy carried away by the emitted neutrino) as opposed to only measuring the kinetic energy of the decay electrons. On the other hand the comparably slow signals produced by calorimetric detectors bring the challenge of using large arrays of small detectors to avoid pileup in the individual calorimeter crystals.
Another approach has been proposed by Jerkins and coworkers  who suggest to perform a full kinematic reconstruction of the decays of trapped tritium atoms. This idea, however, seems to be hampered by conceptual difficulties .
5.1. Project 8
Project 8 is a new effort to measure the neutrino mass lead by groups from MIT, University of California, Santa Barbara and University of Washington. The idea is to use technology from KATRIN's gaseous tritium source combined with a sensitive array of microwave antennae to extract energy spectra of tritium decay electrons from the coherent cyclotron radiation emitted by individual electrons in a magnetic field (see Figure 27). In such a setup, electrons will follow a circular or spiral path with a cyclotron frequency of where is the kinetic energy of the electron, is the magnetic field of the source, is the electron charge, is the electron mass, is the unshifted cyclotron frequency, and is the Lorentz factor. The total power emitted by each electron depends on its relative velocity and the pitch angle between the direction of movement and the magnetic field vector With a suitable set of source parameters, the power of the cyclotron radiation emitted will be large enough to be detected, but not large enough to rapidly change the electron’s momentum and will therefore allow to reconstruct its kinetic energy. The energy resolution of the method depends on the relative uncertainty with which the cyclotron frequency can be determined, which in turn depends on the observation time for an individual decay electron.
Monreal and Formaggio  discuss a reference design with a 1 T magnetic field strength. In this design GHz for decay electrons near the endpoint energy eV. For a resolution of eV we therefore require a relative uncertainty and, via the Nyquist theorem, a minimum observation time of μs. This observation time determines the necessary mean free path of decay electrons in the source before they undergo inelastic scattering from molecules and therefore places constraints on the density of the source. In the reference design the required free path length is given by m corresponding to a maximum density of with a total cross-section cm2 for inelastic scattering of electrons near the endpoint energy from molecular tritium.
Given the long pathlength required for the electrons to achieve the necessary observation time, it becomes clear that the decay electrons need to be trapped within an instrumented volume to realize the actual experiment. This can be achieved using the magnetic mirror effect, where the -field of the source volume is increased by some amount at both ends of the apparatus. The ratio between the field strength in the source volume to that at the edges then determines the minimum pitch angle of the stored electrons. Only storing electrons with pitch angles above has the added benefit that they are also emitting higher signal powers via the dependence of equation 53 on . For stored electrons the emitted signal powers are in the W region while the authors of  estimate the thermal noise contribution of suitable amplifiers for the signal detection to be of the order W. Incoherent noise from non-endpoint and/or low-pitch beta electrons is expected to contribute about 10−24 W/Bq in the signal region which would allow a tritium activity of the source of the order 108–109 Bq. Of course these estimates also depend on how efficient an actual antenna configuration could collect the signal power emitted by the electrons.
The correct reconstruction of electron energies from the detected cyclotron radiation is complicated by the Doppler shift the signals experience due to the motion of the electrons parallel to the magnetic field lines. Two different antenna configurations are proposed to remedy this problem. In the first case, two antennae are placed on axis at both ends of the source volume (see endcap antennae on Figure 27). While one of the antennae will detect a blue-shifted signal due to the parallel velocity component of the electron towards it, the other will detect the red-shifted version of the same signal. Combining these signals, the unshifted electron energy can then be reconstructed.
The second configuration discussed consists of arrays of evenly spaced microwave antennae oriented transverse to the magnetic field. Passing electrons will induce signals that sweep from blue shift to red shift in the individual antennae. When these signals are summed up with an appropriate choice of delay lines, the unshifted cyclotron frequency should sum coherently, while the contributions from Doppler shifts sum incoherently and should thus be suppressed. The resulting output will then appear at the unshifted frequency with an amplitude modulated due to the periodically varying response function of the antennae along the array. This modulation gives rise to two sideband peaks in the signal that can help to discriminate a real high energy electron signal from background signals caused for example by sidebands from lower energy signals.
A simulation performed by Monreal and Formaggio  for tritium decays observed within a μs interval is shown in Figure 28. As discussed above, signals of high energy electrons show up as triplets of lines in the lower part of the frequency spectrum. Because this region is also populated by sideband signals from lower energy electrons, the coincident detection of at least two of the three lines is required to confidently identify an endpoint electron. Extending this requirement to all three lines for the positive identification of a high energy electron, a total source strength of ca. 109 Bq can be allowed to keep background from accidental coincidences on an acceptable level.
The proposed technique presents very different systematic errors than those present in MAC-E filter experiments like KATRIN. Advantages are the ability to correct for fluctuations of the source density, as the complete spectrum is monitored by the experiment. The absence of an electrostatic spectrometer should enable a more scalable experiment as the energy resolution is not limited by the size of a MAC-E filter, but by the observation time for individual electrons and ultimately by the irreducible energy spread due to the excitation of rotational and vibrational final states of the daughter molecule emerging from the tritium decay. On the other hand, care has to be taken to eliminate the effects of Doppler shifts that alter the frequency picked up by the microwave antennae. Inhomogeneities in the -field are a source of line broadening, while drifts of the magnetic field will shift the overall spectrum. Electron- scattering and signal pileup will limit the allowed source density and therefore influence the achievable signal rate in the important endpoint region of the spectrum.
To verify the ability to detect the cyclotron emission of individual electrons and to investigate the achievable resolution, a test experiment  is currently being set up at the University of Washington. In the setup, a small magnetic bottle within a 1 T superconduction magnet is used to capture 17.8 keV conversion electrons from a source. Different antenna designs will be tested in an attempt to observe the cyclotron emission from these electrons.
The “Microcalorimeter Arrays for a Rhenium Experiment” (MARE) collaboration is working to further the development of sensitive microcalorimeters to investigate 187Re -decay. The current activities in the MARE collaboration are organized in two phases : in MARE-1 several groups are working on alternative microcalorimeter concepts which will be tested by setting up neutrino mass experiments with sensitivities in the order of a few eV. Besides the selection of the most sensitive detector technology, this phase will also be used to investigate the use of the EC decay of 163Ho as an alternative to the study of rhenium -decay to determine the neutrino mass  (see also Section 5.3). The consideration of 163Ho was triggered by persisting difficulties with superconducting metallic rhenium absorbers coupled to the sensors . Thermalization of the energy deposited in -decays seems to be hindered by the excitation of long lived states in the rhenium absorber. The nature of these states can presently only be speculated upon . After selecting the most successful technique, a full scale experiment with sub-eV sensitivity to the neutrino mass will then be set up in MARE phase 2.
The technical developments generally aim at two main goals: first to improve the energy resolution of the detectors and secondly to shorten the response time of the signals in order to reduce pileup problems. The various groups in the MARE collaboration work on different techniques to achieve these goals.
Groups from University Milano-Bicocca, NASA/GSFC, and Wisconsin are working together to develop arrays of silicon implanted thermistors coupled to AgReO4 absorbers  (see Figure 29). The experiment that is currently being set up can accommodate up to eight 36 pixel arrays with 0.5 mg AgReO4 absorbers of 0.27 Bq activity each. The absorber crystals used are cut in a regular shape of μm3 from large single crystals and are glued to the thermistors with an intermediate layer of thin pure silicon spacers. For readout of the thermistors a cold buffer stage using JFETs at 120 K is installed as close as possible to the sensors, followed by a main amplifier stage at room temperature. With this system, it is possible to achieve an energy and time resolutions of 25 eV and 250 μs, respectively . After a first test run with 10 AgReO4 crystals on one array and two Sn absorbers on a second array to investigate the environmental background near the rhenium endpoint, a 187Re -decay measurement with 72 channels should start and provide a sensitivity of 4.5 eV at 90% C.L. within three years running time . With all 8 arrays instrumented a sensitivity of 3.3 eV at 90% C.L. would be reached with rhenium decays.
Besides their work on thermistors with AgReO4 absorbers, the group from Milano-Bicocca is also investigating the use of Microwave Kinetic Inductance Detectors (MKIDs) for measurements of 163Ho EC decays. These devices are superconducting resonators in the 1–10 GHz region that exploit the temperature dependence of superconducting films. An advantage with these sensors is the easy readout of a few thousand detectors using frequency multiplexing techniques .
The University of Genoa is working together with groups from University of Miami and University Lisbon/ITN on the development of Transition Edge Sensors (TES) coupled to metallic rhenium absorber crystals . TES sensors exploit the sharp rise of resistance with temperature of a superconductor operated at the phase transition from normal to superconducting behaviour. This allows for a very sensitive detector compared to conventional semiconductor thermistors, but has the drawback of being less stable in operation and having a lower saturation energy. The sensors are based on a Ir-Au bi-layer configured in a double S-shape. The normal resistance of the circuit is 2 Ω and by adjusting the layer thickness, the transition temperature is brought to ca. 80 mK. The superconducting layer is deposited on a silicon substrate onto which, on the other side, a superconducting rhenium absorber of 200–300 mg is glued (see Figure 30). Two Al-Si wires are used for readout of the TES and for connection to a heat bath cooled by a dilution refrigerator. With this design an energy resolution of 11 eV FWHM and a rise time of 160 μs have been achieved . Their planned phase 1 experiment will accommodate 300 TES detectors with about 0.25 Bq activity per element, which would allow a sensitivity on the neutrino mass of 2 eV with three years worth of data taking. As an alternative to rhenium the group at Genoa is also working on the development of TES sensors with 163Ho loaded absorbers.
At the University of Heidelberg, work is ongoing on the development of so-called Metallic Magnetic Calorimeters (MMC) . In contrast to TES sensors or thermistors that exploit the change of resistivity of the sensors with temperature, these detectors are measuring the change in magnetization of a paramagnetic material. Promising results have been obtained with gold absorbers, where resolutions down to 2 eV FWHM and signal rise times of about 90 ns have been achieved for soft X-rays around 6 keV energy . The resolution with rhenium absorbers, however, was found to be 44 eV FWHM at signal rise times below 10 μs. The use of MMC detectors with holmium implanted absorbers produced more encouraging results and led to the formation of the ECHO collaboration to further investigate this technique (see Section 5.4).
All these activities should finally lead to a selection of the most suitable technique for a rhenium or holmium based large-scale neutrino mass experiment with sub-eV sensitivity. According to , a staged approach, where a total of five detector arrays is deployed one per year, would, after 10 years running time, enable a statistical sensitivity better than 0.25 eV. One has to keep in mind however that effects like the beta-environmental fine structure (see Section 3.5) are difficult to estimate and have to be investigated in parallel to the technical developments within the next years.
5.3. Electron Capture on 163Ho
A promising alternative to -decay measurements is the study of electron capture (EC) decays of 163Ho to measure the neutrino mass. The decay process considered is The deexcitation spectrum of the intermediate state is given by a series of lines at energies which correspond to the dissipated binding energy of the electron hole in the fina