eFF computes the energy of a collection of point charge nuclei and spherical Gaussian electrons ( ) as a function of the nuclear and electron positions plus the electron sizes. To obtain ground states and adiabatic dynamics, the electron positions and sizes are varied as a parametric function of the nuclear positions to minimize the overall energy of the system.
To obtain excited state dynamics, electrons are represented as wave packets:
Thus the average position of a wave packet obeys classical dynamics (Ehrenfest's theorem) with the addition that the size of the wave packet obeys classical dynamics as well. The radial effective mass factor depends on the dimension of the wave packet, and becomes for a 2D Gaussian and for a 1D Gaussian. Nuclei move via Ehrenfest dynamics [17,18].
Figure 5 shows a comparison of exact versus semiclassical Gaussian wave packet dynamics propagated in the one-dimensional double well potential . Over a short time interval, the two models match well. However, over longer periods the exact wave packet interferes with itself and eventually delocalizes over both wells, while the Gaussian wave packet bounces back and forth, with no signs of damping.
In systems where electrons are well-localized, semiclassical dynamics should be valid. However, quantum interference and tunneling effects will not be captured by the approximation.
From the equations of motion, it follows that if we define a total kinetic energy as
In eFF, the total energy is the sum of a Hartree product kinetic energy, a Hartree product electrostatic energy, and an antisymmetrization (Pauli) correction:
The eFF potential energy includes electron kinetic energies and electrostatic interactions between the individual electron densities, computed assuming the electrons act as independent particles.
The electron kinetic energies are a manifestation of the Heisenberg uncertainty principle, and vary as the inverse square of each electron's size. The finite size of a hydrogen atom, for example, is a balance between the electrostatic attraction of the electron to the nuclei, and the kinetic energy of the electron, which prevents the electron from collapsing onto the nucleus (Figure 6).
A Pauli exclusion term is needed because electrons are not independent particles - they are indistinguishable fermions, and the overall antisymmetry of the wavefunction upon interchange of any two particles gives rise to a Pauli repulsion between electrons.
The Pauli expression approximates the energy difference between the true antisymmetric wavefunction and the corresponding Hartree product wavefunction (Figure 7). It is derived by considering the difference between the antisymmetric and symmetric combinations of valence bond states, and assuming that kinetic energy differences predominate; analogous approximations have been reported elsewhere [11,12,6,9,4].
In eFF, the parameter corresponds to a dynamic electron mass and governs the overall time scale of electron motions. It does not affect the electronic kinetic energy which appears in , and so has no effect on the electron sizes and bond lengths in ground state molecules.
For electrons in harmonic potentials, the dynamic mass equals the true electron mass. However, in general this need not be the case, as in semiconductors, where the effective electron mass is , which may be much larger than the true electron mass; or in metals, where quasiparticles have an effective mass determined by many-body theory. Currently, we take to be an adjustable parameter.
In our studies on the Auger process in diamondoids, we found that core hole lifetimes were consistent with those obtained from experimental linewidth measurements only when . Aside from an overall time-scaling, changing did not change the basic events observed: core hole filling, ejection of a secondary electron, and excitation of neighboring electrons. The partitioning and distribution of electron energies were also unchanged.
It is also convenient to use a higher when simulating adiabatic dynamics, after Car-Parrinello dynamics, which often use an effective mass of a.u. In the case of the warm dense hydrogen studies , we found that changing did not affect the thermodynamic quantities averaged over time, and so used .
There are cases where it is essential to use a lower electron mass. For example, an ion passing through a solid excites electrons in the region it passes through, which slows the ion down. Here a rigorous separation of nuclear and electron time scales is needed to obtain the correct stopping power.