Search

Recently, the use of polarizable force fields in Molecular Dynamics simulations has been gaining importance, since they allow a better description of heterogeneous systems compared to simple point charges force fields. Among the various techniques developed in the last years the one based on polarizable point dipoles represents one of the most used. In this paper, we review the basic technical issues of the method, illustrating the way to implement intramolecular and intermolecular damping of the electrostatic interactions, either with and without the Ewald summation method. We also show how to reduce the computational overhead for evaluating the dipoles, introducing to the state-of-the-art methods: the extended Lagrangian method and the always stable predictor corrector method. Finally we discuss the importance of screening the electrostatic interactions at short range, defending this technique against simpler approximations usually made. We compare results of density functional theory and classical force field-based Molecular Dynamics simulations of chloride in water. [ABSTRACT FROM AUTHOR]

In this paper, double ionization energies (DIEs) of gas-phase atoms and molecules are calculated by energy difference method with density functional theory. To determine the best functional for double ionization energies, we first study 24 main group atoms in the second, third, and fourth periods. An approximation is used in which the electron density is first obtained from a density functional computation with the exchange-correlation potential Vxc known as statistical average of orbital potentials, after which the energy is computed from that density with 59 different exchange-correlation energy functionals Exc. For the 24 atoms, the two best Exc functional providing DIEs with average absolute deviation (AAD) of only 0.25 eV are the Perdew–Burke–Ernzerhof functional modified by Hammer et al. [Phys. Rev. B 59, 6413 (1999)] and one known as the Krieger–Chen–Iafrate–Savin functional modified by Krieger et al. (unpublished). Surprisingly, none of the 20 available hybrid functionals is among the top 15 functionals for the DIEs of the 24 atoms. A similar procedure is then applied to molecules, with opposite results: Only hybrid functionals are among the top 15 functionals for a selection of 29 molecules. The best Exc functional for the 29 molecules is found to be the Becke 1997 functional modified by Wilson et al. [J. Chem. Phys. 115, 9233 (2001)]. With that functional, the AAD from experiment for DIEs of 29 molecules is just under 0.5 eV. If the two suspected values for C2H2 and Fe(CO)5 are excluded, the AAD improves to 0.32 eV. Many other hybrid functionals perform almost as well. [ABSTRACT FROM AUTHOR]