Next Article in Journal
Approximate Solution of the Thomas–Fermi Equation for Free Positive Ions
Previous Article in Journal
Targeted Cross-Section Calculations for Plasma Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polarizabilities and Rydberg States in the Presence of a Debye Potential

by
Anand K. Bhatia
* and
Richard J. Drachman
NASA/Goddard Space Flight Center, Heliophysics Science Division, Greenbelt, MD 20771, USA
*
Author to whom correspondence should be addressed.
Deceased.
Atoms 2021, 9(4), 86; https://doi.org/10.3390/atoms9040086
Submission received: 11 September 2021 / Revised: 14 October 2021 / Accepted: 16 October 2021 / Published: 22 October 2021

Abstract

:
Polarizabilities and hyperpolarizabilities, α1, β1, γ1, α2, β2, γ2, α3, β3, γ3, δ and ε of hydrogenic systems have been calculated in the presence of a Debye–Huckel potential, using pseudostates for the S, P, D and F states. All of these converge very quickly as the number of terms in the pseudostates is increased and are essentially independent of the nonlinear parameters. All the results are in good agreement with the results obtained for hydrogenic systems obtained by Drachman. The effective potential seen by the outer electron is −α1/x4 + (6β1 − α2)/x6 + higher-order terms, where x is the distance from the outer electron to the nucleus. The exchange and electron–electron correlations are unimportant because the outer electron is far away from the nucleus. This implies that the conventional variational calculations are not necessary. The results agree well with the results of Drachman for the screening parameter equal to zero in the Debye–Huckel potential. We can calculate the energies of Rydberg states by using the polarizabilities and hyperpolarizabilities in the presence of Debye potential seen by the outer electron when the atoms are embedded in a plasma. Most calculations are carried out in the absence of the Debye–Huckel potential. However, it is not possible to carry out experiments when there is a complete absence of plasma at a particular electron temperature and density. The present calculations of polarizabilities and hyperpolarizabilities will provide accurate results for Rydberg states when the measurements for such states are carried out.

1. Introduction

Recently, Qi et al. [1] reported a calculation of the dipole polarizability α1 of the hydrogenic systems in the screened Coulomb potential due to the systems being in hot dense plasmas. The screened potential in this case is the Debye–Huckel potential given by Equation (1).
V(r) = (−Ze2/r) exp(−μr),
where 1/μ = (kBTe/4πe2ne)1/2 is the Debye screening length, kB is the Boltzmann constant, Te and ne are the plasma electron temperature and density, respectively. There are other calculations of energy and polarizability by Zimmermann [2], Bahar et al. [3], Paul and Ho [4], Fowler [5] and Saha et al. [6]. The expression for polarizabilities is given by Equation (1).
Si,k = Σ < 0 │ vi │ N> │ <N│vi │ 0>/(EN − E0)k,
where vi = r1iPi(cosθi) and the various polarizabilities are given by Equation (2).
αi = Si,1, βi = Si,2 and γi = Si,3. ΨN are intermediate states of the appropriate angular momentum and Σ indicates the sum of all such states.

2. Calculations and Results

Qi et al. [1] solved the Schrodinger equation in both the discrete and continuous spectrum of the potential given in Equation (1) by using the symplectic integration scheme. The contribution of the continuum states to the polarizability is particularly important for large μ when the electron binding energy is small and coupling to the continuum is strong. Saha et al. [6] solved the eigenvalue problem by the variational method using a large Slater-type orbitals basis for the discrete states. The contribution of the continuum states was not considered.
In the present calculation, we used pseudostates [7] for helium and the negative hydrogen ion to calculate polarizabilities and hyperpolarizabilties. In [7], the wave functions were of Hylleraas form because the systems consist of two electrons. Definitive results were obtained for the polarizabilities and hyperpolarizabilties. Now we are dealing with single-electron systems and therefore wave functions consist of one electron only and are very simple.
Pseudostates used in this calculation are given by
ΨS = Σi=0Ciri exp(−asr) Y00(Ω),
ΨP = Σi=1Diri exp(−apr) Y10(Ω),
ΨD = Σi=2Eiri exp(−adr) Y20(Ω),
ΨF = Σi=3Firi exp(−afr) Y30(Ω),
where Ci, Di, Ei, and Fi, are the eigenvectors. With a few terms and taking the nonlinear parameters equal to 1.0 for hydrogen atoms, we obtained results for αi, βi, and γi which agree very well for μ = 0 with the known results. We varied all the nonlinear parameters in Equations (3)–(6). However, results were not sensitive to the variation of nonlinear parameters. Therefore, we kept it fixed at 1.0. It is certainly possible to use the hydrogenic functions. The use of pseudostates makes calculations very simple and straightforward. This reduces the efforts involved in the computation. We give in Table 1 the presently calculated results for αi for various values of μ, using 20 terms in the expansion for wave functions (3) and (4). The results are very well converged for shorter expansions as well. We see that the agreement is very good. The present calculation is very easy to carry out compared to that of Qi et al. [1]. They have to calculate each bound state and the wave function corresponding to that state. They also included the continuum of the electron-hydrogenic system. It should be noticed that only those values of μ can be used to calculate α1 for which the ground state of the hydrogen atom remains bound.
In Table 2, we give β1, γ1, α2, β2, γ2, α3, β3, γ3, ε and δ for a few values of μ. The third-order polarizability δ has the following form [8]:
δ = ΣN,M <0 │ vi │ N><N │ vj │ M><M │ <M │ vk │ 0>/(EN − E0)(EM − E0),
The values that (ijk) can take are all the permutations of [1,2].
The fourth-order hyperpolarizability only involving dipole terms has the following form [8]:
ε = ΣN,M,P <0 │ v1 │ N><N │ v1 │ M><M │ v1 │ P │ 0>/(EN –E0)(EM – E0)(EP − E0)
In Equations (7) and (8), potentials v1, v2, and v3 are given reference [8] and are not repeated here. The notation |0>, |N>, |M>. and |P> represent wave functions for the angular momenta S, P, D, and F given in Equations (3)–(6) to be appropriately used to obtain nonzero matrix elements in Equations (7) and (8).
Similarly, by taking the nonlinear parameters equal to 2.0, we obtain results which agree with those obtained by Drachman [7] for μ = 0, using the Dalgarno and Lewis method [9]. The present results are given in Table 3. It can be seen that the polarizabilities and hyperpolarizabilities increase as μ increases. Since the ground state of the He ion is very tightly bound compared to that of the hydrogen atom, it is possible to have much larger values of μ.

3. Rydberg States of He

The long-range potential, in terms of polarizabilities, seen by the outer electron is given by
U(x) = −α1/x4 + (6β1 − α2)/x6 + (δ + 16γ/5)/x 7 + (−α3 + 15β2 − ε − α1β1 − 72γ1 [1 + L(L + 1)/10])/x8,
where x is the distance of the outer electron from the nucleus. In Table 4, we give the expectation value of U(x) in MHz for the wave function of the outer electron in N = 10 and L = 7 and 8 states and compare the present results for μ = 0.0 with those obtained by Drachman [10] without the second-order corrections. We see that the agreement is very good even when very simple wave functions are used.
The expectation values of 1/xn for n = 4, 5, 6, 7, 8, 9 and 10 were calculated for μ = 0.0 and therefore cannot be used when the helium atoms are embedded in the plasma.

4. Transition Rates

Transition rates are given by
A ( n p 1 s ) = 8.032 × 10 9 ( E n p E 1 s ) 3 9 0 R n p ( r ) r 3 R 1 s ( r ) d r
In the above equation, R n p ( r ) and R 1 s ( r ) are the hydrogen functions for np and 1s states. In Table 5, we show how the transition rates from 4p, 3p and 2p states to the 1s state in hydrogen atoms change with the screening parameter. It can be seen that all the rates decrease with the increase in the screening parameter μ. Using the exact wave functions, we find the transition rates for 2p, 3p and 4p states to 1s state equal to 0.626, 0.167 and 0.0682 in units of 109 s−1 when there is no screening and exact hydrogenic functions are used. The last one is not in agreement with the one obtained from pseudostates while the first two are in good agreement with those obtained using pseudostates.

5. Conclusions

We showed that by using pseudostates, it is possible to obtain good results for polarizabilities and hyperpolarizabilties for hydrogenic systems, energies of Rydberg states in the helium atoms and transition rates in hydrogen atoms. The present results obtained using pseudostates for μ = 0 agree with those obtained using elaborate wave functions. Since there is always the presence of a plasma at a particular electron temperature and density, the present results will be useful for comparison with observations. A detailed account of development of the field and possible applications is given in reference [10]. Using the perturbation theory, wave functions, polarizabilities and hyperpolarizabilities given here, the energies of various Rydberg states of any quantum numbers can easily be calculated in the Debye potential.

Author Contributions

A.K.B. and R.J.D. contributed equally. All authors have read and agreed to the published version of the manuscript.

Funding

No funding was received for this research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Qi, Y.Y.; Wang, J.G.; Janev, R.K. Bound-bound transitions in hydrogenlike ions in Debye plasmas. Phys. Rev. A 2009, 80, 032502. [Google Scholar] [CrossRef]
  2. Zimmermann, R. The Green’s function of the Debye potential: Evaluation of the ground-state polarizability. J. Phys. B 1985, 18, 2817. [Google Scholar] [CrossRef]
  3. Bahar, M.K.; Soylu, A.; Poszwa, A. The Hulthen Potential Model for Hydrogen Atoms in Debye Plasma. IEEE Trans. Plasma Sci. 2016, 10, 2297. [Google Scholar] [CrossRef]
  4. Paul, S.; Ho, Y.K. Electric Field on Hydrogen Atom. Phys. Plasmas 2010, 17, 082704. [Google Scholar] [CrossRef]
  5. Fowler, P.W. Polarizability of the Electron in a Debye Potential. J. Phys. B 1986, 19, L1. [Google Scholar] [CrossRef]
  6. Saha, B.; Mukherjee, P.K.; Diercksen, G.H.F. Energy levels and structural properties of compressed hydrogen atoms under Debye screening. Astron. Astrophys. 2002, 396, 337. [Google Scholar] [CrossRef]
  7. Bhatia, A.K.; Drachman, R.J. Polarizabilities of helium and negative hydrogen ion. J. Phys. B Mol. Opt. Phys. 1994, 27, 1299. [Google Scholar] [CrossRef]
  8. Drachman, R.J.; Bhatia, A.K. Rydberg levels of lithium. Phys. Rev. A 1995, 51, 2926. [Google Scholar] [CrossRef]
  9. Dalgarno, A.; Lewis, J.T. The exact calculation of long-range forces between atoms by perturbation theory. Proc. R. Soc. Lond. Ser. A 1955, 233, 70. [Google Scholar]
  10. Drachman, R.J. Long-Range Casimir Forces, Theory and Recent Experiments on Atomic Systems; Levin, F.S., Micha, D.A., Eds.; Plenum Press: New York, NY, USA, 1993; p. 219. [Google Scholar]
Table 1. Comparison of the presently calculated results for αi for hydrogen atom for various values of μ, using 20 terms in the wave functions.
Table 1. Comparison of the presently calculated results for αi for hydrogen atom for various values of μ, using 20 terms in the wave functions.
MPresentQi et al. [1]
0.004.5004.500
0.014.50220A
0.024.508684.50820
0.0254.513464.51299
0.0504.552204.55176
0.06254.580494.58003
0.104.699784.69933
0.205.276375.27661
A: No results are given for this value of μ by the authors.
Table 2. Polarizabilities and hyperpolarizabilities for hydrogen atoms for various values of μ.
Table 2. Polarizabilities and hyperpolarizabilities for hydrogen atoms for various values of μ.
μ = 0.00.020.0250.050.06250.100.20
β15.37505.39175.40095.47515.53055.76366.9384
γ16.64586.67506.69116.82196.91827.33139.5128
α215.00015.054015.083815.324615.500816.249820.0479
β213.375013.445313.484013.797714.028415.019520.3300
γ212.494812.582912.631513.026213.317914.587221.8128
α3131.2500131.989132.396135.6875138.1049148.490204.341
β3102.031102.817103.249106.758109.349120.632185.434
γ383.204484.032084.497488.198090.9547103.142178.390
Ε59.212559.711059.930761.717263.037168.7756101.226
Δ106.500107.036107.331109.724111.484119.044159.391
Table 3. Polarizabilities and hyperpolarizabilities for the helium ion for various values of μ.
Table 3. Polarizabilities and hyperpolarizabilities for the helium ion for various values of μ.
μ = 0.00.020.0250.050.06250.100.20
α10.28125 0.281390.281460.282090.282550.284510.29374
β10.083940.084060.084090.084390.084610.085560.09006
γ10.025960.025990.026010.026140.026230.026650.02864
α20.234380.234590.023470.235680.236400.239450.25390
β20.052250.052320.052350.052670.052910.053890.58670
γ20.012200.012220.012240.012340.012410.012720.14245
α30.512690.513430.513840.517170.519620.530030.58004
β30.099640.099840.099940.100830.101480.104260.11780
γ30.020310.020370.020390.020630.020800.021530.02518
ε0.057920.058020.058080.058530.058860.060270.06716
δ0.416020.416550.416840.419260.421040.428610.46501
Table 4. Comparison of the presently calculated Rydberg states with the results of Drachman.
Table 4. Comparison of the presently calculated Rydberg states with the results of Drachman.
NLPresent
(MHz)
Drachman [6]
(MHz)
107−48.60605124−48.60604738
108−24.17853458−24.17853458
Table 5. Transition rates in sec (−1) 4p, 3p and 2p states of hydrogen atoms to the 1s state. The rates have been multiplied by 10(−9).
Table 5. Transition rates in sec (−1) 4p, 3p and 2p states of hydrogen atoms to the 1s state. The rates have been multiplied by 10(−9).
μ A (2p→1s) A (3p→1s)A (4p→1s)
0.0000.6242930.1666700.071760
0.0200.6183430.1600130.066440
0.0250.6151370.1566130.064150
0.0500.5901800.1318770.053783
0.06250.5726930.1156470.051620
0.10000.5032900.0616070.058523
0.20000.2156530.0349230.070210
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bhatia, A.K.; Drachman, R.J. Polarizabilities and Rydberg States in the Presence of a Debye Potential. Atoms 2021, 9, 86. https://doi.org/10.3390/atoms9040086

AMA Style

Bhatia AK, Drachman RJ. Polarizabilities and Rydberg States in the Presence of a Debye Potential. Atoms. 2021; 9(4):86. https://doi.org/10.3390/atoms9040086

Chicago/Turabian Style

Bhatia, Anand K., and Richard J. Drachman. 2021. "Polarizabilities and Rydberg States in the Presence of a Debye Potential" Atoms 9, no. 4: 86. https://doi.org/10.3390/atoms9040086

APA Style

Bhatia, A. K., & Drachman, R. J. (2021). Polarizabilities and Rydberg States in the Presence of a Debye Potential. Atoms, 9(4), 86. https://doi.org/10.3390/atoms9040086

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop