dft-overview
|
dft-tex,
dft.html,
dft_biol.html,
dft_dft.html,
dft_fof.html,
dft_hamil.html,
dft_hk.html,
dft_intro.html,
dft_perf.html,
dft_refs.html,
dft_scfks.html,
dft_wave.html,
foot_motif.gif,
footnode.html,
image.gif,
img1.gif,
img10.gif,
img100.gif,
img101.gif,
img102.gif,
img103.gif,
img104.gif,
img105.gif,
img106.gif,
img107.gif,
img108.gif,
img109.gif,
img11.gif,
img110.gif,
img111.gif,
img112.gif,
img113.gif,
img114.gif,
img115.gif,
img116.gif,
img117.gif,
img118.gif,
img119.gif,
img12.gif,
img120.gif,
img121.gif,
img122.gif,
img123.gif,
img124.gif,
img125.gif,
img126.gif,
img127.gif,
img128.gif,
img129.gif,
img13.gif,
img130.gif,
img131.gif,
img132.gif,
img133.gif,
img134.gif,
img135.gif,
img136.gif,
img137.gif,
img138.gif,
img139.gif,
img14.gif,
img140.gif,
img141.gif,
img142.gif,
img143.gif,
img144.gif,
img145.gif,
img146.gif,
img147.gif,
img148.gif,
img149.gif,
img15.gif,
img150.gif,
img151.gif,
img152.gif,
img153.gif,
img154.gif,
img155.gif,
img156.gif,
img157.gif,
img158.gif,
img159.gif,
img16.gif,
img160.gif,
img161.gif,
img162.gif,
img163.gif,
img164.gif,
img165.gif,
img166.gif,
img167.gif,
img168.gif,
img169.gif,
img17.gif,
img170.gif,
img171.gif,
img172.gif,
img173.gif,
img174.gif,
img175.gif,
img176.gif,
img177.gif,
img178.gif,
img179.gif,
img18.gif,
img180.gif,
img181.gif,
img182.gif,
img183.gif,
img184.gif,
img185.gif,
img186.gif,
img187.gif,
img188.gif,
img189.gif,
img19.gif,
img190.gif,
img191.gif,
img192.gif,
img193.gif,
img194.gif,
img195.gif,
img196.gif,
img197.gif,
img198.gif,
img199.gif,
img2.gif,
img20.gif,
img200.gif,
img201.gif,
img202.gif,
img203.gif,
img204.gif,
img205.gif,
img206.gif,
img207.gif,
img208.gif,
img209.gif,
img21.gif,
img210.gif,
img211.gif,
img212.gif,
img213.gif,
img214.gif,
img215.gif,
img216.gif,
img217.gif,
img218.gif,
img219.gif,
img22.gif,
img220.gif,
img221.gif,
img222.gif,
img223.gif,
img224.gif,
img225.gif,
img226.gif,
img227.gif,
img228.gif,
img229.gif,
img23.gif,
img230.gif,
img231.gif,
img232.gif,
img233.gif,
img234.gif,
img235.gif,
img236.gif,
img237.gif,
img238.gif,
img239.gif,
img24.gif,
img240.gif,
img241.gif,
img242.gif,
img243.gif,
img244.gif,
img245.gif,
img246.gif,
img247.gif,
img248.gif,
img249.gif,
img25.gif,
img250.gif,
img251.gif,
img252.gif,
img253.gif,
img254.gif,
img255.gif,
img256.gif,
img26.gif,
img27.gif,
img28.gif,
img29.gif,
img3.gif,
img30.gif,
img31.gif,
img32.gif,
img33.gif,
img34.gif,
img35.gif,
img36.gif,
img37.gif,
img38.gif,
img39.gif,
img4.gif,
img40.gif,
img41.gif,
img42.gif,
img43.gif,
img44.gif,
img45.gif,
img46.gif,
img47.gif,
img48.gif,
img49.gif,
img5.gif,
img50.gif,
img51.gif,
img52.gif,
img53.gif,
img54.gif,
img55.gif,
img56.gif,
img57.gif,
img58.gif,
img59.gif,
img6.gif,
img60.gif,
img61.gif,
img62.gif,
img63.gif,
img64.gif,
img65.gif,
img66.gif,
img67.gif,
img68.gif,
img69.gif,
img7.gif,
img70.gif,
img71.gif,
img72.gif,
img73.gif,
img74.gif,
img75.gif,
img76.gif,
img77.gif,
img78.gif,
img79.gif,
img8.gif,
img80.gif,
img81.gif,
img82.gif,
img83.gif,
img84.gif,
img85.gif,
img86.gif,
img87.gif,
img88.gif,
img89.gif,
img9.gif,
img90.gif,
img91.gif,
img92.gif,
img93.gif,
img94.gif,
img95.gif,
img96.gif,
img97.gif,
img98.gif,
img99.gif
|
|
|
Molecular DFT
Introduction to Molecular Approaches of Density Functional Theory
Ohio Supercomputer Center, 1224 Kinnear Rd, Columbus, OH 43221-1153
Return to index
Hohenberg and Kohn theorems
The field of rigorous density functional theory was born in 1964 with
the publication of the Hohenberg and Kohn paper (1964).
They proved the following:
- I.
- Every observable of a stationary quantum mechanical system (including energy), can
be calculated, in principle exactly, from the ground-state density alone, i.e.,
every observable can be written as a functional of the ground-state density.
- II.
- The ground state density can be calculated, in principle exactly,
using the variational method involving only density,
The original theorems refer to the time independent (stationary) ground state, but
are being extended to excited states and time-dependent potentials (for review see, e.g.,
Gross & Kurth, 1994).
How these theorems were derived? By quite an original logic. Within a Born-Oppenheimer
approximation, the ground state of the system of electrons is a result of positions of
nuclei. This is obvious, if we look at the hamiltonian in equation (20).
In this hamiltonian, the kinetic energy of electrons ( ) and the
electron-electron interaction ( ) ``adjust'' themselves to the external
(i.e., coming from nuclei) potential . Once the is in place,
everything else is, including electron density, which simply adjusts itself to give the
lowest possible total energy of the system. The external potential is the
only variable term in this equation, and everything else depends indirectly on it.
Hohenberg and Kohn posed a more interesting question, which is quite opposite to
the traditional logic: Is uniquely determined
from the knowledge of electron density ? Can we find out (in principle,
though it need not be easy) where and what the nuclei are, if we know the density
in the ground state? Is there a precise mapping from to ? The answer to
this question is: Yes!. Actually the mapping is accurate within a constant, which would not
change anything, since Schrödinger equations with and yields
exactly the same eigenfunctions, i.e., states (it is easy to prove based on the linear property
of the hamiltonian), and the energies will be simply elevated by the value of this const. Note that
all energies are known only within some constant, which establishes the frame of reference (e.g.,
we do not include electron-Mars gravitational attraction within most calculations).
Why was this question so important? Because, if this is true, the knowledge
of density provides total information about the system, and formally if we know the density,
we know everything there is to known.
Since determines number of electrons, N:
and determines the , the knowledge of total density is as good, as
knowledge of , i.e., the wave function describing the state of the system.
They proved it through a contradiction:
- 1.
- Assume that we have an exact ground state density
,
- 2.
- Assume that the ground state is nondegenerate (i.e., there is only one wave function
for this ground state (though HK theorems can be easily
extended for degenerate ground
states, see, e.g., Dreizler & Gross, 1990; Kohn, 1985),
- 3.
- Assume that for the density
there are two possible external potentials:
and , which obviously produce two different hamiltonians:
and , respectively. They obviously produce two different
wave functions for the ground state: and , respectively. They correspond to
energies:
and ,
respectively.
- 4.
- Now, let us calculate the expectation value of energy
for the
with the
hamiltonian and use variational theorem:
- 5.
- Now let us calculate the expectation value of energy for the
with the
hamiltonian and use varational theorem:
- 6.
- By adding equations (61) and (62) by sides we obtain:
and it leads to a contradiction.
Since now, we know that determines N and , it also
determines all properties of the ground state, including the kinetic
energy of electrons
and energy of interaction among electrons , i.e., the total ground state energy
is a functional of density with the following
components :
Additionally, HK grouped together all functionals
which are secondary (i.e., which are responses)
to the :
The functional operates only on density and is universal, i.e., its form does not depend on the particular
system under consideration (note that N-representable densities integrate to N, and the information about
the number of electrons can be easily obtained from the density itself).
The second HK theorem provides variational principle in electron density representation .
For a trial density
such that and for which
,
where is the energy functional. In other words, if some density represents the correct number
of electrons N, the total energy calculated from this density cannot be lower than the true energy of the
ground state.
As to the necessary conditions for this theorem, there is still some controversy concerning
the, so called, representability of density. The N-representability, i.e., the fact that the trial has to
sum up to N electrons is easy to achieve by simple rescaling. It is
automatically insured if can be mapped to some wave function
(for further discussion see: Parr & Yang, 1989; Gilbert, 1975; Lieb, 1982;
and Harriman, 1980).
Assuring that the trial density has also -representability
(usually denoted in the literature as v-representability) is not that easy. Levy (1982) and Lieb (1983)
have shown, that there are some ``reasonable'' trial densities, which are not the ground state densities for any
possible potential, i.e., they do not map to any external potential. Such densities do not correspond therefore
to any ground state, and their optimization will not lead to a ground state.
Moreover, during energy minimization, we may take a wrong turn, and get
stuck into some non v-representable density and never be able to converged
to a physically relevant ground state density.
For an interesting discussion, see
Hohenberg et al. (1990). Assuming that we restrict ourselved only to
trial densities which are both N and v representable, the variational principle for density is easly proven, since
each trial density defines a hamiltonian . From the hamiltonian
we can derive the corresponding wave function for the ground state represented by this hamiltonian. And
according to the traditional variational principle, this wave function will not be a ground state
for the hamiltonian of the real system :
where is the true ground state density of the real system.
The condition of minimum for the energy functional:
needs to be constrained
by the N-representability of density which is optimized .
The Lagrange's method of undetermined multipliers is a very
convenient approach for the constrained minimization
problems. In this method
we represent constraints in such a way that their value is
exactly zero when they are satisfied. In our case, the
N representability constraint can be represented as:
These constraints are then multiplied by an undetermined
constants and added to a minimized function or functional.
where is yet undetermined Lagrange multiplier.
Now, we look for the minimum of this expression by requiring that its
differential is equal to zero (a necessary condition of minimum).
Solving this differential equation will provide us with a prescription of
finding a minimum which satisfies the constraint. In our case it leads to:
since and N are constants. Using the definition of the
differential of the functional (see, e.g., Parr & Yang, 1989):
and the fact that differential and integral signs may be interchanged, we obtain
Since integration runs over the same variable and has the same limits, we can
write both expressions under the same integral:
which provides the condition for constrained minimisation and defines the value
of the Lagrange multiplier at minium. It is also expressed here via external potential from equation (65):
Density functional theory
gives a firm definition of the chemical potential , and leads to
several important general conclusions. For review, please refer to Parr & Yang (1989), chapters 4 and 5.
Return to index
Next Section
|