On bound state computations in three- and four-electron atomic systems

A variational approach is developed for bound state calculations in three- and four-electron atomic systems. This approach can be applied to determine, in principle, an arbitrary bound state in three- and four-electron ions and atoms. Our variational wave functions are constructed from four- and five-body Gaussoids that respectively depend on six (r ₁₂, r ₁₃, r ₁₄, r ₂₃, r ₂₄, r ₃₄) and ten (r ₁₂, r ₁₃, r ₁₄, r ₁₅, r ₂₃, r ₂₄, r ₂₅, r ₃₄, r ₃₅ and r ₄₅) relative coordinates. The approach allows operating with the more than one electron spin functions. In particular, the trial wave functions for the ¹ S states in four-electron atomic systems include the two independent spin functions χ₁ = αβαβ + βαβα − βααβ − αββα and χ₂ = 2ααββ + 2ββαα − βααβ − αββα − βαβα − αβαβ. We also discuss the construction of variational wave functions for the excited 2³ S states in four- electron atomic systems.     

原子的解析波函数——Ⅱ.第二周期元素的正常态原子与离子的解析波函数与能量积分的计算

我们在文献[1]中设计了一套五个参数的变分波函数用来计算了周期表中前面十个原子的能量,所得结果比过去一些作者用四参数波函数所算得的结果为好。我们在过去计算经验的基础上,另外找到了一套特别简单的解析波函数,其形式为1s电子:ψ₁(r)=N₁e^-μar,2s电子:ψ₂(r)=N₂[(μr)e^-μr-Ne^-μar],2p电子:ψ₃(r)=N₃(μr)cosθe^-μr,ψ₄(r)=N₄(μr)sinθe^iφ-μr,ψ₅(r)=N₅(μr)sinθe^-iφ-μr,式中的α与μ为变分参数;N₁,N₂,N₃,N₄,N₅为归一化因子;N为正交化系数。μ可用解析法来决定,因而只有一个参数α要由数值法来决定。我们用这样的波函数算出了第二周期元素的正常态原子和离子(共有八十几个原子态)的各电子的各种能量积分值及总能量值,并确定了波函数的最佳参数值。其结果与五参数波函数的计算结果相比,一般相差在万分之一至千分之一的范围内,并比最近有些作者用一种三参数波函数所算的结果还好。根据这些结果,我们还讨论了Slater近似计算法的可靠程度和适用范围。     

原子的解析波函数

我们设计了一套变分波函数,用来计算了周期表中前面十个原子的能量。我们设计的单电子试探波函数具有下列形式:1s:ψ₁(r)=N₁e^-μαr[1+(μbr)²], 2s:ψ₂(r)=N₂[(μr)e^-μr-Ne^-μcr], 2p:ψ₃(r)=N₃(μdr)cosθe^-μdr, ψ₄(r)=N₄(μdr)sinθe^iφ-μdr, ψ₅(r)=N₅(μdr)sinθe^-iφ-μdr。式中的a,b,c,d及μ为五个变分参数。N₁,N₂,N₃,N₄与N₅为归一化因子;N由ψ₁与ψ₂的正交条件来决定。用这种波函数来计算原子的能量,所得的结果比莫尔斯等人(P.M.Morse,L.A.Young and E.S.Haurwitz)用他们设计的四参数波函数所算得的结果为好,更接近实验值,同时也接近于由自洽场所算出的结果。若我们的波函数中固定c等于1不变,这时就变为只有四个参数的波函数,结果仍比莫尔斯等人的好。     

氦原子和类氦离子基态能量的变分计算及相对论修正

采用一个包含坐标伸缩系数的简单有效的变分波函数,同时考虑到核的运动,利用Mathematic a语言开发了一个用变分法计算三体问题的程序,对氦原子和类氦离子(H^-,He,Li ⁺,Be⁺⁺,B³⁺,C⁴⁺,N⁵⁺,O ⁶⁺)的非相对论基态能量和解析波 函数进行了变分计算.在此基础上,对非相对论哈密顿量进行相对论和辐射修正,并考虑到有 限核电荷半径的影响,得到了氦原子和类氦离子高精度的基态能量值. 关键词： 氦原子 类氦离子 变分法 基态能量 相对论修正     

Compact and accurate variational wave functions of three-electron atomic systems constructed from semi-exponential radial basis functions

The semi-exponential basis set of radial functions [A.M. Frolov, Phys. Lett. A 374, 2361 (2010)] is used for variational computations of bound states in three-electron atomic systems. It appears that the semi-exponential basis set has a substantially greater potential for accurate variational computations of bound states in three-electron atomic systems than was originally anticipated. In particular, the 40-term Larson’s wave function improved with the use of semi-exponential radial basis functions now produces the total energy –7.4780581457 a.u. for the ground 1²S-state in the ^￥Li^\infty{\rm Li} atom (only one spin function c₁\chi_1 = aba\alpha\beta\alpha - baa\beta\alpha\alpha was used in these calculations). This variational energy is very close to the exact ground state energy of the ^￥Li^\infty{\rm Li} atom and is substantially lower than the total energy obtained with the original Larson’s 40-term wave function (–7.477944869 a.u.).     

Variational calculations on the hydrogen molecular ion

We present high-precision non-relativistic variational calculations of bound vibrational—rotational state energies for the H⁺ ₂ and D⁺ ₂ molecular ions in each of the lowest electronic states of Σ _g , Σ _u , and Π _u symmetry. The calculations are carried out including coupling between Σ and Π states but without using the Born—Oppenheimer or any adiabatic approximation. Convergence studies are presented which indicate that the resulting energies for low-lying levels are accurate to about 10^?13. Our procedure accounts naturally for the lambda-doubling of the Π _u state.     

Microscopic Analysis of 5d States Splitting and Charge Transfer Energies Dependence on Interionic Distance in Alkaline Earth Fluorides Doped with Light Trivalent Lanthanides

A first principles fully relativistic analysis (K. Ogasawara et al., Phys. Rev. B 2001, 64, 115413) of the dependence of 5d orbitals splitting (10Dq) and charge transfer (CT) energies on interionic distance has been performed for light lanthanides (Ce³⁺, Pr³⁺, Nd³⁺) in CaF₂, SrF₂, BaF₂ crystals. The salient feature of the method is that four‐component molecular orbitals (MO) composed of atomic wave functions are used as the basis set. Without any fitting parameter, the power dependencies for 10Dq and linear dependencies for the CT energies on the distance between rare‐earth (RE) ions and ligands were obtained. A comparison with experimental values is discussed.     

Some perturbation calculations for 1s2p, 1s 22p,and 1s 22s

The Dalgarno interchange theorem is used, together with the U function approach, to evaluate the first order perturbation corrections to <r ₁ ^-1+r ₂ ^-1> and <r ₁ + r ₂> for the two-electron states 1s2p ¹ P and 1s2p ³ P. The results for <r ₁ + r ₂> are extended by using a screening approximation, and are compared with the results of accurate variational calculations. The first order perturbation correction to spin-weighted expectation values of type <ΣV(r_j )s_zj > is given for the three-electron states 1s ²2p ² P and 1s ²2s ² S. The case V(r_j ) = r _j ^-1 is treated in detail.     

Scaling of cross-sections for asymmetric (e, 3e) process on helium-like ions by fast electrons

An approximate simple scaling law is obtained for asymmetric (e, 3e) process on helium-like ions for double ionization by fast electrons. It is based on the equation (Z ^′3π) exp[-Z(r₁ + r₂)],Z′ = Z – (5/16) for ground state wave function of helium-like ions and Z^′2 scaling of energies. The scaling law is found to work very well if the lower energy electron is ejected along the momentum transfer direction and the other one is ejected in the opposite direction. It also works quite well if this electron is ejected within about 90° of the momentum transfer direction with the other electron going in the opposite direction. The scaling law becomes increasingly accurate as the target nuclear charge and the energy increase.     

Variational Calculation of 4He Tetramer Ground and Excited States Using Realistic 4He–4He Potentials

We calculated binding energies and wave functions of the ⁴He tetramer ground and excited states employing various realistic ⁴He?⁴He potentials which includes the currently most accurate one with the adiabatic, relativistic, QED and residual retardation corrections. We used our Gaussian expansion method (GEM) for ab initio variational calculations of few-body systems. We found that precisely the same shape of the short-range correlation (r _ij < 4Å) in the dimer appear in the ground and excited states of trimer and tetramer. The four kinds of the binding energies of the trimer and tetramer ground and excited states, ${B_3^{(0)}, B_3^{(1)}, B_4^{(0)}}$ and ${B_4^{(1)}}$ , for the different potentials exhibit perfect linear correlations over the range of binding energies relevant for ⁴He atoms; namely, six types of the generalized atomic Tjon lines were observed.     

Atoms with vacancies

Methods of orthogonalization of the wave function for an atom with vacancies in the s shell and the energetically lower state for ns(n + 1)s²-type configurations are considered. Non-orthogonal radial orbitals are used. An equation is obtained for the overlap integral for atomic wave functions in terms of single-electron radial orbitals. The parameters of the generalized hydrogenic analytic variational radial orbitals and the energies of F, Ne⁺, and Na²⁺ in the configuration 1s2s²2p⁶ are determined. Results are compared to calculations from other authors and experimental values, demonstrating that the approximation is good for calculating atoms with vacancies.     

A Study of Bound States of Systems of Helium and Lithium Using the Method of Discrete-Variable Representation

The systems of particles He₂, ⁶Li–He, ⁷Li–He, He₃, ⁶Li–He₂, and ⁷Li–He₂, the binding energies of which are small and the bound-state wave functions of which are widely distributed in space, are considered. Because the interaction potential is weak and rather localized compared to the characteristic sizes of wave functions of these systems, the problem of an accurate determination of binding energy and wave functions is complicated. Small changes in input parameters or an inaccuracy of calculations can lead to considerable deviations of calculated results from true values. An essential part of the study is the development and application of the discrete-variable representation method. This method is based on the determination of basis functions and the nodes and weights of a quadrature formula in such way that the values of a function are zero at all these nodes but one. With this representation the time required for calculating the Hamiltonianmatrix elements is reduced several times. The binding energies of several systems consisting of helium and lithium atoms were obtained using the method of discrete-variable representation. Thanks to the application of this approach, the calculation time was significantly reduced without loss in accuracy.     

Spherically compressed helium atom described by perturbative and variational methods

We analyze the energy spectrum of the three lowest-lying S symmetry states for the spherically confined helium atom as a function of the box radius by using an approach based on a time independent perturbation theory and two variational methods. The first treatment depends on exact solutions for confined one-electron atoms, whereas in the latter two methods exponents and linear coefficients are variationally optimized via {s,t,u}-Hylleraas functions and Generalized {r₁,r₂,r₁₂}-Hylleraas basis sets that fulfill appropriate boundary conditions. Although it is found that throughout most of the box radii here analyzed the variational energies for the three states lie below those perturbatively obtained, an opposite trend occurs toward the weak and strong confinement regions for the singlet excited and triplet states, respectively.     

Dynamical coupling in the differential equations approach to atom-diatom exchange reactions

The atomic exchange reaction A + BC → AB + C is investigated quantum mechanically employing a coupled differential equations approach. The relative motion in reactant and product channels is described in the common coordinate R ₃ (the AC nuclear separation) and is developed in three-dimensional space. The total wave functions of the system are expressed as a superposition of valence bond electronic states of the initial (A, BC) and final (AB, C) configurations, with the coefficients describing the relative and internal (vibrational, rotational) nuclear motions. Choosing convenient trial functions with the appropriate boundary conditions and using the Kohn variational principle, a set of differential (rather than the usual integro-differential) equations is obtained for the relative motion wave functions in R ₃. The potential matrix elements turn out to be dynamical in that they depend on the initial k ₁ and final k ₂ wave vectors. Two-state coupled channel calculations of the differential and integral cross sections for the isotopic species D + H₂, H + H₂ and D + D₂ are presented for collision energies up to 0·8 eV.     

Two-time Green's function method in quantum electrodynamics of high-Z few-electron atoms

The two-time Green's function method in quantum electrodynamics of high-Z few-electron atoms is described in detail. This method provides a simple procedure for deriving formulas for the energy shift of a single level and for the energies and wave functions of degenerate and quasi-degenerate states. It also allows one to derive formulas for the transition and scattering amplitudes. Application of the method to resonance scattering processes yields a systematic theory for the spectral line shape. The practical ability of the method is demonstrated by deriving formulas for the QED and interelectronic-interaction corrections to energy levels and transition and scattering amplitudes in one-, two-, and three-electron atoms. Numerical calculations of the Lamb shift, the hyperfine splitting, the bound-electron g factor, and the radiative recombination cross section in heavy ions are also reviewed.     

Wave functions for two- and three-electron atoms and isoelectronic ions

We have developed simple wave functions for two- and three-electron atoms and ions, which have the correct structure when one of the electrons is far away, or when two of the particles are close to each other. These essentially parameter-free wave functions allow us to deduce fairly accurate values for the energies, , for multipolar polarizabilities of two-electron atoms and ions, and for the coefficients of the asymptotic density. Received: 5 August 1998 / Received in final form: 27 November 1998     

Relativistic corrections and very small g-shifts in solids

The relativistic Breit-Margenau correction to the Zeeman-interaction has been calculated for a group of atoms in the periodic table with Herman Skillman wave functions and potentials. The results are applied to the S-state ions Mn²⁺ and Gd³⁺, to the F-center and V_k-center, and to the shallow donors in Si.     

A complex first-principles study of the L2,3-edge XANES spectra and crystal field effects for divalent 3d ions in cubic ZnS

A detailed first-principles analysis of the L_2,3 X-ray absorption near edge structure (XANES) spectra, crystal field strength 10Dq, covalent effects and molecular orbitals (MO) position for all divalent 3d ions from Sc²⁺ to Cu²⁺ in cubic ZnS is performed in the present paper. The calculations were done in the framework of the first-principles fully relativistic discrete variational multi-electron (DVME) based on numerical solution of the Dirac equation with the local density approximation. As a result of the performed calculations, the L_2,3 XANES spectra for all considered ions were calculated and assigned in terms of the electron configurations involved into the absorption transitions; fairly good agreement with available experimental data for Ti²⁺, Mn²⁺, Fe²⁺, Co²⁺ and Ni²⁺ is demonstrated. Experimental XANES spectra for Sc²⁺, V²⁺, Cr²⁺, Cu²⁺ not reported previously were also calculated for the sake of completeness of the present study to enable a systematic analysis of all calculated results for the whole series considered. It was shown that the L₃ and L₂ bands shift to the higher energies on increasing a 3d ion atomic number. In addition, the separation between the L₃ and L₂ bands, the crystal field strength 10Dq and mixture between the 3d ion and sulfur wave functions increase along the considered series, from Sc²⁺ to Cu²⁺. On the other hand, all orbitals of the 3d ions systematically lower down in the same direction. The above formulated trends were confirmed by the experimental data on the crystal field splittings and nephelauxetic effect.     

Simulation of the fben d transitions of lanthanide ions in YPO4 using quantum-chemical calculations

We constructed an effective one-electron Hamiltonian by using the 4f/5d energies and eigenvectors obtained from the first-principles calculation with the relativistic self-consistent discrete variational Slater software package (DV-Xα). From the effective Hamiltonian, we obtained the crystal-field and spin-orbit interaction parameters for the 4f and 5d electrons of lanthanide ions (Ce³⁺, Pr³⁺, Nd³⁺ and Eu³⁺⁾ doped in YPO₄, and these parameters were used to calculate the 4f^N-4f^N-15d transition. Comparison with experiments shows that the obtained parameters are reasonable and the excitation spectra can be well predicted.     

Fully relativistic analysis of the covalence effects for the isoelectronic 3d ions (Cr, Mn, Fe) in SrTiO3

Fully relativistic calculations of the energy levels, absorption spectra, molecular orbitals (MO) compositions, covalence effects and energies of the charge transfer (CT) transitions for three isoelectronic ions Cr³⁺, Mn⁴⁺, Fe⁵⁺ in the SrTiO₃ (STO) crystal have been performed. The recently developed first-principles approach to the analysis of the absorption spectra based on the first principles discrete variational multi-electron method (DV-ME) (K. Ogasawara et al., Phys. Rev. B 64 (2001) 115413) was used in the calculations. As a result, energy levels of the above ions, their absorption spectra and energies of the lowest CT transitions were all calculated. By performing analysis of the MO population, it was shown that the degree of covalence of the chemical bonds between 3d ions and oxygen ions in SrTiO₃ increases in the following order: Cr³⁺→Mn⁴⁺→Fe⁵⁺, whereas the CT energies monotonically decrease in the same order.