Theoretical study on the hyperfine interaction constants and the isotope shift factors for the $3s^2~^1 S_{ 0} ~-~ 3s3p~^{3,1}P^o_1$ transitions in Al$^+$
Tingxian Zhang, Luyou Xie, Jiguang Li, and Zehuang Lu

TL;DR
This theoretical study calculates hyperfine interaction constants and isotope shift factors for specific transitions in Al$^+$ ions, providing precise atomic parameters crucial for nuclear property analysis.
Contribution
It offers detailed multiconfiguration Dirac-Hartree-Fock calculations of hyperfine and isotope shift parameters, including electron correlation and relativistic effects, with uncertainties below 1.5%.
Findings
Hyperfine constants have uncertainties less than 1.5%.
Electron correlations significantly affect isotope shift factors.
Relativistic nuclear recoil corrections are negligible for Al$^+$.
Abstract
We calculated the magnetic dipole and the electric quadrupole hyperfine interaction constants of 3s3p states and the isotope shift, including mass and field shift, factors for transitions from these two states to the ground state 3s in Al ions using the multiconfiguration Dirac-Hartree-Fock method. The effects of the electron correlations and the Breit interaction on these physical quantities were investigated in detail based on the active space approach. It is found that the CC and the higher-order correlations are considerable for evaluating the uncertainties of the atomic parameters concerned. The uncertainties of the hyperfine interaction constants in this work are less than 1.5\%. Although the isotope shift factors are highly sensitive to the electron correlations, reasonable uncertainties were obtained by exploring the effects of the electron…
| Reference configurations | AO() | Model | NCSF | |||
|---|---|---|---|---|---|---|
| {2s22p63s2} | DF | 1 | ||||
| {4s, 3p, 3d, 4f} | CV1s- | 61 | ||||
| {5s, 4p, 4d, 5f, 5g} | CV1s- | 254 | ||||
| {6s, 5p, 5d, 6f, 6g, 6h} | CV1s- | 616 | ||||
| {7s, 6p, 6d, 7f, 7g, 6h} | CV1s- | 1098 | ||||
| {8s, 7p, 7d, 8f, 7g, 6h} | CV1s- | 1603 | ||||
| {9s, 8p, 8d, 9f, 7g, 6h} | CV1s- | 2211 | ||||
| {10s, 9p, 9d, 10f, 7g, 6h} | CV1s- | 2922 | ||||
| {11s, 10p, 10d, 10f, 7g, 6h} | CV1s- | 3509 | ||||
| {12s, 11p, 11d, 10f, 7g, 6h} | CV1s- | 4163 | ||||
| {13s, 12p, 12d, 10f, 7g, 6h} | CV1s- | 4884 | ||||
| {13s, 12p, 12d, 10f, 7g, 6h} | CC2s | 14988 | ||||
| {2s22p63p2} | {8s, 8p, 7d, 8f, 7g, 6h} | MR1- | 23358 | |||
| {2s22p43s24p2; 2s22p43s23d2; | {6s, 7p, 6d, 6f, 6g, 6h} | MR2- | 229659 | |||
| 2s22p53s3p3d } | ||||||
| {2s22p43s23p4p; 2s22p43s23d4d; | {5s, 6p, 5d, 4f, 5g} | MR3- | 357223 | |||
| 2s2p53s24s4p; 2s22p43s24p5p} |
| Reference configurations | AO() | Model | NCSF | |
|---|---|---|---|---|
| {2s22p63s3p} | DF | 2 | ||
| {4s, 4p, 3d, 4f} | CV1s- | 600 | ||
| {5s, 5p, 4d, 5f, 5g} | CV1s- | 2326 | ||
| {6s, 6p, 5d, 6f, 6g, 6h} | CV1s- | 5573 | ||
| {7s, 7p, 6d, 7f, 7g, 6h} | CV1s- | 9860 | ||
| {8s, 8p, 7d, 8f, 7g, 6h} | CV1s- | 14292 | ||
| {9s, 9p, 8d, 9f, 7g, 6h} | CV1s- | 19594 | ||
| {10s, 10p, 9d, 10f, 7g, 6h} | CV1s- | 25766 | ||
| {11s, 11p, 10d, 10f, 7g, 6h} | CV1s- | 30696 | ||
| {12s, 12p, 11d, 10f, 7g, 6h} | CV1s- | 36146 | ||
| {13s, 13p, 12d, 10f, 7g, 6h} | CV1s- | 42116 | ||
| {13s, 13p, 12d, 10f, 7g, 6h} | CC2s | 223468 | ||
| {2s22p63p3d} | {8s, 8p, 8d, 8f, 7g, 6h} | MR1- | 352702 | |
| {2s22p63s5p; 2s22p63s3p; 2s22p63d5p} | {5s, 7p, 5d, 5f, 5g} | MR2- | 400690 | |
| {2s22p43s3p4p2; 2s22p43s3p4d2} | {4s, 6p, 5d, 4f } | MR3- | 1327012 |
| Model | 3s2 | 3s3p | 3s3p | |||||
| Configurations | Configurations | Configurations | ||||||
| MR1 | 2s22p63s2 | 0.9810 | 2s22p63s3p | 0.9896 | 2s22p63s3p | 0.9781 | ||
| 2s22p63p2 | 0.1382 | 2s22p63p3d | 0.0557 | 2s22p63p3d | 0.1439 | |||
| MR2 | 2s22p43s23d2 | 0.0385 | 2s22p63s5p | 0.0175 | 2s22p63s5p | 0.0360 | ||
| 2s22p43s24p2 | 0.0300 | 2s22p63s4p | 0.0209 | 2s22p63s4p | 0.0270 | |||
| 2s22p53s3p3d | 0.0204 | 2s22p63d5p | 0.0085 | 2s22p63d5p | 0.0262 | |||
| MR3 | 2s22p43s23p4p | 0.0179 | 2s22p43s3p4d2 | 0.0194 | 2s22p43s3p4d2 | 0.0188 | ||
| 2s22p43s23d4d | 0.0162 | 2s22p43s3p4p2 | 0.0160 | 2s22p43s3p4p2 | 0.0111 | |||
| 2s2p53s24s4p | 0.0162 | |||||||
| 2s22p43s24p5p | 0.0158 | |||||||
| Model | 3s3p | 3s3p | |||
|---|---|---|---|---|---|
| A | B | A | B | ||
| MR3- | 1327.3(10.2) | 15.1(0.2) | 197.1(2.8) | 22.1(0.3) | |
| Model | |||||||
|---|---|---|---|---|---|---|---|
| DF | 614 | 844 | 0.1944 | 615 | 129 | 0.1945 | |
| CV1s- | 788 | 1646 | 0.1915 | 911 | 732 | 0.1921 | |
| CV1s- | 621 | 967 | 0.1965 | 887 | 56 | 0.1873 | |
| CV1s- | 647 | 943 | 0.2025 | 957 | 86 | 0.1928 | |
| CV1s- | 622 | 972 | 0.2009 | 937 | 46 | 0.1911 | |
| CV1s- | 633 | 975 | 0.2029 | 975 | 49 | 0.1928 | |
| CV1s- | 632 | 975 | 0.2028 | 978 | 48 | 0.1930 | |
| CV1s- | 623 | 974 | 0.2026 | 971 | 49 | 0.1927 | |
| CV1s- | 624 | 974 | 0.2026 | 977 | 50 | 0.1927 | |
| CV1s- | 623 | 974 | 0.2026 | 977 | 51 | 0.1927 | |
| CV1s- | 624 | 974 | 0.2026 | 977 | 51 | 0.1927 | |
| CC2s | 597 | 596 | 0.2021 | 823 | 152 | 0.1930 | |
| MR1- | 617 | 771 | 0.1928 | 841 | 112 | 0.1890 | |
| MR2- | 718 | 1154 | 0.1982 | 1064 | 452 | 0.1916 | |
| MR3- | 608 | 989 | 0.1971 | 931 | 268 | 0.1904 | |
| Transition | |||||
|---|---|---|---|---|---|
| 632(88) | 608(90) | 1003(163) | 989(171) | 0.1971(0.0059) | |
| 958(38) | 931(37) | 279(139) | 268(134) | 0.1904(0.0038) |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Theoretical study on the hyperfine interaction constants and the isotope shift factors for the 3s\bm{{}^{2}}\bm{^{1}S_{0}$$-$${3s3p}}\bm{}transitionsinAl$^+
Tingxian Zhang
College of Physics and Electronic Engineering, Northwest Normal University, Lanzhou 730070 P. R. China
Institute of Applied Physics and Computational Mathematic, Beijing 100088, P. R. China
Luyou Xie
College of Physics and Electronic Engineering, Northwest Normal University, Lanzhou 730070 P. R. China
Jiguang Li
Institute of Applied Physics and Computational Mathematic, Beijing 100088, P. R. China
Zehuang Lu
MOE Key Laboratory of Fundamental Physical Quantities Measurement, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, P. R. China
Abstract
We calculated the magnetic dipole and the electric quadrupole hyperfine interaction constants of 3s3p states and the isotope shift, including mass and field shift, factors for transitions from these two states to the ground state 3s in Al+ ions using the multiconfiguration Dirac-Hartree-Fock method. The effects of the electron correlations and the Breit interaction on these physical quantities were investigated in detail based on the active space approach. It is found that the CC and the higher-order correlations are considerable for evaluating the uncertainties of the atomic parameters concerned. The uncertainties of the hyperfine interaction constants in this work are less than 1.5%. Although the isotope shift factors are highly sensitive to the electron correlations, reasonable uncertainties were obtained by exploring the effects of the electron correlations. Moreover, we found that the relativistic nuclear recoil corrections to the mass shift factors are very small and insensitive to the electron correlations for Al*+*. These atomic parameters present in this work are valuable for extracting the nuclear electric quadrupole moments and the mean-square charge radii of Al isotopes.
I Introduction
The exotic nuclei and the nuclei close to the dripline have some peculiar structures and properties, such as larger nuclear radii and reaction cross section, etc. Investigating the properties of these nuclei, especially for the evolutionary trend of the nuclear properties along the isotope chain, is helpful for understanding the many-body interactions within the nucleus and improving the nuclear structure theory Campbell et al. (2016). To characterize the nuclear properties quantitatively, nuclear electric quadrupole moment and mean-square charge radius \langle$$r^{2}$$\rangle are required. In the chain of aluminum isotopes, the proton-rich isotope 23Al has the proton-halo structure Cai et al. (2002); Nagatomo et al. (2010), the neutron-rich isotopes 31-33Al are in the vicinity of the “island of inversion”, and the 26Al is a self-conjugate nucleus Cooper et al. (1996). Up to date, the and \langle$$r^{2}$$\rangle values are only available for a few isotopes of aluminum, and the accuracies are not high enough except for 27Al Kellö et al. (1999); Sundholm and Olsen (1992); Angeli and Marinova (2013). For example, the value of 26Al, with an uncertainty about 12%, was deduced from the experimental hyperfine structures measured by using the atomic laser spectroscopy Levins et al. (1997) in assistance of the relation ( is the electric quadrupole hyperfine interaction constants). For the 23Al, 28Al, and 31-33Al isotopes, the nuclear electric quadrupole moments were measured by the -ray-detected nuclear magnetic resonance (-NMR) or -ray-detected nuclear quadrupole resonance (-NQR) method Nagatomo et al. (2010); Stöckmann et al. (1978); Nagae et al. (2009); Kameda et al. (2007); Nagatomo et al. (2009); Rydt et al. (2009). However, the error bars, about 20% for 23Al and 31Al, are so large that the nuclear deformation parameter cannot be determined. Compared with the nuclear quadrupole moment, the nuclear charge mean-square-root is scarce. The \langle$$r^{2}$$\rangle value of these interesting nuclei has not been reported. As a consequence, it is indispensable to determine accurately the and \langle$$r^{2}$$\rangle values of aluminium isotopes for exploring their nuclear properties.
The nuclear-model-independent data and \langle$$r^{2}$$\rangle can be extracted by combining the measured hyperfine structures and isotope shifts with theoretical prediction of the electric field gradient and the isotope shift factors. The collinear laser spectroscopy technique has been developed to measure the isotope shifts and the hyperfine structures of the exotic or radioactive nuclei with short lifetime and low production Campbell et al. (2016); Nörtershäuser et al. (2009). In addition, the high-precision measurements of hyperfine structures and isotope shifts of Al+ ions can be obtained based on the Al+ ion optical clock Xu et al. (2016), and thus can be used to verify the calculation. For these reasons, the Al+ ion can be considered as a good candidate for extracting the nuclear electric quadrupole moment and mean-square charge radius \langle$$r^{2}$$\rangle of Al isotopes.
Most of the earlier theoretical works related to the 3s3p 3s2 transitions in Al+ ions focused on transition energies and probabilities Konovalova and Kozlov (2015); Santana and Träbert (2015); M Stanek (1996); Safronova et al. (2000); Safronova (2000); Laughlin and G.A.Victor (1979); Hibbert and Keenan (1987); Brage et al. (1998); Jönsson and Fischer (1997); Zou and Fischer (2000, 2001); Fischer et al. (2006); Kang et al. (2009, 2010); Andersson et al. (2010). For the hyperfine interaction constants we only found theoretical reports by Kang et al. and Andersson et al.. In their works the off-diagonal magnetic dipole and electric quadrupole hyperfine interaction constants of the 3s3p states were calculated for investigating the effect of hyperfine interaction on the lifetime of metastable states 3s3p Kang et al. (2009); Andersson et al. (2010). While the isotope shift factors of the transitions 3s3p 3s for Al*+* ions have not been reported.
In this work, we calculated the hyperfine interaction constants of 3s3p states and the mass shift and field shift factors of the 3s2 3s3p in Al+ using the multi-configuration Dirac-Hartree-Fock (MCDHF) method. The active space approach was adopted to investigate the effects of electron correlations on the hyperfine interaction constants and the isotope shift factors in detail. Based on this, we built a computation model, which can capture electron correlations effectively and allow us to obtain high-precision atomic parameters concerned. In addition, the relativistic nuclear recoil corrections to the isotope shift factors were discussed. We hope this work could support the experimental investigations about the Al+ ion optical clock and the studies about the nuclear properties of Al isotopes.
II Theoretical method
II.1 MCDHF method
In the MCDHF method an atomic state wavefunction (ASF) is a linear combination of configuration state functions (CSFs) with the same parity , total angular momentum and it’s component along the direction P.Grant (2007),
[TABLE]
In the Eq (1) is the mixing coefficient, and represents other appropriate labeling of the CSF. The CSFs are the linear combinations of products of one-electron Dirac orbitals. In the self-consistent field (SCF) calculation both the mixing coefficients and the orbital are optimized. After the virtual orbital set is obtained more electron correlations can be included in the relativistic configuration interaction (RCI) calculations by further expanding the configuration space. In the RCI calculations all orbitals are kept frozen, and only the mixing coefficients are variable. The Breit interaction in the low-frequency approximation,
[TABLE]
is also taken into account in the RCI computation Li et al. (2012).
II.2 Hyperfine interaction
The hyperfine interaction is caused by the interaction between the electrons and the electronmagnetic multipole moments of the nucleus, and its Hamiltonian can be represented using the spherical tensor operators and Lindgren (1984),
[TABLE]
Here, and represent the magnetic dipole and electric quadrupole interactions, respectively, and the higher-order terms are tiny and neglected in this work. For a -electron atom the electronic tensor operators and are the sums of the one-electron operators,
[TABLE]
and
[TABLE]
Here, is the imaginary unit, is the fine-structure constant, and are spherical tensor operators, and is the orbital angular momentum operator. The nuclear tensor operators and are related to magnetic dipole moment and electronic quadrupole moment of the nucleus with spin through Jönsson et al. (1996):
[TABLE]
and
[TABLE]
In the first-order perturbation approximation, the magnetic dipole and the electric quadrupole hyperfine interaction constants and are Jönsson et al. (1996)
[TABLE]
and
[TABLE]
II.3 Isotope shifts
The isotope shift is composed of the mass shift (MS) and the field shift (FS). The former is caused by the motion of nucleus with finite mass and the latter by the nuclear distribution. For a transition the isotope shift is
[TABLE]
In the formulas above, the and are the nuclear masses of the isotopes and , respectively. In addition, . The mass shift (MS) can be separated into the normal mass shift (NMS) and the specific mass shift (SMS). The relativistic nuclear recoil Hamiltonian and correspond to one-body and two-body relatively mass shift operator, respectively Palmer (1987); Torbohm et al. (1985),
[TABLE]
For a given level, the normal mass shift factors and the specific mass shift factors are defined as Li et al. (2016)
[TABLE]
and
[TABLE]
For a transition, the is the difference of the mass shift factors between the upper () and lower () levels, and .
The field shift (FS) factor for a given transition is expressed as Blundell et al. (1987)
[TABLE]
where, is the change of the total electronic probability density at the origin
[TABLE]
III Computational model
Generally the precision of the calculated atomic parameters mainly depend on description of electron correlations in the atomic system. In the frame of the MCDHF method one can consider electron correlations systematically adopting the active space approach. According to the perturbation theory the electron correlations can be divided into the first-order and the higher-order correlations. The first-order correlation effects are captured by the CSFs generated through the single (S) and double (D) excitations from the occupied orbitals in the single reference configuration (SR) sets. In order to capture the electron correlations efficiently the occupied orbitals in the reference configuration are separated into the valence and the core orbitals. Therefore, the first-order correlation is composed of the correlation between the valence electrons (VV correlation), the correlation between valence and core electrons (CV correlation), and the correlation between core electrons (CC correlation). The one beyond the first-order correlation is defined as the higher-order correlation and captured by the configuration space expanding from the multireference (MR) configurations set, i.e. the “MR-SD” model. It was shown that this method is capable of accounting for the electron correlation for the complex ions and atoms Li et al. (2012); Filippin et al. (2016a, b); Carette and Godefroid (2013); Carette et al. (2010); Andersson et al. (2016).
III.1 Capture of the first-order correlations
In our calculations, the atomic state wave function of even 3s2 and odd 3s3p states were optimized separately. We treated the 1s, 2s and 2p orbitals as the core, and the others as the valence orbitals. The VV and CV correlations were taken into account in the SCF calculations. As shown in the second column of Table 1 and Table 2 the active orbitals (AO) were enlarged layer by layer in order to monitor the convergence of the physical concerned. These steps were labeled by CV1s-, where the subscript means the CV correlation between the 1s core and the valence orbitals was involved. and is the maximum principal quantum number and the maximum angular quantum number of the outermost active orbitals at each step, respectively. In our test calculations it was found that the contributions from orbitals with large orbital angular momentum (such as 7i, 7h, 8g orbitals, etc.) are fractional to atomic parameters concerned, so these orbitals were not included in the set of active orbitals. The final set of active orbitals in our calculations was composed of eleven layers of orbitals with , seven layers with , three layers with and one layer with . The number of CSFs for 3s2 and 3s3p states in each step are also presented in Table 1 and Table 2. The SCF calculations started from the Dirac-Hartree-Fock approximation (labeled as DF in the table), in which the occupied orbitals in the reference configurations were optimized as spectroscopic. Subsequently, these orbitals were kept frozen, and only the added orbitals in the active set were variable. At last, the orbitals sets formed in the step CV1s- were fixed in subsequent RCI calculations in which the CC2s and the higher-order correlations were included.
The correlation in the core (labeled as CC2s) was taken into account by allowing the single and double excitations from the 2s and 2p core orbitals to the largest active set. So far, all the first-order correlation were included in our calculations except the correlation between the core orbitals 1s (CC1s correlation). In fact the CC1s correlation is negligible and will be discussed in section below.
III.2 Capture of the higher-order correlations
The higher-order correlation can be captured in two ways. One is to add the CSFs generated by the triple (T) and quadruple (Q) substitutions from the single reference configuration, and the other is to include those produced through the SD substitutions from the multireference configurations. The first way is impracticable for complex atoms, since the configuration space will be expanded too rapidly. Moreover, it is unnecessary to capture all TQ substitutions in practical calculations, because the contributions from most of them are tiny. Actually, the SD excitations from the multireference configurations set is equivalent to the restricted TQ excitations from the single reference configuration, but reserve the important TQ excitations by properly selecting the reference configurations. Therefore, the most important higher-order correlation can be captured by using the MR-SD model. From the physical viewpoint the multi-reference configuration set is composed of the CSFs with large mixing coefficient in the first-order configuration space. In this case, we selected the dominant CSFs in the configuration spaces obtained with the CC2s model to form the multireference configuration set, since the contributions of the CC1s correlations to the mixing coefficient of the dominant CSFs are negligible. The dominant CSFs can be identified according to the weight factor that is defined as
[TABLE]
Here, the sum extend over the CSFs that belong to the same configuration Carette and Godefroid (2013). The Table 3 shows the weight of the selected configurations. It should be noted that the weight factors in this table were obtained by summing the configuration with , because the mixing coefficients only have tiny effects on the calculated weight factor of the configuration.
In order to explore the convergency of the higher-order correlation effect, we expanded the multireference configurations set in terms of the condition , and , i.e. the MR1, MR2 and MR3 sets in Table 3, respectively. It should be emphasized that only the higher-order correlation between valence orbitals (3s and 3p) and the core orbitals (2s and 2p) was considered in our calculations. In addition, in order to ensure the convergency of the atomic parameters concerned for a given MR sets, the orbital in the active set were added layer by layer as mentioned earlier. For example, five layers of correlation orbitals were added for capturing adequately the higher-order correlation under the condition of . This was marked as MR1-, in which is the maximum principal quantum number of the active “s” orbitals for the corresponding MR1 model. The similar regulation was used for the MR2 and MR3 models as well. Moreover, we should point out that in the selection of the multireference configurations we balanced the weights of the given configurations for and states since these levels were optimized on a common orbital basis set. It means that the 2s22p63s5p and 2s22p63d5p were added in the MR2 set, since their weights for state satisfy . At last, the contribution of the Breit interaction was evaluated in the DF model. The calculations in this work were carried out by using the GRASP2K code Jönsson et al. (2007, 2013).
IV Results and discussion
IV.1 Hyperfine interaction constants of 3s3p
Figure 1 shows the magnetic dipole () and the electric quadrupole () hyperfine interaction constants of the 3s3p states in Al+ ions as functions of the configuration spaces. It can be seen that seven layers of correlation orbitals are enough to capture the VV and CV electron correlation effects in the SCF calculations. The extra three layers were augmented for adequately describing the core-core correlation in the subsequent RCI computation. In addition, the relatively large oscillation for the electric quadrupole hyperfine interaction constants of the 3s3p state occurs because the scale in this figure is small. The contributions of CC2s correlations to the and constants are 4% and 9% for the state, respectively, and 9% and 5% for the state. The higher-order correlation on the hyperfine interaction constants were accounted for by the MR model. These contribute to the magnetic dipole and the electric quadrupole hyperfine interaction constants of state by 0.6% and 3%, respectively, and of the state by 8% and 5%. It is worth noting that the higher-order correlation effects counteract the core-core correlation effects. For instance, the effects of the CC2s correlation on the magnetic dipole hyperfine interaction constants of the states are offset by the higher-order correlation, so the final results obtained with the MR model are in good agreement to the one from the CV models. In fact, Engles and Jacek et al. found this phenomenon in their works Bieroń et al. (2009); Engels (1993), that is, the CC correlations always make the agreement worse between the calculated hyperfine interaction constants and the experimental values, and these discrepancies can be offset by the higher-order correlation. Hence, for the calculation asked for high precision, the CC correlation and the higher-order correlation are indispensable, which is also allowed us to evaluate the uncertainties in the calculation.
Kang et al. and Andersson et al.Kang et al. (2009); Andersson et al. (2010) calculated the off-diagonal hyperfine interaction constants using the MCDHF method in order to investigate the influence of hyperfine interaction on the lifetime of metastable states 3s3p . For confirming our computation models are reliable, we made comparisons for the off-diagonal hyperfine interaction constants between our and their results in Table 4. As can be seen from this table, excellent agreement was found between ours and Andersson et al.. The less good agreement for Kang et al. results is caused by the fact that the core-valence correlation between the 1s orbitals and the valence orbitals, the core-core correlation and the higher-order correlation were not considered.
Table 5 shows the calculated hyperfine interaction constants corresponding to the MR3- model and their uncertainties (in the parenthesis) of the 3s3p states in Al+ ions. Generally, the uncertainties result from electron correlations neglected in the computational models and physical effects. In this work the VV, CV1s, CC2s and the main part of higher-order correlations were taken into account systematically. According to the convergence trend illustrated in Figure 2, there are 0.8% uncertainties in the calculations. The CC1s correlation and the higher-order correlation involving the 1s orbital were neglected in our computational models. From our test we found that the effects of the CC1s correlations on the hyperfine interaction constants under investigation are less than . As discussed above, the effects of the core-core correlation related to the shell on the hyperfine interaction constants almost counteract the higher-order correlation effects. We speculated that the cancellation also arises between the CC1s and the related higher-order correlation. Therefore, the uncertainty due to the neglected CC1s and the higher-order effects should be less than 1%. Additionally, the Breit interaction was ignored in the present calculation, which gives rise to 0.05% uncertainty. The total uncertainty of the and is 0.76% and 1.32% for the state, and 1.42% and 1.36% for the state.
IV.2 Isotope shift factors of transitions 3s2 3s3p
Table 6 shows the isotope shift factors, including the relativistic normal and specific mass shifts factors , and the field shift factors , for the transitions 3s2 3s3p as functions of computational models. From this table we can see that the mass shift factors are more sensitive to the electron correlations than the field shift factors, especially for the specific mass shift factors. For instance, from the model DF to CV1s-, the change of the reaches 28% for the transition 3s2 3s3p , and 48% for 3s2 3s3p . While the reduced by a factor of 2 for the transition 3s2 3s3p , and a factor of 6 for the transition 3s2 3s3p , respectively. In contrast, the change for is marginal, less than 2%. The high sensitivity of the SMS factors to the electron correlations is shown again when the CC2s correlations were included in RCI computations. The influences of the CC2s correlation on the is about 16%, and the change is about three times for the of the transition 3s2 3s3p . For the transition 3s2 3s3p the effects of CC2s correlation and related higher-order correlation on the isotope shift factors are opposite, which is similar to the case of the hyperfine interaction constants. However, the CC2s and related higher-order correlations both decrease the of the transition 3s2 3s3p .
We display the relativistic and nonrelativistic NMS and , the SMS and, the FS factors and their uncertainties (in parenthesis) in Table 7 for the transitions 3s2 3s3p . In the non-relativistic approximation, the NMS factors for a given transition are proportional to its transition frequency , i.e. the scaling law Godefroid et al. (2001). Using the experimental transition energies we obtained the non-relativistic normal mass shift factor 615.67 (GHz u) for the 3s2 3s3p transition and (GHz u) for the 3s2 3s3p transition, respectively. Compared with our ab initio calculation, the discrepancies are less than 3%, which due not only to the neglected CC1s and related higher-order correlations but also to the relativistic effect in the atomic state wave functions. The linear correlation between the convergency of the calculated transition energy and the MS factors for a given transition has been found in the B+, C- and Cl- ions Liu et al. (2017); Carette and Godefroid (2013, 2011), which also occurs in the present case. We deduced from the linear correlation that the SMS factors are changed in a right direction with the MR calculations, since the transition energies in MR3 model agree better with the experimental value than those in MR1 and MR2 models. In addition, from Table 7, it was found that for the transition 3s2 3s3p the contribution of the relativistic nuclear recoil corrections to the NMS and SMS factors are 24(GHz u) and 13(GHz u), respectively. For the other transition it has similar value. For the total MS factors of the transitions 3s2 3s3p , the effects of the relativistic nuclear recoil corrections are not more than , which is in our expectation that the relativistic nuclear recoil corrections for ions with are small. In addition, we noticed that the effects of the relativistic nuclear recoil corrections on the NMS and SMS factors are insensitive to the electron correlations, which is illustrated in Figure 2.
For the isotope shift factors, the uncertainties from the VV, CV1s and CC2s correlations reach a satisfactory level. Specifically, for the NMS factors of two transitions it is less than 1%, and for the FS factors less than 0.1%. The uncertainties of the SMS factors resulting from the VV, CV1s and CC2s correlations were controlled less than 2% for the transition 3s2 3s3p and 7% for 3s2 3s3p . Since the limited computing resource, we cannot further expand the configuration space. Therefore, it is difficult to estimate accurately uncertainties of the mass shift factors in the MR computations due to neglected higher-order correlations. Roughly, the uncertainties reach around 9% and 7% for the NMS and SMS factors, respectively, in the 3s2 3s3p transition. For the 3s2 3s3p transition, the errors are about 2% for the NMS and 16% for the SMS. Compared with the mass shift factors, the FS factors are stable with expansion of the configuration space. We estimated the uncertainties for the FS factors in the MR calculations to be around 1% for these two transitions. Additionally, we have tested the effects of the CC1s correlation on the IS factors, and found that for the NMS factors the uncertainties are not more than , but for the SMS factors they are about and for the 3s2 3s3p and 3s2 3s3p transitions, respectively. Nevertheless, the effect of the CC1s correlation on the FS factors is fractional (less than 1%). The Breit interaction corrections to the isotope shift factors, estimated in the DF calculation, are less than 3% for all physical quantities under investigation. To sum up, the total uncertainties are about 15% for the NMS factors, 17% for the SMS factors, and 3% for the FS factor in 3s2 3s3p transition, and 4%, 47% and 2% for the other transition. The calculated NMS, SMS and FS factors together with their uncertainties are listed in Table 7.
V Conclusion
In summery, we calculated the hyperfine interaction constants and the isotope shift factors involving the 3s2 and 3s3p states in Al+ ions using the MCDHF method and the active space approach. We have discussed the effects of the electron correlations and the Breit interaction on the atomic parameters concerned. In this case, we found that for the hyperfine interaction constants the contribution of the higher-order correlation is opposite with CC correlation. Based on the discussion about the contribution of the electron correlations and the Breit interaction, we have obtained reliable uncertainties for our calculated results. For the hyperfine interaction constants, the uncertainties are less than 1.5%. For the isotope shift factors, the uncertainties are 14% and 4% for the NMS factors, 17% and 47% for the SMS factors, and 3% and 2% for the FS factors in transition 3s 3s3p and 3s 3s3p , respectively. Although the higher-order correlations change the level structure slightly, these effects on the hyperfine interaction constants and the isotope shift factors are indispensable. Therefore, it is necessary to include the CC correlation and related higher-order correlation in computational models in order to achieve high accuracy of atom parameters. In addition, for the Al+ ion the effect of the relativistic nuclear recoil correction on the mass shift factor is small (less than 4%) and insensitive to the electron correlations.
With respect to the fact that there are relatively large uncertainties in the calculation of mass shift factors, the partitioned correlation function interaction (PCFI) approach Verdebout et al. (2013) will be a promising method for calculating the isotope shifts more accurately, which can capture the different electron correlations flexible and accurate.
During the review of our article, a paper about hyperfine-mediated electric quadrupole shifts in Al+ and In+ ion clocks was just published Beloy et al. (2017). In this paper, the diagonal magnetic dipole and electric quadrupole hyperfine interaction matrix elements of the 3s3p state for Al+ were calculated by using the method of configuration interaction plus many-body perturbation theory (CI+MBPT). For the diagonal magnetic dipole hyperfine interaction matrix element, our value is in perfect agreement with theirs. While the consistence is less good for the diagonal electric quadrupole hyperfine interaction matrix element, and the difference is around 35%.
Acknowledgements.
This work has been supported by the Nation Natural Science Foundation of China under Grant Nos.11404025, 91536106, U1331122 and U1530142, the China Postdoctoral Science Foundation under Grant No.2014M560061, the Young Teachers Scientific Research Ability Promotion Plan of Northwest Normal University (NWNU-LKQN-15-3).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Campbell et al. (2016) P. Campbell, I. Moore, and M. Pearson, Progress in Particle and Nuclear Physics 86 , 127 (2016) . · doi ↗
- 2Cai et al. (2002) X. Z. Cai, H. Y. Zhang, W. Q. Shen, Z. Z. Ren, J. Feng, D. Q. Fang, Z. Y. Zhu, W. Z. Jiang, Y. G. Ma, C. Zhong, W. L. Zhan, Z. Y. Guo, G. Q. Xiao, J. S. Wang, Y. T. Zhu, J. C. Wang, J. X. Li, M. Wang, J. F. Wang, Z. J. Ning, Q. J. Wang, and Z. Q. Chen, Phys. Rev. C 65 , 024610 (2002) . · doi ↗
- 3Nagatomo et al. (2010) T. Nagatomo, H. Ueno, M. Mihara, K. Matsuta, A. Yoshimi, Y. Ichikawa, K. Yamada, H. Kawamura, A. Ozawa, T. Moriguchi, Y. Ishibashi, K. Asahi, M. Uchida, K. Suzuki, T. Inoue, Y. Hasama, H. Iijima, T. Sumikama, M. Fukuda, and T. Minamisono, Hyperfine Interactions 198 , 103 (2010) . · doi ↗
- 4Cooper et al. (1996) T. G. Cooper, J. Billowes, P. Campbell, and M. R. Pearson, Journal of Physics G: Nuclear and Particle Physics 22 , 99 (1996) .
- 5Kellö et al. (1999) V. Kellö, A. J. Sadlej, P. Pyykkö, D. Sundholm, and M. Tokman, Chemical Physics Letters 304 , 414 (1999) . · doi ↗
- 6Sundholm and Olsen (1992) D. Sundholm and J. Olsen, Phys. Rev. Lett. 68 , 927 (1992) . · doi ↗
- 7Angeli and Marinova (2013) I. Angeli and K. Marinova, Atomic Data and Nuclear Data Tables 99 , 69 (2013) . · doi ↗
- 8Levins et al. (1997) J. M. G. Levins, J. Billowes, P. Campbell, and M. R. Pearson, Journal of Physics G: Nuclear and Particle Physics 23 , 1145 (1997) .
