scispace - formally typeset
Open AccessJournal ArticleDOI

The Wave Mechanics of an Atom with a Non-Coulomb Central Field. Part I. Theory and Methods

D. R. Hartree
- Vol. 24, Iss: 1, pp 89-110
Reads0
Chats0
TLDR
In this article, the characteristic values and functions of the wave equation of Schrodinger for a non-Coulomb central field were determined for the case where the potential is given as a function of the distance r from the nucleus.
Abstract
The paper is concerned with the practical determination of the characteristic values and functions of the wave equation of Schrodinger for a non-Coulomb central field, for which the potential is given as a function of the distance r from the nucleus.The method used is to integrate a modification of the equation outwards from initial conditions corresponding to a solution finite at r = 0, and inwards from initial conditions corresponding to a solution zero at r = ∞, with a trial value of the parameter (the energy) whose characteristic values are to be determined; the values of this parameter for which the two solutions fit at some convenient intermediate radius are the characteristic values required, and the solutions which so fit are the characteristic functions (§§ 2, 10).Modifications of the wave equation suitable for numerical work in different parts of the range of r are given (§§ 2, 3, 5), also exact equations for the variation of a solution with a variation in the potential or of the trial value of the energy (§ 4); the use of these variation equations in preference to a complete new integration of the equation for every trial change of field or of the energy parameter avoids a great deal of numerical work.For the range of r where the deviation from a Coulomb field is inappreciable, recurrence relations between different solutions of the wave equations which are zero at r = ∞, and correspond to terms with different values of the effective and subsidiary quantum numbers, are given and can be used to avoid carrying out the integration in each particular case (§§ 6, 7).Formulae for the calculation of first order perturbations due to the relativity variation of mass and to the spinning electron are given (§ 8).The method used for integrating the equations numerically is outlined (§ 9).

read more

Content maybe subject to copyright    Report

Dr Hariree,
The
wave mechanics
of
an atom, etc.
89
The Wave
Mechanics
of
an Atom tvith
a
Non-Voulomb Ventral
Field.
Part
I.
Theory
and
Methods.
By D. R.
HARTREE, Ph.D.,
St John's College.
[Received
19 November, read 21 November, 1927.]
§
1.
Introduction.
On
the
theory
of
atomic structure proposed
by
Bohr,
in
which
the electrons
are
considered
as
point charges revolving
in
orbits
about
the
nucleus,
the
Orbits being specified
by
quantum
con-
ditions,
it
is well known that both
a
qualitative and an approximate
quantitative explanation
of
many features
of the
simpler optical
spectra
and of
X-ray spectra
of
atoms with many electrons (e.g.
Rydberg sequences
in
optical spectra, term magnitudes
in
both
X-ray
and
optical spectra)
can be
given,
if
the assumption
is
made
that
the
effects
of
the electrons
on one
another can be represented
by supposing each to move in
a
central non-Coulomb field
of
force*;
further,
the
additional concept
of a
spinning electron provides
a
similar explanation
of
other features
of
these specbraf
(e.g.
doublet struoture
of
terms
and
magnitude
of
doublet separation,
anomalous Zeeman effect). This assumption
of a
central field
was
admittedly
a
rough approximation made
in the
absence
of any
detailed ideas about the interaction between the different electrons
in an atom, but in view
of
its success as
a
first approximation
for the
orbital atom model,
the
question arises whether, the same simple
approximations may
not
give useful results when applied
to the
new formulation
of
the quantum theory which
has
been developed
in
the
last two years.
The wave mechanics
of
Schrodinger} appears
to be the
most
suitable form
of the
new quantum theory
to use for
this purpose,
and will
be
adopted throughout. Further,
if i|r
is
a
solution
of
the
wave equation (suitably
normalised),
the suggestion has been made
by Schrodinger,
and
developed
by
Klein§, that |i/r|
2
gives
the
volume density
of
charge
in the
state described
by
this
i/r;
whether this interpretation
is
always applicable
may be
doubtful,
but
for the
wave functions corresponding
to
closed orbits
of
elec-
trons
in an
atom, with which alone this paper will
be
concerned,
it
has the
advantage that
it
gives something
of a
model both
of
the stationary states
(if i/r
only contains
one of the
characteristic
functions)
and of the
process
of
radiation
(if i/r is the sum of
* See,
for
example,
M. Bom,
Vorlesungen Uber
Atommechanik
(or the
English
translation, The
Mechanics
of
tlie
Atom),
Ch. in.
t
For a
general review,
see R. H.
Fowler, Nature, Vol. cxix,
p. 90
(1927)
for
a more detailed treatment,
F.
Hund, Liniewpektren,
Ch. m. '
%
E. Schrodinger, Ann. der Phys., Vol.
LXXIX,
pp.
361,
489; Vol.
LXXX
D 437-
Vol.
LXXXI,
p. 109
(1926); Phys.
Rev.,
Vol. xxrai,
p.
1049 (1926)
§
F.
Klein, Zeit.
f.
Phys., Vol.
su, p. 432
(1927).

90 Dr Hartree, The wave mechanics of an atom
a number of characteristic functions), and gives a simple interpre-
tation of the formula of the perturbation theory, namely, thdt the
change in the energy is the perturbing potential averaged over the
distribution of charge. Also in considering the scattering of radia-
tion of wave length not large compared to atomic dimensions, the
coherent radiation scattered by a hydrogen atom is given by
treating the scattering as classical scattering by the distribution
of charge given by Schrodinger's suggestion (if the wave length
is not too short)*, and this is probably true for any atom. For
the purpose of this paper the suggestion will be adopted literally;
the charge distribution for an atom in a stationary state is then
static (it does not necessarily follow that the charge itself is static).
Further, the distribution of charge for a closed n
k
group of
electrons is centrally symmetricalf, which suggests that on the
wave mechanics the assumption of a central field may give results
more satisfactory in detail than could be expected on the older
form of quantum theory.
Schrodinger's suggestion concerning the interpretation of
yfr
affords a hope that it may be possible to consider the internal field
of the atom as being due to the distribution of charge given by the
characteristic functions for the core electrons; we may, in fact,
attempt to find a field of force such that the total distribution of
charge, given by the characteristic functions in this field (taken in
suitable multiples corresponding to the numbers of electrons in
different n
k
groups), reproduces the field. The solution of this
problem, or rather a refinement of it, is, indeed, the main object of
the quantitative work to be considered here.
One point of contrast between the old ' orbital' mechanics and
the new wave mechanics may be emphasised here at once. On the
orbital mechanics the motion of an electron in an orbit lies wholly
between two radii (the potential energy and so the total energy
depends on the field at distances greater than the maximum
radius, but this does not affect the motion); on the wave
mechanics the solution of the wave equation is different from zero
at all but a finite number of values of the radius, and depends on
the field at all distances (though certainly the solution is very
small, and depends on the field to an extent negligible in practice,
at very large radii and usually at very small radii).
It will be seen later
(§ 2)
that, for a given characteristic function,
it is possible to specify two points which may to a certain extent
be considered to mark the apses of the corresponding orbit in the
orbital atom model, but the characteristic function is not zero
* I. Waller,
Nature,
Vol. oxx, p.
155
(1927);
Phil.
Mag.
Nov.
1927
(Supplement).
f See A. Unsold, Ann. der Phys., Vol.
LXXZII,
p.
355,
§§ 5
and 6. Unsold proves
this by considering the energy of an indefinitely small charge in the
field
of a closed
group; it also follows directly for the charge distribution from Unsdld's formulae
'67),
(69).

with a nan-Coulomb central field 91
outside them; as a consequence, it is not possible to make a strict
separation between 'penetrating' and 'non-penetrating' solutions
of the wave equation as it was possible to divide orbits into two
such classes. This suggests that it would be interesting to enquire
what occurs in the wave mechanics in cases when on the orbital
mechanics an atom may have a ' penetrating' and a ' non-pene-
trating' orbit with the same quantum numbers* «*, and also
whether the comparatively large doublet separations of some terms
corresponding to orbits classed as 'non-penetrating' can be ex-
plained by the non-zero fraction of the total charge which must, on
the wave mechanics, lie inside the core.
Associated with this question is that of the assignment of the
principal quantum number n to a solution of the wave equation in
a non-Coulomb field; in the orbital mechanics this was assigned in
a perfectly definite way, to which there is no direct analogy on the
wave mechanics. In this paper I will be written for the subsidiary
quantum number, taking integer values from zero upwards, which
is less by unity than Bohr's azimuthal quantum number k; this
follows the practice adopted by various writers f; to avoid altering
a notation which has become familiar, k will be retained as a suffix
in referring to the quantum numbers of an electron,
so
that l=k
1.
It seems best to define n such that
n
k
=
n
l
+
l is the number
of values of the radius r for which i|r = 0, excluding r = 0 (if it is
a root) and r = oo; n
I is then the number of values of r for
which |^|
2
is a maximum. This agrees with Bohr's principal
quantum number n for the hydrogen atom.
Both in order to eliminate various universal constants from
the equations and also to avoid high powers of 10 in numerical
work, it is convenient to express quantities in terms of units, which
may be called ' atomic units.' defined as follows:
Unit of
length,
a
B
=
A
2
/4Tr
2
me
2
,
on the orbital mechanics the
radius of the
1-quantum
circular orbit of the H-atom with
fixed nucleus.
Unit of
charge,
e, the magnitude of the charge on the electron.
Unit of
mass,
m, the mass of the electron.
Consistent with these are:
Unit of
action,
hl2ir.
Unit of
energy,
e*/a
= potential energy of charge e at distance a
from an equal charge = 2hcR =
twice
the ionisation energy
of the hydrogen atom with fixed nucleus.
Unit of time,
1/47TC.R.
* For example, on the orbital mechanics, Eb, Cn,
Ag,
Au, have 3
3
X-ray orbits,
and for the neutral atoms of these elements the first d term corresponds to a non-
penetrating
3
3
orbit.
t See, for example, P. Hand, op. dt., passim.

92
Dr
Hartree, The wave mechanics
of
an atom
These units being consistent,
the
ordinary equation?
of
classical
and wave mechanics hold
in
them;
in
particular Schrodinger's
wave equation for the motion of
a
point electron with total energy
E,
in
a
static field
in
which its potential energy
is
V, becomes
V°-^+2(E-V)yJr=*0.
(11)
For
an
attractive field,
V
is always negative;
it
is
convenient
to write
V=-v,
(1-2)
so that
the
quantity
v
calculated
in
numerical work
is
usually
positive.
For terms
of
optical
and
X-ray spectra,
E
is
always negative.
It
is
convenient
to
write
#
=
-H
(1-3)
the factor
\
being introduced
in
order that
e
shall
be
the energy
as
a
multiple
of
the ionisation energy
of
the hydrogen atom;
if a
characteristic value
of
the
solution
of
the
wave equation gives
directly
a
spectral term, of wave number v, then
e = v/R.
(1-4)
The wave equation
in
terms of v and
e is
-e)y}r
^
0.
The present paper
is
divided into
two
parts; this,
the
first,
deals with the methods used
for
solving this equation
for
a
given
non-Coulomb central field, and with the relevant theory; the second
with the question
of
the determination
of
the potential
v,
and with
an account and discussion
of
the results
for
some actual atoms.
§ 2. Theory.
When
the
field
is spherically symmetrical so that v is
a
function
of
the
radius
r
only, and spherical polar coordinates
r,
6,
<f>
are
used,
yfr
separates into
a
product
of a
function
^
(r)
of
r
only, and
a surface spherical harmonic
S
(#,
<f>);
if
I
is
the
order
of
the
spherical harmonic, the function
%
satisfies the equation
i%
+
\
2v
-
e
jj^
x
=o
(2-D
r
dr
\_
r
2
J
A
v
or, writing
'
P=r%
(2"2)
and using dashes
to
denote differentiation with respect
to
r,
P" +
\2v-e-l(l+l)jr*)P=0
(2-3)
This
is
the form
in
which the wave equation
is
used in the greater
part
of
this paper.

with a non-Coulomb central field 93
There are three advantages in working with P rather than
with
•%;
first,, the differential equation is simpler, secondly, P
2
gives the radial density of charge if P is suitably normalised
(i.e.
P*dr I1
P>dr
is the charge lying between radii r and r+dr) *,
and this is the quantity often required in applications, thirdly, as
a consequence of this, -P
2
is the weighting function for the per-
turbing potential at different radii, in the case of a perturbation
which is centrally symmetrical.
If we recall that the wave equation is derived from the classical
Hamiltonian equation of the problem by the substitution
p
x
= i^, etc.
(in atomic units), and also that on the orbital mechanics the
radial momentum p
r
in an orbit of angular momentum k (in
atomic units
hj^-w)
is given by
p
r
*
= 2v
e
--
equation (23) suggests that as far as we can picture an orbit
corresponding to a given solution of the wave equation, its angular
momentum in atomic units is given by
]<?=l(J,+
1) (I integral),
and that its apses are given by the roots of
2v-e-l(l+l)/r*
=
0,
(2"4)
i.e. by the points of inflexion of P other than those which occur at the
points where P =
0.
Usually this expression has two roots, between
which it is positive; between them P has an oscillatory character,
outside them it has an exponential character; it may have one
root only (in the case I = 0) or four (for I > 0 only). On the orbit
model of the atom, when four roots occur they give the apses of
an internal and an external orbit with the same energy f; what
then happens on the wave mechanics will be discussed in the
second part of this paper.
The first requirement is a method for finding values of e for
which, given tiasa function of r, the solution P of (2"3) is zero at
7=0 and
GO
(Schrodinger's condition is that i/r should be finite
* According
to
Schrodinger's interpretation
of
tf/
(§ 1)
the
charge lying
in an
element
of
volume denned
by
drd9d<f>
is
fr2sir\*6drd0d<t>l\\p
2
r
i
a\ti
i
8drd8d<l>,
the
integral being over
all
space,
so
that
the
charge lying between radii
r and
r
+ dr
is
1I
^r'dr
=
P
2
dr
/ I
"'P*dr,
the integration
of
the spherical harmonic factor cancelling
out.
t
In
general
the
internal and external orbits with the same energy will
not
both
be quantum orbits,
but
when they occur
it is
usually possible (always
if
integral
quantum numbers are used)
to
obtain
an
internal
and an
external quantum orbit
with
the
same quantum numbers.

Citations
More filters
Journal ArticleDOI

The density functional formalism, its applications and prospects

TL;DR: In this paper, the authors survey the local density functional formalism and some of its applications and discuss the reasons for the successes and failures of the local-density approximation and some modifications.
Journal ArticleDOI

Electronic excitations: density-functional versus many-body Green's-function approaches

TL;DR: In this paper, the authors compare the theoretical and practical aspects of the two approaches and their specific numerical implementations, and present an overview of accomplishments and work in progress, as well as a comparison of both the Green's functions and the TDDFT approaches.
Journal ArticleDOI

Theory of Solutions of Molecules Containing Widely Separated Charges with Special Application to Zwitterions

TL;DR: In this paper, the electrical contribution to the chemical potential of an ion having an arbitrary charge distribution is calculated with the aid of the Debye-Huckel theory, and the calculation is based upon a general solution in polar coordinates of the approximate Debye Huckel equation, Δψ-κ2ψ=0.
Journal ArticleDOI

Perspective on density functional theory

TL;DR: This perspective reviews some recent progress and ongoing challenges in density functional theory.
Journal ArticleDOI

Time-Dependent Density Functional Theory

TL;DR: An overview of TDDFT from its theoretical foundations to several applications both in the linear and in the nonlinear regime is given.
References
More filters
Journal ArticleDOI

Quantisierung als Eigenwertproblem

Journal ArticleDOI

Anwendung der Quantenmechanik auf das Problem der anomalen Zeemaneffekte

TL;DR: In this paper, the authors untersuch das quantenmechanische Verhalten des durch the Comptonsche Hypothese charakterisierten Atommodells, das die Zeemaneffekte and die Feinstrukturen der Dublettspektren vollstandig erklart werden konnen.
Journal ArticleDOI

Eigenvalues and Whittaker's Function

A. S. Eddington
- 01 Jul 1927 - 
TL;DR: In this paper, it is shown that Schrodinger's differential equation for the problem of determining the eigenvalues and eigenfunctions for the hydrogen atom is fully treated in a standard text-book, Whittaker and Watson's Modern Analysis.
Journal ArticleDOI

On the Evaluation of Certain Integrals Important in the Theory of Quanta

TL;DR: A reduction of the integral (1) to a simple and convenient expression is given, quite indispensable in the case of hyperbolic orbits, and has an interest beyond the special case of the Kepler motions.
Journal ArticleDOI

The Transition from Ordinary Dispersion into Compton Effect

Ivar Waller
- 01 Jul 1927 - 
TL;DR: In this paper, a generalised dispersion formula for very short waves was proposed, holding also for waves short compared with atomic dimensions and making it possible to follow the gradual transformation of ordinary dispersion into the Compton effect.
Related Papers (5)
Frequently Asked Questions (11)
Q1. What are the contributions in this paper?

On the theory of atomic structure proposed by Bohr, in which the electrons are considered as point charges revolving in orbits about the nucleus, the Orbits being specified by quantum conditions, it is well known that both a qualitative and an approximate quantitative explanation of many features of the simpler optical spectra and of X-ray spectra of atoms with many electrons ( e. g. Rydberg sequences in optical spectra, term magnitudes in both X-ray and optical spectra ) can be given, if the assumption is made that the effects of the electrons on one another can be represented by supposing each to move in a central non-Coulomb field of force * ; further, the additional concept of a spinning electron provides a similar explanation of other features of these specbraf ( e. g. doublet struoture of terms and magnitude of doublet separation, anomalous Zeeman effect ). Further, if i|r is a solution of the wave equation ( suitably normalised ), the suggestion has been made by Schrodinger, and developed by Klein§, that |i/r| gives the volume density of charge in the state described by this i/r ; whether this interpretation is always applicable may be doubtful, but for the wave functions corresponding to closed orbits of electrons in an atom, with which alone this paper will be concerned, it has the advantage that it gives something of a model both of the stationary states ( if i/r only contains one of the characteristic functions ) and of the process of radiation ( if i/r is the sum of 

The normalisation integral P'dr is required in calculation ofJ o perturbations; for the series electron the main contribution to this integral is from values of r where the field is effectively that of a point charge C; using some of the results of the previous section an approximate recurrence formula for this integral can be found, and from it an approximate formula for the integral itself. 

For the second term on the right the authors have, on integrating by parts,.\\p- (FT dp = [P>PP f] - Jp ± (P>P>) dp •expansion of the differential coefficient under the integral, followed by substitution for P" from the differential equation for P (5"3) and integration of the remaining term by parts gives finally/p» (PJ dp = [P>PP'] - [PP 2] -ftp* - 2n*p +1 (I +1) - 1] P*dP,so that altogether(»• + 1, I) dp = k [P* PP'] - J [pi*] + \\ [(n*Y -l(l+l) + l]/i* (»•, I) dp - [ i (n» - l - l ) (n* + l)fJP*(n* -1, I)dp,the function P in the integrated terms being P(n*, I). 

The wave mechanics of an atom, etc.at r = 0, and inwards from initial conditions corresponding to a solution zero at r = oo, with a trial value of the parameter (the energy) whose characteristic values are to be determined; the values of this parameter for which the two solutions fit at some convenient intermediate radius are the characteristic values required, and the solutions which so fit are the characteristic functions (§§ 2, 10). 

As explained at the end of the last section, solutions of the equations for M with different values of n* and The authorsatisfying this condition have been calculated, so that from the values of P at a comparatively small radius (only large enough for the deviation from a Coulomb field to be inappreciable) the value of the normalisation integral can be found from (7-6), thus avoiding the numerical evaluation of the integral in each particular case. 

The relativity perturbation term in the classical Hamiltonian is (in ordinary units) — (l/2mc2)(E — V)2f, or in atomic unitsAv = £ a2 (e/2 - v)2 [a2 = (2ire-/cKf = 1/18800],so that Ae = (a2/4>)f (e - 2vfPidrl\\ P>dr. 

Jo /JoIf Z is the effective nuclear charge at any radius (the charge which, placed at the nucleus, would give the same field as the actual field at that radius), and 1 and s are the orbital and spin angular momentum vectors, the spinning electron perturbation term in the classical Hamiltonian is (in ordinary units) £ (eh/27rmcy (Z/r3) ls j , or in atomic units Ja2 (Z/r3) Is, so that(Z/r3) P*dr/f P*dr.f 

as is convenient in the numerical work, the arbitrary constant in P is taken so that, for small r, P is the same for all solutions with the same The author(strictly, so that the limit of P/rl+1 as r-*-0 is the same for all solutions), and the main part of the perturbation arises from small values of r (as is the case forthe two special perturbations considered), then The author2vP2dr will be J o approximately the same for all solutions with the same I, so that approximately§ 

I +1 (corresponding to the ' elliptical orbits') it has been found best to integrate equation (3*4) for £ from r = 0 out to about the first maximum of P, and equation (2-3) for P from there to a point rather beyond the last maximum of P ; . 

These terms vanish at the upper limit p = <x>; if the field were a Coulomb field for all r, they would become infinite at the lower limit except for integral values of n*. 

Doubtless it would be possible to derive formulae for integration over longer intervals with adequate accuracy, but the writer's experience in other similar work is that simple formulae and a large number of intervals are much preferable to complicated formulae and a small number of intervals.