A theoretical and experimental chemist's joint view on hydrogen bonding in ionic liquids and their binary mixtures

Adsorption of Tunable aryl alkyl ionic liquids (TILs) on the graphene and Defective graphene nanosheets: A DFT Study

Optical and electronic characteristics of the graphene nanosheets (GNS) could be altered by some structural defects such as double-vacancy and Stone-Wales ones. The physisorption manner of [MPI][BF4], [MPT1][BF4], [MPT2][BF4], and [MPTT][BF4] ionic liquids on intact and defective GNS surfaces were investigated using M06-2X/cc-pVDZ computational method. Capability for adsorption on the DV and SW graphene surfaces by TILs is increased by about 1.0-4.3 and 0.4-2.0 kcal/mol respectively. The electrostatic potential of the GNS-DV surface is more negative than the GNS-SW one which enables it to interact with cation parts of the adsorbed TILs so extensively. The highest adsorption energy belongs to the [MPI][BF4]/GNS-DV system. Adsorption of the TILs on the GNS surfaces leads to a decrease in the energy of the LUMO molecular orbital as well as their energy gap of them. Results revealed that the electrical conductivity, as well as absorption spectra of the GNS surfaces, are affected by TILs adsorption and defect nature.

Raman Optical Activity of N-Acetyl-L-Cysteine in water and in methanol: the "clusters-in-a-liquid" model and ab initio molecular dynamics simulations

Raman and Raman Optical Activity (ROA) spectra of N-acetyl-L-cysteine (NALC), a flexible chiral molecule, were measured in water and in methanol to evaluate the solvent effects. Two different solvation approaches, i.e. the DFT based clusters-in-a-liquid solvent model and the ab initio molecular dynamics (AIMD) simulations, were applied to simulate the Raman and ROA spectra. Systematic conformational searches were carried out using a recently developed conformational searching tool, CREST, with the inclusion of polarizable continuum model of water and of methanol. The CREST candidates of NALC and the NALC-solvent complexes were re-optimized and their Raman and ROA simulations were done at the B3LYP-D3BJ/def2-TZVP and the B3LYP-aug-cc-pVDZ//cc-pVTZ levels. Also, AIMD simulations , which includes some anharmonic effects and all intermolecular interactions in solution, were performed. By empirically weighting the computed Raman and ROA spectra of each conformer, good agreements with the experimental data were achieved with both approaches, while AIMD offered some improvements in the carbonyl and in the low wavenumber regions over the static DFT approach. The pros and cons of these two different approaches for accounting the solvent effects on Raman and ROA of this flexible chiral system will also be discussed.

Optimized Atomic Partial Charges and Radii Defined by Radical Voronoi Tessellation of Bulk Phase Simulations

We present a novel method for the computation of well-defined optimized atomic partial charges and radii from the total electron density. Our method is based on a two-step radical Voronoi tessellation of the (possibly periodic) system and subsequent integration of the total electron density within each Voronoi cell. First, the total electron density is partitioned into the contributions of each molecule, and subsequently the electron density within each molecule is assigned to the individual atoms using a second set of atomic radii for the radical Voronoi tessellation. The radii are optimized on-the-fly to minimize the fluctuation (variance) of molecular and atomic charges. Therefore, our method is completely free of empirical parameters. As a by-product, two sets of optimized atomic radii are produced in each run, which take into account many specific properties of the system investigated. The application of an on-the-fly interpolation scheme reduces discretization noise in the Voronoi integration. The approach is particularly well suited for the calculation of partial charges in periodic bulk phase systems. We apply the method to five exemplary liquid phase simulations and show how the optimized charges can help to understand the interactions in the systems. Well-known effects such as reduced ion charges below unity in ionic liquid systems are correctly predicted without any tuning, empiricism, or rescaling. We show that the basis set dependence of our method is very small. Only the total electron density is evaluated, and thus, the approach can be combined with any electronic structure method that provides volumetric total electron densities-it is not limited to Hartree-Fock or density functional theory (DFT). We have implemented the method into our open-source software tool TRAVIS.

Probing the Effect of Side Alkyl Chain Length on the Structural and Dynamical Micro-heterogeneities in Dicationic Ionic Liquids

The molecular dynamics simulations and Voronoi tessellation analysis of two dicationic ionic liquids (DILs) including [C₅(mim)₂][NTf₂]₂ and [C₅(mim)₂C₄][NTf₂]₂ have been carried out to investigate the effects of side alkyl chain length on the structural and dynamical micro-heterogeneity of these DILs. Radial distribution functions (RDFs), spatial distribution functions (SDFs), and also neighborhood analysis of ions have been calculated to determine the arrangement of the nearest neighboring ions. To better understand the hydrogen-bonding network, microstructures, inter- and intramolecular orientations of ions in the studied DILs, different kinds of combined distribution functions (CDFs) were computed and analyzed. Also, qualitative and quantitative analyses of the structural heterogeneity were explored through total/partial structure factors, heterogeneity order parameters (HOPs), and domain analysis from Voronoi tessellation. The results showed that the side alkyl chains in DILs have significant effects on their micro-organizations in such a way that [C₅(mim)₂C₄][NTf₂]₂ with longer side chains has more microstructural heterogeneity than [C₅(mim)₂][NTf₂]₂ where the linkage alkyl chain is the same in both of them. Furthermore, to shed light on the dynamical heterogeneity, ion pair, ion cage, and hydrogen-bond stabilities and also the reorientation dynamics of ions have been investigated. Results demonstrated that local dynamics differences originate from local structural heterogeneity.

Aqueous solutions of binary ionic liquids: insight into structure, dynamics, and interface properties by molecular dynamics simulations and DFT methods

The behavior of aqueous solutions of mixtures of ionic liquids (ILs) is of special interest because of their amphiphilic character, from both a fundamental and application viewpoint. In this work, we conducted molecular dynamics (MD) simulations and density functional theory (DFT) calculations to understand the effect of water on the intermolecular interactions in three IL binary mixtures [C4mim]/[Cl]/[BF4], [C4mim]/[Cl]/[PF6] and [C4mim]/[BF4]/[PF6] containing the well-characterized cation, 1-n-butyl-3-methylimidazolium [C4mim]+ and the anions chloride [Cl]-, tetrafluoroborate [BF4]-, and hexafluorophosphate [PF6]-. The perturbation of the structures in the binary IL mixture by water molecules was analyzed in the bulk and at the liquid/vacuum interface using distribution functions, hydrogen-bond statistics, and density profiles. Interactions between anions and cations change drastically when the IL mixtures are dissolved in water. In particular, anion-water interactions are stronger than anion-cation interactions. H-Bonds are the dominant interactions. They are prevalently electrostatic and strong for the two [Cl]-containing systems in both the water-free and the water-containing systems. The very hydrophobic [C4mim]/[BF4]/[PF6] system gains stability from dispersive interactions and consequently segregates water markedly when admixed. The most probable orientations of IL cations in the bulk and at the vicinity of the interface were examined using bivariate distribution calculations and show [PF6]- segregating to the surface in keeping with its highly hydrophobic nature. DFT calculated structures, energies, dipole moments, global hardness and solvation energies using model ion pairs [C4mim][X] or complexes [C4mim]2[X][Y], with [X/Y]- = [Cl]-, [BF4]-, or [PF6]- are completely consistent with the findings for the bulk.

Exploring Free Energy Profiles of Enantioselective Organocatalytic Aldol Reactions under Full Solvent Influence

We present a computational study on the enantioselectivity of organocatalytic proline-catalyzed aldol reactions between aldehydes in dimethylformamide (DMF). To explore the free energy surface of the reaction, we apply two-dimensional metadynamics on top of ab initio molecular dynamics (AIMD) simulations with explicit solvent description on the DFT level of theory. We avoid unwanted side reactions by utilizing our newly developed hybrid AIMD (HyAIMD) simulation scheme, which adds a simple force field to the AIMD simulation to prevent unwanted bond breaking and formation. Our condensed phase simulation results are able to nicely reproduce the experimental findings, including the main stereoisomer that is formed, and give a correct qualitative prediction of the change in syn:anti product ratio with different substituents. Furthermore, we give a microscopic explanation for the selectivity. We show that both the explicit description of the solvent and the inclusion of entropic effects are vital to a good outcome-metadynamics simulations in vacuum and static nudged elastic band (NEB) calculations yield significantly worse predictions when compared to the experiment. The approach described here can be applied to a plethora of other enantioselective or organocatalytic reactions, enabling us to tune the catalyst or determine the solvent with the highest stereoselectivity.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Distance Angle Descriptors of the Interionic and Ion-Solvent Interactions in Imidazolium-Based Ionic Liquid Mixtures with Aprotic Solvents: A Molecular Dynamics Simulation Study

The aim of this paper is to quantify the changes of the interionic and ion-solvent interactions in mixtures of imidazolium-based ionic liquids, having tetrafluoroborate (BmimBF₄), hexafluorophosphate (BmimPF₆), trifluoromethylsulfonate (BmimTFO), or bis(trifluoromethanesulfonyl)imide (BmimTFSI), anions, and polar aprotic molecular solvents, such as acetonitrile (AN), γ-butyrolactone (GBL), and propylene carbonate (PC). For this purpose, we calculate, using the nearest-neighbor approach, the average distance between the imidazolium ring H atom in positions 2, 4, and 5 (H^2,4,5) and the nearest high-electronegativity atom of the solvent or anion (X) as distance descriptors, and the mean angle formed by the C^2,4,5-H^2,4,5 bond and the H^2,4,5···X axis around the H^2,4,5 atom as angular descriptors of the cation-anion and cation-solvent interactions around the ring C-H groups. The behavior of these descriptors as a function of the ionic liquid mole fraction is analyzed in detail. The obtained results show that the extent of the change of these descriptors with respect to their values in the neat ionic liquid depends both on the nature of the anion and on the mixture composition. Thus, in the case of the mixtures of the molecular solvents with BmimBF₄ and BmimTFO, a small change of the distance and a drastic increase of the angular descriptor corresponding to the cation-anion interactions are observed with decreasing mole fraction of the ionic liquid, indicating that the anion moves from the above/below position (with respect to the imidazolium ring plane) to a position that is nearly linearly aligned with the C²-H² bond and hinders the possible interaction between the C²-H² group and the solvent molecules. On the other hand, in the case of mixtures of BmimTFSI and BmimPF₆ with the molecular solvents, both the observed increase of the distance descriptor and the slight change of the angular descriptor with decreasing ionic liquid mole fraction are compatible with the direct interactions of the solvent with the C²-H² group. The behavior of these descriptors is correlated with the experimentally observed ¹H chemical shift of the C²-H² group and the red shift of the C²-H² vibrational mode, particularly at low ionic liquid mole fractions. The present results are thus of great help in interpreting these experimental observations.

Microheterogeneity in Ionic Liquid Mixtures: Hydrogen Bonding, Dispersed Ions, and Dispersed Ion Clusters

Mixtures of ionic liquids (ILs) having a common ion but differing in the identity of the anion or cation represent highly interesting media. By varying the composition, one can successfully modulate specific physicochemical, structural, and biological properties. The molecular interactions (coulombic, hydrogen-bonding, van der Waals, and π–π intermolecular forces) that determine the three-dimensional structure of pure ILs can indeed be modified by the addition of another IL. In this context, we present here a 1H NMR, Fourier transform (FT)-IR, thermogravimetric, and solvatochromic study of the structural features of IL binary mixtures based on a common imidazolium cation ([CnC1im]+) and anions of different size and hydrogen-bond acceptor ability. For each mixture, the analyses were carried out at different molar ratios of the two components.

Noncovalent interactions underlying binary mixtures of amino acid based ionic liquids: insights from theory

Mixtures of ionic liquids formed by blending a common 1-methyl-3-butylimidazolium [Bmim] cation with the dicarboxylic amino acid anions viz., aspartic acid [Asp], asparagine [Asn], glutamic acid [Glu], and glutamine [Gln], have been investigated by employing dispersion corrected density functional theory.

Local Structure in Terms of Nearest-Neighbor Approach in 1-Butyl-3-methylimidazolium-Based Ionic Liquids: MD Simulations

Description of the local microscopic structure in ionic liquids (ILs) is a prerequisite to obtain a comprehensive understanding of the influence of the nature of ions on the properties of ILs. The local structure is mainly determined by the spatial arrangement of the nearest neighboring ions. Therefore, the main interaction patterns in ILs, such as cation-anion H-bond-like motifs, cation-cation alkyl tail aggregation, and ring stacking, were considered within the framework of the nearest-neighbor approach with respect to each particular interaction site. We employed classical molecular dynamics (MD) simulations to study in detail the spatial, radial, and orientational relative distribution of ions in a set of imidazolium-based ILs, in which the 1-butyl-3-methylimidazolium (C4mim(+)) cation is coupled with the acetate (OAc(-)), chloride (Cl(-)), tetrafluoroborate (BF4(-)), hexafluorophosphate (PF6(-)), trifluoromethanesulfonate (TfO(-)), or bis(trifluoromethanesulfonyl)amide (TFSA(-)) anion. It was established that several structural properties are strongly anion-specific, while some can be treated as universally applicable to ILs, regardless of the nature of the anion. Namely, strongly basic anions, such as OAc(-) and Cl(-), prefer to be located in the imidazolium ring plane next to the C-H(2/4-5) sites. By contrast, the other four bulky and weakly coordinating anions tend to occupy positions above/below the plane. Similarly, the H-bond-like interactions involving the H(2) site are found to be particularly enhanced in comparison with the ones at H(4-5) in the case of asymmetric and/or more basic anions (C4mimOAc, C4mimCl, C4mimTfO, and C4mimTFSA), in accordance with recent spectroscopic and theoretical findings. Other IL-specific details related to the multiple H-bond-like binding and cation stacking issues are also discussed in this paper. The secondary H-bonding of anions with the alkyl hydrogen atoms of cations as well as the cation-cation alkyl chain aggregation turned out to be poorly sensitive to the nature of the anion.

Order in the chaos: the secret of the large negative entropy of dissolving 1-alkyl-3-methylimidazolium chloride in trihexyltetradecylphosphonium chloride

A large symmetric ion cluster cage is formed if 1-alkyl-3-methylimidazolium chloride ionic liquids with short alkyl chains are added to trihexyltetradecylphosphonium chloride which results in a large negative entropy of dissolution and an unexpected ion dynamics.

The covalence and infrared spectra of cationic hydrogen bonds and dihydrogen bonds

A theoretical study of hydrogen bonds and dihydrogen bonds formed by ethyl cation, hydrocarbons and magnesium hydride is presented with calculations performed at the BHandHLYP/6-31G(d,p) level of theory. The structural results and IR analyses demonstrated great insights, mainly the strengthening and weakness of the CC bond of the ethyl cation and π or pseudo-π bonds, respectively. The interaction strength was measured through the supermolecule as well as by means of additional approaches. The QTAIM calculations were applied to characterize not only the intermolecular interactions but specifically the covalent character in the H + ⋯ π, H + ⋯ pseudo-π and H + ⋯ H contacts. The NBO calculations were useful to interpret the polarization on the CC bond and whether this effect is related with the bond length reduction as well as increase of charge density and frequency shifts.