A pulsed, mono-energetic and angular-selective UV photo-electron source for the commissioning of the KATRIN experiment

TL;DR

This paper presents a pulsed, mono-energetic, and angular-selective UV photo-electron source designed for calibrating the KATRIN experiment's main spectrometer, with experimental results validated against detailed simulations.

Contribution

It introduces a novel electron source with specific energy and angular properties for KATRIN calibration, and demonstrates its performance through measurements and simulations.

Findings

Successful development of the electron source.

Agreement between measurements and simulations.

Enhanced calibration capabilities for KATRIN.

Abstract

The KATRIN experiment aims to determine the neutrino mass scale with a sensitivity of 200 meV/c^2 (90% C.L.) by a precision measurement of the shape of the tritium $\beta$-spectrum in the endpoint region. The energy analysis of the decay electrons is achieved by a MAC-E filter spectrometer. To determine the transmission properties of the KATRIN main spectrometer, a mono-energetic and angular-selective electron source has been developed. In preparation for the second commissioning phase of the main spectrometer, a measurement phase was carried out at the KATRIN monitor spectrometer where the device was operated in a MAC-E filter setup for testing. The results of these measurements are compared with simulations using the particle-tracking software "Kassiopeia", which was developed in the KATRIN collaboration over recent years.

Tables

Table 1. Table 1: Measured transmission functions at different wavelengths λ 𝜆 \lambda and fixed spectrometer voltage U spec = − 18.6   kV subscript 𝑈 spec times -18.6 kV U_{\mathrm{spec}}=$-18.6\text{\,}\mathrm{k}\mathrm{V}$ (fig.  7 ). The table shows the upper limit of the energy distribution, E max subscript 𝐸 max E_{\mathrm{max}} , and the energy spread, σ E subscript 𝜎 𝐸 \sigma_{E} ; both values are derived from the fit result. The measurement at 266   nm times 266 nm 266\text{\,}\mathrm{n}\mathrm{m} has been performed twice at − 18.6   kV times -18.6 kV -18.6\text{\,}\mathrm{k}\mathrm{V} .

$\lambda$	$E_{\max }$	${\sigma }_E$	${\chi }^2/ndf$
$266.0\mathrm{nm}$	$0.82\left(2\right)\mathrm{eV}$	$0.31\left(5\right)\mathrm{eV}$	$1.39$
$266.0\mathrm{nm}$	$0.82\left(2\right)\mathrm{eV}$	$0.28\left(4\right)\mathrm{eV}$	$1.40$
$272.4\mathrm{nm}$	$0.74\left(1\right)\mathrm{eV}$	$0.22\left(2\right)\mathrm{eV}$	$1.18$
$282.4\mathrm{nm}$	$0.61\left(2\right)\mathrm{eV}$	$0.19\left(3\right)\mathrm{eV}$	$1.23$
$292.4\mathrm{nm}$	$0.47\left(2\right)\mathrm{eV}$	$0.14\left(3\right)\mathrm{eV}$	$3.38$
$302.4\mathrm{nm}$	$0.33\left(2\right)\mathrm{eV}$	$0.09\left(7\right)\mathrm{eV}$	$3.46$

Table 2. Table 2: Measured magnetic reflection curves at different azimuthal directions α az subscript 𝛼 az \alpha_{\mathrm{az}} of the plate angle (figure  8 ). The table shows the reflection angle α ^ = α max ^ 𝛼 subscript 𝛼 max \hat{\alpha}=\alpha_{\mathrm{max}} and the width σ α subscript 𝜎 𝛼 \sigma_{\alpha} (in terms of plate angle) of the reflection curve that was determined by the fit. The angular spread σ θ subscript 𝜎 𝜃 \sigma_{\theta} (in terms of pitch angle) has been computed from the adiabatic transformation eq. ( 15 ) with the known reflection angle.

${\alpha }_{\mathrm{az}}$	${\alpha }_{\max }$	${\sigma }_{\alpha }$	${\sigma }_{\theta }$	${\chi }^2/ndf$
$0°$	$10.06\left(2\right)°$	$0.39\left(3\right)°$	$16.0\left(10\right)°$	$0.71$
$180°$	$11.13\left(3\right)°$	$0.39\left(3\right)°$	$15.2\left(12\right)°$	$0.88$
$90°$	$9.56\left(2\right)°$	$0.40\left(2\right)°$	$16.7\left(9\right)°$	$1.89$
$270°$	$9.73\left(2\right)°$	$0.40\left(3\right)°$	$16.6\left(10\right)°$	$0.82$
weighted average:		$0.396\left(1\right)°$	$16.20\left(7\right)°$

Table 3. Table 3: Measured transmission functions at different plate angles α p subscript 𝛼 p \alpha_{\mathrm{p}} . The table shows the mean angle θ ^ ^ 𝜃 \hat{\theta} and the angular spread σ θ subscript 𝜎 𝜃 \sigma_{\theta} in the spectrometer entrance magnet; the values have been determined by the fit. An expected pitch angle θ ^ ana subscript ^ 𝜃 ana \hat{\theta}_{\mathrm{ana}} is derived analytically from adiabatic transformation eq. ( 5 ). At α p ≥ α max = 10 ​ ° subscript 𝛼 p subscript 𝛼 max 10° \alpha_{\mathrm{p}}\geq\alpha_{\mathrm{max}}=$$ magnetic reflection is observed, which leads to a significant deformation of the transmission function.

${\alpha }_p$	$\hat{\theta }$	${\sigma }_{\theta }$	${\chi }^2/ndf$	${\hat{\theta }}_{\mathrm{ana}}$
$0°$	$1.7\left(13\right)°$	$2.0\left(13\right)°$	$1.09$	$2.0°$
$2°$	$5.7\left(34\right)°$	$9.3\left(26\right)°$	$1.07$	$13.5°$
$4°$	$23.2\left(3\right)°$	$5.8\left(8\right)°$	$1.12$	$25.4°$
$6°$	$38.2\left(2\right)°$	$4.3\left(5\right)°$	$1.31$	$38.6°$
$8°$	$55.2\left(3\right)°$	$5.6\left(4\right)°$	$1.50$	$54.7°$
$10°$	$89.3\left(8\right)°$	$0.8\left(7\right)°$	$10.9$	$85.7°$

Table 4. Table 4: Work functions determined from measured transmission functions at different wavelengths λ 𝜆 \lambda (fig.  7 ). The work functions Φ † superscript Φ † \Phi^{\dagger} are derived from eq. ( 19 ), with E e , max subscript 𝐸 𝑒 max E_{e,\mathrm{max}} the upper limit of the energy distribution (tab.  1 ) and h ​ c / λ ℎ 𝑐 𝜆 hc/\lambda the known photon energy. The results are compared with the work function Φ = 3.78 ​ ( 4 )   eV Φ times 3.78 4 eV \Phi=$3.78(4)\text{\,}\mathrm{e}\mathrm{V}$ that was determined in a Fowler-type measurement (fig.  10 ).

$\lambda$	$hc/\lambda$	${\Phi }^{†}$	${\Phi }^{†}-\Phi$
$266.0\mathrm{nm}$	$4.66\left(2\right)\mathrm{eV}$	$3.84\left(4\right)\mathrm{eV}$	$0.05\left(7\right)\mathrm{eV}$
$266.0\mathrm{nm}$	$4.66\left(2\right)\mathrm{eV}$	$3.84\left(4\right)\mathrm{eV}$	$0.05\left(7\right)\mathrm{eV}$
$272.4\mathrm{nm}$	$4.55\left(4\right)\mathrm{eV}$	$3.81\left(5\right)\mathrm{eV}$	$0.02\left(8\right)\mathrm{eV}$
$282.4\mathrm{nm}$	$4.39\left(4\right)\mathrm{eV}$	$3.79\left(5\right)\mathrm{eV}$	$0.00\left(8\right)\mathrm{eV}$
$292.4\mathrm{nm}$	$4.24\left(3\right)\mathrm{eV}$	$3.77\left(5\right)\mathrm{eV}$	$-0.02\left(8\right)\mathrm{eV}$
$302.4\mathrm{nm}$	$4.10\left(3\right)\mathrm{eV}$	$3.77\left(5\right)\mathrm{eV}$	$-0.02\left(8\right)\mathrm{eV}$
weighted average:		$3.810\left(1\right)\mathrm{eV}$	$0.03\left(3\right)\mathrm{eV}$

Table 5. Table 5: Simulated pitch angles in the spectrometer entrance magnet and derived effective starting angles. The table shows the median pitch angle, θ ^ ^ 𝜃 \hat{\theta} , and the angular spread, σ θ subscript 𝜎 𝜃 \sigma_{\theta} , in the entrance magnet for different plate angles α p subscript 𝛼 p \alpha_{\mathrm{p}} . The simulation results are compared with the pitch angle determined from corresponding measurements, θ ^ meas subscript ^ 𝜃 meas \hat{\theta}_{\mathrm{meas}} . An effective initial pitch angle at the photocathode, θ start ∗ superscript subscript 𝜃 start \theta_{\mathrm{start}}^{*} , can be computed from θ ^ ^ 𝜃 \hat{\theta} using the adiabatic transformation eq. ( 21 ) and the known magnetic fields at the setup ( B start = 27   mT subscript 𝐵 start times 27 mT B_{\mathrm{start}}=$27\text{\,}\mathrm{m}\mathrm{T}$ , B max = 6   T subscript 𝐵 max times 6 T B_{\mathrm{max}}=$6\text{\,}\mathrm{T}$ ).

${\alpha }_p$	$\hat{\theta }$	${\sigma }_{\theta }$	${\hat{\theta }}_{\mathrm{meas}}$	${\theta }_{\mathrm{start}}^{\ast }$
$0°$	$1.5°$	$1.2°$	$1.7\left(13\right)°$	$0.1°$
$2°$	$11.5°$	$1.3°$	$5.7\left(34\right)°$	$0.9°$
$4°$	$23.6°$	$1.3°$	$23.2\left(3\right)°$	$1.6°$
$6°$	$36.8°$	$1.5°$	$38.2\left(2\right)°$	$2.4°$
$8°$	$52.5°$	$1.9°$	$55.2\left(3\right)°$	$3.2°$
$10°$	$78.9°$	$4.2°$	$89.3\left(8\right)°$	$3.9°$
$12°$	— magnetically reflected —			$4.6°$

Equations

$$E_{\perp }^{}=E\cdot {\sin }^2\theta ,E_{\parallel }^{}=E\cdot {\cos }^2\theta .$$

$$\mu =\frac{E_{\perp }^{}}{|{\vec{B}}^{}|}=const.\text{for}\left|\frac{1}{B}\frac{\mathrm{d}B}{\mathrm{d}z}\right|\ll \frac{{\omega }_{\mathrm{c}}}{v_{\parallel }^{}},$$

$$q{U}_{\mathrm{ana}}<E_{\parallel }^{}=E_0\cdot \left(1-{\sin }^2{\theta }_0\cdot \frac{B_{\min }}{B_{\max }}\right).$$

$$\Delta E={E}^{\mathrm{max}}_{\perp}=E_{0}\cdot\frac{B_{\mathrm{min}}}{B_{\mathrm{max}}}=$0.93\text{\,}\mathrm{e}\mathrm{V}$\,.$$

$${\theta }_0\le {\theta }_{\max }=\arcsin \sqrt{\frac{B_{\mathrm{source}}}{B_{\max }}},$$

$$q\Delta U=q(U_{\mathrm{start}}-U_{\mathrm{spec}}),$$

$$0<E_e\le h\nu -\Phi =hc/\lambda -\Phi .$$

$$\vec{F}=q(\vec{E}+\vec{v}\times \vec{B}).$$

$${\vec{v}}_{\mathrm{drift}}=\frac{\vec{E}\times \vec{B}}{B^2},$$

$$T(E,U_{\mathrm{ana}})={\dot{N}}_0\cdot {\int }_E^{\infty }\eta (ϵ){\int }_0^{{\theta }_{\max }(ϵ,U_{\mathrm{ana}})}\zeta (\theta )\mathrm{d}\theta \mathrm{d}ϵ+{\dot{N}}_{\mathrm{b}},$$

$${\theta }_{\max }(E,U_{\mathrm{ana}})=\{\begin{array}{l}\arcsin \left(\sqrt{\frac{E}{U_{\mathrm{ana}}}\cdot \frac{2}{\gamma +1}\cdot \frac{B_{\mathrm{start}}}{B_{\min }}}\right),\\ 0\text{for}\sqrt{\dots }<0.\end{array}$$

$$\eta (E)=\frac{1}{\sqrt{2\pi }}\cdot \{\begin{array}{ll}\frac{1}{{\alpha }_E}\cdot \exp \left(-\frac{1}{2}\frac{(E-\hat{E}{)}^2}{{\alpha }_E^2}\right)& (\kappa =0),\\ \frac{1}{{\alpha }_E-\kappa (E-\hat{E})}\cdot & \\ \cdot \exp \left(-\frac{1}{2{\kappa }^2}\ln {\left[1-\kappa \frac{E-\hat{E}}{{\alpha }_E}\right]}^2\right)& (\kappa \ne 0),\end{array}$$

$${\sigma }_E=\frac{{\alpha }_E}{\kappa }\cdot \sqrt{e^{{\kappa }^2}(e^{{\kappa }^2}-1)}.$$

$$\zeta (\theta )=\frac{1}{\sqrt{2\pi }\sigma }\cdot \left[\exp \left(-\frac{(\theta -\hat{\theta }{)}^2}{2{\sigma }_{\theta }^2}\right)+\exp \left(-\frac{(\theta +\hat{\theta }{)}^2}{2{\sigma }_{\theta }^2}\right)\right],$$

$$\theta \approx \arcsin \left({\alpha }_{\mathrm{p}}\cdot k\cdot \sqrt{\frac{B_{\max }}{B_{\mathrm{start}}}}\right)\approx \arcsin \left(\frac{{\alpha }_{\mathrm{p}}}{{\alpha }_{\max }}\right),$$

$$\Delta E_{\mathrm{ref}}=$18.6\text{\,}\mathrm{k}\mathrm{V}$\cdot\frac{$0.38\text{\,}\mathrm{m}\mathrm{T}$}{$6\text{\,}\mathrm{T}$}=$1.18\text{\,}\mathrm{e}\mathrm{V}$$

$

$$I(\mu )$$

$$\xi (\mu )$$

$$E_{e,\max }=h\nu -\Phi =hc/\lambda -\Phi ,$$

$$\theta =\arcsin \left(\sqrt{{\sin }^2{\theta }_{\mathrm{start}}^{\ast }\cdot \frac{B_{\max }}{B_{\mathrm{start}}}}\right).$$

$$\theta =\arcsin \left({\alpha }_{\mathrm{p}}\cdot k\cdot \sqrt{\frac{B_{\max }}{B_{\mathrm{start}}}}\right).$$

Full text

††journal: Eur. Phys. J. C

1 Introduction

The KArlsruhe Tritium Neutrino experiment KATRIN DesignReport2005 aims to measure an ‘effective mass’ of the electron anti-neutrino, given by an incoherent sum over the mass eigenstates NuMass2013 . It performs kinematic measurements of tritium $\beta$-decay to achieve a neutrino mass sensitivity down to $200\text{\,}\mathrm{m}\mathrm{e}\mathrm{V}\mathrm{/}\mathrm{c}^{2}$ at 90% C. L., improving the results of the predecessor experiments in Mainz NuMassMainz2005 and Troitsk NuMassTroitsk2011 by one order of magnitude. As the evolution of the neutrino mass results of these experiments showed, the study of systematic effects is of major importance: Underestimated or unknown “energy loss” processes caused too positive or even negative values for the square of the neutrino mass NuMass2008 . A detailed understanding of systematic uncertainties at the KATRIN experiment is crucial to achieve its target sensitivity.

The outline of the KATRIN experiment is depicted in figure 1 DesignReport2005 ; NuMass2013 . Molecular tritium is fed into the $10\text{\,}\mathrm{m}$ long beam tube of the windowless gaseous tritium source (WGTS WGTS2015 ). Superconducting magnets along the beam line create an adiabatic guiding field in a $191\text{\,}\mathrm{T}\,\mathrm{c}\mathrm{m}^{2}$ magnetic flux tube, and $\beta$-decay electrons emitted in forward direction propagate towards the spectrometer section. The electrons then enter the transport and pumping section that reduces the tritium flow by a factor of ${10}^{14}$ in total DPS2006 , using a combination of a differential pumping section (DPS DPS2012 ) with turbo-molecular pumps and a cryogenic pumping section (CPS CPS2010 ) where tritium is adsorbed by an argon frost layer. The kinetic energy of the decay electrons is analyzed in a tandem of MAC-E filter111Magnetic Adiabatic Collimation with Electrostatic filter. spectrometers MACE1992 ; MACE1985 ; MACE1981 . The main spectrometer achieves an energy resolution of $0.93\text{\,}\mathrm{e}\mathrm{V}$ at the tritium endpoint of $E_{0}(\text{T}_{2})=$18,571.8(12)\text{,}\mathrm{e}\mathrm{V}$$ NuMass2008 ; TritiumEndpoint2015 by a combination of an electrostatic retarding potential and a magnetic guiding field. Electrons with sufficient kinetic energy pass the retarding potential and are counted at the focal-plane detector (FPD FPD2015 ). An integral energy spectrum is measured by varying the filter energy close to the tritium endpoint. The effective neutrino mass is determined by fitting the convolution of the theoretical $\beta$-spectrum with the response function of the spectrometer to the data, taking into account important parameters such as the final states distribution and the energy loss spectrum and other systematic corrections DesignReport2005 ; WGTS2012 . The spectrometer high-voltage is monitored by a pair of precision high-voltage dividers HVdivider2013 ; HVdivider2009 that support voltages up to $35\text{\,}\mathrm{k}\mathrm{V}$ and $65\text{\,}\mathrm{k}\mathrm{V}$, respectively. An absolute voltage calibration is achieved by measuring the divider’s output voltage with $\mathrm{p}\mathrm{p}\mathrm{m}$ precision using a digital voltmeter. Additionally, the stability of the retarding potential is monitored continuously at the monitor spectrometer MonSpec2014 . This MAC-E filter is connected to the main spectrometer high voltage system and measures natural conversion lines of ${}^{\text{83m}}\text{Kr}$, where changes in the retarding potential are observed as shifts in the measured line position.

A precise knowledge of the transmission properties of the KATRIN main spectrometer is crucial to limit systematic uncertainties and reach the desired neutrino mass sensitivity. The transmission properties are affected by inhomogeneities of the electromagnetic fields in the main spectrometer. In addition to simulations, dedicated measurements are necessary to determine the spectrometer transmission function over the complete magnetic flux tube. Such measurements require a mono-energetic and angular-selective electron source, which we present in this work. A pulsed electron beam allows us to access additional information from the electron time-of-flight (ToF) ToF2013 .

This article is structured as follows: Section 2 discusses the revised technical design of the photoelectron source that was developed at WWU Münster over the recent years Egun2009 ; Egun2011 ; Egun2014 . The design underwent many improvements for the second commissioning phase of the KATRIN main spectrometer. In section 3 we show results from test measurements at the KATRIN monitor spectrometer. We determine important source characteristics such as the energy and angular spread of the produced electrons and the effective work function of the photocathode. Section 4 discusses simulation results that were produced by Kassiopeia, a particle-tracking software that has been developed as a joint effort in the KATRIN collaboration over recent years Kassiopeia2016 . These simulations allow us to gain a detailed understanding of the electron acceleration and transport processes inside the electron source.

2 Setup and Design

2.1 Principle of the MAC-E filter

The principal design of the MAC-E filter is based on the combination of an electric retarding potential with a spatially inhomogeneous magnetic field MACE1992 . In the following we describe this principle on the basis of the technical implementation at the KATRIN experiment. Two solenoids located at the entrance and exit regions produce a strong magnetic field $B_{\mathrm{max}}$, which drops to a minimal value $B_{\mathrm{min}}$ at the central plane of the spectrometer. The value $B_{\mathrm{min}}$ can be adjusted by a system of air coils, which are placed around the spectrometer. The beam tube and the electrodes at the spectrometer entrance and exit are on ground potential, while the spectrometer vessel and the central electrodes are operated at high voltage. The absolute value of the retarding potential increases towards the central spectrometer plane and reaches a maximum of $U_{\mathrm{ana}}\approx$-18.6\text{,}\mathrm{k}\mathrm{V}$$ at the position of the magnetic field minimum. This point lies on the so-called analyzing plane. The electromagnetic conditions in the analyzing plane define the transmission function for electrons that propagate through the spectrometer. Inside the MAC-E filter, electrons follow a cyclotron motion around the magnetic field lines. The kinetic energy $E$ can be split into a longitudinal component ${E}^{\mathrm{}}_{\parallel}$ into the direction of the field line and a transversal component ${E}^{\mathrm{}}_{\perp}$, wich corresponds to the gyration around the field line. Both components of the electron’s kinetic energy can be described by the polar angle of the electron momentum relative to the magnetic field line, the pitch angle $\theta=\angle{(\vec{p},\vec{B})}$:

[TABLE]

The adiabatic motion of the electrons is one of the key features of the MAC-E filter. When the relative change of the magnetic field over one cyclotron turn is small, the magnetic moment $\mu$ is conserved (here written non-relativistically):

[TABLE]

where $\omega_{\mathrm{c}}$ denotes the cyclotron frequency of the electron. The reduction of the magnetic field towards the analyzing plane leads to a decrease in transversal energy ${E}^{\mathrm{}}_{\perp}$. The longitudinal component ${E}^{\mathrm{}}_{\parallel}$ increases accordingly in this process because of energy conservation. This behavior results in a momentum collimation of the electron beam, and electrons that enter the spectrometer at a strong magnetic field $B_{\mathrm{max}}$ reach a minimal transversal energy ${E}^{\mathrm{}}_{\perp}$ at $B_{\mathrm{min}}$ in the analyzing plane. Because the retarding potential $U_{\mathrm{ana}}$ only analyzes the longitudinal energy component, this principle of adiabatic collimation allows measuring the energy of electrons from an isotropic source with high precision. The transmission condition for an electron with charge $q$ that enters the spectrometer with energy $E_{0}$ and pitch angle $\theta_{0}$ is

[TABLE]

At nominal conditions, the KATRIN main spectrometer achieves a minimal magnetic field of $B_{\mathrm{min}}=$0.3\text{,}\mathrm{m}\mathrm{T} in the analyzing plane and a maximal magnetic field of $B_{\mathrm{max}}=$6\text{\,}\mathrm{T} at the pinch magnet, which is positioned at the exit of the main spectrometer.

The energy resolution (more precisely: filter width) $\Delta E$ of the MAC-E filter for an isotropic source is derived by postulating an electron that enters the spectrometer with maximum pitch angle, \theta_{0}=$$$ or {E}^{\mathrm{start}}{\perp}=E{0}$. The energy resolution corresponds to the remaining transversal kinetic energy in the analyzing plane after adiabatic collimation in the MAC-E filter:

[TABLE]

At KATRIN, electrons start in the source at $B_{\mathrm{source}}<B_{\mathrm{max}}=$6\text{,}\mathrm{T}$$. This leads to magnetic reflection of electrons with large pitch angles at the pinch magnet. The magnetic mirror effect occurs independently of the spectrometer transmission condition and reduces the acceptance angle of the MAC-E filter,

[TABLE]

and electrons created with larger pitch angles do not contribute to the measurement. At KATRIN, the acceptance angle is \theta_{\mathrm{max}}=$$$ for electrons starting in the WGTS with B_{\mathrm{source}}=$3.6\text{,}\mathrm{T}$$. This excludes electrons with excessive path lengths as a result of their cyclotron motion, and thereby reduces systematic uncertainties caused by energy losses.

The KATRIN beam line transports a maximum magnetic flux of $191\text{\,}\mathrm{T}\,\mathrm{c}\mathrm{m}^{2}$ from the source to the detector. Electrons that are created at the source follow different magnetic field lines, depending on their initial radial and azimuthal position. The transmission function for electrons is affected by inhomogeneities in the analyzing plane of the electric potential ($\Delta U_{\mathrm{ana}}<$1.2\text{,}\mathrm{V}) and the magnetic field ($\Delta B_{\mathrm{min}}<$50\text{\,}\mathrm{\SIUnitSymbolMicro}\mathrm{T}). Because these variations are too large to be neglected, the detector features a pixelated wafer that can adequately resolve the position in the analyzing plane. This allows us to consider the electromagnetic inhomogeneities by determining transmission functions for individual detector pixels. The exact value of the retarding potential $U_{\mathrm{ana}}$ and the magnetic field $B_{\mathrm{min}}$ can be accessed through measurements with an electron source that generates electrons at defined kinetic energy and pitch angle. A source that fulfills these requirements has been developed at WWU Münster for the commissioning of the KATRIN main spectrometer.

2.2 Principle of the electron source

It was demonstrated in Egun2011 that angular selectivity can be achieved by a combination of non-parallel electric and magnetic fields. An earlier design that used a gold-plated quartz tip, which was illuminated by UV light from optical fibers on the inside of the tip was able to produce electrons with non-zero pitch angles. This setup achieved an insufficiently large angular spread of the electrons. The source was therefore not usable as a calibration source for a MAC-E filter. The design was further refined in Egun2014 and the setup now resembles a plate capacitor that introduces a homogeneous electric acceleration field. The setup can be tilted against the magnetic field lines to imprint a well-defined pitch angle on the generated electrons. This design uses a planar photocathode, which is back-illuminated by UV light from a single optical fiber.

This setup is shown in figure 2. The emission spot $p_{e}$ is located on the back plate (red), which is put on a negative potential $U_{\mathrm{start}}$ and thus defines the kinetic energy of the generated electrons, $E_{\mathrm{kin}}=qU_{\mathrm{start}}$. The surplus energy of the electrons in the analyzing plane,

[TABLE]

then amounts to the remaining kinetic energy that is available to overcome the retarding potential $U_{\mathrm{ana}}$ of the spectrometer. The inhomogeneity $\Delta U_{\mathrm{ana}}$ of the retarding potential caused by the finite dimensions of the spectrometer is called potential depression. It results in an effective retarding potential $U_{\mathrm{ana}}=U_{\mathrm{spec}}+\Delta U_{\mathrm{ana}}$ that is more positive than the spectrometer voltage $U_{\mathrm{spec}}$. The value $U_{\mathrm{ana}}$ is affected by further inhomogeneities of the electromagnetic conditions in the spectrometer (e. g. work function fluctuations), which can be resolved by transmission function measurements with our electron source.

The front plate (blue) with an aperture for electrons is mounted parallel to the back plate and placed in front of the emission spot. A potential difference $U_{\mathrm{acc}}=U_{\mathrm{front}}-U_{\mathrm{start}}\leq$5\text{,}\mathrm{k}\mathrm{V}$$ is applied between the plates to create an electric field perpendicular to the photocathode surface. The plates are mounted inside a grounded cage (yellow) to shield the electric field at the photocathode against outside influences. The whole setup can be mechanically tilted against the direction of the magnetic field. After passing the front plate, the electrons are accelerated adiabatically towards the ground potential at the spectrometer entrance where they achieve their maximum kinetic energy.

The electrons are emitted from a photocathode that consists of a thin gold (or silver) layer. The photocathode is back-illuminated via an optical fiber by UV light with a variable wavelength $\lambda$, which can be tuned to match the work function $\Phi$ of the photocathode material, $\lambda\lesssim hc/\Phi$. The energy distribution of the emitted electrons is defined by the photon energy $h\nu=hc/\lambda$ and the work function $\Phi$:

[TABLE]

For metallic surfaces used in our electron source, the work function is $\Phi<$5\text{,}\mathrm{e}\mathrm{V}$$. Factors such as surface roughness or impurities caused by adsorbed gas molecules modify the Fermi level at the photocathode, which typically reduces the observed work function Photoemission2012 . It is possible to perform an in situ measurement of the work function using the well-known approach by Fowler Photoemission1931 . In section 3.6 we present results from applying this technique.

In our electron source setup with a planar photocathode, the initial emission angle of the emitted electrons is expected to follow a $\cos\theta$-distribution Photoemission2002 . A well-defined pitch angle $\theta$ is created by collimating the electron beam with the electric field $\vec{E}$ between the two plates. The electrostatic acceleration dominates the electron propagation because of their low kinetic energy after emission, according to the Lorentz equation

[TABLE]

The magnetic guiding field $\vec{B}$ takes over as the electrons gain more kinetic energy, and eventually the electrons enter an adiabatic cyclotron motion around the magnetic field line. The resulting pitch angle of the electrons in the spectrometer entrance magnet is minimal if the plate setup is aligned with the magnetic field, \alpha_{\mathrm{p}}=$$$. By tilting the source against the magnetic field by the *plate angle* \alpha_{\mathrm{p}}>, the non-adiabatic acceleration by the electric field works against the magnetic guiding field. This increases the transversal kinetic energy of the electrons, thereby creating an angular distribution of gaussian shape with a defined mean pitch angle $\theta>. Because the plate setup is located inside a grounded cage, the electric acceleration field at the photocathode is constant for different plate angles.

The pitch angle of the emitted electrons transforms adiabatically during propagation towards the spectrometer entrance magnet, where the electrons enter a higher magnetic field. According to eq. (2) the transformation depends on the ratio of the magnetic fields at the emission spot, $B_{\mathrm{start}}$, and the magnetic field in the entrance magnet, $B_{\mathrm{mag}}\leq B_{\mathrm{max}}$. The pitch angle increases because $B_{\mathrm{start}}\ll B_{\mathrm{mag}}$. The electron source we present here produces defined pitch angles that cover the full range of $\theta=$0\text{,}\mathrm{\SIUnitSymbolDegree}90\text{,}\mathrm{\SIUnitSymbolDegree}$$ in the entrance magnet with plate angles \alpha_{\mathrm{p}}\leq$$$ (section [3.5](#S3.SS5)). The energy spread in the entrance magnet is defined by the initial energy distribution of the emitted electrons, because the acceleration by the electric field does not deform the energy distribution. The kinetic energy is merely shifted by qU_{\mathrm{start}}$, while the spectral shape of the initial distribution is unaffected. A spectroscopic measurement of the electron energy, e. g. with a MAC-E filter, therefore allows us to determine the initial energy distribution of the emitted electrons. The method is also suited to determine the photocathode work function, which is discussed in section 3.3.

2.3 Technical design

The technical design of the electron source is based on the plate-capacitor setup depicted in figure 2. We use two stainless steel disks with radius $r_{\mathrm{p}}=$30\text{,}\mathrm{m}\mathrm{m} for the front and back plate, which are placed at a distance $d=$10\text{\,}\mathrm{m}\mathrm{m}. Both plates were electro-polished before installation. The front plate has a thickness of $d_{\mathrm{fp}}=$2\text{,}\mathrm{m}\mathrm{m} and features an aperture with a radius $r_{\mathrm{afp}}=$3\text{\,}\mathrm{m}\mathrm{m} for the emitted electrons. The back plate has a thickness of $d_{\mathrm{bp}}=$3\text{,}\mathrm{m}\mathrm{m} and allows mounting a photocathode holder at its center. The holder has an aperture to glue-in an optical fiber with diameter $200\text{\,}\mathrm{\SIUnitSymbolMicro}\mathrm{m}$. The holder with the optical fiber is manually polished to create a flat surface, and the photocathode material is deposited on the surface by electron beam physical vapor deposition (EBPVD). For the measurements presented here we used a gold photocathode with a layer thickness of $20\text{\,}\mathrm{n}\mathrm{m}$; we also used silver with a thickness of $40\text{\,}\mathrm{n}\mathrm{m}$ in other measurements. The plates are isolated against each other and the grounded cage by polyether ether ketone (PEEK) insulators. The grounded cage has an inner radius of $r_{\mathrm{c}}=$50\text{\,}\mathrm{m}\mathrm{m} with an aperture $r_{\mathrm{ac}}=$35\text{,}\mathrm{m}\mathrm{m}$$ at the front.

The grounded cage is gimbal-mounted to allow tilting against two axes. The center of rotation is aligned with the emission spot on the back plate. This design ensures that the magnetic field line that the electron is following does not change when tilting the source cage. A precise readout of the plate angle is achieved by rotating piezo-electric motors (Attocube ANR240) that are installed at the pivot joints of the gimbal mount. These motors do not provide sufficient torque to tilt the electron source under vacuum conditions, but allow the relative tilt angle to be measured with a precision of [math]. To actuate the gimbal mount under vacuum conditions, our design uses two air-pressure linear motors (Bibus Tesla 1620) that are mounted outside the vacuum chamber. The linear motion of the motors is transferred onto the chamber by Bowden cables that are attached to each axis of the gimbal mount. By operating the motors, each axis can be tilted separately. The motors are controlled with a LabView software, which also takes care of the transformation between the two-axial and polar/azimuthal coordinate system for the plate angles.

When electrons are reflected by the electric retarding potential at the analyzing plane or by the magnetic field at the spectrometer entrance, they may become stored between the spectrometer and the electron source in a setup similar to a Penning trap. This can lead to a discharge, which has disastrous consequences for the photocathode. To avoid such storing conditions and the subsequent discharge, a dipole electrode is placed in the beamline, between the source cage and the entrance magnet. The electric field $\vec{E}$ induced by the electrode results in a drift of the stored electrons,

[TABLE]

with $\vec{B}$ the magnetic field at the dipole electrode. In our setup, the half-shell dipole electrode spans [math] at a radius of $r_{\mathrm{dip}}=$30\text{,}\mathrm{m}\mathrm{m} and is operated at a voltage $U_{\mathrm{dip}}\leq$4\text{\,}\mathrm{k}\mathrm{V}. Measurements confirmed that this electrode removes trapped electrons efficiently and prevents Penning discharges; this is discussed in section 3.2.

The optical system to provide the UV light for the photocathode allows choosing between two light sources. A frequency-quadrupled $\text{Nd:YVO}_{4}$ laser (InnoLas mosquito-266-0.1-V) provides UV light at a wavelength of $266\text{\,}\mathrm{n}\mathrm{m}$ ($1\text{\,}\mathrm{n}\mathrm{m}$ FWHM) at high intensity ($<$10\text{,}\mathrm{m}\mathrm{W} output power). The intensity can be adjusted by an internal attenuator ($\lambda/2$-plate with polarizing filter) and by a neutral density (ND) filter, which is placed in the laser beam. Behind the ND filter, a fraction of approximately 0.5% of the UV light is coupled out by a beam splitter to measure the UV light intensity with a photodiode. The laser light is focused by an aspheric lens into a $\varnothing$200\text{\,}\mathrm{\SIUnitSymbolMicro}\mathrm{m} optical fiber and guided into the source chamber. The laser is operated in pulsed mode with frequencies of $40\text{\,}\mathrm{k}\mathrm{H}\mathrm{z}100\text{\,}\mathrm{k}\mathrm{H}\mathrm{z}$ at a pump diode current of $6\text{\,}\mathrm{A}8\text{\,}\mathrm{A}$. The current and frequency setting determines the output power, which can be tuned to produce a desired electron rate of several $\mathrm{k}\mathrm{c}\mathrm{p}\mathrm{s}$ ($\mathrm{c}\mathrm{p}\mathrm{s}$: counts per second) at the detector. The pulse width of $<$20\text{,}\mathrm{n}\mathrm{s}$$ allows for time-of-flight measurements with a precisely known starting time of the electrons.

Alternatively, an array of LEDs can be used as light source to provide UV light with $\lambda=$260\text{,}\mathrm{n}\mathrm{m}320\text{,}\mathrm{n}\mathrm{m}$$. Six ball-lens UV LEDs (Roithner UVTOP260–310) with peak wavelengths of $265\text{\,}\mathrm{n}\mathrm{m}275\text{\,}\mathrm{n}\mathrm{m}$ etc. on are mounted on a revolver that is moved by a stepper motor. This allows us to automatically place the desired LED on the optical axis without manual adjustments. To achieve a sharp line width, a UV monochromator with $4\text{\,}\mathrm{n}\mathrm{m}$ FWHM is used. The monochromator is operated by another stepper motor. The LED revolver in combination with the monochromator allows selecting arbitrary wavelengths in the available range. Like in the laser setup, a beam splitter with photodiode is used to monitor the light intensity. The divergent light beam of the LEDs is focused by an optical telescope consisting of two convex lenses, and guided into the electron source through an optical fiber. The current to operate the LEDs is provided by a function generator in pulse mode, using the internal $50\text{\,}\mathrm{\SIUnitSymbolOhm}$ resistor with an output voltage of $8.5\text{\,}\mathrm{V}$. With this setting, the LEDs are driven by a peak current of $200\text{\,}\mathrm{m}\mathrm{A}$, which corresponds to a mean current of $20\text{\,}\mathrm{m}\mathrm{A}$ at 10% duty cycle. Under nominal conditions, a pulse frequency of $100\text{\,}\mathrm{k}\mathrm{H}\mathrm{z}$ is used with a pulse length of $1\text{\,}\mathrm{\SIUnitSymbolMicro}\mathrm{s}$. Time-of-flight measurements are thus also possible with LEDs as a light source. Depending on the LED and the monochromator setting, electron rates in the $\mathrm{k}\mathrm{c}\mathrm{p}\mathrm{s}$ range can be achieved.

The optical system (laser device, the stepper motors of the LED system and the two photodiodes), the actuation of the plate angle and the power supply for the dipole electrode are controlled and monitored by a LabView software that has been developed for use with the electron source. The photodiode read-out allows us to monitor the stability of the UV light source, where intensity changes (e. g. because of warm-up effects) could result in fluctuations of the observed electron rate.

2.4 Analytical transmission function

The observed transmission functions from measurements with the electron source can be modeled by an analytical description of the MAC-E filter groh:phd . The conditions for transmission eq. (3) and magnetic reflection eq. (5) are applied to the theoretical energy distribution $\eta(E)$ and angular distribution $\zeta(\theta)$. The analytical transmission function $T(E)$ is given by the integrated energy distribution, which is modified by the range of pitch angles that are transmitted through the spectrometer:

[TABLE]

where $\dot{N}_{\mathrm{0}}$ is the amplitude of the electron signal and $\dot{N}_{\mathrm{b}}$ the observed background. The term $\theta_{\mathrm{max}}$ describes the largest pitch angle that can be transmitted according to eq. (3):

[TABLE]

This analytical method includes all relevant effects into the model (e. g. the transformation of the pitch angle resulting from adiabatic collimation), and allows us to determine the underlying distributions independently erhard:phd ; behrens:phd .

The asymmetric energy distribution is described by a generalized normal distribution GeneralizedGaussian1997 ,

[TABLE]

with the mean energy $\hat{E}$ and the energy width $\alpha_{E}$. For the transmission model, the energy distribution is evaluated in the range $E=[0;\infty)$. The asymmetry is described by the skewness parameter $\kappa$; at $\kappa=0$ the distribution is equivalent to a symmetric normal distribution. For $\kappa>0$, the function is limited to $E=[0;\hat{E}+\frac{\alpha_{E}}{\kappa})$. The width $\alpha_{E}$ can be converted into an energy spread $\sigma_{E}$, which can be compared independently of the skewness:

[TABLE]

The angular distribution is modeled by the sum of two normal distributions that are placed around $\theta=$$$,

[TABLE]

with the mean angle $\hat{\theta}$ and the angular spread $\sigma_{\theta}$. For the transmission model, the angular distribution is evaluated in the range \theta=[$$;$$]. The summing takes into account that the distribution is deformed for $\theta\rightarrow$$$ because the pitch angle is only defined for positive values.

The measured transmission functions presented in this paper have been fitted by a Markov-Chain Monte Carlo (MCMC) method of minimizing the $\chi^{2}$ value, using a code that was implemented in Python. It utilizes emcee emcee2013 for the MCMC fit process behrens:phd .

3 Measurements

The photoelectron source presented in this work has been commissioned successfully at the KATRIN monitor spectrometer. The corresponding measurements were carried out in the summer of 2014 and allowed us to verify the two key features of the electron source – angular selectivity and a small energy spread – and to study other important characteristics of the device. The monitor spectrometer was chosen because it could be operated independently of the main spectrometer during hardware preparations for its second commissioning phase. The electron source was subsequently mounted at the main spectrometer for the commissioning measurements of the spectrometer and detector section erhard:phd ; behrens:phd ; kraus:phd ; barrett:phd ; wierman:phd .

3.1 Experimental setup

In contrast to the main spectrometer, the monitor spectrometer features a symmetric magnetic field setup with $B_{\mathrm{max}}=$6\text{,}\mathrm{T} at the spectrometer entrance and exit. In our measurements, the spectrometer was operated at voltages $U_{\mathrm{ana}}\approx$-18.6\text{\,}\mathrm{k}\mathrm{V} with a minimal magnetic field $B_{\mathrm{min}}=$0.38\text{,}\mathrm{m}\mathrm{T} in the analyzing plane. The electron source was mounted in a vacuum chamber in front of the spectrometer beamline. The full setup is shown in figs. [3](#S3.F3) and [4](#S3.F4). The source was installed at a fixed position where the electron beam is always central to the spectrometer axis, as it was not necessary to perform measurements on different magnetic field lines. The emission spot at the back plate is located at $z_{\mathrm{es}}=$2.635\text{\,}\mathrm{m} in the spectrometer coordinate system, where the analyzing plane is located at $z_{\mathrm{ana}}=$0\text{,}\mathrm{m}. The two solenoid magnets are placed symmetrically around the analyzing plane at an axial distance of $z_{\mathrm{mag}}=\pm$2.01\text{\,}\mathrm{m}; the reference point $z_{0}=$0\text{,}\mathrm{m} refers to the analyzing plane at the spectrometer center. The magnetic field at the emission spot was measured with a Hall probe, yielding a value $B_{\mathrm{start}}=$21\text{\,}\mathrm{m}\mathrm{T}. To achieve electromagnetic conditions that are comparable to the main spectrometer setup, it is important to adjust $B_{\mathrm{start}}$ to achieve the magnetic field the electron source will be exposed to. Using an additional air-cooled coil that is placed close to the source at $z_{\mathrm{coil}}=$-2.640\text{,}\mathrm{m} and operated at a coil current $I_{\mathrm{coil}}=$35\text{\,}\mathrm{A}, the field at the emission spot increases to $B^{+}_{\mathrm{start}}=$27\text{,}\mathrm{m}\mathrm{T}$$. This is comparable to the main spectrometer setup where $29\text{\,}\mathrm{m}\mathrm{T}$ are reached.

Figure 5 shows the high voltage scheme of the monitor spectrometer setup. The electron source is connected with a small difference voltage to the high voltage of the spectrometer in order to cancel out voltage fluctuations that would occur if two independent power supplies were used. The back plate voltage, $U_{\mathrm{start}}$, can be varied against the spectrometer voltage $U_{\mathrm{spec}}$ by combining a power supply that operates at $0\text{\,}\mathrm{k}\mathrm{V}-1.25\text{\,}\mathrm{k}\mathrm{V}$ with a battery that delivers a voltage of about $90\text{\,}\mathrm{V}$. By putting the two voltage sources in series, it is possible to vary the starting voltage to achieve a surplus energy $q\Delta U=q(U_{\mathrm{start}}-U_{\mathrm{spec}})=$-90\text{,}\mathrm{e}\mathrm{V}1160\text{,}\mathrm{e}\mathrm{V}$$ without requiring a polarity-switching power supply. The voltage difference between electron source and spectrometer, $U_{\mathrm{start}}-U_{\mathrm{spec}}$, is measured by a difference voltmeter (DVM) to monitor the electron surplus energy. Transmission functions can be measured by varying the starting voltage within a few $\mathrm{V}$ around zero while observing the electron rate at the detector. The high voltage system is mainly located inside a Faraday cage, which is operated on the spectrometer high voltage. This cage is put inside another grounded HV cabinet to allow safe operation. The acceleration voltage for the front plate, $U_{\mathrm{acc}}$, is provided by an additional power supply that generates up to $5\text{\,}\mathrm{k}\mathrm{V}$ w. r. t. the back plate voltage. The acceleration voltage is thus kept constant while varying $U_{\mathrm{start}}$. This power supply is isolated for voltages up to $35\text{\,}\mathrm{k}\mathrm{V}$ and can be placed outside the HV cabinet.

3.2 Electron rate

As a first test that the electron source is operating as expected, the achieved electron rate at the detector was determined. The electron rate depends on the UV light setting. For the test measurements at the monitor spectrometer, we used a nominal laser setting of $f_{\mathrm{las}}=$100\text{,}\mathrm{k}\mathrm{H}\mathrm{z} and $I_{\mathrm{las}}=$6\text{\,}\mathrm{A} with a 1% ND filter in the optical beamline. This yields an electron rate of about $\dot{N}=$1500\text{,}\mathrm{c}\mathrm{p}\mathrm{s} at the detector in full transmission. The statistical uncertainty of a $10\text{\,}\mathrm{s}$ measurement is $<$1\text{\,}\% in this case, which is sufficient for our investigations and allows us to measure a typical transmission function in less than 10 minutes. For some measurements, the light source was switched to LEDs to allow wavelength-dependent measurements. Because the electron emission is influenced by the photon energy (UV wavelength) and the work function of the photocathode material, the electron rate that is achieved with LEDs varies in a typical range of $200\text{\,}\mathrm{c}\mathrm{p}\mathrm{s}1000\text{\,}\mathrm{c}\mathrm{p}\mathrm{s}$. The rate can be tuned over a wide range by changing the duty cycle of the function generator (pulser) that drives the LEDs. The LEDs were typically operated at a pulse frequency $f_{\mathrm{LED}}=$100\text{,}\mathrm{k}\mathrm{H}\mathrm{z}, a pulse width $\tau_{\mathrm{LED}}=$1000\text{\,}\mathrm{n}\mathrm{s} (10% duty cycle), and a forward voltage $U_{\mathrm{LED}}=$8.5\text{,}\mathrm{V}$$.

The dipole electrode in front of the electron source is intended to remove stored electrons from the beamline between source and spectrometer. The removal efficiency depends on the strength of the induced $E\times B$ drift eq. (9), and thus increases with a larger dipole voltage222The alternative solution of reducing the global magnetic field would change the magnetic fields at the MAC-E filter and is thus disfavored.. In our setup, the magnetic field at the electrode is $B_{\mathrm{dip}}=$78\text{,}\mathrm{m}\mathrm{T} with $E_{\mathrm{dip}}\approx$40\text{\,}\mathrm{k}\mathrm{V}\mathrm{/}\mathrm{m} according to simulations. The removal efficiency of the dipole electrode was investigated by measuring transmission functions in direction of increasing and decreasing electron surplus energy.

Figure 6 shows that the observed transmission function is affected by a hysteresis effect that depends on the dipole voltage, which allows investigating the removal efficiency of the dipole electrode. The observed transmission functions show a similar behavior, except for the nominal electron rate that is reached at full transmission. The small rate drift that can be observed in the upper panel can be explained by fluctuations in UV light intensity. The hysteresis effect can be explained by the continuous filling of the trap from the beginning of the measurement when measuring in direction of increasing surplus energy, because the surplus energy at the beginning is too small for electrons to be transmitted. Electrons with a given energy stay trapped until they lost kinetic energy (e. g. through synchrotron radiation) or are removed by the dipole field. Scattering processes with electrons of higher kinetic energy that are generated at a later time during the measurement cause some of the trapped electrons to gain kinetic energy, thereby increasing transmission probability towards the detector. The effect does not occur when the measurement is performed in inverse direction, where the higher-energetic electrons are transmitted at the beginning of the measurement wierman:phd . This leads to a hysteresis effect in the electron rate between the two scanning directions, which becomes smaller when the dipole voltage is increased and more electrons are removed from the trap. The observed rate difference is therefore a direct measure for the dipole efficiency. Our measurement indicates that a dipole voltage of $U_{\mathrm{dip}}=$2\text{,}\mathrm{k}\mathrm{V}$$ is sufficient to avoid the hysteresis effect. With lower dipole voltages, the observed rate difference between the two scanning directions increases, indicating an insufficient removal of stored electrons.

3.3 Energy spread

The energy resolution of a MAC-E filter eq. (4) depends on the retarding potential $U_{\mathrm{ana}}$. At low voltages $|U_{\mathrm{ana}}|\ll$18.6\text{,}\mathrm{k}\mathrm{V} and low electron energies $E\approx qU_{ana}$, the energy resolution improves because of the smaller amount of transversal energy left in the analyzing plane. A low voltage measurement with $U_{\mathrm{ana}}\approx$-200\text{\,}\mathrm{V} allows us to directly determine the energy distribution of the produced electrons. Unfortunately, at the monitor spectrometer it is not possible to detect electrons with $E\ll$10\text{,}\mathrm{k}\mathrm{e}\mathrm{V}$$ due to the energy threshold of the detector. Fortunately, the energy distribution can also be determined from a measurement performed at nominal high voltage.

If the source is operated at the so-called zero angle setting, it produces the smallest possible pitch angle in the spectrometer entrance magnet. The zero angle position has to be found manually by varying the plate angle $\alpha_{\mathrm{p}}$ around [math] independently for the vertical and horizontal axis. At an electron surplus energy $E=q\Delta U\approx$0\text{,}\mathrm{e}\mathrm{V} eq. ([6](#S2.E6)), the transmission probability is entirely dominated by the pitch angle of the emitted electrons. The observed electron rate is thus sensitive to small changes of the produced pitch angle, and the rate dependency w. r. t. the plate angle shows a maximum at zero angle $\alpha_{\mathrm{p}}=\alpha_{0}$. At the monitor spectrometer, the zero angle offset was found to be $\alpha_{\mathrm{hor}}=$0.04(1)\text{\,}\mathrm{\SIUnitSymbolDegree} and $\alpha_{\mathrm{ver}}=$1.13(1)\text{,}\mathrm{\SIUnitSymbolDegree} at $U_{\mathrm{dip}}=$2\text{\,}\mathrm{k}\mathrm{V}. This offset is caused by mechanical imperfections, which result in a minor misalignment that can be easily corrected by such a measurement. The impact of the angular spread on the observed transmission function is marginal when the zero angle is applied. In this case the actual mean and width of the angular distribution are not relevant to the analytical transmission model as long as \theta<$$$, and the energy spread dominates the shape of the resulting transmission function. It is thus possible to fit an (integrated) energy distribution to the measured transmission function while assuming a fixed angular distribution at a small pitch angle. For the case discussed here, an angular distribution with mean angle \hat{\theta}= and angular spread $\sigma_{\theta}= was used. These values are consistent with particle-tracking simulations (section 4.4) and complementary measurements of the angular distribution that were performed at the monitor spectrometer (section 3.5).

The measurements discussed in this section use the analytical transmission model eq. (10) with five free parameters: the amplitude and background of the electron signal, as well as the mean, width and shape of the energy distribution. The statistical uncertainty at each data point is derived from the measured rate fluctuations by computing the median- and $1\sigma$-percentiles of the rate taken at $2\text{\,}\mathrm{s}$ intervals for each data run at a fixed value of $U_{\mathrm{start}}$ (constant surplus energy). In most cases, the uncertainty determined by this method matches the $\sqrt{N}$ expectation from Poisson statistics. However, the percentile method is believed to be more robust against asymmetric rate fluctuations, and is thus preferred. For the transmission function measurements, an uncertainty of $\pm$60\text{,}\mathrm{m}\mathrm{e}\mathrm{V}$$ is assumed for the surplus energy, which is included in the fit as an additional term in the uncertainty of each data point. The value has been estimated from the contributions of the individual power supplies that are used in the setup erhard:phd .

Figure 7 shows transmission functions that were measured using the UV laser ($266\text{\,}\mathrm{n}\mathrm{m}$) and UV LEDs ($272\text{\,}\mathrm{n}\mathrm{m}302\text{\,}\mathrm{n}\mathrm{m}$) at $U_{\mathrm{ana}}=$-18.6\text{,}\mathrm{k}\mathrm{V} and $B_{\mathrm{min}}=$0.38\text{\,}\mathrm{m}\mathrm{T}. The laser measurement was performed twice and produced consistent results in terms of the corresponding energy distribution. All measurements used the $U_{dip}=$2\text{,}\mathrm{k}\mathrm{V} setting that was determined earlier. The observed transmission functions have been normalized to show a transmission probability with an average background of [math] and an average nominal amplitude of $1$ in the plot. Statistical fluctuations can yield negative amplitudes, as seen in the $302\text{\,}\mathrm{n}\mathrm{m}$ measurement. This normalization procedure was also applied in the subsequent measurements. When the UV-LEDs are used as light sources, the observed width of the transmission function decreases for larger wavelengths. This matches the expectation that the photoelectrons emitted from the photocathode material have a lower energy spread when the wavelength is closer to the work function of the photocathode, according to eq. ([7](#S2.E7)). The energy distributions that were determined from the observed transmission functions at different wavelengths are shown in the lower panel of fig. [7](#S3.F7). The reduced energy spread for increasing UV wavelengths is clearly visible. The asymmetric shape of the energy distribution that is observed here is expected from the theoretical model of the photoeffect [Photoemission1964a ](#bib.bib35). At lower photon energies, the low-energy fraction of the underlying energy distribution is cut off at $E=$0\text{\,}\mathrm{e}\mathrm{V}, which results in a more symmetric shape of the observed distribution.

Table 1 lists the parameters of the energy distribution, which are derived from the measured transmission functions. The upper limit of the energy distribution, $E_{\mathrm{max}}$, is used as an indicator for the achieved energy spread. It corresponds to the energy where the distribution drops to 25% of its maximum. Using this definition, the range $[0;E_{\mathrm{max}}]$ includes at least 90% of the distribution’s integral. It is also possible to determine the work function directly from a transmission measurement by relating the measured value of $E_{\mathrm{max}}$ to the known wavelength $\lambda$. This approach is discussed in section 3.6 below. The value $\sigma_{E}$ refers to the width of a symmetric normal distribution, which can be derived from the generalized normal distribution. The transformation to $\sigma_{E}$ takes into account the asymmetry of the distribution and allows comparing distributions with different asymmetry.

The results indicate that owing to the small angular spread in this setting, the width of the measured transmission function is fully dominated by the energy distribution of the electrons. This is true especially for measurements with zero angle and small wavelengths, where the angular distribution has only a minor effect on the transmission function and the energy spread is comparably large.

3.4 Magnetic reflection

Magnetic reflection occurs when the electron pitch angle reaches [math] and the total kinetic energy is in the transversal component. The pitch angle increases from the source towards $B_{\mathrm{max}}=$6\text{,}\mathrm{T} at the spectrometer entrance magnet as a result of adiabatic transformation eq. ([2](#S2.E2)). Magnetic reflection can be investigated by increasing the plate angle, $\alpha_{\mathrm{p}}$, until a rate decrease is observed at the detector. To ensure that electrons are reflected only magnetically and not because of an insufficient surplus energy eq. ([3](#S2.E3)), the measurement is performed at large surplus energies $q\Delta U\geq$10\text{\,}\mathrm{e}\mathrm{V}. The rate gradually decreases with increasing $\alpha_{\mathrm{p}}$ as more electrons are reflected as a result of the angular distribution in the magnet. The rate dependency can be modeled by a symmetric error function, which allows us to investigate the angular distribution at large pitch angles \theta\rightarrow$$$. The center position of the error function is referred to as *reflection angle* \alpha_{\mathrm{max}}$; it corresponds to the plate angle where 50% of electrons are reflected.

At the monitor spectrometer, this measurement was performed at four different azimuthal directions of the plate angle to investigate possible asymmetries, $\alpha_{az}=$0\text{,}\mathrm{\SIUnitSymbolDegree}90\text{,}\mathrm{\SIUnitSymbolDegree}180\text{,}\mathrm{\SIUnitSymbolDegree}270\text{,}\mathrm{\SIUnitSymbolDegree}$$. The results are shown in fig. 8 (solid lines). The underlying gauss-curves shown at the bottom of the figure allows a better comparison of the angular distributions.

As will be shown later in section 4.4, the produced pitch angle $\theta$ increases non-linearly with the plate angle $\alpha_{\mathrm{p}}$. This results from the adiabatic transformation towards the spectrometer entrance magnet:

[TABLE]

where $k$ is a scaling factor that depends on the non-adiabatic acceleration of the emitted electrons, and $B_{\mathrm{start}}$, $B_{\mathrm{max}}$ are the magnetic fields at the electron source and the spectrometer entrance, respectively.

Table 2 shows the fit results of these measurements. For nominal magnetic field at the electron source ($B_{\mathrm{start}}=$27\text{,}\mathrm{m}\mathrm{T}$$), reflection occurs at a plate angle \alpha_{\mathrm{max}}\approx$$$. The width of the angular distribution is consistent over the four measurements, yielding a width of \sigma_{\alpha}= for the underlying Gaussian distribution. The adiabatic transformation eq. ([15](#S3.E15)) converts the value $\sigma_{\alpha}$ to an effective angular spread $\sigma_{\theta}$ in the magnet. The conversion employs the constraint that magnetic reflection occurs at $\alpha_{\mathrm{p}}=\alpha_{\mathrm{max}}$ with $\theta=. This yields an average angular spread of $\sigma_{\theta}=$$$ at the maximal pitch angle of [math]. Note that the angular spread close to magnetic reflection increases because of the non-linearity of eq. (15), and is significantly lower at smaller pitch angles.

The discrepancy between the measurements in four azimuthal directions can be explained in two ways. Firstly, particle-tracking simulations indicate that misalignments of the emission spot relative to the plate setup of the electron source result in significant offsets of the produced pitch angles. Such misalignments can result from mechanical imperfections of the setup and are likely the explanation for the observed asymmetry behrens:phd . Secondly, phase effects can affect the electron acceleration processes in the source. The cyclotron phase of the emitted electrons differs depending on the azimuthal direction into which the electron beam is collimated. This results in slight variations of the produced pitch angle, which depend on the azimuthal plate angle $\alpha_{\mathrm{az}}$. The asymmetry in vertical direction ($\alpha_{\mathrm{az}}=$0\text{,}\mathrm{\SIUnitSymbolDegree}180\text{,}\mathrm{\SIUnitSymbolDegree}$$) is further increased by the electric field of the dipole electrode.

The magnetic reflection measurements have been fitted with MINUIT2 MINUIT1975 , using a normal distribution in integral form (scaled error function) to model the shape of the reflection curve. The uncertainty of the electron rate has been determined like explained above. An uncertainty of [math] is assumed for the plate angle and included in the fit; this value corresponds to the uncertainty of the plate angle read-out at the source (section 2.3).

3.5 Angular selectivity

When the plate angle at the source is increased, a larger pitch angle relative to the magnetic field vector is imprinted on the emitted electrons. The pitch angle in the analyzing plane of the spectrometer is a result of the adiabatic transformation eq. (2) and thus depends on the produced pitch angle at the source and the magnetic field variation between source and analyzing plane. The increased pitch angle leads to a shift of the measured transmission function to higher surplus energies, as the transversal component of the kinetic energy is larger in this case and needs to be compensated for. Such a shift is only observable if the electron source can produce large-pitch angles with a small angular spread, referred to as angular selectivity. The measured shift between the minimal pitch angle (zero angle \theta\approx$$$) and the maximal pitch angle (\theta=$$$ in the entrance magnet) allows us to determine the energy resolution eq. (4) of the spectrometer.

Figure 9 shows measured transmission functions at different plate angles $\alpha_{\mathrm{p}}$. The zero angle setting \alpha_{\mathrm{p}}=$$$ is used as a reference for the other measurements. The transmission functions are clearly separated and the expected shift to larger surplus energies is observed when increasing the plate angle. Table [3](#S3.T3) shows the corresponding parameters of the derived angular distribution. Magnetic reflection occurs at \alpha_{\mathrm{p}}\geq (as expected from the magnetic reflection measurement, which yields $\alpha_{\mathrm{max}}=; cmp. section 3.4). This results in a significantly deformed angular distribution, because reflected electrons are missing from the observed transmission function. Because the fit is based on a reference measurement at \alpha_{\mathrm{p}}=$$$ to obtain the corresponding energy distribution, the deformation affects the fit result and explains the large \chi^{2}$ value.

The transmission functions were fitted as explained above. However in this case, free parameters were the amplitude and background of the electron signal and the mean angle and the angular spread. This allows us to determine the produced pitch angle directly from the measurement, while assuming a known energy distribution of the electron source. In this case, a reference measurement at nominal settings (\alpha_{\mathrm{p}}=$$$, \lambda=266\text{\,}\mathrm{n}\mathrm{m}$$, U_{\mathrm{ana}}=-18.6\text{\,}\mathrm{k}\mathrm{V}$$) was used for the energy distribution (section [3.3](#S3.SS3)). The fit using the analytical model of the transmission function is not very sensitive to the actual shape of the angular distribution for \theta\rightarrow and $\theta\rightarrow. The angular distribution determined from the measurements at \alpha_{\mathrm{p}}=$$$ and \alpha_{\mathrm{p}}=\alpha_{\mathrm{max}}$ thus yield large uncertainties, and the angular spread is significantly smaller than at intermediate pitch angles. However, the fit results match expectations from an analytical calculation of the pitch angle based on the magnetic reflection limit discussed in section 3.4). The observed pitch angles are also confirmed by simulation results (section 4.4).

The total shift between \alpha_{\mathrm{p}}=$$$ and \alpha_{\mathrm{p}}=$$$ corresponds to the maximal difference of pitch angles in the spectrometer entrance magnet, $\theta=$0\text{,}\mathrm{\SIUnitSymbolDegree}90\text{,}\mathrm{\SIUnitSymbolDegree}, and is thus equivalent to the energy resolution of the spectrometer eq. ([4](#S2.E4)). The observed shift of $\Delta E_{\theta}=$1.20(6)\text{\,}\mathrm{e}\mathrm{V} corresponds to the expected energy resolution of

[TABLE]

for the monitor spectrometer operating at $B_{\mathrm{min}}=$0.38\text{,}\mathrm{m}\mathrm{T} at the spectrometer’s center and $B_{\mathrm{max}}=$6\text{\,}\mathrm{T}.

3.6 Work function

The measured energy spread (section 3.3) depends on the UV wavelength (photon energy) and the work function of the photocathode material eq. (7). The upper limit of the energy distribution is given by $E_{\mathrm{max}}=h\nu-\Phi$, and thus the work function can be determined from an energy distribution measurement. In addition, a direct measurement of the work function is possible by the method conceived by Fowler Photoemission1931 . Here the electron yield $I$ is measured at varying UV wavelengths $\lambda$, and the work function $\Phi$ can be determined by fitting the Fowler function to the data,

[TABLE]

with $\mu=(h\nu-\Phi)/(k_{B}T)$. Here $k_{B}$ is Boltzmann’s constant and $T$ the temperature of the photocathode. In comparison to alternate methods such as using a Kelvin probe KelvinProbe1991 , this in situ measurement allows us to determine the actual work function of the photocathode under nominal conditions at the experimental site. The determined work function thus can be compared with the measured energy distributions of the electron source.

Figure 10 shows the result of a Fowler-type measurement that has been performed at the monitor spectrometer setup. Six different UV LEDs have been used in combination with a monochromator to scan a wavelength range between $261\text{\,}\mathrm{n}\mathrm{m}$ and $321\text{\,}\mathrm{n}\mathrm{m}$. The emission profiles of the used LEDs can be approximated by a Gaussian as shown at the bottom of the plot. The achieved electron rate drops significantly when moving away from the LED’s peak wavelength, and the lower rate in these cases is compensated by an increased measurement time to reach a similar statistical uncertainty. For wavelengths $\lambda\gtrsim$270\text{,}\mathrm{n}\mathrm{m}, the observed electron yield matches the expectation from the Fowler equation eq. ([17](#S3.E17)). For smaller wavelengths a deviation from the expected behavior is observed, which gets more emphasized for decreasing wavelengths. At $\lambda<$270\text{\,}\mathrm{n}\mathrm{m}, the electron yield reduces again, which is in contrast to the well-known three-step model of photoemission developed by Berglund and Spicer Photoemission1964a ; Photoemission1964b . This observation can be explained by considering effects that become more dominant at higher photon energies Photoemission1995 . Close to the work function threshold the photoemission is dominated by excitons, which are created from incident light and move towards the surface where they release their electron (exciton dissociation). At photon energies well above the threshold (here $\gtrsim$4.5\text{,}\mathrm{e}\mathrm{V}$$), the photoabsorption becomes dominated by electron-excitation into unbound states. These electrons are not emitted from the surface, and the total electron yield is thus reduced at smaller wavelengths.

The fit results in a work function of $\Phi=$3.78\text{,}\mathrm{e}\mathrm{V}$\pm$0.03\text{,}\mathrm{e}\mathrm{V}$\mathrm{(sys)}\pm$0.01\text{,}\mathrm{e}\mathrm{V}$\mathrm{(stat)}$ for the gold photocathode with $20\text{\,}\mathrm{n}\mathrm{m}$ layer thickness at $T=$300\text{,}\mathrm{K}. The systematic uncertainty is estimated to $0.03\text{\,}\mathrm{e}\mathrm{V}$ from the uncertainty of the wavelength caused by to the filter width of the monochromator (0.01\text{,}\mathrm{e}\mathrm{V}$\hat{=}$1\text{,}\mathrm{n}\mathrm{m}$$), and the uncertainty of the LED peak wavelength ($2\text{\,}\mathrm{n}\mathrm{m}$). The monochromator was calibrated beforehand using the known wavelength of the UV laser ($266\text{\,}\mathrm{n}\mathrm{m}$). The determined work function is equivalent to a wavelength of $\lambda_{\mathrm{opt}}=$328.2\text{,}\mathrm{n}\mathrm{m}$\pm$2.3\text{,}\mathrm{n}\mathrm{m}$\mathrm{(sys)}\pm$0.7\text{,}\mathrm{n}\mathrm{m}$\mathrm{(stat)}$. The energy spread of the electron source can be minimized by matching the UV wavelength to this value. Unfortunately, the available LEDs limit the usable wavelength range to about $320\text{\,}\mathrm{n}\mathrm{m}$, as the very low rate at larger wavelengths would require unfeasibly long measurement times. However, even at wavelengths well above $266\text{\,}\mathrm{n}\mathrm{m}$, the optimum for maximal intensity, the resulting energy spread of $0.3\text{\,}\mathrm{e}\mathrm{V}$ or less is sufficiently small to determine the transmission properties of the spectrometer (section 3.3).

The work function determined from this measurement can be compared with the result from investigating the energy distribution. The maximal kinetic energy $E_{e,\mathrm{max}}$ of the photo-electrons emitted by the electron source is given by the relation

[TABLE]

with Planck’s constant $h$, the speed of light $c$ and the UV wavelength $\lambda$. The work function $\Phi$ can thus be determined from the upper limit of the energy distribution of the photo-electrons, which is given in table 1 (cmp. figure 7). The upper limit shifts to lower values when the UV wavelength is increased and the incident photons have less energy (i. e. the distribution gets narrower). The resulting work functions from this method should be consistent for measurements performed at different wavelengths. Table 4 shows the results from using this approach, and compares the determined work function $\Phi^{\dagger}$ with the work function $\Phi$ yielded by the Fowler-type measurement above. A combined analysis of the resulting work functions yields $\overline{\Phi^{\dagger}}=$3.810(1)\text{,}\mathrm{e}\mathrm{V}$$, using a weighted average that takes into account the uncertainties of $E_{e,\mathrm{max}}$. This result is consistent with the value determined by the Fowler-type measurement. It is thus verified that both methods produce consistent results, and that the determined work function is applicable to the measured transmission functions.

Our determined value for the work function of a gold surface is far below the theoretical expectation of $4.2\text{\,}\mathrm{e}\mathrm{V}5.1\text{\,}\mathrm{e}\mathrm{V}$ WorkFunction1970 ; WorkFunction1995 . However, the work function of our photocathode is affected by various effects like the surface roughness, impurities in the material and electric fields at the surface WorkFunctionEffects1973 ; WorkFunctionEffects2005 . These effects are largely eliminated when work functions are determined under ideal laboratory conditions, which makes these results incomparable to our in situ measurement. This observation is confirmed by other measurements performed at the KATRIN main spectrometer behrens:phd and in a test setup at WWU Münster zacher:phd ; winzen:diploma .

4 Simulations

The particle-tracking software Kassiopeia was developed as a joint effort from members of the KATRIN collaboration to simulate trajectories of charged particles such as electrons or ions in complex electromagnetic fields with very high precision Kassiopeia2016 . Kassiopeia is embedded in the so-called KASPER framework, the overall KATRIN software package. The software is used to study the transmission properties of the KATRIN spectrometers and to investigate background processes, among other simulation tasks. For the development of the electron source presented in this paper, Kassiopeia simulations provided substantial input for optimizations of the existing design. Detailed simulations were performed to investigate the electron acceleration processes within the source and to understand how the well-defined pitch angles are produced.

4.1 Implementation into Kassiopeia

Kassiopeia performs tracking of charged particles in electromagnetic fields based upon a given simulation geometry. Electric fields are computed by the boundary element method (BEM) from a set of charge densities at the electrode surfaces. The charge densities are pre-computed from the given electrode potentials with the iterative Robin Hood method RobinHood2012 . For axially symmetric electric fields, an approximation method known as zonal harmonic expansion can be used to speed up the field computations with negligible loss of accuracy Elcd2011 . To accurately model the electron source with all relevant components (e. g. the half-shell dipole electrode) it is necessary to use geometric shapes that break axial symmetry, thus no such approximation can be used. KEMField supports OpenGL-based graphics processing unit (GPU) acceleration, a feature that was utilized to considerably reduce the required computation time of such complex geometric structures. Magnetic fields are computed from a given set of coil geometries (solenoids and air coils) via elliptic integration; it is possible to apply zonal harmonic expansions here as well Magfield2011 . The simulations of the electron source use a detailed model of the magnet system at the monitor spectrometer. The particle-tracking in Kassiopeia is carried out by discretizing the trajectory into a finite number of steps. At each step the electromagnetic fields $\vec{E}(\vec{x}),\vec{B}(\vec{x})$ are evaluated and the equation of motion is solved by integration furse:phd ; groh:phd ; Kassiopeia2016 , after which the particle propagates to the next step. For charged particles, the Lorentz force eq. (8) defines the equation of motion.

The electrode geometry of the electron source was implemented in Kassiopeia based on CAD drawings of the electron source design. The position of the electrodes w. r. t. the spectrometer setup was determined from measurements at the experimental site and from comparisons of simulated with measured magnetic fields (section 3.1). Figure 11 shows the simulated magnetic field and electric potential between the photocathode of the electron source and the entrance magnet of the monitor spectrometer.

4.2 Energy and angular distributions

The simulations allow us to investigate the electron acceleration mechanisms inside the source. An important question is the effect of the electromagnetic fields on the energy and angular distributions achieved. Electrons were started from the emission spot on the back plate (radius $100\text{\,}\mathrm{\SIUnitSymbolMicro}\mathrm{m}$, according to the dimensions of the optical fiber in the experimental setup), where the starting voltage $U_{\mathrm{start}}\approx$-18.6\text{,}\mathrm{k}\mathrm{V} is applied. The initial energy is normal-distributed in the range $0\text{\,}\mathrm{e}\mathrm{V}0.6\text{\,}\mathrm{e}\mathrm{V}$ ($\mu=\sigma=$0.2\text{\,}\mathrm{e}\mathrm{V}). The initial polar angle w. r. t. the back plate follows a $\cos\theta$-distribution in the range $0\text{\,}\mathrm{\SIUnitSymbolDegree}90\text{\,}\mathrm{\SIUnitSymbolDegree}$. The parameters of the energy distribution were chosen according to measurement results, which yield an energy spread of up to $0.3\text{\,}\mathrm{e}\mathrm{V}$ (section 3.3), while the angular distribution matches the results from Photoemission2002 .

In the simulations presented here, $1000$ electrons were created at the back plate for each setting. The electrons were tracked up to the spectrometer entrance magnet in order to determine the energy and angular distributions. Both distributions are key parameters for the analysis of transmission function measurements.

Figure 12 shows the correlation between the initial and the final energy distributions w. r. t. the initial pitch angle. The simulations used the zero-angle setting (\alpha_{\mathrm{p}}=$$$) at U_{\mathrm{acc}}=5\text{\,}\mathrm{k}\mathrm{V}$$ and U_{\mathrm{dip}}=3\text{\,}\mathrm{k}\mathrm{V}$$333Note that the simulations use a 3\text{,}\mathrm{k}\mathrm{V}$dipole setting instead of the$2\text{,}\mathrm{k}\mathrm{V}$setting used in the measurements discussed above, because they are intended to be comparable with later measurements carried out at the KATRIN main spectrometer.. The distributions are characterized by their median and the$1\sigma$-width, which are both computed using quantiles. The energy distribution in the magnet yields a median energy of$\hat{E}=0.24\text{\,}\mathrm{e}\mathrm{V}$$ with an asymmetric width of \sigma_{E}^{-}=0.14\text{\,}\mathrm{e}\mathrm{V}$$ and \sigma_{E}^{+}=0.16\text{\,}\mathrm{e}\mathrm{V}$$, which is equivalent to the initial energy distribution. The observed asymmetry results from excluding negative energies from the underlying normal distribution (E\geq0\text{\,}\mathrm{e}\mathrm{V}$$). It can be seen that the energy distribution is completely unaffected by the acceleration processes inside the electron source. The resulting distribution is shifted to larger energies by the electrostatic acceleration, E=E_{0}+qU_{\mathrm{start}}$, but consistent in width and shape. Further simulations showed that this is also true for different values of$U_{\mathrm{start}}\neq-18.6\text{\,}\mathrm{k}\mathrm{V}$$ and non-zero plate angles \alpha_{\mathrm{p}}>$$$. The measured energy distribution in the magnet (section 3.3) is therefore fully equivalent to the initial energy distribution at the photocathode. This allows investigating the energy spread of the generated electrons by transmission function measurements, and to determine the work function of the photocathode according to eq. (7) (section 3.6).

Figure 13 shows results of the same simulation, but here the correlation between the initial and final angular distributions is investigated w. r. t. the initial kinetic energy. The electron pitch angle is changed by the non-adiabatic acceleration and the subsequent adiabatic transport to the entrance magnet. While the initial pitch angles follow a $\cos\theta$-distribution with angles up to [math], the pitch angles in the magnet are narrowly distributed. As above, the distributions were analyzed by their median and width. The distribution in the magnet has a median pitch angle of \hat{\theta}=$$$ with an asymmetric width of \sigma_{\theta}^{-}= and $\sigma_{\theta}^{+}=. Here the asymmetry is caused by the fact that the pitch angle is limited to the range $0\text{\,}\mathrm{\SIUnitSymbolDegree}90\text{\,}\mathrm{\SIUnitSymbolDegree}$ by definition. Whenever the pitch angle would assume negative values resulting from adiabatic transformation, it is instead mirrored to a positive value. The observed distribution is therefore “wrapped” into the positive regime at \theta=$$$, and thus becomes asymmetric when this effect occurs. As indicated by the coloring in the figure, the kinetic energy of the electrons influences the produced angular distribution as well. Electrons with higher kinetic energies contribute more to the observed angular spread than low-energetic electrons. The same effect is also observed at larger plate angles \alpha_{\mathrm{p}}>$$$. This is explained by the efficiency of the non-adiabatic acceleration in the plate setup of the electron source, which is responsible for imprinting a well-defined pitch angle on the electrons. According to the Lorentz equation eq. (8), the electrostatic acceleration becomes less effective as the electron energy increases (cmp. section 2.2). Low-energetic electrons are therefore more strongly collimated, while for electrons with higher initial energies the observed angular spread increases. It is thus possible to further reduce the angular spread by tuning the electron source to produce a small energy spread, which can be achieved by matching the UV wavelength to the photocathode work function.

4.3 Electron acceleration and transport

The performance of the electron selectivity in the electron source can be assessed by an investigation of the pitch angle transformation for \alpha_{\mathrm{p}}>$$$. Figure [14](#S4.F14) shows the evolution of the pitch angles along the electron trajectory between the photocathode and the spectrometer entrance magnet. The produced pitch angle along the electron trajectory depends on the plate angle \alpha_{\mathrm{p}}$, as indicated by the color scheme. The initial pitch angles are quickly collimated into a narrow distribution. Already at a distance$d\lesssim1\text{\,}\mathrm{m}\mathrm{m}$$ from the photocathode, the electron beam reaches an angular spread of less than [math] for any given setting of \alpha_{\mathrm{p}}$.

Electrons that pass the front plate are further accelerated to their full kinetic energy $E=qU_{\mathrm{start}}$ inside the source cage and transported adiabatically towards the spectrometer magnet. Because transmitted electrons pass the dipole electrode only once, the electric dipole field has no significant influence on the pitch angle transformation. However, the stray electric field of the dipole electrode affects the electron acceleration process itself: because of the asymmetric dipole field, a vertical electric field gradient is generated inside the source cage. Depending on the cyclotron phase of the electrons (and thus, depending on the azimuthal plate angle $\alpha_{\mathrm{az}}$) the electrons are accelerated differently and the pitch angle changes accordingly. Simulations show that dipole voltages $U_{\mathrm{dip}}=$0\text{,}\mathrm{k}\mathrm{V}4\text{,}\mathrm{k}\mathrm{V}$$ lead to deviations of the pitch angle up to [math], an observation also made by corresponding measurements. The deviations can be corrected by an empirical determination of the zero angle (cmp. section 3). The pitch angle increases towards $B_{\mathrm{max}}$ as a result of adiabatic transformation. When the pitch angle exceeds $\theta_{\mathrm{max}}=$$$, electrons are magnetically reflected. These electrons can get stored between the photocathode and the entrance magnet and need to be removed by the dipole electrode to avoid a possible Penning discharge.

4.4 Production of well-defined pitch angles

Figure 15 shows the simulated angular distributions in the spectrometer entrance magnet that are produced by the angular-selective electron source. When $\alpha_{\mathrm{p}}$ is increased, the angular distribution is shifted towards larger pitch angles, while the angular spread and shape is not affected; for plate angles \alpha_{\mathrm{p}}\gtrsim$$$, a broadening is observed as \theta$approaches the [math] limit. In case of the$\alpha_{\mathrm{p}}= setting, only a low number of electrons reaches the magnet, as the majority is magnetically reflected. At $\alpha_{\mathrm{p}}=, the observed distribution is asymmetric because the pitch angle cannot reach negative values; hence, the negative part of the distribution is mirrored at $\theta=$$$.

The dependency of $\theta$ in the entrance magnet on the plate angle $\alpha_{\mathrm{p}}$ can be described by the adiabatic transformation eq. (2) of an electron that propagates towards the spectrometer entrance magnet where a magnetic field $B_{\mathrm{max}}$ is achieved. The electron starts from the photocathode with an effective starting angle $\theta_{\mathrm{start}}^{*}$ (cmp. table 5) at the initial magnetic field $B_{\mathrm{start}}$:

[TABLE]

One can assume a strictly linear dependency of the small effective starting angle $\theta_{\mathrm{start}}^{*}$ and the mechanical plate angle $\alpha_{\mathrm{p}}$, $\theta_{\mathrm{start}}^{*}=f(\alpha_{\mathrm{p}})\approx k\cdot\alpha_{\mathrm{p}}$, which is indicated in fig. 16. A relation between $\alpha_{\mathrm{p}}$ and the produced pitch angle $\theta$ can be derived by employing the approximation $\sin^{2}x\approx x^{2}$ for small $x$:

[TABLE]

Because the transformation is fully adiabatic, it depends only on the ratio of magnetic fields, $B_{\mathrm{start}}$ and $B_{\mathrm{max}}$. This follows from eq. (2) and eq. (1) with the kinetic energy $E_{0}=qU_{\mathrm{start}}=\operatorname{const.}$ The factor $k$ then describes the effect of the non-adiabatic acceleration in the electron source, which produces the effective initial pitch angle $\theta_{\mathrm{start}}^{*}$ at the end of the grounded source cage.

Table 5 lists the corresponding pitch angles and angular spreads, which correspond to the median and the $1\sigma$-width of the angular distributions. Again, the values were computed using percentiles. It should be noted that the pitch angle at \alpha_{\mathrm{p}}=$$$ is systematically larger because of the asymmetric shape of the distribution, which shifts the median to larger values. Similarly, the median at \alpha_{\mathrm{p}}= is systematically smaller due to the deformation of the angular distribution, which is caused by magnetic reflection. The angular spread is comparable over a wide range of plate angles with $\sigma_{\theta}\approx. The spread becomes significantly larger for \theta\rightarrow$$$ as a result of adiabatic transformation eq. ([21](#S4.E21)). The simulated pitch angles and the angular spread are in good agreement with the corresponding measurements (section [3.5](#S3.SS5)). Table [5](#S4.T5) also lists the measured pitch angles \theta_{\mathrm{meas}}(cmp. tab. [3](#S3.T3)), and shows that both results are typically in agreement. An effective starting angle\theta_{\mathrm{start}}^{*}has been computed via eq. ([21](#S4.E21)), showing a strictly linear relation to the plate angle\alpha_{\mathrm{p}}$.

The measurements and simulations discussed in this work clearly show that the electron source achieves angular selectivity and can produce well-defined pitch angles with small angular spread. Figure 16 shows the produced pitch angle in the spectrometer entrance magnet (solid red line) and the pitch angle at the end of the source chamber (dashed blue line) according to simulations. At \alpha_{\mathrm{max}}=$$$, the pitch angle reaches [math] and magnetic reflection occurs. Resulting from the finite angular spread, at \alpha_{\mathrm{p}}= a fraction of the produced electrons is already reflected and cut off from the observed angular distribution. The simulated reflection angle is in excellent agreement with the magnetic reflection measurement (section [3.4](#S3.SS4)), where $\alpha_{\mathrm{max}}=$10.06(3)\text{\,}\mathrm{\SIUnitSymbolDegree}$$ was observed for $\alpha_{az}=. The effective initial pitch angle shows a strictly linear dependency to the plate angle with a factor $k$ and a constant angular spread $\sigma_{\theta,\mathrm{start}}=$$$.

5 Conclusion

An angular-selective electron source has been developed for the commissioning measurements of the KATRIN main spectrometer. In the first major measurement campaign at the KATRIN main spectrometer, several design improvements could be identified. After their implementation, the electron source was tested successfully at the monitor spectrometer in 2014. These preparation measurements demonstrated that the design requirements are completely fulfilled and that the electron source achieves all key features:

• •

Angular selectivity: The source produces well-defined electron pitch angles in the spectrometer entrance magnet. Magnetic reflection occurs when the pitch angle exceeds [math], which was observed at a plate angle of [math] in measurements. This value is in excellent agreement with the corresponding simulations, which also yield a reflection angle of [math].

• •

Small energy spread: Depending on the wavelength of the used UV light source, an energy spread between $0.09(7)\text{\,}\mathrm{e}\mathrm{V}$ at $302\text{\,}\mathrm{n}\mathrm{m}$ and $0.031(5)\text{\,}\mathrm{e}\mathrm{V}$ at $266\text{\,}\mathrm{n}\mathrm{m}$ was observed in transmission function measurements at $U_{\mathrm{start}}=$-18.6\text{,}\mathrm{k}\mathrm{V}. A measurement at low voltage $U_{\mathrm{spec}}\approx$-200\text{\,}\mathrm{V} allows us to determine the energy spread with much higher precision because of the improved energy resolution of the spectrometer. While this feature cannot be applied at the monitor spectrometer, it is of great use for the commissioning of the main spectrometer.

• •

Small angular spread: At the monitor spectrometer setup with $B_{\mathrm{max}}=$6\text{,}\mathrm{T} and $B_{\mathrm{start}}=$27\text{\,}\mathrm{m}\mathrm{T}, an angular spread of [math] or less was observed in transmission function measurements at different plate angles. Simulations indicate that the angular spread is typically even smaller (about [math]) for pitch angles $\theta\lesssim$$$.

• •

Electron rate: The electron source achieves a stable electron rate at the detector of $1500\text{\,}\mathrm{c}\mathrm{p}\mathrm{s}$ with the laser, and up to $400\text{\,}\mathrm{c}\mathrm{p}\mathrm{s}$ with the LEDs as light source. It is possible to regulate the rate by tuning the intensity of the UV photon system, e. g. varying the pulse width of the LED pulser or by adjusting the laser diode current.

• •

Pulsed mode: The light sources were operated in pulsed mode during the monitor spectrometer measurements. The pulsed mode allows time-of-flight (ToF) measurements to characterize several properties of the MAC-E filter. The ToF mode plays an important role in the commissioning measurements of the main spectrometer.

The energy spread of the generated electrons depends on the work function of the photocathode and can be minimized by adjustung the UV wavelength to the properties of the utilized material. For our gold photocathode a work function of 3.78\text{\,}\mathrm{e}\mathrm{V}$\pm$0.03\text{\,}\mathrm{e}\mathrm{V}$\mathrm{(sys)}\pm$0.01\text{\,}\mathrm{e}\mathrm{V}$\mathrm{(stat)}$ was found. Our value was determined *in situ* and is considerably lower than the literature value for a clean gold surface that has been prepared under ultra-high vacuum conditions. This observation is explained by surface impurities from the continuous operation at $p\approx${10}^{-7}\text{\,}\mathrm{m}\mathrm{b}\mathrm{a}\mathrm{r}, where residual gas can be adsorbed onto the photocathode surface. Other effects, such as the unavoidable surface roughness and strong electric fields at the photocathode, can additionally lead to a reduction of the work function. We used two different methods to directly determine the work function of the photocathode. A Fowler-type measurement, which investigates the wavelength-dependent electron yield, and a direct investigation of the energy distribution of the photo-electrons, which is derived from transmission function measurements. One advantage of the latter method is that the work function can be determined without requiring a dedicated wavelength scan. We demonstrated that this alternate method produces comparable results.

Particle-tracking simulations were performed with the Kassiopeia software, providing vital input for the analysis of the measurements, and allow us to get a precise understanding of the electron acceleration processes in the electron source. The simulation results are typically in good agreement with the measurements. We showed that the energy distribution of the electrons in the spectrometer entrance magnet corresponds to to the initial energy distribution, while both distributions show the same width and shape in the simulations. It is thus possible to fully determine the electron energy spectrum by performing transmission function measurements with a MAC-E filter. The angular distribution in the spectrometer magnet results from the non-adiabatic acceleration of the emitted electrons in the plate setup of the electron source and the subsequent adiabatic transport towards the spectrometer entrance. The electron beam is collimated by the strong electric acceleration field at the photocathode and reaches an effective angular spread of roughly [math] when leaving the non-adiabatic acceleration region. According to simulations, an angular spread of [math] (up to [math] for $\theta\rightarrow$$$) is reached in the spectrometer magnet. This value is slightly lower than the measured angular spread of [math]. The produced pitch angle and the angular spread in the magnet strongly depend on the magnetic fields at the setup. The differences between measurements and simulations can therefore be explained by undetected misalignments of the setup and entailing inaccuracies of the computed fields.

Our electron source allows us to investigate major characteristics of a MAC-E filter, such as the transmission properties and the effective energy resolution of the spectrometer. We studied key features of the electron source in measurements at the KATRIN monitor spectrometer and in a suite of accompanying simulations. We fully characterized our electron source and demonstrated a reliable operation in a MAC-E filter setup. The electron source can be utilized as a vital tool for the commissioning of the KATRIN main spectrometer and in preparation of the upcoming neutrino mass measurements.

Acknowledgements.

This work has been supported by the Helmholtz Association (HGF) and the Bundesministerium für Bildung und Forschung (BMBF) with project numbers 05A14PMA and 05A14VK2.