`Journal of Atomic, Molecular, and Optical PhysicsVolume 2012 (2012), Article ID 569876, 6 pageshttp://dx.doi.org/10.1155/2012/569876`
Research Article

## Energies, Fine Structures, and Hyperfine Structures of the States for the Beryllium Atom

School of Physics, Beijing Institute of Technology, Beijing 100081, China

Received 12 July 2012; Revised 15 August 2012; Accepted 27 August 2012

Academic Editor: Derrick S. F. Crothers

Copyright © 2012 Chao Chen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Energies and wave functions of the states for the beryllium atom are calculated with the full-core plus correlation wave functions. Fine structures and hyperfine structures are calculated with the first-order perturbation theory. For the state, the calculated energies, fine structure, and hyperfine structure parameters are in good agreement with the latest theoretical and experimental data in the literature; it is shown that atomic parameters of the low-lying excited states for the beryllium atom can be calculated accurately using this theoretical method. For the () states, the present calculations may provide valuable reference data for future theoretical calculations and experimental measurements.

#### 1. Introduction

In recent years, studies of energies, fine structures, and hyperfine structures of the low-lying excited states for the beryllium atom [110] have been of great interest to spectroscopists because there are many strong optical transitions suitable for spectral and hyperfine structure measurements. On the other hand, studies of the low-lying excited states for the beryllium atom play an important role in developing the excited state theory of multielectron atoms and better understanding the complicated correlation effects between electrons. The fine structure comes from the spin-orbit, spin-other-orbit, and spin-spin interactions. The hyperfine structure of atomic energy levels is caused by the interaction between the electrons and the electromagnetic multipole moments of the nucleus. The leading terms of this interaction are the magnetic dipole and electric-quadrupole moments. The fine and hyperfine structure is sensitive to the correlation effects among electrons. Experimentally, some properties of the atomic nucleus can be obtained by investigating the hyperfine structure of the atomic energy levels. The nuclear electric-quadrupole moment, which is difficult to measure directly with nuclear physics techniques, can be determined using the measured hyperfine structure and the accurate theoretical results.

The state of the beryllium atom is of interest since it is the lowest excited state in which hyperfine effects can occur, and the ground state has no hyperfine splitting because it is . It is generally a very demanding task to calculate hyperfine structure accurately. Polarization of the closed shells in the   core, due to the Coulomb interaction with open shells, can have a large effect on the hyperfine structure. Up till now, the most sophisticated theoretical calculations of the hyperfine structure parameters for the state of the Be atom have been carried out using linked-cluster many-body perturbation (LC MBPT) theory [5, 6], Hartree-Fock and CI allowing all SD excitations to correlation orbitals of Slater type by Beck and Nicolaides [7], as well as multiconfiguration Hartree-Fock (MCHF) method [8, 9]. Experimentally, the magnetic dipole and electric-quadrupole hyperfine constants have been determined very accurately with the atomic-beam magnetic-resonance technique [10] for the state in beryllium. To the best of our knowledge, few results on energies, fine structures, and hyperfine structures have been investigated for the states of the beryllium atom due to the restriction of resolution from experiments and the numerical unsteadiness in theoretical calculations.

An elegant and complete variation approach, namely, the full core plus correlation (FCPC) method, has been developed by Chung [11, 12]. This method has been successfully applied to three- and four-electron systems, with the core. Many elaborate calculations, especially for the dipole polarizabilities [13], quadrupole and octupole polarizabilities [14], and total atomic scattering factors [15], show that FCPC wave functions have a reasonable behavior over the whole configuration space for three-electron systems. This method has also been used to calculate the hyperfine structure of the and states () for the lithium isoelectronic sequence; the results are in good agreement with the Hylleraas calculations and with the experiment data [16]. As is well known, theoretical calculations of the hyperfine structure parameters depend sensitively on the behavior of the wave function in the proximity of the nucleus. In addition, core polarization effects for the low states need to be included in the nonrelativistic wave function. It would be interesting to find out whether the FCPC wave function can also be successful for calculating hyperfine structure parameters of low-lying excited states for the beryllium atom. In this work, the FCPC wave functions are carried out on the () states of the beryllium atom. The energies, fine structures, and hyperfine structures are calculated and compared with the data available in the literature. The purpose of this work is to explore the capacity of the FCPC wave function to calculate the atomic parameters of the low-lying excited states for the beryllium atom and provide more reliable theoretical data to stimulate further experimental measurements.

#### 2. Theory

According to the FCPC method [11, 12], the wave function for the four-electron state can be written as where is an antisymmetrization operator. is a predetermined -core wave function which is represented by a CI basis set, the angular part is (1,2) is a two-electron singlet spin function. The linear and nonlinear parameters in (2) are determined by optimizing the energy of the two-electron core. The factor represents the wave function of the two outer electrons which is given by the angular part is The latter wave function of (1) describes the core relaxation and the intrashell electron correlation in the four-electron system. It is given by where A different set of   is used for each . The angular part is

To simplify notation, this angular function is simply denoted as with the understanding that and couple into , the total orbital angular momentum. There are three possible spin functions for the state, namely,

For the radial basis functions of each angular-spin component, a set of linear and nonlinear parameters is chosen. These parameters are determined in the energy optimization process. For each set of  , , , and , we try all possible and and keep the ones which make significant contribution to the energy in (1).

The fine structure perturbation operators [1, 2] are given by where the spin-orbit, spin-other-orbit, and spin-spin operators are

To calculate the fine structure splitting, the coupling scheme is used:

The fine structure energy levels are calculated by first-order perturbation theory

For an -electron system, the hyperfine interaction Hamiltonian can be represented as follows [17, 18]: where and are spherical tensor operators of rank in the electronic and nuclear spaces, respectively. The term represents the magnetic-dipole interaction between the magnetic field generated by the electrons and nuclear magnetic dipole moments, and the term the electric quadrupole interaction between the electric-field gradient from the electrons and the nonspherical charge distribution of the nucleus. The contributions from higher-order terms are much smaller and can often be neglected.

In the nonrelativistic framework, the electronic tensor operators, in atomic units, can be written as where is the orbital electron factor, and is the electron spin factor. is the nuclear mass. The tensor is connected to the spherical harmonics by

The hyperfine interaction couples the electronic angular momenta J and the nuclear angular momenta I to a total angular momentum F = I + J. The uncoupling and coupling hyperfine constants are defined in atomic units as [17, 18]: where and . In these expressions, is the nuclear magnetic moment and is the nuclear electric quadrupole moment. is the nuclear spin, and is the atomic electronic angular moment.

#### 3. Results and Discussions

In order to achieve accurate calculation results for various properties of the low-lying excited states for the beryllium atom, the choice of basis function with sufficiently high quality is critical and it is our major concern. The seven components (0,0), (1,1), (2,2), (3,3), (4,4), (5,5), (6,6) altogether 159 terms are used for the core. The in (1) has four angular components, is summed from 0 to 3 with the angular components (0,1), (1,2), (2,3), and (3,4), and the number of terms in ranges from 36 to 15. Most of the other correlation effects are included in (6), which accounts for the intershell as well as the intrashell correlations. Many relevant angular and spin couplings are important for the energy, these basis functions are tried to include in (1, 2, 3, 4) with significant energy contribution. For each set of orbital angular momenta , , , and , there could be several ways to couple this set into the desired total orbital angular momentum. In this work, for states, the important angular series (, , , ) are , , , and so forth. In both cases, the value of is from 0 to 6, as the energy contribution from set with is small and negligible. In order to get the high-quality wave function, the number of angular-spin components in the wave functions ranges from 15 to 66, and the number of terms in the of (6) is about 790. The linear and nonlinear parameters are individually optimized in the energy minimization process. Using the Rayleigh-Ritz variational method, the basic wave function and the corresponding eigenvalue are determined.

Nonrelativistic energies of the () states for the beryllium atom are given in Table 1. As Table 1 shows, for the () states, the nonrelativistic energies in this work are lower and better than those of Hibbert and Weiss [1, 2], the improvement ranging from 0.0479 a.u. to 0.0473 a.u. Hibbert and Weiss reported a set of large-scale configuration interaction (CI) calculations for the () states, which can give an accurate approximation for each state, but it may tend to obscure the global picture of the spectrum which is so transparent in the other approach. The work of Hibbert and of Weiss did not include any intrashell correlation in the shell, as the calculations were of transitions in the outer subshells. The correlation energy of the shell is almost independent of the nuclear charge and also of the number of additional electrons outside the shell. For Be, it is about 0.0457 a.u. and this accounts for the main difference between earlier work and the present; more accurate results are presented in Table 1. Of course, for the calculation of hyperfine parameters, correlation within the shell is crucial in obtaining accurate hyperfine parameters, and this has been achieved in the present work. For the state, the present calculation from the FCPC method is also lower than the result of Weiss [2].

Table 1: Nonrelativistic energies of the states for the beryllium atom (in a.u.).

If including the effects of the spin-orbit, spin-other-orbit, and spin-spin interactions, the energies of the fine structure resolved levels are obtained. In this work, the fine structure splittings of the triplet states are calculated with the , , and operators using the first-order perturbation theory. Table 2 gives the fine-structure splittings of the () states for the beryllium atom. The experimental Be splittings are 2.35 () and 0.64 cm−1 [3]. They agree with our prediction 2.36 and 0.64 cm−1. Although many theoretical studies have been done on the BeI excited systems, the published theoretical fine structure results are scarce. One exception is Laughlin, Constantinides, and Victor [4]. They use a model potential calculation and predict the splittings to be 2.53 and 0.71 cm−1 for the state, which should be considered as quite good in view of the simplicity in their computation and fall in experimental uncertainties. Present calculations for this state are more accurate due to correlation effect well described in this method. For the experiment, the splitting of is not resolved. But the splitting from the state to the , 0 is determined to be 0.37 cm−1. In this work, the calculated splittings are 0.35 and 0.092 cm−1. This implies that the predicted splitting from to the center of gravity of and 0 should be 0.373 cm−1. It agrees with the experiment. The good agreement with experiment could be used as the indication of the accuracy of the wave function constructed here. For the state, our calculated splittings are hoped to offer reference for further experimental measurements.

Table 2: Fine structure splittings of the states for the beryllium atom (in cm−1).

The hyperfine structure parameters of the states for the beryllium atom are calculated in this work: Fermi contact , the spin dipolar , the orbital , and the electric quadrupole . In the present calculation, , = −1.177492 nm, for Be are taken from [19]. The hyperfine interaction in the state for the beryllium atom is of interest since it is the lowest excited state in which hyperfine effects can occur, which has been studied over the past four decades [510]. Table 3 gives the hyperfine structure parameters of the state for the beryllium atom through the FCPC wave function to compare with data in the literature. As can be seen from Table 3, the present results for hyperfine structure parameters are better than the earlier theoretical results [57] in whole. The present calculations also agree with the results by FE MCHF (finite-element multiconfiguration Hartree-Fock) method [8] to two significant figures. The calculated Fermi contact term in this work differs from the results from the latest calculation through MCHF method [9] by only 0.07%, and the differences for the other terms are on the order of a few parts in a thousand. This means that the wave function used in the present work is reasonable and accurate in the full configuration space. The hyperfine coupling constants are also listed in Table 3 to compare with results from other calculations and experiments. Our calculated hyperfine coupling constants agree perfectly with the experimental value [10] to four significant figures. That is also true for the MCHF calculation of [9]. It is shown that hyperfine structure parameters of the low-lying excited states for the beryllium atom can be calculated accurately using the present FCPC wave function. For the ( states, to the best of our knowledge, there is no report on hyperfine structure parameters in the literature. The present predictions for the hyperfine structure parameters and coupling constants are listed in Table 4, which may provide valuable reference data for other theoretical calculations and experimental measurements in future.

Table 3: The hyperfine structure parameters (in a.u.) and coupling constants (in MHz) of the state for the beryllium atom.
Table 4: The hyperfine structure parameters (in a.u.) and coupling constants (in MHz) of the states for the beryllium atom.

#### 4. Summary

In this work, energies, fine-structure splittings, and hyperfine structure parameters of the () states for the beryllium atom are calculated with the FCPC wave functions. The obtained nonrelativistic energies are much lower than the previous published theoretical values. The calculated fine structure splittings are in good agreement with experiment. For the state, the calculated hyperfine structure parameters are in good agreement with the latest theoretical and experimental data in the literature; it is shown that hyperfine constants of the low-lying excited states for the beryllium atom can be calculated accurately using this kind of wave function. For other states, the present predicted hyperfine structure parameters may provide valuable reference data for future theoretical calculations and experimental measurements.

#### Acknowledgments

The author is grateful to Dr. Kwong T. Chung for his computer code. The work is supported by National Natural Science Foundation of China and the Basic Research Foundation of Beijing Institute of Technology.

#### References

1. A. Hibbert, “Oscillator strengths of transitions involving $2s3{l}^{3}L$ states in the beryllium sequence,” Journal of Physics B, vol. 9, no. 16, pp. 2805–2811, 1976.
2. A. W. Weiss, “Calculations of the $2sn{s}^{1}S$ and $2p3{p}^{3,1}P$ Levels of Be I,” Physical Review A, vol. 6, no. 4, pp. 1261–1266, 1972.
3. L. Johansson, “The spectrum of the neutral beryllium atom,” Arkiv For Fysik, vol. 23, pp. 119–128, 1962.
4. C. Laughlin, E. R. Constantinides, and G. A. Victor, “Two-valence-electron model-potential studies of the Be I isoelectronic sequence,” Journal of Physics B, vol. 11, no. 13, pp. 2243–2250, 1978.
5. S. N. Ray, T. Lee, and T. P. Das, “Many-body theory of the magnetic hyperfine interaction in the excited state ($1{s}^{2}2s2p$${}_{}{}^{3}{P}_{}^{}$) of the beryllium atom,” Physical Review A, vol. 7, no. 5, pp. 1469–1479, 1973.
6. S. N. Ray, T. Lee, and T. P. Das, “Study of the nuclear quadrupole interaction in the excited ($2{}_{}{}^{3}{P}_{}^{}$) state of the beryllium atom by many-body perturbation theory,” Physical Review A, vol. 8, no. 4, pp. 1748–1752, 1973.
7. D. R. Beck and C. A. Nicolaides, “Fine and hyperfine structure of the two lowest bound states of Be- and their first two ionization thresholds,” International Journal of Quantum Chemistry, vol. 26, supplement 18, pp. 467–481, 1984.
8. D. Sundholm and J. Olsen, “Large MCHF calculations on the hyperfine structure of Be(3PO): the nuclear quadrupole moment of 9Be,” Chemical Physics Letters, vol. 177, no. 1, pp. 91–97, 1991.
9. P. Jönsson and C. F. Fischer, “Large-scale multiconfiguration Hartree-Fock calculations of hyperfine-interaction constants for low-lying states in beryllium, boron, and carbon,” Physical Review A, vol. 48, no. 6, pp. 4113–4123, 1993.
10. A. G. Blachman and A. Lurio, “Hyperfine structure of the metastable ($1{s}^{2}2s2p$)${}_{}{}^{3}P$ states of ${}_{4}{}^{}{\text{Be}}_{}^{9}$ and the nuclear electric quadrupole moment,” Physical Review, vol. 153, no. 1, pp. 164–176, 1967.
11. K. T. Chung, “Ionization potential of the lithiumlike $1{s}^{2}2s$ states from lithium to neon,” Physical Review A, vol. 44, no. 9, pp. 5421–5433, 1991.
12. K. T. Chung, X. W. Zhu, and Z. W. Wang, “Ionization potential for ground states of berylliumlike systems,” Physical Review A, vol. 47, no. 3, pp. 1740–1751, 1993.
13. Z. W. Wang and K. T. Chung, “Dipole polarizabilities for the ground states of lithium-like systems from $Z=3$ to 50,” Journal of Physics B, vol. 27, no. 5, pp. 855–864, 1994.
14. C. Chen and Z. W. Wang, “Quadrupole and octupole polarizabilities for the ground states of lithiumlike systems from $Z=3$ to 20,” The Journal of Chemical Physics, vol. 121, no. 9, pp. 4171–4174, 2004.
15. C. Chen and Z. W. Wang, “Total atomic scattering factors for the ground states of the lithium isoelectronic sequence from Na8+ to Ca17+,” The Journal of Chemical Physics, vol. 122, no. 2, Article ID 024305, 5 pages, 2005.
16. X. X. Guan and Z. W. Wang, “The hyperfine structure of the $1{s}^{2}{\text{ns}}^{2}\text{S}$ and $1{s}^{2}\text{np}$${}_{}{}^{2}{P}_{}^{}$ states ($n=2,3,4$, and 5) for the lithium isoelectronic sequence,” The European Physical Journal D, vol. 2, no. 1, pp. 21–27, 1998.
17. J. Carlsson, P. Jönsson, and C. Froese Fischer, “Large multiconfigurational Hartree-Fock calculations on the hyperfine-structure constants of the 7Li $2s$${}_{}{}^{2}{S}_{}^{}$ and $2p$${}_{}{}^{2}{P}_{}^{}$ states,” Physical Review A, vol. 46, no. 5, pp. 2420–2425, 1992.
18. A. Hibbert, “Developments in atomic structure calculations,” Reports on Progress in Physics, vol. 38, no. 11, pp. 1217–1338, 1975.
19. P. Raghavan, “Table of nuclear moments,” Atomic Data and Nuclear Data Tables, vol. 42, no. 2, pp. 189–291, 1989.