A multilevel approach is presented to assess the ability of several popular dispersion corrected density functionals (M06-2X, CAM-B3LYP-D3, BLYP-D3, and B3LYP-D3) to reliably describe two-body interaction potential energy surfaces (IPESs). To this end, the automated Picky procedure ( Cacelli et al. J. Comput. Chem. 2012 , 33 , 1055 ) was exploited, which consists in parametrizing specific intermolecular force fields through an iterative approach, based on the comparison with quantum mechanical data. For each of the tested functionals, the resulting force field was employed in classical Monte Carlo and Molecular Dynamics simulations, performed on systems of up to 1000 molecules in ambient conditions, to calculate a number of condensed phase properties. The comparison of the resulting structural and dynamic properties with experimental data allows us to assess the quality of each IPES and, consequently, even the quality of the DFT functionals. The methodology is tested on the benzene dimer, commonly used as a benchmark molecule, a prototype of aromatic interactions. The best results were obtained with the CAM-B3LYP-D3 functional. Besides assessing the reliability of DFT functionals in describing aromatic IPESs, this work provides a further step toward a robust protocol for the derivation of sound force field parameters from quantum mechanical data. This method can be relevant in all those cases where standard force fields fail in giving accurate predictions.

An automated protocol is proposed and validated, which integrates accurate quantum mechanical calculations with classical numerical simulations. Intermolecular force fields, (FF) suitable for molecular dynamics (MD) and Monte Carlo simulations, are parameterized through a novel iterative approach, fully based on quantum mechanical data, which has been automated and coded into the PICKY software, here presented. The whole procedure is tested and validated for pyridine, whose bulk phase, described through MD simulations performed with the specifically parameterized FF, is characterized by computing several of its thermodynamic, structural, and transport properties, comparing them with their experimental counterparts.

Comput., 2018, 14, 543-556) to evaluate a large number of dimer interaction energies. The resulting quantum mechanically derived FFs are then used in extensive molecular dynamics simulations, in order to evaluate a number of thermodynamic, structural, and dynamic properties of the heterocycle's gas and liquid phases. The comparison with the available experimental data is good and furnishes a validation of the presented approach, which can be confidently exploited for the design of novel and more complex materials.

Noncovalent interactions between homodimers of several aromatic heterocycles (pyrrole, furan, thiophene, pyridine, pyridazine, pyrimidine, and pyrazine) are investigated at the ab initio level, employing the Möller-Plesset second-order perturbation theory, coupled with small Gaussian basis sets (6-31G* and 6-31G**) with specifically tuned polarization exponents. The latter are modified using a systematic and automated procedure, the MP2 approach, based on a comparison with high level CCSD(T) calculations extrapolated to a complete basis set. The MP2 results achieved with the modified 6-31G** basis set show an excellent agreement with CCSD(T)/CBS reference energies, with a standard deviation less than 0.3 kcal/mol. Exploiting its low computational cost, the MP2 approach is then used to explore sections of the intermolecular energy of the considered homodimers, with the aim of rationalizing the results. It is found that the direct electrostatic interaction between the monomers electron clouds is at the origin of some observed features, and in many cases multipoles higher than dipole play a relevant role, although often the interplay with other contributions to the noncovalent forces (as for instance induction, π-π or XH-π interactions) makes a simple rationalization rather difficult.

A robust and automated protocol for the derivation of sound force field parameters, suitable for condensed-phase classical simulations, is here tested and validated on several halogenated hydrocarbons, a class of compounds for which standard force fields have often been reported to deliver rather inaccurate performances. The major strength of the proposed protocol is that all of the parameters are derived only from first principles because all of the information required is retrieved from quantum mechanical data, purposely computed for the investigated molecule. This a priori parametrization is carried out separately for the intra- and intermolecular contributions to the force fields, respectively exploiting the Joyce and Picky programs, previously developed in our group. To avoid high computational costs, all quantum mechanical calculations were performed exploiting the density functional theory. Because the choice of the functional is known to be crucial for the description of the intermolecular interactions, a specific procedure is proposed, which allows for a reliable benchmark of different functionals against higher-level data. The intramolecular and intermolecular contribution are eventually joined together, and the resulting quantum mechanically derived force field is thereafter employed in lengthy molecular dynamics simulations to compute several thermodynamic properties that characterize the resulting bulk phase. The accuracy of the proposed parametrization protocol is finally validated by comparing the computed macroscopic observables with the available experimental counterparts. It is found that, on average, the proposed approach is capable of yielding a consistent description of the investigated set, often outperforming the literature standard force fields, or at least delivering results of similar accuracy.

1-n-Butyl-2,3-dimethylimidazolium (BMMI) ionic liquids (ILs) associated with different anions undergo H/D exchange preferentially at 2-Me group of the imidazolium in deuterated solvents. This process is mainly related to the existence of ion pairs rather than the anion basicity. The H/D exchange occurs in solvents (CDCl3 and MeCN for instance) in which intimate contact ion pairs are present and the anion possesses a labile H in its structure, such as hydrogen carbonate and prolinate. In D2 O, separated ion pairs are formed and the H/D exchange does not occur. A plausible catalytic cycle is that the IL behaves as a neutral base in the course of all H/D exchange processes. NMR experiments, density functional calculations, and molecular dynamics simulations corroborate these hypotheses.

The reliability of molecular dynamics (MD) simulations in predicting macroscopic properties of complex fluids and soft materials, such as liquid crystals, colloidal suspensions, or polymers, relies on the accuracy of the adopted force field (FF). We present an automated protocol to derive specific and accurate FFs, fully based on ab initio quantum mechanical (QM) data. The integration of the JOYCE and PICKY procedures, recently proposed by our group to provide an accurate description of simple liquids, is here extended to larger molecules, capable of exhibiting more complex fluid phases. While the standard JOYCE protocol is employed to parameterize the intramolecular FF term, a new automated procedure is here proposed to handle the computational cost of the QM calculations required for the parameterization of the intermolecular FF term. The latter is thus obtained by integrating the old PICKY procedure with a fragmentation reconstruction method (FRM) that allows for a reliable, yet computationally feasible sampling of the intermolecular energy surface at the QM level. The whole FF parameterization protocol is tested on a benchmark liquid crystal, and the performances of the resulting quantum mechanically derived (QMD) FF were compared with those delivered by a general-purpose, transferable one, and by the third, "hybrid" FF, where only the bonded terms were refined against QM data. Lengthy atomistic MD simulations are carried out with each FF on extended 5CB systems in both isotropic and nematic phases, eventually validating the proposed protocol by comparing the resulting macroscopic properties with other computational models and with experiments. The QMD-FF yields the best performances, reproducing both phases in the correct range of temperatures and well describing their structure, dynamics, and thermodynamic properties, thus providing a clear protocol that may be explored to predict such properties on other complex fluids or soft materials.

It is well known that the macroscopic physico-chemical properties of ionic liquids (ILs) are influenced by the presence of water that strongly interferes with the supramolecular organization of these fluids. However, little is known about the function of water traces within this confined space and restricted ionic environments, i.e. between cations and anions. Using specially designed ILs namely 1,2,3-trimethyl-1H-imidazol-3-ium imidazol-1-ide (MMMI·Im) and 3-n-butyl-1,2-dimethyl-1H-imidazol-3-ium imidazol-1-ide (BMMI·Im), the structure and function of water have been determined in condensed, solution and gas phases by X-ray diffraction studies, NMR, molecular dynamics simulations (MDS) and DFT calculations. In the solid state the water molecule is trapped inside the ionic network (constituted of contact ion pairs formed by π(+)-π(-) interaction) through strong H-bonds involving the water hydrogens and the nitrogens of two imidazolate anions forming a guest@host supramolecular structure. A similar structural arrangement was corroborated by DFT calculations and MDS. The presence of a guest@host species (H2O@ILpair) is maintained to a great extent even in solution as detected by (1)H-(1)H NOESY-experiments of the ILs dissolved in solvents with low and high dielectric constants. This confined water catalyses the H/D exchange with other substrates containing acidic-H such as chloroform.