Penning-trap measurement of the Q value of electron capture in 163Ho for the determination of the electron neutrino mass

Similar content being viewed by others

Main

The observation of neutrino flavour oscillations proves that neutrinos are massive particles, establishing that the weak neutrino flavour eigenstates are a superposition of three neutrino mass eigenstates, in contradiction to the standard model of particle physics1,2. Oscillation experiments can only investigate the differences between the squared neutrino mass eigenvalues, leaving the absolute scale of the neutrino mass an open question. Thus, the absolute scale of the neutrino mass remains one of the most sought-after quantities in nuclear and particle physics, cosmology and theories beyond the standard model that could potentially explain the origin of the neutrino rest mass3,4,5,6.

Neutrinos are produced in weak nuclear decays; a model-independent measurement of their rest mass can be performed in a kinematic study of the decay products, where the neutrino itself is not directly detected. Relying on energy and momentum conservation, this is currently the most model-independent approach for neutrino mass determinations. Kinematic investigations constrain the effective rest mass of the electron neutrino (or antineutrino) \({m}_{{\nu }_{e}}^{2}=\mathop{\sum }\nolimits_{i = 1}^{3}| {U}_{ei}{| }^{2}{m}_{i}^{2}\), where Ufi are the elements of the Pontecorvo-Maki-Nakagawa-Sakata matrix, which describes the superposition of mass eigenstates mi (i ∈ {1, 2, 3}) in the flavour eigenstates νf (f ∈ {e, μ, τ}). The individual mass eigenstates are not resolved in these experiments as the squared mass differences are well below current instrumental resolutions—the largest such difference is \(\Delta {m}_{32}^{2}=(2.453\pm 0.033)\times 1{0}^{-3}\,{{{{\rm{eV}}}}}^{2}\,{c}^{-4}\) (ref. 7).

The most stringent constraint on the neutrino mass scale comes from the analysis of the matter distribution in the Universe, which provides a limit on the sum of the neutrino masses of <120 meV c−2 (ref. 8). The most stringent direct limit of 0.8 eV c−2 (90% confidence level) from a kinematic study of tritium β-decay was set by the KATRIN collaboration9,10.

Complementary to this approach, there are several experiments using calorimetric techniques to investigate the neutrino rest mass directly. The first calorimetric approaches were the MANU and MIBETA experiments investigating 187Re β-decay, yielding upper limits of 19 and 15 eV c−2 (90% confidence level), respectively11. Two current experiments, ECHo12,13 and HOLMES14,15, are investigating the electron capture in 163Ho → 163Dy + νe + Ecal, where Ecal is the energy detected in a calorimeter. The current upper limit of the electron neutrino rest mass is on a level of <150 eV c−2 (ref. 13) and the ECHo and HOLMES collaborations aim to achieve sensitivities well below <1 eV c−2 (ref. 12).

Within the ECHo collaboration, metallic magnetic calorimeters are used for the measurement of the energy of all emitted radiation, except for the energy carried away by the neutrino. This is obtained by implanting 163Ho ions directly into the absorber material of the detector. The calorimetrically measured decay spectrum is subsequently analysed by fitting it to a theoretical spectral shape from which the Q value, as well as \({m}_{{\nu }_{e}}\), can be determined. To quantitatively investigate systematic effects in the interpretation of the calorimetrically measured spectra that might arise due to the 163Ho ions being implanted into a metallic material, this Q value is best compared to one obtained from an independent direct measurement. The required accuracy of ~1 eV c−2 can only be reached at present using high-precision Penning-trap mass spectrometry. In Penning-trap mass spectrometry, the Q value is addressed directly through a measurement of the mass difference of the mother and daughter nuclides, 163Ho and 163Dy, respectively12,16, by measuring the free-space cyclotron frequency ratio of the two species in a strong homogeneous magnetic field B. Within a magnetic field, an ion with charge-to-mass ratio q/m is forced into a circular orbit, in which it revolves with the free-space cyclotron frequency \({\nu }_{c}=\frac{1}{2\uppi }\frac{q}{m}B\). In a Penning trap, a superimposed weak quadrupolar electrostatic potential confines the ion along the magnetic field lines and modifies the ion's radial motion: the free-space cyclotron motion splits into the magnetron motion with a frequency ν− and the modified cyclotron motion with a frequency ν+. The quadrupolar electrostatic potential also induces a harmonic oscillatory motion with frequency νz along the magnetic field lines. From a measurement of all three motional eigenfrequencies, the free-space cyclotron frequency can be reconstructed using the invariance theorem \({\nu }_{c}^{2}={\nu }_{+}^{2}+{\nu }_{z}^{2}+{\nu }_{-}^{2}\) (ref. 17). From subsequent measurements of the free-space cyclotron frequency, the ratio \({R}_{q+}={\nu }_{c}\left({\,}\!^{163}{\rm{Dy}}^{q+}\right)/{\nu }_{c}\left({\,}\!^{163}{\rm{Ho}}^{q+}\right)\) is determined, which finally allows the Q value to be determined by including atomic physics calculations of the binding energy difference \(\Delta {E}_{{{{\rm{B}}}}}^{\;q+}\) of the removed electrons:

$$Q={m}_{{{{\rm{Dy}}}}}^{q+}\left({R}_{q+}-1\right)+\Delta {E}_{{{{\rm{B}}}}}^{\;q+}.$$

(1)

\(\Delta {E}_{{{{\rm{B}}}}}^{\;q+}\) is given by the difference in the sum of the binding energies of the n missing electrons in the highly charged ions (HCIs) of both nuclides, and \({m}_{{{{\rm{Dy}}}}}^{q+}\) is the 'reference' mass of the HCI of dysprosium. q = n ⋅ e is the charge of the ions, where e is the elementary charge and n the number of removed electrons (charge state). To enhance the readability of formulae, sometimes q also denotes the number of missing electrons n.

The Penning-trap experiment PENTATRAP

Experimental set-up

Measurements of R of the two HCIs (163Hoq+ and 163Dyq+) were carried out with the high-precision Penning-trap mass spectrometer PENTATRAP located at the Max-Planck-Institute for Nuclear Physics in Heidelberg, Germany18,19. An overview of the apparatus is given in Fig. 1.

The upper horizontal part of the beamline is located on the ground floor, while the superconducting magnet is located in a dedicated laboratory in the basement. The TIP-EBIT is an electron beam ion trap specifically designed for very small samples sizes21. Following the TIP-EBIT in the horizontal beamline, a Bradbury-Nielsen gate is used to separate a single charge state. HCIs produced in the TIP-EBIT are guided through the electrostatic beamline to the stack of five identical Penning traps in the superconducting magnet. Deceleration electrodes with appropriately timed voltage pulses are used to capture the HCIs in the Penning traps. A more detailed view of the Penning-trap stack is shown on the right.

HCIs of the synthetic radioisotope 163Ho, which was produced by neutron irradiation of stable 162Er (ref. 20), and HCIs of the stable 163Dy were produced in a compact room-temperature electron beam ion trap (EBIT) that was specifically designed and constructed for the production of HCIs from samples available only in limited quantities (TIP-EBIT)21. For the measurements reported here only 2 × 1015 atoms of 163Ho were used, with a typical sample containing about 1014 atoms of 163Ho. HCIs of the two species were extracted with a kinetic energy of 4.4 keV q−1 from the EBIT and transported through an electrostatic beamline towards the Penning traps. Individual charge states n = {38, 39, 40} were selected using a Bradbury-Nielsen Gate and a fast switching electronic circuit22,23 located about 1.5 m from the EBIT. Just before reaching the mass spectrometer, the HCIs were decelerated to a few electronvolts per q by appropriately timed voltage pulses on two cylindrical drift tubes.

The mass spectrometer consisted of a stack of five identical cylindrical Penning traps located in the cold bore of a 7 T actively shielded superconducting magnet18,24. The voltages applied to the Penning-trap electrodes were supplied by an ultrastable voltage source25. The Penning traps and the detection system were located inside a vacuum chamber immersed in liquid helium at a temperature of about 4 K. Two (trap 2 and trap 3; Fig. 2a) of the five Penning traps were equipped with a non-destructive image current detection system18,26,27,28 and were used to measure the motional frequencies of the ions. Trap 1 and trap 4 served as storage traps, while trap 1 was also used as a capture trap when a new set of ions was loaded into the trap stack.

a, Rendering of the stack of five identical cylindrical Penning traps of the PENTATRAP experiment. Traps 2 and 3, with labels marked in red, are used as measurement traps and are equipped with a detection system. Shuttling the ions from Position 1 into Position 2 effectively swaps the ion species in traps 2 and 3, resulting in the alternating data structure as shown in b. Traps 1 and 4 are used as storage traps, while trap 5 is not used in this measurement. b, Exemplary dataset of the measured νc of 163Ho (orange) and 163Dy (blue) in traps 2 (upper panel) and 3 (lower panel) for one measurement run in charge state q = 38 ⋅ e. For traps 2 and 3 frequency offsets of 25,081,589 Hz and 25,081,620 Hz were subtracted. The linear drift of the free-space cyclotron frequency can be attributed to the slow decay of the magnetic field of the superconducting magnet due to the flux creep effect30,31. Note that the y axis is broken for illustrative purposes and that the error bars are smaller than the data point markers. c, Ri of νc of 163Dy and 163Ho in traps 2 (filled circles) and 3 (empty circles) determined from the full dataset of two runs for the charge state n = 38. The data for each run are shown in separate panels; the ratios from b are shown in the left panel; the right panel shows ratios from the second dataset. The horizontal black line indicates the weighted average of all measured ratios for this charge state with the light red band marking the 1σ uncertainty band. Error bars correspond to the 1σ statistical uncertainty that is propagated by Gaussian uncertainty propagation.

Environmental parameters affecting the magnetic field in the traps were stabilized (such as the temperature in the laboratory to 0.1 K/day and the liquid helium level and pressure of helium gas inside the cold bore of the magnet). In these conditions, the magnetic field exhibited a relative drift of a few 10−10 per hour29. Frequency measurements were performed overnight and on weekends, when external perturbations were minimal.

The measurement started with loading a set of three ions in the order 163Dy, 163Ho and 163Dy into traps 2, 3 and 4, respectively (Fig. 2a). The motional frequencies of the HCIs in traps 2 and 3 were measured simultaneously, starting with the ions in Position 1. Subsequently, the ions were shuttled to Position 2, which effectively swapped the ion species in traps 2 and 3 (Fig. 2a) and the measurement was repeated. The resulting data structure is shown in Fig. 2b where νc is plotted as a function of the measurement time. Alternating data points for 163Dy and 163Ho result from the swapping of the ion species in traps 2 and 3. More details of the ion preparation and the measurement sequence are given in the Methods.

Data analysis

To extract frequency ratios Ri from νc, the magnetic field behaviour is interpolated between the individual frequency measurement data points from one species to the time when the other species' frequencies were measured.

Figure 2b shows exemplary the free-space cyclotron frequencies from one measurement run performed on ions with the charge state q = 38 ⋅ e. The linear slope of the data points can be attributed to the slow decay of the magnetic field of the superconducting magnet due to the flux creep effect30,31 and is on the order of a few 10−10 per hour relative to the absolute magnetic field of ~7 T.

In the data analysis, the frequency of 163Dy was linearly interpolated between two data points to the time at which 163Ho was measured. Ri was determined from this interpolated data point, as illustrated in Fig. 3a. This procedure was followed for the full dataset. Residual nonlinear behaviour of the cyclotron frequency drift, originating from physical effects that alter the temperature and position of magnetic materials that surround the Penning traps and change the magnetic field within the traps, was taken into account in the uncertainty of the interpolated Ri. For this, the frequency data points were interpolated back to themselves (Fig. 3b), and the sum of the residuals divided by the number of residuals was included as an additional uncertainty in the ratio. The resulting ratios \({R}_{i}={\nu }_{c,i}\left({\,}\!^{163}{{{{\rm{Dy}}}}}^{38+}\right)/{\nu }_{c,i}\left({\,}\!^{163}{{{{\rm{Ho}}}}}^{38+}\right)\) for the two measurement runs are shown in Fig. 2c for both traps. The ratios for the individual traps were consistent, therefore the final ratio was calculated as the weighted average and shown as a red line including the 1σ uncertainty band. For the calculation of the uncertainty of the final ratio, the inner error \({\sigma }_{{{{\rm{int}}}}}^{2}\) and the outer error \({\sigma }_{{{{\rm{ext}}}}}^{2}\) were calculated, and the larger of the two was used as the final uncertainty32,33:

$${\sigma }_{{{{\rm{int}}}}}^{2}=\frac{1}{{\sum }_{i}\frac{1}{{\sigma }_{i}^{2}}}$$

(2)

$${\sigma }_{{{{\rm{ext}}}}}^{2}=\frac{{\sum }_{i}\frac{1}{{\sigma }_{i}^{2}}{({R}_{i}-\tilde{R})}^{2}}{(N-1){\sum }_{i}\frac{1}{{\sigma }_{i}^{2}}}.$$

(3)

Here, Ri and σi are individual cyclotron frequency ratios and their corresponding 1σ uncertainties, \(\tilde{R}\) is the weighted average and N is the total number of ratios.

Detailed plot of the first few datapoints of the cyclotron frequency νc from Fig. 2b to illustrate the data analysis procedure. From the frequency values, an offset of 25,081,589 Hz is subtracted. a, Linear interpolation between two 163Dy datapoints to the time at which 163Ho was measured for the determination of Ri. Note that the vertical axis was broken for illustration purposes. b, Exemplary estimation of the nonlinearity by interpolation of the data onto itself. Here we linearly interpolate between the first and third datapoints and determine the difference between the measured datapoint between and the interpolated one. The sum of these 'residuals' divided by the number of residuals in the full dataset is taken into account as an additional uncertainty of the ratio. For details of the analysis procedure see the main text. Error bars correspond to the 1σ statistical uncertainty that is propagated by Gaussian uncertainty propagation.

The total systematic uncertainty (for example, field anharmonicities and inhomogeneity, image charge shift and relativistic shift) was strongly suppressed due to the fact that 163Ho and 163Dy in the same charge state form a unique mass doublet with a sufficiently small mass difference of about 2.8 keV. With a difference in the mass-to-charge ratio of only about 10−8, the same trapping potential was used for both 163Ho and 163Dy and the magnetron and axial frequencies were sufficiently equal. Thus, all systematic uncertainties in the free-space cyclotron frequency measurement cancelled to a large extent in the determination of R and were smaller than 10−12. Extended Data Table 1 summarizes the considered systematic shifts. An additional systematic uncertainty can arise from the fact that HCIs might have long-lived low-energy atomic metastable states, as observed in previous measurements34. This is undesirable as it would shift the determined Q value by the energy of the metastable state. In the 'Q value determination' section we compare the Q values resulting from the measurements of three different charge states, which allowed us to exclude potential shifts in the Q value due to long-lived electronic metastable states that would influence each charge state differently. The final ratios of free-space cyclotron frequencies of the ions in the different charge states are summarized in Table 2.

Calculation of the binding energy differences

Theoretical calculations provided the binding energies of the electrons removed from the neutral Ho and Dy atoms. The Dy and Ho atoms are in the [Xe]4f106s2 5I8 and [Xe]4f116s2 4I15/2 electronic states, respectively. For better control of systematic effects, several HCI of Dy and Ho were considered in the experiment, namely, Dy38+,39+,40+, with the ground states [Ar]3d10,9,8, respectively, and Ho38+,39+,40+ with ground states [Ar]3d104s, [Ar]3d10 and [Ar]3d9, respectively.

Configuration interaction method

In a first set of calculations, the binding energies were calculated in Quanty35,36,37 using the configuration interaction method. The starting point was a fully relativistic density functional theory calculation with the full-potential local-orbital minimum-basis code FPLO38,39,40. The density functional theory calculation determined the ground-state density of the ion around which a configuration interaction expansion was made. The corresponding Kohn-Sham orbitals were used as single-particle bases to construct the Slater determinants that span a configuration space. The Hamiltonian comprises Coulomb and static Breit interactions between the electrons, as well as their relativistic kinetic energies and potential energies due to Coulomb attraction of the ion's nucleus. Diagonalization of this Hamiltonian on a given configuration space using the Lánczos algorithm determined the ground-state energy of an ion.

At first, only the space of the ground-state configuration was considered. The configuration space was then iteratively expanded to include single, double and triple excitations of electrons into orbitals with higher principal quantum numbers. Details of these calculations are given in the Methods. We obtained the calculated binding energy differences given in Table 1.

MCDHF method

In the second set of calculations, we used the MCDHF method41 combined with Brillouin-Wigner many-body perturbation theory42,43.

In the MCDHF method, the atomic state function is modelled as a superposition of configuration state functions (CSFs) with fixed angular momentum, magnetic and parity quantum numbers. The CSFs are built as Slater determinants of Dirac orbitals in the jj coupling scheme. Using the parallel GRASP2018 codes43, we expanded the space of virtual orbitals used for the construction of CSFs by single- and double-electron exchanges in a systematic manner. The convergence of the energies with respect to the maximal principal quantum number of virtual orbitals was monitored, and the spread of values resulting from different correlation models was used as a measure of the leading contribution (90%) of the theoretical uncertainties. In case of the HCI, the set of CSFs was generated with exchanges including all occupied orbitals from 1s onwards, and with virtual orbitals up to typically 10h. Virtual orbitals were optimized in a layer-by-layer fashion43,44. The effects of the Breit interaction, recoil and approximate quantum electrodynamic corrections were accounted for by the configuration interaction method using orbitals from the MCDHF procedure43. More details are given in the Methods. We obtained the theoretical binding energy differences listed in Table 1.

In a third set of calculations, we used the MCDFGME code45 to check the previous results. The calculation was performed in the optimized level mode, where all correlation orbitals are fully relaxed instead of the layer-by-layer method. Convergence was much more difficult in this case and limited the number of extra orbitals that could be added in the evaluation of correlation. In this calculation, the magnetic and retardation interaction at the Breit level were included in the Dirac-Fock equations on the same footing as the Coulomb interaction, meaning that the Breit interaction was included to all orders in the correlation energy46. The Uehling potential was also evaluated to all orders47. Finally, self-energy screening was calculated using both the Welton method48 and the model operator method49. For the HCIs, energies obtained by exciting occupied orbitals from 3s or 3d to open shells (4f, 6p, 5d, 7s, 7p and 5g) were compared. For neutrals, values obtained by exciting the core from 3d and 4s were compared. Calculations included only single and double excitations, as triple excitations led to unmanageably large numbers of magnetic and retardation integrals. All possible single excitations were included, even those obeying Brillouin's theorem45. The results are given in Table 1 and are in good agreement with the GRASP2018 evaluation. Both sets of values are in agreement with the uncorrelated values50, confirming the good compensation of correlation between the two ions.

Final values for the binding energy difference

The final binding energy \(\Delta {E}_{{{{\rm{B}}}}}^{\;q+}\) for each charge state q was calculated as the weighted average of the values from the configuration interaction and MCDHF calculations (Table 1). The uncertainty was determined by comparing the inner and outer errors and using the larger as the final uncertainty on \(\Delta {E}_{{{{\rm{B}}}}}^{\;q+}\). For the charge states q = {39, 40} ⋅ e, the larger of the two uncertainties was averaged with the uncertainty assuming correlations between the configuration interaction and MCDHF methods; that is, with the uncertainty of 0.8 eV of the MCDHF method. The resulting \(\Delta {E}_{{{{\rm{B}}}}}^{\;q+}\) are consistent with the MCDFGME calculations described above, as well as with the calculations recently published in ref. 51.

Q value determination

The Q value of the electron capture in 163Ho was determined from the measured R (see 'The Penning-trap experiment') and the theoretically calculated binding energy differences \(\Delta {E}_{{{{\rm{B}}}}}^{\;q+}\) (see 'Calculation of the binding energy differences' and Table 1) for each charge state q = {38, 39, 40} ⋅ e according to equation (1). The (reference) mass \({m}_{{{{\rm{Dy}}}}}^{q+}\) of 163Dyq+ was calculated starting from the mass of atomic 163Dy, mDy (ref. 52) and subtracting the masses of the n missing electrons53,54 and their binding energies55. Table 2 lists R for the three measured charge states, as well as the 1σ uncertainty δR, which was computed using standard Gaussian uncertainty propagation.

Using equation (1) and the binding energy differences, the final Q values were calculated for the three charge states and are summarized in Table 2.

The resulting Q values agree within their 1σ uncertainties. Resulting from this very good agreement, systematic deviations from either the free-space cyclotron ratio measurement or from the calculation of the binding energy difference can largely be excluded. Furthermore, the influence of unknown metastable electronic states can also largely be ruled out as it is very unlikely that an electronic metastable state would have exactly the same excitation energy in all three of the measured charge states.

The final Q value was calculated as the weighted average of the Q values obtained for the three charge states resulting in:

$$Q=2{,}863.2\pm0.6\,{{{\rm{eV}}}}\,{c}^{-2}.$$

(4)

This value is in 1σ agreement with the previously measured value at SHIPTRAP of 2,833 ± 34 eV c−2 (ref. 56), but 50 times more precise. In Fig. 4 the most recent measurements of the Q value of 163Ho from cryogenic microcalorimetry, Penning-trap mass spectrometry and the Atomic Mass Evaluation (AME) 2020 (ref. 52) are shown. The value from AME 2020 is an average from three different microcalorimetric measurements. Our value is slightly higher than the current AME adjustment, and agrees within 1.2σ.

Data from Penning-trap mass spectrometry (SHIPTRAP56 and this work), microcalorimetry (ECHo 201757 and ECHo 201913) and the most recent AME adjustment (AME 2020)52 are shown. Error bars correspond to 1σ combined uncertainties.

The 163Ho electron capture Q value was obtained by combining a high-precision measurement of the free-space cyclotron frequency of HCIs of the mother and daughter nuclides in a Penning trap and precise atomic physics calculations of the electronic binding energies of the missing electrons. Experiments investigating the electron neutrino mass by microcalorimetric measurements of the decay spectrum of 163Ho such as those performed by the ECHo and HOLMES collaborations are now provided with an independently measured Q value with a precision of 0.6 eV, which allows the assessment of systematic uncertainties in the neutrino mass determination using cryogenic microcalorimetry on a level of <1 eV.

Methods

Measurement preparation and sequence

The measurement preparation started with loading a set of three ions, in the order 163Dy, 163Ho and 163Dy into traps 2, 3 and 4, respectively (Fig. 2a). Each HCI was first loaded into trap 2, where its motional amplitudes were reduced by resistive cooling17,58. Great care was taken to ensure that only a single HCI was captured in a trap and cooled. For this also the 'magnetron cleaning' technique was applied59. After being prepared in this way, the ions were moved to one of the following traps and stored until the set of ions for a given measurement was complete.

In both measurement traps (traps 2 and 3), the motional frequencies of the HCIs were measured using the single-dip, double-dip and pulse-and-phase techniques58,60. The magnetron frequency was small compared with the other two motional frequencies and depends only very weakly on the ion's mass; it was therefore measured only once a day using the double-dip technique before the main measurement sequence. Thus, the main measurement sequence reduced to a measurement of the modified cyclotron frequency (pulse-and-phase technique) and the axial frequency (double-dip technique), which were performed simultaneously in traps 2 and 3. Compared with a single-trap measurement, this effectively doubled the statistics and furthermore allowed different analysis methods to be employed, as well as systematic checks through comparing the results obtained in both traps.

In the pulse-and-phase cycle, the starting phase of the cyclotron motion was set by exciting it using a dipolar pulse with the frequency determined with the double-dip method during the preparation. The modified cyclotron motion then evolved freely during the phase evolution time Tevol (about 40 s) while the axial frequency was determined using a dip measurement19,29,61,62,63,64. Following the phase evolution time, the phase information that accumulated in the modified cyclotron motion was coupled to the axial motion using a π pulse on the sideband frequency and the final phase was measured with the image current detection system60. This is done in traps 2 and 3, starting with the ions in Position 1. Subsequently, the ions were shuttled to Position 2, which effectively swapped the ion species in traps 2 and 3 and the measurement was repeated (Fig. 2a). This sequence was repeated 24 times in one main measurement loop and can be continued, in principle, infinitely. Typically, the measurement is stopped due to either external magnetic field perturbations or charge exchange of the HCIs. Lifetimes of the HCIs until a charge exchange process happens were up to 36 hours. Reloading ions is beneficial as it allows one to compare different sets of ions, and therefore also systematic checks for contaminant ions that might be present in the Penning traps during the measurement or for possible metastable electronically excited states in the HCIs34.

Convergence studies with the configuration interaction and multiconfiguration Dirac-Fock methods

In the configuration interaction calculations made with the Quanty code, we iteratively expanded the configuration space with single, double and triple excitations into single-electron states with higher principal quantum numbers. For the ions this implies iterative inclusion of excitations into orbitals with n = 4, 5, 6, 7, 8, and for the neutral atom n = 5, 6.

The evolution of the ground-state energy with expanding configuration space was monitored and showed to good approximation 1/n behaviour that allowed extrapolation of the ground-state energy to estimate its uncertainty due to a truncated configuration space. Considering further uncertainties due to numerical accuracy, the choice of single-particle basis sets and triple excitations, we arrived at a total uncertainty of 1 eV for the estimation of the differences in binding energies of Hoq+ and Dyq+.

As consistency check, for every step where the configuration space was increased the binding energy difference between Hoq+ and Dyq+ was calculated. It showed an approximate 1/n2 behaviour, which again allowed extrapolation. Within our uncertainties, we obtained the same results as in Table 1.

In case of the multiconfiguration Dirac-Fock calculations with the GRASP2018 package, as described in the main text, we found that for neutral atoms, the inclusion of all spectroscopic orbitals into the active space would lead to several tens of millions of CSFs, which is not tractable at present. For these systems, we included exchanges from the 3s orbital up to typically 8h. To bridge the different models used for the neutral atoms and the HCIs, we also studied the intermediate Pd-like HCIs Dy20+ and Ho21+, with excitations from the 2s orbital to typically 10h. We observed that the correlation terms largely cancelled in energy differences such as [E(Ho) − E(Ho21+)] − [E(Dy) − E(Dy20+)] and [E(Ho21+) − E(Ho40+)] − [E(Dy20+) − E(Dy40+)] due to the structural similarities of nearby charge states. These differences converged more quickly when the set of virtual orbitals was extended than with the individual energies E(Ho) and E(Dy) of the neutral atoms. We note that this inclusion of an intermediate system was useful because of the high charge states 38+, 39+ and 40+ in the experiment, and allowed uncertainties to be reduced.

Data availability

Datasets generated during the current study are available from the corresponding author on reasonable request.

Code availability

Code used in the analysis of the data is available from the corresponding author on reasonable request.

References

1. Fukuda, Y. et al. Evidence for oscillation of atmospheric neutrinos. Phys. Rev. Lett. 81, 1562-1567 (1998).

   Article  ADS  Google Scholar 

2. Ahmad, Q. R. et al. Direct evidence for neutrino flavor transformation from neutral-current interactions in the Sudbury Neutrino Observatory. Phys. Rev. Lett. 89, 011301 (2002).

   Article  ADS  Google Scholar 

3. King, S. F. Neutrino mass models. Rep. Progr. Phys. 67, 107-157 (2003).

   Article  ADS  Google Scholar 

4. Drexlin, G., Hannen, V., Mertens, S. & Weinheimer, C. Current direct neutrino mass experiments. Adv. High Energy Phys. 2013, 293986 (2013).

   Article  Google Scholar 

5. Formaggio, J. A., de Gouvêa, A. L. C. & Robertson, R. G. H. Direct measurements of neutrino mass. Phys. Rep. 914, 1-54 (2021).

   Article  ADS  Google Scholar 

6. de Gouvêa, A. Neutrino mass models. Annu. Rev. Nucl. Part. Sci. 66, 197-217 (2016).

   Article  ADS  Google Scholar 

7. Group, P. D. et al. Review of particle physics. Progr. Theor. Exp. Phys. 2022, 083C01 (2022).

   Article  Google Scholar 

8. Planck Collaboration Planck 2018 results—VI. Cosmological parameters. Astron. Astrophys. 641, A6 (2020).

   Article  Google Scholar 

9. Aker, M. et al. The design, construction, and commissioning of the KATRIN experiment. J. Instrum. 16, T08015 (2021).

   Article  Google Scholar 

10. Aker, M. et al. Direct neutrino-mass measurement with sub-electronvolt sensitivity. Nat. Phys. 18, 160-166 (2022).

    Article  Google Scholar 

11. Nucciotti, A. The use of low temperature detectors for direct measurements of the mass of the electron neutrino. Adv. High Energy Phys. 2016, 9153024 (2016).

    Article  Google Scholar 

12. Gastaldo, L. et al. The electron capture in 163Ho experiment-ECHo. Euro. Phys. J. Spec. Top. 226, 1623-1694 (2017).

    Article  ADS  Google Scholar 

13. Velte, C. et al. High-resolution and low-background 163Ho spectrum: interpretation of the resonance tails. Euro. Phys. J. C 79, 1026 (2019).

    Article  ADS  Google Scholar 

14. Faverzani, M. et al. The HOLMES experiment. J. Low Temp. Phys. 184, 922-929 (2016).

    Article  ADS  Google Scholar 

15. Nucciotti, A. et al. Status of the HOLMES experiment to directly measure the neutrino mass. J. Low Temp. Phys. 193, 1137-1145 (2018).

    Article  ADS  Google Scholar 

16. Eliseev, S., Novikov, Y. N. & Blaum, K. Penning-trap mass spectrometry and neutrino physics. Annal. Physik 525, 707-719 (2013).

    Article  ADS  Google Scholar 

17. Brown, L. S. & Gabrielse, G. Geonium theory: physics of a single electron or ion in a Penning trap. Rev. Mod. Phys. 58, 233-311 (1986).

    Article  ADS  Google Scholar 

18. Repp, J. et al. PENTATRAP: a novel cryogenic multi-Penning-trap experiment for high-precision mass measurements on highly charged ions. Appl. Phys. B 107, 983-996 (2012).

    Article  ADS  Google Scholar 

19. Filianin, P. et al. Direct Q-value determination of the β− decay of 187Re. Phys. Rev. Lett. 127, 072502 (2021).

    Article  ADS  Google Scholar 

20. Dorrer, H. et al. Production, isolation and characterization of radiochemically pure 163Ho samples for the ECHo-project. Radiochim. Acta 106, 535-547 (2018).

    Article  Google Scholar 

21. Schweiger, C. et al. Production of highly charged ions of rare species by laser-induced desorption inside an electron beam ion trap. Rev. Sci. Instrum. 90, 123201 (2019).

    Article  ADS  Google Scholar 

22. Bradbury, N. E. & Nielsen, R. A. Absolute values of the electron mobility in hydrogen. Phys. Rev. 49, 388-393 (1936).

    Article  ADS  Google Scholar 

23. Schweiger, C. et al. Fast silicon carbide MOSFET based high-voltage push-pull switch for charge state separation of highly charged ions with a Bradbury-Nielsen gate. Rev. Sci. Instrum. 93, 094702 (2022).

    Article  ADS  Google Scholar 

24. Roux, C. et al. The trap design of PENTATRAP. Appl. Phys. B 107, 997-1005 (2012).

    Article  ADS  Google Scholar 

25. Böhm, C. et al. An ultra-stable voltage source for precision Penning-trap experiments. Nucl. Instrum. Methods Phys. Res. Sect. A 828, 125-131 (2016).

    Article  ADS  Google Scholar 

26. Wineland, D. J. & Dehmelt, H. G. Principles of the stored ion calorimeter. J. Appl. Phys. 46, 919-930 (1975).

    Article  ADS  Google Scholar 

27. Feng, X., Charlton, M., Holzscheiter, M., Lewis, R. A. & Yamazaki, Y. Tank circuit model applied to particles in a Penning trap. J. Appl. Phys. 79, 8-13 (1996).

    Article  ADS  Google Scholar 

28. Nagahama, H. et al. Highly sensitive superconducting circuits at ~700 kHhz with tunable quality factors for image-current detection of single trapped antiprotons. Rev. Sci. Instrum. 87, 113305 (2016).

    Article  ADS  Google Scholar 

29. Kromer, K. et al. High-precision mass measurement of doubly magic 208Pb. Euro. Phys. J. A 58, 202 (2022).

    Article  ADS  Google Scholar 

30. Anderson, P. W. Theory of flux creep in hard superconductors. Phys. Rev. Lett. 9, 309-311 (1962).

    Article  ADS  Google Scholar 

31. Anderson, P. W. & Kim, Y. B. Hard superconductivity: theory of the motion of Abrikosov flux lines. Rev. Mod. Phys. 36, 39-43 (1964).

    Article  ADS  Google Scholar 

32. Nagy, S., Blaum, K. & Schuch, R. in Lecture Notes in Physics: Trapped Charged Particles and Fundamental Interactions Vol. 749 (eds Herfurth, F. & Blaum, K.) 1-36 (Springer, 2008).

33. Birge, R. T. The calculation of errors by the method of least squares. Phys. Rev. 40, 207-227 (1932).

    Article  ADS  Google Scholar 

34. Schüssler, R. X. et al. Detection of metastable electronic states by Penning trap mass spectrometry. Nature 581, 42-46 (2020).

    Article  ADS  Google Scholar 

35. Haverkort, M. W., Zwierzycki, M. & Andersen, O. K. Multiplet ligand-field theory using Wannier orbitals. Phys. Rev. B 85, 165113 (2012).

    Article  ADS  Google Scholar 

36. Haverkort, M. W. et al. Quanty (Quanty, 2022); http://www.quanty.org

37. Braβ, M. & Haverkort, M. W. Ab initio calculation of the electron capture spectrum of 163Ho: Auger-Meitner decay into continuum states. New J. Phys. 22, 093018 (2020).

    Article  ADS  Google Scholar 

38. Koepernik, K. & Eschrig, H. Full-potential nonorthogonal local-orbital minimum-basis band-structure scheme. Phys. Rev. B 59, 1743-1757 (1999).

    Article  ADS  Google Scholar 

39. Opahle, I., Koepernik, K. & Eschrig, H. Full-potential band-structure calculation of iron pyrite. Phys. Rev. B 60, 14035-14041 (1999).

    Article  ADS  Google Scholar 

40. Eschrig, H., Richter, M. & Opahle, I. in Relativistic Electronic Structure Theory Theoretical and Computational Chemistry Vol. 14 (ed. Schwerdtfeger, P.) 723-776 (Elsevier, 2004).

41. Grant, I. P. Relativistic calculation of atomic structures. Adv. Phys. 19, 747-811 (1970).

    Article  ADS  Google Scholar 

42. Kotochigova, S., Kirby, K. P. & Tupitsyn, I. I. Ab initio fully relativistic calculations of X-ray spectra of highly charged ions. Phys. Rev. A 76, 052513 (2007).

    Article  ADS  Google Scholar 

43. Fischer, C. F., Gaigalas, G., Jönsson, P. & Bieroń, J. GRASP2018—a Fortran 95 version of the general relativistic atomic structure package. Comput. Phys. Commun. 237, 184-187 (2019).

    Article  ADS  Google Scholar 

44. Fischer, C. F., Godefroid, M., Brage, T., Jönsson, P. & Gaigalas, G. Advanced multiconfiguration methods for complex atoms: I. Energies and wave functions. J. Phys. B 49, 182004 (2016).

    Article  ADS  Google Scholar 

45. Indelicato, P., Lindroth, E. & Desclaux, J. P. Nonrelativistic limit of Dirac-Fock codes: the role of Brillouin configurations. Phys. Rev. Lett. 94, 013002 (2005).

    Article  ADS  Google Scholar 

46. Indelicato, P. Projection operators in multiconfiguration Dirac-Fock calculations: application to the ground state of helium-like ions. Phys. Rev. A 51, 1132-1145 (1995).

    Article  ADS  Google Scholar 

47. Indelicato, P. Nonperturbative evaluation of some QED contributions to the muonic hydrogen n = 2 Lamb shift and hyperfine structure. Phys. Rev. A 87, 022501 (2013).

    Article  ADS  Google Scholar 

48. Indelicato, P., Gorveix, O. & Desclaux, J. P. Multiconfigurational Dirac-Fock studies of two-electron ions. II. Radiative corrections and comparison with experiment. J. Phys. B 20, 651-663 (1987).

    Article  ADS  Google Scholar 

49. Shabaev, V. M., Tupitsyn, I. I. & Yerokhin, V. A. Model operator approach to the Lamb shift calculations in relativistic many-electron atoms. Phys. Rev A. 88, 012513 (2013).

    Article  ADS  Google Scholar 

50. Rodrigues, G. C., Indelicato, P., Santos, J. P., Patté, P. & Parente, F. Systematic calculation of total atomic energies of ground state configurations. Atom. Data Nucl. Data Tables 86, 117-233 (2004).

    Article  ADS  Google Scholar 

51. Savelyev, I. M., Kaygorodov, M. Y., Kozhedub, Y. S., Tupitsyn, I. I. & Shabaev, V. M. Calculations of the binding-energy differences for highly-charged Ho and Dy ions. JEPT Lett. 118, 87-91 (2023).

    ADS  Google Scholar 

52. Wang, M., Huang, W. J., Kondev, F. G., Audi, G. & Naimi, S. The AME 2020 atomic mass evaluation (II). Tables, graphs and references. Chinese Phys. C 45, 030003 (2021).

    Article  ADS  Google Scholar 

53. Tiesinga, E., Mohr, P. J., Newell, D. B. & Taylor, B. N. CODATA recommended values of the fundamental physical constants: 2018. Rev. Mod. Phys. 93, 025010 (2021).

    Article  ADS  Google Scholar 

54. Sturm, S. et al. High-precision measurement of the atomic mass of the electron. Nature 506, 467-470 (2014).

    Article  ADS  Google Scholar 

55. Kramida, A. et al. NIST Atomic Spectra Database v.5.9 (NIST, 2021); https://physics.nist.gov/PhysRefData/ASD/ionEnergy.html

56. Eliseev, S. et al. Direct measurement of the mass difference of 163Ho and 163Dy solves the Q-value puzzle for the neutrino mass determination. Phys. Rev. Lett. 115, 062501 (2015).

    Article  ADS  Google Scholar 

57. Ranitzsch, P. C.-O. et al. Characterization of the 163Ho electron capture spectrum: a step towards the electron neutrino mass determination. Phys. Rev. Lett. 119, 122501 (2017).

    Article  ADS  Google Scholar 

58. Cornell, E. A., Weisskoff, R. M., Boyce, K. R. & Pritchard, D. E. Mode coupling in a Penning trap: π pulses and a classical avoided crossing. Phys. Rev. A 41, 312-315 (1990).

    Article  ADS  Google Scholar 

59. Heiße, F. et al. High-precision mass spectrometer for light ions. Phys. Rev. A 100, 022518 (2019).

    Article  ADS  Google Scholar 

60. Cornell, E. A. et al. Single-ion cyclotron resonance measurement of M(CO+)/M(N2+). Phys. Rev. Lett. 63, 1674-1677 (1989).

    Article  ADS  Google Scholar 

61. Rischka, A. et al. Mass-difference measurements on heavy nuclides with an eV/c2 accuracy in the PENTATRAP spectrometer. Phys. Rev. Lett. 124, 113001 (2020).

    Article  ADS  Google Scholar 

62. Ketter, J. et al. Classical calculation of relativistic frequency-shifts in an ideal Penning trap. Int. J. Mass Spectrom. 361, 34-40 (2014).

    Article  Google Scholar 

63. Ketter, J., Eronen, T., Höcker, M., Streubel, S. & Blaum, K. First-order perturbative calculation of the frequency-shifts caused by static cylindrically-symmetric electric and magnetic imperfections of a Penning trap. Int. J. Mass Spectrom. 358, 1-16 (2014).

    Article  Google Scholar 

64. Schuh, M. et al. Image charge shift in high-precision Penning traps. Phys. Rev. A 100, 023411 (2019).

    Article  ADS  Google Scholar 

Download references

Acknowledgements

We acknowledge funding and support from the Max-Planck-Gesellschaft (C.S., V.D., M.D., P.F., Z.H., J.H., C.H.K., K.K., D.L., Y.N.N., A.R., R.X.S., S.E., K.B.) and the International Max-Planck Research School for precision tests of fundamental symmetries (IMPRS-PTFS) (C.S., M.D., K.K.). This project was also funded by the European Research Council (ERC) under the European Union's Horizon 2020 Research and Innovation Programme under grant agreement numbers 832848-FunI (K.K, K.B.) and 824109-EMP (C.E.), the Deutsche Forschungsgemeinschaft (DFG) under Project-ID 273811115-SFB1225-ISOQUANT (M.B., V.D., M.W.H., C.S.), the Deutsche Forschungsgemeinschaft through grant number INST 40/575-1 FUGG (JUSTUS 2 cluster) (M.B., M.W.H.), the Deutsche Forschungsgemeinschaft Research Unit FOR2202 Neutrino Mass Determination by Electron Capture in 163Ho, ECHo (funding under grant numbers HA 6108/2-1 (M.B., M.W.H.), GA 2219/2-2 (L.G.), EN 299/7-2 (C.E.), EN299/8-2 (C.E.), BL 981/5-1 (K.B.) and DU 1334/1-2 (C.E.D., H.D., D.R.)), the Max-Planck-RIKEN-PTB Center for Time, Constants and Fundamental Symmetries (K.B.) and the state of Baden-Württemberg through bwHPC (M.B., M.W.H.). K.B., P.I. and Y.N.N. are members of the Allianz Program of the Helmholtz Association, contract number EMMI HA-216 'Extremes of Density and Temperature: Cosmic Matter in the Laboratory'. The results presented in this paper are based on work performed before 24 February 2022. This work comprises parts of the PhD thesis work of C.S. to be submitted to Heidelberg University, Germany.

Funding

Open access funding provided by Max Planck Society.

Author information

Author notes

1. Martin Braß

   Present address: Institute of Solid State Physics, TU Wien, Vienna, Austria

2. Rima X. Schüssler

   Present address: Van der Waals-Zeeman Institute, Institute of Physics, University of Amsterdam, Amsterdam, the Netherlands

Authors and Affiliations

1. Max-Planck-Institut für Kernphysik, Heidelberg, Germany

   Christoph Schweiger, Vincent Debierre, Menno Door, Pavel Filianin, Zoltán Harman, Jost Herkenhoff, Christoph H. Keitel, Kathrin Kromer, Daniel Lange, Alexander Rischka, Rima X. Schüssler, Sergey Eliseev & Klaus Blaum

2. Institute for Theoretical Physics, Heidelberg University, Heidelberg, Germany

   Martin Braß & Maurits W. Haverkort

3. Department of Chemistry, Johannes Gutenberg-Universität Mainz, Mainz, Germany

   Holger Dorrer, Christoph E. Düllmann & Dennis Renisch

4. Helmholtz-Institut Mainz, Mainz, Germany

   Christoph E. Düllmann & Dennis Renisch

5. GSI Helmholtzzentrum für Schwerionenforschung GmbH, Darmstadt, Germany

   Christoph E. Düllmann

6. Kirchhoff Institute for Physics, Heidelberg University, Heidelberg, Germany

   Christian Enss & Loredana Gastaldo

7. Laboratoire Kastler Brossel, Sorbonne Université, CNRS, ENS-PSL Research University, Collège de France, Campus Pierre et Marie Curie, Paris, France

   Paul Indelicato

8. Heidelberg University, Heidelberg, Germany

   Daniel Lange

9. NRC 'Kurchatov Institute'-Petersburg Nuclear Physics Institute, Gatchina, Russia

   Yuri N. Novikov

Contributions

C.S., M.D., C.E.D., C.E., P.F., L.G., Y.N.N., A.R., R.X.S., S.E. and K.B. conceived and designed the experiments. C.S., M.D., P.F., J.H., K.K. and S.E. performed the experiments. C.S. and S.E. analysed the data. M.B., V.D., Z.H., M.W.H. and P.I. performed and evaluated the theoretical calculations. M.B., V.D., Z.H., M.W.H., C.H.K. and P.I. contributed to supervising, interpreting and discussing the theoretical calculations. C.S., M.D., H.D., C.E.D., K.K., D.L., D.R., S.E. and K.B. contributed materials/analysis tools. C.S., M.B., V.D., Z.H., M.W.H., P.I. and S.E. wrote the paper. All authors took part in the critical review of the paper.

Corresponding author

Correspondence to Christoph Schweiger.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Physics thanks Pascal Quinet and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Additional information

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Schweiger, C., Braß, M., Debierre, V. et al. Penning-trap measurement of the Q value of electron capture in 163Ho for the determination of the electron neutrino mass. Nat. Phys. (2024). https://doi.org/10.1038/s41567-024-02461-9

Download citation

• Received: 14 June 2023

• Accepted: 01 March 2024

• Published: 19 April 2024

• DOI: https://doi.org/10.1038/s41567-024-02461-9