Dataset for "Comparing the Performances of Force Fields in Conformational Searching of Hydrogen Bond-Donating Catalysts"

Data collection method:

Conformational searches were performed with the force fields available in Schrödinger’s MacroModel v12.6 (OPLS3e, OPLS-2005, MMFF, MMFFs, AMBER*, OPLS, MM2* and MM3*) on each of 20 molecules. Aside from the force field potentials, the other settings used in the conformational searches were as follows: solvent was set to “None” (gas-phase); the energy minimization method was set to “Polak-Ribiere Conjugate Gradient” (PRCG), converging on the gradient with convergence threshold set to 0.001; the search method was set to “Mixed torsional/Low-mode sampling” with the energy window for saving structures set to 50 kJ mol-1 and the maximum atom deviation for structures to be considered different conformers was set to 0.25 Å. All DFT calculations were performed using the Gaussian 16, Revision A.03 software. Geometry optimizations were performed on all conformers of all molecules at the M06-2X/6-31G(d) level of theory in the IEFPCM(benzene) solvent model. Once all of the conformer structures had reached minima and were stationary points, single-point energies were calculated at the M06-2X/def2-TZVPP IEFPCM(benzene) level of theory. This level of theory was chosen since it has previously been used with success in previous reaction modelling studies.

Technical details and requirements:

Schrödinger MacroModel v12.6; Gaussian 16, Revision A.03

Additional information:

Each of the 20 molecules in the dataset was given a nickname that is used for the calculation files. The molecule numbers (used in the publication) and 2D molecular structures that these nicknames correspond to can be found from the images "molecule_names_and_structures_1.jpg" and "molecule_names_and_structures_2.jpg". Within each molecule directory, there is one directory for each force field that was able to conformationally search that molecule (format: <mol_name>/<mol_name_FF>). Within the <mol_name_FF> directories are the .mol2 files that correspond to the force field structures from the conformational searches. As well as the force field structures, the <mol_name_FF> directories each contain a "mo62x_6-31Gd" directory that contain the M06-2X/6-31G(d) optimised minima as Gaussian .out files and the M06-2X/def2-TZVPP//M06-2X/6-31G(d) single-point energy calculations as "_SPEopt.out" files. The "csvs" directory contains one .csv file for each molecule, each containing all of the force field and DFT energetic data for all conformers from each force field. The columns of each .csv data file are: conf: The conformer numbers from the searches FFE: The energies of the conformers accourding to the force fields Structure: The filenames of the DFT minima E_SPC: The single-point corrected DFT energy E: The uncorrected M06-2X/6-31G(d) energy ZPE: The zero-point energy H_SPC: The single-point corrected enthalpy T.S: The product of the temperature (298.15 K) and the entropy T.qh-S: The product of the temperature and the quasi-harmonic entropy G(T)_SPC: The single-point corrected Gibbs free energy qh-G(T)_SPC: The single-point corrected quasi-harmonic Gibbs free energy