Journal
PHYSICAL CHEMISTRY CHEMICAL PHYSICS
Volume 22, Issue 7, Pages 3855-3866Publisher
ROYAL SOC CHEMISTRY
DOI: 10.1039/c9cp06792b
Keywords
-
Funding
- Australian Research Council [DE160100807]
- Japan Society for the Promotion of Science Program for Advancing Strategic International Networks
- National Institutes of Health [R01 GM135392]
- Oklahoma Center for the Advancement of Science and Technology [HR18-130]
Ask authors/readers for more resources
In this paper, the performance of a wide range of DFT methods is assessed for the calculation of interaction energies of thermal clusters of a solute in water. Three different charge states (neutral, proton transfer transition state and zwitterion) of glycine were solvated by 1 to 40 water molecules as sampled from molecular dynamics simulations. While some ab initio composite methods that employ insufficiently large basis sets incurred significant errors even for a cluster containing only 5 water molecules relative to the W1X-2 benchmark, the DLPNO-CCSD(T)/CBS and DSD-PBEP86 (triple zeta basis set) levels of theory predicted very accurate interaction energies. These levels of theory were used to benchmark the performance of 16 density functionals from different rungs of Jacob's Ladder. Of the Rung 4 functionals examined, the omega B97M-V and omega B97X-V functionals stood out for predicting absolute interaction energies in 40-water clusters with mean absolute deviations (MAD) similar to 4 kJ mol(-1). The B3LYP-D3(BJ) functional performed exceptionally well with a MAD similar to 1.7 kJ mol(-1) and is the overall best performing method. Calculations of relative interaction energies allow for cancellation of systematic errors, including basis set truncation and superposition errors, and the omega B97M-V and B3LYP-D3(BJ) double zeta basis set calculations yielded relative interaction energies that are within similar to 3 kJ mol(-1) of the benchmark. The ONIOM approximation provides another strategy for accelerating the calculation of accurate absolute interaction energies provided that the calculations have converged with respect to the size of the high-level-layer.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available