Detection of the 5p - 4f orbital crossing and its optical clock transition in Pr9.

Recent theoretical works have proposed atomic clocks based on narrow optical transitions in highly charged ions. The most interesting candidates for searches of physics beyond the Standard Model are those which occur at rare orbital crossings where the shell structure of the periodic table is reordered. There are only three such crossings expected to be accessible in highly charged ions, and hitherto none have been observed as both experiment and theory have proven difficult. In this work we observe an orbital crossing in a system chosen to be tractable from both sides: Pr[Formula: see text]. We present electron beam ion trap measurements of its spectra, including the inter-configuration lines that reveal the sought-after crossing. With state-of-the-art calculations we show that the proposed nHz-wide clock line has a very high sensitivity to variation of the fine-structure constant, [Formula: see text], and violation of local Lorentz invariance; and has extremely low sensitivity to external perturbations.

Introduction

Current single ion clocks reach fractional frequency uncertainties around , and enable sensitive tests of relativity and searches for potential variations in fundamental constants[1-5]. These experiments could be improved by exploiting clock transitions with a much reduced sensitivity to frequency shifts caused by external perturbations. Transitions in highly charged ions (HCI) naturally meet this requirement due to their spatially compact wave functions[6,7]. Unfortunately, in most cases this raises the transition frequencies beyond the range of current precision lasers. While some forbidden fine-structure transitions remain in the optical range, these are generally insensitive to physics beyond the Standard Model. More interestingly, at configuration crossings due to re-orderings of electronic orbital binding energies along an isoelectronic sequence, many optical transitions between the nearly degenerate configurations can exist[7,8]. For the and orbitals, this was predicted to occur at Sn-like Pr, see Fig. 1. Here, the P – G magnetic octupole (M3) transition seems ideally suited for an ultra-precise atomic clock and searches for physics beyond the Standard Model with HCI[8-10], recently reviewed in ref. [11]. It is highly sensitive to potential variation of the fine-structure constant, , and to violation of local Lorentz invariance (LLI).

Fig. 1

The orbital crossing and level structures of selected ions along the isoelectronic sequence.

a Configuration-averaged binding energies for relevant configurations of the Sn-like (50 electron) isoelectronic sequence, as a function of atomic number . b–d Grotrian diagrams for low-lying levels of the , , and configurations for , 59, and 62, respectively. Pr is situated close to the configuration crossing point, and the corresponding diagram shows that inter-configuration optical transitions are allowed. Fig. 1c shows lines that were measured (blue) in the electron beam ion trap (EBIT). Strongly mixed levels are indicated with multiple colors. We identified potential clock transitions (dashed magenta lines), which are not observable in the EBIT.

With two electrons above closed shells, Pr has a less complex electronic structure than the open -shell systems studied in previous works[12-14]. Nonetheless, predictions do not reach the accuracy needed for finding the clock transition in a precision laser spectroscopy experiment.

In the following, we present measurements of all the optical magnetic-dipole (M1) transitions with rates of at least order 100 s taking place between the fine-structure states of the and configurations. Since these configurations have both even parity, strongly mixed levels exist, allowing for relatively strong M1 transitions between them. By measuring these and applying the Rydberg-Ritz combination principle, the wavelength of the extremely weak clock transition is inferred to be 452.334(1) nm. Our calculations show that the proposed nHz-wide clock line, proceeds through hyperfine admixture of its upper state with an E2-decaying level. Moreover, we find that it has competitive sensitivities to physics beyond the Standard Model. These favorable properties can be exploited with the detailed experimental scheme that we propose.

Results

Measurement and identification of the lines

The Heidelberg electron beam ion trap (HD-EBIT) was employed to produce and trap Pr ions[15]. In this setup, a magnetically focused electron beam traverses a tenuous beam of CHOPr molecules (CAS number 15492-48-5), which are disassociated by electron impact; further impacts sequentially raise the charge state of the ions until the electron beam energy cannot overcome the binding energy of the outermost electron. The combination of negative space-charge potential of the electron beam and voltages applied to the set of hollow electrodes (called drift tubes) trap the HCI inside the central drift tube. By suitably lowering the longitudinal trapping potential caused by the drift tubes, lower ionization states are preferentially evaporated, so that predominantly Pr ions remain trapped. Several million of these form a cylindrical cloud with a length of ~5 cm and radius of 200 μm. Electron-impact excitation of the HCI steadily populates states which then decay along many different fluorescent channels. Spectra in the range from 220 to 550 nm were recorded using a 2-m focal length Czerny Turner type spectrometer equipped with a cooled CCD camera[16]. Exploratory searches with a broad entrance slit detected weak lines at reduced resolution. By monitoring the line intensities while scanning the electron beam energy we determined their respective charge state, finding 22 Pr lines in total; see Fig. 1c. The charge state identification was made on the basis of a comparison between the estimated electron beam energy at maximum intensity of the lines (135(10) eV), and the predicted ionization energy of Pr (147 eV). Here, lines from neighboring charge states appear approximately an order of magnitude weaker compared to their respective maximal intensity. Tentative line identifications were based on wavelengths and line strengths predicted from ab initio calculations. Our Fock-space coupled-cluster[17] calculations were found to reproduce the spectra with average difference (theory experiment) meV, while our AMBiT[18] calculations (implementing particle-hole configuration interaction with many-body perturbation theory) were accurate to meV (see Methods for further details). Intensity ratios of the lines were compared to predictions by taking into account the wavelength-dependent efficiency of the spectrometer setup and the Pr population distribution in the EBIT. The latter was determined from collisional radiative modeling using the Flexible Atomic Code (FAC)[19].

For confirmation and determination of level energies at the part-per-million level, high-resolution measurements were carried out. This revealed their characteristic Zeeman splittings in the  = 8.000(5) T magnetic field at the trap center. Since the only stable Pr isotope ( = 141) has a rather large nuclear magnetic moment and a nuclear spin of , hyperfine structure (HFS) had to be taken into account to fit the line shapes. We extended the previously employed Zeeman model[12,16], to include HFS in the Paschen-Back regime because for the involved states. We performed a global fit of the complete dataset to ensure consistent factors extracted from lines connecting to the same fine-structure levels. The good agreement of these with AMBiT predictions conclusively confirmed our line identifications, see Fig. 2, Table 1, and Table 2. Level energies with respect to the P ground state were determined from the wavelengths using the LOPT program[20], yielding those necessary to address the F and G clock states, see Table 1. A prominent example of configuration crossing is given by the D and the P levels: their separation of 19 is comparable to Zeeman splitting in the strong magnetic field. This leads to a pronounced quantum interference of the magnetic substates and to an asymmetry of the emission spectra of the above fine-structure levels when decaying to the F state (see Fig. 2), also seen in e.g. the D lines of alkali metals[21]. The non-diagonal Zeeman matrix element between the near-degenerate fine-structure states, characterizing the magnitude of interference, was extracted by fitting the experimental line shapes, and was found to be in good agreement with AMBiT  predictions. The magnetic field strength acts in this setting as a control parameter of quantum interference. In future laser-based high-precision measurements, a much weaker magnetic field will be used. This will minimize interference and enable a better determination of field-free transition energies, and consequently of the clock lines.

Fig. 2

Spectra of four distinctive Pr lines.

The measured spectra (black) were fitted with the hyperfine Paschen-Back model shown in red; blue arrows indicate the positions and relative strengths of Zeeman components (collapsing the hyperfine substructure for a clearer presentation). The orange fit in the lower panel includes the effect of Zeeman mixing of the D and P levels.

Table 1

Overview and comparison of values relevant to the atomic structure of Pr.

Level	Energy (cm^−1)							g_J		A_HFS	q
	Expt.	AMBiT	ΔE	FSCC	ΔE	4-val[10]	ΔE	Expt.	AMBiT	(GHz)	(cm^−1)
5p^2^3P_0	0	0	0	0	0	0	0	0	0	0	0
5p4f^3G_3	22,101.36(5)	21,368	−733	22,248	147	21,895(450)	−206	0.875(2)	0.853	7.771	69,918
5p4f^3F_2	24,494.00(5)	23,845	−649	24,525	31	24,199(370)	−295	0.889(5)	0.883	−1.688	64,699
5p4f^3D_3	27,287.09(5)	26,372	−915	27,575	288	27,002(570)	−285	1.136(4)	1.145	−3.857	74,073
5p^2^3P_1	28,561.063(6)	27,789	−772	28,526	−35	28,436(320)	−125	1.487(3)	1.5	−3.203	39,097
5p4f^3G_4a	29,230.87(6)	28,367	−864	29,482	251	29,343(590)	112	1.130(3)	1.115	5.692	74,358
5p^2^3D_2	36,407.48(6)	35,550	−857	35,980	−427	36,217(380)	−190	1.19(1)	1.139	11.004	51,620
5p4f^3F_3	55,662.43(5)	54,852	−810	55,737	75	55,220(710)	−442	0.940(2)	0.943	2.568	110,266
5p4f^3F_4	59,184.84(5)	58,469	−716	59,393	208			1.158(2)	1.161	2.31	112,108
5p^2^3F_2	62,182.14(2)	61,325	−857	62,380	198			1.028(5)	1.054	2.224	101,716
5p4f^3G_5	63,924.17(6)	62,788	−1136	64,214	290			1.202(2)	1.2	2.347	113,269
5p4f^1F_3	63,963.57(6)	62,721	−1243	64,379	415			1.197(6)	1.226	0.485	112,004
5p^2^3P_2	67,290.97(5)	66,350	−941	67,343	52			1.210(4)	1.207	2.331	98,759
5p4f^3D_1	67,309.3(1)	66,429	−880	67,925	616			0.54(1)	0.5	−2.432	110,679
5p4f^3G_4b	69,861.70(8)	68,528	−1334	70,193	331			1.039(4)	1.023	2.62	111,833
4f^2^3F_2		79,693		81,801					0.813	1.026	118,508
5p4f^1D_2		80,569		82,657					0.907	1.555	123,702

Measured values are shown in the Expt. columns. Calculated energies and Landé factors of the Pr states stem from AMBiT, Fock-space coupled-cluster (FSCC), and CI+all-order calculations with four valence electrons (4-val). Magnetic-dipole hyperfine structure constants and sensitivities to -variation were calculated using the AMBiT code

Table 2

Energies and rates of Pr lines. The measured values (Expt.) include a 1 standard deviation uncertainty given in brackets, they are compared to ab initio theory predictions.

Lower	Upper	Energy (eV)			Rate (s^−1)
		Expt.	AMBiT	FSCC
5p^2^3D_2	5p^2^3F_2	3.19563(1)	3.1957	3.2732
5p4f^3D_3	5p4f^3F_3	3.51807(2)	3.5311	3.4916
5p^2^3P_0	5p^2^3P_1	3.5411202(8)	3.4454	3.5368	1021
5p4f^3G_4a	5p4f^3F_4	3.713810(3)	3.7322	3.7085	1086
5p^2^3D_2	5p^2^3P_2	3.829068(5)	3.8187	3.8885	374
5p4f^3F_2	5p4f^3F_3	3.864391(3)	3.8444	3.8698	823
5p4f^3D_3	5p4f^3F_4	3.954826(3)	3.9795	3.9449	675
5p4f^3G_3	5p4f^3F_3	4.161043(2)	4.1515	4.1521	1342
5p^2^3P_1	5p^2^3F_2	4.168481(3)	4.1579	4.1974	493
5p4f^3G_4a	5p4f^3G_5	4.301421(1)	4.2677	4.3062	4334
5p4f^3G_4a	5p4f^1F_3	4.306314(3)	4.2594	4.3267	640
5p4f^3D_3	5p4f^1F_3	4.547299(3)	4.5067	4.5631	1991
5p4f^3G_3	5p4f^3F_4	4.597753(6)	4.5999	4.6054	1660
5p4f^3F_2	5p^2^3F_2	4.672740(8)	4.6469	4.6934	483
5p^2^3P_1	5p^2^3P_2	4.80191(1)	4.7810	4.8127	1060
5p4f^3D_3	5p^2^3P_2	4.959843(4)	4.9566	4.9306	762
5p4f^3G_3	5p^2^3F_2	4.969391(7)	4.9540	4.9757	386
5p4f^3G_4a	5p4f^3G_4b	5.037586(9)	4.9793	5.0475	876
5p4f^3G_3	5p4f^1F_3	5.19021(2)	5.1271	5.2236	252
5p4f^3D_3	5p4f^3G_4b	5.27855(2)	5.2267	5.2840	878
5p4f^3F_2	5p^2^3P_2	5.30616(2)	5.2699	5.3088	1022
5p4f^3F_2	5p4f^3D_1	5.30842(2)	5.2797	5.3809	1199

The rate is calculated using theoretical matrix elements from AMBiT and experimental transition frequencies

Properties of the clock transition

After discovering this orbital crossing and the interesting clock transition, in the following we discuss its properties. Without HFS, the G state would decay through a hugely suppressed M3 transition with a lifetime of order 10 million years – a 3 fHz linewidth. However, admixture with the F state by hyperfine coupling induces much faster E2 transitions (lifetime years) with widths on the order of nHz, see Fig. 3. In state-of-the-art optical clocks such transitions have been probed[4] and become accessible in HCI using well-established quantum-logic spectroscopy techniques[22-24]. The state is decoupled from the ground state but decays to the component with a lifetime of 400 years. For comparison, the F E2 transition to the ground state has a much broader linewidth of 6.4 mHz, similar to that of the Al clock[5,25].

Fig. 3

Schematic level diagram.

Shown are the ground and lowest two fine-structure states in Pr between which the clock transitions take place, energy splittings not to scale. Levels with the same total angular momentum are admixed by the magnetic hyperfine interaction, allowing the G clock states to decay via E2 transition rather than M3, which would take of order 10 million years to decay.

Blackbody-radiation (BBR) shift is a dominant source of systematic uncertainty in some atomic clocks such as Yb[4]. However, the static dipole polarisability (to which the BBR shift is proportional) is strongly suppressed in HCI by both the reduced size of the valence electron wave functions and the typically large separations between mixed levels of opposite parity. Together, these lead to a scaling where is the effective screened charge that the valence electron experiences[8]. Our calculations of the static polarisability for Pr yield for the clock state of a.u. and confirm the expected suppression. Furthermore, the ground state polarisability is rather similar, so the differential polarisability for the clock transition is a.u., ten times smaller than that of the excellent Al clock transition with a.u.[5]. An atomic clock based on Pr would therefore be extremely resilient to BBR even at room temperature.

Beyond their favourable metrological properties, HCI have been suggested for probing variation of fundamental constants due to their high sensitivity to this[6,7]. Sensitivity to variation of a transition with frequency is usually characterized by the parameter , defined by the equationwhere is the frequency at the present-day value of the fine-structure constant and . Calculated values for Pr levels are presented in Table 1, and compared to other proposed transitions in Table 3. The Pr M3 clock transition has a sensitivity similar to that of the 467 nm E3 clock transition in Yb, SF ( cm[26]), but with opposite sign. The sign change can be understood in the single-particle model: the Yb transition is , while in Pr it is , leading to opposite sign. Comparison of these two clocks would therefore lead to improved limits on  variation and allow control of systematics.

Table 3

Predicted sensitivities to physics beyond the Standard Model.

Ion	Ref.	Level	K	⟨J∣∣T^(2)∣∣J⟩
Pr^9+		5p4f^3G_3	6.32	74.2
		5p4f^3F_2	5.28	57.8
Ca^+	[29]	3d^2D_3∕2		7.09 (12)
		3d^2D_5∕2		9.25 (15)
Yb^+	[64,65]	4f^145d^2D_3∕2	1.00	9.96
		4f^145d^2D_5∕2	1.03	12.08
		4f^136s^2^2F_7∕2	−5.95	−135.2
Dy	[28,65]	4f^105d6sJ=10	0.77	69.48
		4f^95d^26sJ=10	2.55	49.73
Hg^+	[1,65]	5d^96s^2^2D_5∕2	−2.94

The relative sensitivity to variation of the fine-structure constant is given with respect to the ground state, but note that in dysprosium the fractional sensitivity of transitions between the two upper levels is many orders of magnitude higher since these levels are almost degenerate. Sensitivity to local Lorentz invariance is represented with the reduced matrix elements of the operator, given in atomic units

We also investigate the sensitivity of the clock transitions to invariance under local Lorentz transformations. LLI is a fundamental feature of the Standard Model and has been tested in all sectors of physics[27]. While Michelson-Morley experiments verify the isotropy of the speed of light, recent atomic experiments have placed strong limits on LLI-breaking parameters in the electron-photon sector[28,29]. The sensitivity of transitions for such studies is given by the reduced matrix element of , defined bywhere is the speed of light, (, ) are Dirac matrices, and is the momentum of a bound electron. We find a.u. for the G state, similar in magnitude to the most sensitive Dy and Yb clock transitions, see Table 3. Again, the sign is opposite to Yb E3, making their comparison more powerful and improving the control of potential systematic effects. Furthermore, the value compares well with other HCI[30].

Proposal for precision laser spectroscopy

Future precision spectroscopy of the clock transitions will require that the internal Pr state be prepared and detected using quantum-logic protocols in HCI that are sympathetically cooled in a cryogenic Paul trap[31,32]. Populations calculated using FAC show that a Pr ion ejected from the EBIT ends up after a few minutes in either the P (25%) or G state (75%). We propose to employ state-dependent oscillating optical dipole forces (ODF) formed by two counter-propagating laser beams detuned by one of the trapping frequencies of the two-ion crystal with respect to each other. Electronic- and hyperfine-state selectivity may be achieved by tuning the ODF near one of the HCI resonances. If the HCI is in the target state, the ODF exerts an oscillating force onto the HCI. This displaces its motional state, which can be detected efficiently on the co-trapped Be ion[23,33], further enhanced by employing non-classical states of motion[34]. We estimate an achievable displacement rate of tens of kHz exciting coherent states of motion for realistic parameters of laser radiation at 408 nm detuned by 1 MHz from the P – F transition. Similarly, motional excitation rates of a few kHz can be achieved by detuning laser radiation at 452 nm by 10 Hz from the nHz-wide P – G clock transition. Since the ion is long-lived in both states, we distinguish between the two clock states by detuning the beams forming the ODF to additionally change the magnetic substate[24]. State selectivity can be provided by the global detuning of the ODF and the unique -factor of the electronic and hyperfine states. In case the HCI is in one of the excited hyperfine states of G, the lowest hyperfine state can be prepared deterministically by driving appropriate microwave -pulses between states, followed by detection if the target state has been reached. The substate is then prepared by driving states using appropriately detuned Raman laser beams on the red sideband, followed by ground state cooling[22,24].

Despite the high resolution of the measured lines, the P – G clock transition with expected nHz linewidth is only known to within 1.5 GHz. We can improve this by measuring the P – F and G – F transitions and performing a Ritz combination. The transitions can be found by employing again counter-propagating Raman beams that excite motion inversely proportional to the detuning of the Raman resonance to the electronic transition[23]. Using realistic laser parameters similar to above, tuned near the electronic resonance, we estimate that a scan with 10 MHz steps using 10 ms probe pulses can be performed without missing the transition. Assuming that each scan point requires 100 repetitions we can scan the range in 10 min. By reducing the Raman laser power while extending the probe time, the resolution can be enhanced to the sub-MHz level. Once identified, quantum-logic spectroscopy following[1,22] on the clock transition can commence. Since the linewidth of the transition is significantly narrower compared to the currently best available lasers[35], we expect laser-induced ac Stark shifts that need mitigation using hyper-Ramsey spectroscopy or a variant thereof[36-45].

Discussion

We have measured optical inter-configuration lines of Pr, finding the – orbital crossing, and thereby determined the frequency of the proposed P – G clock transition with an accuracy sufficient for quantum-logic spectroscopy at ultra-high resolution. Our state-of-the art calculations agree well with the measurements, thus we used the obtained wave functions to predict the polarizabilities of the levels and their sensitivities physics beyond the Standard Model. These are crucial steps towards future precision laser spectroscopy of the clock transition, for which we have also proposed a detailed experimental scheme.

Methods

Line-shape model

In the hyperfine Paschen-Back regime, the energy shift of a fine-structure state’s magnetic sublevel in an external magnetic field is given byHence, a transition between two fine-structure states has multiple components with energies with the transition energy without an external field, and . Here, and in the following, primed symbols differentiate upper states from lower states. Taking into account the Gaussian shape of individual components, the line-shape function is defined aswhere the sum is taken overall combinations of upper and lower magnetic sublevels. The factors take into account the known efficiencies of the setup for the two perpendicular linear polarizations of the light, and is the common linewidth, which is determined by the apparatus response and Doppler broadening. The magnetic-dipole matrix elements are given byThe large parentheses denotes a Wigner -symbol. It follows that . In case of the asymmetric lines, the above initial fine-structure states were transformed to the eigenstates of the Zeeman Hamiltonian also taking into account non-diagonal coupling.

Equation (4) was fitted to each measured line with the following free parameters: , , a constant baseline value (noise floor), and overall amplitude. The results for the extracted transition energies are presented in Table 2. The values of all involved states were kept as global fit parameters for the whole dataset of 22 lines. Shot noise and read-out noise were taken into account for the weighting of the data points. The final uncertainties on the transition energies for each line were taken as the square root of its fit uncertainty and calibration uncertainty added in quadrature.

Fock-space coupled-cluster calculations

The Fock-space coupled-cluster (FSCC) calculations of the transition energies were performed within the framework of the projected Dirac-Coulomb-Breit Hamiltonian[46]. In atomic units (),Here, is the one-electron Dirac Hamiltonian, and are the four-dimensional Dirac matrices and . The nuclear potential takes into account the finite size of the nucleus, modeled by a uniformly charged sphere[47]. The two-electron term includes the nonrelativistic electron repulsion and the frequency-independent Breit operator,and is correct to second order in the fine-structure constant .

The calculations started from the closed-shell reference [Kr]4 configuration of Pr. In the first stage the relativistic Hartree-Fock equations were solved for this closed-shell reference state, which was subsequently correlated by solving the coupled-cluster equations. We then proceeded to add two electrons, one at at time, recorrelating at each stage, to reach the desired valence state of Pr. We were primarily interested in the and the configurations of Pr; however, to achieve optimal accuracy we used a large model space, comprised of 4 , 5 , 4 , 5 , 3 , 2 , and 1 orbitals. The intermediate Hamiltonian method[48] was employed to facilitate convergence.

The uncontracted universal basis set[49] was used, composed of even-tempered Gaussian type orbitals, with exponents given by

The basis set consisted of 37 s, 31 p, 26 d, 21 f, 16 g, 11 h, and 6 i functions; the convergence of the obtained transition energies with respect to the size of the basis set was verified. All the electrons were correlated.

The energy calculations were performed using the Tel-Aviv Relativistic Atomic Fock-Space coupled-cluster code (TRAFS-3C), written by E. Eliav, U. Kaldor and Y. Ishikawa. The final FSCC transition energies were also corrected for the QED contribution, calculated using the AMBiT program.

Polarizabilities

The polarizabilities were also calculated using the Fock-space coupled-cluster method within the finite-field approach[50,51]. We used the DIRAC17 program package[52], as the Tel-Aviv program does not allow for addition of external fields. The v3z basis set of Dyall was used[53]; 20 electrons were correlated and the model space consisted of 5 and 4 orbitals. These calculations were carried out in the framework of the Dirac-Coulomb Hamiltonian, as the Breit term is not yet implemented in the DIRAC program.

CI+MBPT calculations

Further calculations of energies and transition properties were obtained using the atomic code  AMBiT[18]. This code implements the particle-hole CI+MBPT formalism[54], which builds on the combination of configuration interaction and many-body perturbation theory described in ref. [55] (see also ref. [56]). This method also seeks to solve Eqs. (5)–(7), but treats electron correlations in a very different way. Full details may be found in ref. [18]; below we present salient points for the case of Pr.

For the current calculations, we start from relativistic Hartree-Fock using the same closed-shell reference configuration used in the FSCC calculations: [Kr]. We then create a B-spline basis set in this potential, including virtual orbitals up to and . In the particle-hole CI+MBPT formalism the orbitals are divided into filled shells belonging to a frozen core, valence shells both below and above the Fermi level, and virtual orbitals.

The CI space includes single and double excitations from the , , and (“leading configurations”) up to , including allowing for particle-hole excitations from the and shells. This gives an extremely large number of configuration state functions (CSFs) for each symmetry, for example the matrix has size . To make this problem tractable we use the emu CI method[57] where the interactions between highly-excited configurations with holes are ignored.

Correlations with the frozen core orbitals (including , , shells and those below) as well as the remaining virtual orbitals ( or ) are treated using second-order MBPT by modifying the one and two-body radial matrix elements[55]. The effective three-body operator is applied to each matrix element separately; to reduce computational cost it is included only when at least one of the configurations involved is a leading configuration[54]. Finally, for the energy calculations we include an extrapolation to higher in the MBPT basis[10] and Lamb shift (QED) corrections[58-60].

Diagonalisation of the CI matrix gives energies and many-body wave functions for the low-lying levels in Pr. Using these wave functions we have calculated electromagnetic transition matrix elements (and hence transition rates), hyperfine structure, and matrix elements of (Eq. (2)). For all of these we have included the effects of core polarisation using the relativistic random-phase approximation (see ref. [61] for relevant formulas). By contrast values (Eq. (1)) were obtained in the finite-field approach by directly varying in the code and repeating the energy calculation. In Table 3, the predicted matrix elements and sensitivities are compared to those of other proposed systems and those already under investigation in relevant high-precision experiments.

Lifetime calculations

Direct decay of the G clock state to the ground state proceeds as an M3 transition, which is hugely suppressed and would indicate a lifetime of order 10 million years. However, the hyperfine components of the Pr clock state have a small admixture of levels, allowing for decay via a much faster E2 transition. The rate of the hyperfine-interaction-induced decay can be expressed as a generalized E2 transitionwhere is the quantum number of total angular momentum of the upper state () and the amplitude can be expressed asHere and are the operators of the hyperfine dipole interaction and the electric quadrupole amplitude, respectively. For the clock transition the sum over intermediate states in Eq. 10 is dominated by the lowest states with : F and P.

MCDF calculations

Transition energies, hyperfine structure coefficients, and factors were also evaluated in the framework of the MCDF and relativistic CI methods, as implemented in the GRASP2K atomic structure package[62], and were found to be in reasonable agreement with the experiment on the level of the CI+MBPT results. As a first step, an MCDF calculation was performed, with the active space of one- and two-electron exchanges ranging from the , , and spectroscopic orbitals up to . In the second step, the active space was extended to also include the orbitals for a better account of core polarisation effects, and CI calculations were performed with the optimized orbitals obtained in the first step. The extension of the active space in the second step has lead to 946k -coupled configurations. For a more detailed modeling of the spectral line shapes, the non-diagonal matrix elements of the hyperfine and Zeeman interactions[63] and mixing coefficients for sublevels of equal magnetic quantum numbers were also evaluated.