Finding the best density functional approximation to describe interaction energies and structures of ionic liquids in molecular dynamics studies

Ionic liquids raise interesting but complicated questions for theoretical investigations due to the fact that a number of different inter-molecular interactions, e.g., hydrogen bonding, long-range Coulomb interactions, and dispersion interactions, need to be described properly. Here, we present a detailed study on the ionic liquids ethylammonium nitrate and 1-ethyl-3-methylimidazolium acetate, in which we compare different dispersion corrected density functional approximations to accurate local coupled cluster data in static calculations on ionic liquid clusters. The efficient new composite method B97-3c is tested and has been implemented in CP2K for future studies. Furthermore, tight-binding based approaches which may be used in large scale simulations are assessed. Subsequently, ab initio as well as classical molecular dynamics simulations are conducted and structural analyses are presented in order to shed light on the different short- and long-range structural patterns depending on the method and the system size considered in the simulation. Our results indicate the presence of strong hydrogen bonds in ionic liquids as well as the aggregation of alkyl side chains due to dispersion interactions.

Lewis Acidity and Basicity of Mixed Chlorometallate Ionic Liquids: Investigations from Surface Analysis and Fukui Function

Mixed chlorometallate ionic liquids (ILs) have been regarded as potential solvents, catalysts, and reagents for many organic processes. The acidity and basicity of these ILs were correlated with theoretically estimated parameters such as electrostatic surface potential maxima and minima, average local surface ionization energy, and Fukui and dual descriptor functions. The introduction of metal chloride into the anions would influence the acidity/basicity of ILs by withdrawing the electron density from the cationic counterpart. For the [C₄mim]-based ILs with the mixed-metal anions, the acidity tends to attenuate while the basicity becomes stronger, as compared to the corresponding chloroaluminate ILs. However, the acidity of [(C₂H₅)₃NH]-based ILs with the mixed-metal anions are greater than that of the net chloroaluminate ILs. The Fukui function values showed that most of the mixed chlorometallate ILs belong to bifunctional distribution. The mixed chlorometallate ILs both have electrophilic and nucleophilic sites, which would be beneficial for their applications.

Prediction of 1H NMR chemical shifts for ionic liquids: strategy and application of a relative reference standard

The computational predictions of ¹H NMR chemical shifts for ionic liquids were investigated. To calculate the chemical shifts more accurately, the approach of relative reference standard (RRS) was proposed. This straightforward computational technique uses organic compounds and ionic liquids that are similar to the studied ionic liquids as standards. The calculated chemical shifts of single ion pairs were strongly influenced by the anion type and the local environment. Using the RRS methodology, the calculated results agreed well with the experimental chemical shifts due to the cancellation of errors caused by the anion. Ionic clusters consisting of 4 ion pairs were also employed to model the ionic liquids with strongly coordinating anions. Large clusters slightly improve the accuracy by reducing systematic errors. Although the experimental ¹H NMR data of the reference ionic liquid should be used, the RRS methodology has been shown to predict ¹H NMR chemical shifts efficiently at different levels of theory.

Using an RRS method to calculate the ¹H NMR chemical shifts of ionic liquid agreed well with the experimental value.

Accurate prediction of the structure and vibrational spectra of ionic liquid clusters with the generalized energy-based fragmentation approach: critical role of ion-pair-based fragmentation

The GEBF method with the ion-pair-based fragmentation has been developed to facilitate ab initio calculations of general ionic liquid clusters.

Mapping the Free Energy of Lithium Solvation in the Protic Ionic Liquid Ethylammonuim Nitrate: A Metadynamics Study

Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 eV. Other less probable conformations having different solvation pattern are also investigated.

Quantum Chemical Modeling of Hydrogen Bonding in Ionic Liquids

Hydrogen bonding (H-bonding) is an important and very general phenomenon. H-bonding is part of the basis of life in DNA, key in controlling the properties of water and ice, and critical to modern applications such as crystal engineering, catalysis applications, pharmaceutical and agrochemical development. H-bonding also plays a significant role for many ionic liquids (IL), determining the secondary structuring and affecting key physical parameters. ILs exhibit a particularly diverse and wide range of traditional as well as non-standard forms of H-bonding, in particular the doubly ionic H-bond is important. Understanding the fundamental nature of the H-bonds that form within ILs is critical, and one way of accessing this information, that cannot be recovered by any other computational method, is through quantum chemical electronic structure calculations. However, an appropriate method and basis set must be employed, and a robust procedure for determining key structures is essential. Modern generalised solvation models have recently been extended to ILs, bringing both advantages and disadvantages. QC can provide a range of information on geometry, IR and Raman spectra, NMR spectra and at a more fundamental level through analysis of the electronic structure.

The role of errors related to DFT methods in calculations involving ion pairs of ionic liquids

Ionic liquids (ILs) play a key role in many chemical applications. As regards the theoretical approach, ILs show added difficulties in calculations due to the composition of the ion pair and to the fact that they are liquids. Although density functional theory (DFT) can treat this kind of systems to predict physico-chemical properties, common versions of these methods fail to perform accurate predictions of geometries, interaction energies, dipole moments, and other properties related to the molecular structure. In these cases, dispersion and self-interaction error (SIE) corrections need to be introduced to improve DFT calculations involving ILs. We show that the inclusion of dispersion is needed to obtain good geometries and accurate interaction energies. SIE needs to be corrected to describe the charges and dipoles in the ion pair correctly. The use of range-separated functionals allows us to obtain interaction energies close to the CCSD(T) level. © 2017 Wiley Periodicals, Inc.

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.

Controlling Electronic Transitions in Fullerene van der Waals Aggregates via Supramolecular Assembly

Morphologies crucially determine the optoelectronic properties of organic semiconductors. Therefore, hierarchical and supramolecular approaches have been developed for targeted design of supramolecular ensembles of organic semiconducting molecules and performance improvement of, e.g., organic solar cells (OSCs), organic light emitting diodes (OLEDs), and organic field-effect transistors (OFETs). We demonstrate how the photonic properties of fullerenes change with the formation of van der Waals aggregates. We identified supramolecular structures with broadly tunable absorption in the visible spectral range and demonstrated how to form aggregates with targeted visible (vis) absorption. To control supramolecular structure formation, we functionalized the C60-backbone with polar (bis-polyethylene glycol malonate-MPEG) tails, thus yielding an amphiphilic fullerene derivative that self-assembles at interfaces. Aggregates of systematically tuned size were obtained from concentrating MPEGC60 in stearic acid matrices, while different supramolecular geometries were provoked via different thin film preparation methods, namely spin-casting and Langmuir-Blodgett (LB) deposition from an air-water interface. We demonstrated that differences in molecular orientation in LB films (C2v type point group aggregates) and spin-casting (stochastic aggregates) lead to huge changes in electronic absorption spectra due to symmetry and orientation reasons. These differences in the supramolecular structures, causing the different photonic properties of spin-cast and LB films, could be identified by means of quantum chemical calculations. Employing supramolecular assembly, we propounded that molecular symmetry in fullerene aggregates is extremely important in controlling vis absorption to harvest photons efficiently, when mixed with a donor molecule, thus improving active layer design and performance of OSCs.

Vibrational analysis and formation mechanism of typical deep eutectic solvents: An experimental and theoretical study

Deep eutectic solvents (DESs), as ionic liquid analogues for green solvents, have gained increasing attentions in chemistry. In this work, three typical kinds of DESs (ChCl/Gly, ChCl/AcOH and ChCl/Urea) were successfully synthesized and characterized by Fourier transform infrared spectroscopy (FTIR) and Raman. Then comprehensive and systematical analyses were performed by the methods of density functional theory (DFT). Two methods (B3LYP/6-311++G(2d,p) and dispersion-corrected B3LYP-D3/6-311++G(2d,p)) were employed to investigate the structures, vibrational frequencies and assign their ownership of vibrational modes for the DESs, respectively. Nearly all the experimental characteristic peaks of IR and Raman were identified according to the calculated results. By linear fitting of the combined calculated vs experimental vibration frequencies, it can be found that both of the two methods are excellent to reproduce the experimental results. Besides, hydrogen bonds were proved to exist in DESs by IR spectrum, structure analysis and RDG analysis. This work was aimed at predicting and understanding the vibrational spectra of the three typical DESs based on DFT methods. Moreover, by comparing experimental and theoretical results, it provides us a deep understanding of the formation mechanisms of DESs.

Structures and Electronic Properties of Lithium Chelate-Based Ionic Liquids

The conformations, electronic properties, and interaction energies of four chelate-based ionic liquids [Li(EA)][Tf2N], [Li(HDA)][Tf2N], [Li(DEA)][Tf2N], and [Li(DOBA)][Tf2N] have been theoretically explored. The reliability of the located conformers has been confirmed via the comparison between the simulated and experimental infrared spectra. Our results show that the N-Li and O-Li coordinate bonds in cation are elongated as the numbers of coordinate heteroatoms of alkanolamine ligands to Li(+) increased. Also the binding energies between Li(+) and ligands are increased and the interaction energies between cations and Tf2N anion are decreased. The cation-anion interaction energies follow the order of [Li(DOBA)][Tf2N] < [Li(HDA)][Tf2N] < [Li(DEA)][Tf2N] < [Li(EA)][Tf2N], which fall within the energetic ranges of conventional ionic liquids. Interestingly, the strongest stabilization orbital interactions in these ionic liquids and their cations revealed by the natural bond orbital analysis lie in the interaction between the lone pair (LP) of the coordinate heteroatoms in ligands or anion as donors and the vacant valence shell nonbonding orbital (LP*) of Li(+) as acceptors, which are very different from that of conventional ionic liquids. Moreover, the charges transferred from cations to anion are quite similar, and the charge of Li(+) is proposed for possibly predicting the order of the interaction energies of ionic liquids in series. The present study allows for the deeper understanding the differences between chelate-based ionic liquids and conventional ionic liquids.

Self-interaction error in DFT-based modelling of ionic liquids

The modern computer simulations of potential green solvents of the future, involving the room temperature ionic liquids, heavily rely on density functional theory (DFT). In order to verify the appropriateness of the common DFT methods, we have investigated the effect of the self-interaction error (SIE) on the results of DFT calculations for 24 ionic pairs and 48 ionic associates. The magnitude of the SIE is up to 40 kJ mol(-1) depending on the anion choice. Most strongly the SIE influences the calculation results of ionic associates that contain halide anions. For these associates, the range-separated density functionals suppress the SIE; for other cases, the revPBE density functional with dispersion correction and triple-ζ Slater-type basis is suitable for computationally inexpensive and reasonably accurate DFT calculations.

Interaction between ionic liquid cation and water: infrared predissociation study of [bmim]+·(H2O)n clusters

The infrared predissociation spectra of [bmim]⁺·(H₂O)_n, n = 1–8, in the 2800–3800 cm⁻¹ region are presented and analyzed with the help of electronic structure calculations.

Range-Separated meta-GGA Functional Designed for Noncovalent Interactions

The accuracy of applying density functional theory to noncovalent interactions is hindered by errors arising from low-density regions of interaction-induced change in the density gradient, error compensation between correlation and exchange functionals, and dispersion double counting. A new exchange-correlation functional designed for noncovalent interactions is proposed to address these problems. The functional consists of the range-separated PBEsol exchange considered in two variants, pure and hybrid, and the semilocal correlation functional of Modrzejewski et al. (J. Chem. Phys. 2012, 137, 204121) designed with the constraint satisfaction technique to smoothly connect with a dispersion term. Two variants of dispersion correction are appended to the correlation functional: the atom-atom pairwise additive DFT-D3 model and the density-dependent many-body dispersion with self-consistent screening (MBD-rsSCS). From these building blocks, a set of four functionals is created to systematically examine the role of pure versus hybrid exchange and the underlying models for dispersion. The new functional is extensively tested on benchmark sets with diverse nature and size. Truly outstanding performance is demonstrated for water clusters of varying size, ionic hydrogen bonds, and thermochemistry of isodesmic n-alkane fragmentation reactions. The merits of each component of the new functional are discussed.

Consistent structures and interactions by density functional theory with small atomic orbital basis sets

A density functional theory (DFT) based composite electronic structure approach is proposed to efficiently compute structures and interaction energies in large chemical systems. It is based on the well-known and numerically robust Perdew-Burke-Ernzerhoff (PBE) generalized-gradient-approximation in a modified global hybrid functional with a relatively large amount of non-local Fock-exchange. The orbitals are expanded in Ahlrichs-type valence-double zeta atomic orbital (AO) Gaussian basis sets, which are available for many elements. In order to correct for the basis set superposition error (BSSE) and to account for the important long-range London dispersion effects, our well-established atom-pairwise potentials are used. In the design of the new method, particular attention has been paid to an accurate description of structural parameters in various covalent and non-covalent bonding situations as well as in periodic systems. Together with the recently proposed three-fold corrected (3c) Hartree-Fock method, the new composite scheme (termed PBEh-3c) represents the next member in a hierarchy of "low-cost" electronic structure approaches. They are mainly free of BSSE and account for most interactions in a physically sound and asymptotically correct manner. PBEh-3c yields good results for thermochemical properties in the huge GMTKN30 energy database. Furthermore, the method shows excellent performance for non-covalent interaction energies in small and large complexes. For evaluating its performance on equilibrium structures, a new compilation of standard test sets is suggested. These consist of small (light) molecules, partially flexible, medium-sized organic molecules, molecules comprising heavy main group elements, larger systems with long bonds, 3d-transition metal systems, non-covalently bound complexes (S22 and S66×8 sets), and peptide conformations. For these sets, overall deviations from accurate reference data are smaller than for various other tested DFT methods and reach that of triple-zeta AO basis set second-order perturbation theory (MP2/TZ) level at a tiny fraction of computational effort. Periodic calculations conducted for molecular crystals to test structures (including cell volumes) and sublimation enthalpies indicate very good accuracy competitive to computationally more involved plane-wave based calculations. PBEh-3c can be applied routinely to several hundreds of atoms on a single processor and it is suggested as a robust "high-speed" computational tool in theoretical chemistry and physics.

Ion pairing in ionic liquids

In the present article we briefly review the extensive discussion in literature about the presence or absence of ion pair-like aggregates in ionic liquids. While some experimental studies point towards the presence of neutral subunits in ionic liquids, many other experiments cannot confirm or even contradict their existence. Ion pairs can be detected directly in the gas phase, but no direct method is available to observe such association behavior in the liquid, and the corresponding indirect experimental proofs are based on such assumptions as unity charges at the ions. However, we have shown by calculating ionic liquid clusters of different sizes that assuming unity charges for ILs is erroneous, because a substantial charge transfer is taking place between the ionic liquid ions that reduce their total charge. Considering these effects might establish a bridge between the contradicting experimental results on this matter. Beside these results, according to molecular dynamics simulations the lifetimes of ion-ion contacts and their joint motions are far too short to verify the existence of neutral units in these materials.

Structure, dynamics and stability of water/scCO2/mineral interfaces from ab initio molecular dynamics simulations

The boundary layer at solid-liquid interfaces is a unique reaction environment that poses significant scientific challenges to characterize and understand by experimentation alone. Using ab initio molecular dynamics (AIMD) methods, we report on the structure and dynamics of boundary layer formation, cation mobilization and carbonation under geologic carbon sequestration scenarios (T = 323 K and P = 90 bar) on a prototypical anorthite (001) surface. At low coverage, water film formation is enthalpically favored, but entropically hindered. Simulated adsorption isotherms show that a water monolayer will form even at the low water concentrations of water-saturated scCO2. Carbonation reactions readily occur at electron-rich terminal Oxygen sites adjacent to cation vacancies that readily form in the presence of a water monolayer. These results point to a carbonation mechanism that does not require prior carbonic acid formation in the bulk liquid. This work also highlights the modern capabilities of theoretical methods to address structure and reactivity at interfaces of high chemical complexity.

Multiple-Site SO2 -Capture Ionic Liquids with Nearly Uniform Site Performance

We propose a series of azolium poly(azolyl)borate ionic liquids (ILs) for reversible SO₂ capture. Density functional calculations demonstrate that the designed borate anions can strongly bond to SO₂ at multiple sites with nearly uniform binding energies. Thus, as well as high overall uptakes, the ILs can achieve much higher effective uptakes (the uptake difference between absorption and desorption conditions) than existing SO₂ -capture reagents. The larger size of the borate anions, the evenly distributed negative charge among the azolyl rings, and the blocking of the conjugation by the tetrahedral boron concertedly reduce absorbate-absorbate repulsion, which leads to a large disparity among binding sites in other multiple-site SO₂ sorbents.

Estimation of Lattice Enthalpies of Ionic Liquids Supported by Hirshfeld Analysis

New measurements of vaporization enthalpies for 15 1:1 ionic liquids are performed by using a quartz-crystal microbalance. Collection and analysis of 33 available crystal structures of organic salts, which comprise 13 different cations and 12 anions, is performed. Their dissociation lattice enthalpies are calculated by a combination of experimental and quantum chemical quantities and are divided into the relaxation and Coulomb components to give an insight into elusive short-range interaction enthalpies. An empirical equation is developed, based on interaction-specific Hirshfeld surfaces and solvation enthalpies, which enables the estimation of the lattice enthalpy by using only the crystal-structure data.

Carboxylate-assisted C-H activation of phenylpyridines with copper, palladium and ruthenium: a mass spectrometry and DFT study

The transition state of metal carboxylate mediated C–H activation is associated with carbon–metal bond formation supported by electron-poor carboxylates.

The C–H activation of 2-phenylpyridine, catalyzed by copper(ii), palladium(ii) and ruthenium(ii) carboxylates, was studied in the gas phase. ESI-MS, infrared multiphoton dissociation spectroscopy and quantum chemical calculations were combined to investigate the intermediate species in the reaction. Collision induced dissociation (CID) experiments and DFT calculations allowed estimation of the energy required for this C–H activation step and the subsequent acetic acid loss. Hammett plots constructed from the CID experiments using different copper carboxylates as catalysts revealed that the use of stronger acids accelerates the C–H activation step. The reasoning can be traced from the associated transition structures that suggest a concerted mechanism and the key effect of the carbon–metal bond pre-formation. Carboxylates derived from stronger acids make the metal atom more electrophilic and therefore shift the reaction towards the formation of C–H activated products.

Hydrogen bonding in ionic liquids

Ionic liquids (IL) and hydrogen bonding (H-bonding) are two diverse fields for which there is a developing recognition of significant overlap. Doubly ionic H-bonds occur when a H-bond forms between a cation and anion, and are a key feature of ILs. Doubly ionic H-bonds represent a wide area of H-bonding which has yet to be fully recognised, characterised or explored. H-bonds in ILs (both protic and aprotic) are bifurcated and chelating, and unlike many molecular liquids a significant variety of distinct H-bonds are formed between different types and numbers of donor and acceptor sites within a given IL. Traditional more neutral H-bonds can also be formed in functionalised ILs, adding a further level of complexity. Ab initio computed parameters; association energies, partial charges, density descriptors as encompassed by the QTAIM methodology (ρBCP), qualitative molecular orbital theory and NBO analysis provide established and robust mechanisms for understanding and interpreting traditional neutral and ionic H-bonds. In this review the applicability and extension of these parameters to describe and quantify the doubly ionic H-bond has been explored. Estimating the H-bonding energy is difficult because at a fundamental level the H-bond and ionic interaction are coupled. The NBO and QTAIM methodologies, unlike the total energy, are local descriptors and therefore can be used to directly compare neutral, ionic and doubly ionic H-bonds. The charged nature of the ions influences the ionic characteristics of the H-bond and vice versa, in addition the close association of the ions leads to enhanced orbital overlap and covalent contributions. The charge on the ions raises the energy of the Ylp and lowers the energy of the X-H σ* NBOs resulting in greater charge transfer, strengthening the H-bond. Using this range of parameters and comparing doubly ionic H-bonds to more traditional neutral and ionic H-bonds it is clear that doubly ionic H-bonds cover the full range of weak through to very strong H-bonds.

Simulating the vibrational spectra of ionic liquid systems: 1-ethyl-3-methylimidazolium acetate and its mixtures

The vibrational spectra of the ionic liquid 1-ethyl-3-methylimidazolium acetate and its mixtures with water and carbon dioxide are calculated using ab initio molecular dynamics simulations, and the results are compared to experimental data. The new implementation of a normal coordinate analysis in the trajectory analyzer TRAVIS is used to assign the experimentally observed bands to specific molecular vibrations. The applied computational approaches prove to be particularly suitable for the modeling of bulk phase effects on vibrational spectra, which are highly important for the discussion of the microscopic structure in systems with a strong dynamic network of intermolecular interactions, such as ionic liquids.

On the origin of ionicity in ionic liquids. Ion pairing versus charge transfer

In this paper we show by using static DFT calculations and classical molecular dynamics simulations that the charge transfer between ionic liquid ions plays a major role in the observed discrepancies between the overall mobility of the ions and the observed conductivities of the corresponding ionic liquids, while it also directly suppresses the association of oppositely charged ions, thus the ion pairing. Accordingly, in electrochemical applications of these materials it is important to consider this reduction of the total charges on the ions, which can greatly affect the performance of the given process or device in which the ionic liquid is used. By slightly shifting from the salt-like to a molecular liquid-like system via the decreased charges, the charge transfer also fluidizes the ionic liquid. We believe that this vital information on the molecular level structure of ionic liquids offers a better understanding of these materials, and allows us to improve the a priori design of ionic liquids for any given purpose.

Theoretical investigation of the interaction between aromatic sulfur compounds and [BMIM](+)[FeCl4](-) ionic liquid in desulfurization: A novel charge transfer mechanism

In this work, interaction nature between a group of aromatic sulfur compounds and [BMIM](+)[FeCl4](-) have been investigated by density functional theory (DFT). A coordination structure is found to be critical to the mechanism of extractive desulfurization. Interaction energy and extractive selectivity follow the order: thiophene (TH)<dibenzothiophene (DBT)≈benzothiophene (BT). Alkylation of TH or BT (e.g. 3-methylthiophene, and 3-methylbenzothiophene) leads to a stronger interaction with ionic liquid, but steric hindrance effects of some alkylic derivatives (e.g. 2,7-dimethylbenzothiophene) lead to a weaker interaction with ionic liquid. The mechanism of extractive desulfurization is attributed to the charge transfer effect. During extractive desulfurization, electrons on aromatic sulfur compounds transfer into the Lewis part of ionic liquid, namely, [FeCl4](-). Furthermore, it is better to consider the Lewis acidity of Fe-containing ionic liquid by the whole unit (such as [FeCl4](-) and aromatic sulfur compounds (X)) rather than only Fe or S atom.

Molecular dynamics study of the effect of alkyl chain length on melting points of [CnMIM][PF6] ionic liquids

Based on molecular dynamics simulations, the melting points Tm of a series of 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids [CnMIM][PF6] with n = 2, 4, 10, 12, and 14 were studied using the free energy-based pseudosupercritical path (PSCP) method. The experimental trend that the Tm decreases with increasing alkyl chain length for ILs with short alkyl chains and increases for the ones with long alkyl chains was correctly captured. Further analysis revealed that the different trends are the results of the balance between fusion enthalpy and fusion entropy. For the ILs with short alkyl chains (ethyl and butyl groups), fusion entropy plays the dominant role so that [C4MIM][PF6], which has a larger fusion entropy due to its higher liquid phase entropy has the lower melting temperature. As for the ILs with long alkyl chains, due to the enhanced van der Waals interactions brought about by the long non-polar alkyl chains, enthalpy becomes the deciding factor and the melting points increase when the alkyl chain goes from C10 to C14. While the melting points for [C2MIM][PF6] and [C4MIM][PF6] were quantitatively predicted and the trends for the long chain ILs were captured correctly, the absolute melting points for [C10MIM][PF6], [C12MIM][PF6] and [C14MIM][PF6] were systematically overestimated in the simulations. Three possible reasons for the overestimation were studied but all ruled out. Further simulation or experimental studies are needed to explain the difference.