Английская Википедия:Helium atom

Содержание

Текст

Текст

Текст

Текст

Текст

Текст

Текст

Текст

Текст

A helium atom is an atom of the chemical element helium. Helium is composed of two electrons bound by the electromagnetic force to a nucleus containing two protons along with two neutrons, depending on the isotope, held together by the strong force. Unlike for hydrogen, a closed-form solution to the Schrödinger equation for the helium atom has not been found. However, various approximations, such as the Hartree–Fock method, can be used to estimate the ground state energy and wavefunction of the atom. Historically, the first such helium spectrum calculation was done by Albrecht Unsöld in 1927. Its success was considered to be one of the earliest signs of validity of Schrödinger's wave mechanics.

Introduction

Файл:Helium-term-scheme.svg

Schematic termscheme for Para- and Orthohelium with one electron in ground state 1s and one excited electron.

The quantum mechanical description of the helium atom is of special interest, because it is the simplest multi-electron system and can be used to understand the concept of quantum entanglement. The Hamiltonian of helium, considered as a three-body system of two electrons and a nucleus and after separating out the centre-of-mass motion, can be written asШаблон:Citation needed <math display="block"> H(\mathbf{r}_1,\, \mathbf{r}_2) = \sum_{i=1,2}\left(-\frac{\hbar^2}{2\mu} \nabla^2_{r_i} -\frac{Ze^2}{4\pi\varepsilon_0 r_i}\right) - \frac{\hbar^2}{M} \nabla_{r_1} \cdot \nabla_{r_2} + \frac{e^2}{4\pi\varepsilon_0 r_{12}} </math>

where <math> \mu = \frac{mM}{m+M}</math> is the reduced mass of an electron with respect to the nucleus, <math> \mathbf{r}_1 </math> and <math> \mathbf{r}_2 </math> are the electron-nucleus distance vectors and <math> r_{12} = |\mathbf{r}_1 - \mathbf{r}_2| </math>. It is important to note, that it operates not in normal space, but in a 6-dimensional configuration space <math>(\mathbf{r}_1,\, \mathbf{r}_2)</math>. The nuclear charge, <math> Z </math> is 2 for helium. In the approximation of an infinitely heavy nucleus, <math> M = \infty </math> we have <math> \mu = m </math> and the mass polarization term <math display="inline">\frac{\hbar^2}{M} \nabla_{r_1} \cdot \nabla_{r_2}</math> disappears, so that in operator language, the Hamiltonian simplifies to:

<math display="block"> H = -\frac{1}{2m} \mathbf p_1^2 - \frac{1}{2m} \mathbf p_2^2 - \frac{kZe^2}{ \mathbf r_1} - \frac{kZe^2}{ \mathbf r_2}+\frac{k e^2}{|\mathbf r_1-\mathbf r_2|} . </math>

The wavefunction belongs to tensor product of combined spin states and combined spatial wavefunctions, and since Hamiltonian only acts on spatial wavefunctions, we can neglect spin states until after solving the spatial wavefunction. This is possible since for any general vector: <math display="inline"> |\Phi \rangle = \sum_{ij} c_{ij} |\varphi_i \rangle |\chi_j \rangle </math> where <math display="inline"> |\varphi_i \rangle </math> is combined spatial wavefunction and <math display="inline"> |\chi_i \rangle </math> is the combined spin component, the Hamiltonian operator since it only acts on the spatial component gives eigenvector equation:

<math> H |\Phi \rangle = \sum_{ij} c_{ij} H |\varphi_i \rangle |\chi_j \rangle = \sum_{ij} c_{ij} (H| \varphi_i \rangle ) |\chi_j \rangle = E \sum_{ij} c_{ij} |\varphi_i \rangle |\chi_j \rangle = \sum_{ij} c_{ij} E |\varphi_i \rangle |\chi_j \rangle </math>

which implies finding solutions for <math> H |\psi_j \rangle = E | \psi_j \rangle </math> where <math display="inline"> | \psi_j \rangle =\sum_i c_{ij} | \varphi_i \rangle </math> is a general combined spatial wavefunction. This energy, however, is not degenerate by the dimension of combined spin state because of symmetrization postulate which requires that physical solutions for identical fermions should be totally antisymmetric, imposing restriction on choice of <math display="inline"> |\chi_i \rangle </math> based on solutions <math display="inline"> | \psi_i \rangle </math>. Hence the solutions are of form: <math> |\psi_i \rangle| \chi_i \rangle </math> where <math display="inline"> | \psi_i \rangle </math> is the energy eigenket spatial wavefunction and <math display="inline"> |\chi_i \rangle </math> is spin wavefunction such that <math> |\psi_i \rangle| \chi_i \rangle </math> is antisymmetric and <math display="inline"> |\Phi \rangle </math> is merely some superpostion of these states.

Since the Hamiltonian is independent of spin, it commutes with all spin operators. Since it is also rotationally invariant, the total x, y or z component of angular momentum operator also commutes with the Hamiltonian. From these commutation relations, <math display="inline"> S_{\pm} </math> and <math display="inline"> L_{\pm} </math> also commutes with the Hamiltonian which implies that energy is independent of <math display="inline"> m_l </math> and <math display="inline"> m_s </math>. Although, the purely spatial form of the Hamiltonian implies that the energy is independent of <math display="inline"> s </math>, this would only be true in the absence of symmetrization postulate. Due to symmetrization postulate, the choice of <math display="inline"> s </math> will influence the type of wavefunction required by symmetrization postulate which would in turn influence the energy of the state.

Other operators that commute with Hamiltonian are spatial exchange operator and parity operator. However a good combination of mutually commuting operators are: <math display="inline"> L^2 </math>, <math display="inline"> L_z </math>, <math display="inline"> S^2 </math>, and <math display="inline"> S_z </math>. Hence the final solutions are given as:

<math display="block"> H | N, L, m_l, S, m_s \rangle = E_{nls} | N, L, m_l, S, m_s \rangle </math>

Where the energy is <math display="inline"> (2l+1) ( 2s + 1 ) </math> fold degenerate. For electrons, the total spin can have values of 0 or 1. A state with the quantum numbers: principal quantum number <math>n</math>, total spin <math>S</math>, angular quantum number <math>L</math> and total angular momentum <math>J = |L-S|, \dots, L+S</math> is denoted by <math>n^{2S+1}L_J</math>.

States corresponding to <math> S = 0</math>, are called parahelium (singlet state, so called as there exist <math display="inline"> 2s+1=1</math> state) and <math> S = 1</math> are called orthohelium (triplet state, so called as there exist <math display="inline"> 2s+1=3</math> states). Since spin exchange operator can be expressed in terms of dot product of spin vectors, eigenkets of spin exchange operators are also eigenkets of <math> S = (S_1 + S_2)^2 </math>. Hence parahelium can also be said to be the spin anti-symmetric state (singlet state) or orthohelium to be spin symmetric state (triplet state).

The singlet state is given as:<math display="block"> |\chi_{s=0,m=0} \rangle = \frac 1{\sqrt 2} (|\uparrow\rangle |\downarrow \rangle- |\downarrow\rangle|\uparrow \rangle) </math><math display="block"> \psi(\mathbf r_1, \mathbf r_2) =\psi(\mathbf r_2\,,\mathbf r_1) </math>

and triplet states are given as:<math display="block"> \begin{align} &|\chi_{s=1,m=0} \rangle = \frac 1{\sqrt 2} (|\uparrow\rangle|\downarrow \rangle+ |\downarrow\rangle|\uparrow\rangle)\; ; \\[4pt] &|\chi_{s=1,m=1} \rangle = \;|\uparrow\rangle|\uparrow\rangle\; ; \;\; |\chi_{s=1,m=-1} \rangle = \;|\downarrow\rangle|\downarrow \rangle \; . \end{align}</math><math display="block"> \psi(\mathbf r_1, \mathbf r_2) = -\psi(\mathbf r_2 \,, \mathbf r_1) </math>

as per symmetrization and total spin number requirement. It is observed that triplet states are symmetric and singlet states are antisymmetric. Since the total wavefunction is antisymmetric, a symmetric spatial wavefunction can only be paired with antisymmetric wavefunction and vice versa. Hence orthohelium (triplet state) has symmetric spin wavefunction but antisymmetric spatial wavefunction and parahelium (singlet state) has antisymmetric spin wavefunction but symmetric spatial wavefunction. Hence the type of wavefunction of each state is given above. The <math display="inline"> (2l+1) </math> degeneracy solely comes from this spatial wavefunction. Note that for <math display="inline"> l=0 </math>, there is no degeneracy in spatial wavefunction.

Alternatively, a more generalized representation of the above can be provided without considering the spatial and spin parts separately. This method is useful in situations where such manipulation is not possible, however, it can be applied wherever needed. Since the spin part is tensor product of spin Hilbert vector spaces, its basis can be represented by tensor product of each of the set, <math> \{ |\uparrow\rangle_1 \,, |\downarrow\rangle_1 \} </math> with each of the set, <math> \{ |\uparrow\rangle_2 \,, |\downarrow\rangle_2 \} </math>. Note that here <math> |\uparrow\rangle_1|\downarrow\rangle_2 \ne |\downarrow\rangle_1|\uparrow\rangle_2 </math> but are infact orthogonal. In the considered approximation (Pauli approximation), the wave function can be represented as a second order spinor with 4 components <math> \psi_{ij}(\mathbf{r}_1,\, \mathbf{r}_2)</math>, where the indices <math> i,j = \,\uparrow,\downarrow </math> describe the spin projection of both electrons in this coordinate system.Шаблон:Better source needed The usual normalization condition, <math display="inline"> \sum_{ij}\int d\mathbf{r}_1 d\mathbf{r}_2|\psi_{ij}|^2 = 1</math>, follows from the orthogonality of all <math> |i\rangle_1|j\rangle_2 </math> elements. This general spinor can be written as 2×2 matrix:

<math display="block"> \boldsymbol \psi = \begin{pmatrix}\psi_{\uparrow\uparrow} & \psi_{\uparrow\downarrow} \\ \psi_{\downarrow\uparrow} & \psi_{\downarrow\downarrow} \end{pmatrix}</math>

If the Hamiltonian had been spin dependent, we would not have been able to treat each these components independently as shown previously since the Hamiltonian need not act in the same manner for all four components.

The matrix can also be represented as a linear combination of any given basis of four orthogonal (in the vector-space of 2×2 matrices) constant matrices <math> \boldsymbol \sigma^i_k </math> with scalar function coefficients <math> \phi^k_i (\mathbf{r}_1,\, \mathbf{r}_2) </math> as <math display="inline"> \boldsymbol \psi = \sum_{ik}\phi^k_i (\mathbf{r}_1,\, \mathbf{r}_2) \boldsymbol \sigma^i_k </math>.

A convenient basis consists of one anti-symmetric matrix (with total spin <math> S = 0</math>, corresponding to a singlet state)<math display="block"> \boldsymbol \sigma^0_0 = \frac 1{\sqrt 2} \begin{pmatrix}0 & 1\\ -1 & 0 \end{pmatrix} = \frac 1{\sqrt 2} (\uparrow\downarrow - \downarrow\uparrow) </math>
and three symmetric matrices (with total spin <math> S = 1</math>, corresponding to a triplet state)

<math display="block"> \begin{align} &\boldsymbol \sigma^1_0 = \frac 1{\sqrt 2} \begin{pmatrix}0 & 1\\ 1 & 0 \end{pmatrix} = \frac 1{\sqrt 2} (\uparrow\downarrow + \downarrow\uparrow)\; ; \\[4pt] &\boldsymbol \sigma^1_1 = \begin{pmatrix}1 & 0\\ 0 & 0 \end{pmatrix} = \;\uparrow\uparrow\; ; \;\; \boldsymbol \sigma^1_{-1} = \begin{pmatrix}0 & 0\\ 0 & 1 \end{pmatrix} = \;\downarrow\downarrow \; . \end{align}</math>

It is easy to show, that the singlet state is invariant under all rotations (a scalar entity), while the triplet are spherical vector tensor representations of an ordinary space vector <math>(\sigma_x,\sigma_y,\sigma_z) </math>, with the three components: <math display="block"> \sigma_x = \frac 1{\sqrt 2} \begin{pmatrix}1 & 0\\ 0 & -1 \end{pmatrix} , \quad \sigma_y = \frac i{\sqrt 2} \begin{pmatrix}1 & 0\\ 0 & 1\end{pmatrix}, \quad \sigma_z = \frac 1{\sqrt 2} \begin{pmatrix}0 & 1\\ 1 & 0 \end{pmatrix} .</math> Since all spin interaction terms between the four components of <math> \boldsymbol \psi </math> in the above (scalar) Hamiltonian are neglected (e.g. an external magnetic field, or relativistic effects, like angular momentum coupling), the four Schrödinger equations can be solved independently. This is identical to the previously discussed method of finding spatial wavefunction eigenstates independently of the spin states, here spatial wavefunctions of different spin states correspond to the different components of the matrix.

The spin here only comes into play through the Pauli exclusion principle, which for fermions (like electrons) requires antisymmetry under simultaneous exchange of spin and coordinates (totally antisymmetric wavefunction condition) <math display="block">\boldsymbol \psi_{ij} (\mathbf{r}_1,\, \mathbf{r}_2) = - \boldsymbol \psi_{ji} (\mathbf{r}_2,\, \mathbf{r}_1).</math>

Parahelium is then the singlet state <math> \boldsymbol \phi = \psi_+ (\mathbf{r}_1,\, \mathbf{r}_2) \boldsymbol \sigma^0_0</math> with a symmetric spatial function <math>\psi_+ (\mathbf{r}_1,\, \mathbf{r}_2) = \psi_+ (\mathbf{r}_2,\, \mathbf{r}_1)</math> and orthohelium is the triplet state <math> \boldsymbol \phi_m = \psi_- (\mathbf{r}_1,\, \mathbf{r}_2)\boldsymbol \sigma^1_m,\; m =-1,0,1 </math> with an antisymmetric spatial function <math>\psi_-(\mathbf{r}_1,\, \mathbf{r}_2) = -\psi_- (\mathbf{r}_2,\, \mathbf{r}_1)</math>.

Approximation methods

Following from the above approximation, effectively reducing three body problem to two body problem, we have:

<math display="block"> H = -\frac{1}{2m} \mathbf p_1^2 - \frac{1}{2m} \mathbf p_2^2 - \frac{kZe^2}{ \mathbf r_1} - \frac{kZe^2}{ \mathbf r_2}+\frac{k e^2}{|\mathbf r_1-\mathbf r_2|} . </math>

This Hamiltonian for helium with two electrons can be written as a sum of two terms:<math display="block"> H = H_0 + H' </math>

where the zero-order unperturbed Hamiltonian is

<math display="block"> H_0 = -\frac{1}{2m} \mathbf p_1^2 -\frac{1}{2m} \mathbf p_2^2 - \frac{kZe^2}{r_1} - \frac{kZe^2}{r_2} </math>

while the perturbation term:<math display="block"> H' = \frac{k e^2}{|\mathbf r_1-\mathbf r_2|} </math>

is the electron-electron interaction. Шаблон:Math is just the sum of the two hydrogenic Hamiltonians: <math display="block"> H_0 = -\frac{1}{2m} \mathbf p_1^2 -\frac{1}{2m} \mathbf p_2^2 - \frac{kZe^2}{r_1} - \frac{kZe^2}{r_2} = H_1 + H_2 </math> where<math display="block"> H_i = -\frac{1}{2m} \mathbf p_i^2 - \frac{kZe^2}{r_i} </math>are independent Coulomb field Hamiltonian of each electron. Since the unperturbed Hamiltonian is a sum of two independent Hamiltonians (ie. are separable), the wavefunction must be of form <math display="inline"> | \psi \rangle = |\psi_1 \rangle |\psi_2 \rangle </math> where <math display="inline"> | \psi_1 \rangle </math> and <math display="inline"> |\psi_2 \rangle </math> are eigenkets of <math display="inline"> H_1 </math> and <math display="inline"> H_2 </math> respectively. However, the spatial wavefunction of the form <math display="inline"> \psi(\mathbf r_1, \mathbf r_2) = \psi_a(\mathbf r_1) \psi_b(\mathbf r_2) </math> need not correspond to physical states of identical electrons as per the symmetrization postulate. Thus, to obtain physical solutions symmetrization of the wavefunctions <math display="inline"> | \psi_1 \rangle </math> and <math display="inline"> |\psi_2 \rangle </math> is carried out.

The proper wave function then must be composed of the symmetric (+) and antisymmetric(−) linear combinations:<math display="block">\psi_\pm^{(0)}(\mathbf r_1 \,, \mathbf r_2) = \frac 1 {\sqrt 2}(\psi_a(\mathbf{r}_1) \psi_b(\mathbf{r}_2) \pm \psi_a(\mathbf{r}_2) \psi_b(\mathbf{r}_1)) </math>or for the special cases of <math> \psi_a = \psi_b</math> (both electrons have identical quantum numbers, parahelium only): <math>\psi_+^{(0)} = \psi_a(\mathbf{r}_1) \psi_a(\mathbf{r}_2) </math>.

This explains the absence of the <math>1^3 S_1</math> state (with <math> \psi_a = \psi_b = \psi_{1s} </math>) for orthohelium, where consequently <math>2^3 S_1</math> (with <math> \psi_a = \psi_{1s}, \psi_b = \psi_{2s} </math>) is the metastable ground state.

Note that all wavefunction obtained thus far cannot be separated into wavefunctions of each particle (even for electrons with identical <math display="inline"> (n_1\,, l_1 \,, m_1 ) </math> and <math display="inline"> (n_2\,, l_2 \,, m_2 ) </math> where wavefunction is <math display="inline"> \psi^{(0)}(\mathbf{r}_1, \mathbf{r}_2) = \psi_{n_1,\ell_1,m_1}(\mathbf{r}_1) \psi_{n_2,\ell_2,m_2}(\mathbf{r}_2) </math> because then, the spin of the electrons are in a superposition of different spin states: <math> |\uparrow\rangle|\downarrow\rangle </math> and <math> |\downarrow\rangle|\uparrow\rangle </math> from <math display="inline"> \sigma^0_0 </math>) ie. the wavefunctions are always in superposition of some kind. In other words, one cannot completely determine states <math display="inline"> (n\,, l \,, m_l \,, m_s) </math> of particle 1 and 2, or measurements of all details, <math display="inline"> (n\,, l \,, m_l \,, m_s) </math> of each electrons cannot be made on one particle without affecting the other. This follows since the wavefunction is always a superposition of different states where each electron has unique <math display="inline"> (n\,, l \,, m_l \,, m_s) </math>. This is in agreement with Pauli exclusion principle.

We can infer from these wavefunctions that <math display="inline"> E= E_a+E_b </math>.

The corresponding energies are:<math display="block"> E^{(0)}_{n_1,n_2} = E_{n_1} + E_{n_2} = - \frac{kZ^2e^2}{2a_0} \left[\frac{1}{n_1^2} + \frac{1}{n_2^2} \right] </math>

A good theoretical descriptions of helium including the perturbation term can be obtained within the Hartree–Fock and Thomas–Fermi approximations (see below).

The Hartree–Fock method is used for a variety of atomic systems. However it is just an approximation, and there are more accurate and efficient methods used today to solve atomic systems. The "many-body problem" for helium and other few electron systems can be solved quite accurately. For example, the ground state of helium is known to fifteen digits. In Hartree–Fock theory, the electrons are assumed to move in a potential created by the nucleus and the other electrons.

Ground state of Helium: Perturbation method

Since ground state corresponds to (1,0,0) state, there can only be one representation of such wavefunction whose spatial wavefunction is:<math display="block"> \psi_0(r_1\,, r_2) =\psi_{1,0,0}(r_1) \psi_{1,0,0}(r_2) = \frac{Z^3}{\pi a_0^3} e^{-\frac{Z(r_1 + r_2)}{a_0}} </math>

We note that the ground state energy of unperturbed helium atom as:<math display="block"> E_0 = - 4 \frac{k e^2}{a_0} = -108.8 \text{ eV} </math>Which is 30% larger than experimental data.

We can find the first order correction in energy due to electron repulsion in Hamiltonian <math display="inline"> \frac{k e^2}{|\mathbf r_1-\mathbf r_2|} </math>:<math display="block"> \Delta^{(1)}= \int \int \frac{Z^6}{\pi^2 a_0^6}\frac{ke^2}{|\mathbf r_1-\mathbf r_2|} e^{-\frac{2Z(r_1+r_2)}{a_0}} d^3\mathbf r_1 d^3\mathbf r_2 = \frac{5}{4} \frac{k e^2}{a_0} </math>

The energy for ground state of helium in first order becomes <math display="inline"> E \approx -74.8 \text{ eV} </math> compared to its experimental value of Шаблон:Val. A better approximation for ground state energy is obtained by choosing better trial wavefunction in variational method.

Screening effect

The energy that we obtained is too low because the repulsion term between the electrons was ignored, whose effect is to raise the energy levels. As Шаблон:Mvar gets bigger, our approach should yield better results, since the electron-electron repulsion term will get smaller.<math display="block"> \bar{H'} = \frac{ke^2}{r_{12}} - \frac{kZe^2}{r_1} -V(r_1) - \frac{kZe^2}{r_2} - V(r_2) </math> Шаблон:Math is a central potential which is chosen so that the effect of the perturbation <math> \bar{H'} </math> is small. The net effect of each electron on the motion of the other one is to screen somewhat the charge of the nucleus, so a simple guess for Шаблон:Math is <math display="block"> V(r) = -\frac{k(Z-S)e^2}{r} = - \frac{kZ_ee^2}{r} </math> where Шаблон:Mvar is a screening constant and the quantity Шаблон:Math is the effective charge. The potential is a Coulomb interaction, so the corresponding individual electron energies are given by <math display="block"> E_0 = -\frac{k(Z-S)^2e^2}{2a_0n^2} = - \frac{kZ_e^2e^2}{2a_0n^2} </math> and the corresponding spatial wave function is given by <math display="block"> \psi_0(r_1\,, r_2) =\sqrt{\frac{Z_e^3}{\pi a_0^3}} e^{-Z_e(r_1 + r_2)} </math>If Z_e was 1.70, that would make the expression above for the ground state energy agree with the experimental value E₀ = −2.903 a.u. of the ground state energy of helium. Since Шаблон:Math in this case, the screening constant is S = 0.30. For the ground state of helium, for the average shielding approximation, the screening effect of each electron on the other one is equivalent to about <math display="inline"> \frac{1}{3} </math> of the electric charge.

Ground state of Helium: The variational method

To obtain a more accurate energy the variational principle can be applied to the electron-electron potential Шаблон:Math using the wave function <math display="block"> \psi_0(\mathbf{r}_1,\, \mathbf{r}_2) = \frac{8}{\pi a^3} e^{-2(r_1+r_2)/a} </math> <math display="block"> \langle H \rangle = 8E_1 + \langle V_{ee} \rangle = 8E_1 + \left(\frac{e^2}{4\pi\varepsilon_0}\right) \left(\frac{8}{\pi a^3}\right)^2 \int \frac{e^{-4(r_1 + r_2)/a}} {|\mathbf{r}_1 - \mathbf{r}_2|}\, d^3\mathbf{r}_1 \, d^3\mathbf{r}_2 </math>

After integrating this, the result is: <math display="block"> \langle H \rangle = 8E_1 + \frac{5}{4a} \left(\frac{e^2}{4\pi\epsilon_0}\right) = 8E_1 - \frac{5}{2}E_1 = -109 + 34 = -75 \text{ eV} </math>

This is closer to the experimental value, but if a better trial wave function is used, an even more accurate answer could be obtained. An ideal wave function would be one that doesn't ignore the influence of the other electron. In other words, each electron represents a cloud of negative charge which somewhat shields the nucleus so that the other electron actually sees an effective nuclear charge Z that is less than 2. A wave function of this type is given by: <math display="block"> \psi(\mathbf{r}_1, \mathbf{r}_2) = \frac{Z^3}{\pi a^3} e^{-Z(r_1+r_2)/a} </math>

Treating Z as a variational parameter to minimize H. The Hamiltonian using the wave function above is given by: <math display="block"> \langle H \rangle = 2 Z^2 E_1 + 2(Z-2) \left(\frac{e^2}{4\pi\varepsilon_0}\right) \left\langle \frac{1}{r} \right\rangle + \left\langle V_{ee} \right\rangle </math>

After calculating the expectation value of <math display="inline"> \frac{1}{r} </math> and V_ee the expectation value of the Hamiltonian becomes: <math display="block"> \langle H \rangle = \left[-2Z^2 + \frac{27}{4}Z\right]E_1 </math>

The minimum value of Z needs to be calculated, so taking a derivative with respect to Z and setting the equation to 0 will give the minimum value of Z: <math display="block"> \frac{d}{dZ} \left(\left[-2Z^2 + \frac{27}{4}Z\right] E_1\right) = 0 </math> <math display="block"> Z = \frac{27}{16} \sim 1.69 </math>

This shows that the other electron somewhat shields the nucleus reducing the effective charge from 2 to 1.69. This result matches with experimental results closely as well as the calculations of effective Z in screening effect. Hence, we obtain the most accurate result yet: <math display="block"> \frac{1}{2} \left(\frac{3}{2}\right)^6 E_1 = -77.5 \text{ eV} </math>

Where again, Шаблон:Math represents the ionization energy of hydrogen.

By using more complicated/accurate wave functions, the ground state energy of helium has been calculated closer and closer to the experimental value Шаблон:Val. The variational approach has been refined to very high accuracy for a comprehensive regime of quantum states by G.W.F. Drake and co-workers as well as J.D. Morgan III, Jonathan Baker and Robert Hill using Hylleraas or Frankowski-Pekeris basis functions. One needs to include relativistic and quantum electrodynamic corrections to get full agreement with experiment to spectroscopic accuracy.

Perturbation theory for Helium

Consider the same setting where unperturbed Hamiltonian is:<math display="block"> H_0 = -\frac{1}{2} \nabla_{r_1}^2 - \frac{1}{2} \nabla_{r_2}^2 - \frac{kZe^2}{r_1} - \frac{kZe^2}{r_2} </math>and perturbation is electron repulsion: <math display="inline"> \frac{k e^2}{|\mathbf r_1-\mathbf r_2|} </math>.

In general, for (1s)(nl) state, in first order perturbation theory:<math display="block"> E = E_{100} + E_{nlm} + \Delta^{(1)} </math>with:<math display="block"> \Delta^{(1)}= I \pm J </math>where I is known as direct integral and J is known as exchange integral or exchange energy. If the combined spatial wavefunction is symmetric, its energy level has the + symbol in <math display="inline"> \Delta^{(1)} </math>, whereas for the antisymmetric combined spatial wavefunction, <math display="inline"> \Delta^{(1)} </math> has a minus symbol. Since due to the symmetrization postulate, the combined spatial wavefunctions differ on symmetric or antisymmetric nature, the J term is responsible for the splitting of energy levels between ortho and para helium states.

They are calculated as:

<math display="block"> I = \int \int |\psi_{100}(\mathbf r_1)|^2 |\psi_{nlm}(\mathbf r_2)|^2 \frac{ke^2}{|\mathbf r_1 - \mathbf r_2|} d^3\mathbf r_1 d^3\mathbf r_2 </math><math display="block"> J = \int \int \psi_{100}(\mathbf r_1)\psi_{nlm}(\mathbf r_2) \frac{ke^2}{|\mathbf r_1 - \mathbf r_2|} \psi_{nlm}^{\text{*}}(\mathbf r_1)\psi_{100}^{\text{*}}(\mathbf r_2) d^3\mathbf r_1 d^3\mathbf r_2 </math>The first integral is said to be analogous to classical potential due to Coulomb interaction, where the squares of wavefunctions are interpreted as electron density. However, no such classical analog exists for the J term. Using Green’s theorem, one can show that the J terms are always positive. From these, the diagram for energy level splitting can be roughly sketched. It also follows that for these states of helium, the energy of parallel spins cannot be more than that of antiparallel spins.

Experimental value of ionization energy

Helium's first ionization energy is Шаблон:Val. This value was measured experimentally. The theoretic value of Helium atom's second ionization energy is Шаблон:Val. The total ground state energy of the helium atom is Шаблон:Val, or Шаблон:Val, which equals Шаблон:Val.

See also

• Isotopes of helium
• Araki–Sucher correction
• Hydrogen molecular ion
• Lithium atom
• List of quantum-mechanical systems with analytical solutions
• Quantum field theory
• Quantum mechanics
• Quantum states
• Theoretical and experimental justification for the Schrödinger equation
• "Helium atom" on Wikiversity