Tips - Variational Principles and Methods in Theoretical PDF
Tips - Variational Principles and Methods in Theoretical PDF
This book brings together the essential ideas and methods behind current applica-
tions of variational theory in theoretical physics and chemistry. The emphasis is on
understanding physical and computational applications of variational methodology
rather than on rigorous mathematical formalism.
The text begins with an historical survey of familiar variational principles in
classical mechanics and optimization theory, then proceeds to develop the vari-
ational principles and formalism behind current computational methodology for
bound and continuum quantum states of interacting electrons in atoms, molecules,
and condensed matter. It covers multiple scattering theory, as applied to electrons
in condensed matter and in large molecules. The specific variational principles
developed for electron scattering are then extended to include a detailed presenta-
tion of contemporary methodology for electron-impact rotational and vibrational
excitation of molecules. The book also provides an introduction to the variational
theory of relativistic fields, including a detailed treatment of Lorentz and gauge
invariance for the nonabelian gauge field of modern electroweak theory.
Ideal for graduate students and researchers in any field that uses variational
methodology, this book is particularly suitable as a backup reference for lecture
courses in mathematical methods in physics and theoretical chemistry.
ROBERT K. NESBET
IBM Almaden Research Center
The Pitt Building, Trumpington Street, Cambridge, United Kingdom
http://www.cambridge.org
vii
viii Contents
As theoretical physics and chemistry have developed since the great quantum rev-
olution of the 1920s, there has been an explosive speciation of subfields, perhaps
comparable to the late Precambrian period in biological evolution. The result is
that these life-forms not only fail to interbreed, but can fail to find common ground
even when placed in proximity on a university campus. And yet, the underlying
intellectual DNA remains remarkably similar, in analogy to the findings of recent
research in biology. The purpose of this present text is to identify common strands
in the substrate of variational theory and to express them in a form that is intelligi-
ble to participants in these subfields. The goal is to make hard-won insights from
each line of development accessible to others, across the barriers that separate these
specialized intellectual niches.
Another great revolution was initiated in the last midcentury, with the introduction
of digital computers. In many subfields, there has been a fundamental change
in the attitude of practicing theoreticians toward their theory, primarily a change
of practical goals. There is no longer a well-defined barrier between theory for
the sake of understanding and theory for the sake of predicting quantitative data.
Given modern resources of computational power and the coevolving development
of efficient algorithms and widely accessible computer program tools, a formal
theoretical insight can often be exploited very rapidly, and verified by quantitative
implications for experiment. A growing archive records experimental controversies
that have been resolved by quantitative computational theory.
It has been said that mathematics is queen of the sciences. The variational branch
of mathematics is essential both for understanding and predicting the huge body
of observed data in physics and chemistry. Variational principles and methods lie
in the bedrock of theory as explanation, and theory as a quantitative computational
tool. Quite simply, this is the mathematical foundation of quantum theory, and
quantum theory is the foundation of all practical and empirical physics and chem-
istry, short of a unified theory of gravitation. With this in mind, the present text is
xiii
xiv Preface
subdivided into four parts. The first reviews the variational concepts and formalism
that developed over a long history prior to the discovery of quantum mechanics,
subdivided into chapters on history, on classical mechanics, and on applied math-
ematics (severely truncated out of respect for the vast literature already devoted
to this subject). The second part covers variational formalism and methodology in
subfields concerned with bound states in quantum mechanics. There are separate
chapters on time-independent quantum mechanics, on independent-electron mod-
els, which may at some point be extended to independent-fermion models as the
formalism of the Standard Model evolves, and on time-dependent theory and linear
response. The third part develops the variational theory of continuum states, includ-
ing chapters on multiple scattering theory (the essential formalism for electronic
structure calculations in condensed matter), on scattering theory relevant to the
true continuum state of a quantum target and an external fermion (with emphasis
on methodology for electron scattering by atoms and molecules), continuing to a
separate chapter on the currently developing theory of electron-impact rotational
and vibrational excitation of molecules. The fourth part develops variational theory
relevant to relativistic Lagrangian field theories. The single chapter in this part de-
rives the nonquantized field theory that underlies the quantized theory of the current
Standard Model of elementary particles.
This book grew out of review articles in specialized subfields, published by
the author over nearly fifty years, including a treatise on variational methods in
electron–atom scattering published in 1980. Currently relevant topics have been
extracted and brought up to date. References that go more deeply into each of the
topics treated here are included in the extensive bibliography. The purpose is to
set out the common basis of variational formalism, then to open up channels for
further exploration by any reader with specialized interests. The most recent source
of this text is a course of lectures given at the Scuola Normale Superiore, Pisa, Italy
in 1999. These lectures were presented under the present title, but concentrated
on the material in Parts I and II here. The author is indebted to Professor Renato
Colle, of Bologna and the Scuola Normale, for making arrangements that made
these lectures possible, and to the Scuola Normale Superiore for sponsoring the
lecture series.
I
Classical mathematics and physics
[5] Akhiezer, N.I. (1962). The Calculus of Variations (Blaisdell, New York).
[26] Blanchard, P. and Brüning, E. (1992). Variational Methods in Mathematical Physics
(Springer-Verlag, Berlin).
[78] Dieudonné, J. (1981). History of Functional Analysis (North-Holland, Amsterdam).
[147] Goldstine, H.H. (1980). A History of the Calculus of Variations from the 17th
through the 19th Century (Springer-Verlag, Berlin).
[210] Lanczos, C. (1966). Variational Principles of Mechanics (University of Toronto
Press, Toronto).
[322] Pars, L.A. (1962). An Introduction to the Calculus of Variations (Wiley, New York).
[436] Yourgrau, W. and Mandelstam, S. (1968). Variational Principles in Dynamics and
Quantum Theory, 3rd edition (Dover, New York).
The idea that laws of nature should satisfy a principle of simplicity goes back
at least to the Greek philosophers [436]. The anthropomorphic concept that the
engineering skill of a supreme creator should result in rules of least effort or of most
efficient use of resources leads directly to principles characterized by mathematical
extrema. For example, Aristotle (De Caelo) concluded that planetary orbits must be
perfect circles, because geometrical perfection is embodied in these curves: “. . . of
lines that return upon themselves the line which bounds the circle is the shortest.
That movement is swiftest which follows the shortest line”. Hero of Alexandria
(Catoptrics) proved perhaps the first scientific minimum principle, showing that
the path of a reflected ray of light is shortest if the angles of incidence and reflection
are equal.
The superiority of circular planetary orbits became almost a religious dogma
in the Christian era, intimately tied to the idea of the perfection of God and of
His creations. It was replaced by modern celestial mechanics only after centuries in
which the concept of esthetic perfection of the universe was gradually superseded by
a concept of esthetic perfection of a mathematical theory that could account for the
3
4 1 History of variational theory
where ds is a path element, and v is the phase velocity. Fermat’s principle is that
the value of the integral T should be stationary with respect to any infinitesimal
deviation of the path x(t) from its physical value. This is valid for geometrical optics
as a limiting case of wave optics. The mathematical statement is that δT = 0 for
all variations induced by displacements δx(t). In this and subsequent variational
formulas, differentials defined by the notation δ · · · are small increments evaluated
in the limit that quadratic infinitesimals can be neglected. Thus for sufficiently small
displacements δx(t), the integral T varies quadratically about its physical value. For
planar reflection consider a ray path from P : (−d, −h) to the observation point
Q : (−d, h) via an intermediate point (0, y) in the reflection plane x = 0. Elapsed
time in a uniform medium is
T (y) = d + (h + y) + d + (h − y)
2 2 2 2 v, (1.2)
The problem posed by Bernoulli is that of the brachistochrone. If two points are
connected by a wire whose shape is given by an unknown function y(x) in a vertical
plane, what shape function minimizes the time of descent of a bead sliding without
friction from the higher to the lower point? The mass of a bead moving under gravity
is not relevant. It can easily be verified by trial and error that a straight line does not
give the minimum time of passage. Always in such problems, conditions appropriate
to physically meaningful solution functions must be specified. Although this is a
vital issue in any mathematically rigorous variational calculus, such conditions
will be stated as simply as possible here, strongly dependent on each particular
application of the theory. Clearly the assumed wire in the brachistochrone problem
must have the physical properties of a wire. This requires y(x) to be continuous,
but does not exclude a vertical drop. Since no physical wire can have an exact
discontinuity of slope, it is reasonable to require velocity of motion along the wire
to be conserved at any such discontinuity, so that the hypothetical sliding bead does
not come to an abrupt stop or bounce with undetermined loss of momentum. It can
easily be verified that a vertical drop followed by a horizontal return to the smooth
brachistochrone curve always increases the time of passage. Thus such deviations
from continuity of the derivative function do not affect the optimal solution.
The calculus of variations [5, 322] is concerned with problems in which a function
is determined by a stationary variational principle. In its simplest form, the problem
is to find a function
x1 y(x) with specified values at end-points x0 , x1 such that the
integral J = x0 f (x, y, y )d x is stationary. The variational solution is derived
from
∂f ∂f
δJ = δy + δy dx = 0
∂y ∂y
after integrating by parts to eliminate δy (x). Because
∂f ∂ f x1 d ∂f
δy
d x = δy − δy d x,
∂y ∂ y x0 dx ∂y
δ J = 0 for fixed end-points δy(x0 ) = δy(x1 ) = 0 if
∂f d ∂f
− = 0. (1.5)
∂y dx ∂y
This is a simple example of the general form of Euler’s equation (1744), derived
directly from a variational expression.
Blanchard and Brüning [26] bring the history of the calculus of variations into
the twentieth century, as the source of contemporary developments in pure math-
ematics. A search for existence and uniqueness theorems for variational problems
engendered deep studies of the continuity and compactness of mathematical entities
1.2 The variational calculus 7
that generalize the simple intuitive definitions assumed by Euler and Lagrange. The
seemingly self-evident statement that, for free variations of the function y(x),
∂f d ∂f
− δyd x = 0
∂y dx ∂y
implies Euler’s equation, was first proven rigorously by Du Bois-Reymond in 1879.
With carefully stated conditions on the functions f and y, this made it possible to
prove the fundamental theorem of the variational calculus [26], on the existence of
extremal solutions of variational problems.
where
∂f ∂f ∂f y
= 0, = 0, = .
∂x ∂y ∂y 1 + y 2
In this example, Euler’s equation takes the form of the geodesic equation
d y
= 0.
d x 1 + y 2
on integrating by parts and using dxa = ẋa dt. The final term here can be replaced,
using the energy constraint condition. Then, using d ẋa = ẍa dt,
Q
Q
∂V
δA = m a ẋa · δxa | P − m a ẍa + · δxa dt.
a a P ∂xa
If the end-points are fixed, the integrated term vanishes, and A is stationary if
and only if the final integral vanishes. Since δxa is arbitrary, the integrand must
vanish, which is Newton’s law of motion. Hence Lagrange’s derivation proves that
the principle of least action is equivalent to Newtonian mechanics if energy is
conserved and end-point coordinates are specified.
2
Classical mechanics
becomes pk = ∂∂q̇T when kinetic energy is expressed as T ({q, q̇}). The equations
k
of motion ṗa = X transform into ṗ k = Q k . Although this can be shown by direct
11
12 2 Classical mechanics
Generalized momenta are defined by pk = ∂∂q̇Lk and generalized forces are defined
∂L
by Q k = ∂q k
. When applied using t as the independent variable, Euler’s variational
equation for the action integral I takes the form of Lagrange’s equations of motion
d ∂L ∂L
− = 0. (2.3)
dt ∂ q̇ ∂q
Hamilton considered variations δq(t) that define continuous generalized trajectories
with fixed end-points, q(t0 ) = q0 and q(t1 ) = q1 . The variational expression is
t1
∂L ∂L
δI = dt δ q̇(t) + δq(t) = 0.
t0 ∂ q̇ ∂q
Following the logic of Euler, after integrating by parts to replace the term in δ q̇ by
one in δq, this implies the Lagrangian equations of motion.
The dissipative term takes the form of a nonconservative force, Q̂ k , acting on the
generalized coordinate qk .
X j ({qk }, t) = 0,
ak j q̇ k + akt = 0,
k
14 2 Classical mechanics
M ẍ − Mg sin φ + λ = 0,
Mr 2 θ¨ − λr = 0,
r θ̇ − ẋ = 0.
derivatives are
∂G ∂G
= v(x, y), = −x,
∂y ∂u
expressed as functions of the original variables.
This transformation is used to define thermodynamic functions that depend
on easily measurable variables such as pressure and temperature [57, 403]. For
example, the thermodynamic function X :enthalpy is determined by V :volume,
S:entropy, T :temperature, and P:pressure through the differential form d X =
T dS + V dP. A Legendre transformation from the abstractly defined entropy S to
the directly measurable temperature T defines the Gibbs function G = X − T S.
Changes dG = −S dT + V dP of the Gibbs function are determined by temperature
and pressure changes. In an isothermal process, ∂G∂P
= V is the system volume. In
∂G
an isobaric process, ∂ T = −S directly measures the entropy.
motion are
∂ L ∂H
= q̇ k − wk = q̇ k − = 0;
∂ pk ∂ pk
d ∂ L ∂ L ∂H
− = ṗ k + = 0;
dt ∂ q̇ k ∂qk ∂qk
∂ L ∂L ∂H
= − pk = − = 0.
∂wk ∂wk ∂wk
The first and second of these equations are Hamilton’s equations of motion. The
∂L
first and third establish wk = q̇ k and pk = ∂w k
as dynamical conditions, equivalent
∂L
to pk = ∂ q̇ .
k
1 2
T = m(ṙ 2 + r 2 θ̇ ).
2
p2 p2
The Hamiltonian is H = 2mr + 2mrθ 2 + V (r ), where pr = ∂∂Tṙ = m ṙ , and pθ = ∂∂Tθ̇ =
mr 2 θ̇ . Using the relevant partial derivatives, Hamilton’s equations of motion are
pθ
ṗ θ = 0; θ̇ = ;
mr 2
2 ∂V pr
ṗr = mr θ̇ − ; ṙ = .
∂r m
1
m ṙ = pr = [2m(E − V (r )) − 2 /r 2 ] 2 .
dθ
= 2 [2m(E − V (r )) − 2 /r 2 ]− 2 .
1
dr r
2.3 Conservation laws 17
∂ ∂ H ∂ ∂ H
=
˙ q̇ + ṗ = − ,
k
∂qk k ∂ pk k k
∂qk ∂ pk ∂ pk ∂qk
xi = xi + α(t), α → 0.
= Lt − pt t = (T − V + E )t = 2T t ,
such that
∂ ∂
= (Lt ) − pt = 0.
∂t ∂t
Since kinetic energy is positive, the action integral A = 2T t dτ is nondecreas-
ing. This suggests using a global parameter τ = s defined by the Riemannian line
element
ds 2 = m i j dqi dq j ,
i, j
Because variations are restricted to those that conserve energy, 2T = 2(E − V ) and
1
ds = (2(E − V )) 2 dt along the varied trajectories. The action integral becomes
A = 2(E − V )dt = m i j q̇ i q̇ j dt = pi dqi .
i, j i
20 2 Classical mechanics
δ A = 0 subject to the energy constraint restates the principle of least action. When
the external potential function is constant, the definition of ds as a path element
implies that the system trajectory is a geodesic in the Riemann space defined by
the mass tensor m i j . This anticipates the profound geometrization of dynamics
introduced by Einstein in the general theory of relativity.
ds 2 = c2 dt 2 − d x 2 − dy 2 − dz 2 ,
where c is the constant speed of light. The space and time intervals can be measured
in any reference frame in uniform linear motion, referred to as an inertial frame.
Space and time intervals measured in two such inertial frames are related by a
linear Lorentz transformation, parametrized by the velocity and direction of relative
motion. It is convenient to introduce Minkowski space-time coordinates
d ∂ ∂
− = 0.
dτ ∂u µ ∂ xµ
2.5 Special relativity 21
[77] Dennis, J.E. and Schnabel, R.B. (1983). Numerical Methods for Unconstrained
Optimization and Nonlinear Equations (Prentice-Hall, Englewood Cliffs,
New Jersey).
[125] Fletcher, R. (1987). Practical Methods of Optimization, 2nd edition (Wiley,
New York).
[168] Hestenes, M.R. (1966). Calculus of Variations and Optimal Control Theory
(Wiley, New York).
[332] Pulay, P. (1987). Analytical derivative methods in quantum chemistry, Adv. Chem.
Phys. 69, 241–286.
[439] Zerner, M.C. (1989). Analytic derivative methods and geometry optimization, in
Modern Quantum Chemistry, eds. A. Szabo and N.S. Ostlund (McGraw-Hill,
New York), pp. 437–458.
25
26 3 Applied mathematics
I current displacement
M inductance reciprocal mass tensor
R resistance viscous force
C capacitance inverse spring constant
E external emf driving force
f 0 f 0 λ0 + f 0 f 1 λ1 − µ = f 0 F
f 1 f 0 λ0 + f 1 f 1 λ1 − µ = f 1 F
λ0 + λ1 = 1.
On solving these equations for the coefficients {λ j }, the solution of minimum norm
is the interpolated gradient vector P, such that µ = |P|2 0, at the interpolated
coordinate vector Q. The Lagrange multiplier µ in this method provides an estimate
of the residual error.
At each step of an iterative method for nonlinear optimization, the subsequent
coordinate step q must be estimated. The vector Q, interpolated in a selected
m-simplex of prior coordinate vectors, must be combined with an iterative esti-
mate of the component of q orthogonal to the hyperplane of the simplex.
28 3 Applied mathematics
A Hessian matrix
∂2W
Fi j = ,
∂qi ∂q j
is defined by small displacements about a given point, such that p = Fq. A
generalization of Newton’s method can be used to estimate the location of a zero
value of a vector from its given value and derivative. If an approximation to F −1 is
maintained and updated at each iterative step, then
q = q⊥ + qi λi ,
i
where
q⊥ = −F −1 p + (F −1 p · qi )[(q · q)]i−1
j qj.
i, j≤m
Alternatively, the interpolation and extrapolation steps can be combined, using the
formula [75]
q = Q − F −1 P = (qi − F −1 pi )λi .
i
At a general point, displaced from the extremum, gradients pi = ∂∂qWi do not vanish.
Newton’s formula estimates a displacement toward the extremum q = −Gp0 ,
where G is the inverse of the Hessian matrix F. Since F is not constant in a nonlinear
problem, if it cannot be computed directly it must be deduced from an initial
estimate followed by iterative updates. In many circumstances, the gradients at a
general point can be computed directly or estimated with useful accuracy. Then
each successive coordinate increment q is associated with a gradient increment
p. If the Hessian matrix is known, Fq = p for sufficiently small increments.
Given an estimated F 0 , it must in general be updated to be consistent with the
computed gradient increment p. This implies a linear formula for the increment of
F, Fq = p − F 0 q. Since this provides only n equations for the n 2 elements
of F, the incremental matrix cannot be determined from data obtained in a single
iterative step.
The practical problem is to find a way to use the information obtained in each
iterative step without making unjustified changes of the Hessian. This could be
3.3 Iterative update of the Hessian matrix 29
accomplished if the Hessian were projected onto two orthogonal subspaces, one of
which corresponds to a subset of m linearly independent increments q and the
corresponding p. A rank-m (Rm) update in the vector space spanned by these
coordinate increments is defined such that q† Fq = q† (p − F 0 q). This
condition is satisfied by
F = Fmm + F 0 − Pm F 0 Pm† ,
where Fmm q = p. Replacing this by a form that leaves the complementary
projection of F 0 unchanged,
F = Fmm + (I − Pm )F 0 (I − Pm† ),
d p † dq = dq † d p = dq † F dq
30 3 Applied mathematics
In test calculations [75] this algorithm was found to produce rapid convergence.
When combined with single-step (m = 1) BFGS update of the inverse Hessian, this
is a very efficient algorithm.
The GDIIS algorithm can locate saddle points (transition states), because it
searches specifically for a point at which all gradients vanish, independently of
the sign of the second derivative. A reaction path can be followed by selecting the
eigenvector of G that belongs to the eigenvalue of greatest magnitude [332].
k=1
i, j
2 2
II
Bound states in quantum mechanics
In 1926, Schrödinger [365] recognized that the variational theory of elliptical differ-
ential equations with fixed boundary conditions could produce a discrete eigenvalue
spectrum in agreement with the energy levels of Bohr’s model of the hydrogen atom.
This conceptually startling amalgam of classical ideas of particle and field turned
out to be correct. Within a few years, the new wave mechanics almost completely
replaced the ad hoc quantization of classical mechanics that characterized the “old”
quantum theory initiated by Bohr. Although the matrix mechanics of Heisenberg
was soon shown to be logically equivalent, the variational wave theory became the
standard basis of methodology in the physics of electrons.
The nonrelativistic Schrödinger theory is readily extended to systems of N in-
teracting electrons. The variational theory of finite N-electron systems (atoms and
molecules) is presented here. In this context, several important theorems that follow
from the variational formalism are also derived.
Hartree atomic units will be used here. In these units, the unit of action h, the
mass m of the electron, and the magnitude e of the electronic charge −e are all set
equal to unity. The velocity of light c is 1/α, where α is the fine-structure constant.
35
36 4 Time-independent quantum mechanics
The unit of length is a0 , the first Bohr radius of atomic hydrogen. The Hartree unit
of energy is e2 /a0 , approximately 27.212 electron volts.
{H − }ψ = 0,
which implies
( f 1 (x)L[ f 2 (x)] − f 2 (x)L[ f 1 (x)]) d x = p( f 1 f 2 − f 1 f 2 )|xx10 .
[u(x)] + λu(x) = 0.
[G(x, ξ )] = −δ(x, ξ ),
where the Lagrange multiplier µ = λ−1 . The maximized value (u|G|u) is µ(u|u) =
λ−1
0 , selecting the lowest eigenvalue. Given u 0 and λ0 , consider the modified Green
function G 1 (s, t) = G(s, t) − u 0 (s)λ−1
0 u 0 (t). If (u|G 1 |u) is maximized, subject to
the orthogonality constraint (u|u 0 ) = 0 and to normalization (u|u) = 0, the same
procedure obtains u 1 and the eigenvalue λ1 . Equivalently, u n+1 maximizes (u|G|u)
subject to orthogonality to all eigenfunctions u i with i ≤ n. Thus an entire countable
sequence of eigenvalues and orthonormal eigenfunctions can be constructed.
If the full set of eigenfunctions is complete, any function χ in the Hilbert space
that is orthogonal to all eigenfunctions must vanish identically. To prove this, sup-
pose that some function χ exists such that (u i |χ ) = 0, i ≤ n, but (χ|χ ) = 1 and
λ(χ|G|χ ) = 1 for finite positive λ. The construction given above develops the
expansion
n
G(s, t) = lim u i (s)u i (t)/λi,
n→∞
i=0
H = t̂ + v(r) = t̂ − Z /r.
ˆ · Y
ˆ m = ( + 1)Ym
ˆz Ym = mYm .
40 4 Time-independent quantum mechanics
the surface term vanishes, leaving the volume term as defined in subsequent liter-
ature, T = (ψ|t̂|ψ). The surface term becomes important in scattering theory and
in solid-state theory using local atomic cells.
Defining the functional V = (ψ|v|ψ), where v(r) is a local potential for which
E = T + V is bounded below, the Schrödinger variational principle requires E
to be stationary subject to normalization (ψ|ψ) = 1. The variation δ E induced
δE
by δψ is expressed using the functional derivative δψ ∗ = Hψ, defined such that
3 ∗ δE
δ E = δ(ψ|H|ψ) = d r{δψ δψ ∗ + cc}. The variational condition is
∗ δE
δ(E − [(ψ|ψ) − 1]) = d r δψ 3
− ψ + cc = 0,
δψ ∗
for unconstrained δψ in the orbital Hilbert space, which implies the Schrödinger
equation
δE
= Hψ = {t̂ + v}ψ = ψ.
δψ ∗
Use of the Rayleigh quotient [ψ] = (ψ|H|ψ)/(ψ|ψ) as a variational func-
tional is an alternative to using a Lagrange multiplier to ensure normalization.
Assuming that the eigenfunctions {ψi } of H are complete and orthonormal in
the relevant Hilbert space, (ψ|ψ) = i |(i|ψ)|2 and (ψ|H|ψ) = i i |(i|ψ)|2 ≥
0 (ψ|ψ). This implies that [ψ] ≥ 0 , a variational principle for the ground
state. For other eigenstates, {H − i }ψi = 0, so that [φi ] = (i|i)/(i|i) = i .
Let ψ = ψi + δψ. Then δ(ψ|ψ) = (δψ|i) + (i|δψ), and δ(ψ|H|ψ) = (δψ|H|i) +
(i|H|δψ) = i δ(ψ|ψ). It then follows that
δ(ψ|H|ψ)(ψ|ψ) − (ψ|H|ψ)δ(ψ|ψ)
δ|ψ→ψi =
(ψ|ψ)2 ψ→ψi
δ(ψ|H|ψ) − i δ(ψ|ψ)
= = 0.
(ψi |ψi )
By construction, this stationary principle holds for any eigenstate.
which implies
{(a|H|b) − (a|b)}xb = 0.
b
In terms of the matrices (a|H|b) and (a|b), the set of coefficients xb is a null vector
of the matrix (a|H|b) − (a|b) when is determined so that the determinant of this
matrix vanishes. If the basis functions are orthonormalized, this condition becomes
det(h ab − δab ) = 0 and the coefficient vector is an eigenvector of the Hermitian
matrix h ab .
Suppose that h (n−1) is diagonalized in a basis of dimension n − 1, and this basis
is extended by adding an orthonormalized function ηn . The diagonalized matrix
is augmented by a final row and column, with elements h ni , h in respectively, for
i < n. The added diagonal element is h nn . Modified eigenvalues are determined
by the condition that the bordered determinant of the augmented matrix h (n) −
should vanish. This is expressed by
(n) |h ni |2
det h − = i<n (i − ) (h nn − ) − = 0.
−
i<n i
If the matrix elements h ni do not vanish, the final factor may vanish at values of
different from the original eigenvalues, giving the equation
|h ni |2
h nn − = .
−
i<n i
The function on the right here has poles at each of the original eigenvalues. It has
positive slope everywhere, because
d |h ni |2 |h ni |2
=
d i − (i − )2
for each term in the sum. This function increases from −∞ to ∞ between each
adjacent pair of unequal original eigenvalues, and must cross the straight line h nn −
once in each such interval.
Hence the new set of eigenvalues interleaves with the the original set. One new
eigenvalue lies below all of the originals, one lies above, and the rest each occur
in one of the intermediate intervals. This construction shows that the mth old eigen-
value is an upper bound for the mth new eigenvalue, for m < n. It follows that each
of the n eigenvalues of a variational matrix eigenvalue problem is an upper bound
for the corresponding eigenvalue of the Hamiltonian operator. If the eigenvalues
are indexed in nondecreasing sequence, and h nn falls between two adjacent nonde-
generate eigenvalues k−1 and k , the new eigenvalues i(n) for i = k are displaced
away from h nn . Thus i(n) ≤ i (i < k), and i+1
(n)
≥ i (i ≥ k). This displacement
4.2 Hellmann–Feynman and virial theorems 43
rule is evident from the 2 × 2 submatrix for ηn and any ψi . The secular equation is
(h nn − i(n) )(i − i(n) ) = |h ni |2 , which requires both factors to have the same sign.
lies in the Hilbert space of variational trial functions. In the limit of an infinitesimal
transformation, ξ → 0, (i Ĝ!|!) + (!|i Ĝ!) = 0, and
[H, Ĝ] = 0.
hypervirial theorem is
i
∂ H ∂H
[H, V] = · pi − ri · .
h i
∂pi ∂ri
4.3 The N-electron problem 45
If the Hilbert space of the variational basis set is invariant under scale transforma-
tion, then the hypervirial theorem implies
∂ H ∂H
· pi − ri · = 0. (4.1)
i
∂pi ∂ri
∂H
The force acting on nucleus a is Fa = − ∂Ra
. From Eqs. (4.1) and (4.2), the general
virial theorem for a molecule is
2T + V − Fa · Ra = 0.
a
For a diatomic molecule with internuclear distance R,
∂E
2T + V + R = 0.
∂R
coordinates,
(H − E)! = 0,
where H = T̂ + Û + V̂ . The individual operators, for kinetic energy, interelec-
tronic Coulomb interaction, and the external potential, respectively, are
1 2
T̂ = − ∇i ,
2 i
1 1 1
Û = = u(ri j ),
2 i, j ri j 2 i, j
V̂ = v(ri ).
i
Spin indices and summations will be suppressed in the notation here, but are to
be understood in connection with coordinate indices and integrations over coor-
dinates. Thus d 3 r · · · is intended to imply an implicit summation over a spin
index. The energy
3 functional 3 for an antisymmetric wave function !(1, . . . , N ) is
∗
(!|H
3 |!) = d r
13 1 ∗ N· · · d r N ! H ! and the normalization integral is (!|!) =
1 d r1 · · · N d r N ! !.
The Schrödinger variational principle requires (!|H |!) to be stationary subject
to constant normalization (!|!). Introducing a Lagrange multiplier, the variational
condition is
δ(!|H − E|!) = (δ!|H − E|!) + (!|H − E|δ!) = 0,
for variations in the N-electron Hilbert space. The Hermitian character of H implies
that (!|H − E|δ!) = (δ!|H − E|!)∗ . Since δ! is arbitrary, it can be multiplied
by i. Then (!|H − E|δ!) = −(δ!|H − E|!)∗ . Together, these equations imply
that (!|H − E|δ!) = 0 if and only if (δ!|H − E|!) = 0. Hence if δ! is unre-
stricted in the Hilbert space, the variational condition implies the time-independent
Schrödinger equation (H − E)! = 0.
from the residual unoccupied set. The general form is "µ = det{φ1µ (1) · · · φ Nµ (N )}.
If the orbital functions are indexed by i, j, . . . ≤ N for the occupied set, and by
N < a, b, . . . for the residual unoccupied set, a consistent notation for virtual ex-
citations is given by
Taking spin indices into account, all two-electron integrals are of the form
(i j|ū|kl) = (i j|u|kl) − (i j|u|lk), with the convention that orbitals with different
spin indices are orthogonal. It is convenient to truncate summations by the use of
occupation numbers n i , which are in principle determined by Fermi–Dirac statis-
tics. At zero temperature, occupation numbers are determined by the structure of the
reference state. Then n i = 1, n a = 0 for i ≤ N < a. A convention used for double
summation indices is i j : i < j ≤ N , ab : N < a < b.
With these conventions of notation, simple formulas exist for all matrix elements
of H in the basis of Slater determinants generated by virtual excitations from
a reference state [72]. Denoting the latter by "0 , and defining an effective one-
electron operator H = h + j n j (i j|ū|i j) for "0 , the sequence of nondiagonal
elements is
a
|H |0 = (a|H|i),
i
ab
ij |H |0 = n i n j (1 − n a )(1 − n b )(ab|ū|i j),
ij ab
abc
i jk |H |0 = 0, etc.
All such matrix elements vanish if more than two occupied orbitals are replaced,
because H contains only a two-electron interaction operator. Given (0|H |0) =
1
i n i (i|h|i) + 2 i, j n i n j (i j|ū|i j), diagonal elements follow a simple rule:
a
i |H |ia − (0|H |0) = (a|H|a) − (i|H|i) − (ai|ū|ai),
ab
i j |H |i j − (0|H |0) = (a|H|a) + (b|H|b) − (i|H|i) − ( j|H| j)
ab
The Hamiltonian is
1 1 1
H = − ∇12 − ∇22 + + O(1)
2 2 r12
1 1
= − ∇r2 − ∇q2 + + O(1).
4 q
Given relative angular momentum , the singular part of the Schrödinger equation is
∂2 2 ∂ ( + 1) 1
Hq = − − + + + O(1).
∂q 2 q ∂q q2 q
To avoid a singularity that cannot be cancelled by any one-electron potential, the co-
efficient of q −1 must vanish when = 0. This implies the Coulomb cusp condition
f 0 (q) = f 00 (1 + 12 q + · · ·). A similar expansion is valid for any > 0. Because the
4.4 Symmetry-adapted functions 49
1
mA
nB
(!µA |H |!ν B ) = (kµ /kν ) 2 xνi (" j A |H |"i B )
A∗ B
aµj
j=1 i=1
1
nA
mB
= (kν /kµ ) 2 aν j ("i A |H |" j B ).
A∗ B
xµi (4.3)
i=1 j=1
(!µA |H |!ν B ) =
! 12
A B A
B
xνi (" j A |H |"i B )
A∗ B
aµj aν j ("i A |H |" j B )
A∗ B
xµi . (4.4)
j i i j
Equations
i =2 3 4
1 −1 1
Coefficients kµ
− 12 1∗ 1 2
∗
1 1
2
− 12 3
2
"2 = det(aα bα cβ ),
"3 = det(aα cα bβ ),
"4 = det(bα cα aβ ).
H22 = E 0 − (ab|u|ba),
H23 = (bc|u|cb),
52 4 Time-independent quantum mechanics
H33 = E 0 − (ac|u|ca),
H24 = −(ac|u|ca),
H34 = (ab|u|ba),
H44 = E 0 − (bc|u|cb).
Matrix elements for the two orthonormal projected doublet (S = 12 ) functions,
computed using Eq. (4.4), are
1 1
(!2 |H |!2 ) = H22 + H23 − H24
2 2
1 1
= E 0 − (ab|u|ba) + (ac|u|ca) + (bc|u|cb),
2 2
1 1 1
(!3 |H |!2 ) = − H22 + H23 − H24
2 2 2
12
1 1
+ H32 + H33 − H34 [(H23 + H24 )]
2 2
12
3
= [−(ac|u|ca) + (bc|u|cb)][−(ac|u|ca) + (bc|u|cb)]
4
12
3
=− [(ac|u|ca) − (bc|u|cb)],
4
1
(!3 |H |!3 ) = − (H23 + H24 ) + (H33 + H34 )
2
1 1
= E 0 + (ab|u|ba) − (ac|u|ca) − (bc|u|cb).
2 2
5
Independent-electron models
Despite the simple and universal structure of the nonrelativistic Hamiltonian for
N interacting electrons, it produces a broad spectrum of physical and chemical
phenomena that are difficult to conceptualize within the full N-electron theory.
Starting with the work of Hartree [162] in the early years of quantum mechanics,
it was found to be very rewarding to develop a model of electrons that interact
only indirectly with each other, through a self-consistent mean field. A deeper
motivation lies in the fact that the relativistic quantum field theory of electrons is
53
54 5 Independent-electron models
{G − i }φi = 0,
occupied orbitals are “relaxed” to satisfy the OEL equations, Janak’s theorem takes
the form
∂E
dE = dn i = i dn i .
i
∂n i i
This theorem justifies the procedure of Slater and Wood for any version of orbital
functional theory.
δ Ex
= v̂ x φi = − n j ( j|u|i)φ j .
n i δφi∗ j
58 5 Independent-electron models
i
n i δφi∗ j
δE 3
+ δn i φi∗ d r = 0. (5.1)
i
n i δφi∗
δE
= Gφi = n j φ j λ ji . (5.2)
n i δφi∗ j
{G − i }φi = 0,
The coefficients ciabj can be obtained only by constructing the wave function !. An
exact formal expression can be derived by partitioning the N-electron Hilbert space
using the projection operator P = ""† and its orthogonal complement Q = I −
""† . Using the definition given above, the correlation energy is E c = ("|H |Q!).
A more explicit expression is obtained from the partitioned Schrödinger equation
P(H − E)P! + P(H − E)Q! = 0,
Q(H − E)P! + Q(H − E)Q! = 0.
Here QH Q defines a reduced Hamiltonian operator Ĥ . On substituting a formal
solution of the second set of equations into the first, it follows for η → 0+ that
E = ("|H − H ( Ĥ − E − iη)−1 H |"),
where ( Ĥ − E − iη)−1 denotes an inverse operator in the Q-space. This expresses
the correlation energy as an exact implicit functional of " and hence of its occupied
orbital functions [290],
E c = −("|H ( Ĥ − E 0 − E c − iη)−1 H |"). (5.3)
The coefficients ciabj are given by an expression that depends only on the reference
state and on a parametric value of E c ,
ciabj = "iabj |Q!
= − "iabj |( Ĥ − E 0 − E c − iη)−1 H |" .
Since H is specified, Eq. (5.3) defines E c as a functional of the occupied orbitals
of ", although it cannot be expressed as an explicit closed formula. The orbital
functional derivative nδi δφ
Ec
∗ = v̂c φi defines a nonlocal correlation potential v̂c in
i
60 5 Independent-electron models
exact OEL equations. Iterative solution for E c combines a search for a root of
the partitioned N-electron secular determinant with the Brueckner condition that
defines ". This would be a legitimate variational method if the inverse of matrix
Q(H − E)Q could be evaluated exactly. In practice this is not possible, but standard
approximations of many-body theory are available for particular applications of the
formalism.
An infinitesimal unitary transformation of the orbital basis that modifies " must
mix occupied and unoccupied orbital functions. For a typical orbital variation,
δ" = "ia δcia and δφi = φa δcia . Unitarity induces δφa = −φi δcia∗ = −φi (δφi |φa ).
In terms of functional derivatives δ/δ̃φ defined for independent variations of occu-
pied and unoccupied orbitals,
δ Ec
δ Ec = n i δφi
i n i δ̃φi∗
δ Ec
+ (1 − n a ) δφa + cc. (5.4)
a (1 − n a )δ̃φa
Treating the coefficients ciabj as constants in the current cycle of an iteration loop,
and using Eq. (5.5), (a|δ E c /n i δφi∗ ) defines the matrix element
1
(a|v̂c |i) = nj (1 − n c )(1 − n b )(a j|ū|cb)(cb|c̄|i j)
2 j b,c
1
− nk n j (1 − n b )(k j|ū|ib)(ab|c̄|k j), (5.6)
2 j,k b
in agreement with matrix elements deduced from the double virtual excitation terms
in Q! [275]. The effective correlation “potential” is a linear operator v̂c with the
kernel
1
vc (r, r ) = nj (1 − n c )(1 − n b )( j|ū|b)φc (r)φc∗ (r )(b|c̄| j)
2 j b,c
1
− nk n j (1 − n b )(b|c̄| j)φk (r)φk∗ (r )( j|ū|b). (5.7)
2 j,k b
5.3 Hartree–Fock theory 61
This subject will be treated in a later chapter on continuum states and scattering
theory.
o
1
+ ( j|u| j)χic − ( j|u|i)χ oj = Hc χic ,
j
2
δ("|H |")
c
= hχ o
+ 2( j|u| j)χio − ( j|u|i)χ cj
n i δχi
o o∗ i
j
o
+ ( j|u| j)χio − ( j|u|i)χ oj = Ho χio .
j
For comparison, the UHF effective Hamiltonians, indexed by spin α and β, are
δ("|H |")
α
β
α
α α∗ = hφ + ( j|u| j)φiα − ( j|u|i)φ αj + ( j|u| j)φiα = Hα φiα ,
n i δφi i
j j
δ("|H |") β
α
β
β
β β β
β β∗
= hφi + ( j|u| j)φi + ( j|u| j)φi − ( j|u|i)φ j = Hβ φi .
n i δφi j j
For the 1s 2 2s 2 S ground state of atomic Li, using the common notation Jn χx =
(n|u|n)χx , K n χx = (n|u|x)χn , the RHF Hamiltonians are
1
Hc = h + 2J1s + J2s − K 1s − K 2s ,
2
Ho = h + 2J1s + J2s − K 1s − K 2s ,
Hα = h + 2J1s + J2s − K 1s − K 2s ,
Hβ = h + 2J1s + J2s − K 1s .
This is the stationary condition for the energy functional, when orbital functions
are expanded in the specified basis set. RHF equations are related to UHF equations
exactly as they are in the theory based on integrodifferential OEL equations.
giving in general a different effective Hamiltonian for each orbital. The Euler–
Lagrange equations are to solved subject to the orthonormality constraints
(i| j) = δi j . Introducing Lagrange multipliers λ ji = ( j|Hi |i) = λi j = (i|H j | j),
these equations take the general form
Hi φi = φ j λ ji .
j
Orbital construction and diagonalization of the CI matrix are alternated until the
calculations converge.
or
where
gi (r, r ) = (1 − n a )φa (r)(a − i )−1 φa∗ (r ).
a
This defines δφi (r)/δvxc (r ) = −gi (r, r )φi (r ) as the kernel of a linear operator.
The variational equation induced by δvxc → δφi is
ni (1 − n a )(i|δvxc |a)(a − i )−1 (a|G|i) = 0.
i a
5.4 The optimized effective potential (OEP) 67
On substituting (a|GOEP |i) = 0 from the OEP orbital equation, and using (a|GOEP −
G|i) = (a|vxc − v̂xc |i), this gives the OEP integral equation
n i φi∗ (r) (1 − n a )φa (r)(a − i )−1 (a|vxc − v̂xc |i) = 0.
i a
1
I = ni (1 − n a )(i|vxc (r) − v̂xc |a)(a − i )−1 (a|vxc (r) − v̂xc |i),
2 i a
evaluated for solutions of the OEP orbital equations, is stationary for variations
of v . If this local potential function is expanded in a basis set πq , vxc (r) =
xc
q πq (r)Lq , the coefficient vector Lq is determined by linear algebraic equations
[69, 151]
(i| p|a)(a − i )−1 (a|q|i)Lq = (i| p|a)(a − i )−1 (a|v̂xc |i),
q i,a i,a
where (a|q|i) = φa∗ (r)πq (r)φi (r)d 3 r. The variational formalism [151, 181] pro-
vides a practical methodology for molecular wave functions by introducing a basis-
set expansion of the optimized local potential. Because only off-diagonal elements
(a|vxc |i) occur in the OEP equations, the local effective potential is determined
only up to an additive constant. The physical boundary condition that potentials
for an isolated system must vanish at infinite distance fixes a specific value of this
constant.
The variational energy E OEP = ("OEP |H |"OEP ) is minimized by the OEP equa-
tions subject only to the constraint of locality. Hence E OEP ≥ E UHF , the Hartree–
Fock ground-state energy. If an exact local exchange potential existed for UHF
ground states, it would be determined by the OEP, implying E OEP = E UHF . Calcu-
lations for the closed-shell atoms He, Be, and Ne [1, 97] obtain OEP energies E OEP
and UHF energies E UHF that agree for He, but E OEP , in Hartree units, is −14.5724
for Be and −128.5455 for Ne, above E UHF , −14.5730 for Be and −128.5471 for Ne,
by amounts greater than the residual computational errors. When E OEP > E UHF ,
an exact local exchange potential cannot exist.
68 5 Independent-electron models
then
where E[v] is the ground-state energy for external potential v(r). The minimizing
wave function is !v and the density is ρv . If v = vρ is chosen such that ρv = ρ,
5.5 Density functional theory (DFT) 69
this defines
F[ρ] = Fv [ρv ] = E[v] − vρv d 3 r,
δF
= fˆφi = {tˆ + vh (r) + v̂xc }φi .
n i δφi∗
5.5 Density functional theory (DFT) 71
if a functional derivative δ δρ
E xc
= vxc (r) exists in the form of a local potential function.
Assuming the existence of such a Fréchet functional derivative [26, 102] constitutes
the locality hypothesis. If this hypothesis were valid, the OEL and Kohn–Sham
equations would be equivalent, determining the same model or reference state.
It might appear paradoxical to invoke the locality hypothesis for vxc (r) while
using the Schrödinger operator tˆ rather than an assumed Fréchet derivative δT δρ
=
vT (r) for the kinetic energy. Applying Hohenberg–Kohn theory to a noninteracting
N-electron system, Kohn and Sham [205] show that the kinetic energy can be
expressed as a functional Ts [ρ] of the same ground-state density used to define
Fs [ρ]. Following the constrained-search logic of Levy [222], the density constraint
can be enforced by a Lagrange multiplier field w(r) that acts as a local effective
potential in the Kohn–Sham orbital equations. Thus it would appear that ground-
state theory for an interacting N-electron system can be replaced by a noninteracting
model constructed to obtain the same density, and that the theory requires only a
local potential w(r). Such a conclusion depends on the locality hypothesis, that
Fréchet functional derivatives exist for density functionals. This question is most
easily examined for a noninteracting system in the context of Thomas–Fermi theory,
where only the kinetic energy is relevant. Before examining this issue, consistency
conditions are derived for orbital and density functional derivatives of the same
functional.
If this becomes a density functional Fs [ρ] for ground states, the density functional
variation is
δ Fs
∗
δ Fs = d r 3
n i δφi (r) φi (r) + cc .
i
δρ(r)
72 5 Independent-electron models
An orbital functional derivative in general defines a linear operator such that niδδφ F
∗ =
δ Fs i
v̂ F φi . If the locality hypothesis were valid, then δρ(r) = v F (r), and the implied local
potential function could be computed directly from the sum rule,
v F (r)ρ(r) = n i φi∗ v̂ F φi .
i
This formula was used by Slater [385] to define an effective local exchange potential.
The generally unsatisfactory results obtained in calculations with this potential
indicate that the locality hypothesis fails for the density functional derivative of the
exchange energy E x [294].
" #
δ Fs [ρ] + V [ρ] − µ ρ d 3 r − N = 0,
or, if a Fréchet functional derivative δδρFs exists,
δ Fs
δρ + v(r) − µ d 3 r = 0.
δρ
If Fs is defined for unrestricted variations of ρ in any infinitesimal function neigh-
borhood of a solution, this implies the Thomas–Fermi (TF) equation
δ Fs
= µ − v(r).
δρ
The Lagrange multiplier µ, determined by normalization, is the chemical potential
[232], such that µ = ∂ E/∂ N when the indicated derivative is defined. This deriva-
tion requires the locality hypothesis, that a Fréchet derivative of Fs [ρ] exists as a
local function v F (r).
The locality hypothesis can be tested in a noninteracting model, in which the func-
tional Fs is replaced by Ts . The kinetic energy orbital functional is T = i n i (i|tˆ|i)
and the OEL equations are just orbital Schrödinger equations
tˆφi = {i − v(r)}φi . (5.8)
If the locality hypothesis is valid, then ∂∂ρTs = vT (r), and the Thomas–Fermi equation
is
vT (r) = µ − v(r). (5.9)
The trace sum i n i (i|vT − tˆ|i) must vanish if the OEL and TF equations are
equivalent. Equations (5.8) and (5.9), multiplied by appropriate factors, summed
over orbitals, and integrated, imply the sum rule [288]
n i i = N µ.
i
Since i n i = N and all i ≤ µ, this implies that all i are equal, in violation of
the exclusion principle for any compact system with more than two electrons. The
failure of this sum rule implies that in general the assumed Fréchet derivative of
Ts [ρ] cannot exist for more than two electrons, and there can be no exact Thomas–
Fermi theory.
If the hypothetical Fréchet derivative vT (r) could be replaced by the operator t̂
when acting on occupied orbital functions φi , there would be no contradiction. It
74 5 Independent-electron models
can be shown that this is the correct implication of variational theory. The functional
derivative in question can be shown to be a Gâteaux derivative [26, 102], the analog
in functional analysis of a partial derivative, equivalent to a linear operator that acts
on orbital functions. Consider variations of Ts induced by unconstrained orbital
variations about ground-state solutions of the OEL equations
δTs = d 3 r n i (δφi∗ t̂φi + cc) = d 3 r n i (δφi∗ {i − v}φi + cc)
i i
= d 3r n i {i − v(r)}δρi (r).
i
Since ρ = i n i ρi , this equation determines partial functional derivatives
δTs
= i − v(r). (5.10)
n i δρi
This defines a Gâteaux functional derivative [26, 102], whose value depends on a
“direction” in the function space, reducing to a Fréchet derivative only if all i are
equal. Defining H = t̂ + v, an explicit orbital index is not needed if Eq. (5.10) is
interpreted to define a linear operator acting on orbital wave functions, H − v =
t̂. The elementary chain rule nδT s
i δρi
= ni∂ρ∂ρi δT
δρ
s
= δT
δρ
s
is valid when the functional
derivatives are interpreted as linear operators. This confirms the chain rule, δT s
δρ i
φ =
δT
n i δφi∗
= t̂φi .
This argument shows that the locality hypothesis fails for more than two electrons
because the assumed Fréchet derivative must be generalized to a Gâteaux derivative,
equivalent in the context of OEL equations to a linear operator that acts on orbital
wave functions. The conclusion is that the use by Kohn and Sham of Schrödinger’s
operator t̂ is variationally correct, but no equivalent Thomas–Fermi theory exists
for more than two electrons. Empirical evidence (atomic shell structure, chemical
binding) supports the Kohn–Sham choice of the nonlocal kinetic energy operator,
in comparison with Thomas–Fermi theory [288]. A further implication is that if
an explicit approximate local density functional E xc is postulated, as in the local-
density approximation (LDA) [205], the resulting Kohn–Sham theory is variation-
ally correct. Typically, for E xc = exc (ρ)ρ d r, the density functional derivative
3
is a Fréchet derivative, the local potential function vxc = exc + ρ dexc /dρ.
If an exact local exchange potential does not exist, there is no reason for UHF,
OEP, and KSC results to be the same in the UHF model. That the OEP density is
not exactly equal to that of the UHF ground state is indicated by an analysis of OEP
results [412], and by recent test calculations [69]. This would imply that the density
constraint in KSC is a true variational constraint for more than two electrons, so that
E KSC > E OEP . The calculations considered here may not have sufficient numerical
accuracy to establish this evidently small energy difference.
6
Time-dependent theory and linear response
[79] Dirac, P.A.M. (1930). Note on exchange phenomena in the Thomas atom, Proc.
Camb. Phil. Soc. 26, 376–385.
[94] Ehrenreich, H. and Cohen, M.H. (1959). Self-consistent field approach to the
many-electron problem, Phys. Rev. 115, 786–790.
[155] Gross, E.K.U., Dobson, J.F. and Petersilka, M. (1996). Density functional theory of
time-dependent phenomena, in: Density Functional Theory, ed. R.F. Nalewajski
(Springer, Berlin).
[238] McLachlan, A.D. and Ball, M.A. (1964). Time-dependent Hartree-Fock theory for
molecules, Rev. Mod. Phys. 36, 844–855.
[292] Nesbet, R.K. (2001). Orbital functional theory of linear response and excitation, Int.
J. Quantum Chem. 86, 342–346.
[329] Petersilka, M., Gossmann, U.J. and Gross, E.K.U. (1996). Excitation energies from
time-dependent density-functional theory, Phys. Rev. Lett. 76, 1212–1215.
[351] Runge, E. and Gross, E.K.U. (1984). Density-functional theory for time-dependent
systems, Phys. Rev. Lett. 52, 997–1000.
[407] Thouless, D.J. (1961). The Quantum Mechanics of Many-body Systems (Academic
Press, New York).
77
78 6 Time-dependent theory and linear response
for variations of trial functions ψ that belong to the usual Hilbert space and satisfy
specified boundary conditions on the surface of the volume of integration. 3 The
Lagrange multiplier is determined so that the normalization integral d xψ ∗ ψ
remains constant. The resulting Schrödinger equation
{H − }ψ = 0,
is modified in time-dependent theory to
∂
ih ψ(xt) = H(xt)ψ(xt).
∂t
The simplified notation (xt) is used here to denote (x, t). If H is independent of
time, these two equations are equivalent. The wave function ψ(x) is modified by a
time-dependent phase factor exp(t/i h), which has no physical consequences.
Following Hamilton’s principle in classical mechanics, the required time depen-
dence can be derived from a variational principle based on a seemingly artificial
Lagrangian density, integrated over both space and time to define the functional
t1
∂
A[ψ] = dt d 3 xψ ∗ i h − H ψ.
t0 ∂t
Treating variations of ψ and ψ ∗ as independent, because they lead to equivalent
Euler–Lagrange equations, and integrating the time integral by parts as in Euler’s
theory, the resulting variational expression is
t1
∗ ∂
δA = dt d x δψ i h − H ψ + cc = 0,
3
t0 ∂t
3 ∗
if the normalization integral d xψ ψ is held constant. This implies the time-
dependent Schrödinger equation in the time interval t0 ≤ t ≤ t1 .
t integral A is not changed if the trial function ψ is multiplied by a phase
The action
factor exp( γ (t )dt /i h), while H is increased by a time-dependent but spatially
uniform potential γ (t). This is an example of gauge invariance, taken out of the
usual context of electromagnetic theory. Indicating the modified wave function by
ψγ , the modified action integral is
t1
∗ ∂
Aγ = dt d x δψγ i h − H − γ (t) ψγ + cc .
3
t0 ∂t
6.2 The independent-electron model as a quantum field theory 79
†
Occupation numbers are defined by n p = "|a p a p |", such that n i = 1, n a = 0
for i ≤ N < a. An (N−1)-electron basis state is defined such that |"i = ai |" if
i ≤ N . The orthonormality condition is verified by
†
"i |" j = "|ai a j |" = n i δi j .
†
An (N+1)-electron basis state is defined such that |"a = aa |" if a > N . Here
80 6 Time-dependent theory and linear response
†
"a |"b = "|aa ab |" = (1 − n a )δab .
∗
For cia (t) = X ia e−iωt + Yia∗ eiω t , treating frequency ω as a complex number,
(a − i − hω)X ia + (n i − n a ) (a j|ū|ib)X bj + (ab|ū|i j)Y jb
j b
= − (n i − n a )(a|w|i),
(a − i + hω)Yia + (n i − n a ) (i j|ū|ab)X bj + (ib|ū|a j)Y jb
j b
= − (n i − n a )(i|w|a).
∂
ih ρ̂ = [H, ρ̂].
∂t
Equivalently, time-dependent canonical Hartree–Fock equations are assumed to
take the same form as the time-dependent Schrödinger equation.
These concepts, inherent in the TDHF formalism, generalize immediately to
orbital functional theory, when electronic correlation energy is included in the
model. Given some definition ! → " that determines a reference state " for any
N-electron state !, correlation energy can be defined for any stationary state by
E c = E − E 0 , where E 0 = ("| Ĥ |") and E = ("| Ĥ |!). Conventional normaliza-
tion ("|!) = ("|") = 1 is assumed. A formally exact functional E c ["] exists for
stationary states, for which a mapping " → ! is established by the Schrödinger
equation [292]. Because both " and ρ̂ are defined by the occupied orbital func-
tions {φi }, for fixed occupation numbers n i , E["], E[ρ̂] and E[{φi }] are equivalent
functionals. Since E 0 is an explicit orbital functional, any approximation to E c as
an orbital functional defines a TOFT theory. Because a formally exact functional
E c exists for stationary states, linear response of such a state can also be described
by a formally exact TOFT theory. In nonperturbative time-dependent theory, total
energy is defined only as a mean value E(t), which lies outside the range of def-
inition of the exact orbital functional E c [{φi }] for stationary states. Although this
may preclude a formally exact TOFT theory, the formalism remains valid for any
model based on an approximate functional E c .
Any postulated orbital functional E[{φi }] defines an action integral,
t1
∗ ∂φi
A[{φi }] = dt n i d x φi i h
3
− E[{φi }] ,
t0 i
∂t
that is stationary for t0 ≤ t ≤ t1 if and only if the occupied orbital functions satisfy
the time-dependent Euler–Lagrange equations
∂φi (xt) δE
ih = = G(xt)φi (xt).
∂t n i δφi∗ (xt)
While E["] = E[!], the mean values of the time derivative are not equal. Even if
Ĥ is independent of time, these values are constants i n i i and E, respectively.
These constants can be reconciled by inserting a gauge potential γ = E − i n i i
into A, determined by equating the time-dependent phase factors. This also pre-
serves the normalization ("|!) = ("|") = 1 for the time-dependent wave func-
tions, so that the definition of correlation energy E c = ("| Ĥ |! − ") remains un-
changed. When the Hamiltonian is time-dependent, γ (t) must be determined so
that ("γ |i h ∂t∂ |"γ ) = i h ∂t∂ , where "γ incorporates the gauge-dependent phase
factor. For any rule ! → " applied t for all t0 ≤ t ≤ t1 , the time derivative terms in
the action integrals are #s = t01 dt("|i h ∂t∂ |") and # = t01 dti h ∂t∂ . The present
t
for variations of ! within the N-electron Hilbert space. Given a particular field v(xt),
the stationary condition determines !v (t) through the time-dependent Schrödinger
equations, but in general the density ρv (xt) computed from this wave function dif-
fers from the specified ρ(xt). A Lagrange-multiplier
t1 3 field vρ (xt) is to be chosen such
that ρv (xt) = ρ(xt).
Since vρ = t0 dt d x v(xt)ρ(xt) is constant for given v,
Av [ρ] = B[ρ] + vρ is defined as a functional of ρ, where B[ρ] = B[!v ], eval-
uated for v = vρ . This construction must succeed whenever the time-dependent
Schrödinger equation has a solution in the interval t0 ≤ t ≤ t1 with !(t0 ) = !0 ,
and if the Lagrange multiplier field vρ (xt) exists. This latter condition is “v-
representability”. These conditions seem very plausible for variations about well-
defined physical solutions, but there is no rigorous theory giving general uniqueness
and existence conditions. When this construction succeeds, it determines a func-
tional Av [ρ] that has a stationary value when ρ(xt) is determined by a solution of
the time-dependent Schrödinger equation defined by the external potential function
vρ (xt). Uniqueness follows from uniqueness of the solution of this time-dependent
equation [351], within a class of equivalent solutions corresponding to gauge trans-
formations. The trial functional B is “universal” in the sense that it depends only
on Ĥ − V̂ , the same for all external potentials. However, the variational condition
implies different results for different t initial conditions, !(t0 ) = !0 .
Defining the functional # = t01 dti h ∂t∂ , then Av = # − , where =
t1
t0 dt E v (t). Using E v (t) = Ĥ = f + V̂ = F + V , stationary Av implies
ˆ
88 6 Time-dependent theory and linear response
determines a density equal to the specified ρ(xt) when w(xt) = wρ (xt). “Nonin-
teracting v-representability” is established by any solution of these equations. This
defines a “universal” density functional Bs [ρ] when evaluated for ρ and wρ . For
fixed w(xt) but variable ρ, the functional Asw [ρ] = Bs [ρ] + wρ is stationary
when ρ = ρw , the density obtained by solving the noninteracting time-dependent
Schrödinger equation in which the time-dependent potential is wρ = w.
Defining the functional #s = t01 dt("|i h ∂t∂ |"), then Asw = #s − s , where
t
t1
s = t0 dt("|T̂ + Ŵ |"). Stationary Asw implies Euler–Lagrange equations of
the form
δ#s δTs
= + w(xt),
δρ δρ
if the indicated density functional derivatives exist as Fréchet derivatives.
In the noninteracting problem both required functional derivatives can be eval-
uated from explicit orbital functional derivatives. In detail, nδ# s
i δφi
∗ = ih
∂
φ and
∂t i
δTs ∗
n i δφi∗
= t φi . For variations of the partial densities ρi = φi φi , this implies Gâteaux
ˆ
functional derivatives [26, 102] nδ# s
i δρi
φi = i h ∂t∂ φi and nδT s
i δρi
φi = tˆφi . The implied or-
bital Euler–Lagrange equations are the time-dependent Schrödinger equations for
the occupied orbitals i ≤ N ,
∂
ih φi (xt) = (tˆ + w(xt))φi (xt).
∂t
Thus for noninteracting electrons, density functional analysis reproduces the usual
Schrödinger equations, but a theory based solely on densities (a time-dependent
6.7 Excitation energies and energy gaps 89
Thomas–Fermi theory) requires a Fréchet derivative for the kinetic energy, and
cannot exist for more than two electrons [288].
For interacting electrons, the required functional derivatives cannot be evaluated
explicitly. In the TDFT derivation [155], the density functional derivative of the
kinetic energy term in is treated correctly as equivalent to the Schrödinger operator
tˆ. Nonetheless, the functional derivative of the exchange-correlation functional is
assumed without proof to be a Fréchet derivative, equal to a time-dependent local
potential function vxc (xt). It is argued above that the difference functional # − #s
can consistently be set to zero by a gauge transformation. Taking this into account,
and assuming that vxc exists, both interacting and noninteracting equations have
the same form, and by construction produce the same density function ρ(xt). It
can then by inferred that the local potentials in these equations are equal, up to a
gauge potential γ (t). Then the Kohn–Sham potential could be decomposed into
w(xt) = v(xt) + vh (xt) + vxc (xt). Hence if a Fréchet derivative vxc did exist, the
correct interacting ρ(xt) would be obtained by solving the implied time-dependent
Kohn–Sham equation.
This proposition has been tested in the exact-exchange limit of the implied
linear-response theory [329]. The TDFT exchange response kernel disagrees qual-
itatively with the corresponding expression in Dirac’s TDHF theory [79, 289]. This
can be taken as evidence that an exact local exchange potential does not exist
in the form of a Fréchet derivative of the exchange energy functional in TDFT
theory.
93
94 7 Multiple scattering theory for molecules and solids
be subdivided into space-filling cells, each containing no more than one atomic
nucleus. A spherical polar coordinate system in each cell is centered about this
Coulombic singularity if it exists. The closed surface of cell τµ is denoted by σµ .
7.1.1 Definitions
µ µ
Solid-harmonic solutions JL and N L of the homogeneous Helmholtz equation in
cell µ are products of spherical Bessel functions and spherical harmonics. Specific
functional forms for the regular and irregular solid harmonics, respectively, are [188]
1
JL = κ 2 j (κr )Y,m (θ, φ),
1
N L = −κ 2 n (κr )Y,m (θ, φ).
The standard angular momentum indices , m, and an implicit spin index, are here
summarized by a single collective index L. Relative normalization is such that
(N L |Wσ |JL ) = δ L ,L , expressed in terms of Wronskian integrals over surface σ
defined by
(φ1 |Wσ |φ2 ) = [φ1∗ ∇n φ2 − (∇n φ1 )∗ φ2 ] dσ.
σ
Here ∇n denotes the outward normal gradient on σ . It can easily be verified from
the defining differential equation that the Helmholtz Green function is given by the
expansion
G 0 (r, r ) = − JL (r)N L∗ (r ), r < r ,
L
in any specified atomic cell [188]. The definition of irregular solid harmonics N L
here establishes boundary conditions appropriate to the principal value Green func-
tion in scattering theory.
form
µ
νµ
N L (r + Xν − Xµ ) = − JLν (r )g L ,L ,
L
µν
for r, r < |Xµ − Xν |. The Hermitian matrix g L ,L defines structure constants for
the two-center expansion of the Green function in locally regular solid harmonics.
Through the Lippmann–Schwinger equation, matrix g propagates a given wave
function in one cell into incremental wave functions in the other cell. MST collects
the total incremental wave function in each cell as a multicenter expansion, to
establish a linear consistency condition for a molecule or solid. For a regular periodic
lattice, contributions from outside an elementary polyatomic translational cell are
multiplied by translational phase factors and summed over an infinite space-lattice.
expansion of χ in cell τµ ,
χ (r) = G 0 vψ.
ν=µ τν
From the two-center expansion of the Green function, this equation implies
µν
χ (r) = − JL (r) g L ,L SLν ,
L L ν=µ
where SLν = − τν JL∗ vψ. On equating both expressions for χ , and matching coef-
ficients in the regular expansion about the coordinate origin in cell τµ ,
µ
µν
CL + g L ,L SLν = 0.
L ν=µ
This derivation implies that these equations are valid under the condition that no
adjacent nucleus should lie within the enclosing sphere of a given local cell.
The wave function ψ can be expanded as a linear combination of basis functions
defined separately in each atomic cell. It is convenient to define specific basis func-
tions φ L (r) constructed by solving the local Schrödinger equation at a specified
orbital energy and matching to JL (r) at the cell origin. Given these primitive basis
functions, in one-to-one correspondence with locally defined regular solid harmon-
ics, the local expansion is ψ = L φ L γ L . If the local potential is nonspherical,
each basis function φ L becomes a sum of spherical-harmonic components as it is
integrated outward from the local cell origin. Nevertheless, it is completely char-
acterized by the single index L, and by a spin index that is assumed in the notation
and discussion here. In this basis,
µν µν
δ L ,L C Lν ,L + g L ,L SLν ,L γ Lν = 0,
ν,L ,L
where
C L ,L = δ L ,L + N L∗ vφ L , SL ,L = − JL∗ vφ L . (7.1)
τ τ
If the matrix C is not singular, which requires the number of basis functions to
match the number of solid harmonics used to expand the Green function, a local
t-matrix is defined by t = −SC −1 . The consistency condition expressed above in
terms of C and S matrices then reduces to the simple matrix expression
(t −1 − g)Sγ = 0.
This is the fundamental equation of multiple scattering theory. It has the re-
markable property of concentrating effects of the local potential function into the
7.1 Full-potential multiple scattering theory 99
t-matrix, separately determined in each local atomic cell, while the geometry of
the polyatomic superstructure is entirely characterized by the matrix g of structure
constants. For a regular periodic solid, the matrices here are summed over equiv-
alent translational cells, with coefficient phase factors that depend on an effective
momentum vector k, to give
{t −1 () − g(; k)}S()γ (; k) = 0.
if the sums converge. The notation here denotes matching of both function value
and normal gradient. Even if the separate series in functions JL and N L do not
converge, it can be shown that well-defined functions exist, providing summation
formulas for the possibly divergent series. This justifies formal operations in MST
that depend on Wronskian integrals of this expansion.
100 7 Multiple scattering theory for molecules and solids
An interior function χ (r) = −(G 0 (r, rσ )|Wσ |ψ) is defined that is regular at the
cell origin. The one-center expansion of the Green function converges in the largest
sphere enclosed by surface σ , sphere S0 of radius r0 . This determines a series ex-
pansion χ(r) = L JL (N L |Wσ |ψ), valid for r < r0 . Similarly, an exterior func-
tion η(r) = (G 0 (r, rσ )|Wσ |ψ) is defined that is bounded at infinite separation.
Again using the convergent expansion of the Green function, but in the exterior
region outside an enclosing sphere S1 of radius r1 , η(r) = − L N L (JL |Wσ |ψ)
for r > r1 .
The auxiliary functions χ and η have finite values on σ because the surface area
element removes the singularity of the Green function in their defining surface inte-
grals. Both functions are defined throughout the region r0 ≤ r ≤ r1 . The boundary
conditions determine unique solutions of the Helmholtz equation for integration of
χ outward from S0 and of η inward from S1 . Similarly, the function ψ given on σ
can be extended throughout this region by integrating the Helmholtz equation both
inwards to S0 and outwards to S1 . Green’s theorem implies for such solutions of
the Helmholtz equation that Wronskian surface integrals are conserved on nested
closed surfaces. Hence the difference function ψ − χ − η has vanishing Wronskian
integrals with all regular and irregular solid harmonics on both spherical surfaces
S0 and S1 . It follows from Green’s theorem that this difference function must vanish
identically in the region between these two surfaces, and in particular on surface σ .
This establishes the surface matching theorem [280],
The series expansion of χ may diverge in the interior moon region [148], between
the enclosed sphere S0 and σ , while the corresponding series for η may diverge
between σ and S1 . Nonetheless, these functions are well defined by Wronskian
integrals of the Green function, independently of such expansions. The implied
surface integrals are properties of the integrated functions, and justify formal use
of the surface series expansion in boundary matching conditions based on these
Wronskian integrals. If applied to basis functions that are solutions of the same
Schrödinger equation in τ , for which (φi |Wσ |φ j ) = 0, surface expansion of either
function implies
(φi |Wσ |JL )(N L |Wσ |φ j ) − (φi |Wσ |N L )(JL |Wσ |φ j ) = 0.
L L
This equation determines the local t-matrix. Any local basis function φi has the
formal expansion on σ
φi ≡σ [JL (N L |Wσ |φi ) − N L (JL |Wσ |φi )].
L
Canonical basis orbitals are defined here such that φ L → JL at the local cell origin.
Introducing a matrix notation for surface Wronskian integrals,
C L ,L = (N L |Wσ |φ L ), SL ,L = −(JL |Wσ |φ L ), (7.2)
the formal expansion of a primitive basis orbital φ L on σ is
φ L ≡σ (JL C L ,L + N L SL ,L ).
L
The surface Wronskian integrals C and S defined here are identical with the volume
integrals defined in Eqs. (7.1). Using (∇ 2 + κ 2 )φ L = vφ L , integration by parts of
the MST volume integrals, Eqs. (7.1), gives equivalent surface integrals, verifying
Eqs. (7.2),
SL ,L = − JL∗ (∇ 2 + κ 2 )φ L = −(JL |Wσ |φ L ),
τ
C L ,L = δ L ,L + N L∗ (∇ 2 + κ 2 )φ L = (N L |Wσ |φ L ).
τ
the enclosed sphere of an atomic cell, and η outside. Here φ is a local solution of
the Schrödinger equation, while χ and η are sums of regular and irregular solid
harmonics, respectively. The matching on S0 is exact in the muffin-tin model,
so that these functions provide a basis set suitable for variational calculations.
In practice, because the long-range η functions are not orthogonal to inner shell
functions in displaced atoms [193], they are expanded in the energy derivatives φ̇ of
the basis functions φ. The resulting linearized muffin-tin orbital (LMTO) method
[13, 9, 384] solves the Schrödinger secular equation in an MTO basis, constructed
using KKR/MST structure constants.
The surface matching theorem makes it possible to generalize the idea of muffin-
tin orbitals to a nonspherical Wigner–Seitz cell τ . Each local basis orbital is rep-
resented as φ ≡σ χ + η on the cell surface σ , where χ and η are the auxiliary
functions defined by the surface matching theorem. An atomic-cell orbital (ACO)
is defined as the function φ − χ , regular inside τ . By construction, the smooth
continuation of this ACO outside τ is the function η. The specific functional
forms are
χL = JL C L ,L , ηL = N L SL ,L .
L L
Either MTO or ACO functions are valid as basis functions for expanding a global
wave function ψ in all atomic cells. By construction, they are regular in τ , smooth
at σ , and bounded outside. When the matrix C is nonsingular, modified canonical
basis functions can be defined such that
φ̂ L = φ L (C −1 ) L ,L ,
L
where β Lν= −(JLν |Wν |ψ).If the sum here diverges, specific integrated functions
ην should be used.
All Wronskian integrals of ψ − ξ should vanish on the honeycomb lattice. Given
µ µ µ
the local expansion ψ = L φ L γ L in any particular cell τµ , (JL |Wµ |ψ − ξ ) = 0
µ µ µ µ
determines β L = L SL ,L γ L . (N L |Wµ |ψ − ξ ) = 0 implies that
µν µν
δ L ,L C Lν ,L + g L ,L SLν ,L γ Lν = 0. (7.4)
ν,L ,L
These are the standard MST equations, in the form (t −1 − g)Sγ = 0, where
t = −SC −1 = (J |Wσ |φ)(N |Wσ |φ)−1 . (7.5)
(∇ 2 + κ 2 )G 0 (x, x ) = δ(x, x ).
where, integrating the finite volume integrals by parts, and using the coefficient
matrices defined in MST,
µ
χ =ψ−
in
G 0 vψ = JL C L ,
τµ L
µν
χ out = G 0 vψ = − JL g L ,L SLν .
ν=µ τν ν=µ L ,L
µ
In a finite basis of primitive functions φ L , the stationary condition is
µ†
µ†
µν µν
−δ = δγ L SL L 1 δ L 1 L 2 C Lν 2 L + g L 1 L 2 SLν 2 L γ Lν + hc = 0,
L ,L L 1 ,L 2 ν
which implies the standard MST equations. Energy eigenvalues are determined by
adjusting the trial energy so that the secular determinant of the matrix of MST
equations is singular. Continuity conditions at cell boundaries are not addressed
by this variational principle. It is usually assumed that ψ lies in the orbital Hilbert
space (continuous function value and normal gradient at any cell interface). These
conditions are required for consistency with the Schrödinger equation, but are not
required by the KR variational principle and the MST equations.
106 7 Multiple scattering theory for molecules and solids
which can be evaluated in closed form. The contracted matrix product gS can
also be evaluated in closed form by defining a generalized C-matrix, for displaced
µν µ µ
irregular solid harmonics N , C̃ L 1 ,L = (N L 1 |Wν |φ Lν ). The function N L 1 defined in
cell τµ is evaluated on the surface σν of a displaced cell. This is justified because
N -functions are valid outside their cell of definition. Using the Green function to
µν µν µν
propagate N µ , this becomes C̃ L 1 ,L = − L 2 g L 1 L 2 (JLν2 |Wν |φ Lν ) = L 2 g L 1 L 2 SLν 2 ,L .
In a matrix notation omitting L indices and summations, the KKR/MST equation is
S µ† (δ µν C ν + g µν S ν )γ ν = µν γ ν = 0,
ν ν
where all matrix products are contracted over the long dimension. In the Hermitian
µν †
matrix the triple product S µ† g µν S ν = S µ† C̃ = (C̃ )µν S ν . Multiplying on the
†
left by S(S † C)−1 , and using (C̃ )µν = S µ† g µν , this takes the form β = 0, defining
a modified secular matrix
†
µν = δ µν + S µ (S µ† C µ )−1 (C̃ )µν .
The coefficients β = Sγ are the expansion coefficients of the global matching
µ µ
function ξ (σ ) = µ,L N L (rσ )β L . Defining t = −S(S † C)−1 S † , which reduces to
t = −SC for square matrices [281], µν takes the standard MST form I − tg,
−1
without implying inversion of a rectangular matrix. Using these closed forms, the
coefficients β and γ are related by
†
γ µ = −(S µ† C µ )−1 (C̃ )µν β ν .
ν=µ
In alloy theory, this equation determines the variational wave function for an atomic
cell τµ embedded in a statistical medium defined by the vector of coefficients β for
all other cells [157, 281].
108 7 Multiple scattering theory for molecules and solids
where σ is the enclosing boundary of τin . The sign of the first term here is reversed
because the notation Wσ implies an outward normal derivative in → out, while
by convention, ∇n ψ denotes the outward normal gradient across the boundary that
encloses the region of definition of ψ. The notation Wµ will be used for a Wronskian
integral with outward normal gradients over the full bounding surface of atomic cell
τµ . The SM variational functional adds a Wronskian surface integral Z σ = 12 (ψ in +
ψ out |Wσ |ψ in − ψ out ) to the Schrödinger functional Z τ = 3 ψ ∗ (H − )ψ. This
interface surface integral vanishes if ψ in ≡σ ψ out for independent trial functions
ψ in and ψ out . Then the global trial function ψ is in the usual Hilbert space and
the variational condition δ Z = 0 is just the Schrödinger condition for stationary
energy .
Schlosser and Marcus [359] showed that for variations about such a continuous
trial function, the induced first-order variations of Z τ and Z σ exactly cancel, even if
the orbital variations are discontinuous at σ or have discontinuous normal gradients.
After integration by parts, the variation of Z τ about an exact solution is a surface
Wronskian integral
The variation of Z σ is
for each atomic cell surface σµ , summed over all adjacent cells. Although two
surface integral conditions are needed in general to eliminate both value and gradient
discontinuities, these equations can be justified because each interface σµν occurs
twice, for cells µ and ν, respectively. Expansion in a local orbital basis in each cell
gives
µ µ µ
ν ν
φ L |Wµν | φL γL − φ L γ L = 0.
ν=µ L L
For an exact solution, the external function ψ ν would be identical to the global
matching function ξ on each interface. An alternative algorithm can be based
on fitting a linear combination of basis functions in each cell to ξ , which is
uniquely determined. The VCM equations for an orbital basis at fixed energy reduce
to [282]
µ
µ µ µ
φ L |Wµ | φ L γ L − ξ = − φ L |Wµ |ξ = 0.
L
omitting L-indices and summations to simplify the notation. The first term is derived
using the general formula
(φi |Wσ |JL )(N L |Wσ |φ j ) = (φi |Wσ |N L )(JL |Wσ |φ j ),
L L
from the surface matching theorem. 2δ Z agrees with the KR variational expression
if intermediate sums over angular indices are complete [148]. Because the SM
variational principle implies stationary energy eigenvalues when the continuity
conditions are exactly satisfied, this is also true for the KKR/MST equations when
variational trial functions are continuous with continuous gradients, as for muffin-
tin and atomic-cell orbitals.
The Green function G 0 and solid-harmonic functions JL and N L are required only
for the fixed energy κ02 , usually taken to be zero. This amounts to adding a linear
energy term to the potential function in the interstitial region, justified in full-
potential theory or the ASA because the net interstitial volume is reduced to zero.
The practical effect is slower convergence of the angular momentum expansions. Al-
though the ASA violates the original KKR conditions that ensure an exact solution
of the muffin-tin model problem, the ASA model of overlapping atomic spheres is
in many ways a more satisfactory approximation to full-potential MST. In particu-
lar, in the surface integral formalism, approximating Wronskian integrals over the
surface of a Wigner–Seitz cell by the corresponding integrals over an equivalent-
volume sphere can be considered as a first approximation to a surface quadrature
scheme for polyhedral cells. The full DFT local potential function can be used within
the atomic sphere. The muffin-tin model, which assumes a constant potential in the
near-field region, has no special justification as a computational approximation for
atomic potentials.
In the usual ASA, the t-matrix is real and diagonal, but the elements have both
zeroes and poles as functions of energy. The original proposal was to simplify KKR
calculations by using constant-energy structure constants while fitting diagonal ele-
ments of t −1 in the standard secular equation, (t −1 − g)β = 0, to rational functions
of energy [8]. This also simplifies the energy dependence of muffin-tin orbitals,
which can be used as basis functions in the standard Schrödinger variational prin-
ciple for the energy eigenvalues. The resulting LMTO formalism (called LCMTO
in the original presentation) produces a linear eigenvalue problem that includes
a “combined correction” matrix representing the difference of the overlap matrix
evaluated in atomic spheres versus atomic polyhedra together with a correction for
basis functions of higher quantum numbers. LMTO energy-band calculations ob-
tain all energy eigenvalues at a given k-point from a single matrix diagonalization.
These eigenvalues are valid within an energy panel in which the rational fit to the
t −1 -matrix is sufficiently accurate.
In considering nonspherical local potentials, a simple rational fit to diagonal ele-
ments of t −1 is no longer valid. From the general definition of the t-matrix as −SC −1
in terms of the standard MST matrices, it is clear that for energy-independent solid-
harmonic functions, linear energy expansion of the local basis functions φ L (; r) is
equivalent to fitting elements of the t-matrix by a simple rational formula. In the
ASA this produces the parametrization proposed by Andersen. It was recognized at
the same time [13] that a linear energy expansion of the basis functions also solves
the practical problem in MTO theory that in the local representation χ = φ − JC
of primitive MTOs the regular solid-harmonic functions are not orthogonal to inner-
shell occupied orbitals. The practical solution is to replace the solid harmonics by
linear combinations of the energy-derivative functions denoted here by φ̇ L . In the
7.3 Energy-linearized methods 115
standard LMTO method, the basis functions are energy-independent MTOs ex-
pressed as linear combinations of functions φ L and φ̇ L computed at a fixed panel
energy [9, 384]. Basis functions of this kind are used in alternative energy-linearized
methods [9, 430] which will not be considered here.
This agrees with the form of an MST null-vector expanded in the energy-dependent
basis at k = 0 + k . Because of numerical approximations inherent in practi-
cal calculations, it is generally inconvenient to diagonalize matrix ω directly. The
116 7 Multiple scattering theory for molecules and solids
eigenvalues are complex numbers if there are any residual numerical errors. For
this reason, a standard variational principle based on Hermitian matrices is used
in the second stage of linearized calculations. The φ, φ̇ expansion must be used
with caution unless all elements of the ω-matrix are small. A large energy shift
may imply that basis functions fall outside the energy panel used in a particular
calculation, so that the interpolated functions are not accurate. So-called “ghost”
bands can occur, especially for higher angular quantum numbers, when functions
φ and φ̇ tend to become linearly dependent within a local cell [11].
where L-indices and sums are suppressed to simplify the notation. In these equa-
tions, the internal summations in Eqs. (7.10) are replaced by direct matching
across adjacent cell interfaces. If the generally complex unsymmetric matrix ω
µ
is diagonalized, with eigenvector elements c Lk , Eqs. (7.11) reduce to the VCM
7.3 Energy-linearized methods 117
µ µ
equations for ψ µ = L φ L (k )c Lk ,
µ
φ L (k )|Wµν |ψ µ − ψ ν = 0,
ν=µ
∇ 2 v = −8πρ, (7.13)
gives the well-known one-center expansion of 1/r12 [188]. It is assumed that the
local density function is subdivided into spherical harmonic components as ρ =
L ρ L . Each of these components defines a particular solution within the enclosing
sphere of the local cell. A very efficient numerical algorithm is available for solving
the radial Poisson equation [229], generalized to arbitrary -values [278]. Denoting
these primitive solutions by v̂ L , the ACO construction can be used to define a
particular solution v0L for each ρ L such that v0L = v̂ L − J ĉ in the local cell τ .
v0L matches onto an external function N ŝ on the cell surface σ . The coefficient
matrices here are defined by
Exploiting the high efficiency of the integration algorithm for the radial Poisson
equation [229], MST surface integral formalism [274, 278] reduces the computation
of multipole moments of cellular charge densities to the evaluation of Wronskian
surface integrals. Electrostatic multipole moments Q L are defined by the asymp-
totic external potential for each cell τ . For the particular choice of solid-harmonic
functions given above,
lim r →∞ v0L = N L ŝ L ,L = −8π NL Q L . (7.20)
L
A factor −2 included in the last term here compensates for the use of Rydberg units
and for the omission of the negative electronic charge in potential functions derived
from Eq. (7.14). Hence the electrostatic multipole moments of atomic cell τµ are
µ 1 µ 1 µ
µ
QL = − ŝ = JL |Wσµ | v̂ L . (7.21)
8π L L ,L 8π L
7.5.1 Definitions
For energy κ in Rydberg units, the Green function of the Helmholtz equation, with
2
This equation can be applied to a wide class of problems including the electronic
structure of vacancies, impurities, and of other localized perturbations of solids and
atomic clusters. The methods considered here make it possible to construct G(x, x )
for any potential function v(x) defined throughout the coordinate space 3 .
It is assumed that 3 is subdivided into space-filling cells. A typical cell τµ has
closed surface σµ . An external empty cell can be included to fill space for a finite
system. Solutions of the homogeneous Helmholtz equation about the origin of cell µ
define regular and irregular generalized solid-harmonic functions, denoted here by
µ µ
JL and N L , respectively. Each generalized solid-harmonic function is the product
of a radial factor and a spherical-harmonic function, defined here by a combined
index L that denotes the quantum numbers (m). Normalization is assumed to be
consistent with the Kronecker-delta relation,
Expressed in local coordinates for each atomic cell, the Helmholtz Green function
can be expanded in solid-harmonic functions,
G 0 (r, r ) = − JL (r)N L∗ (r ), r < r . (7.27)
L
To verify this expansion, consider G 0 (r, r ) for a fixed point r . For r < r , this
function is a regular solution of the Helmholtz equation and must have an expansion
122 7 Multiple scattering theory for molecules and solids
of the form
G 0 (r, r ) = JL (r)(N L |W0 |G 0 ), (7.28)
L
where the Wronskian integral is over the surface of any sphere centered at the cell
origin with radius r0 < r . For r > r , this function is a bounded solution of the
Helmholtz equation and must have an expansion of the form
G 0 (r, r ) = − N L (r)(JL |W1 |G 0 ), (7.29)
L
where the Wronskian integral is over the surface of any sphere centered at the
cell origin with radius r1 > r . If Eq. (7.22) is multiplied by N L∗ (r), and the cor-
responding null term with the Helmholtz operator applied to N L∗ is subtracted, on
integrating over the volume defined by r0 < r < r1 , between radii r0 and r1 , this
integral reduces to Wronskian surface integrals. This implies
Since r1 > r , Eq. (7.29) implies that the first integral on the right-hand side van-
ishes. The second integral gives the coefficient in Eq. (7.28), which reduces to
Eq. (7.27). Similarly, if Eq. (7.22) is multiplied by JL∗ (r) and integrated over the
volume between radii r0 and r1 ,
In this case, the second integral on the right-hand side vanishes and the first integral
gives the coefficient in Eq. (7.29). This provides a formula valid for r > r , com-
plementary to Eq. (7.27). Comparison of these two formulas shows that the linear
operator G 0 is Hermitian for real energies.
When coordinate x = Xµ + r lies in cell τµ and x = Xν + r lies in a different
µ
cell τν , a two-center expansion of the Green function is needed. Since N L (x − Xµ )
is a regular solution of the Helmholtz equation in a sphere of radius |Xµ − Xν |
about the origin of a displaced cell ν,
µ µ
νµ
N L (r + Xν − Xµ ) = N L (r )δµν − JLν (r )g L ,L , (7.32)
L
νµ
for r < |Xµ − Xν |. The expansion coefficients g L ,L are structure constants at the
given energy. This matrix is Hermitian for real energies. Using the above expansion,
the Helmholtz Green function can be written in the form,
µ µν
µ
G 0 (Xµ + r, Xν + r ) = JL (r)g L L JLν∗ (r ) − JL (r)N Lν∗ (r )δµν , (7.33)
L L L
7.5 Green functions 123
where r < r within any single cell µ. The two-center expansion converges in gen-
eral for muffin-tin geometry, when the two coordinate values lie within nonover-
lapping spheres. It is assumed here that cell coordinate origins are chosen to satisfy
this geometric condition for some neighborhood of each cell origin.
These definitions can be generalized to give an expansion of the Green function
G(x, x ) in terms of local solutions of the Schrödinger equation in each cell τ . Local
regular solutions φ L and local irregular solutions ζ L can be defined by matching to
functions JL and N L , respectively, on an infinitesimal sphere about the cell origin,
or, alternatively, by requiring similar conditions on a local cell boundary σ . This
construction imposes the Wronskian condition
(ζ L |Wσ |φ L ) = δ L ,L . (7.34)
These integrals vanish between any two φ or two ζ functions. Green’s theorem
implies that this Kronecker-delta formula is valid for Wronskian integrals over
any closed surface that encloses the local cell center and excludes all other singular
points of the potential function. If singularities of the potential function occur only at
cell centers, these functions can be assumed to be have regular extensions out to the
nearest neighbor cell center, but in general no further. For this reason, the primitive
function ζ must be modified by the addition of a sum of regular functions φ in its
cell of origin to define a modified irregular function ζ̃ that is regular everywhere
except at the defining cell center. With these definitions [53, 441, 444], by analogy
to the Helmholtz Green function,
G(r, r ) = − φ L (r)ζ̃ ∗L (r ), r < r , (7.35)
L
µ µ
νµ
ζ̃ L (r + Xν − Xµ ) = ζ L (r )δµν − φ Lν (r )G L ,L , (7.36)
L
valid for r < r within a single cell, or for coordinates in two nonoverlapping
µν
spheres. The coefficients G L ,L define the structural matrix of the Schrödinger Green
function. For ν = µ, ζ̃ µ − ζ µ is a sum of regular functions φ µ with coefficients
µµ
−G L ,L .
The definition of irregular functions ζ depends on boundary conditions imposed
on the Green function. Any variation of the set of functions ζ defined by adding
linear combinations of the regular functions φ with coefficients that constitute an
Hermitian matrix simply moves the corresponding term between the two summa-
tions in Eq. (7.37), leaving the net sum invariant, and preserving the Kronecker-delta
124 7 Multiple scattering theory for molecules and solids
conditions of Eq. (7.34). The present derivations are covariant with respect to the
group of such transformations. For particular applications, the representation can be
chosen for the greatest convenience, subject only to the Kronecker-delta condition.
which defines G(z) = (z − H)−1 as the resolvent operator. G(z) is an analytic op-
erator function in the upper half of the complex z-plane. The singularities are poles
or a continuum branch cut corresponding to energy eigenvalues. If H is Hermitian,
the eigenvalues lie on the real energy axis. The corresponding orthonormal set of
eigenfunctions ψi defines a spectral representation
G(z; x, x ) = ψi (x)(z − i )−1 ψi∗ (x ). (7.39)
i
This formula indicates that the residue of G(z) at an eigenvalue pole is just the
density matrix, in the coordinate representation, summed over all eigenfunctions
with this eigenvalue. This residue can be extracted conveniently by using Dirac’s
formula after displacing each eigenvalue into the lower complex plane,
1 1
lim δ→0 =P − iπ δ(z − i ). (7.40)
z − i + iδ z − i
Specializing to the electronic density for electrons of one spin in the local coordi-
nates of a particular cell, the local density per unit energy is
1
n(; r) = − G(; r, r). (7.41)
π
Summing over all eigenenergies below a chemical potential or Fermi level µ, the
local density function is defined by a contour integral passing above all poles on
the real axis,
1 µ
ρ(r) = − G(z; r, r) dz. (7.42)
π −∞
In methodology based on computation of the Green function, this formula re-
places the usual sum over eigenfunctions. When integrated over coordinate space,
7.5 Green functions 125
These formulas must be summed over a spin index to give the corresponding total
densities. In multiple scattering theory, the Green function as given by Eq. (7.37) is
subdivided into terms valid in separate atomic cells. The elements with r = r are
single sums over the cells. Hence the density of states of given spin in a particular
cell τµ is
µ 1
n () = − G µ (; r, r), (7.44)
π τµ
again omitting L-indices and summations. This is a Dyson equation in the form
of a matrix representation G = g + gt G of the Lippmann–Schwinger equation,
Eq. (7.24). This equation can also be derived by using the idea of atomic-cell
orbitals (ACO), generalized to irregular functions. An irregular ACO is defined in
a reference cell τλ as the canonical function ζ λ . In all other cells τµ it is the local
expansion of N λ , expressed in terms of the structure constants g µλ . Similarly, a
regular ACO takes the form φ − J in its indexed cell and Jgt in all other cells. The
ACO form of ζ̃ λ in cell τµ is
ζ̃ λ = ζ µ δ µλ − φ µ G µλ + J µ (δ µν − g µν t ν )G νλ − g µλ , (7.48)
ν
the form,
(N µ |Wµν |(N µ δ µλ − J µ g µλ ) − (N ν δ νλ − J ν g νλ )) = 0, (7.49)
ν=µ
While these equations are consistent with the GFCM, their relationship to the
VCM variational equations suggests that the alternative formalism defined by the
JVCM might be more appropriate. In the global solution sought in the present case,
the irregular term is fixed, and all variational increments are linear combinations
of regular functions JL . Hence variation of coefficients gives the JVCM Eqs. (7.7).
A similar situation occurs for the Poisson equation, since a particular solution of
the inhomogeneous equation is fixed, and all incremental functions are regular
solid harmonics. Recent calculations comparing these methods for the Poisson
equation [443] show improved convergence in the JVCM formalism. It was shown
that extension of the Schlosser–Marcus variational principle to the Poisson equation
implies the JVCM equations [443]. In the present case, Eqs. (7.7) take the form
(J µ |Wµν |(N µ δ µλ − J µ g µλ ) − (N ν δ νλ − J ν g νλ )) = 0, (7.51)
ν=µ
for all values of µ and λ. These equations determine the structural matrix G for
any set of space-filling cells.
Equations (7.53) simplify if the system has point-group or translational sym-
metry. Then local basis functions of equivalent cells are related to those of a smaller
number of generating cells by phase factors and elementary rotations. Considering
only translational symmetry, it is convenient to index cells in a reference transla-
tional cell by indices µ, etc. and to index translated equivalent cells by the corre-
sponding indices µ̃, etc., such that the displacement of cell τµ̃ relative to cell τµ
128 7 Multiple scattering theory for molecules and solids
is a lattice translational vector dµ̃µ . Using this notation, Eqs. (7.53) imply separate
equations for each k-vector in the reduced Brillouin zone, in the form
µ ν̃ νλ
ν̃=µ (φ |Wµν̃ |φ ) exp(ik · dν̃ν )G (k) =
µ ν̃ νλ
ν̃=µ (φ |Wµν̃ |ζ ) exp(ik · dν̃ν )δ + δ µλ , (7.54)
for all values of µ and λ that index cells in the reference translational cell. The
corresponding equations derived from Eq. (7.50) have been verified by comput-
ing structure constants for an fcc lattice along the ( − K line in the reduced
Brillouin zone, checked against a published structure-constant program STR [384].
Equations (7.53) have not been tested, but are expected to show improved conver-
gence because they are variationally correct for this problem.
8
Variational methods for continuum states
129
130 8 Variational methods for continuum states
When substituted into the first equation, this gives an effective multichannel
Schrödinger equation of the form
{ H̃ − E}! P = {P(H − E)P − PHQ{Q(H − E)Q}−1 QHP}! P = 0. (8.3)
132 8 Variational methods for continuum states
Equation (8.3) provides a common basis for the major computational methods used
in low-energy electron scattering theory. The terms containing {Q(H − E)Q}−1 in
this equation describe polarization response and other correlation effects in terms
of a matrix optical potential. Singularities due to this energy denominator cancel
out exactly [264, 265, 270].
Effective one-electron equations for the channel orbital functions can be obtained
either by evaluating orbital functional derivatives of the variational functional
or more directly by projecting Eq. (8.3) onto the individual target states # p . With
appropriate normalizing factors, (# p |!s ) = ψ ps . Equations for the radial channel
functions f ps (r ) are obtained by projecting onto spherical harmonics and elemen-
tary spin functions. The matrix operator acting on channel orbitals is
m̂ pq = (# p | H̃ − E|A#q ). (8.4)
The projection of this operator onto spherical harmonics and spin functions defines
a radial operator m pq = ĝ pq − p δ pq , where p = E − E p .
α0 = I, α1 = K ,
which defines the reactance matrix K = α1 (α0 )−1 . For exact solutions of the
projected Schrödinger equation, the open-channel K-matrix is real and sym-
metric [257, 302]. It can be diagonalized by an orthogonal transformation to
define eigenchannels, defined by the eigenvectors. The eigenvalues are tan ησ , for
σ = 1, . . . , n o , which defines the eigenphases ησ up to a multiple of π radians. The
scattering matrix S is defined by [257]
S = (I + i K )(I − i K )−1 ,
The S-matrix is unitary and symmetric, while the T-matrix is symmetric. This
particular definition of the T-matrix reduces for scattering by a central potential to
the phase-shift factor in the scattering amplitude,
1 ∞
f (θ ) = (2 + 1) exp(iη ) sin η P (cos θ ),
k =0
The total cross section for unpolarized scattering is obtained by summing σqp over
degenerate final states and by averaging over initial states.
satisfactory way has been found to include the continuum of unbound target states
required for completeness of this expansion.
Since only a relatively small number of coupled equations can be solved in
practice, the target states # p must be selected carefully. Virtual target excitation
into the ionization continuum must be approximated by inclusion of closed-channel
pseudostates that cannot be target eigenstates but have the character of wave packets
in the continuum. Target polarization response is treated by including polarization
pseudostates #γ ( p) .
If only a single target state is included, itself represented as a single-determinant
model state " p , the formalism reduces to the static-exchange (SE) model. Because
of the complexity of more detailed scattering models, much work in fixed-nuclei
electron scattering by molecules has been carried out at this level of approximation
[215]. This model is generally unsatisfactory for low-energy scattering, and must
be augmented by a valid approximation to the dipole polarization response of the
target system [270, 135]. Models that augment static exchange with polarization
response are denoted by the acronym SEP. It will be shown below that a model with
this computational structure can be derived from formally exact orbital functional
theory (OFT), if the correlation potential operator v̂c is approximated in terms of
polarization pseudostates [290].
f ps (r ) = F0 p α0 ps + F1 p α1 ps (8.6)
the variational approximation to the (N +1)-electron wave function takes the form
!s = A# p φi p + "µ cµ αi ps .
ip
ip µ
ip
In the matrix variational method, the coefficients cµ are determined separately
for each set of indices i, p from the matrix equations, for all µ,
"µ |H − E|A#q φ jq + "ν cνjq = 0, (8.7)
ν
which imply the variational condition ("µ |H − E|!s ) = 0. For arbitrary values of
the coefficients αi ps these equations follow from the variational condition
∂st
i p∗
= 0,
∂cµ
for all µ, i, p. The variational functional st = (!s |H − E|!t ) here is an n o × n o
matrix.
When Eq. (8.7) is satisfied, the variational functional becomes an explicit
quadratic function of the coefficients αi ps , which can be assumed to be real numbers.
Thus
pq
st = αi ps m i j α jqt ,
ip jq
for all indices i, p. Each such solution defines a 2n o element column vector of
the rectangular matrix of coefficients α. In a matrix notation, suppressing channel
indices, the variational functional is
† †
= α † mα = α0 (m 00 α0 + m 01 α1 ) + α1 (m 10 α0 + m 11 α1 ).
In this notation, Eq. (8.8) is
m 00 m 01 α0
mα = = 0.
m 10 m 11 α1
matrices to give stationary expressions for these matrices. This logic produces
multichannel versions of the variational principle originally derived by Hulthén
[175, 176], Kohn [202], and Kato [194, 195]. As applied to the K-matrix, Kohn’s
variational principle can be derived by considering variations about an exact solution
that maintain the canonical form α0 = I, α1 = K . Then δα0 = 0 and δα1 = δ K ,
which implies δ = 12 δ K . This implies that the Kohn functional
[K ] = K − 2
1
δ = δ K † (m 10 + m 11 K t ) + (m 10 + m 11 K t )† δ K + δ K .
2
This is a real symmetric matrix by construction, and remains valid for an exact
variational solution, when = 0. Hence this derivation proves that an exact open-
channel K-matrix is real and symmetric.
The difficulty with this formula is that there is no guarantee that the auxiliary
matrix m 11 is not singular [264, 265]. At an energy for which an eigenvalue of m 11
vanishes, K t is not defined unless all columns of the matrix m 10 are orthogonal
to the null eigenvector. Specific examples [265] show that m 10 does not have this
property in general for otherwise valid variational trial wave functions. The number
of null values of m 11 increases with the number of basis functions, leading to energy
and basis-dependent anomalies in the Kohn formalism [369, 370]. In an alternative
version of the theory [176, 350], variations are restricted to α0 = K −1 and α1 = I .
This implies that anomalies occur at null points of m 00 , which in general do not co-
incide with those of m 11 . This fact can be exploited to generate various anomaly-free
variants of the Kohn theory [265]. Especially in the matrix variational formalism,
this method has been applied to many calculations of electron–atom scattering of
high accuracy [270]. This method was used to compute e− –He scattering cross
sections that established a calibration standard for subsequent experimental work
[268, 269].
8.2 Kohn variational theory 139
limitation of the matrix variational method [270], in that analytic formulas are
available only for specialized forms of the long-range potentials that can occur,
especially for molecules.
(h − )wi (r ) = 0
Acting on any function F(r ) that vanishes for r → ∞ and is regular at the origin
r → 0,
∞ r
G F(r ) = 2 w0 (r ) w1 (r ) + w1 (r ) w0 (r ) F(r ) dr .
r 0
2 2
are used to construct the principal value Green function. This Green function im-
poses the asymptotic form f (r ) ∼ w0 (r ) + w1 (r ) tan η, where
∞
tan η = −2 w0 (r )v(r ) f (r ) dr .
0
Other asymptotic forms consistent with unit Wronskian define different but equally
valid Green functions, with different values of the asymptotic coefficient of w1 . In
particular, if w1 ∼ k − 2 exp i(kr − 12 π ), this determines the outgoing-wave Green
1
Coulomb or static dipole potentials should be included in the model potential whose
solutions w0 , w1 are used to construct the Green function.
Specializing the present derivation to the principal value Green function, the
unsymmetrical expression tan η = −2(w0 |v| f ) is exact for an exact solution of
the Lippmann–Schwinger equation, but it is not stationary with respect to infinites-
imal variations about such a solution. Since w0 = f + Gv f for such a solution,
this can be substituted into the unsymmetrical formula to give an alternative, sym-
metrical expression tan η = −2( f |v + vGv| f ), which is also not stationary.
However, these expressions can be combined to define the Schwinger functional
which is stationary for variations of the trial function f about an exact solution.
The variation δ[tan η] induced by an infinitesimal variation of f about an exact
solution is
(δ f |v|w0 ) (δ f |v + vGv| f )
− + hc [tan η]
( f |v|w0 ) ( f |v + vGv| f )
= −2(δ f v|w0 − f − Gv f ) + hc.
This implies that [tan η] is stationary if and only if f , in the range of v, satisfies
the Lippmann–Schwinger equation.
The Schwinger functional has several remarkable properties. It combines three
formulas that may have different numerical values:
For [tan η] to be stationary with respect to variations of the coefficients (a| f ), this
8.3 Schwinger variational theory 143
Combining these equations, all three expressions for tan η give the same result for
the stationary functional,
[tan η] = (w0 |a)vac (v + vGv)−1 cd vdb (b|w0 ).
a,b c,d
where
m pq = (h − ) pq + v pq .
It is assumed that all closed-channel components have been eliminated by the
partitioning transformations that determine Eq. (8.4), from which the radial operator
m pq is derived. Equations of this form can be solved using the Green function of
the model problem
(h − ) ps wisq (r ) = 0.
s
144 8 Variational methods for continuum states
The multichannel model functions w0sq and w1sq are real-valued continuum solu-
tions of the model problem, in a K-matrix formalism. They are normalized to the
matrix Wronskian condition
(wisp w jsq − wisp w jsq ) = δ pq (δi0 δ j1 − δi1 δ j0 ).
s
q q r 0
such that
(h − ) pp G σp q Fq (r ) = F p (r ).
σ p,q
Using the matrix Wronskian condition, the proof is the same as for the single-
channel problem.
This multichannel matrix Green function determines a multichannel Lippmann–
Schwinger equation
u pτ = w0 pτ − G σpp v p q u qτ
σ p ,q
which implies that the incremental K-matrix in the model eigenchannel representa-
tion is K = −2(w0 |v|u), in a matrix notation. The Hermitian conjugate relation
is K = −2(u|v|w0 ).
8.3 Schwinger variational theory 145
K = −2(u|v + vGv|u).
[K ]σ τ =
−2 (w0 pσ |v|η p a )[(η|v + vGv|η)]−1
p a,q b (ηq b |v|w0qτ ).
a,b p, p q,q
The same result is obtained by substituting a formal solution for u into any of the
three approximate expressions for K combined in the Schwinger functional matrix.
which become
m pq u qs (r ) = η pa λas
p
. (8.12)
q a
The indices refer to both open and closed channels, which do not have to be dis-
tinguished at the finite radius r1 . As defined for solutions of radial close-coupling
equations, R pq is a real symmetric matrix. Extension of the theory and definitions to
more general geometry will be considered below. For electron–molecule scattering,
a molecular center is chosen somewhat arbitrarily such that the electronic charge
is approximately enclosed in some smaller sphere of radius r0 , while bound–free
exchange can be neglected outside the larger R-matrix radius r1 .
The theory of the R-matrix was developed in nuclear physics. As usually pre-
sented, the theory makes use of a Green function to relate value and slope of the
radial channel orbitals at r1 , expanding these functions for r < r1 as linear com-
binations of basis functions that satisfy fixed boundary conditions at r1 . The true
logarithmic derivative (or reciprocal of the R-matrix in multichannel formalism)
148 8 Variational methods for continuum states
is computed from Green’s theorem, despite the use of basis functions whose log-
arithmic derivatives at r1 have a fixed but arbitrary value. Because of the inherent
discontinuity of the boundary derivative this expansion tends to converge slowly,
and requires correction by an approximate method due to Buttle [54].
In nuclear physics, the specifically nuclear interaction is of short range, so that
full scattering information is obtained by matching the R-matrix at r1 to asymptotic
external wave functions. The method has been extended to electron–atom scatter-
ing, where long-range potentials are important, by combining basis expansion of
the channel orbital functions within r1 with explicit numerical solution of close-
coupling equations outside [50, 47]. This method makes it possible to process
algebraic equations defined by matrix elements of nonlocal operators within r1 ,
while solving the simple asymptotic close-coupling equations, without exchange,
outside r1 . This requires r1 to be large enough that exchange can be neglected
and that the nonlocal optical potential can be approximated by an asymptotic local
potential.
The R-matrix can be matched at r1 to external channel orbitals, solutions in
principle of external close-coupling equations, to determine scattering matrices.
Radial channel orbital vectors, of standard asymptotic form for the K-matrix,
− 12 1
w0 pq (r ) ∼ k p sin k p r − p π δ pq ,
2
− 12 1
w1 pq (r ) ∼ k p cos k p r − p π δ pq ,
2
are defined by integrating inwards for r > r1 . The asymptotic forms must be modi-
fied appropriately for Coulomb or fixed-dipole scattering. The function to be fitted
at r1 is
f ps (r1 ) = [w0 pq (r1 )δqs + w1 pq (r1 )K qs ].
q
The R-matrix relation between function value and gradient can be solved for the
K-matrix,
−1
K st = − w1 ps (r1 ) − r1 R pq (r1 )w1qs (r1 )
p q
sp
× w0 pt (r1 ) − r1 R pq (r1 )w0qt (r1 ) .
q pt
As written, these equations refer to open channels only. When external closed
channels are considered, an external closed channel orbital that vanishes as r → ∞
must be included for each such channel. The indices p, q, s run over all channels,
8.4 Variational R-matrix theory 149
but the index t refers to open channels only. K st is a rectangular matrix, whose
open-channel submatrix can be shown to be real and symmetric.
The theory of the R-matrix can be understood most clearly in a variational for-
mulation. The essential derivation for a single channel was given by Kohn [202],
as a variational principle for the radial logarithmic derivative. If h is the radial
Hamiltonian operator, the Schrödinger variational functional is
∞
= f (h − ) f dr.
0
If f (r ) is an exact solution for r > r1 , this reduces the functional to the finite integral
considered by Kohn,
r1 r1
1 2 1 2 1
λ = f (h − ) f dr = ( f ) + v(r ) − k f 2 dr − λ f 2 (r1 ),
0 0 2 2 2
obtained after integration by parts. The parameter λ is the logarithmic derivative
f / f = (r1 R)−1 , defining the parameter R as the single-channel R-matrix.
For an infinitesimal variation δ f , with fixed λ, the variation of λ is
r1
δλ = 2 δ f (h − ) f dr + δ f (r1 )[ f (r1 ) − λ f (r1 )].
0
where
r1
1 1
Aab = ηa ηb + ηa v(r ) − k 2 ηb dr. (8.14)
0 2 2
The parameter λ is uniquely determined as a consistency condition [202]. From
Eq. (8.13),
1 −1
ca = λ (A )ab ηb (r1 ) f (r1 ),
2 b
so that
1
f (r1 ) = ηa (r1 )ca = λ ηa (r1 )(A−1 )ab ηb (r1 ) f (r1 ).
a 2 a,b
150 8 Variational methods for continuum states
Then, by construction,
1
−1
f ps (r1 ) = R ps r12 , f ps (r1 ) = δ ps r1 2 ,
integrating the kinetic energy term by parts. ∇n ψ here is the outward normal gradient
on σ of the trial function ψ. The residual volume integral A, incorporating a Bloch-
modified [30] kinetic energy term, is
1 ∗ ∗
A= ∇ψ · ∇ψ + ψ (v − )ψ .
τ 2
This will be symbolized here by ψ(σ ) = R∇n ψ(σ ) = Rξ (σ ). If the normal gradi-
ent ξ is specified, this defines a classical Neumann boundary condition on σ , which
determines a unique solution of the Schrödinger equation in the enclosed volume
τ . The value of the boundary integral is
1
A1 = ψ ∗ ξ.
2 σ
The functional must vanish for an exact solution of the Schrödinger equation
in τ . This implies that A1 and A are equal, and because A is real,
A = A1 = A∗1 .
and
1
δ A1 = δψ ∗ ξ.
2 σ
which is real for any trial function. From the variations of A1 and A,
−1 1 ∗ −1 1 ∗ ∗
δ[A] = [A] A1 δψ ξ − A δψ ∇n ψ + δψ (h − )ψ + cc .
2 σ 2 σ τ
For variations about an exact solution, for which A = A1 = A∗1 , this reduces to
1 ∗ ∗
δ[A] = δψ (ξ − ∇n ψ) − δψ (h − )ψ + cc .
2 σ τ
For variations of ψ that are unconstrained throughout τ and on σ , this implies that
[A] is stationary if and only if (h − )ψ = 0 in τ and ∇n ψ = ξ on σ .
When ψ(x) = a φa (x)ca is expressed as a sum of basis functions, variations
δψ are driven by variations of the coefficients ca . If δ[A] vanishes, this implies
1
∗ ∗
φa (h − )φb cb = φ ξ− ∇n φb cb .
b τ 2 σ a b
Explicitly,
1 −1
cb = [A ]ba φa∗ ξ, all b.
2 a σ
When these values are substituted into the definitions of A, A1 , and [A], all of these
quantities are equal to
1
[A] = ξ ∗ (1)φa (1)[A−1 ]ab φ ∗ (2)b ξ (2).
4 a,b σ1 σ2
Because [A] is stationary and the surface function ξ is arbitrary, the operator
R must itself be stationary. It is evidently real and symmetric by construction.
The basis functions φa (x) are required only to be linearly independent in τ . En-
forcing orthonormality through a fixed boundary condition on σ imposes a lack
of completeness with respect to determining the R-operator. The factor 1/r1 is
omitted from the definition given above, making the operator R dimensionally a
length. This cannot be avoided in general geometry, since no unique radius is de-
fined for a nonspherical surface. The usual R-matrix in spherical geometry is the
surface-harmonic representation of the R-operator as defined here, divided by the
dimensional scale factor r1 .
Bound–free correlation
For an N-electron target state, E c is a sum of pair-correlation energies. In a two-
electron system, such as atomic He, a major contribution to the correlation energy
8.4 Variational R-matrix theory 159
where all integrals are evaluated for r ≤ r1 . The second line of this general formula
is a correction to target correlation energy due to removal of φκ from the unoccupied
set. The first line can be shown to include a multipole polarization potential.
To model the effect of close-coupling equations including a closed pseudostate
channel, suppose that a pseudostate orbital φ p j represents the first-order pertur-
bation of a particular target orbital φ j in a multipole field of order λ, and that
p
the pseudostate excitation energy is E j − E 0 . Because φ p j is selected to interact
strongly within the inner radius r0 , the residual orthogonalized set of unoccupied
orbitals φq remain approximately complete in r0 ≤ r ≤ r1 for functions orthogonal
to occupied orbitals. The close-coupling equations imply that a closed-channel or-
bital φqκ coupled to φ p j is approximately of the form of a slowly decaying function
times the open-channel orbital φκ . Hence the average energy of the two-electron
p
virtual excitations j, κ → p j , q can be approximated by E j − E 0 . This argument
justifies a closure approximation, and the first line of Eq. (8.16) becomes
(κ|v̂c |κ) nj (1 − n q )(κ j|ū|q p j )(q p j |c̄|κ j)
j q
p −1
− n j κ|( j|u| p j ) E j − E 0 ( p j |u| j)|κ . (8.17)
j
160 8 Variational methods for continuum states
SE SEP
k (au) Present Ref. [383] Present Ref. [383]
=0
0.3 2.7049 2.7037 2.7424 2.7459
0.5 2.4357 2.4325 2.4812 2.4836
0.7 2.1943 2.1963 2.2433 2.2518
0.9 1.9836 1.9808 2.0339 2.0426
=1
0.3 0.0140 0.0108 0.0340 0.0297
0.5 0.0447 0.0426 0.0866 0.0847
0.7 0.0980 0.0947 0.1629 0.1567
0.9 0.1570 0.1552 0.2366 0.2311
161
162 9 Electron-impact rovibrational excitation of molecules
transfer of kinetic energy between electronic and nuclear motion can be neglected.
1
In the expansion parameter (m/M) 4 , where M is a typical nuclear mass, elec-
tronic energy is of zeroeth order, vibrational energy of second order, and rotational
energy of fourth order. Far from neglecting these Born–Oppenheimer corrections,
a valid theory of electron-impact rovibrational excitation must compute the scat-
tering effects of such terms with quantitative accuracy, because the phenomenon
of interest is precisely such an energy transfer between electronic and nuclear
motion.
Rotational level spacing is in fact so small for most molecules that it can be
resolved only by precise spectroscopic techniques. Typical electron scattering data
does not have sufficient energy resolution, and most published data is for rotationally
averaged scattering cross sections. Special methods based on an adiabatic theory
of slowly moving nuclei are valid under these circumstances, not dependent on
variational theory and methods. However, near excitation or scattering thresholds,
rotational analysis may play an essential role in interpreting complex scattering
phenomena due to dynamical long-range potentials, such as the dipole potential
for heteropolar diatomic molecules, quenched by molecular rotation. For these
reasons, the present discussion of computational methodology will concentrate on
threshold effects characteristic of low-energy electron scattering and on vibrational
excitation. The phenomenon of dissociative attachment, when an incident electron is
captured by a target molecule, which then dissociates to produce an ionic fragment,
is a particular example of such energy transfer, and will be considered here as an
application of rovibrational threshold theory.
An electron scattering resonance, in abstract but general mathematical terms, is
characterized by a singularity of the scattering matrix when analytically continued to
a complex energy value res − 2i γ for relatively small width parameter γ > 0 [270].
In simple effective-potential models, a so-called shape resonance occurs when an
electronic energy level in a potential well is degenerate with an energy continuum
outside a potential barrier that confines the well. An electron initially bound in
the state corresponding to this energy level leaks out through the potential barrier
with a time constant h/γ . Such poles of the analytically continued scattering matrix
can also cause significant energy-dependent Wigner cusp or rounded step structures
[427, 32] in scattering cross sections just below or just above the energy threshold at
which a continuum becomes energetically accessible. This behavior is characteristic
of long-range potentials and does not require a well-defined potential barrier [270].
A typical phenomenon is a virtual state, which occurs in simple potential models
at the threshold for a continuum with orbital angular momentum = 0, hence no
centrifugal barrier.
Resonance and threshold structures in fixed-nuclei electron-scattering cross sec-
tions are replicated, with characteristic energy shifts, in rovibrational excitation
9.1 The local complex-potential (LCP) model 163
cross sections. One of the most striking examples is provided by the prominent
multiple peaks observed in both the electronically elastic and vibrational excita-
tion cross sections for electron scattering by the N2 molecule [368, 143, 95, 215].
These peaks are associated with a fixed-nuclei resonance at approximately 2 eV,
but their width arises from electron–vibrational coupling and their spacing does
not correspond precisely to the vibrational level structure of either the neutral
target molecule or of a vibrating transient negative ion. A quantitative theory of
these excitation structures has been developed through several advances in both
formalism and computational technique, ultimately based on variational theory of
the interacting electron–vibrational system. Another striking example is the ob-
servation of excitation peaks associated with each successive vibrational excita-
tion threshold in electron scattering by dipolar molecules such as HF, HCl, and
HBr [346]. Because of the electric dipole moment of such target molecules, if
low-energy scattering theory is to be relevant it must provide a detailed analy-
sis of rotational excitation and of rotational screening of the long-range dipole
potential.
potential for large internuclear separation. In this case, an outgoing wave packet
becomes a true bound state of the negative ion after crossing outwards, and disso-
ciates without returning. Excitation and decay processes inside the crossing point
are the same as in the boomerang model. An incident electron, captured into ex-
cited ionic vibrational levels, creates a vibronic wave packet that moves outwards
while decaying into vibrational excitations, but eventually simply crosses outwards
as a molecular negative ion that dissociates into stable fragments. Applied to the
dissociative attachment of H2 [16, 19, 415], this parametrized model accurately
reproduces available experimental data.
A very significant practical problem with the LCP model is the lack of a unique
and well-defined theory of the required resonance state when its decay width be-
comes large. In the case of H−2 , serious discrepancies between various parametrized
or computed values of the energy and width of the 2 'u+ shape resonance in the LCP
model of dissociative attachment indicate that a more fundamental, first-principles
theory is needed to give fully convincing results. The difficulty is that as the
decay width becomes large, the concept of a well-defined LCP model becomes
questionable, since the scattering resonance fades into the background scattering
continuum [271]. Analytic theory implies an energy shift which may qualitatively
change the character of the assumed resonance potential curve [59, 82].
The coupled Schrödinger equations can be projected onto the φd · · · φd∗ subspace
by Feshbach partitioning, giving an equation for the coefficient function χd (q)
in the component φd χd (q) of the total wave function. The effective Hamiltonian
in this equation is T̂ n + Vd (q) + V̂opt , which contains an optical potential that is
nonlocal in the q-space. This operator is defined by its kernel in the φd · · · φd∗
subspace,
−1
1 ∗
vopt (q, q ; E) = k dk d)k Vdk (q) E − k 2 − T̂ n − V0 (q) + iη Vdk (q ).
2
Here kv2 = 2(E − E v ) for the bound or continuum vibrational state indexed by v.
Thus the Feshbach formalism implies energy-dependent, nonlocal energy-shift and
width functions for a resonance.
Neglecting nonresonant scattering, the resonant contribution to the transition
matrix is [82]
∗
T (k f , v f ; ki , vi ) = v f Vdk f
G(E)Vdki
vi ,
since the projected equation for the coefficient function χd (q) implies χd (q) =
G(E)Vdki χvi . This projected equation is
This matrix can be computed from the general variational formula derived in
Chapter 8, using a complete set of vibronic basis functions
It is assumed that target states # p are indexed for each value of q such that a smooth
diabatic energy function E p (q) is defined. This requires careful analysis of avoided
crossings. The functions χα should be a complete set of vibrational functions for
the target potential V p = E p , including functions that represent the vibrational con-
tinuum. All vibrational basis functions are truncated at q = qd , without restricting
their boundary values. The radial functions ψα should be complete for r ≤ a.
Free boundary conditions are not allowed in the formulation of the theory
by Schneider et al. [363], based on the nonvariational theory of Wigner and
Eisenbud [428]. Specific boundary conditions are imposed using a Bloch opera-
tor. This determines boundary conditions correctly at energy poles of the R-matrix
determinant, but requires a Buttle correction [54] for energy values between such
poles [363]. This becomes problematic for the internuclear coordinate, because the
physical model of the dissociating state is a complex potential function for q ≤ qd ,
so that fixed boundary conditions imply complex energy eigenvalues. Nevertheless,
in the usual case that R-matrix poles are associated with homogeneous Neumann
boundary conditions on the R-matrix boundary, the Wigner–Eisenbud theory and
variational R-matrix theory derive the same equations for the vibronic R-matrix.
The projection integrals on the electronic wall are
qd
( pLv|α) = (A# p χv |!α )r =a = χv∗ (q)ψα (q; a)χα (q) dq.
0
This is the reciprocal of the logarithmic derivative of the wave function χd (q) in the
dissociation channel, for q = qd . At given total energy E these R-matrix elements
are matched to external scattering wave functions by linear equations that determine
the full scattering matrix for all direct and inverse processes involving nuclear
motion and vibrational excitation. Because the vibronic R-matrix is Hermitian by
construction (real and symmetric by appropriate choice of basis functions), the
vibronic S-matrix is unitary.
Schneider et al. [363] use Born–Oppenheimer vibronic basis functions as indi-
cated above, and neglect Born–Oppenheimer corrections determined by the inter-
nuclear momentum operator acting on the electronic wave function. Radial basis
functions ψk (q; r ) correspond to R-matrix pole states, whose energy values E k (q)
define an indexed vibrational potential. Vibrational basis functions χkµ (q) are com-
puted as eigenfunctions of the corresponding Hamiltonian. Since resonance states
are not treated by projection, the method depends on effective completeness of the
double expansion in electronic and vibrational eigenfunctions. The first application
of this method, to the multipeaked vibrational excitation structure observed in e−N2
scattering, was remarkably successful [364], in much closer agreement with exper-
iment than were comparable calculations using the hybrid ADN close-coupling
formalism [61]. Electronic wave functions with fixed boundary conditions at r = a
were used for the four lowest R-matrix pole states, and vibrational wave func-
tions were computed without considering the vibrational continuum. Boundary
172 9 Electron-impact rovibrational excitation of molecules
values on the cap surface are not relevant since dissociation is not involved. A
Buttle correction for the electronic basis was computed using adiabatic theory.
These calculations were later extended up to 30 eV [138] scattering energy, in-
cluding differential cross sections for elastic scattering and vibrational excitation.
The original method was modified to use a fixed electronic basis set for all inter-
nuclear distances, in order to mitigate problems arising from avoided crossings of
R-matrix pole state potential curves. Vibrational wave functions were represented
in a basis of shifted Legendre polynomials. These procedures were used for cal-
culations of integral and differential cross sections in e−HF scattering [250] and
for similar calculations on HCl [251], both examples of threshold excitation peaks
due to the molecular dipole moment. Because results computed with this method,
based on a fixed boundary condition at q = qd , depend strongly on the choice of qd ,
a modified theory has been proposed in which the evidently successful R-matrix
theory of vibrational excitation is combined with resonance-state theory for nu-
clear motion [103]. The approximation of neglecting derivatives of electronic wave
functions with respect to internuclear coordinates appears to be satisfactory in all
of these applications [362].
Evaluated at specified r = a, the phase matrix " of the R-matrix has analogous
properties. For real energies, " has a monotonically increasing trace, which can be
made continuous by suitable choice of the branch of each multivalued eigenphase
function, adding or subtracting integral multiples of π at each energy value. A
point of maximum slope defines a “precursor resonance” [284] independently of
any specific physical model, which corresponds to a pole of the analytically con-
tinued S-matrix as the R-matrix boundary is increased. The proposed methodology
uses time-delay analysis [389] to separate a given fixed-nuclei phase matrix into
a rapidly varying resonant part and a slowly varying background part. The rovi-
bronic phase matrix is constructed by applying resonance theory to the resonant
phase matrix only, treating the nonresonant background part by energy-modified
adiabatic theory. By separating out rapid variations of " this methodology reduces
the completeness requirement from the double expansion inherent in the method of
Schneider et al. [363] to a single expansion for a well-defined precursor resonance
state. Vibrational completeness is obtained by using a special basis of “spline-delta”
functions that make no distinction between bound and continuum vibrational wave
functions [295].
continued U -matrix as defined above. This implies an analytic formula for the
resonant phase matrix "1 , which will be considered in more detail below. Given
the phase matrix ", this construction of "1 defines a background matrix "0 by
subtraction, such that " = "0 + "1 . This procedure can be repeated for several
neighboring precursor resonances if necessary. Since " is unchanged, no informa-
tion is lost.
1
where vv = [(E − E v )(E − E v )] 2 . The geometric mean is appropriate to the
dimensionless vibronic phase matrix for general short-range potentials [236], as
in the earlier EMA theory [267].
The EMAP method has been used in ab initio calculations of near-threshold rota-
tional and vibrational excitation in electron scattering by polar molecules [409, 410].
Computed differential cross sections are in quantitative agreement with available
experimental data for e−HF scattering. Because of the computational simplifica-
tions resulting from combining R-matrix theory with the adiabatic approximation,
this methodology was able to obtain results equivalent to converged close-coupling
calculations including both vibrational and rotational degrees of freedom. Specific
treatment of rotational structure is essential for such molecules because of dynam-
ical rotational screening of the long-range dipole potential. These calculations pro-
vide a detailed analysis of the striking threshold peak structures observed for such
dipolar molecules. For e−H2 vibrational excitation, in a direct comparison with
the FONDA (first-order nondegenerate adiabatic) approximation [252, 2, 254], and
with benchmark vibrational close-coupling results [411], the EMAP method was
found to be computationally efficient and reasonably accurate at energies somewhat
above threshold and away from a scattering resonance [234].
The EMAP method has been used to compute elastic scattering and symmetric-
stretch vibrational excitation cross sections for electron scattering by CO2 [235].
This is one of the first ab initio calculations of vibrational excitation for a poly-
atomic molecule. The results are in good agreement with experiment, which shows
unusually large low-energy cross sections. The theory identifies a near-threshold
9.3 Vibronic R-matrix theory 175
1
y(q)γ 2 (q)[res (q) − (q)]−1 γ 2 (q)y † (q).
1 1
R1FN (q; ) = (9.2)
2
Using the basic rationale of EMA theory [267], the parametric function (q) be-
comes op (q) = E−Hn when the kinetic energy of nuclear motion cannot be neg-
lected. However, the operator (res (q) − op (q)) has a well-defined c-number value
in vibrational eigenstates determined by the eigenvalue equation
vibrational excitation is
1 1 1 †
( pLv|R1 | p L v ) = v y pL γ 2 s (E s − E)−1 s γ 2 y p L v . (9.4)
2 s
where "d is a postulated discrete state, defined for q ≤ qd , that interacts with
the background electronic continuum. For this precursor state, the projection
integrals on the electronic wall of the vibronic hypercylinder are ( pLv|s) =
(A# p Y L χv |!s )r =a , and the corresponding projection integral on the dissociation
cap is (d|s) = ("d |!s )q=qd = χs (qd ), a normalized eigenfunction of Eq. (9.3). The
vibronic R-matrix for this precursor resonance state is
1
( pLv|Rres | p L v ) = ( pLv|s)(E s − E)−1 (s| p L v ), (9.5)
2 s
were valid. The NADP formalism postulates this to be true, implying that the R-
matrix connecting wall and cap surfaces is
1 1 †
(d|R1 | pLv) = χs (qd )(E s − E)−1 s|γ 2 y pL |v ,
2 s
9.3 Vibronic R-matrix theory 177
[73] Cottingham, W.N. and Greenwood, D.A. (1998). Introduction to the Standard
Model of Particle Physics (Cambridge University Press, New York).
[93] Edelen, D.G.B. (1969). Nonlocal Variations and Local Invariance of Fields
(American Elsevier, New York).
[121] Feynman, R.P. (1961). Quantum Electrodynamics (Benjamin, New York).
[217] Leader, E. and Predazzi, E. (1996). An Introduction to Gauge Theories and Modern
Particle Physics, Vols. 1 and 2 (Cambridge University Press, New York).
[336] Ramond, P. (1989). Field Theory: A Modern Primer, 2nd edition (Addison-Wesley,
Redwood City, California).
[340] Renton, P. (1990). Electroweak Interactions (Cambridge University Press,
New York).
[352] Ryder, L.H. (1985). Quantum Field Theory (Cambridge University Press,
New York).
[373] Schwinger, J. (1958). Quantum Electrodynamics (Dover, New York).
[391] Swanson, M.S. (1992). Path Integrals and Quantum Processes (Academic Press,
New York).
[421] Weinberg, S. (1995). The Quantum Theory of Fields, Vol. I (Cambridge University
Press, New York).
[422] Weinberg, S. (1996). The Quantum Theory of Fields, Vol. II (Cambridge University
Press, New York).
181
182 10 Relativistic Lagrangian theories
empirical data to which it is applicable. Classical and quantized QED will be used
here to introduce the basic formalism of field theory, including the variational the-
ory of invariance properties. This theory, especially gauge invariance, is central to
recent developments of electroweak theory (EWT) and quantum chromodynamics
(QCD).
Because both the Maxwell and the Dirac fields appear to be truly elementary, not
simply approximate models of inaccessible underlying structures, QED serves as a
model for a more general field theory that might ultimately describe all phenomena.
In recent years, some of the mystery has been removed from nuclear and high-energy
physics by the discovery (or invention) of quarks and gauge boson fields, as analogs
and generalizations of electrons and the electromagnetic field, respectively. Starting
from QED as a generic gauge theory, the formalism of electroweak theory (EWT)
is developed here, following the implications of nonabelian gauge symmetries in
field theory. The same principles apply to the quantum chromodynamics (QCD) of
quarks and gluons, which will not be considered in detail.
Fµν = ∂µ Aν − ∂ν Aµ .
The Dirac and Maxwell fields are coupled through the covariant derivative
Dµ = ∂µ + i(−e/hc)Aµ ,
For zero rest mass, using Eqs. (10.1), the 2 × 2 operator that acts on 2-component
184 10 Relativistic Lagrangian theories
H L = σ · (i hc∇ − eA) − e A0
e
= −cσ · p + A − eφ.
c
This is an exact result of the theory, and must retain its validity for quantized fields. If
there is no external field, any scalar field resulting from coupling u ∗p u q must vanish.
This is required in the vacuum state of the quantized theory, and a valid classical
model should omit such terms. Because the functions u p (x) are elementary Pauli
spinors, the nonscalar field A/(x) has triplet-odd character, describing a vector boson
in quantized theory. Then A/(x) is an odd-parity field, whose matrix elements vanish
except for transitions between negative and positive energies, because massless
fermions have definite helicity. Chirality-breaking virtual transitions are required
in order to produce a nonvanishing mass from the electromagnetic interaction.
In relativistic theory, it is desirable to retain covariant forms in the equations of
motion. Instead of Eq. (10.3), Eq. (10.2) can be rewritten as
defining δ m̂c2 = −e A/. This is an operator or matrix in the Pauli spinor representa-
tion, since its nonvanishing matrix elements describe transitions between helicity
10.1 Classical relativistic electrodynamics 185
states. If there were a bare mass m 0 , a term m 0 c2 ψ would be added to the right-hand
side of this equation, defining m̂ = m 0 + δ m̂. Since the bare mass m 0 is an unde-
termined parameter of the theory, the only nonarbitrary value that can be assigned
to it is zero, unless it can be attributed to a physical mechanism omitted from the
field theory. The dynamical mass δm is in principle determined by the theory, but
is found to be infinite in the absence of a renormalization procedure equivalent
to introducing an arbitrary cutoff for divergent sums or integrals over relativistic
momentum transfer [121].
that vanishes when t1 < t0 . If G 0 (x1 , x0 ) is defined in the same way, without δ m̂,
this implies a Lippmann–Schwinger equation
G m (x1 , x0 ) = G 0 (x1 , x0 ) − i d 4 yG 0 (x1 , y)δ m̂(y)G m (y, x0 ).
δφ(x) = 0; δx µ = δa µ .
where
∂L
E= d x 3
φ̇ − L
∂ φ̇
∂L
P = − d 3 x ∇φ.
∂ φ̇
These formulas are valid for the nonrelativistic one-electron Schrödinger equa-
tion. The Lagrangian density is
L = ψ ∗ (i hψ̇ − Hψ).
P = d 3 x ψ ∗ (−i h∇)ψ.
For the Dirac bispinor, the irreducible representation matrix Dab for each helicity
component is a Pauli spin matrix σ multiplied by h/2. Then
1
δ Dab ψb = hδθk σab
k
ψb .
2
Given
∂L/∂(∂µ ψ) = i hcψ † γ 0 γ µ ,
∂µ F µν = (4π/c) j ν .
on the assumption that the current density j µ is not affected by this transformation.
However, the Lagrangian density L A is not invariant, since the coupling term is
modified such that
L A → L A − h j µ ∂µ χ (x).
Thus the equations of motion are gauge invariant, but the action integral is not.
For the Dirac field in an externally determined Maxwell field, the Lagrangian density
including a renormalized mass term is
L D = i hcψ † γ 0 {γ µ Dµ + im}ψ.
Dµ = ∂µ + i(−e/hc)Aµ ,
i hcγ µ Dµ ψ = mc2 ψ.
10.3 Gauge invariance 191
Because
−ieχ (x) ie
e Dµ ψ → ∂µ + ie(∂µ χ ) − (Aµ + hc∂µ χ ) ψ = Dµ ψ,
hc
where τ0 is the unit matrix. This is a product of U (1) factors and 2 × 2 SU(2)
matrices. The fermion field ψ becomes an array consisting of two spinors, each
carrying an additional “weak isospin” index I3 . In EWT this index distinguishes
members of isospin doublets: electrons and neutrinos or down quarks and up quarks.
A local unitary transformation of a fermion field ψ is defined by
ψ → U ψ; ψ † → ψ †U −1 ,
where U commutes with the Dirac matrices γ µ . Interaction with a gauge field E µ
is described by a covariant derivative Dµ in the fermion Lagrangian density
E µ → U E µU −1 − (∂µU )U −1 , (10.6)
For the abelian phase factor of the U (1) group, U = exp{ieχ(x)} and E µ =
−(ie/hc)Aµ , which produces the usual gauge transformation Aµ → Aµ +
hc∂µ χ (x). In electroweak theory, a U (1) gauge field Bµ is defined such that
U = exp{− 12 ig yχ(x)} and E µ = (ig y/2hc)Bµ . The gauge transformation is
Bµ → Bµ + hc∂µ χ (x).
10.3 Gauge invariance 193
in any inertial frame, where Tν0 = Tν0 (ψ) + Tν0 (A) such that
using the antisymmetry of F µλ and the field equations. For a noninteracting field the
residual term vanishes and the resulting symmetrical energy–momentum tensor is
1 µλ 1 σλ
#µν (A) = − F Fνλ + gνµ F Fσ λ .
4π 16π
The electric and magnetic field vectors in vacuo are, respectively,
1 ∂A
E=− − ∇φ, B = ∇ × A.
c ∂t
Using italic indices for 3-vectors, with Aµ = (φ, −A), ∂µ = (∂/c∂t, ∇),
Ei = ∂0 Ai − ∂i A0 = F0i ,
Bi = −i jk (∂ j Ak − ∂k A j ) = −i jk F jk .
In terms of these vector fields, the Maxwell field Lagrangian density is (1/8π )(E 2 −
B2 ), and the field energy and momentum are
E 2 + B2 E ×B
E(A) = 3
d x , P(A) = d3x .
8π 4π c
For interacting fields, the Maxwell field energy is not separately conserved.
A gauge-covariant derivation follows from the inhomogeneous field equations
(Maxwell equations in vacuo),
∂
∇×E + B = 0, ∇ · E = 4πρ,
c∂t
∂ 4π
∇×B − E= j, ∇ · B = 0.
c∂t c
The local conservation law follows immediately [318] from
1 ∂ 4π
B · (∇ × E) − E · (∇ × B) = ∇ · (E × B) = − (E 2 + B2 ) − E · j.
2 c∂t c
Expressed in terms of the field energy and momentum densities defined above, this
is
∂ E 2 + B2 E ×B
+ c2 ∇ · = −E · j.
∂t 8π 4π c
When integrated over a finite volume, the divergence term becomes a surface
integral of the Poynting vector [318]. Expressed in covariant notation, ∂ ν #0ν =
− 1c F 0ν jν , showing the explicit modification of energy–momentum conservation
due to the final dissipative term here.
10.4 Energy and momentum of the coupled fields 199
E 2 + B2 E ×B
(x) = , π (x) = ,
8π 4π c
the relativistic mass density of the field is defined by
without invoking a specific requirement for a symmetric stress tensor. For a plane-
wave radiation field in vacuo,
(E 2 − B2 )2 + 4(E · B)2
µ2 c 4 = = 0,
64π 2
because the field vectors are orthogonal and of equal magnitude (in mixed Gaussian
units).
Using italic indices for spatial 3-vectors, but retaining vector notation for the abstract
3-vector index of the matrices τ and the gauge field W, the fermion “Hamiltonian”
operator is
0 λ g 1
HW = −i h(cγ γ Dλ − ∂t ) = −cσi pi − τ · Wi + gτ · W0 .
2c 2
The SU(2) energy–momentum tensor can be symmetrized and made gauge invariant
by adding an incremental tensor
µ 1 µλ g
Tν = W · ∂λ Wν + Wν × Wλ .
4π hc
Because ∂µ Tνµ does not reduce to terms that vanish even for a noninteracting
field, this construction must be verified. The energy and 3-momentum of the gauge
field derived from the resulting symmetric energy–momentum tensor #µν are
1
E(W ) = − d 3 x(W0i · W0i − W jk · W jk | j<k ),
8π
1
Pi (W ) = d 3 x W0 j · Wi j .
4π c
The local conservation law for the interacting gauge field can be derived from
the covariant field equations, as was done above for the Maxwell field. Using the
SU(2) field equations and expanding (∂ ν W0λ ) · Wνλ as
1 0 νλ g 0 ν
∂ (W · Wνλ ) + [∂ (W × Wλ ) − 2∂ ν (W0 × Wλ )] · Wνλ ,
4 2hc
there are two distinct expansions of the 4-divergence
1 ν 0λ 1 1 0 0j
− ∂ (W · Wνλ ) = − ∂ i (W0 j · Wi j ) − ∂ (W · W0 j )
4π 4π 4π
1 0 νλ 1
=− ∂ (W · Wνλ ) − W0λ · JλW
16π c
g
− [∂ 0 (Wν × Wλ ) − 2∂ ν (W0 × Wλ )] · Wνλ .
8π hc
The final term here can be expressed in the form
1 0λ g ν 1
W · W × Wνλ + = W0λ · JλW − jλW + .
c 4π h c
10.5 The Standard Model 201
The effect of this term is to remove the self-interaction current density JλW − jλW
from the dissipative term in the local conservation law. The residual invariant
can be shown to vanish identically. Using the identity
1 g
(∂ ν Wλ ) × Wνλ = (∂ ν Wλ − ∂ λ Wν ) × Wνλ ≡ (Wν × Wλ ) × Wνλ ,
2 2hc
the residual invariant is
g
= [2∂ ν (W0 × Wλ ) − ∂ 0 (Wν × Wλ ) − 2W0λ × Wν ] · Wνλ
8π hc
g2
= W0 · [2Wλ × (Wν × Wνλ ) + (Wν × Wλ ) × Wνλ ]
8π h2 c2
g2
= W0 · [Wν (Wλ · Wνλ ) + Wλ (Wν · Wνλ )] ≡ 0.
8π h2 c2
The implied local conservation law for the interacting SU(2) field is
1
∂ ν #0ν = − W0λ · jλW .
c
Further analysis is required to conclude that the local mass density must necessarily
vanish.
Dµ = ∂µ + i(e0 /hc)Uµ .
Here e0 is the magnitude of the electronic charge. The equation of motion for the
massless fermion field is
i hcγ µ Dµ ψ = 0.
10.5 The Standard Model 203
Defining Wµ± = √1 (W 1
2 µ ± Wµ2 ) and τ± = 12 (τ1 ± iτ2 ),
1 1
Uµ = (y + τ3 )Aµ + (−y tan θW + τ3 cot θW )Z µ0
2 2
1
+√ (τ+ Wµ− + τ− Wµ+ ).
2 sin θW
The parameter y is −1 for leptons and 1/3 for quarks, while τ3 = ∓1.
[1] Aashamar, K., Luke, T.M. and Talman, J.D. (1978). Optimized central potentials
for atomic ground-state wavefunctions, At. Data Nucl. Data Tables 22,
443–472.
[2] Abdolsalami, M. and Morrison, M.A. (1987). Calculating vibrational-excitation
cross sections off the energy shell: A first-order adiabatic theory, Phys. Rev. A 36,
5474–5477.
[3] Adhikari, S.K. (1998). Variational Principles and the Numerical Solution of
Scattering Problems (Wiley, New York).
[4] Aissing, G. and Monkhorst, H.J. (1993). On the removal of the exchange singularity
in extended systems, Int. J. Quantum Chem. S27, 81–89.
[5] Akhiezer, N.I. (1962). The Calculus of Variations (Blaisdell, New York).
[6] Amaldi, U. et al. (1987). Comprehensive analysis of data pertaining to the weak
neutral current and the intermediate vector meson masses, Phys. Rev. D 36,
1385–1407.
[7] Andersen, O.K. (1971). Comments on the KKR wavefunctions; extension of the
spherical wave expansion beyond the muffin tins. In Computational Methods in
Band Theory, eds. P.M. Marcus, J.F. Janak, A.R. Williams (Plenum, New York),
178–182.
[8] Andersen, O.K. (1973). Simple approach to the band-structure problem, Solid State
Commun. 13, 133–136.
[9] Andersen, O.K. (1975). Linear methods in band theory, Phys. Rev. B 12, 3060–3083.
[10] Andersen, O.K. and Jepsen, O. (1977). Advances in the theory of one-electron
states, Physica B 91, 317–328.
[11] Andersen, O.K., Jepsen, O. and Glötzel, D. (1985). Canonical description of the
band structures of metals. In Highlights of Condensed-Matter Theory, Soc. Ital. Fis.
Corso 89, 59–176.
[12] Andersen, O.K. and Kasowski, R.V. (1971). Electronic states as linear combinations
of muffin-tin orbitals, Phys. Rev. B 4, 1064–1069.
[13] Andersen, O.K. and Woolley, R.G. (1973). Muffin-tin orbitals and molecular
calculations: General formalism, Mol. Phys. 26, 905–927.
[14] Baluja, K., Burke, P.G. and Morgan, L.A. (1982). R-matrix propagation program
for solving coupled second-order differential equations, Comput. Phys. Commun.
27, 299–307.
[15] Bardsley, J.N. (1968). The theory of dissociative recombination, J. Phys. B 1,
365–380.
205
206 References and bibliography
[16] Bardsley, J.N., Herzenberg, A. and Mandl, F. (1966). Electron resonances of the
H−2 ion, Proc. Phys. Soc. (London) 89, 305–320.
[17] Bardsley, J.N., Herzenberg, A. and Mandl, F. (1966). Vibrational excitation and
dissociative attachment in the scattering of electrons by hydrogen molecules, Proc.
Phys. Soc. (London) 89, 321–340.
[18] Bardsley, J.N. and Mandl, F. (1968). Resonant scattering of electrons by molecules,
Rep. Prog. Phys. 31, 471–531.
[19] Bardsley, J.N. and Wadehra, J.M. (1979). Dissociative attachment and vibrational
excitation in low-energy collisions of electrons with H2 and D2 , Phys. Rev. A 20,
1398–1405.
[20] Berman, M. and Domcke, W. (1984). Projection-operator calculations for shape
resonances: A new method based on the many-body optical approach, Phys. Rev. A
29, 2485–2496.
[21] Berman, M., Estrada, H., Cederbaum, L.S. and Domcke, W. (1983). Nuclear
dynamics in resonant electron–molecule scattering beyond the local approximation:
The 2.3-eV shape resonance in N2 , Phys. Rev. A 28, 1363–1381.
[22] Berman, M., Mündel, C. and Domcke, W. (1985). Projection-operator calculations
for molecular shape resonances: The 2 'u+ resonance in electron-hydrogen
scattering, Phys. Rev. A 31, 641–651.
[23] Bingel, W.A. (1967). A physical interpretation of the cusp conditions for molecular
wave functions, Theor. Chim. Acta 8, 54–61.
[24] Birtwistle, D.T. and Herzenberg, A. (1971). Vibrational excitation of N2 by
resonance scattering of electrons, J. Phys. B 4, 53–70.
[25] Bjorken, J.D. and Drell, S.C. (1965). Relativistic Quantum Fields (McGraw-Hill,
New York).
[26] Blanchard, P. and Brüning, E. (1992). Variational Methods in Mathematical
Physics: A Unified Approach (Springer-Verlag, Berlin).
[27] Blatt, J.M. and Biedenharn, L.C. (1952). The angular distribution of scattering and
reaction cross sections, Rev. Mod. Phys. 24, 258–272.
[28] Blatt, J.M. and Weisskopf, V.F. (1979). Theoretical Nuclear Physics
(Springer-Verlag, New York).
[29] Bleecker, D. (1985). Gauge Theory and Variational Principles (Benjamin,
New York).
[30] Bloch, C. (1957). A unified formulation of the theory of nuclear reactions, Nucl.
Phys. 4, 503–528.
[31] Born, M. and Oppenheimer, J.R. (1927). Ann. Phys. 84, 457.
[32] Brandsden, B.H. (1970). Atomic Collision Theory (Plenum, New York).
[33] Breit, G. (1959). Handbuch der Physik 41, No.1 (Springer, Berlin).
[34] Breit, G. and Wigner, E.P. (1936). Capture of slow neutrons, Phys. Rev. 49,
519–531.
[35] Brenig, W. (1957). Zweiteilchennäherungen des Mehrkörperproblems I, Nucl. Phys.
4, 363–374.
[36] Brescansin, L.M., Leite, J.R. and Ferreira, L.G. (1979). A study of the ground states
and ionization energies of H2 , C2 , N2 , F2 , and CO molecules by the variational
cellular method, J. Chem. Phys. 71, 4923–4930.
[37] Brillouin, L. (1934). Actual. Sci. Ind. 159, 28.
[38] Brown, L.M. (1958). Two-component fermion theory, Phys. Rev. 111, 957–964.
[39] Brown, R.G. and Ciftan, M. (1983). Generalized non-muffin-tin band theory, Phys.
Rev. B 27, 4564–4579.
[40] Brown, R.G. and Ciftan, M. (1983). Phys. Rev. B 28, 5992–6007.
References and bibliography 207
[41] Broyden, C.G. (1967). Quasi-Newton methods and their application to function
minimisation, Math. Comput. 21, 368–381.
[42] Broyden, C.G. (1970). The convergence of single-rank quasi-Newton methods,
Math. Comput. 24, 365–382.
[43] Brueckner, K.A. and Wada, W. (1956). Nuclear saturation and two-body forces:
Self-consistent solutions and the effects of the exclusion principle, Phys. Rev. 103,
1008–1016.
[44] Burke, P.G. (1973). The R-matrix method in atomic physics, Comput. Phys.
Commun. 6, 288–302.
[45] Burke, P.G. (1977). Potential Scattering in Atomic Physics (Plenum, New York).
[46] Burke, P.G. and Robb, W.D. (1972). Elastic scattering of electrons by hydrogen and
helium atoms, J. Phys. B 5, 44–54.
[47] Burke, P.G. and Robb, W.D. (1975). The R-matrix theory of atomic processes, Adv.
At. Mol. Phys. 11, 143–214.
[48] Burke, P.G. and Seaton, M.J. (1971). Numerical solutions of the integro-differential
equations of electron–atom collision theory, Methods Comput. Phys. 10, 1–80.
[49] Burke, P.G. and West, J.B., eds. (1987). Electron–Molecule Scattering and
Photoionization (Plenum, New York).
[50] Burke, P.G, Hibbert, A. and Robb, W.D. (1971). Electron scattering by complex
atoms, J. Phys. B 4, 153–161.
[51] Butler, W.H. and Nesbet, R.K. (1990). Validity, accuracy, and efficiency of
multiple-scattering theory for space-filling scatterers, Phys. Rev. B 42,
1518–1525.
[52] Butler, W.H. and Zhang, X.-G. (1991). Accuracy and convergence properties of
multiple-scattering theory in three dimensions, Phys. Rev. B 44, 969–983.
[53] Butler, W.H., Zhang, X.-G. and Gonis, A. (1992). The Green function cellular
method and its relation to multiple scattering theory, Mat. Res. Symp. Proc. 253,
205–210.
[54] Buttle, P.J.A. (1967). Solution of coupled equations by R-matrix techniques,
Phys. Rev. 160, 719–729.
[55] Cairo, L. (1965). Variational Techniques in Electromagnetism, tr. G.D. Sims
(Gordon and Breach, New York).
[56] Callaway, J. and March, N.H. (1985). Density functional methods: theory and
applications, Solid State Phys. 38, 135–221.
[57] Callen, H.B. (1963). Thermodynamics (Wiley, New York).
[58] Carathéodory, C. (1967). Calculus of Variations and Partial Differential Equations
of the First Order, Part II, tr. R.B. Dean and J.J. Brandstatter (Holden-Day,
San Francisco).
[59] Cederbaum, L.S. and Domcke, W. (1981). Local against nonlocal complex potential
in resonant electron–molecule scattering, J. Phys. B 14, 4665–4689.
[60] Ceperley, D.M. and Alder, B.J. (1980). Ground state of the electron gas by a
stochastic method, Phys. Rev. Lett. 45, 566–569.
[61] Chandra, N. and Temkin, A. (1976). Hybrid theory and calculation of e – N2
scattering, Phys. Rev. A 13, 188–215.
[62] Chang, E.S. and Fano, U. (1972). Theory of electron–molecule collisions by frame
transformation, Phys. Rev. A 6, 173–185.
[63] Chang, E.S. and Temkin, A. (1969). Rotational excitation of diatomic molecules by
electron impact, Phys. Rev. Lett. 23, 399–403.
[64] Chang, E.S. and Temkin, A. (1970). Rotational excitation of diatomic molecular
systems. II. H+2 , J. Phys. Soc. Japan 29, 172–179.
208 References and bibliography
[65] Chase, D.M. (1956). Adiabatic approximation for scattering processes, Phys. Rev.
104, 838–842.
[66] Chatwin, R.A. and Purcell, J.E. (1971). Approximate solution of a Sturm–Liouville
system using nonorthogonal expansions: Application to α–α nuclear scattering,
J. Math. Phys. 12, 2024–2030.
[67] Chen, J.C.Y. (1966). Dissociative attachment in rearrangement electron collision
with molecules, Phys. Rev. 148, 66–73.
[68] Clementi, E. and Roetti, C. (1974). Roothaan–Hartree–Fock atomic wavefunctions,
At. Data Nucl. Data Tables 14, 177–478.
[69] Colle, R. and Nesbet, R.K. (2001). Optimized effective potential in finite-basis-set
treatment, J. Phys. B 34, 2475–2480.
[70] Colle, R. and Salvetti, O. (1975). Approximate calculation of the correlation energy
for the closed shells, Theoret. Chim. Acta 37, 329–334.
[71] Collins, J.C. (1984). Renormalization (Cambridge University Press, Cambridge).
[72] Condon, E.U. and Shortley, G.H. (1935). The Theory of Atomic Spectra (Cambridge
University Press, Cambridge).
[73] Cottingham, W.N. and Greenwood, D.A. (1998). An Introduction to the Standard
Model of Particle Physics (Cambridge University Press, New York).
[74] Courant, R. and Hilbert, D. (1953). Methods of Mathematical Physics (Interscience,
New York).
[75] Csaszar, P. and Pulay, P. (1984). Geometry optimization by direct inversion in the
iterative subspace, J. Mol. Struct. 114, 31–34.
[76] Damburg, R.J. and Karule, E. (1967). A modification of the close-coupling
approximation for e – H scattering allowing for the long-range interaction, Proc.
Phys. Soc. (London) 90, 637–640.
[77] Dennis, J.E. and Schnabel, R.B. (1983). Numerical Methods for Unconstrained
Optimization and Nonlinear Equations (Prentice-Hall, Englewood Cliffs,
New Jersey).
[78] Dieudonné, J. (1981). History of Functional Analysis (North-Holland, Amsterdam).
[79] Dirac, P.A.M. (1930). Note on exchange phenomena in the Thomas atom, Proc.
Camb. Phil. Soc. 26, 376–385.
[80] Dirac, P.A.M. (1933). Physikalische Z. der Sowjetunion 3, 64.
[81] Dirac, P.A.M. (1958). Quantum Mechanics, 4th edition (Oxford University Press,
London).
[82] Domcke, W. (1991). Theory of resonance and threshold effects in electron–molecule
collisions: the projection-operator approach, Phys. Rep. 208, 97–188.
[83] Domcke, W. and Cederbaum, L.S. (1977). Theory of the vibrational structure of
resonances in electron–atom scattering, Phys. Rev. A 16, 1465–1482.
[84] Domcke, W. and Cederbaum, L.S. (1981). On the interpretation of low-energy
electron-HCl scattering phenomena, J. Phys. B 14, 149–173.
[85] Domcke, W., Cederbaum, L.S. and Kaspar, F. (1979). Threshold phenomena in
electron–molecule scattering: a nonadiabatic theory, J. Phys. B 12, L359–L364.
[86] Domcke, W. and Mündel, C. (1985). Calculation of cross sections for vibrational
excitation and dissociative attachment in HCl and DCl beyond the local complex
potential approximation, J. Phys. B 18, 4491–4509.
[87] Domps, A., Reinhard, P.-G. and Suraud, E. (1998). Time-dependent Thomas–
Fermi approach to electron dynamics in metal clusters, Phys. Rev. Lett. 80,
5520–5523.
[88] Dovesi, R., Saunders, V.R., Roetti, C., Causà, M., Harrison, N.M., Orlando, R. and
Aprà, E. (1996). CRYSTAL95 User’s Manual (University of Torino, Torino, Italy).
References and bibliography 209
[112] Ferreira, L.G. and Leite, J.R. (1978). Variational cellular model of the molecular
and crystal electronic structure, Phys. Rev. Lett. 40, 49–52.
[113] Ferreira, L.G. and Leite, J.R. (1978). General formulation of the variational cellular
method for molecules and crystals, Phys. Rev. A 18, 335–343.
[114] Ferreira, L.G. and Leite, J.R. (1979). Self-consistent calculation of the electronic
structure of N2 and CO by the variational cellular method, Phys. Rev. A 20,
689–699.
[115] Feshbach, H. (1958). Unified theory of nuclear reactions, Ann. Phys. (New York) 5,
357–390.
[116] Feshbach, H. (1962). A unified theory of nuclear reactions, Ann. Phys. (New York)
19, 287–313.
[117] Feynman, R.P. (1939). Forces in molecules, Phys. Rev. 56, 340–343.
[118] Feynman, R.P. (1948). A relativistic cut-off for classical electrodynamics, Phys.
Rev. 74, 939–946.
[119] Feynman, R.P. (1948). Space-time approach to nonrelativistic quantum mechanics,
Rev. Mod. Phys. 20, 367–387.
[120] Feynman, R.P. (1949). Space-time approach to quantum electrodynamics, Phys.
Rev. 76, 769–789.
[121] Feynman, R.P. (1961). Quantum Electrodynamics (Benjamin, New York).
[122] Feynman, R.P. and Gell-Mann, M. (1958). Theory of the Fermi interaction, Phys.
Rev. 109, 193–198.
[123] Feynman, R.P. and A.R. Hibbs (1965). Quantum Mechanics and Path Integrals
(McGraw-Hill, New York).
[124] Fletcher, R. (1970). A new approach to variable metric algorithms, Comput. J. 13,
317–322.
[125] Fletcher, R. (1987). Practical Methods of Optimization, 2nd edition (Wiley,
New York).
[126] Fliflet, A.W. and McKoy, V. (1978). Discrete basis set method for
electron–molecule continuum wave functions, Phys. Rev. A 18, 2107–2114.
[127] Fock, V. (1930). Z. Phys. 61, 126.
[128] Fock, V. (1957). The Theory of Space, Time and Gravitation, tr. N. Kemmer
(Pergamon Press, London).
[129] Fried, H. (1972). Functional Methods and Models in Quantum Field Theory (MIT
Press, Cambridge, Massachusetts).
[130] Froese Fischer, C. (1977). The Hartree–Fock Method for Atoms (Wiley, New York).
[131] Garrod, G. and Percus, J.K. (1964). Reduction of the N -particle variational
problem, J. Math. Phys. 5, 1756–1776.
[132] Gerjuoy, E., Rau, A.R.P. and Spruch, L. (1983). A unified formulation of the
construction of variational principles, Rev. Mod. Phys. 55, 725–774.
[133] Gianturco, F.A. and Stoecklin, T. (1996). The elastic scattering of electrons from
CO2 molecules: I. Close coupling calculations of integral and differential cross
sections, J. Phys. B 29, 3933–3954.
[134] Gianturco, F.A., Jain, A. and Pantano, L. (1987). Electron-methane scattering via a
parameter-free model interaction, J. Phys. B 20, 571–586.
[135] Gianturco, F.A., Thompson, D.G. and Jain, A. (1995). Electron scattering from
polyatomic molecules using a single-center expansion formulation, in
Computational Methods for Electron–Molecule Collisions, eds. W.M. Huo and
F.A. Gianturco (Plenum, New York), pp. 75–118.
[136] Gilbert, T.L. (1975). Hohenberg–Kohn theorem for nonlocal external potentials,
Phys. Rev. B 12, 2111–2120.
References and bibliography 211
[137] Gill, P.E., Murray, W. and Wright, M.H. (1981). Practical Optimization (Academic
Press, New York).
[138] Gillan, C.J., Nagy, O., Burke, P.G., Morgan, L.A. and Noble, C.J. (1987). Electron
scattering by nitrogen molecules, J. Phys. B 20, 4585–4603.
[139] Gillan, C.J., Tennyson, J. and Burke, P.G. (1995). The UK molecular R-matrix
scattering package: a computational perspective, in Computational Methods for
Electron–Molecule Collisions, eds. W.M. Huo and F.A. Gianturco (Plenum,
New York), pp. 239–254.
[140] Glashaw, S.L. (1961). Partial symmetries of weak interactions, Nucl. Phys. 22,
579–588.
[141] Glimm, A. and Jaffe, A. (1981). Quantum Mechanics - A Functional Integral Point
of View (Springer-Verlag, Berlin).
[142] Glowinski, R. (1983). Numerical Methods for Nonlinear Variational Problems
(Springer-Verlag, Berlin).
[143] Golden, D.E. (1966). Low-energy resonances in e− −N2 total scattering cross
sections: the temporary formation of N− 2 , Phys. Rev. Lett. 17, 847–851.
[144] Golden, D.E., Lane, N.F., Temkin, A. and Gerjuoy, E. (1971). Low energy
electron–molecule scattering experiments and the theory of rotational excitation,
Rev. Mod. Phys. 43, 642–678.
[145] Goldfarb, D. (1970). A family of variable-metric methods derived by variational
means, Math. Comput. 24, 23–26.
[146] Goldstein, H. (1983). Classical Mechanics, 2nd edition (Wiley, New York).
[147] Goldstine, H.H. (1980). A History of the Calculus of Variations from the 17th
Through the 19th Century (Springer-Verlag, Berlin).
[148] Gonis, A. and Butler, W.H. (2000). Multiple Scattering in Solids (Springer,
New York).
[149] Gonze, X. and Scheffler, M. (1999). Exchange and correlation kernels at the
resonance frequency: implications for excitation energies in density-functional
theory, Phys. Rev. Lett. 82, 4416–4419.
[150] Gordon, R.G. and Kim, Y.S. (1972). Theory for the forces between closed-shell
atoms and molecules, J. Chem. Phys. 56, 3122–3133.
[151] Görling, A. (1999). New KS method for molecules based on an exchange charge
density generating the exact KS exchange potential, Phys. Rev. Lett. 83,
5459–5462.
[152] Görling, A. and Ernzerhof, M. (1995). Energy differences between Kohn–Sham and
Hartree–Fock wave functions yielding the same electron density, Phys. Rev. A 51,
4501–4513.
[153] Grimm-Bosbach, T., Thümmel, H.T., Nesbet, R.K. and Peyerimhoff, S.D. (1996).
Calculation of cross sections for rovibrational excitation of N2 by electron impact,
J. Phys. B 29, L105–L112.
[154] Gross, E.K.U. and Kohn, W. (1985). Local density-functional theory of
frequency-dependent linear response, Phys. Rev. Lett. 55, 2850–2852.
[155] Gross, E.K.U., Dobson, J.F. and Petersilka, M. (1996). Density functional theory of
time-dependent phenomena, in Density Functional Theory, ed. R.F. Nalewajski,
Series “Topics in Current Chemistry” (Springer, Berlin).
[156] Guralnik, G.S., Hagen, C.R. and Kibble, T.W.B. (1964). Global conservation laws
and massless particles, Phys. Rev. Lett. 13, 585–587.
[157] Gyorffy, B.L. (1972). Coherent-potential approximation for a nonoverlapping
muffin-tin potential model of random substitutional alloys, Phys. Rev. B 5,
2382–2384.
212 References and bibliography
[158] Ham, F. and Segal, B. (1961). Energy bands in periodic lattices – Green’s function
method, Phys. Rev. 124, 1786–1796.
[159] Harris, J. (1979). The role of occupation numbers in HKS theory, Int. J. Quantum
Chem. S13, 189–192.
[160] Harris, J. (1984). Adiabatic-connection approach to Kohn–Sham theory, Phys. Rev.
A 29, 1648–1659.
[161] Harris, J. and Jones, R.O. (1974). The surface energy of a bounded electron gas,
J. Phys. F 4, 1170–1186.
[162] Hartree, D.H. (1928). The wave mechanics of an atom with a non-Coulomb central
field. Part II. Some results and discussion, Proc. Camb. Phil. Soc. 24, 111–132.
[163] Hartree, D.H. (1957). The Calculation of Atomic Structures (Wiley, New York).
[164] Hazi, A.U., Orel, A.E. and Rescigno, T.N. (1981). Ab initio study of dissociative
attachment of low-energy electrons to F2 , Phys. Rev. Lett. 46, 918–922.
[165] Head, J.D., Weiner, B. and Zerner, M.C. (1988). A survey of optimization
procedures for stable structures and transition states, Int. J. Quantum Chem. 33,
177–186.
[166] Head, J.D. and Zerner, M.C. (1988). Newton-based optimization methods for
obtaining molecular conformation, Adv. Quantum Chem. 20, 241–290.
[167] Herzenberg, A. (1968). Oscillatory energy dependence of resonant
electron–molecule scattering, J. Phys. B 1, 548–558.
[168] Hestenes, M.R. (1966). Calculus of Variations and Optimal Control Theory (Wiley,
New York).
[169] Higgs, P.W. (1964). Broken symmetries and the masses of gauge bosons, Phys. Rev.
Lett. 13, 508–509.
[170] Higgs, P.W. (1966). Spontaneous symmetry breakdown without massless bosons,
Phys. Rev. 145, 1156–1163.
[171] Hill, E.L. (1951). Hamilton’s principle and the conservation theorems of
mathematical physics, Rev. Mod. Phys. 23, 253–260.
[172] Höck, A. and Engel, E. (1998). Pseudopotentials from orbital-dependent
exchange-correlation functionals, Phys. Rev. A 58, 3578–3581.
[173] Hohenberg, P. and Kohn, W. (1964). Inhomogeneous electron gas, Phys. Rev. 136,
B864–B871.
[174] Hoyle, F.M. and Narlikar, J.V. (1974). Action at a Distance in Physics and
Cosmology (Freeman, San Francisco).
[175] Hulthén, L. (1944). Kgl. Fysiogr. Sällsk. Lund Förh. 14, No.21.
[176] Hulthén, L. (1948). Arkiv. Mat. Astron. Fysik 35A, No.25.
[177] Huo, W.M. (1995). The Schwinger variational method, in Computational Methods
for Electron–Molecule Collisions, eds. W.M. Huo and F.A. Gianturco (Plenum,
New York), pp. 326–355.
[178] Huo, W.M. and Gianturco, F.A. (1995). Computational Methods for Electron–
Molecule Collisions (Plenum, New York).
[179] Hurley, A.C. (1954). The electrostatic calculation of molecular energies II.
Approximate wave functions and the electrostatic method, Proc. Roy. Soc. A 226,
179–192.
[180] Hurley, A.C. (1976). Electron Correlation in Small Molecules (Academic Press,
New York).
[181] Ivanov, S., Hirata, S. and Bartlett, R.J. (1999). Exact exchange treatment for
molecules in finite-basis-set Kohn–Sham theory, Phys. Rev. Lett. 83,
5455–5458.
[182] Jackson, J.D. (1975). Classical Electrodynamics, 2nd edition (Wiley, New York).
References and bibliography 213
[183] Jackson, J.L. (1951). A variational appeoach to nuclear reactions, Phys. Rev. 83,
301–304.
[184] Jacob, M. and Wick, G.C. (1959). On the general theory of collisions for particles
with spin, Ann. Phys. (New York) 7, 404–428.
[185] Janak, J.F. (1978). Proof that ∂ E/∂n i = i in density-functional theory, Phys. Rev.
B 18, 7165–7168.
[186] Jauch, J.M. and Rohrlich, F. (1985). The Theory of Photons and Electrons
(Springer-Verlag, Berlin).
[187] Jeffreys, H. and Jeffreys, B.S. (1956). Methods of Mathematical Physics, 3rd
edition (Cambridge University Press, New York).
[188] Joachain, C.J. (1975). Quantum Collision Theory, Vol.1 (North-Holland,
Amsterdam).
[189] Jones, R.O. and Gunnarsson, O. (1989). The density functional formalism, its
applications and prospects, Rev. Mod. Phys. 61, 689–746.
[190] Itzykson, C. and Zuber, J.-B. (1980). Quantum Field Theory (McGraw-Hill,
New York).
[191] Kamimura, M. (1977). A coupled channel variational method for microscopic study
of reactions between complex nuclei, Prog. Theor. Phys. Suppl. 62, 236–294.
[192] Kapur, P.L. and Peierls, R. (1938). The dispersion formula for nuclear reactions,
Proc. Roy. Soc. (London) A166, 277–295.
[193] Kasowski, R.V. and Andersen, O.K. (1972). Muffin tin orbitals in open structures,
Solid State Commun. 11, 799–802.
[194] Kato, T. (1950). Variational methods in collision problems, Phys. Rev. 80, 475.
[195] Kato, T. (1951). Upper and lower bounds of scattering phases, Prog. Theor. Phys. 6,
394–407.
[196] Kato, T. (1957). On the eigenfunctions of many-particle systems in quantum
mechanics, Commun. Pure Appl. Math. 10, 151–177.
[197] Kaufmann, K., Baumeister, W. and Jungen, M. (1989). Universal Gaussian basis
sets for an optimum representation of Rydberg and continuum wavefunctions,
J. Phys. B 22, 2223–2240.
[198] Kibble, T.W.B. (1967). Symmetry breaking in non-abelian gauge theories, Phys.
Rev. 155, 1554–1561.
[199] Klein, F. (1918). Nachr. Ges. Wiss. Göttingen, 171.
[200] Kline, M. (1972). Mathematical Thought from Ancient to Modern Times, 3 vols.
(Oxford University Press, New York).
[201] Kohl, H. and Dreizler, R.M. (1986). Time-dependent density-functional theory:
conceptual and practical aspects, Phys. Rev. Lett. 56, 1993–1995.
[202] Kohn, W. (1948). Variational methods in nuclear collision problems, Phys. Rev. 74,
1763–1772.
[203] Kohn, W. (1952). Variational methods for periodic lattices, Phys. Rev. 87, 472–481.
[204] Kohn, W. and Rostoker, N. (1954). Solution of the Schrödinger equation in periodic
lattices with an application to metallic lithium, Phys. Rev. 94, 1111–1120.
[205] Kohn, W. and Sham, L.J. (1965). Self consistent equations including exchange and
correlation effects, Phys. Rev. 140, A1133–A1138.
[206] Koopmans, T. (1933). Über die Zuordnung von Wellenfunktionen und Eigenwerten
zu den einzelnen Elektronen eines Atoms, Physica 1,104–113.
[207] Korringa, J. (1947). On the calculation of the energy of a Bloch wave in a metal,
Physica 13, 392–400.
[208] Kuperschmidt, B.A. (1990). The Variational Principles of Dynamics (World
Scientific, New York).
214 References and bibliography
[209] Lai, C.H. (1981). Gauge Theory of Electromagnetic and Weak Interactions, ed.
C.H. Lai (World Scientific, Singapore).
[210] Lanczos, C. (1966). Variational Principles of Mechanics (University of Toronto
Press, Toronto).
[211] Landau, L.D. (1956). The theory of a Fermi liquid, Zh. Eksp. Teor. Fiz. 30,
1058–1064. [Sov. Phys. JETP 3, 920–925 (1956)]
[212] Landau, L.D. (1957). Oscillations in a Fermi liquid, Zh. Eksp. Teor. Fiz. 32, 59–66.
[Sov. Phys. JETP 5, 101–108 (1957)]
[213] Landau, L.D. and Lifshitz, E.M. (1958). Quantum Mechanics, Non-Relativistic
Theory, translated by J.B. Sykes and J.S. Bell (Pergamon, London).
[214] Lane, A.M. and Thomas, R.G. (1958). R-matrix theory of nuclear reactions, Rev.
Mod. Phys. 30, 257–353.
[215] Lane, N.F. (1980). The theory of electron–molecule collisions, Rev. Mod. Phys. 52,
29–119.
[216] Langreth, D.C. and Mehl, M.J. (1981). Easily implementable nonlocal
exchange-correlation energy functional, Phys. Rev. Lett. 47, 446–450.
[217] Leader, E. and Predazzi, E. (1996). An Introduction to Gauge Theories and Modern
Particle Physics, Vols. 1 and 2 (Cambridge University Press, New York).
[218] Lee, B.W. and Zinn-Justin, J. (1972). Spontaneously broken gauge symmetries.
I. Preliminaries, Phys. Rev. D 5, 3121–3137.
[219] Lee, C., Yang, W. and Parr, R.G. (1988). Development of the Colle-Salvetti
correlation energy formula into a functional of the electron density, Phys. Rev. B 37,
785–789.
[220] Lee, C.M. (1974). Spectroscopy and collision theory. III. Atomic eigenchannel
calculation by a Hartree–Fock-Roothaan method, Phys. Rev. A 10, 584–600.
[221] Leibbrandt, G. (1987). Introduction to noncovariant gauges, Rev. Mod. Phys. 59,
1067–1119.
[222] Levy, M. (1979). Universal variational functionals of electron densities, first-order
density matrices, and natural spin-orbitals and solution of the v-representability
problem, Proc. Natl. Acad. Sci. 76, 6062–6065.
[223] Lieb, E.H. (1983). Density functionals for Coulomb systems, Int. J. Quantum
Chem. 24, 243–277.
[224] Lieb, E.H. and Simon, B. (1977). The Thomas–Fermi theory of atoms, molecules
and solids, Adv. Math. 23, 22–216.
[225] Light, J.C. and Walker, R.B. (1976). An R-matrix approach to the solution of
coupled equations for atom–molecule scattering, J. Chem. Phys. 65, 4272–4282.
[226] Lima, M.A.P., Leite, J.R. and Fazio, A. (1971). Theoretical study of the F2 molecule
using the variational cellular method, J. Phys. B 14, L533–L535.
[227] Lino, A.T., Leite, J.R., Ferraz, A.C. and Takahashi, E.K. (1987). Self-consistent
formulation of the variational cellular method applied to periodic structures –
results for sodium and silicon, J. Phys. Chem. Solids 48, 911–919.
[228] Lippmann, B.A. and Schwinger, J. (1950). Variational principles for scattering
processes. I, Phys. Rev. 79, 469–480.
[229] Loucks, T.L. (1967). Augmented Plane Wave Method (Benjamin, New York),
pp. 98–103.
[230] Lucchese, R.R., Takatsuka, K. and McKoy, V. (1986). Applications of the
Schwinger variational principle to electron–molecule collisions and molecular
photoionization, Phys. Rep. 131, 147–221.
[231] March, N.H. (1957). The Thomas–Fermi approximation in quantum mechanics,
Adv. Phys. 6, 1–101.
References and bibliography 215
[232] March, N.H. (1975). Self-consistent Fields in Atoms (Pergamon Press, Oxford).
[233] Mazevet, S., Morrison, M.A., Boydstun, O. and Nesbet, R.K. (1999). Inclusion of
nonadiabatic effects in calculations on vibrational excitation of molecular hydrogen
by low-energy electron impact, Phys. Rev. A 59, 477–489.
[234] Mazevet, S., Morrison, M.A., Boydstun, O. and Nesbet, R.K. (1999). Adiabatic
treatments of vibrational dynamics in low-energy electron–molecule scattering,
J. Phys. B 32, 1269–1294.
[235] Mazevet, S., Morrison, M.A., Morgan, L.A. and Nesbet, R.K. (2001). Virtual-state
effects on elastic scattering and vibrational excitation of CO2 by electron impact.
Phys. Rev. A 64, 040701.
[236] Mazevet, S., Morrison, M.A. and Nesbet, R.K. (1998). Application of the
nonadiabatic phase matrix method to vibrational excitation near a short-lived
resonance: the case of e-H2 scattering, J. Phys. B 31, 4437–4448.
[237] McCurdy, C.W., Rescigno, T.N. and Schneider, B.I. (1987). Interrelation between
variational principles for scattering amplitudes and generalized R-matrix theory,
Phys. Rev. A 36, 2061–2066.
[238] McLachlan, A.D. and Ball, M.A. (1964). Time-dependent Hartree–Fock theory for
molecules, Rev. Mod. Phys. 36, 844–855.
[239] McWeeny, R. (1989). Methods of Molecular Quantum Mechanics (Academic Press,
New York).
[240] Mercier, A. (1959). Analytical and Canonical Formalism in Physics (Interscience,
New York).
[241] Mermin, N.D. (1965). Thermal properties of the inhomogeneous electron gas, Phys.
Rev. 137, A1441–A1443.
[242] Merzbacher, E. (1961). Quantum Mechanics (Wiley, New York).
[243] Messiah, A. (1961). Quantum Mechanics, Vol. I (North Holland, Amsterdam).
[244] Miller, W.H. and Jansen op de Haar, B.M.D.D. (1987). A new basis set for quantum
scattering calculations, J. Chem. Phys. 86, 6213–6220.
[245] Mito, Y. and Kamimura, M. (1976). The generator coordinate method for
composite-particle scattering based on the Kohn-Hulthén variational principle,
Prog. Theor. Phys. 56, 583–598.
[246] Moiseiwitsch, B.L. (1966). Variational Principles (Interscience, New York).
[247] Morgan, J. van W. (1977). Integration of Poisson’s equation for a complex system
with arbitrary geometry, J. Phys. C 10, 1181–1202.
[248] Morgan, L.A. (1984). A generalized R-matrix propagation program for solving
coupled second-order differential equations, Comput. Phys. Commun. 31, 419–422.
[249] Morgan, L.A. (1998). Virtual states and resonances in electron scattering by CO2 ,
Phys. Rev. Lett. 80, 1873–1875.
[250] Morgan, L.A. and Burke, P.G. (1988). Low-energy electron scattering by HF,
J. Phys. B 21, 2091–2105.
[251] Morgan, L.A., Burke, P.G. and Gillan, C.J. (1990). Low-energy electron scattering
by HCl, J. Phys. B 23, 99–113.
[252] Morrison, M.A. (1986). A first-order nondegenerate adiabatic theory for calculating
near-threshold cross sections for rovibrational excitation of molecules by electron
impact, J. Phys. B 19, L707–L715.
[253] Morrison, M.A. (1988). Near-threshold electron–molecule scattering, Adv. At. Mol.
Phys. 24, 51–156.
[254] Morrison, M.A., Abdolsalami, M. and Elza, B.K. (1991). Improved accuracy in
adiabatic cross sections for low-energy rotational and vibrational excitation, Phys.
Rev. A 43, 3440–3459.
216 References and bibliography
[255] Morse, P.M. and Feshbach, H. (1953). Methods of Theoretical Physics, Vols. I
and II (McGraw-Hill, New York).
[256] Moruzzi, V.L., Janak, J.F. and Williams, A.R. (1978). Calculated Electronic
Properties of Metals (Pergamon, New York).
[257] Mott, N.F. and Massey, H.S.W. (1965). The Theory of Atomic Collisions (Oxford
University Press, New York).
[258] Mündel, C. and Domcke, W. (1984). Nuclear dynamics in resonant electron–
molecule scattering beyond the local approximation: model calculations on
dissociative attachment and vibrational excitation, J. Phys. B 17, 3593–3616.
[259] Mündel, C., Berman, M. and Domcke, W. (1985). Nuclear dynamics in
resonant electron–molecule scattering beyond the local approximation: Vibrational
excitation and dissociative attachment in H2 and D2 , Phys. Rev. A 32,
181–193.
[260] Nagy, A. (1995). Coordinate scaling and adiabatic connection formula for
ensembles of fractionally occupied excited states, Int. J. Quantum Chem. 56,
225–228.
[261] Nesbet, R.K. (1955). Configuration interaction in orbital theories, Proc. Roy. Soc. A
230, 312–321.
[262] Nesbet, R.K. (1961). Approximate methods in the quantum theory of
many-Fermion systems, Rev. Mod. Phys. 33, 28–36.
[263] Nesbet, R.K. (1961). Construction of symmetry-adapted functions in the
many-body problem, J. Math. Phys. 2, 701–709.
[264] Nesbet, R.K. (1968). Analysis of the Harris variational method in scattering theory,
Phys. Rev. 175, 134–142.
[265] Nesbet, R.K. (1969). Anomaly-free variational method for inelastic scattering,
Phys. Rev. 179, 60–70.
[266] Nesbet, R.K. (1971). Where semiclassical radiation theory fails, Phys. Rev. Lett. 27,
553–556.
[267] Nesbet, R.K. (1979). Energy-modified adiabatic approximation for scattering
theory, Phys. Rev. A 19, 551–556.
[268] Nesbet, R.K. (1979). Accurate e− −He cross sections below 19 eV, J. Phys. B 12,
L243–L248.
[269] Nesbet, R.K. (1979). Variational calculations of accurate e− −He cross sections
below 19 eV, Phys. Rev. A 20, 58–70.
[270] Nesbet, R.K. (1980). Variational Methods in Electron–Atom Scattering Theory
(Plenum, New York).
[271] Nesbet, R.K. (1981). The concept of a local complex potential for nuclear motion in
electron–molecule collisions, Comments At. Mol. Phys. 11. 25–35.
[272] Nesbet, R.K. (1984). R-matrix formalism for local cells of arbitrary geometry.
Phys. Rev. B 30, 4230–4234.
[273] Nesbet, R.K. (1984). Asymptotic distorted-wave approximation for
electron–molecule scattering. J. Phys. B 17, L897–L900.
[274] Nesbet, R.K. (1986). Linearized atomic-cell orbital method for energy-band
calculations, Phys. Rev. B 33, 8027–8034.
[275] Nesbet, R.K. (1986). Nonperturbative theory of exchange and correlation in
one-electron quasiparticle states, Phys. Rev. B 34, 1526–1538.
[276] Nesbet, R.K. (1988). Variational methods for cellular models, Phys. Rev. A 38,
4955–4960.
[277] Nesbet, R.K. (1990). Full-potential multiple scattering theory, Phys. Rev. B 41,
4948–4952.
References and bibliography 217
[278] Nesbet, R.K. (1990). Atomic cell method for total energy calculations, Bull. Am.
Phys. Soc. 35, 418.
[279] Nesbet, R.K. (1992). Variational principles for full-potential multiple scattering
theory, Mat. Res. Symp. Proc. 253, 153–158.
[280] Nesbet, R.K. (1992). Internal sums in full-potential multiple scattering theory, Phys.
Rev. B 45, 11491–11495.
[281] Nesbet, R.K. (1992). Full-potential revision of coherent-potential-approximation
alloy theory, Phys. Rev. B 45, 13234–13238.
[282] Nesbet, R.K. (1992). Full-potential multiple scattering theory without structure
constants, Phys. Rev. B 46, 9935–9939.
[283] Nesbet, R.K. (1996). Alternative density functional theory for atoms and molecules,
J. Phys. B 29, L173–L179.
[284] Nesbet, R.K. (1996). Nonadiabatic phase-matrix method for vibrational excitation
and dissociative attachment in electron–molecule scattering, Phys. Rev. A 54,
2899–2905.
[285] Nesbet, R.K. (1997). Fractional occupation numbers in density-functional theory,
Phys. Rev. A 56. 2665–2669.
[286] Nesbet, R.K. (1997). Local response model of the generalized polarization
potential, Phys. Rev. A 56. 2778–2783.
[287] Nesbet, R.K. (1997). Recent developments in multiple scattering theory and density
functional theory for molecules and solids, in Conceptual Perspectives in Quantum
Chemistry, eds. J.-L. Calais and E. Kryachko (Kluwer, Dordrecht), pp. 1–58.
[288] Nesbet, R.K. (1998). Kinetic energy in density-functional theory, Phys. Rev. A 58,
R12–R15.
[289] Nesbet, R.K. (1999). Exact exchange in linear-response theory, Phys. Rev. A 60,
R3343–R3346.
[290] Nesbet, R.K. (2000). Bound-free correlation in electron scattering by atoms and
molecules, Phys. Rev. A 62, 040701(R).
[291] Nesbet, R.K. (2001). Local potentials in independent-electron models, Int. J.
Quantum Chem. 81, 384–388.
[292] Nesbet, R.K. (2002). Orbital functional theory of linear response and excitation,
Int. J. Quantum Chem. 86, 342–346.
[293] Nesbet, R.K. and Colle, R. (1999). Does an exact local exchange potential exist?
J. Math. Chem. 26, 233–242.
[294] Nesbet, R.K. and Colle, R. (2000). Tests of the locality of exact Kohn–Sham
exchange potentials, Phys. Rev. A 61, 012503.
[295] Nesbet, R.K. and Grimm-Bosbach, T. (1993). Use of a discrete complete set of
vibrational basis functions in nonadiabatic theories of electron–molecule scattering,
J. Phys. B 26, L423–L426.
[296] Nesbet, R.K. and Sun, T. (1987). Self-consistent calculations using canonical
scaling in the linearized atomic-cell orbital method: Energy bands of fcc Cu, Phys.
Rev. B 36, 6351–6355.
[297] Nesbet, R.K., Mazevet, S. and Morrison, M.A. (2001). Procedure for correcting
variational R-matrix calculations for polarization response, Phys. Rev. A 64,
034702.
[298] Nesbet, R.K., Noble, C.J. and Morgan, L.A. (1986). Calculations of elastic electron
scattering by H2 for fixed nuclei, Phys. Rev. A 34, 2798–2807.
[299] Nesbet, R.K., Noble, C.J., Morgan, L.A. and Weatherford, C.A. (1984). Variational
R-matrix calculations of e− + H2 scattering using numerical asymptotic basis
functions, J. Phys. B 17, L891–L895.
218 References and bibliography
[300] Nestmann, B.M. and Peyerimhoff, S.D. (1990). Optimized Gaussian basis sets for
representation of continuum wave functions, J. Phys. B 22, L773–L777.
[301] Nestmann, B.M., Nesbet, R.K. and Peyerimhoff, S.D. (1991). A concept for
improving the efficiency of R-matrix calculations for electron–molecule scattering,
J. Phys. B 24, 5133–5149.
[302] Newton, R.G. (1966). Scattering Theory of Waves and Particles (Plenum,
New York).
[303] Noble, C.J. and Nesbet, R.K. (1984). CFASYM, a program for the calculation of the
asymptotic solutions of the coupled equations of electron collision theory, Comput.
Phys. Commun. 33, 399–411.
[304] Noether, E. (1918). Nachr. Ges. Wiss. Göttingen, 235.
[305] Norcross, D.W. and Padial, N.T. (1982). The multipole-extracted adiabatic nuclei
approximation for electron–molecule collisions, Phys. Rev. A 25, 226–238.
[306] Nordholm, S. and Bacskay, G. (1978). Generalizations of the finite-element and
R-matrix methods, J. Phys. B 11, 193–207.
[307] Norman, M.R. and Koelling, D.D. (1984). Towards a Kohn–Sham potential via the
optimized effective potential method, Phys. Rev. B 30, 5530–5540.
[308] Nozières, P. (1964). Theory of Fermi liquids, in 1962 Cargèse Lectures in
Theoretical Physics, ed. M. Lévy (Benjamin, New York), IV, 1–37,
[309] Nozières, P. (1964). Theory of Interacting Fermi Systems (Benjamin, New York).
[310] Oberoi, R.S. and Nesbet, R.K. (1973). Variational formulation of the R-matrix
method for multichannel scattering, Phys. Rev. A 8, 215–219.
[311] Oberoi, R.S. and Nesbet, R.K. (1973). Numerical asymptotic functions in
variational scattering theory, J. Comput. Phys. 12, 526–533.
[312] Oberoi, R.S. and Nesbet, R.K. (1974). Addendum to “Variational formulation of the
R matrix method for multichannel scattering”, Phys. Rev. A 9, 2804–2805.
[313] Oksyuk, Yu.D. (1966). Excitation of the rotational levels of diatomic molecules by
electron impact in the adiabatic approximation, Zh. Eksp. Teor. Fiz. 49, 1261–1273.
[Sov. Phys. JETP 22, 873–881 (1966)]
[314] O’Malley, T.F. (1966). Theory of dissociative attachment, Phys. Rev. 150, 14–29.
[315] O’Malley, T.F. (1967). Theory of dissociative attachment (Erratum), Phys. Rev.
156, 230.
[316] Ortiz, J.V. (1999). Toward an exact one-electron picture of chemical bonding, Adv.
Quantum Chem. 35, 33–52.
[317] Padial, N.T. and Norcross, D.W. (1984). Ro-vibrational excitation of HCl by
electron impact, Phys. Rev. A 29, 1590–1593.
[318] Panofsky, W.K.H. and Phillips, M. (1969). Classical Electricity and Magnetism
(Addison-Wesley, Reading, Mass.).
[319] Parr, R.G. (1983). Density functional theory, Ann. Rev. Phys. Chem. 34, 631–656.
[320] Parr, R.G. and Bartolotti, L.J. (1983). Some remarks on the density functional
theory of few-electron systems, J. Phys. Chem. 87, 2810–2815.
[321] Parr, R.G. and Yang, W. (1989). Density-Functional Theory of Atoms and
Molecules (Oxford University Press, New York).
[322] Pars, L.A. (1962). An Introduction to the Calculus of Variations (Wiley,
New York).
[323] Pauli, W. (1958). Theory of Relativity, tr. G. Field (Pergamon Press, New York).
[324] Payne, P.W. (1979). Density functionals in unrestricted Hartree–Fock theory,
J. Chem. Phys. 71, 490–496.
[325] Percival, I.C. and Seaton, M.J. (1957). The partial wave theory of
electron-hydrogen atom collisions, Proc. Camb. Phil. Soc. 53, 654–662.
References and bibliography 219
[326] Perdew, J.P. and Norman, M.R. (1982). Electron removal energies in Kohn–Sham
density-functional theory, Phys. Rev. B 26, 5445–5450.
[327] Perdew, J.P. and Norman, M.R. (1984). Reply to Comment on ‘Electron removal
energies in Kohn–Sham density-functional theory’, Phys. Rev. B 30, 3525–3526.
[328] Perdew, J.P., Parr, R.G., Levy, M. and Balduz, J.L. (1982). Density-functional
theory for fractional particle number: derivative discontinuities of the energy, Phys.
Rev. Lett. 49, 1691–1694.
[329] Petersilka, M., Gossmann, U.J. and Gross, E.K.U. (1996). Excitation energies from
time-dependent density-functional theory, Phys. Rev. Lett. 76, 1212–1215.
[330] Pfingst, K., Nestmann, B.M. and Peyerimhoff, S.D. (1995). Tailoring the R-matrix
approach for application to polyatomic molecules, in Computational Methods for
Electron–Molecule Collisions, eds. W.M. Huo and F. Gianturco (Plenum,
New York), pp. 293–308.
[331] Pokorski, S. (1987) Gauge Field Theories (Cambridge University Press,
Cambridge).
[332] Pulay, P. (1987). Analytical derivative methods in quantum chemistry, Adv. Chem.
Phys. 69, 241–286.
[333] Pulay, P. and Fogarasi, G. (1992). Geometry optimization in redundant internal
coordinates, J. Chem. Phys. 96, 2857–2860.
[334] Purcell, J.E. (1969). Nuclear structure calculations in the continuum-application to
neutron-carbon scattering, Phys. Rev. 185, 1279–1285.
[335] Quigg, C. (1983). Gauge Theories of the Strong, Weak, and Electromagnetic
Interactions (Benjamin-Cummings, New York).
[336] Ramond, P. (1989). Field Theory: A Modern Primer, 2nd edition (Addison-Wesley,
Redwood City, California).
[337] Lord Rayleigh (1892). Phil. Mag. 34, 481.
[338] Lord Rayleigh (1892). The Theory of Sound, 2nd edition, Vol. 2, p. 328, reprinted
by Dover Publications, New York (1976).
[339] Reed, M. and Simon, B. (1980). Methods of Modern Mathematical Physics I:
Functional Analysis (Academic Press, New York).
[340] Renton, P. (1990). Electroweak Interactions (Cambridge University Press, New
York).
[341] Rescigno, T.N., McCurdy, C.W., Orel, A.E. and Lengsfield, B.H., III (1995). The
complex Kohn variational method, in Computational Methods for
Electron–Molecule Collisions, eds. W.M. Huo and F.A. Gianturco (Plenum,
New York), pp. 1–44.
[342] Rescigno, T.N., McCurdy, C.W. and McKoy, V. (1974). Discrete basis set approach
to nonspherical scattering, Chem. Phys. Lett. 27, 401–404.
[343] Rescigno, T.N., McCurdy, C.W. and McKoy, V. (1974). Discrete basis set approach
to nonspherical scattering. II., Phys. Rev. A 10, 2240–2245.
[344] Rescigno, T.N., McCurdy, C.W. and McKoy, V. (1975). Low energy e− −H2 elastic
cross sections using discrete basis functions, Phys. Rev. A 11, 825–829.
[345] Rivers, R.J. (1987). Path Integral Methods in Quantum Field Theory (Cambridge
University Press, New York).
[346] Rohr, K. and Linder, F. (1976). Vibrational excitation of polar molecules by
electron impact I. Threshold resonances in HF and HCl, J. Phys. B 9, 2521–2537.
[347] Roothaan, C.C.J. (1951). New developments in molecular orbital theory, Rev. Mod.
Phys. 23, 69–89.
[348] Roothaan, C.C.J. (1960). Self-consistent field theory for open shells of electronic
systems, Rev. Mod. Phys. 32, 179–185.
220 References and bibliography
[349] Roothaan, C.C.J. and Bagus, P.S. (1963). Atomic self-consistent field calculations
by the expansion method, Methods Comput. Phys. 2, 47–94.
[350] Rubinow, S.I. (1955). Variational principle for scattering with tensor forces, Phys.
Rev. 98, 183–187.
[351] Runge, E. and Gross, E.K.U. (1984). Density-functional theory for time-dependent
systems, Phys. Rev. Lett. 52, 997–1000.
[352] Ryder, L.H. (1985). Quantum Field Theory (Cambridge University Press,
Cambridge).
[353] Salam, A. (1968). In Proceedings of the Eighth Nobel Symposium, ed. N. Svartholm
(Almqvist and Wiksell, Stockholm).
[354] Sarpal, B.J., Tennyson, J. and Morgan, L.A. (1991). Vibrationally resolved electron
HeH+ collisions using the nonadiabatic R-matrix method, J. Phys. B 24,
1851–1866.
[355] Sarpal, B.J., Tennyson, J. and Morgan, L.A. (1994). Dissociative recombination
without curve crossing: study of HeH+ , J. Phys. B 27, 5943–5953.
[356] Schlegel, H.B. (1982). Optimization of equilibrium geometries and transition
structures, J. Comput. Chem. 3, 214–218.
[357] Schlegel, H.B. (1987). Optimization of equilibrium geometries and transition
structures, Adv. Chem. Phys. 67, 249–286.
[358] Schlessinger, L. and Payne, G.L. (1974). Procedure for finding the scattering
solutions of the Schrödinger equation, Phys. Rev. A 10, 1559–1567.
[359] Schlosser, H. and Marcus, P. (1963). Composite wave variational method for
solution of the energy-band problem in solids, Phys. Rev. 131, 2529–2546.
[360] Schneider, B.I. (1975). R-matrix theory for electron–atom and electron–molecule
collisions using analytic basis set expansions. Chem. Phys. Lett. 31, 237–241.
[361] Schneider, B.I. (1976). Role of the Born–Oppenheimer approximation in the
vibrational excitation of molecules by electrons, Phys. Rev. A 14, 1923–1925.
[362] Schneider, B.I. (1995). An R-matrix approach to electron–molecule collisions, in
Computational Methods for Electron–Molecule Collisions, eds. W.M. Huo and
F. Gianturco (Plenum, New York), pp. 213–226.
[363] Schneider, B.I., LeDourneuf, M. and Burke, P.G. (1979). Theory of vibrational
excitation and dissociative attachment: an R-matrix approach, J. Phys. B 12,
L365–L369.
[364] Schneider, B.I., LeDourneuf, M. and Vo Ky Lan (1979). Resonant vibrational
excitation of N2 by low-energy electrons: an ab initio R-matrix calculation, Phys.
Rev. Lett. 43, 1926–1929.
[365] Schrödinger, E. (1926). Quantisiering als Eigenwertproblem, Ann. der Physik 81,
109–139.
[366] Schulman, L.S. (1981). Techniques and Applications of Path Integration (Wiley,
New York).
[367] Schulz, G.J. (1962). Vibrational excitation of nitrogen by electron impact, Phys.
Rev. 125, 229–232.
[368] Schulz, G.J. (1964). Vibrational excitation of N2 , CO, and H2 by electron impact,
Phys. Rev. 135, A988–A994.
[369] Schwartz, C. (1961). Variational calculations of scattering, Ann. Phys. (N.Y.) 16,
36–50.
[370] Schwartz, C. (1961). Electron scattering from hydrogen, Phys. Rev. 124,
1468–1471.
[371] Schwinger, J. (1947). A variational principle for scattering problems, Phys. Rev. 72,
742.
References and bibliography 221
[372] Schwinger, J. (1951). On the Green’s functions of quantized fields, Proc. Nat. Acad.
Sci. 37, 452–455.
[373] Schwinger, J. (1958). Quantum Electrodynamics (Dover, New York).
[374] Seaton, M.J. (1953). The Hartree–Fock equations for continuous states with
applications to electron excitation of the ground configuration terms of O I, Phil.
Trans. Roy. Soc. (London) A 245, 469–499.
[375] Seaton, M.J. (1966). A variational principle for the scattering phase matrix, Proc.
Phys. Soc. (London) 89, 469–470.
[376] Seaton, M.J. (1973). Close coupling, Comput. Phys. Commun. 6, 247–256.
[377] Seaton, M.J. (1983). Quantum defect theory, Rep. Prog. Phys. 46, 167–257.
[378] Seitz, L.I. (1955). Quantum Mechanics, 2nd edition (McGraw-Hill, New York).
[379] Shanno, D.F. (1970). Conditioning of quasi-Newton methods for function
minimisation, Math. Comput. 24, 647–656.
[380] Sharp, R.T. and Horton, G.K. (1953). A variational approach to the unipotential
many-electron problem, Phys. Rev. 90, 317.
[381] Shimamura, I. (1977). In The Physics of Electronic and Atomic Collisions, invited
papers, X ICPEAC, Paris (1977), (North-Holland, Amsterdam).
[382] Siegert, A.J.F. (1939). On the derivation of the dispersion formula for nuclear
reactions, Phys. Rev. 56, 750–752.
[383] Sinfailam, A.L. and Nesbet, R.K. (1972). Variational calculations of
electron–helium scattering, Phys. Rev. A 6, 2118–2125.
[384] Skriver, H.L. (1984). The LMTO Method (Springer-Verlag, Berlin).
[385] Slater, J.C. (1951). A simplification of the Hartree–Fock method, Phys. Rev. 81,
385–390.
[386] Slater, J.C. (1972). Statistical exchange-correlation in the self-consistent field,
Adv. Quantum Chem. 6, 1–92.
[387] Slater, J.C. (1974). The Self-Consistent Field for Molecules and Solids
(McGraw-Hill, New York).
[388] Slater, J.C. and Wood, J.H. (1971). Statistical exchange and the total energy of a
crystal, Int. J. Quantum Chem. Suppl. 4, 3–34.
[389] Smith, F.T. (1960). Lifetime matrix in collision theory, Phys. Rev. 118, 349–356.
[390] Sundermeyer, K. (1983). Constrained Dynamics (Springer-Verlag, Berlin).
[391] Swanson, M.S. (1992). Path Integrals and Quantum Processes (Academic Press,
New York).
[392] Symanzik, K. (1954). Schwinger’s functional in field theory, Z. Natürforschung 9a,
809–824.
[393] Synge, J.L. (1956). Relativity: the Special Theory (Interscience, New York).
[394] Szabo, A. and Ostlund, N. S. (1982). Modern Quantum Chemistry; Introduction to
Advanced Electronic Structure Theory (McGraw-Hill, New York).
[395] Takahashi, Y. (1957). On the generalized Ward identity, Nuovo Cimento 6, 371–375.
[396] Takatsuka, K. and McKoy, V. (1981). Extension of the Schwinger variational
principle beyond the static exchange approximation, Phys. Rev. A 24, 2473–2480.
[397] Takatsuka, K., Lucchese, R.R. and McKoy, V. (1981). Relationship between the
Schwinger and Kohn-type variational principles in scattering theory, Phys. Rev. A
24, 1812–1816.
[398] Talman, J.D. and Shadwick, W.F. (1976). Optimized effective atomic central
potential, Phys. Rev. A 14, 36–40.
[399] Taylor, J.T. (1983). Scattering Theory (Krieger, FL).
[400] Teller, E. (1962). On the stability of molecules in the Thomas–Fermi theory, Rev.
Mod. Phys. 34, 627–631.
222 References and bibliography
[426] Wheeler, J.A. and Feynman, R.P. (1945). Interaction with the absorber as the
mechanism of radiation, Rev. Mod. Phys. 17, 157–179.
[427] Wigner, E.P. (1948). On the behavior of cross sections near thresholds, Phys. Rev.
73, 1002–1009.
[428] Wigner, E.P. and Eisenbud, L. (1947). Higher angular momenta and long range
interaction in resonance reactions, Phys. Rev. 72, 29–41.
[429] Wigner, E.P. and Seitz, F. (1933). On the constitution of metallic sodium, Phys. Rev.
43, 804–810.
[430] Williams, A.R., Kübler, K. and Gelatt, C.D. (1979). Cohesive properties of metallic
compounds: Augmented-spherical-wave calculations, Phys. Rev. B 19,
6094–6118.
[431] Williams, A.R. and Morgan, J. van W. (1972). Multiple scattering by non-muffin-tin
potentials, J. Phys. C 5, L293–L298.
[432] Williams, A.R. and Morgan, J. van W. (1974). Multiple scattering by non-muffin-tin
potentials: general formulation, J. Phys. C 7, 37–60.
[433] Wolfenstein, L. (1999). Neutrino physics, Rev. Mod. Phys. 71, S140–S144.
[434] Yang, C.N. (1977). Magnetic monopoles, fiber bundles, and gauge fields, Ann. N.Y.
Acad. Sci. 294, 86–97.
[435] Yang, C.N. and Mills, R.L. (1954). Conservation of isotopic spin and isotopic gauge
invariance, Phys. Rev. 96, 191–195.
[436] Yourgrau, W. and Mandelstam, S. (1968). Variational Principles in Dynamics and
Quantum Theory, 3rd edition (Dover, New York).
[437] Zeller, R. (1987). Multiple-scattering solution of Schrödinger’s equation for
potentials of general shape, J. Phys. C 20, 2347–2360.
[438] Zeller, R. (1988). Empty-lattice test for non-muffin-tin multiple-scattering
equations, Phys. Rev. B 38, 5993–6002.
[439] Zerner, M.C. (1989). Analytic derivative methods and geometry optimization, in
Modern Quantum Chemistry, Introduction to Advanced Electronic Structure
Theory, eds. A. Szabo and N.S. Ostlund (McGraw-Hill, New York), pp. 437–458.
[440] Zhang, J.Z.H. and Miller, W.H. (1989). Quantum reactive scattering via the
S-matrix version of the Kohn variational principle: Differential and integral cross
sections for D + H2 → HD + H, J. Chem. Phys. 91, 1528–1547.
[441] Zhang, X.-G. and Butler, W.H. (1992). Simple cellular method for the exact
solution of the one-electron Schrödinger equation, Phys. Rev. Lett. 68, 3753–3756.
[442] Zhang, X.-G. and Butler, W.H. (1992). Multiple-scattering theory with a truncated
basis set, Phys. Rev. B 46, 7433–7447.
[443] Zhang, X.-G., Butler, W.H., MacLaren, J.M. and van Ek, J. (1994). Cellular
solutions for the Poisson equation in extended systems, Phys. Rev. B 49,
13383–13393.
[444] Zhang, X.-G., Butler, W.H., Nicholson, D.M. and Nesbet, R.K. (1992).
Green-function cellular method for the electronic structure of molecules and solids,
Phys. Rev. B 46, 15031–15039.
[445] Zhao, Q. and Parr, R.G. (1993). Constrained-search method to determine wave
functions from electronic densities, J. Chem. Phys. 98, 543–548.
[446] Zhao, Q., Morrison, R.C. and Parr, R.G. (1994). From electron densities to
Kohn–Sham kinetic energies, orbital energies, exchange-correlation potentials, and
exchange-correlation energies, Phys. Rev. A 50, 2138–2142.
[447] Zumbach, G. and Maschke, K. (1983). New approach to the calculation of density
functionals, Phys. Rev. A 28, 544–554.
[448] Zvijac, D.J., Heller, E.J. and Light, J.C. (1975). Variational correction to Wigner
R-matrix theory of scattering, J. Phys. B 8, 1016–1033.
Index
225
226 Index