Publication date: Mar 16, 2022
Non-covalent bonding patterns are commonly harvested as a design principle in the field of catalysis, supramolecular chemistry, and functional materials to name a few. Yet, their computational description generally neglects finite temperature and environment effects, which promote competing interactions and alter their static gas-phase properties. Recently, neural network potentials (NNPs) trained on Density Functional Theory (DFT) data have become increasingly popular to simulate molecular phenomena in condensed phase with an accuracy comparable to ab initio methods. To date, most applications have centered on solid-state materials or fairly simple molecules made of a limited number of elements. Herein, we focus on the persistence and strength of chalcogen bonds involving benzotelluradiazole in condensed phase. While the tellurium-containing heteroaromatic molecules are known to exhibit pronounced interactions with anions and lone pairs of different atoms, the relevance of competing intermolecular interactions, notably with the solvent, is complicated to monitor experimentally but also challenging to model at an accurate electronic structure level. Here, we train direct and baselined NNPs to reproduce hybrid DFT energies and forces in order to identify what are the most prevalent non-covalent interactions occurring in a solute-Cl-THF mixture. The simulations in explicit solvent highlight the clear competition with chalcogen bonds formed with the solvent and the short-range directionality of the interaction with direct consequences for the molecular properties in the solution. The comparison with other potentials (e.g., AMOEBA, direct NNP, and continuum solvent model) also demonstrates that baselined NNPs offer a reliable picture of the non--covalent interaction interplay occurring in solution.
No Explore or Discover sections associated with this archive record.
File name | Size | Description |
---|---|---|
input.xml
MD5md5:91eb0f8ec9be129a246f6270a5d1fd8a
|
3.4 KiB | Example input to i-PI to run MD using MTS approach with baselined and direct NNP |
lmp1.in
MD5md5:9a2d4d5485b2d5cc78481da3895b7026
|
2.3 KiB | Example LAMMPS input to run MD with NNP using i-PI |
model_baselined.tar.gz
MD5md5:54c66b0695a147dadc8a4d7748c7723e
|
117.2 MiB | Parameters for the baselined NNPs in n2p2 format |
model_direct.tar.gz
MD5md5:2df14da4885d249d681a62f2b2336926
|
779.0 KiB | Parameters for the direct NNPs in n2p2 format |
pbe.cp2k
MD5md5:6575bad7f35df5bc7dd59e1ba1c106d8
|
1.7 KiB | Input for PBE computations in cp2k |
pbe0.cp2k
MD5md5:88396e7112a32ce8b98d3aaab4be5677
|
3.3 KiB | Input for PBE0 computations in cp2k |
plumed.in
MD5md5:ffd2bebe722010fcc6ce47c14a95aaa8
|
157 Bytes | Example PLUMED input for umbrella sampling |
2022.42 (version v1) [This version] | Mar 16, 2022 | DOI10.24435/materialscloud:90-vd |