If one admits that the linear response to an external electric field E is: the energy associated with this external field is then:
that is: .
This gives the correct expression for the polarization potential, a short range attractive potential deeper than the Coulomb potential, for short radial distance r.
Using experimental spectroscopy and neglecting fine structure splittings, it is possible to measure the wavelengths of the transitions of these heavier atoms. These transitions depart from those of the non-relativistic of the hydrogen atom:
The regular spectrum of hydrogen is modified to deal with the experimental alkaline atoms spectra, such as Na, Li, or Mg, with an analytical approach introducing the quantum defects: .
The modification for an electron jump is given by:
It is possible to model the atoms with the so-called “optical electron” scheme, meaning that far from the cloud of the remaining electrons, appears a dipolar force. The compound cloud with the Z protons of the nucleus and the remaining electrons interacts with the outer “optical” electron, acts as a dipolar interaction, thus repelling the compound cloud.
Using quantum defect theory, the quantum defects can be calculated two fold if the transition in an alkaline atom is observed experimentally two equations, suffice to specify the values of and .
Thus for an transition,
The theory of these defects is explained when one introduces the polarization potential, which acts in the same sign as the Coulomb potential. Its expression is: atomic units and is the static dipole polarizability. The polarization effect on the outer electron, assuming nearly hydrogenic states, can be expressed as:
Using known results  of the expectation value of: , the correction to the Coulomb potential is:
Solving this equation using a limited development, thus assuming that , the theorical quantum defect is fixed. Once the theorical is set, a simple code furnishes the numerical quantum defects associated with a spectroscopically identified transition wavelength of an alkaline atom. The comparison of these two values enables us an estimate the static polarizability for the considered atom species.
2. Experimental Quantum Defects for MgI and Ions MgII MgIII, NaI, NaII, NaII
For alkaline atoms, such as Na, Ca, Li, Mg, K, and their ions, it is possible to measure the wavelengths associated with identified transitions. These can be found in tables Moore (1949) . It is easy to use symbolic software like Mathematica, and solve two equations, if one wants to have two quantum defects: and , that requires two identified transitions of any alkaline, for instance for the Na (sodium) atom and Mg (magnesium). These data are listed in Table 1 and Table 2.
Using Mathematica syntax gives the way to proceed with the experimental determination of quantum defects and and is quite simple. Neglecting
Table 1. Values of NaI quantum defects l = 0, 1, 2.
Table 2. Values of MgI quantum defects l = 0, 1.
fine structure, one can find in tables of data , two transitions, observed in experimental spectroscopy such as for MgI that is neutral Magnesium:
MgI and that is in atomic units:
and assuming that the quantum defects depend on the value of the momenta and that is an transition.
One can perform the calculation of the two quantum defects by using the Solve function common in many symbolic softwares (Mathematica, Mapple, Matlab, Python). It is necessary to use two independent equations for two distinct quantum defects.
The two variables are: . A lot of solutions appear, it is easy two isolate the good set of variables.
The numerical results are the following:
The complete set of solutions contains imaginary terms for the defects: thus ignored, is out of the range. It remains the correct and .
This solve function can be used for any couple of identified transitions, for example charged ions: MgII, or MgIII MgIV, that is Mg+, Z = 1, Mg++, Z = 2, Mg+++, Z = 3.
Using the first order of the development, requiring a small quantum defect (a D state), and setting atomic units and .
Some simple Mathematica instructions, reproduce quickly the relation written in Equation (10) .
The results obtained with a small code are in accordance with Equation (10) in :
Jungen formula (16)
solution is: (21)
my solution (22)
Jungen value (23)
Assuming a small quantum defect: and for Na+ atom, this equation enables us to give an estimate of the static polarizability .
3. Some Remarks on the Laguerre Polynomial Mathematica Function LaguerreL
Table 3 below shows how the quantum defects vary for an atom like Magnesium and its ions. The greater is, the greater the quantum defects are. These data are taken from the Topbase on line database.
To deal with quantum wave functions, one uses the following form called the generalized Laguerre polynomial, which satisfies the differential equation:
The equation below shows the hydrogenic wave functions; it is given because it serves to show how the quantum defect wave functions are modified from .
The LaguerreL polynomial is written as given in the Mathematica software
Table 3. Values of Mg ions quantum defects l = 0, 1 as function of the .
a# is the charge that moves the optical electron. The building of the wave function for low quantum numbers: breaks down. A relativistic theory is needed and is an open research problem.
where the serves to represent the hydrogen wave functions. If one deals with an hydrogenic ion with a nuclear charge Z the radial r has to be changed to . When solved with correct physical parameters using atomic units: , this equation for the motion of a bound electron gives the discrete spectrum for the bound states:
The total energy H is and the full solution for the radial part is obtained using spherical coordinates . This is called the Schrodinger wave equation and the bound states are obtained by solving:
is easily solved, using two Mathematica functions. The Solve function, when correctly written, gives the correct solution for the attractive Coulomb potential
with the centrifugal term .
In the complete Schrodinger equation, the angles and separate, from the variable r, it is well known that the eigen function for the angles is the spherical harmonics: .
4. Building of the Wave Function with a Principal Quantum Number
The Case of the Alkaline Optical Electron Radial Wave Functions
For what concerns the construction of a wave function using the non integer principal quantum number, defined as: , we need to change the arguments, of the well built LaguerreL polynomial. In any version of the Mathematica software, the l momentum has to be changed to together with for the projection of the momentum . With this software, it possible is to obtain a correct result, for non integer arguments: for instance the
gives good numbers.
I will simply change the arguments of the Laguerre polynomial:
One may notice that some arguments do not change with the transformation that is .
The factorial function with Mathematica, works for any number: and the Laguerre polynomial: that is the generalized function n,l, fractional arguments.
If one deals with an ionic core of charge Z, one has to change the r radial variable into in the .
For what concerns the calculation of the quantum defect wave function it is possible to extend continuously the transformation to .
The author shows how the quantum theory, with the use of wave functions taking into account the modification of the quantum numbers, as explained by Kostelecky and Nieto  provides a simple framework for computing atomic and ionic wave functions using current symbolic software.
This is a consequence of the extreme refinement of symbolic softwares now available: Mathematica, Mapple, MatLab, Python. The old quantum theory exposed by Max Born  explains that the effect of the polarization potential results in a precession of the elliptical trajectory of the electron around its perihelion. This work shows how to simply produce, with the help of a modified Schrodinger equation, the probability distribution for the optical electron near the nucleus of alkaline atoms and their ions. Our approach to this problem is self-consistent that is: a simple code is given, this enables to calculate the quantum defects associated with a spectroscopically identified transition of an alkaline atom such as MgI Equation (2), then another code is used to construct the quantum wave functions using these set of quantum defects Equation (29), as shown in the figures.
Our approach is coherent: it conciliates the calculation of the quantum defects coming from spectroscopical data, for which exist databases, with the building of the wave functions. These change from those of hydrogen with the the set of quantum numbers. The change means that the set quantum numbers used are . It is the simple and easy way to determine the static polarizability of many atoms, these quantities are still an open research problem, beyond the scope of this paper.
Figure 1 shows the hydrogen wave function , (a S state), and two alkaline wave functions and .
Figure 2 concerns P state, and three wave functions (with different ).
polarization potential effect dominating the Coulomb potential .
Figure 1. Red curve is the n = 4 l = 0 , green curve corresponds to , with , the scarlet curve corresponds to , with , this curve shows clearly that the probability for the electron to be near the nucleus is greater than the hydrogen electron.
Figure 2. Red curve is the n = 4 l = 1 hydrogen wave function , green, blue and scarlet curves correspond to the wave function , , r) with , the higher is , the higher is the contribution of the alkaline wave function near the nucleus, here obtained near the origin r = 0.
Finally it is possible to build the difference function, here called: and to give a correct illustration of the effect inherent to the polarization potential . These are the main purposes of this article.
6. Comments on the Figures
Figure 3 contains S and P states representations of the hydrogen wave functions
Figure 3. This shows how the introduction of the quantum defects in the wave function differs from the hydrogen wave functions red curves. The first plot and the second represent respectively a S state and a P state for the hydrogen electron (red curve) and the other color curves describe the model alkaline electron. The third plot superposes all these.
and those of alkaline wave functions modified by the introduction of the quantum defects .
Figure 4 represents the square of the wave functions and , these squared wave functions are the true probabilities for the electron of the hydrogen atom, and the optical electron to be near the nucleus, a proton or an ionic core with .
Figure 5 shows that the optical electron with its squared wave function , with and gets nearer to the nucleus than the electron of the hydrogen atom .
This is a consequencce of the attractive effect of the polarisation potential. It is not difficult to write a small code, taking into account the modification of the wave functions, using non integer and , the limits of validity of the modified wave function:
Figure 4. The red curve is with the blue curve with . These curves are everywhere positive, showing the real probabilty functions for the electron of the hydrogen atom, and the optical electron of a model alkaline atom.
Figure 5. The difference function: is plotted, that is n = 3 S l = 0 and P l = 1 states of the hydrogen squared wave functions subtracted from the squared wave functions modified by different quantum defects ( ). The greater the quantum defects are, the higher are the contributions of the polarization potential near the origin r = 0, that is near the core nucleus.
Figure 6. and . These plots illustrate and isolate the quantum effect of the polarization potential once subtracted the Coulomb potential for hydrogen.
are and , this is well explained by two authors Kostelecky and Nieto with their factor: where the integer factor or , this transformation enables to obtain a continuous wave function: even when one deals with low quantum numbers n and high quantum defect . Figure 6 illustrates the quantum strength of the polarization potential alone, it can be visualized in the old quantum theory as a precession motion added to the elliptical trajectory described by the outer electron.
The author thanks Dr Régis Courtin, for improving the english of this paper, and for his carefull reading of the text.
1At this stage, when the quantum defect is defined one can equalize the quantity: . This exactly represents the difference between the hydrogen energy levels from the deeper energy levels associated with the effect of the polarization potential. This leads to an estimate of the static polarizability by solving the following set of equations. Dealing with charged ions one has to change the factor .