Appropriate clusterset selection for the prediction of thermodynamic properties of liquid water with QCE theory.

Evident in many physical and chemical phenomena, thermodynamics is the study of how energy is stored, transformed and transferred in a molecule or material. However, prediction of these properties with simulation techniques is a non-trivial task as several factors such as composition and intermolecular interactions come into play. While molecular dynamics and ab initio molecular dynamics are the most common techniques for the prediction of thermodynamic properties, there exists many shortcomings associated with their use. Therefore, in this work we instead apply QCE theory to predict the thermodynamic properties of liquid water. This theory assumes that a condensed phase system can be represented as a 'mixture' of varying sized clusters rather than as a continuum. As QCE theory relies on first-principle simulations and statistical thermodynamics to determine the thermodynamic behavior of a system, appropriate selection of clusters is a crucial step towards achieving accurate predictions. In this study, we use molecular dynamics and ab initio calculations to obtain initial configurations of 400 water clusters, Wn where n = 3 to 10 and contrast their stability using two different criteria. The role of entropy towards cluster stabilization is investigated by comparing the binding (ΔEBIND/mol) and Gibbs free binding energy per molecule (ΔGBIND/mol) of various Wn at 298.15 K. Initial clustersets are constructed by exploring two-, three-, four and five-combinations of clustersets using the minimum ΔGBIND/mol structures of Wn. We also expand the ΔGBIND/mol criteria for Wn of sizes 3 to 7 to include values larger than 0.0 kJ mol-1 and smaller than 3.0 kJ mol-1 as a means of improving thermodynamic predictions. 37 of the 459 resulting clustersets predicted the correct boiling point of water and its thermodynamic properties with an accuracy of 95%. A scaled population-weighted infrared spectrum was compared to experimental results to validate the composition of the top 5 clustersets. The predicted spectra showed an exact match within the low frequency range (<1000 cm-1) with some discrepancy at the high frequency range (>3400 cm-1). This work highlights that ΔGBIND/mol is so far the best criteria to apply when determining an appropriate clusterset for QCE theory.

Structural elucidation of polydopamine facilitated by ionic liquid solvation.

Minimal understanding of the formation mechanism and structure of polydopamine (pDA) and its natural analogue, eumelanin, impedes the practical application of these versatile polymers and limits our knowledge of the origin of melanoma. The lack of conclusive structural evidence stems from the insolubility of these materials, which has spawned significantly diverse suggestions of pDA's structure in the literature. We discovered that pDA is soluble in certain ionic liquids. Using these ionic liquids (ILs) as solvents, we present an experimental methodology to solvate pDA, enabling us to identify pDA's chemical structure. The resolved pDA structure consists of self-assembled supramolecular aggregates that contribute to the increasing complexity of the polymer. The underlying molecular energetics of pDA solvation and a macroscopic picture of the disruption of the aggregates using IL solvents have been investigated, along with studies of the aggregation mechanism in water.

Toward an extreme-scale electronic structure system.

Electronic structure calculations have the potential to predict key matter transformations for applications of strategic technological importance, from drug discovery to material science and catalysis. However, a predictive physicochemical characterization of these processes often requires accurate quantum chemical modeling of complex molecular systems with hundreds to thousands of atoms. Due to the computationally demanding nature of electronic structure calculations and the complexity of modern high-performance computing hardware, quantum chemistry software has historically failed to operate at such large molecular scales with accuracy and speed that are useful in practice. In this paper, novel algorithms and software are presented that enable extreme-scale quantum chemistry capabilities with particular emphasis on exascale calculations. This includes the development and application of the multi-Graphics Processing Unit (GPU) library LibCChem 2.0 as part of the General Atomic and Molecular Electronic Structure System package and of the standalone Extreme-scale Electronic Structure System (EXESS), designed from the ground up for scaling on thousands of GPUs to perform high-performance accurate quantum chemistry calculations at unprecedented speed and molecular scales. Among various results, we report that the EXESS implementation enables Hartree-Fock/cc-pVDZ plus RI-MP2/cc-pVDZ/cc-pVDZ-RIFIT calculations on an ionic liquid system with 623 016 electrons and 146 592 atoms in less than 45 min using 27 600 GPUs on the Summit supercomputer with a 94.6% parallel efficiency.

Behavior of counterpoise correction in many-body molecular clusters of organic compounds: Hartree-Fock interaction energy perspective.

The counterpoise (CP) correction by Boys and Bernardi has been well accepted as a reliable strategy to account for basis set superposition error (BSSE) in intermolecular complexes. The behavior of the CP correction was thoroughly studied in individual molecules of molecular complexes. This work studies the performance of the CP correction in many-body clusters including three-body clusters of organic compounds in the 3B-69 dataset. Additionally, we used crystal structures of polymorphs of benzene, aspirin, and oxalyl dihydrazide (ODH) to construct a many-body cluster dataset, abbreviated as the MBC-36 dataset, consisting of two, four and eight molecules, and 16 molecules in the case of benzene. A series of Dunning's basis sets-cc-pXZ and aug-cc-pXZ (X = D and T)-were used to predict CP-corrected Hartree-Fock (HF) interaction energies of the 3B-69 and MBC-36 datasets. The CP-corrected interaction energies were found to be basis-set independent, whereas the non-CP corrected interaction energies were found not to a follow a smooth exponential fitting as previously found for electronic energies of individual molecules. This observation was attributed to the presence of non-additive induction forces in some clusters. Two 2 × 2 × 2 supercells of benzene polymorphs were constructed to explore the local nature of BSSE effects. A cut-off radius of 10 Å was demonstrated to be sufficient to fully recover these effects. Although the behavior of CP correction was found to be non-conventional in many-body clusters of organic compounds, the use of a small basis set such as cc-pVDZ showed excellent performance in the prediction of HF interaction energies.

A DLPNO-CCSD(T) benchmarking study of intermolecular interactions of ionic liquids.

The accuracy of correlation energy recovered by coupled cluster single-, double-, and perturbative triple-excitations, CCSD(T), has led to the method being considered the gold standard of computational chemistry. The application of CCSD(T) has been limited to medium-sized molecular systems due to its steep scaling with molecular size. The recent development of alternative domain-based local pair natural orbital coupled-cluster method, DLPNO-CCSD(T), has significantly broadened the range of chemical systems to which CCSD(T) level calculations can be applied. Condensed systems such as ionic liquids (ILs) have a large contribution from London dispersion forces of up to 150 kJ mol-1 in large-scale clusters. Ionic liquids show appreciable charge transfer effects that result in the increased valence orbital delocalization over the entire ionic network, raising the question whether the application of methods based on localized orbitals is reliable for these semi-Coulombic materials. Here the performance of DLPNO-CCSD(T) is validated for the prediction of correlation interaction energies of two data sets incorporating single-ion pairs of protic and aprotic ILs. DLPNO-CCSD(T) produced results within chemical accuracy with tight parameter settings and a non-iterative treatment of triple excitations. To achieve spectroscopic accuracy of 1 kJ mol-1 , especially for hydrogen-bonded ILs and those containing halides, the DLPNO settings had to be increased by two orders of magnitude and include the iterative treatment of triple excitations, resulting in a 2.5-fold increase in computational cost. Two new sets of parameters are put forward to produce the performance of DLPNO-CCSD(T) within chemical and spectroscopic accuracy.

Explainable Solvation Free Energy Prediction Combining Graph Neural Networks with Chemical Intuition.

The prediction of a molecule's solvation Gibbs free (ΔGsolv) energy in a given solvent is an important task which has traditionally been carried out via quantum chemical continuum methods or force field-based molecular simulations. Machine learning (ML) and graph neural networks in particular have emerged as powerful techniques for elucidating structure-property relationships. This work presents a graph neural network (GNN) for the prediction of ΔGsolv which, in addition to encoding typical atom and bond-level features, incorporates chemically intuitive, solvation-relevant parameters into the featurization process: semiempirical partial atomic charges and solvent dielectric constant. Solute-solvent interactions are included via an interaction map layer which can be visualized to examine solubility-enhancing or -decreasing interactions learnt by the model. On a test set of small organic molecules, our GNN predicts ΔGsolv in water and cyclohexane with an accuracy comparable to polarizable and ab initio generated force field methods [mean absolute error (MAE) = 0.4 and 0.2 kcal mol-1, respectively], without the need for any molecular simulation. For the FreeSolv data set of hydration free energies, the test MAE is 0.7 kcal mol-1. Interpretability and applicability of the model is highlighted through several examples including rationalizing the increased solubility of modified diaminoanthraquinones in organic solvents. The clear explanations afforded by our GNN allow for easy understanding of the model's predictions, giving the experimental chemist confidence in employing ML models toward more optimized synthetic routes.

Experimental Evidence of Long-Lived Electric Fields of Ionic Liquid Bilayers.

Herein we demonstrate that ionic liquids can form long-lived double layers, generating electric fields detectable by straightforward open circuit potential (OCP) measurements. In imidazolium-based ionic liquids an external negative voltage pulse leads to an exceedingly stable near-surface dipolar layer, whose field manifests as long-lived (∼1-100 h) discrete plateaus in OCP versus time traces. These plateaus occur within an ionic liquid-specific and sharp potential window, defining a simple experimental method to probe the onset of interfacial ordering phenomena, such as overscreening and crowding. Molecular dynamics modeling reveals that the OCP arises from the alignment of the individual ion dipoles to the external electric field pulse, with the magnitude of the resulting OCP correlating with the product of the projected dipole moment of the cation and the ratio between the cation diffusion coefficient and its volume. Our findings also reveal that a stable overscreened structure is more likely to form if the interface is first forced through crowding, possibly accounting for the scattered literature data on relaxation kinetics of near-surface structures in ionic liquids.

Prediction of lattice energy of benzene crystals: A robust theoretical approach.

We present an inexpensive and robust theoretical approach based on the fragment molecular orbital methodology and the spin-ratio scaled second-order Møller-Plesset perturbation theory to predict the lattice energy of benzene crystals within 2 kJ⋅mol-1 . Inspired by the Harrison method to estimate the Madelung constant, the proposed approach calculates the lattice energy as a sum of two- and three-body interaction energies between a reference molecule and the surrounding molecules arranged in a sphere. The lattice energy converges rapidly at a radius of 13 Å. Adding the corrections to account for a higher correlated level of theory and basis set superposition for the Hartree Fock (HF) level produced a lattice energy of -57.5 kJ⋅mol-1 for the benzene crystal structure at 138 K. This estimate is within 1.6 kJ⋅mol-1 off the best theoretical prediction of -55.9 kJ⋅mol-1 . We applied this approach to calculate lattice energies of the crystal structures of phase I and phase II-polymorphs of benzene-observed at a higher temperature of 295 K. The stability of these polymorphs was correctly predicted, with phase II being energetically preferred by 3.7 kJ⋅mol-1 over phase I. The proposed approach gives a tremendous potential to predict stability of other molecular crystal polymorphs.

Establishing the accuracy of density functional approaches for the description of noncovalent interactions in ionic liquids.

We test a number of dispersion corrected versatile Generalized Gradient Approximation (GGA) and meta-GGA functionals for their ability to predict the interactions of ionic liquids, and show that most can achieve energies within 1 kcal mol-1 of benchmarks. This compares favorably with an accurate dispersion corrected hybrid, ωB97X-V. Our tests also reveal that PBE (Perdew-Burke-Ernzerhof GGA) calculations using the plane-wave projector augmented wave method and Gaussian Type Orbitals (GTOs) differ by less than 0.6 kJ mol-1 for ionic liquids, despite ions being difficult to evaluate in periodic cells - thus revealing that GTO benchmarks may be used also for plane-wave codes. Finally, the relatively high success of explicit van der Waals density functionals, compared to elemental and ionic dispersion models, suggests that improvements are required for low-cost dispersion correction models of ions.

Electrochemical characterization and thermodynamic analysis of TEMPO derivatives in ionic liquids.

In this study we investigate the reversibility of the reduction process of three TEMPO derivatives - TEMPOL, 4-cyano-TEMPO, and 4-oxo-TEMPO. The [C2mim][BF4] and [C4mpyr][OTf] ionic liquids (ILs) were used to perform cyclic voltammetry (CV) to analyse the redox potentials of the TEMPO derivatives. The former was previously shown to quench the aminoxy anion of TEMPO through a proton transfer reaction with the cation, whereas the latter supported the irreversibility of the TEMPO reduction process. In CV results on TEMPO derivatives, it was shown that [C4mpyr][OTf] could allow for a high degree of reversibility in the reduction of 4-cyano-TEMPO and a moderate degree of reversibility in the reduction of TEMPOL. In comparison, reduction of 4-cyano-TEMPO was predominantly irreversible in [C2mim][BF4], whilst TEMPOL showed complete irreversibility. 4-Oxo-TEMPO did not show any notable reduction reversibility in either IL tested. Reduction potentials showed little variation between the derivatives and 0.2 V variation between the ILs, with the most negative reduction potential being observed at -1.43 V vs. Fc/Fc+ for TEMPOL in [C4mpyr][OTf]. To explain the varying degrees of reversibility of the reduction process, four types of side reactions involving proton transfer to the aminoxy anion were studied using highly correlated quantum chemical methods. Proton transfer from the IL cation was shown to have the ability to quench all three aminoxy anions depending on the IL used. On average, TEMPOL was shown to be the most susceptible to proton transfer from the IL cation, having an average Gibbs free energy (GFE) of 10.5 kJ mol-1 more negative than that of 4-cyano-TEMPO, which was shown to have the highest GFE of proton transfer. Side reactions between water and aminoxy anions were also seen to have the potential to contribute to degradation of the aminoxy anions tested, with 4-oxo-TEMPO being shown to be the most reactive to degradation with water with a GFE of -12.6 kJ mol-1. 4-Oxo-TEMPO was found to be highly susceptible to self-quenching by its aminoxy anion and radical form with highly negative proton transfer GFEs of -47.9 kJ mol-1 and -57.7 kJ mol-1, respectively. Overall, 4-cyano-TEMPO is recommended as being the most stable of the aminoxy anions tested with TEMPOL, thus providing a viable alternative to improve solubility should the IL be tuned to maximize its stability.

Ordered Solvents and Ionic Liquids Can Be Harnessed for Electrostatic Catalysis.

Herein, we employ classical molecular dynamics simulations using the Drude oscillator-based polarizable force field, quantum chemical calculations, and ONIOM multiscale calculations to study (a) how an external field orders the solvent environment in a chemical reaction and then (b) whether in the absence of this same applied field the ordered solvent environment alone can electrostatically catalyze a chemical reaction when compared with the corresponding disordered solvent. Our results show that a 0.2 V/Å external electric field, which is below the threshold for bond breaking of solvent molecules, leads to significant ordering of bulk methanol solvent and the ionic liquid [EMIM][BF4]. Importantly, in the absence of this same field, the ordered solvent lowers the activation energy of the hydrogen-transfer reaction of o-alkylphenyl ketones in excess of 20 kcal/mol when the solvent is methanol and by over 30 kcal/mol for [EMIM][BF4]. Even a 0.1 V/Å external field has effects of ca. 10 and 20 kcal/mol, respectively. This work suggests a possible strategy for scaling electrostatic catalysis by applying a pulsed external field to the reaction medium to maintain solvent ordering while allowing the reaction to proceed largely in the absence of an external field.

The effect of descriptor choice in machine learning models for ionic liquid melting point prediction.

The characterization of an ionic liquid's properties based on structural information is a longstanding goal of computational chemistry, which has received much focus from ab initio and molecular dynamics calculations. This work examines kernel ridge regression models built from an experimental dataset of 2212 ionic liquid melting points consisting of diverse ion types. Structural descriptors, which have been shown to predict quantum mechanical properties of small neutral molecules within chemical accuracy, benefit from the addition of first-principles data related to the target property (molecular orbital energy, charge density profile, and interaction energy based on the geometry of a single ion pair) when predicting the melting point of ionic liquids. Out of the two chosen structural descriptors, ECFP4 circular fingerprints and the Coulomb matrix, the addition of molecular orbital energies and all quantum mechanical data to each descriptor, respectively, increases the accuracy of surrogate models for melting point prediction compared to using the structural descriptors alone. The best model, based on ECFP4 and molecular orbital energies, predicts ionic liquid melting points with an average mean absolute error of 29 K and, unlike group contribution methods, which have achieved similar results, is applicable to any type of ionic liquid.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids.

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Increased stability of nitroxide radicals in ionic liquids: more than a viscosity effect.

Radical stability has been subject to continuous research due to its importance in polymerization as well as in all-organic batteries. Recently, the SOMO-HOMO conversion was identified as the main factor in controlling the stability of distonic radicals, for which the negative charge resides on the same molecule. Based on this finding, the idea of ionic liquids stabilizing radicals was hypothesized in this study. A series of ionic liquids were tested in EPR measurements of the 3-carboxy-2,2,5,5-tetramethyl-pyrroline-1-oxyl. Unusually high rotational diffusion constants (τR), 4 times larger compared to conventional media such as dichloromethane (DCM), were recorded at room temperature. This finding could only be explained by a strong interaction existing between the radical and ionic liquid ions, which was confirmed with quantum chemical calculations, with interaction energies falling between -17.1 kJ mol-1 for tetramethylphosphonium tetrafluoroborate and -85.6 kJ mol-1 for 1,3-dimethylimidazolium triflate. Elevated temperature measurements performed at 80 °C reduced the viscosity of the ionic liquids to that of DCM, while the τR values remained relatively high, thus further confirming that the rotational hindrance occurred due to radical-ionic liquid interactions. The calculated interaction energies between the radical and ionic liquids ions were also found to correlate well with experimental rotational diffusion constants, thus offering us a valuable tool in tailoring ionic liquids for enhanced stability of nitroxide radicals. The findings of this study showcase the ability of ionic liquids to reduce reactivity of nitroxides without the need for any chemical modification of the radical.

Contrasting Synergistic Heterobimetallic (Na-Mg) and Homometallic (Na or Mg) Bases in Metallation Reactions of Dialkylphenylphosphines and Dialkylanilines: Lateral versus Ring Selectivities.

A series of dialkylphenylphosphines and their analogous aniline substrates have been metallated with the synergistic mixed-metal base [(TMEDA)Na(TMP)(CH2 SiMe3 )Mg(TMP)] 1. Different metallation regioselectivities for the substrates were observed, with predominately lateral or meta-magnesiated products isolated from solution. Three novel heterobimetallic complexes [(TMEDA)Na(TMP)(CH2 PCH3 Ph)Mg(TMP)] 2, [(TMEDA)Na(TMP)(m-C6 H4 PiPr2 )Mg(TMP)] 3 and [(TMEDA)Na(TMP)(m-C6 H4 NEt2 )Mg(TMP)]  4 and two homometallic complexes [{(TMEDA)Na(EtNC6 H5 )}2 ] 5 and [(TMEDA)Na2 (TMP)(C6 H5 PEt)]2  6 derived from homometallic metallation have been crystallographically characterised. Complex 6 is an unprecedented sodium-amide, sodium-phosphide hybrid with a rare (NaNNaP)2 ladder motif. These products reveal contrasting heterobimetallic deprotonation with homometallic induced ethene elimination reactivity. Solution studies of metallation mixtures and electrophilic iodine quenching reactions confirmed the metallation sites. In an attempt to rationalise the regioselectivity of the magnesiation reactions the C-H acidities of the six substrates were determined in THF solution using DFT calculations employing the M06-2X functional and cc-pVTZ Dunning's basis set.

Anisotropic Thermal and Guest-Induced Responses of an Ultramicroporous Framework with Rigid Linkers.

The interdependent effects of temperature and guest uptake on the structure of the ultramicroporous metal-organic framework [Cu3 (cdm)4 ] (cdm=C(CN)2 (CONH2 )- ) were explored in detail by using in situ neutron scattering and density functional theory calculations. The tetragonal lattice displays an anisotropic thermal response related to a hinged "lattice-fence" mechanism, unusual for this topology, which is facilitated by pivoting of the rigid cdm anion about the Cu nodes. Calculated pore-size metrics clearly illustrate the potential for temperature-mediated adsorption in ultramicroporous frameworks due to thermal fluctuations of the pore diameter near the value of the target guest kinetic diameter, though in [Cu3 (cdm)4 ] this is counteracted by a competing contraction of the pore with increasing temperature as a result of the anisotropic lattice response.

Spectroscopic Studies on Photoinduced Reactions of the Anticancer Prodrug, trans,trans,trans-[Pt(N3 )2 (OH)2 (py)2 ].

The photodecomposition mechanism of trans,trans,trans-[Pt(N3 )2 (OH)2 (py)2 ] (1, py=pyridine), an anticancer prodrug candidate, was probed using complementary Attenuated Total Reflection Fourier Transform Infrared (ATR-FTIR), transient electronic absorption, and UV/Vis spectroscopy. Data fitting using Principal Component Analysis (PCA) and Multi-Curve Resolution Alternating Least Squares, suggests the formation of a trans-[Pt(N3 )(py)2 (OH/H2 O)] intermediate and trans-[Pt(py)2 (OH/H2 O)2 ] as the final product upon 420 nm irradiation of 1 in water. Rapid disappearance of the hydroxido ligand stretching vibration upon irradiation is correlated with a -10 cm-1 shift to the antisymmetric azido vibration, suggesting a possible second intermediate. Experimental proof of subsequent dissociation of azido ligands from platinum is presented, in which at least one hydroxyl radical is formed in the reduction of PtIV to PtII . Additionally, the photoinduced reaction of 1 with the nucleotide 5'-guanosine monophosphate (5'-GMP) was comprehensively studied, and the identity of key photoproducts was assigned with the help of ATR-FTIR spectroscopy, mass spectrometry, and density functional theory calculations. The identification of marker bands for some of these photoproducts (e.g., trans-[Pt(N3 )(py)2 (5'-GMP)] and trans-[Pt(py)2 (5'-GMP)2 ]) will aid elucidation of the chemical and biological mechanism of anticancer action of 1. In general, these studies demonstrate the potential of vibrational spectroscopic techniques as promising tools for studying such metal complexes.

Bulk properties of aqueous graphene oxide and reduced graphene oxide with surfactants and polymers: adsorption and stability.

A diverse range of molecular surfactants and polymers have been incorporated into aqueous graphene oxide (GO) and reduced graphene oxide (rGO) dispersions in order to understand the complex relationship between surface chemistry, surface forces and interfacial thermodynamics of these materials with typical amphiphiles. Surfactant additives were systematically varied in terms of their charge and hydrophobicity to reveal important structure-function relationships affecting adsorption and interaction with GO and rGO surfaces. Small-angle (and ultra small-angle) neutron scattering was employed to examine and monitor the interactions and self-assembly in each system. Charge was found to be the overriding factor driving adsorption, as cationic surfactants very readily adsorbed to both GO and rGO, whereas anionic surfactants gave little to no evidence of adsorption despite possessing hydrophobic tail-groups. Molecules of neutral charge such as nonionic and zwitterionic surfactants as well as neutral polymers also showed strong affinities for GO and rGO, indicating that dispersion and dipole (induction polarisation) interactions also play a significant role in adsorption with these materials. Modelling the neutron data revealed in many cases a q-2 slope in the low q and ultra low q regions, indicating that scattering was occurring from large, flat surfaces (lamellae or bilayers), suggesting an effective flattening of the sheets in dispersion. The results presented thus help to form a roadmap for the behaviour of GO and rGO with surfactants and polymers, relevant to adsorption, stabilisation, formulation and coating in aqueous environments as adsorbent and functional materials.

Cluster approach to the prediction of thermodynamic and transport properties of ionic liquids.

The prediction of physicochemical properties of ionic liquids such as conductivity and melting point would substantially aid the targeted design of ionic liquids for specific applications ranging from solvents for extraction of valuable chemicals to biowaste to electrolytes in alternative energy devices. The previously published study connecting the interaction energies of single ion pairs (1 IP) of ionic liquids to their thermodynamic and transport properties has been extended to larger systems consisting of two ion pairs (2 IPs), in which many-body and same-ion interactions are included. Routinely used cations, of the imidazolium and pyrrolidinium families, were selected in the study coupled with chloride, tetrafluoroborate, and dicyanamide. Their two ion pair clusters were subjected to extensive configuration screening to establish most stable structures. Interaction energies of these clusters were calculated at the spin-ratio scaled MP2 (SRS-MP2) level for the correlation interaction energy, and a newly developed scaled Hartree-Fock method for the rest of energetic contributions to interaction energy. A full geometry screening for each cation-anion combination resulted in 192 unique structures, whose stability was assessed using two criteria-widely used interaction energy and total electronic energy. Furthermore, the ratio of interaction energy to its dispersion component was correlated with experimentally observed melting points in 64 energetically favourable structures. These systems were also used to test the correlation of the dispersion contribution to interaction energy with measured conductivity.

Prediction of 1H NMR chemical shifts for clusters of imidazolium-based ionic liquids.

Nuclear magnetic resonance (NMR) has been widely used to elucidate the bulk structure of ionic liquids. In this work, we calculated 1H NMR chemical shifts of 1-ethyl-3-methylimidazolium (C2mim+) ionic liquids combined with various anions such as chloride (Cl), tetrafluoroborate (BF4), hexafluorophosphate (PF6), acetate (OAc), trifluoroacetate (TFA), and dicyanamide (DCA). The previously established level of theory, HF/6-311G+(3df,2p), was used for the accurate prediction of NMR chemical shifts both in gas phase and in solvents with varying dielectric constant such as CHCl3 and ethanol. The following factors affecting the predicted proton chemical shifts were considered. Firstly, ionic clusters consisting of 2, 8 and 16 ion pairs were optimized to model interionic interactions present in the bulk of ionic liquids. In larger clusters the distribution of the calculated chemical shifts of individual protons in the C2mim+ cation was examined with respect to the position of the cation in the cluster. We further confirmed that electronic properties of ionic liquids such as magnetic shielding had local nature, thus allowing us to accurately predict proton NMR chemical shifts of ionic liquids from relatively small-sized clusters. Secondly, solvent effects in single ion pairs as well as larger ionic clusters were accounted through a Conductor-like Polarisable Continuum Model (CPCM). Solvent effects generated through a dielectric constant of either chloroform or ethanol were found to be important in single ion pairs due to improved description of interionic distances. With increasing cluster size the difference between gas-phase and CPCM optimized structures became minimal, thus resulting in similar values for calculated 1H NMR chemical shifts. We also established that the model size that produced the best results for imidazolium ionic liquids strongly depended on the anion type. Strongly coordinating anions such as chloride and acetate require calculations of clusters consisting of at least 8 ion pairs, whereas weakly coordinating anions produce excellent accuracy for single ion pairs optimized in the presence of solvent. The polarity of the solvent was found to play a minor role.