Measuring the nuclear magnetic quadrupole moment of optically trapped ytterbium atoms in the metastable state

We propose a scheme to measure a nuclear magnetic quadrupole moment (MQM), a CP -violating electromagnetic moment that appears in the nuclear sector, using the long-lived 3P2 metastable state in neutral 173Yb atoms. Laser-cooling and trapping techniques enable us to prepare ultracold 173Yb atoms in the 3P2 state trapped in an optical lattice or an optical tweezer array, providing an ideal experimental platform with long spin coherence time. In addition, our relativistic configuration interaction calculation for the 3P2 electronic wavefunction reveals a large magnetic field gradient generated by the atomic electrons in this state, which amplifies the measurable effect of an MQM. Our scheme could lead to an improvement of more than one order of magnitude in MQM sensitivity compared to the best previous measurement (Murthy et al 1989 Phys. Rev. Lett. 63, 965).


Introduction
While the standard model (SM) of particle physics is a well-established model [1], the violation of the symmetries of charge conjugation and parity (CP) that emerges in the Cabibbo-Kobayashi-Maskawa matrix in the SM is too small to explain the baryon asymmetry of the universe [2,3,4], one of the unsolved questions in modern physics.This motivates the search for sources of larger CP violations that can be described by physics beyond the SM.Under the assumption of CPT invariance (T denoting time reversal), CP violation implies the violation of T symmetry.
The effect of CP violation in elementary particles can appear in low-energy-highresolution measurements using atoms and molecules [5,6,7].An electric dipole moment of an electron (eEDM) probes CP-violation in the leptonic sector, whereas the nuclear Schiff moment and magnetic quadrupole moment (MQM) probe CP-violation in the hadronic sector.Measurements from these two sectors are complementary since they are sensitive to different types of new physics that cause CP-violation [5,7].Despite the excellent precision and accuracy of Schiff moment measurements using Hg [8,9] and Xe [10,11], there remains some ambiguity in relating the measured energy shifts to microscopic CP violation parameters.Nuclear many-body calculations are required to determine the upper limit of CP-violating parameters, such as proton's EDM.However, the discrepancy between the calculations reported is considerable as even the signs of the values may not agree [12], which affects the determination of the upper limit of CPviolation parameters [9,11].Broadly, this difficulty arises because the Schiff moment is a collective nuclear property, which becomes measurable only because of the polarization of the electron cloud in an atom by the finite-sized nucleus [13,14].In contrast, the contributions of CP-violating sources in the MQM can be estimated more accurately because single nucleon P-and T -violating effects directly contributes to the MQM [15,16].We note that the MQM has been investigated not only at the nucleon level , but also at the quark level [17].Importantly, the MQM is enhanced in some heavily deformed nuclei, as has been confirmed using spherical basis calculations [15,18] and Nilsson model calculations [16].While there are several nuclei that can be used for MQM measurements, including some radioactive nuclei, there has been considerable interest in the stable isotope 173 Yb, which has I = 5/2 and large quadrupole deformation that leads to huge enhancement of the MQM [19,20].This fact has led to experimental proposals to measure the MQM using molecules like YbF [20] and YbOH [21,22], exploiting the internal magnetic field gradient induced by the electric field in polarized molecules.Recently, MQM measurements using solids containing deformed nuclei have also been proposed [23,24,25].
In this work, we present a scheme to measure the MQM of 173 Yb using ultracold 173 Yb atoms in the long-lived metastable state 3 P 2 , and accurately compute the MQM sensitivity of this system.Our relativistic configuration interaction theory calculation for the 3 P 2 electronic wavefunction reveals a large magnetic field gradient, a key ingredient for the magnetic quadrupole moment measurement.The preparation of ultracold Yb atoms in the 3 P 2 state in an optical lattice [26] and an optical tweezer array [27] has already been demonstrated, providing an ideal experimental platform with long spin coherence time for the MQM measurement.
Figure 1 shows the overview of the proposed experiment.Figure 2 summarizes the 3 P 2 (4f 14 6s6p) main results of this work.Energy shifts, produced by an MQM of M = 1 µ N fm under an electric field of 50 kV/cm and a magnetic field of 1 G, are shown for all hyperfine levels in the 3 P 2 state.Since magnetic field drifts and noise can affect the measurement of the small shifts due to an MQM, we have also indicated the magnetic moment of each state in the figure .From the information in Fig. 2, we estimate the MQM sensitivity obtained from this scheme based on a typical setup of the optical lattice or optical tweezer array as δM = 5.2 × 10 −8 µ N fm, indicating sufficient precision to improve upon the current experimental bound obtained with 133 Cs atoms [28].We also discuss a way to mitigate the deleterious effects of fluctuating magnetic fields during the experiment, by tuning the experimental parameters such that a chosen pair of hyperfine states has large differential MQM sensitivity but nearly identical magnetic moments.This paper is organized as follows: In Sec. 2 we explain the relevant theory for the atomic energy level shifts due to the nuclear MQM.The computational details for the internal magnetic field gradient are described, along with the details of the calculation of MQM sensitivities for the hyperfine states [29].The proposed experimental scheme for the MQM measurement is presented in Sec. 3. Section 4 is devoted to conclusions.Energy shifts of the hyperfine states in the 3 P 2 state of 173 Yb due to the MQM interaction, ∆E MQM (F, M F ), in units of Hz (above the states).These energy shifts are for an MQM magnitude M = 1 µ N fm and an applied electric field E z = 50 kV/cm.As a measure of the sensitivity of a hyperfine state to magnetic field noise and systematic errors, the magnetic moments, µ F,M F , in units of µ B given by Eq. 24 (below the states in parentheses) are shown for an applied magnetic field B z = 1 G.

Nuclear magnetic quadrupole moment
Similar to the nuclear electric quadrupole interaction [30], the Hamiltonian for the nuclear magnetic quadrupole interaction is given by where i and j are Cartesian coordinate indices, M ij is the i, j component of the magnetic quadrupole moment tensor, and B is the magnetic field.As the magnetic quadrupole moment is extremely small, it is sufficient to consider just the diagonal matrix elements of this Hamiltonian which lead to first-order perturbations.Thus the effective Hamiltonian for a nuclear spin, I, quantized along the z-axis can be expressed (with i = j = z) as Here M is the magnitude of the MQM, M I is the projection of I along the quantization axis z.Inside an atom, the gradient of the magnetic field (∇B) ij is produced by the spins and orbital motion of the atomic electrons.In the four-component relativistic framework, the magnetic field gradient in atomic units is α is the Dirac matrix, and r is the distance of an electron from the center of the nucleus.

Atomic energy level shifts induced by the nuclear MQM
While the nuclear MQM itself has no direct interaction with an electric field, MQM measurements are carried out by measuring energy shifts in the presence of a static electric field, as shown in Fig. 1 and 2. This is because the magnetic field gradient is a P-odd quantity, and a nonzero magnetic field gradient in an atom is only generated after parity mixing between the electronic orbitals in an externally-applied electric field.
In this section, we formally show how the energy shifts due to the MQM arise in the electric field.The discussion is along similar lines to the case of eEDM-induced energy shifts [31,32,33].

The case with no hyperfine coupling
First, for simplicity, we ignore the hyperfine coupling between I and J, and thus the product states |JM J ⟩ |IM I ⟩ are the eigenstates.We express the contribution from ĤMQM in the perturbation theory.The eigenvalue problem of the total Hamiltonian ( Ĥe ) can be written as where E e m and |Ψ e m ⟩ are the energy eigenvalue and eigenstate, respectively.We focus on the 3 P 2 state for |Ψ e m ⟩ state in this study.The total Hamiltonian Ĥe can be separated into the unperturbed ( Ĥe 0 ) and perturbed ( Ĥ′ ) Hamiltonians Here, Ĥe 0 is the atomic electronic Hamiltonian where c is the speed of light, α and β are Dirac matrices, p is the momentum operator, V nuc (r i ) is the electron-nucleus interaction potential, and g(r i , r j ) is the two-electron interaction operator.When V nuc is a spherical potential §, Ĥ0 is P-even, and the zerothorder wavefunction m ⟩ is P-definite.Ĥ′ consists of the contributions from the MQM and the lowest-order electric-dipole interaction of electrons with an externally applied static electric field, Note that the interaction term between the MQM and external fields does not exist.As a result, the first-order contribution for the energy shift E m is purely zero, where m ⟩ is the zeroth-order wavefunction of the unperturbed Hamiltonian Ĥ0 given by, Ĥ0 Ψ (0) This is due to the fact that the MQM term ĤMQM and electric dipole r i are P-odd, while Next, we consider the second-order energy shift.The first-order wavefunction can be expressed by When we retain only the terms that are proportional to M and ignoring terms proportional to M 2 , then the second-order energy shift in an atom is given by, § The case where V nuc is not spherical corresponds to molecules where the atomic states (atomic orbitals) with opposite parity can mix without external electric fields in the molecule-fixed frame.Instead, a (smaller) electric field has to be applied to mix the rotational or Λ(Ω)-doublet states with opposite parity, in order to orient the molecule axis.
Here, we consider E ext = E z ẑ, so that Now, the energy shift can be expressed as E (2) where the atomic EDM ⟨D m ⟩ induced by the MQM can be expressed as follows: The perturbative expression derived above is useful to clarify the physical background.However, our relativistic many-body calculations include the term of the electric dipole interaction with an external electric field in the unperturbed Hamiltonian, instead of the perturbation, Ĥe = Ĥe In this study, its eigenfunction |Ψ e m ⟩ is computed based on the relativistic configuration interaction (CI) theory.The energy shift due to the MQM is Here we replaced m with M J , which is the projection of the total electronic angular momentum J along the quantization axis, to distinguish the ones in the latter section, where the hyperfine coupling is taken into account.With current technology, the maximum external electric field that can be applied in the MQM experiment is about 200 kV/cm ∼ 4 × 10 −4 a.u.[34], which is well within the perturbative regime.Therefore the energy shift is proportional to the electric field E z applied in the experiment to a very good approximation.This fact allows us to define the quantity which only depends on atomic properties and is independent of experimental parameters such as the applied electric field.T M J is thus a measure of the intrinsic sensitivity of an atomic state to the MQM, analogous to the enhancement factor R [35] (or K [36,37]) used in atomic eEDM calculations.Eq. 15 can be expressed in terms of T M J as Here, ĤJz MQM is the effective MQM operator that provides the matrix element described by Eq. 17 in the |JM J ⟩ |IM I ⟩ basis.ĤA 0 includes the hyperfine interactions, the Zeeman interaction of the electronic angular momentum and nuclear spin with a magnetic field as well as the direct current (DC) Stark shift due to an external static electric field, and is given as, Here A hfc and B hfc are the usual hyperfine structure coefficients; g J (g I ) is the g factor of the total electron angular momentum (nuclear spin); µ N and µ B are the nuclear magneton and Bohr magneton, respectively.The atomic parameters in Eq. ( 19) for 173 Yb in the 3 P 2 state are summarized in Table 1.The corresponding eigenstates and eigenvalues are ĤA Since the hyperfine interaction is dominant in the ĤA 0 , F and M F are good quantum numbers, and the energy shift from the effective MQM Hamiltonian can be obtained from the first-order perturbation theory, as follows: The MQM energy shifts for all hyperfine states are shown in Fig. 2. The presented values are obtained by carrying out the calculation described in the following subsections.The key feature is that the sign of the energy shift depends on the signs of M F and E z , reflecting the T -odd and P-odd nature of the MQM Hamiltonian.

Computational details
The electronic structure calculations were carried out using a development version (git hash f49c9b1) of the DIRAC program package [41,42].All calculations were based on the Dirac-Coulomb-Gaunt Hamiltonian including the external electric field (1.0 × 10 −5 a.u ∼ 50 kV/cm) in the form of the electric-dipole approximation to the quantization axis.The ⟨SS|SS⟩-type integrals are explicitly computed.The Dyall basis set with valence triple-ζ quality (dyall v3z) [43] was employed in the uncontracted form, and the small components of the basis sets were generated based on the kinetic balance [44].Diffuse functions are added to s, p, d, and f orbitals in an even-tempered manner.We employed the KR-CI module [45,46,47] for the relativistic configuration-interaction calculations.The Gaussian-type nuclear charge distribution [48] was employed for the electron-nuclear interaction potential.
The atomic spinors, which are the reference state functions for the CI wavefunction, were generated using the average-of-configuration Hartree-Fock (AOC-HF) method [49] at the electronic configuration [Xe]4f 14 6s 1 6p 1 .We employed the generalized-activespace (GAS) method to obtain the CI wavefunctions, where the number of allowed holes can be specified for each shell.We employ three models shown in Tables 2 and 3. Here, two holes are allowed in the 4d, 5s, 5p, 6s, 6p shells, while the truncation of the virtual spinors and the number of holes in the 4p shells are different between each model.To obtain the target electronic configuration [Xe]4f 14 6s 1 6p 1 with certainty, we specified 1 as the maximum occupation number in 6s shell, which provides a hole forcibly in this shell.The remaining electron occupies the 6p shell.The final value for the M J = 1 subspace of 3 P 2 is T M J = −2.43(8)a.u., and the value in the M J = 2 subspace is T M J = −2.85(21)a.u.
Note that Yb has another open-shell state 4f 13 6s 2 5d (J = 2) that has been recently observed [50].This state has a high sensitivity to various phenomena in new physics [51,52,53].According to our preliminary calculations, however, it is not sufficiently sensitive for an MQM search.The parameters T M J =1 = 0.13 a.u. and T M J =2 = 0.16 a.u.for this state, obtained with the relativistic CI method, are much smaller than for the 3 P 2 state.

MQM sensitivities
From the calculated values of T M J , we can evaluate the energy shift due to the MQM.In Fig. 2, we show the energy shift ∆E MQM (F, M F ) (Eq. 21) for a magnetic quadrupole moment magnitude M = 1 µ N fm.Here, the eigenstate Ψ A 0,F,M F is expanded by the basis |IM I ⟩ |JM J ⟩.The magnetic moment µ F,M F [54], a property that parametrizes the sensitivity of |F, M F ⟩ state to external magnetic field fluctuations, is given by Here, g F is the g-factor.The computed values of µ F,M F are shown in Fig. 2.
We find that the maximum MQM sensitivity is not obtained in the state M F = F , except for the F = 9/2 states.A qualitative explanation of this effect is as follows: E e MQM,M J given by Eq. 15 can be separated into nuclear (I and M I ) and electronic (M J ) contributions.From Table 4, the electronic contribution is approximately common for M J = +1, +2 states, and for M J = −1, −2 states.The nuclear part (3M 2 I − I(I + 1)) is −8, −2, 10 for |M I | = 1/2, 3/2, 5/2, respectively.Therefore, the largest MQM shift appears in states where only |M I | = 5/2 contributes, such as |F = 9/2, M F = ±9/2⟩.In the other states, the sensitivity depends on the Clebsch-Gordan (CG) coefficient ⟨IM I JM J |IJF M F ⟩, and therefore the most sensitive nuclear spin sublevels |IM I = ±5/2, JM J ̸ = 0 > does not always have a large contribution to a given |F, M F = ±F ⟩ state.

Experimental scheme
From an experimental perspective, it is desirable to perform spectroscopic measurements between pairs of hyperfine states that have large differential values of ∆E MQM (F, M F ) and small differential µ F,M F .Further, among the many possible combinations of available hyperfine state pairs, transitions between states with the the same |M F | values are useful, because of the insensitivity of their energy differences to the tensor light shift from the optical trapping field, as well as tensor shifts due to the applied static electric field.In particular, a large difference ∆E MQM (F, M F ) and the smallest differential magnetic moment can be obtained using the transition between |F = 3/2, M F = 1/2⟩ → |F = 3/2, M F = −1/2⟩ (see Fig. 2).The reason for the large differential MQM shift between these states is that the CG coefficient ⟨IM I = ±5/2, JM J = ∓2|IJF M F = ±1/2⟩ is largest when F = 3/2 and smallest when F = 9/2.This pair of states is also preferred because they can be connected by a single-photon radio frequency (RF) transition.These states (F = 3/2) can also be conveniently populated from the ground electronic 1 S 0 (F = 5/2) state in Yb using a hyperfineinduced E1 transition [55].
An important source of technical noise and fluctuations in an MQM measurement is the instability of the magnetic field applied to the atoms.In this respect, small sensitivity of the transition frequency to magnetic field fluctuations is advantageous.The energy of a hyperfine state is not simply proportional to the magnetic field, due to the competition between the hyperfine and Zeeman interactions.This fact gives rise to an interesting situation where the resonance frequency of particular transition can be made insensitive to magnetic fields by finding the condition ∂E A 0,F,M F /∂B ext = ∂E A 0,F ′ ,M ′ F /∂B ext .We refer to such a condition as a "magic magnetic field" (MMF).  5.
Similar ideas have been explored in the context of eEDM search with nearly-degenerate states in polyatomic molecules [56] and in the more general context of molecular CPviolation searches with transitions between non-degenerate states, and in the latter case have shown their general existence among many molecular species [21].The possibility of MMFs are investigated in the experimentally convenient region B < 20 G, and the MMF values are summarized in Table 5, where we find several MMFs for E ext values of 50 kV/cm and 100 kV/cm.Two examples of MMFs at E ext = 50 kV/cm are shown in Table 5. Summary of the magic magnetic fields B ext (G) that provide the same gradient of the energy with respect to the magnetic field ∂E A 0,F,M F /∂B ext (MHz/G) for the two states, |F, M F ⟩ and |F ′ , M ′ F ⟩. ∆E (GHz) is the energy difference between the two states.∆E MQM (Hz/µ N fm) is the expected differential MQM shift between the two states assuming M = 1 µ N fm, obtained from Eq. 21.Two different static electric fields E ext = 50, 100 (kV/cm) are considered.δM is the sensitivity to the MQM (×10 −8 µ N fm) for a shot-noise-limited measurement.Bold labels indicate magneticdipole-allowed single-photon transitions between states.Fig. 3, where the magnetic field sensitivity is plotted as a function of magnetic field.
In very small external magnetic fields, the quantization axis could potentially be disturbed by light shifts from the trapping light used in the tweezer array or optical lattice.For our expected experimental parameters, the Stark shift from the static electric field is significantly larger than the light shift, and thus the electric field direction, parallel to the magnetic field, robustly defines the quantization axis.
The proposed MQM measurement can be carried out using ultracold atoms in an optical lattice [57] or an optical tweezer array [58], which has several practical advantages for MQM measurements compared to experiments with beams or vapour cells.One advantage of the optical lattice and optical tweezer systems is the small volume occupied by the trapped-atom ensemble, as a result of which the inhomogeneities of magnetic and electric fields are usually very small.Trapping atoms in ultra-high vacuum conditions easily allows the application of electric fields up to 100 kV/cm.Cold atoms in traps are also ideal platforms to apply advanced quantum metrology methods such as spin squeezing, enabling more precise measurements than the standard quantum limit.Spin squeezing methods for cold atoms, including Yb, have been demonstrated [59,60,61,62,63,64].In addition, a Heisenberg-limited measurement could be possible by using the Greenberger-Horne-Zeilinger state [65], which has been successfully created for an optical tweezer array [66].This yields measurement uncertainty proportional to 1/N , rather than 1/ √ N in the case of the standard quantum limit, where N is the number of atoms.
Our group has demonstrated the preparation of more than 10 4 atoms of 173 Yb in the ground state 1 S 0 at several tens of nano-kelvin, and successful loading of the atoms into a 3D optical lattice, resulting in the formation of a Mott-insulating state where typically one atom is occupied in one site [67,68].The efficient excitation of Yb atoms to the 3 P 2 state in a 3D optical lattice has been demonstrated [26].When we use the transitions between the different |M F | states or different F states for the MQM measurements, as in Table 5, tensor light shifts are important sources of systematic and statistical uncertainties.However, we can work with a magic wavelength for trapping laser light where the tensor light shift vanishes.Our preliminary calculations based on reported values [69,70,71] suggest the existence of such a wavelength around 900 nm for the 3 P 2 state.
Recent experiments from our group have also demonstrated the preparation of Yb atoms in the ground state 1 S 0 in each of 50 sites of an optical tweezer array, followed by excitation to the 3 P 2 state [27].We anticipate the preparation of 10 3 atoms of 173 Yb in the 3 P 2 state is possible.This is justified by the estimation of necessary tweezer light power of about 10 mW per site and the fact that the laser power of more than 10 W at 532 nm is commercially available.Since the radiative lifetime of the 3 P 2 state is 10 s [39], we expect the coherence time between hyperfine states can be made comparable, by employing magnetic field stabilization and dynamical decoupling pulses to stabilize the field at the micro-gauss level.As discussed above, working at a magic magnetic field value can further relax the constraints on magnetic field stability required to achieve 10 s coherence time.
We can use the energy level shift coefficients computed in this paper to estimate the experimental sensitivity to the MQM.For an experiment limited by quantum projection noise from the atomic ensemble ("shot-noise limit"), the statistical uncertainty for the transition frequency is given by δf = 1/ √ N τ T , where N is the number of atoms, τ is the coherence time, and T is the total measurement time.Using conservative experimental parameters, N = 10 3 , τ = 1 s and, T = 2.6 × 10 6 s (thirty days), the frequency precision is δf ≈ 2.0 × 10 −5 Hz.The corresponding sensitivity to the MQM is: The MQM sensitivity values are listed in Table 5.The best sensitivity in this table is δM = 3.7 × 10 −8 µ N fm, under the condition B ext = 16.2G and E ext = 100 kV/cm.In the case of the transition |F = 3/2, M F = 1/2⟩ → |F = 3/2, M F = −1/2⟩ in the condition of Fig. 2, δM = 5.2 × 10 −8 µ N fm is expected.These estimates indicate that an experiment using 173 Yb in the 3 P 2 state has sufficient precision to improve upon the current experimental bound M < 1.6 × 10 −6 µ N fm obtained using Cs [28].

Summary
We have described a competitive method to search for nuclear MQMs using 173 Yb atoms in the 3 P 2 state trapped in an optical lattice or optical tweezer array, and computed relevant atomic properties such as the MQM enhancement, hyperfine magnetic moments and magic magnetic field conditions.Our estimated sensitivity δM ≤ 3.7 × 10 −8 µ N fm has the potential to improve upon the current experimental limit by more than one order of magnitude.This proposal focuses on 173 Yb atoms in the metastable 3 P 2 state, but the methods and ideas described here are generally applicable to new physics searches with other atoms and molecules.

Figure 1 .
Figure 1.Overview of experiment to search for the MQM using an 173 Yb (I = 5/2) atom in the 3 P 2 state.The energy shifts in the hyperfine levels are induced by the MQM as a result of the application of a DC electric field (E ext ) parallel and anti-parallel to the magnetic field (B ext ).The difference in energy shifts between the two cases signals the MQM.This figure only describes the cases of the |F = 1/2, M F = +1/2⟩ and |F = 1/2, M F = −1/2⟩ states, but others can also be employed for the MQM search.

Figure 2 .
Figure 2. Energy shifts of the hyperfine states in the 3 P 2 state of 173 Yb due to the MQM interaction, ∆E MQM (F, M F ), in units of Hz (above the states).These energy shifts are for an MQM magnitude M = 1 µ N fm and an applied electric field E z = 50 kV/cm.As a measure of the sensitivity of a hyperfine state to magnetic field noise and systematic errors, the magnetic moments, µ F,M F , in units of µ B given by Eq. 24 (below the states in parentheses) are shown for an applied magnetic field B z = 1 G.

Table 1 .
(19)ic parameters in Eq(19)that are used in our study.A hfc (MHz) B hfc (MHz) g I (µ N ) g J (µ B )The case of nonzero hyperfine coupling Taking the hyperfine coupling into account, atomic states should be described by the total angular momentum F = Î + Ĵ.The total atomic Hamiltonian to obtain the MQM-energy shift in the |F M F ⟩ states is given by, ĤA = ĤA 0 + ĤJz MQM .

Table 2 .
Generalized active space models for the CI wave functions for the 3 P 2 state.The parameters k and N determine the constraints of the occupation for each GAS shell.The space with N virtual Kramers pairs consists of the canonical AOC-HF spinors.The cutoff energies corresponding to N are provided in Table3.

Table 3 .
Definition of the GASCI wavefunction models.The value of L in the parentheses (L a.u.) is the cutoff energy corresponding to the number of virtual spinors N .e1(e2)meansthateven-tempered diffuse functions 1s1p1d1f (2s2p2d2f ) are added to the dyall v3z basis sets.cvrefers to the core-valence correlation.Electronic structure calculations The results of the computation of T M J are summarized in Table4.Although dyall.v3zincludes the primitive functions for 6p orbitals, the effects of the diffuse functions are significant, especially for M J = 2 state.This is because T M J is a response property to an external electric field, where diffuse functions are required to correctly describe the polarization of the wavefunction.The size of v3z+e1 basis sets is sufficient because the difference of T M J between the v3z+e1 and v3z+e2 is negligible.The M J = 2 state is more sensitive to correlation effects than the M J = 1 state.The final values and uncertainties in Table4are obtained as follows:

Table 4 .
The values of T M J (a.u.) of the Yb atom in 3 P 2 state with various basis sets and different wavefunction models.The main text describes how the final values and uncertainty were estimated.