Premium
Electrostatically embedded molecules‐in‐molecules approach and its application to molecular clusters
Author(s) -
Tripathy Vikrant,
Saha Arjun,
Raghavachari Krishnan
Publication year - 2021
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.26492
Subject(s) - fragmentation (computing) , embedding , electrostatics , molecule , chemistry , range (aeronautics) , charge (physics) , electrostatic interaction , atomic physics , quantum chemistry , computational chemistry , molecular physics , physics , chemical physics , materials science , quantum mechanics , computer science , operating system , organic chemistry , artificial intelligence , composite material , supramolecular chemistry
We report the application of our fragment‐based quantum chemistry model MIM ( M olecules‐ I n‐ M olecules) with electrostatic embedding. The method is termed “EE‐MIM ( E lectrostatically E mbedded M olecules‐ I n‐ M olecules)” and accounts for the missing electrostatic interactions in the subsystems resulting from fragmentation. Point charges placed at the atomic positions are used to represent the interaction of each subsystem with the rest of the molecule with minimal increase in the computational cost. We have carefully calibrated this model on a range of different sizes of clusters containing up to 57 water molecules. The fragmentation methods have been applied with the goal of reproducing the unfragmented total energy at the MP2/6‐311G(d,p) level. Comparative analysis has been carried out between MIM and EE‐MIM to gauge the impact of electrostatic embedding. Performance of several different parameters such as the type of charge and levels of fragmentation are analyzed for the prediction of absolute energies. The use of background charges in subsystem calculations improves the performance of both one‐ and two‐layer MIM while it is noticeably important in the case of one‐layer MIM. Embedded charges for two‐layer MIM are obtained from a full system calculation at the low‐level. For one‐layer MIM, in the absence of a full system calculation, two different types of embedded charges, namely, Geometry dependent (GD) and geometry independent (GI) charges, are used. A self‐consistent procedure is employed to obtain GD charges. We have further tested our method on challenging charged systems with stronger intermolecular interactions, namely, protonated ammonia clusters (containing up to 30 ammonia molecules). The observations are similar to water clusters with improved performance using embedded charges. Overall, the performance of NPA charges as embedded charges is found to be the best.