Force Feynman

In this minimal END approximation, the electronic basis functions are centered on the average nuclear positions, which are dynamical variables. In the limit of classical nuclei, these are conventional basis functions used in moleculai electronic structure theoiy, and they follow the dynamically changing nuclear positions. As can be seen from the equations of motion discussed above the evolution of the nuclear positions and momenta is governed by Newton-like equations with Hellman-Feynman forces, while the electronic dynamical variables are complex molecular orbital coefficients that follow equations that look like those of the time-dependent Hartree-Fock (TDHF) approximation [24]. The coupling terms in the dynamical metric are the well-known nonadiabatic terms due to the fact that the basis moves with the dynamically changing nuclear positions. 

Note that the exact adiabatic functions are used on the right-hand side, which in practical calculations must be evaluated by the full derivative on the left of Eq. (24) rather than the Hellmann-Feynman forces. This forai has the advantage that the R dependence of the coefficients, c, does not have to be considered. Using the relationship Eq. (78) for the off-diagonal matrix elements of the right-hand side then leads directly to... 

Importantly for direct dynamics calculations, analytic gradients for MCSCF methods [124-126] are available in many standard quantum chemistiy packages. This is a big advantage as numerical gradients require many evaluations of the wave function. The evaluation of the non-Hellmann-Feynman forces is the major effort, and requires the solution of what are termed the coupled-perturbed MCSCF (CP-MCSCF) equations. The large memory requirements of these equations can be bypassed if a direct method is used [233]. Modem computer architectures and codes then make the evaluation of first and second derivatives relatively straightforward in this theoretical framework. 

The intriguing point about the second set of equations is that q is now kept constant. Thus the vector ip evolves according to a time-dependent Schrddinger equation with time-independent Hamilton operator H[q) and the update of the classical momentum p is obtained by integrating the Hellmann-Feynman forces [3] acting on the classical particles along the computed ip t) (plus a constant update due to the purely classical force field). 

Unfortunately, this only holds for the exact wavefunction and certain other types ol leavefuiiction (such as at the Hartree-Fock limit). Moreover, even though the Hellmarm-Feynman forces are much easier to calculate they are very unreliable, even for accurate wavefunctions, giving rise to spurious forces (often referred to as Pulay forces [Pulay l )S7]). 

The first term is the Hellmann-Feynman force and the second is the wave function response. The latter now contains contributions both from a change in basis functions and MO coefficients. 

The central of these is recognized as the Hellmann-Feynman force. The two-electron... 

The central term is again the Hellmann-Feynman force, which vanishes since the two-electron operator g is independent of the nuclear positions. 

Figure 1 Convergence of the total energy and of the Hellmann-Feynman forces for ensembles of paramagnetic Fe atoms with 4 to 32 atoms. Part (a) shows the results of non-selfconsistent calculations performed with a fixed potential, part (b) the results of selfconsistent calculations. Full lines represent the RMM-DIIS (iterative diagonal-ization) results, broken lines the CGa (total-energy minimization) calculations. (4. text.

The Hellmann-Feynman forces are also sensitive to the effect of moving ions on the basis set (pj) of the electronic wave function (, = This... 

The electrostatic Hellmann-Feynman theorem states that for an exact electron wave function, and also of the Hartree-Fock wave function, the total quantum-mechanical force on an atomic nucleus is the same as that exerted classically by the electron density and the other nuclei in the system (Feynman 1939, Levine 1983). The theorem thus implies that the forces on the nuclei are fully determined once the charge distribution is known. As the forces on the nuclei must vanish for a nuclear configuration which is in equilibrium, a constraint may be introduced in the X-ray refinement procedure to ensure that the Hellmann-Feynman force balance is obeyed (Schwarzenbach and Lewis 1982). 

Here, the one-electron operator Zj (rj- R /lq-Ral3 is referred to as the Hellmann-Feynman force operator it is the derivative of the Hamiltonian with respect to displacement of center-a in the x, y, or z direction. 

Here ojaa>(R) = (Ea(R) — Ea/(R))/h and iLaa/ is the Liouville operator that describes classical evolution determined by the mean of the Hellmann-Feynman forces corresponding to adiabatic states a and a, ... 

Grochowski P, Lesyng B (2003) Extended Hellmann-Feynman Forces, Canonical Representations, and Exponential Propagators in the Mixed Quantum-Classical Molecular Dynamics. J. Chem. Phys. 119 11541-11555... 

The time-step of 0.5 fs is used to simulate the dynamic system to 4.0 ps. The temperature of 300 K is used throughout the simulations. The MD simulations are performed using the Nose-Hoover thermostat for temperature control. The Hellmann-Feynman forces acting on the atoms are calculated from the ground-state electronic energies at each time step and are subsequently used in the integration of Newton s equation of motion. 

Fhfb is the Hellman-Feynman force plus the correction for the orbital basis set dependence on the nuclear coordinate, x. The term is specific to the present LCGTO-LSDF implementation. Equation 10 is equivalent to... 

It should be pointed out that Schwarz (20),using double perturbation theory,has demonstrated that it is possible to rationalize the relativistic bond length contraction in terms of the attractive Hellmann-Feynman force due to the relativistic change in electron density.In such an approach it would be necessary to analyze and get a physical picture of the relevant density changes... 

---