Next: Volume vs. energy, volume Up: Theoretical Background Previous: Improved functional form for   Contents   Index

N.B. This document is no longer maintained, please visit our wiki.

## Forces

Within the finite temperature LDA forces are defined as the derivative of the generalized free energy. This quantity can be evaluated easily. The functional depends on the wavefunctions , the partial occupancies , and the positions of the ions . In this section we will shortly discuss the variational properties of the free energy and we will explain why we calculate the forces as a derivative of the free energy. The formulas given are very symbolic and we do not take into account any constraints on the occupation numbers or the wavefunctions. We denote the whole set of wavefunctions as and the set of partial occupancies as .

The electronic groundstate is determined by the variational property of the free energy i.e.

for arbitrary variations of and . We can rewrite the right hand side of this equation as

For arbitrary variations this quantity is zero only if and , leading to a system of equations which determines and at the electronic groundstate. We define the forces as derivatives of the free energy with respect to the ionic positions i.e.

force

At the groundstate the first two terms are zero and we can write

force

i.e. we can keep and fixed at their respective groundstate values and we have to calculate the partial derivative of the free energy with respect to the ionic positions only. This is relatively easy task.

Previously we have mentioned that the only physical quantity is the energy for . It is in principle possible to evaluate the derivatives of E( ) with respect to the ionic coordinates but this is not easy and requires additional computer time.

N.B. Requests for support are to be addressed to: vasp.materialphysik@univie.ac.at