Avoiding fractional electrons in subsystem DFT based ab-initio molecular dynamics yields accurate models for liquid water and solvated OH radical
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102 (United States)
- CNR-ISTM: Institute of Molecular Sciences and Technologies, Milano (Italy)
In this work we achieve three milestones: (1) we present a subsystem DFT method capable of running ab-initio molecular dynamics simulations accurately and efficiently. (2) In order to rid the simulations of inter-molecular self-interaction error, we exploit the ability of semilocal frozen density embedding formulation of subsystem DFT to represent the total electron density as a sum of localized subsystem electron densities that are constrained to integrate to a preset, constant number of electrons; the success of the method relies on the fact that employed semilocal nonadditive kinetic energy functionals effectively cancel out errors in semilocal exchange–correlation potentials that are linked to static correlation effects and self-interaction. (3) We demonstrate this concept by simulating liquid water and solvated OH{sup •} radical. While the bulk of our simulations have been performed on a periodic box containing 64 independent water molecules for 52 ps, we also simulated a box containing 256 water molecules for 22 ps. The results show that, provided one employs an accurate nonadditive kinetic energy functional, the dynamics of liquid water and OH{sup •} radical are in semiquantitative agreement with experimental results or higher-level electronic structure calculations. Our assessments are based upon comparisons of radial and angular distribution functions as well as the diffusion coefficient of the liquid.
- OSTI ID:
- 22660741
- Journal Information:
- Journal of Chemical Physics, Journal Name: Journal of Chemical Physics Journal Issue: 23 Vol. 144; ISSN JCPSA6; ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Spin densities from subsystem density-functional theory: Assessment and application to a photosynthetic reaction center complex model
FDE-vdW: A van der Waals inclusive subsystem density-functional theory