Variational method for ground-state energy of helium atom in N dimensions

IL NUOVO CIMENTO Vol. 25 B, N. 9 Settembre 2 DOI.393/ncb/i2-97-2 Basic topics: Fundamental Physics Variational method for ground-state energy of helium atom in N dimensions Sami M. AL-Jaber( ) Department of Physics, An-Najah National University - P.O. Box 7, Nablus, Palestine (ricevuto il 22 Maggio 2; approvato il 29 Giugno 2; pubblicato online il 6 Settembre 2) Summary. The ground-state energy of the N-dimensional helium atom is presented by applying the variational principle. The calculations are made for the unscreened and screened cases. It is shown that, in both cases, the magnitude of the ground-state energy decreases (less negative) as the spatial dimension N increases. For the unscreened case, the relative contribution of the electron-electron interaction term to the ground-state energy is calculated for different dimensions, and it is found that this ratio approaches one half as N. For the screened case, the effective nuclear charge is computed for different dimensions and its limiting value is found to be 3/2 in the infinite-dimensional space. In addition, the relative contribution of screening to ground-state energy is calculated in different dimensions and it is shown that it reaches /8 as N. Furthermore, the asymptotic behavior of ground-state energy, contribution of electron-electron interaction and contribution of screening effect are presented. PACS 3.65.-w Quantum mechanics. PACS 3.65.Ca Formalism. PACS 3.65.Ta Foundations of quantum mechanics; measurement theory.. Introduction The helium ground-state energy has been a subject of interest for a long time [-7] and continues to attract the attention of many investigators [8-]. The variational technique is considered as a powerful tool for estimating the ground-state energy of a quantum system. Over the last few years, this method has been applied to helium atom in different settings: Komasa and Rychlewski [2] considered correlated Gaussian functions in variational calculations to estimate the ground-state energy of helium dimmer, Banerjee [3] and Flores-Riveros et al. [4] studied energy spectrum of spherically confined helium atom, Theodorakis et al. [5] considered optimization method to find ( ) E-mail: jaber@najah.edu; http:www.najah.edu c Società Italiana di Fisica 99

SAMI M. AL-JABER approximate solutions of field equations, Lobanova et al. [6] and Biaye et al. [7] used variational method in calculating energies of excited states of helium to illustrate the importance of electron correlation effects. Hydrogen-like systems are considered a prototype of physical problems in N dimensions [8-2] and have an important role in different areas of physics: For example, in quantum field theory [2], in quantum chemistry [22], in quantum computation [23] and in nanotechnology [24]. In addition, the generalization of physical problems to higher space dimensions is useful in random walks [25], in Casimir effect [26], in harmonic oscillators [27-29] and in statistical mechanics [3, 3]. Furthermore, the N-dimensional radial Schrödinger equation has been examined for different kind of potentials [32-36]. In the present paper, we will derive an analytic expression for the ground-state energy for helium atom in N dimensions via the variational method. We do this for two cases: In the first case, we assume the nuclear charge to be unshielded by the electrons, and in the second case we consider screening of the nuclear charge by the electrons. The paper is organized as follows: in sect. 2, we present some aspects of the N-dimensional hydrogenic system which are relevant to our work. In sect. 3, we calculate the ground-state energy of helium atom in N dimensions by applying the variational method assuming no screening of the nuclear charge. In sect. 4, we include the effect of screening on the ground-state energy. Section 5 is devoted for conclusions. 2. The helium atom in N dimensions: preliminaries We consider a helium atom which consists of two electrons and a nucleus with a charge Ze (with Z = 2). The Hamiltonian for this system (assuming infinite nuclear mass and ignoring fine structure and small corrections) is given by () H = H + V ee, where H is the sum of two Hamiltonians for one-electron hydrogenic-like atoms, given by (2) H = h2 ( 2 2m + 2 e 2) 2 ( Z + ), 4πε r 2 and V ee represents the electron-electron interaction which is given by (3) V ee = e2 4πε 2, with 2 =/ r 2.IntheN-dimensional space, the Laplacian operator in the polar coordinates (r, θ,θ 2,...,θ,ϕ)ofR N has the form (4) 2 = 2 (N ) + r2 r r + r 2 Ω2, where Ω 2 is a partial differential operator on the unit hyper-sphere R N and whose eigenvalues are l(l + N 2). The angular solutions are the hyper-spherical harmonics Y {m} l which depend on (N ) angular coordinates with ranges θ j π and ϕ 2π. Each of these spherical harmonics is determined by (N ) integers l, m,m 2,...,m that can assume all values consistent with the inequalities

VARIATIONAL METHOD FOR GROUND-STATE ENERGY OF HELIUM ATOM ETC. l m m 2... m. The eigenfunctions of H are the product of the eigenfunctions of the one-electron Schrödinger equation for and r 2 coordinates, namely (5) Φ () (, r 2 )=ψ( )ψ( r 2 ), where (6) ψ( r) =AR Nnl (r)y m l (θ j,ϕ). j =, 2, 3,...,N 2. The functions R Nnl (r) are the solutions for the radial Schrödinger equation and are given by [37] (7) R Nnl (r) =e ρ/2 ρ l L (2l+) (n+l+n 3) (ρ), where L(ρ) is the associated Laguerre polynomial and ρ =(8m( E)/ h 2 ) /2. The energy eigenvalues of H are (8) EnN =(EnN ) +(EnN 2Z 2 ) 2 = ( ) n + N 3 2 E, 2 where E is the three-dimensional ground-state energy of a one hydrogen atom that is given by (9) E = 2a(4πε ), whose value equals 3.6eV and a =4πε h 2/me 2 is the Bohr radius. Since we are interested in the ground-state energy of the helium atom in N dimensions, eq. (8) gives the ground-state energy of H by putting Z =2andn =, with the result () E N = e 2 32 (N ) 2 E. 3. Variational method without screening In order to find the contribution of the electron-electron interaction term V ee to the ground-state energy, we apply the variational principle by adopting the trial wave function () φ(,r 2 )=A exp [ ] 2Z(r + r 2 ), a(n ) where A is a normalization constant. Normalization requirement of the above trial wave function over the two sets of the N-dimensional electron coordinates gives (2) = A 2 [ exp 4Z(r ] + r 2 ) r N r2 N d dr 2. a(n )

2 SAMI M. AL-JABER The integration over the r coordinates can be easily carried out by letting x =4Zr/a(N ) and using (3) x N e x dx =Γ(N), which immediately gives the value of the constant A, namely (4) A = [ ] N 4Z. Γ(N) a(n ) We apply the variational principle, using the trial wave function given in eq. (), to find the contribution of V ee to the ground-state energy of the helium atom. This is given by (5) V ee = e2 A 2 4πε It is helpful to expand / r 2 as [38] ( r 2 exp 4Z(r ) + r 2 ) r N r2 N d dr 2 dω dω 2. a(n ) (6) r 2 = l= r< l r> l+ C ()/2 l (cos θ), where θ is the angle between and r 2 and C ()/2 l are the Gegenbauer polynomials. The addition theorem for the M linearly independent hyperspherical harmonics Yl m of degree l on the S N sphere has the form [39] (7) C ()/2 l (cos θ) = M N 2 (2l + N 2) Ω(N) m= Yl m (Ω )Yl m (Ω 2 ), where Ω(N) =2π N/2 /Γ(N/2) is the surface area of the N-dimensional unit sphere, and M =(2l + N 2)(l + N 3)!/(N 3)!. Therefore eq. (5) becomes (8) V ee = e2 4πε A 2 exp = e2 4πε A 2 = e2 4πε A 2 l= M Ω(N)(N 2) (2l + N 2) m [ 4Z( + r 2 ) a(n ) d Yl m (Ω )Yl m (Ω 2 )dω dω 2 ] 2 d dr 2 2 exp e λr d [ r [ 4Z( + r 2 ) a(n ) 2 dr 2 e λr2 + ] r > r 2 dr 2 e λr2 ], where we put λ = 4Z/a(N ) and used the orthonormality of hyperspherical harmonics [39] (9) Y m l (Ω)Y m l (Ω)dΩ = δ ll δ mm.

VARIATIONAL METHOD FOR GROUND-STATE ENERGY OF HELIUM ATOM ETC. 3 Using tables of integrals [4], (2) gives us (2) (22) and r 2 dr 2 e λr2 = x n e ax dx = e ax e λr N j= r2 dr 2 e λr2 = e λr n j= ( ) j x n j a j+ n! (n j)!, ( ) j (N )! λ j+ (N j)! rn j ( )N (N )! λ N, j= ( ) j (N 2)! λ j+ (N 2 j)! r j, (23) r (N ) e λr d = ( )N (N )! λ N. Therefore, after some algebra, we get (24) r 2 e λr2 dr 2 + r 2 e λr2 dr 2 = ( )N (N )! λ N (e λr ) +e λr j=o ( ) j j(n 2)! λ j+ (N j)! r j. The integration over, given in eq. (8), is now straightforward with the use of eq. (23) and the result is [ ( V ee = e2 A 2 ( )2N 3 (N 2)! 4πε λ 2N (N )! ) (25) 2 N ] j(2n 3 j)! 2 2 j, (N j)! j= j= which upon the substitution of A and λ becomes ( ) [ ( V ee = e2 4Z 4πε a (N ) 3 (N )! ) (26) (N 2)! 2 N ] j(2n 3 j)! 2 2 j. (N j)! The use of eq. (9) and after little manipulations, the above equation takes the form (27) V ee = 8E Z ( (N ) 2 2 N j(2n 3 j)! ) 2 N (N )! 2 j. (N j)! j=

4 SAMI M. AL-JABER Table I. Ground-state energies in different space dimensions. N EN V ee E N V ee /EN 3 8.8 34 74.8.325 4 48.3555 6.6222 3.7333.34374 5 27.2 9.8825 7.387.36328 6 7.48 6.562.846.37695 7 2.889 4.689 7.48.3872 8 8.8863 3.558 5.374.39526 9 6.8 2.7323 4.6769.48 5.37284 2.8868 3.8467.47264 4.352.79259 2.5594.49 2 3.5967.49588 2.8.459 3 3.222.26755.7546.494 4 2.575.882.487.4225 5 2.224.94428.276.42527 Our result clearly shows the dependence of the contribution of the electron-electron interaction on the space dimension N. Therefore, using eq. (), the ground-state energy, which corresponds to the Hamiltonian H, in N-dimensional space is (28) E N = E N + V ee = 32 (N ) 2 E + V ee. For the three-dimensional case (N = 3), eq. (27) gives V ee = 5E Z/4 and thus E N =2Z 2 E 5E Z/4, which is an expected result [4]. It is constructive and illuminating to compute EN, V ee and E N for different values of N. We show these results in table I, where we used E = 3.6eV and Z = 2. It is interesting to calculate the relative contribution of the electron-electron interaction to the ground-state energy, namely V ee /EN, which we include in the last column of table I. For comparison purposes, we include in the second column of table I the ground-state energy without V ee, which decreases as the space dimension N increases. This is due to the extra repulsive term in the effective potential which arises in dimensions greater than three, whose role has been emphasized by the present author more than a decade ago [37]. Furthermore, the results in the third column show that the contribution of the electron-electron interaction decreases as the space dimension increases. This is so since, as N increases, the repulsive extra term in the effective Hamiltonian tries to repel the electron away from the nucleus and hence the average of /r decreases. The fourth column gives the ground-state energies including the contribution V ee. The last column is interesting: It shows that the relative contribution of the electron-electron interaction, V ee /EN, increases as the space dimension increases. In fact, the terms in the curly bracket in eq. (27) goes to unity as N and thus, eq. (27), with Z = 2, gives (29) V ee = 6E /(N ) 2.

VARIATIONAL METHOD FOR GROUND-STATE ENERGY OF HELIUM ATOM ETC. 5 Therefore, with the use of eq. (), the relative contribution of V ee for large N is (3) lim N V ee /E N =/2. One must notice that though each of V ee and EN limiting ratio equals exactly one half. 4. Variational method with screening vanishes as N, yet their In this section, we consider the effect of electron screening of the nuclear charge. On the average, each electron represents a cloud of negative charge which partially shields the nucleus, so that the second electron sees an effective nuclear charge, Z. Therefore, we set Z a variational parameter in our variational calculations for the ground-state energies in different dimensions. We may write the Hamiltonian given in eq. (2) as (3) H = h2 2m ( 2 + 2 2) e 2 4πε ( Z + Z r 2 ) + e2 4πε The expectation value of the above Hamiltonian is clearly ( ) e (32) H =2Z 2 EN 2 +2(Z 2) + V ee, 4πε r ( Z 2 + Z 2 ) +. r 2 r 2 where EN =4E /(N ) 2 is the ground-state energy of a one hydrogen atom in N dimensions. The expectation value of /r in the N-dimensional one-particle hydrogenic atom with nuclear charge Z is given by [37] (33) = r which gives for the ground state (n =) Z a ( n + N 3 2 ) 2, (34) e 2 4πε = e2 r 4πε a 4Z (N ) 2 = 2ZE N. V ee is given by eq. (27) which can be written as (35) V ee = 2ZEN ( 2 N j(2n 3 j)! ) 2 N (N )! 2 j. (N j)! j= The substitution of eqs. (34) and (35) into eq. (32) yields (36) ( [ H = 2Z 2 +8Z 2Z ( 2 N ) 2 N (N )! j= ]) j(2n 3 j)! 2 j EN. (N j)!

6 SAMI M. AL-JABER Table II. Effective charge Z, ground-state energy with screening, H and its deviation and percentage from that without screening. N Z H H E N H E N E N % 3.6875 78.29 3.49 4.665 4.65625 33.68.4285 4.5 5.63672 8.267.89747 5.82 6.623.46373.677 5.695 7.6279 7.86.453 6.6 8.6474 5.7795.3469 6.459 9.5989 4.3426.2744 6.747.59273 3.4746.2228 6.996.5889 2.74395.8459 7.225 2.5849 2.25634.555 7.44 3.5859.88758.329 7.578 4.57749.625.49 7.7297 5.57472.37652.4 7.8687 We minimize H with respect to Z by setting Z H = which immediately gives (37) Z =2 ( 2 N j(2n 3 j)! ) 2 2 N (N )! 2 j. (N j)! j= Clearly, the above equation shows that the effective nuclear charge which is seen by the electron depends on the space dimension N, and for the three-dimensional space (N = 3) it yields an effective charge Z =27/6 =.6875. With the help of eq. (37), it is tempting to calculate the effective nuclear charge Z, which minimizes H, in different space dimensions. Then eq. (36) can be used to calculate the minimum of H, which is just the ground-state energy that includes the screening effect. We also calculate the contribution of screening, H E N, and its percentage, H EN E N %, to the ground-state energy. These results are presented in table II. The results in table II show that as the space dimension N increases, the effective charge Z decreases (i.e. screening increases) and the corresponding ground-state energy, H, decreases in magnitude (becomes less bound). In addition, comparison of column three in table II with that in table I shows that the ground-state energy with screening is more negative than the ground-state energy, E N, without screening. Furthermore, it is worth noticing that the ratio H EN E N increases as N increases. It is interesting to examine the very large N limit: As mentioned earlier, the curly bracket in eq. (27) goes to unity and thus eq. (37) yields an effective charge Z =3/2, which upon its substitution into eq. (36), we get H =8E /(N ) 2. As was shown in the previous section, the asymptotic behavior of the ground-state energy is E N =6E /(N ) 2 and therefore the ratio H EN E N =/8 asn, which represents the relative contribution of screening to the ground-state energy of the helium atom in the infinite-dimensional space.

VARIATIONAL METHOD FOR GROUND-STATE ENERGY OF HELIUM ATOM ETC. 7 5. Conclusions The variational method has been used to estimate the ground-state energy of an N-dimensional helium atom. This has been done for two cases: the first without screening and the second with screening. In the first case, the contribution of the electron-electron interaction term, V ee, to the ground-state energy, was shown to depend on the space dimension. Thus the ground-state energy decreases in magnitude (less negative) as N increases and it asymptotically behaves as 6E /N 2. Our results also show that the relative contribution of V ee to the ground-state energy E N increases as N increases and it approaches a value of one half in the infinite-dimensional space (N ) despite the vanishing limiting value of each of V ee and E N in that limit. In the second case, the effective nuclear charge, due to electron screening, and its effect on the ground-state energy, were calculated in different dimensions. The results show that as N increases, the effective nuclear charge Z decreases (i.e. screening increases) and it approaches the value 3/2 as N, and the ground-state energy H decreases in magnitude and it behaves asymptotically as 8E /N 2. It was also shown that the relative contribution of screening increases with increasing N and it approaches the value /8 as N. I would like to thank An-Najah National University for granting me a one year sabbatical leave, and I also would like to thank Palestine Technical University-Kaddorie (PTUK) for the kind hospitality I received during this leave. REFERENCES [] Hylleraas E. and Midtdal J., Phys. Rev., 3 (956) 829. [2] Perkeris C. L., Phys. Rev., 2 (958) 649. [3] Perkeris C. L., Phys. Rev., 5 (959) 26. [4] Schwartz C., Phys. Rev., 28 (962) 46. [5] Kinoshita T., Phys. Rev., 5 (959) 366. [6] Rau A. R. P., Phys. Rev. A, 4 (97) 27. [7] Leopold J. G. and Percival I. C., J. Phys. B: At. Mol. Phys., 3 (98) 37. [8] Korobov V. I., Phys. Rev. A, 66 (22) 245. [9] Korobov V. I., Phys. Rev. A, 6 (2) 6453. [] Sims J. S. and Hagstrom S. A., Phys. Rev. A, 8 (29) 5257. [] Sims J. S. and Hagstrom S. A., J. Chem. Phys, 24 (26) 7. [2] Komasa J. and Rychlewski J., Mol. Phys., 9 (997) 99. [3] Banerjee A., Kamal C. and Chowdhury A., Phys. Lett. A, 35 (26) 2. [4] Flores-Riveros A., Aquino N. and Montgomery H. E., Phys. Lett. A, 374 (2) 246. [5] Theodorakis S. and Kindyni N., Am. J. Phys., 78 (2) 244. [6] Lobanova E. M., Sheinerman S. A. and Gerchikov L. J., J. Exp. Theor. Phys., 5 (27) 486. [7] Biaye M., Knoté A., Faye N. A. B. and WaquéA., Eur. Phys. J. D, 3 (2) 2. [8] Herschbach D., Avery J. and Goscinski O., Dimensional Scaling in Chemical Physics (Kluwer, Dordrecht) 993. [9] Avery J., Hyperspherical Harmonics and Generalized Sturmians (Kluwer, Dordrecht) 2. [2] Avery J. and Herschbach D., Int. J. Quantum Chem., 4 (992) 673. [2] Witten E., Phys. Today, 33 (98) 38. [22] Aquilanti V., Cavalli S. and Colleti C., Chem. Phys., 24 (997).

8 SAMI M. AL-JABER [23] Dylman M. I., Paltzman P. M. and Seddigard P. M. P., Phys. Rev. B, 67 (23) 5542. [24] Li S. S. and Xia J. B., Phys. Lett. A, 366 (27) 2. [25] Mackay T. D., Barlett S. D., Stephensen L. T. and Sanders B. C., J. Phys. A: Math. Gen., 35 (22) 2745. [26] Bender C. M., Boethtcher S. and Lipatov L., Phys. Rev. D, 46 (992) 5557. [27] Oyewumi K. J., Akimpilu F. O. and Ayboola A. D., Int. J. Theor. Phys., 47 (28) 39. [28] Al-Jaber S. M., Int. J. Theor. Phys., 47 (28) 853. [29] Rothos V. M., Nistazkis H. E., Kevyebidis P. G. and Frantzeskakis D. J., J. Phys. A: Math. Theor., 42 (29) 2527. [3] Al-Jaber S. M., Int. J. Theor. Phys., 38 (999) 98. [3] Al-Jaber S. M., An-Najah J. Res., 22 (28) 67. [32] Khan G. R., Eur. Phys. J. D, 53 (29) 23. [33] Calogero F., Eur. Phys. J. D, 57 (2). [34] Al-Jaber S. M., Nuovo Cimento B, 2 (997) 76. [35] Land M., Found. Phys., 37 (27) 455. [36] Chen G., Z. Naturforsch. A, 59 (24) 875. [37] Al-Jaber S. M., Int. J. Theor. Phys., 37 (998) 289. [38] Mehta M. L. and Normand J. M., J. Phys. A, 3 (997) 867. [39] Crosche C. and Steiner F., J. Math. Phys., 36 (995) 2354. [4] Dwight H. B., Tables of Integrals and other Mathematical Data (Macmillan Publishing Co, New York) 96. [4] Griffiths D. J., Introduction to Quantum Mechanics (Prentice Hall, New Jersey) 995.