High-precision mass measurement of doubly magic \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{208}$$\end{document}208Pb

The absolute atomic mass of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{208}$$\end{document}208Pb has been determined with a fractional uncertainty of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$7\times 10^{-11}$$\end{document}7×10-11 by measuring the cyclotron-frequency ratio R of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{208}$$\end{document}208Pb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{41+}$$\end{document}41+ to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{132}$$\end{document}132Xe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{26+}$$\end{document}26+ with the high-precision Penning-trap mass spectrometer Pentatrap and computing the binding energies \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\text {Pb}}$$\end{document}EPb and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\text {Xe}}$$\end{document}EXe of the missing 41 and 26 atomic electrons, respectively, with the ab initio fully relativistic multi-configuration Dirac–Hartree–Fock (MCDHF) method. R has been measured with a relative precision of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$9\times 10^{-12}$$\end{document}9×10-12. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\text {Pb}}$$\end{document}EPb and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\text {Xe}}$$\end{document}EXe have been computed with an uncertainty of 9.1 eV and 2.1 eV, respectively, yielding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$207.976\,650\,571(14)$$\end{document}207.976650571(14) u (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {u}=9.314\,941\,024\,2(28)\times 10^{8}$$\end{document}u=9.3149410242(28)×108 eV/c\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^2$$\end{document}2) for the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{208}$$\end{document}208Pb neutral atomic mass. This result agrees within \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.2\sigma $$\end{document}1.2σ with that from the Atomic-Mass Evaluation (AME) 2020, while improving the precision by almost two orders of magnitude. The new mass value directly improves the mass precision of 14 nuclides in the region of Z = 81–84 and is the most precise mass value with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A>200$$\end{document}A>200. Thus, the measurement establishes a new region of reference mass values which can be used e.g. for precision mass determination of transuranium nuclides, including the superheavies.

The absolute atomic mass of 208 Pb has been determined with a fractional uncertainty of 7 × 10 −11 by measuring the cyclotron-frequency ratio R of 208 Pb 41+ to 132 Xe 26+ with the high-precision Penning-trap mass spectrometer Pentatrap and computing the binding energies E Pb and EXe of the missing 41 and 26 atomic electrons, respectively, with the ab initio fully relativistic multiconfiguration Dirac-Hartree-Fock (MCDHF) method.R has been measured with a relative precision of 9 × 10 −12 .E Pb and EXe have been computed with an uncertainty of 9.1 eV and 2.1 eV, respectively, yielding 207.976 650 571 (14) u (u = 9.314 941 024 2(28) × 10 8 eV/c 2 ) for the 208 Pb neutral atomic mass.This result agrees within 1.2σ with that from the Atomic-Mass Evaluation (AME) 2020, while improving the precision by almost two orders of magnitude.The new mass value directly improves the mass precision of 14 nuclides in the region of Z = 81 − 84 and is the most precise mass value with A > 200.Thus, the measurement establishes a new region of reference mass values which can be used e.g. for precision mass determination of transuranium nuclides, including the superheavies.

I. INTRODUCTION
Heavy and superheavy nuclides beyond the doubly magic nucleus of 208 Pb can only exist due to nuclear shell effects holding them together by counteracting the rapidly increasing Coulomb repulsion with growing proton number Z [1].Insight into these quantum-mechanical nuclear structure effects can be derived from the masses of such nuclides.In addition to some direct heavy mass measurements [2][3][4][5], a network of nuclear transitions and relative mass measurements, i.e. the Atomic-Mass Evaluation (AME), provides mass values for most heavy and superheavy nuclides by tracing them back to a few wellknown masses of uranium isotopes [6].However, no nuclide beyond Z = 70 can be found whose mass is known to a relative precision of better than 2 × 10 −9 to act as a precise reference point for these heavy elements.This directly limits the achievable precision in the heavier mass regions and can possibly lead to tensions or shifts of the relative measured masses due to their referencing to only one reference point.The limitations by mass dependent shifts can be reduced significantly once there is a reference mass with similar mass known to high precision [7].The need for new anchor points for the AME arose during recent mass measurements with TRIGA-TRAP [5,8] at the research reactor TRIGA in Mainz, specifically, an improved absolute mass of 208 Pb [9].Measuring this mass will also directly improve the masses of several Pb isotopes and other nuclides in this mass region [6].
In addition to the impact as a mass reference for other mass measurements, the mass of 208 Pb will soon be * kromer@mpi-hd.mpg.deneeded when the magnetic moment, or the g-factor, of the bound electron of hydrogen-like 208 Pb is planned to be determined by the Penning-trap experiments Alphatrap at the MPIK in Heidelberg [10] and Artemis at GSI Darmstadt [11].This measurement could be the most stringent test of bound-state quantum electrodynamics in strong fields.The error of the mass of the nucleus, however, enters the error budget and therefore needs to be known to high precision [12].With the results of this paper, the error of the mass of 208 Pb will be negligible in future g-factor determinations.
Based on the accurate absolute mass of 132 Xe [13,14], in this paper, we present a determination of the absolute atomic mass of 208 Pb with a fractional uncertainty of 7 × 10 −11 .This is the result of measuring the cyclotron-frequency ratio of 208 Pb 41+ and 132 Xe 26+ with the high-precision Penning-trap mass spectrometer Pentatrap [15,16] in combination with a computation of the binding energies of the missing 41 and 26 atomic electrons, respectively, using the ab initio fully relativistic multi-configuration Dirac-Hartree-Fock (MCDHF) method.The masses of 132 Xe 26+ and 208 Pb 41+ can be related to their neutral counterparts via with m e = 5.485 799 090 65(16) × 10 −4 u being the electron rest mass [17] and m 132 Xe = 131.904155 086 (10) u being the mass of a neutral 132 Xe atom [13,14], each has a relative accuracy of 2.9 × 10

II. EXPERIMENTAL AND THEORETICAL METHODS
If one introduces a charged particle into a magnetic field B, it will describe a free space cyclotron motion with the frequency ω c = q m B, with q/m being the charge-tomass ratio.The working principle of a Penning trap is based on a strong homogeneous magnetic field in combination with an electrostatic quadrupole potential.While the electrostatic potential prevents the ion from escaping in axial direction, forcing it onto an oscillatory axial motion with frequency ω z , the magnetic field forces the ion in radial direction onto a circular orbit with a modified cyclotron frequency ω + .The cross product of the two fields in the Lorentz equation leads to an additional slow drift around the trap center called magnetron motion with frequency ω − .When comparing these three Penning-trap eigenfrequencies to the movement of a free charged particle in a purely magnetic field, it holds [18]: From this equation we can see that the determination of eigenfrequencies of an ion in a Penning trap can be used to determine its mass, if the magnetic field inside the trap is known.However, a determination of a magnetic field of B ≈ 7 T inside a volume of just a few 10 µm 3 to sufficient precision is not possible.Therefore, a relative measurement is chosen at Pentatrap, using a reference ion and a sequential measurement scheme to determine mass ratios [15].Highly charged ions are used due to the advantage that with higher q/m the modified cyclotron frequency increases and can therefore be measured to a higher relative precision.For each mass determination a reference nuclide and charge states have to be chosen that form a q/m doublet with the nuclide of interest in order to largely suppress systematic effects in the cyclotron-frequency ratio determination [15,16].The advantage being, that with q/m doublets the same trapping voltage can be used to match the axial frequency to the detection tank circuit's resonance frequency.Using the same trapping voltage reduces systematic shifts due to trap anharmonicities.In addition, the absolute mass of the reference nuclide has to be known better than the aimed uncertainty of the mass of the nuclide of interest.More technical restrictions are posed by the production of the reference ion, limited by binding energies and the availability of probe material.For these reasons, the near q/m doublet 208 Pb 41+ (q/m = 0.197 138 e/u) and 132 Xe 26+ (q/m = 0.197 113 e/u) [13,14] was chosen.The 132 Xe 26+ ion was created from a gaseous natural source inside a commercial Dresden electron beam ion trap (DREEBIT) [19,20].The DREEBIT is connected to a beamline with a large bender magnet for q/m selection, see Fig. 1a) upper beamline.The 208 Pb 41+ ion was produced in a Heidelberg Compact electron beam ion trap (compact EBIT) [21] equipped with an in-trap laser-desorption target of monoisotopic 208 Pb [22].After ion breeding, the q/m selection was achieved using the time-of-flight separation technique with fast high-voltage switches recently developed at the MPIK [23], supplying the voltages to a Bradbury-Nielson gate [24], see Fig. 1a) lower beamline.Once the ions were selected and decelerated by two pulsed drift tubes, they were consecutively trapped in the first of Pentatrap's five traps and transported down to their individual traps.Due to the five stacked Penning traps available, see Fig. 1b), a simultaneous measurement in two traps is possible, increasing the measurement speed by higher statistics and offering up the opportunity for cross checks between the traps and several analysis methods.Out of the other three traps, two are needed for ion storage and one trap is planned for monitoring, however, currently not in use.
The ion's frequencies depend on the magnetic field and the electrostatic potential.All environmental influences on these quantities need to be stabilized over the duration of the measurement.For this, the Pentatrap laboratory is temperature-stabilized to δT < 50 mK/h and the height of the liquid helium level z lHe used for cooling the superconducting magnet, Penning traps, and the detection system is stabilized to δz lHe < 1 mm/h along with the pressure of helium gas inside the magnet's bore to δp < 10 µbar/h [15].
We employ the Fourier-transform ion-cyclotronresonance detection technique [25] using cryogenic tank circuits connected to the Penning traps to pick up the small image current induced in the trap electrodes by the ion.The largest frequency ω + , and therefore the frequency with the highest contribution to the overall error, is measured using the phase-sensitive pulse and phase (PnP) method [26,27].This method, described in more detail below, sets an initial phase of the reduced cyclotron frequency, then the motion is left decoupled for a variable phase accumulation time t acc during which the phase can evolve freely, before reading out the final phase φ meas .The other two frequencies ω z and ω − are measured with the Fast-Fourier-Transform (FFT) dip and the double-dip technique, respectively [28].
The measurements of 208 Pb 41+ versus 132 Xe 26+ were carried out with the measurement scheme shown in Fig. 2.After a rough estimate of all three frequencies of both ions in both positions, shown in Fig. 1b), using the dip and double-dip technique, the measurement run starts with an N determination, with N being the integer number of full turns of the reduced cyclotron motion during the phase accumulation time.This preparatory measurement is necessary before the actual PnP measurement because the cyclotron phase of the ion increases linearly with time with the increment of the frequency ω + and will thereby pass a full turn of 2π many times during the phase accumulation time.The integer N needs to be known in order to determine the modified cyclotron frequency through the total accumulated phase with φ meas being the measured phase at the end of the accumulation time t acc .The N determination utilizes 9 different phase accumulation times between 0.1 and 40.05 s and finds an ω + for which each N for all accumulation times is integer.
A constant phase offset, unavoidable due to the phase readout, is cancelled out by subtracting a short reference phase with an accumulation time of 0.1 s from each long measurement phase.After concluding the N determination in both traps for lead and xenon, the actual PnP loop is started, see Fig. 2 lower half.Here, a starting phase is imprinted on the modified cyclotron motion by an ω + dipole excitation pulse.The phase can then evolve freely for t acc before the final modified cyclotron phase is imprinted on the axial motion by a coupling π-pulse, where it can be detected as an axial phase [26,27].All excitation and coupling pulses are shaped with a Tukey window [29] in order to avoid systematic phase shifts during the phase readout.During the phase evolution time of the PnP sequence, an FFT axial-dip measurement is performed.This simultaneous phase determination and dip detection leads to a reduction of systematic effects associated with the temporal variation of the trap potential and the magnetic field, since they cancel out when calcu-FIG.2. The figure depicts a flowchart of the measurement scheme and the relevant measurement times during the lead mass campaign.In blue at the top one can see the intial N determination followed by the main measurement PnP loop in orange with the simultaneous measurement of axial frequency and modified cyclotron phase.The ions are frequently swapped in order to minimize the magnetic field drift between position 1 and 2.
lating the free cyclotron frequency, using the invariance theorem in Eq. ( 5).After repeating the measurement of the two ions in trap 2 and 3 ten times, the ions are swapped.If 132 Xe 26+ was in the trap, then 208 Pb 41+ is swapped in and vice versa.This is repeated for around 12 hours before restarting the whole measurement scheme again with the N determination.The magnetron frequency, being a factor ≈ 1, 600 smaller than the modified cyclotron frequency, does not need to be measured repeatedly since the double-dip determination during the preparation phase is sufficiently precise.After one measurement run a relative statistical uncertainty of ≈ 10 −11 is reached.
To determine a neutral mass of 208 Pb from the ions' free cyclotron frequency ratio we need to include the mass of the missing electrons in combination with their binding energies (in the following we always refer to the absolute binding energies).We employ the ab initio fully relativistic multiconfiguration Dirac-Hartree-Fock (MCDHF) method [30][31][32] to compute the binding energies E Xe and E Pb of the outermost 26 and 41 electrons, respectively, in neutral Xe and Pb atoms.First, for the case of Xe, the ionization energy of the outermost 8 electrons has been experimentally determined to be of 424.7 (7) eV [33].Thus, one only needs to calculate the binding-energy difference between the ground states of Xe 26+ ([Ar]3d 10 1 S 0 ) and Xe 8+ ([Kr]4d 10 1 S 0 ).Similarly, since the ionization energy of the outermost 4 electrons in neutral Pb has been measured to be 96.71904(61) eV [33], only the corresponding bindingenergy difference between the ground states of Pb 41+ ([Kr]4d 5 4 P 5/2 ) and Pb 4+ ([Xe]4f 14 5d 10 1 S 0 ) needs to be determined theoretically.In the following, we use E Xe 08−26 and E Pb 04−41 to represent these two terms.
Within the MCDHF scheme, the many-electron atomic state function (ASF) is constructed as a linear combination of configuration state functions (CSFs) with common total angular momentum (J), magnetic (M ), and parity (P ) quantum numbers: |ΓP JM = k c k |γ k P JM .The CSFs |γ k P JM are given as jjcoupled Slater determinants of one-electron orbitals, and γ k summarizes all the information needed to fully define the CSF, i.e. the orbital occupation and coupling of single-electron angular momenta.Γ collectively denotes all the γ k included in the representation of the ASF.The set of CSF basis is generated by the GRASP2018 code [32] via single and double (SD) excitation of electrons from the reference configurations to high-lying virtual orbitals.After solving the self-consistent MCDHF equations for the radial wavefunctions and the mixing coefficients c k , the relativistic configuration interaction (RCI) method is applied to account for the corrections arising from the quantum electrodynamic terms and Breit interactions.We systematically expand the size of the basis set by adding and optimizing virtual orbitals layer by layer up to n = 10 (n is the principle quantum number), with the final correlation energies being derived by extrapolating to n = ∞.
As the ground states of Xe 26+ and Xe 8+ are both in closed-shell configurations, the CSF basis sets used for the calculations can be generated by allowing SD excitations from all the core electrons starting with the 1s orbitals.This gives a contribution from the SD electron correlation energy of 25.6 eV.The contributions from the Breit interactions and QED effects are 4.0 and 0.5 eV, respectively.
For the calculation of Pb 41+ , however, due to its open 4d 5 configuration, the number of CSFs for J = 5/2 easily grows above 4 million even for the SD exchanges of the 4s 2 4p 6 4d 5 electrons and becomes not tractable.Therefore, to access E Pb 04−41 , one needs to construct an ion chain in the calculation to reduce the errors.We first calculate the binding-energy difference between Pb 4+ and Pb 22+ via SD excitations from core electrons down to the 3p subshell, and then calculate the binding-energy difference between Pb 22+ and Pb 36+ by allowing SD excitations of all the core electrons.Finally, the connection between Pb 36+ and Pb 41+ is bridged over Pb 42+ via SD excitations from the 4s orbitals.In total, the SD electron correlation effects contribute 58.1 eV to E Pb 04−41 .The Breit interactions and QED terms give rise to corrections of 9.3 and −0.6 eV, respectively, to the binding energy difference.
Until now, only the SD correlation energies are included in the calculations.Considering that the uncertainties in Breit and QED terms are small, the neglected higher-order correlations will account for the systematic errors.To estimate these errors, we make use of the accurate ionization data of the outermost 8 and 4 electrons in Xe and Pb, respectively.The estimations are based on three observations.First, as a self-consistent theory, the MCDHF always approaches the real groundstate energy from above.Thus, the MCDHF binding energy of a given ionic ground state is always smaller than its real value.Second, for a given element, the contributions from higher-order correlation terms scale with the number of bound electrons.Therefore, the differences between the MCDHF binding energy and its real value is more likely to be smaller in highly charged ions.Lastly, within the same isoelectronic sequence, the contributions from higher-order correlations are always smaller for highly charged ions.This is because, perturbatively, in the denominator of each perturbation term, the energy differences between atomic states in highly charged ions are much larger than those in lower charged ions.As a result, for closed-shell ions, the calculated binding-energy differences based on the SD excitations are always smaller than their real values, but their deviations become narrower when the charge states become larger.
For the case of Xe, we find that the calculated bindingenergy difference between Xe 8+ and Xe is 3.5 eV smaller than the experimentally measured value of 424.7 (7) eV, with the single-electron ionization energies of Xe, Xe 7+ , and Xe 8+ being 0.32 eV, 0.22 eV and 0.10 eV, respectively, smaller than their experimental measurements.Though the deviations of the single-electron ionization energies for some open-shell ions between Xe 8+ and Xe 26+ may be larger than the 0.22 eV deviation of that in Xe 7+ , one can still conservatively expect that the average deviation of the 16 electrons will not be larger than 0.22 eV.This indicates that the systematic shift of E Xe 08−26 shall be within 4.0 eV.To cover the range between 0 to 4.0 eV, one can add a systematic correction of 2.0(2.0)eV, with both systematic shift and uncertainty being 2.0 eV.This leads to E Xe 08−26 = 8, 546.5(2.0)eV and E Xe = 8, 971.2(2.1)eV.
With a similar procedure, the ionization energy of Pb 2+ and Pb 3+ are found to be 1.26 and 0.95 eV smaller than their experimental measured values when SD excitation from the 3p subshells are considered.This indicates an average deviation of < 1.0 eV for the singleelectron ionization energies for Pb 4+ to Pb 22+ , and a systematic correction of 9.0(9.0)eV to the binding energy of the corresponding 18 electrons.For the ions between Pb 22+ and Pb 41+ , since they are close to the isoelectronic systems of Xe ions, one can conservatively assume a < 0.22 eV average deviation of the corresponding single-electron ionization energies.After the summation, we obtain E Pb 04−41 = 28, 633.9(9.1)eV and E Pb = 28, 730.6(9.1)eV.

III. RESULTS
After calculating the free cyclotron frequencies during each PnP loop, the interpolation method [34] is applied to calculate the frequency ratios, see Fig. 3.This method uses two consecutive cyclotron frequencies of one trap in position 1 and interpolates them to the time the cyclotron frequency of the position 2 in the same trap was measured.Then the frequency ratio of the interpolated value of position 1 and the matching value of position 2 can be formed cancelling out in first order the magnetic field drift over time.The linear drift of the magnetic field is ∆B/B = −2.3× 10 −10 /h.The impact of the non-linear drifts of the B-field was thoroughly investigated and found insignificant on the level of the achieved statistical uncertainty.The final measured ion frequency ratio is with a relative statistical uncertainty of 9 × 10 −12 .
The measured ratio is then corrected for known systematic shifts and their respective uncertainties.An overview of the relevant systematic effects and their size is listed in Tab.I.
The largest systematic shift comes from the image charge shift (ICS) [35].The highly charged ions induce an oscillating image charge in the trap electrodes.While this is necessary for detection, it causes a shift of the ions' frequencies by generating a counteracting electric field.The image charge shift depends strongly on the mass difference of the ions and on the radius of the trap, the latter being in the case of Pentatrap 5.000(5) mm.The ICS was determined to be R meas − R = ∆(R meas ) ICS = 1.85(9) × 10 −10 , with R being the corrected ratio.In addition to this, the relativistic shift due to relativistic mass increase [36] leads to another systematic uncertainty related to the size of the excited radii during the PnP measurement scheme: ∆(R meas ) rel = 0(4) × 10 −12 .
All other known systematic effects, due to e.g.trap potential anharmonicity, are on the order of 10 −13 and below and are therefore neglected.Thus, the final ω cratio is R − 1 = 1.252 192 39(9)(10)(13) × 10 −4 , where the number in the first, second, and third brackets indicate the statistical, systematic, and total uncertainty, respectively.
Combining the binding energies of the missing electrons calculated by theory, the experimentally determined mass ratio, and the mass excess of the reference isotope of 132 Xe [13,14] as listed in the AME2020, the mass excess of 208 Pb is determined to be −21, 749.855 (13)

IV. DISCUSSION AND CONCLUSION
In addition to the improvement of precision of the mass of 208 Pb itself, our measurement also improves the masses of a series of lead isotopes, connected by (n, γ) reactions.So far, their mass precision was limited by the precision of 208 Pb, but is now limited by the precision of the energy of the (n, γ) reactions.Furthermore, Fig. 4 shows the improvement in precision of, in total, 14 neighbouring nuclides' masses connected to the mass of 208 Pb via different decays, e.g.α decay, for which the energy is well known.Since the new value of the mass excess of 208 Pb is shifted downward, all these connected nuclides will be shifted down by 1.4 keV.Tab.II lists the new mass values for all nuclides which were significantly improved.With the reported measurement we have established a new region in the nuclear chart with reference masses for experiments on heavy and superheavy nuclides.When measuring masses around m = 200 u the error due to the reference mass will be as low as a few 10 −10 and therefore negligible for mass determinations on radionuclides.
Furthermore, with the new mass precision of 208 Pb of 7×10 −11 the g-factor of the bound electron can be determined to the same level of precision.It therefore allows to carry out the experiments on 208 Pb 81+ at Alphatrap and Artemis without having a large mass dependent error limiting their g-factor determination.

FIG. 1
FIG.1.a) Schematic illustration of the ion production section, the two beamlines, and the combined deceleration region ending in the trap chamber.The ion bunches with their respective energies are shown after they have been q/m selected by the bender magnet or the Bradbury-Nielson gate (BNG).b) Schematic drawing of the Penning-trap tower with two measurement ion configurations.The ions are moved from position 1 to 2 or back every ≈ 15 min.The frequency measurements are carried out in traps 2 and 3. Traps 1 and 4 are utilized as storage traps and trap 5 is currently not in use but is planned to be used to monitor magnetic field fluctuations in the future.

FIG. 3 .
FIG.3.The plot shows the frequency ratios for the different PnP loops.Each color represents a different measurement run.The 1σ error band of the averaged ratio Rmeas is visualized in red.
keV, which amounts to a neutral atomic mass of 207.976 650 571(14) u.The new value improves the mass uncertainty of neutral 208 Pb by a factor of 88 to a relative uncertainty of δm/m = 7 × 10 −11 and shifts the mass excess value by −1.4(1.1)keV.

FIG. 4 .
FIG.4.This colormap depicts a of the nuclide chart.The color corresponds to the improvement factor (with 1 being no improvement) in mass precision after including the new mass value of 208 Pb in the AME[6].The labeled arrows show the relevant connections for the mass determination of the different nuclides via different forms of decays from which the energy is known.

TABLE I .
Systematic shifts and their uncertainties on the measured modified-cyclotron frequency ratio Rmeas.For more details see text.

TABLE II .
[6] mass values of affected nuclides, when including the new mass value of 208 Pb in the AME2020[6].