public subroutine grad_en_hellman_feynman(basis, coord, zq, denab, de, logtol)
@brief Hellmann-Feynman force
@details Compute derivative contributions due to the Hamiltonian
operator change w.r.t. shifts of nuclei. The contribution
of the form is evaluated by Gauss-Rys quadrature.
This version handles spdfg and L shells.
Note
No relativistic methods available
@author Vladimir Mironov
@date Sep, 2018 Initial release
@param[in,out] denab density matrix in packed format, remains unchanged on return
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.