New Bohr model calculates Helium ground state energy correctly

Y. Tsubono
303 5-16 Daita6-chome, Setagaya-ku Tokyo 155-0033, Japan
URL: http://www7b.biglobe.ne.jp/~kcy05t/index.html

PDF version

Abstract :

The quantum mechanical methods can calculate the behavior of atoms and molecules almost correctly using the methods such as perturvation, variation, and density functional theories. But these methods can not express concretely the motion of the electron such as spin and circular movement. So in this study, we try to describe the motion of the two electrons of the helium atom using the classical methods based on the Bohr model. and using this new model, we could theoretically calculate the ground state energy of the helium, which is consistent with the experimental value (-79.005eV). So it is possible that this method would be a new tool for calculating the electronic behavior.

Introduction :

In 1913, by the study of the emission spectroscopy of atomic hydrogen, Niels Bohr postulates the Bohr model by using classical mechanics and electromagnetic theory. In the Bohr model, it was hypothesized that an electron is circularly moving around the nucleus, and the orbital angular momentum is quantized in units of hbar= h / 2 π (h is Plank's constant = 6.62606896 x 10^(-34) Js). So we have

(1-1)

where m and v are the mass (= 9.1093826 x 10^(-31) kg) and the velocity (m/sec, this "m" means meter, not mass ) of the electron, and n is a positive integer. Equating the centrifugal force to the Coulomb force, we have

(1-2)

where r is the circular orbital radius (m), e is the electron charge (= 1.60217653 x 10^(-19) C), and ε is the permittivity of vacuum (= 8.854187817 x 10^(-12) C^(2) / Nm^(2)). From the equation 1-1 and 1-2, the total energy E (J) is

(1-3)

This results agreed with the observed hydrogen spectrum, and this model also applies to any one-electron atom. In 1923 Louis de Broglie suggested that electrons might have a wave aspect and its wavelength λ is equal to h / mv. In 1927, Davisson and Germer experimentally confirmed de Broglie's hypothesis by reflecting electrons from metals and observing diffraction effects [1]. Recently the results of the two-slit experiment of an electron showed its wavelike properties [2].

If the circular orbital length is equal to a integer times the wavelength of the electron, we have, 2 π r = n × (h / mv), mvr = n × (h / 2 π ), so it is consistent with the equation 1-1 of the Bohr theory. But the Bohr model could not explain about the spin of the electron and the two-electron atoms such as the helium. Because of such problems, the Bohr theory was replaced by the quantum mechanical theory based on the Schroedinger equation in 1920's. The absolute square of the wavefunction Ψ gives probability density of electrons. But in quantum mechanics, it is impossible to specify both the position and the momentum of an electron at some instant of time. The solution of the Schroedinger wave equation showed that the orbital angular momentum of the electron in the ground state of the hydrogen atom is zero. So the Coulomb potential may be infinitely negative, when the electron is close to the nucleus. And the solution also showed the relation between L (the total angular momentum quantum number) and M (the quantum number of the z component of the angular momentum),

M = -L, -L+1, -L+2, .......... L-2, L-1, L (2)

The relation 2 and the results of the Stern-Gerlach experiment indicated that an electron has ± 1/2 hbar spin angular momentum. But we can not visualize the spin rotational motion. By equating the angular momentum of the spinnig sphere of the electron to ± 1/2 hbar, the sphere speed leads to about one hundred times the speed of light. Experimentally, the spin magnetic moment μ is known to be

(3)

where g, the spin g-factor, is 2 [3], and β is the Bohr magneton (= 9.27400949 x 10^(-24) J/T). To solve the above problem, we try forgetting quantum mechanical methods and going back to the Bohr model. We suppose that the orbital angular momentum (S2) of the electron in the ground state of the hydrogen atom is ± 1 × hbar (here, spin is zero) and its g factor (g2) is 1. Then the magnetic moment of the electron in that state μ2 = - (g2)β(S2) / hbar is equal to μ of the equation 3. In this state, the problems about infinitely negative Coulomb potential and spinning sphere speed do not occur. As the magnetic moment of the electron is not changed, the splitting pattern of the energy levels in the magnetic field (Paschen-Back Effect) is also not changed.

The helium atom has two electrons and a nucleus of charge +2e. The Schroedinger equation for the herium can not be precisely solved. So the perturbation theory such as the variation method is applied to the calculation. Using the perturbation method, we can get the ground state energy of the helium atom, which is very close to the experimental value -79.005eV. But the calculation of the high order correction is very difficult.
In this paper, we try to calculate the ground state energy of the helium atom using the new theory based on the Bohr model, and check if the calculation value is equal to the experimental value -79.005eV.
First, suppose we have the following model (Fig.1). In that model, two electrons of the helium are on the opposite sides of the nucleus and moving on the same circular orbital.

Fig.1

Fig.1 The two electrons are moving on the opposite sides.

Equating the centrifugal force to the Coulomb force (the magnetic force is too small to be considered), we have

(4-1)

The circular orbital length is supposed to be an integer times the wavelength of the electron, we have

(4-2)

The total energy E of the electrons is the sum of the kinetic energy and the Coulomb potential energy, so

(4-3)

Solving the above three equations (4-1, 4-2, 4-3), the ground state energy (n=1) is -83.33eV. The value is lower than the experimental value -79.0eV. In this model, as the two electrons are on the same orbital, it is possible that the two electronic waves interfere with each other and their motions are affected, so the equations (4-1, 4-2, 4-3) may not be satisfied.
To avoid such problems, we suppose another model as shown in Fig.2 and Fig.3. In that model, the planes of the two same-shaped orbitals are perpendicular to each other. Any point on the electron 1 orbital is at the same distance from the points on the both-side( ± z) electron 2 orbital. On the electron 1 orbital, the wave of electron 2 may interfere with itself and vanish, so the motion of electron 1 may not be affected by the wave of electron 2.

Fig.2

Fig.2 The two same-shaped orbital planes are perpendicular to each other.
This figure shows one quarter of the orbitals

Fig.3

Fig.3 The two electrons have moved one quarter of their orbitals.

Methods and results :

The computer program (class filename: MathMethod) written in the JAVA language (version 1.5.0) to calculate the electronic orbital of the herium is shown in Fig.4.
As shown in Fig.2 and Fig.3, the helium nucleus is at the origin, the orbital plane of the electron 1 is the xy-plane, and the orbital plane of the electron 2 is xz-plane. The electron 1 initially at (r1, 0, 0) (Fig.2) moves one quarter of its orbital to (0, r2, 0) (Fig.3), while the electron 2 initially at (-r1, 0, 0) (Fig.2) moves to (0, 0, r2) (Fig.3). Meter and second are rather large units for measurement of atomic behavior, here we use new convenient units as follows,

unit

Fig.4 The program to calculate the electronic orbital of the helium

sample JAVA program

In this program, we first input the initial x-coordinate r1 (in MM) of the electron 1, and the absolute value of the total energy E (in eV) of the helium. From the inputted value, we calculate the initial velocity of the electron 1. And at intervals of 1 SS we calculate the Coulomb force among the two electrons and the nucleus. When the electron 1 is at (xx, yy, 0), the electron 2 is at (-xx, 0, yy) (in MM). Change MM to meter as follows; x (m) = xx × 10^(-14). y (m) = yy × 10^(-14). So the x component of the acceleration (m/sec^(2)) of the electron 1 is as follows;

(5-1)

where the first term is by the Coulomb force between the nucleus and the electron 1, and the second term is by the force between the two electrons.
. In the program, we use the reduced mass instead of the electron mass. In the same way, the y component of the acceleration (m/sec^(2)) is as follows;

(5-2)

Change m/sec^(2) to MM/SS^(2) using the next relation;

(6)

Based on that calculation value we change the velocity vector and the position of the electrons. We suppose electron 1 moves only on the xy-plane, so the z component of the acceleration of the electron 1 is neglected. We also calculate the wavelength of the electron from the velocity ( λ = h / mv) at intervals of 1 SS. The number of that wave ( λ in length) contained in that short movement section (The sum of them is WN) is,

(7)

where (VX, VY) are the velocity of the electron 1 (in MM/SS), the denominator is the wavelength (in m (meter)), and the numerator is the movement distance (in meter) for 1 SS.
When the electron 1 has moved one quarter of the orbital and its x-coordinate is zero (Fig.3), this program checked the y-component of the electron 1 velocity (last VY). As we want the last VY value which is close to zero in Fig.3, only when -0.0001 < VY < 0.0001 (MM/SS) is satisfied, the program displays the following values on the screen, r1, VY, preVY (VY 1SS ago), and (mid)WN (the total number of the electronic wave contained in one quarter of the orbital). The initial inputted x-coordinate of the electron 1 is automatically increased by 1 MM per above calculation to +100 MM.

Fig.5 shows the results in which the last VY is the closest to zero. Fig.6 graphs the results in Fig.5. As shown in Fig.5 and Fig.6, when the total energy of the helium (E) is -79.005 eV (the experimental value), WN x 4 is the closest to 1.

E (eV) r1 (MM) WN WN x 4
-77.000 3154.0 0.25324 1.01296
-77.500 3134.0 0.25242 1.00968
-78.000 3114.0 0.25161 1.00644
-78.500 3094.0 0.25080 1.00320
-79.000 3074.5 0.25001 1.00004
-79.005 3074.0 0.25000 1.00000
-79.500 3055.0 0.24922 0.99688
-80.000 3036.0 0.24844 0.99376
-80.500 3017.0 0.24767 0.99068
-81.000 2998.5 0.24690 0.98760

Fig.5 The results in which VY in Fig.3 is the closest to zero

graph

Fig.6 Plots of the number of the electronic wave contained in one orbital at various energy levels.

This results demonstrate that two electrons in the ground state of the helium may be actually moving around the nucleus on the orbitals of just one electronic wave length as shown in Fig.2 and Fig.3.

In this study, we do not use the quantum mechanical method such as the Schroedinger equation. The two main forces that affect the electronic motion are supposed to be the Coulomb force and the force by the de Broglie's wave. By a simple calculation, the magnetic force between the two electrons is too small in comparison to the Coulomb force. In the standard model, it is said that as the helium atom has spin-up and spin-down electrons, the magnetism does not occur. But strictly speaking, except in the part at the same distance from the two electrons, magnetic fields are produced by the electrons even in the standard helium model. So as the electrons moves, they lose energy by emitting electromagnetic waves. The same thing can be said about this new helium model. Probably the de Broglie's wave is related to the stabilization of the electronic motions. It is a very meaningful and astonishing thing that the calculation value by this new method based on the Bohr model is about the same as the experimental value. In the future, this method will be possibly to be a new useful tool for calculating electronic behavior.

References :

[1] C. Davisson, and L. H. Germer, "Scattering of Electrons by a Single Crystal of Nickel", Nature 119, 558, 1927

[2] A. Tonomura, J. Endo, T. Matsuda, et al., "Demonstration of Single-Electron Buildup of an Interference Pattern", Am. J. Phys., 57, 117, 1989

[3] G. Gabrielse, D. Hanneke, T. Kinoshita, et al., "New Determination of the Fine Structure Constant from the Electron g Value and QED", Phys. Rev. Lett., 97, 030802, 2006

to

2009/2/11 This site is link free.