7.4 Density Functional Molecular Dynamics (DFMD)

Born-Oppenheimer approximation:

- the heavy atomic cores (``ions'') consisting of nucleus and inner electronic shells are treated as classical particles

- the valence and conduction electrons are represented by wave functions which assume the configuration of least energy in the momentary field created by the ions and all other valence and conduction electrons.

Let denote the wave functions of the electrons; the electron density at position is then

(7.60) |

Total energy of the electrons:

(7.62) | |||

(7.63) | |||

(7.64) | |||

(7.65) |

Approximation for : usually

- Expand
in plane waves,

(7.66) - Find that set of expansion coefficients
which minimizes the energy functional
7.61; as an additional requirement,
maintain the orthonormality condition

Apply variational calculus:

Let denote an abstract ``mass''
assigned to each electronic wave function .
The ``kinetic energy'' due to the change of in time is
then

(7.68) |

Now introduce a *Lagrangian function*

(7.69) |

describing the ``motion''of . Here is the ionic mass, and the Lagrange multipliers allow for the conditions 7.67. From the Lagrangian the ``equations of motion''are derived:

The second equation describes the classical dynamics of the ions. The first equation represents the abstract ``motion'' in the space of the electronic degrees of freedom.

By keeping the ``kinetic energy'' of this motion small we can constrain the -subsystem always close to the momentary minimum of the energy surface defined by the slowly varying ionic configuration.

To control the temperature of the

- rescale all from time to time, or

- use one of the thermostats available from statistical-mechanical simulation; see Chapter 6.

The electronic degrees of freedom should adjust quickly to the varying energy landscape mass small compared to the ionic masses; e. g. atomic mass unit.

- amorphous silicon [CAR 85]

- lithium [WENTZCOVICH 91]

microclusters of alkali metals [VITEK 89]

- molten carbon [GALLI 90B]

- ionic melts [GALLI 90A]

For current applications, see the web or reviews such as [VITEK 89,MAKRI 99,MAZZONE 99,OHNO 99].

See also: