Introducing KICK-MEP: exploring potential energy surfaces in systems with significant non-covalent interactions

Context

Exploring potential energy surfaces (PES) is fundamental in computational chemistry, as it provides insights into the relationship between molecular energy, geometry, and chemical reactivity. We introduce Kick-MEP, a hybrid method for exploring the PES of atomic and molecular clusters, particularly those dominated by non-covalent interactions. Kick-MEP computes the Coulomb integral between the maximum and minimum electrostatic potential values on a 0.001 a.u. electron density isosurface for two interacting fragments. This approach efficiently estimates interaction energies and selects low-energy configurations at reduced computational cost. Kick-MEP was evaluated on silicon-lithium clusters, water clusters, and thymol encapsulated within Cucurbit[7]uril, consistently identifying the lowest energy structures, including global minima and relevant local minima.

Methods

Kick-MEP generates an initial population of molecular structures using the stochastic Kick algorithm, which combines two molecular fragments (A and B). The molecular electrostatic potential (MEP) values on a 0.001 a.u. electron density isosurface for each fragment are used to compute the Coulomb integral between them. Structures with the lowest Coulomb integral are selected and refined through gradient-based optimization and DFT calculations at the PBE0-D3/Def2-TZVP level. Molecular docking simulations for the thymol-Cucurbit[7]uril complex using AutoDock Vina were performed for benchmarking. Kick-MEP was validated across different molecular systems, demonstrating its effectiveness in identifying the lowest energy structures, including global minima and relevant local minima, while maintaining a low computational cost.