Simulations of ionic liquids, solutions, and surfaces

Geminal Dicationic Ionic Liquids (GDILs) and Their Adsorption on Graphene Nanoflakes

In this work, the configuration and stability of 15 geminal dicationic ionic liquids (GDILs) and their adsorption mechanism on the graphene nanoflake (GNF) are investigated using the density functional theory (DFT) method. We find that the interactions of dications ([DAm]+, [DIm]+, [DImDm]+, [DPy]+, and [DPyrr]+)) are stabilized near the anions ([BF4]-, [PF6]-, and [Tf2N]-) in the most stable configurations of GDILs through electrostatic interactions, van der Waals (vdW) interactions, and hydrogen bonding (H-bonding). Our calculations show that the adsorption of the GDILs on the GNF is consistent with the charge transfer and occurs via X···π (X = N, O, F), C-H···π, and π···π noncovalent interactions, leading to a decrease in the strength of the intermolecular interactions between the dications and anions in the GDILs. The thermochemistry calculations reveal that the formation of GDIL@GNF complexes is an exothermic and favorable reaction. The adsorption energy (Eads) calculations show that the highest Eads values for the interaction of GDILs containing [BF4]-, [PF6]-, and [Tf2N]- anions with the GNF are observed for the [DPy][BF4]@GNF (-23.56 kcal/mol), [DPy][PF6]@GNF (-29.29 kcal/mol), and [DPyrr][Tf2N]@GNF (-24.74 kcal/mol) complexes, respectively. Our results show that the adsorption of the GDILs on the GNF leads to the decrease of the chemical potential (μ), chemical hardness (η), and HOMO-LUMO energy gap (Eg) values and an increase in the electrophilicity index (ω) value of the GNF. In addition, the effect of GDIL adsorption on the UV-vis absorption spectrum was studied at the TD-M06-2X/cc-pVDZ level of theory. We find that the adsorption of GDILs results in minimal change in the shape of the main absorption peak (at λ = 363 nm) in the GNF spectrum and only shifts it to higher wavelengths. On the other hand, a new peak appears in the GNF spectrum upon adsorption of [DPy][Y] (Y = [BF4]-, [PF6]-, and [Tf2N]-) due to the relatively strong π···π interactions between the [DPy]+ dication and GNF. Finally, the transition density matrix (TDM) heat maps show that electron transfers related to the excitation states in the GDIL@GNF complexes occur mainly through π(C=C) → π*(C=C) transitions in the GNF and the transitions from [DPy]+ dication to the GNF.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Structural and Dynamical Properties of Liquids in Confinements: A Review of Molecular Dynamics Simulation Studies

Molecular dynamics (MD) simulations are a powerful tool for detailed studies of altered properties of liquids in confinement, in particular, of changed structures and dynamics. They allow, on one hand, for perfect control and systematic variation of the geometries and interactions inherent in confinement situations and, on the other hand, for type-selective and position-resolved analyses of a huge variety of structural and dynamical parameters. Here, we review MD simulation studies on various types of liquids and confinements. The main focus is confined aqueous systems, but also ionic liquids and polymer and silica melts are discussed. Results for confinements featuring different interactions, sizes, shapes, and rigidity will be presented. Special attention will be given to situations in which the confined liquid and the confining matrix consist of the same type of particles and, hence, disparate liquid-matrix interactions are absent. Findings for the magnitude and the range of wall effects on molecular positions and orientations and on molecular dynamics, including vibrational motion and structural relaxation, are reviewed. Moreover, their dependence on the parameters of the confinement and their relevance to theoretical approaches to the glass transition are addressed.

Insights into Ionic Liquids: From Z-Bonds to Quasi-Liquids

Ionic liquids (ILs) hold great promise in the fields of green chemistry, environmental science, and sustainable technology due to their unique properties, such as a tailorable structure, the various types available, and their environmentally friendly features. On the basis of multiscale simulations and experimental characterizations, two unique features of ILs are as follows: (1) strong coupling interactions between the electrostatic forces and hydrogen bonds, namely in the Z-bond, and (2) the unique semiordered structure and properties of ultrathin films, specifically regarding the quasi-liquid. In accordance with the aforementioned theoretical findings, many cutting-edge applications have been proposed: for example, CO₂ capture and conversion, biomass conversion and utilization, and energy storage materials. Although substantial progress has been made recently in the field of ILs, considerable challenges remain in understanding the nature of and devising applications for ILs, especially in terms of e.g. in situ/real-time observation and highly precise multiscale simulations of the Z-bond and quasi-liquid. In this Perspective, we review recent developments and challenges for the IL research community and provide insights into the nature and function of ILs, which will facilitate future applications.

Carbon-carbon supercapacitors: Beyond the average pore size or how electrolyte confinement and inaccessible pores affect the capacitance

Carbon-carbon supercapacitors are high power electrochemical energy storage systems, which store energy through reversible ion adsorption at the electrode-electrolyte interface. Due to the complex structure of the porous carbons used as electrodes, extracting structure-property relationships in these systems remains a challenge. In this work, we conduct molecular simulations of two model supercapacitors based on nanoporous electrodes with the same average pore size, a property often used when comparing porous materials, but different morphologies. We show that the carbon with the more ordered structure, and a well defined pore size, has a much higher capacitance than the carbon with the more disordered structure and a broader pore size distribution. We analyze the structure of the confined electrolyte and show that the ions adsorbed in the ordered carbon are present in larger quantities and are also more confined than for the disordered carbon. Both aspects favor a better charge separation and thus a larger capacitance. In addition, the disordered electrodes contain a significant amount of carbon atoms, which are never in contact with the electrolyte, carry a close to zero charge, and are thus not involved in the charge storage. The total quantities of adsorbed ions and degrees of confinement do not change much with the applied potential, and as such, this work opens the door to computationally tractable screening strategies.

Dynamical properties across different coarse-grained models for ionic liquids

Room-temperature ionic liquids (RTILs) stand out among molecular liquids for their rich physicochemical characteristics, including structural and dynamic heterogeneity. The significance of electrostatic interactions in RTILs results in long characteristic length- and timescales, and has motivated the development of a number of coarse-grained (CG) simulation models. In this study, we aim to better understand the connection between certain CG parameterization strategies and the dynamical properties and transferability of the resulting models. We systematically compare five CG models: a model largely parameterized from experimental thermodynamic observables; a refinement of this model to increase its structural accuracy; and three models that reproduce a given set of structural distribution functions by construction, with varying intramolecular parameterizations and reference temperatures. All five CG models display limited structural transferability over temperature, and also result in various effective dynamical speedup factors, relative to a reference atomistic model. On the other hand, the structure-based CG models tend to result in more consistent cation-anion relative diffusion than the thermodynamic-based models, for a single thermodynamic state point. By linking short- and long-timescale dynamical behaviors, we demonstrate that the varying dynamical properties of the different CG models can be largely collapsed onto a single curve, which provides evidence for a route to constructing dynamically-consistent CG models of RTILs.

Structural transitions at electrodes, immersed in simple ionic liquid models

We used a recently developed classical Density Functional Theory (DFT) method to study the structures, phase transitions, and electrochemical behaviours of two coarse-grained ionic fluid models, in the presence of a perfectly conducting model electrode. Common to both is that the charge of the cationic component is able to approach the electrode interface more closely than the anion charge. This means that the cations are specifically attracted to the electrode, due to surface polarization effects. Hence, for a positively charged electrode, there is competition at the surface between cations and anions, where the latter are attracted by the positive electrode charge. This generates demixing, for a range of positive voltages, where the two phases are structurally quite different. The surface charge density is also different between the two phases, even at the same potential. The DFT formulation contains an approximate treatment of ion correlations, and surface polarization, where the latter is modelled via screened image interactions. Using a mean-field DFT, where ion correlations are neglected, causes the phase transition to vanish for both models, but there is still a dramatic drop in the differential capacitance as proximal cations are replaced by anions, for increasing surface potentials. While these findings were obtained for relatively crude coarse-grained models, we argue that the findings can also be relevant in "real" systems, where we note that many ionic liquids are composed of a spherically symmetric anion, and a cation that is asymmetric both from a steric and a charge distribution point of view.

Structure, Molecular Interactions, and Dynamics of Aqueous [BMIM][BF4] Mixtures: A Molecular Dynamics Study

Molecular dynamics simulations with many-body polarizable force fields were carried out to investigate the thermodynamic, structural, and dynamic properties of aqueous solutions of 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF₄]). The radial distribution functions exhibit well-defined features, revealing favored structural correlations between [bmim]⁺, [BF₄]^-, and H₂O. The addition of water is shown to alter ionic liquid structural organizations by replacing counterions in the coordination shells and disrupt the cation-anion network. At low water concentration, the majority of water molecules are isolated from each other and have lower average dipole moment than that in pure water. With increasing hydration level, while [bmim][BF₄] ionic network breaks up and becomes isolated ion pairs or free ions in the dilute limit, water begins to form clusters of increasing sizes and eventually forms a percolating network. As a result, the average water dipole moment increases and approaches its bulk value. Water is also observed to have a substantial influence on the dynamics of ionic liquids. At low water content, the cation and anion have similar diffusion coefficients due to the correlated ionic motion of long-lived ion pairs. As the water concentration increases, both ions exhibit greater mobility and faster rotations from the breakup of ionic network. Consequently, the ionic conductivity of [bmim][BF₄] aqueous solutions rises with increasing water composition.

Hysteresis in the MD Simulations of Differential Capacitance at the Ionic Liquid-Au Interface

In this Letter, we report the first observation of the capacitance-potential hysteresis at the ionic liquid | electrode interface in atomistic molecular dynamics simulations. While modeling the differential capacitance dependence on the potential scan direction, we detected two long-living types of interfacial structure for the BMImPF₆ ionic liquid at specific charge densities of the gold Au(111) surface. These structures differ in how counterions overscreen the surface charge. The high barrier for the transition from one structure to another slows down the interfacial restructuring process and leads to the marked capacitance-potential hysteresis.

The structure and interaction properties of two task-specific ionic liquids and acetonitrile mixtures: A combined FTIR and DFT study

The mixtures of ionic liquid (IL) and acetonitrile (CH₃CN) can be used as reaction media, supercapacitors and thermally stable electrolytes. The macroscopic properties of ILs-CH₃CN mixtures have been extensively studied. However, some fundamental questions regarding the microscopic properties of ILs-CH₃CN mixtures still remain to be answered. In this work, the structure properties and hydrogen-bond interactions of two task-specific ILs, i.e., 1-propylnitrile-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([PCNMIM][Tf₂N]) and 1-(2'-hydroxylethyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([C2OHMIM][Tf₂N]), and CH₃CN were studied using the combination of Fourier transform infrared spectroscopy (FTIR) and density functional theory (DFT) calculations. The aromatic C‒H stretching vibration region of the cation was an area of special focus. Excess infrared spectroscopy with enhanced resolution was applied to analyse the original infrared spectra. It is found that: (1) The two ILs form stable hydrogen-bonds with CH₃CN. (2) Ion cluster, ion cluster-acetonitrile, and ion pair-acetonitrile are identified in the mixture. Acetonitrile cannot break apart the strong electronic interaction between the cation and anion in the examined concentration range. (3) The hydrogen-bonds are weak strength, closed shell and electrostatic dominant interactions. (4) The preferred interaction site of [PCNMIM]⁺ cation is the hydrogen atom at the C2 site, while that of [C2OHMIM]⁺ cation is the hydrogen atom in the hydroxyl group.

A combination of FTIR and DFT to study the microscopic structure and hydrogen-bonding interaction properties of the [BMIM][BF4] and water

The structure and hydrogen-bond interaction property of water and a model ionic liquid (IL): 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF₄]) were studied using the combination of Fourier transform infrared spectroscopy (FTIR) and density functional theory (DFT) calculations. The O‒D stretching vibration region of the deuterated water was an area of special focus. Excess infrared spectroscopy with enhanced resolution was applied to analyse the original infrared spectra of v(O‒D). It is found that: (1) [BMIM][BF₄] forms stable hydrogen-bonds with water in the mixture. (2) The hydrogen-bonds are weak strength, closed shell and electrostatic dominant interactions. The preferred interaction site of [BMIM]⁺ cation is the hydrogen atom at the C2. (3) Cage hexamer water, cyclic tetramer water, cyclic trimer water, ion cluster-water complex, ion pair-water, and anion-water complexes are identified in the mixture. When the mole fraction of D₂O (x(D₂O)) is larger than 0.9, ion cluster and ion pair were broken apart into individual cations and anions. The cage hexamer water, cyclic tetramer water, and cyclic trimer water disappear at x(D₂O) < 0.8, 0.5, and 0.3, respectively. HDO formed by H/D isotope exchange was detected when x(D₂O) is less than 0.3.

Tuning Water Networks via Ionic Liquid/Water Mixtures

Water in nanoconfinement is ubiquitous in biological systems and membrane materials, with altered properties that significantly influence the surrounding system. In this work, we show how ionic liquid (IL)/water mixtures can be tuned to create water environments that resemble nanoconfined systems. We utilize molecular dynamics simulations employing ab initio force fields to extensively characterize the water structure within five different IL/water mixtures: [BMIM+][BF4−], [BMIM+][PF6−], [BMIM+][OTf−], [BMIM+][NO3−] and [BMIM+][TFSI−] ILs at varying water fraction. We characterize water clustering, hydrogen bonding, water orientation, pairwise correlation functions and percolation networks as a function of water content and IL type. The nature of the water nanostructure is significantly tuned by changing the hydrophobicity of the IL and sensitively depends on water content. In hydrophobic ILs such as [BMIM+][PF6−], significant water clustering leads to dynamic formation of water pockets that can appear similar to those formed within reverse micelles. Furthermore, rotational relaxation times of water molecules in supersaturated hydrophobic IL/water mixtures indicate the close-connection with nanoconfined systems, as they are quantitatively similar to water relaxation in previously characterized lyotropic liquid crystals. We expect that this physical insight will lead to better design principles for incorporation of ILs into membrane materials to tune water nanostructure.

Screening of highly charged ions in an ionic liquid; when will ion pairs form?

The properties of pairs of doubly charged solute ions are studied as a function of their separation in the ionic liquid, dimethylimidazolium chloride ([dmim][Cl]).

Formation of Aqueous Biphasic Systems with an Ionic Liquid Induced by Metallic Salts: Nanoscopic Views from Molecular Dynamics Simulations

The formation of aqueous biphasic systems (ABSs) based on aqueous ionic liquid (IL)/salt mixtures has been investigated via molecular dynamics simulations (with IL butyl-methyl-imidazolium triflate; salts NaCl, CsCl, SrCl₂, and EuCl₃). The analysis of ion distributions, solvation, and mutual interactions during the dynamics reveals the heterogeneity of all solutions due to ion segregation into mutually exclusive IL and salt domains, even in monophasic solutions ("ionic sociology"). Ion segregation and ABS formation are found to increase with (i) the salt content and (ii) the IL content, (iii) in the order Na⁺ < Sr²⁺ < Eu³⁺, and (iv) when the IL ion "polarity" is diminished, following experimental trends. The structuration of the solution is rationalized as a synergistic water transfer from the best donating ion pair (first hydration shell of hydrophobic moieties of IL ions) to the best accepting pair (M ⁿ⁺ and Cl^- ions, beyond their first shell). In ABSs, the IL- and salt-containing phases are linked by a well-defined "interface" that decreases in width when MCl _n becomes more hydrophilic and/or more concentrated. In the IL-rich phase of ABSs, the hydration of IL ions and their mutual interactions are shown to be similar to those displayed at aqueous interfaces.

Effect of Concentration on the Interfacial and Bulk Structure of Ionic Liquids in Aqueous Solution

Bio and aqueous applications of ionic liquids (IL) such as catalysis in micelles formed in aqueous IL solutions or extraction of chemicals from biologic materials rely on surface-active and self-assembly properties of ILs. Here, we discuss qualitative relations of the interfacial and bulk structuring of a water-soluble surface-active IL ([C₈MIm][Cl]) on chemically controlled surfaces over a wide range of water concentrations using both force probe and X-ray scattering experiments. Our data indicate that IL structuring evolves from surfactant-like surface adsorption at low IL concentrations, to micellar bulk structure adsorption above the critical micelle concentration, to planar bilayer formation in ILs with <1 wt % of water and at high charging of the surface. Interfacial structuring is controlled by mesoscopic bulk structuring at high water concentrations. Surface chemistry and surface charges decisively steer interfacial ordering of ions if the water concentration is low and/or the surface charge is high. We also demonstrate that controlling the interfacial forces by using self-assembled monolayer chemistry allows tuning of interfacial structures. Both the ratio of the head group size to the hydrophobic tail volume as well as the surface charging trigger the bulk structure and offer a tool for predicting interfacial structures. Based on the applied techniques and analyses, a qualitative prediction of molecular layering of ILs in aqueous systems is possible.

Ionic liquid interface at an electrode: simulations of electrochemical properties using an asymmetric restricted primitive model

We use Monte Carlo simulations of a coarse-grained model to investigate structure and electrochemical behaviours at an electrode immersed in room temperature ionic liquids (RTILs). The simple RTIL model, which we denote the asymmetric restricted primitive model (ARPM), is composed of monovalent hard-sphere ions, all of the same size, in which the charge is asymmetrically placed. Not only the hard-sphere size (d), but also the charge displacement (b), is identical for all species, i.e. the monovalent RTIL ions are fully described by only two parameters (d, b). In earlier work, it was demonstrated that the ARPM can capture typical static RTIL properties in bulk solutions with remarkable accuracy. Here, we investigate its behaviour at an electrode surface. The electrode is assumed to be a perfect conductor and image charge methods are utilized to handle polarization effects. We find that the ARPM of the ionic liquid reproduces typical (static) electrochemical properties of RTILs. Our model predicts a declining differential capacitance with increasing temperature, which is expected from simple physical arguments. We also compare our ARPM, with the corresponding RPM description, at an elevated temperature (1000 K). We conclude that, even though ion pairing occurs in the ARPM system, reducing the concentration of 'free' ions, it is still better able to screen charge than a corresponding RPM melt. Finally, we evaluate the option to coarse-grain the model even further, by treating the fraction of the ions that form ion pairs implicitly, only through the contribution to the dielectric constant of the corresponding dipolar (ion pair) fluid. We conclude that this primitive representation of ion pairing is not able to reproduce the structures and differential capacitances of the system with explicit ion pairs. The main problem seems to be due to a limited dielectric screening in a layer near the electrode surface, resulting from a combination of orientational restrictions and a depleted dipole density.

Surface induced smectic order in ionic liquids - an X-ray reflectivity study of [C22C1im]+[NTf2]. Phys Chem Chem Phys 2018;19:26651-26661. [PMID: 28960006 DOI: 10.1039/c7cp04852a]

Surface induced smectic order was found for the ionic liquid 1-methyl-3-docosylimidazolium bis(trifluoromethlysulfonyl)imide by X-ray reflectivity and grazing incidence scattering experiments. Near the free liquid surface, an ordered structure of alternating layers composed of polar and non-polar moieties is observed. This leads to an oscillatory interfacial profile perpendicular to the liquid surface with a periodicity of 3.7 nm. Small angle X-ray scattering and polarized light microscopy measurements suggest that the observed surface structure is related to fluctuations into a metastable liquid crystalline SmA₂ phase that was found by supercooling the bulk liquid. The observed surface ordering persists up to 157 °C, i.e. more than 88 K above the bulk melting temperature of 68.1 °C. Close to the bulk melting point, we find a thickness of the ordered layer of L = 30 nm. The dependency of L(τ) = Λ ln(τ/τ₁) vs. reduced temperature τ follows a logarithmic growth law. In agreement with theory, the pre-factor Λ is governed by the correlation length of the isotropic bulk phase.

XPS enables visualization of electrode potential screening in an ionic liquid medium with temporal- and lateral-resolution

We present an X-ray photoelectron spectroscopic (XPS) investigation of potential screening across two gold electrodes fabricated on a porous polymer surface which is impregnated with the ionic liquid (IL) N-N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide [DEME-TFSI]. The IL provides a sheet of conducting layers to the insulating polymer film, and allows monitoring charging and screening dynamics at the polymer + IL/vacuum interface in a laterally resolved fashion across the electrodes. Time-resolved measurements are also implemented by recording F1s peaks of the IL, while imposing 10 mHz square-wave (SQW) pulses across the two electrodes in a source-drain geometry. Variations in the F1s binding energy reflect directly the transient local electrical potential, and allow us to visualize screening of the otherwise built-in local voltage drop on and across the metal electrodes in the range of millimeters. Accordingly, the device is partitioned into two oppositely polarized regions, each following polarization of one electrode through the IL medium. On the other extreme, upon imposing relatively fast 1 kHz SQW pulses the charge screening is prevented and the device is brought to assume a simple resistor role. A simple equivalent circuit model also reproduces the observed voltage transients qualitatively. The presented structure and variants of XPS measurements, enabling us to record voltage transients in unexpectedly large lateral distances away from the electrodes, can impact the understanding of various electrochemical concepts.

Molecular scale structure and dynamics at an ionic liquid/electrode interface

The structural arrangement and dynamics of ions near the IL/electrode interface during charging and discharging was studied by a combination of time resolved X-ray reflectivity and impedance spectroscopy.

On the thickness of the double layer in ionic liquids

The Au(111)|BF₄⁻interface model in which BF₄⁻reorients and spontaneously dissociates at surface coverageθ= 1/3.

Molecular Dynamics Studies on Liquid/Vapor Interface Properties and Structures of 1-Ethyl-3-methylimidazolium Dimethylphosphate-Water

1-Ethyl-3-methylimidazolium dimethylphosphate ([Emim][Dmp])-water binary solution is one of the promising new working-pairs for absorption heat pump and absorption chillers, which are widely used to recover industrial waste heat. In the absorption process, the mass and heat transfer at the interface greatly depend on interface microscopic structure. Therefore, in order to understand the absorption process, it is very important to study the interface microscopic structure. The liquid-vapor interface properties, as well as the orientation of [Emim]⁺, [Dmp]^-, water at the interface and its aqueous solution with different water mole fraction, were studied using classical all-atom force field by molecular dynamic simulations. The simulated bulk mass density fitted by hyperbolic tangent function for each system was in good agreement with the experiment data, with the relative deviation between simulated and experimental value within 2%. The simulated results indicate that anion is always distributed at the outmost layer of the interface, followed by cation and water molecule. In [Emim][Dmp], the tilt angle of imidazolium rings to the surface normal is in the range of 0° < θ < 12°; for most cation, their ethyl and methyl tilted toward gas phase and bulk, respectively, but for a few cation, their ethyl and the methyl take the opposite orientation. For anion, one methyl prefers to turn toward gas phase and another methyl (PC vector from P atom to C atom) lie nearly parallel to the surface, while one PO vector (from P atom to O atom) turns toward liquid bulk and another PO vector is nearly parallel to the surface. In aqueous solution of [Emim][Dmp], the tilt angle of the imidazolium ring to the surface normal becomes larger (0° < θ < 37°) at the interface, but almost all ethyl intend to tilt toward gas phase and the methyl tilt toward liquid bulk compared with pure [Emim][Dmp]. Two methyl in anion prefer to turn toward gas phase and its two PO vectors toward liquid bulk. This orientation indicates that pure [Emim][Dmp] absorb water in gas phase more easily than [Emim][Dmp]+H₂O system does. Water molecules are distributed in the inner layer of the interface with two OH vectors (from O atom to H atom) tilting toward external surface.

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.

Initial dissolution of D2O at the gas-liquid interface of the ionic liquid [C4min][NTf2] associated with hydrogen-bond network formation

We have studied the initial dissolution of D₂O at the interfacial surface of the flowing jet sheet beam of the ionic liquid (IL) [C₄min][NTf2] using the King and Wells method as a function of both the temperature and collision energy of the IL. The initial dissolution probability of D₂O into the IL [C₄min][NTf2] was found to follow the general propensity that the solubility of gases into a liquid decreases with temperature. However, a large partial molar enthalpy and entropy for the initial dissolution of D₂O in the IL [C₄min][NTf2] were observed from the temperature dependence of the initial dissolution probability: ΔH_l = -53 ± 8 kJ mol^-1, ΔS_l = -210 ± 30 J mol^-1 K^-1. In addition, it was found that the collision energy significantly reduced the initial dissolution probability. We propose that the associated D₂O molecules at the interface of the IL [C₄min][NTf2] make a hydrogen-bond network around the [NTf2]^- anion before dissolution into the deeper portion of the interface layer.

Structures of Ionic Liquids Having Both Anionic and Cationic Octyl Tails: Lamellar Vacuum Interface vs Sponge-Like Bulk Order

Numerous experimental and computational studies have shown that the structure of ionic liquids is significantly influenced by confinement and by interactions with interfaces. The nature of the interface can affect the immediate ordering of cations and anions, changing important rheological characteristics relevant to lubrication. Most studies suggest that such changes are local or short-ranged and that bulk properties are reestablished on a length scale of a few nanometers. The current study focuses on the 1-methyl-3-octylimidazolium octylsulfate ionic liquid for which both the cation and anion have moderate length linear alkyl tails. For this system, we find that the bulk phase is dominated by the very common sponge-like morphology characteristic of many ionic liquids. However, at the vacuum interface, a lamellar structure is observed that is not restricted to the vicinity of the surface but instead extends across the full 9 nm slab of our simulation. We suspect that in reality it could extend significantly beyond this.

Computational and Experimental Study of Li-Doped Ionic Liquids at Electrified Interfaces

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Ψ). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Ψ. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Ψ-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Ψ and suggest a source of impedance for Li+ to reach the electrode. Finally, we make predictions of the specific energy at ideal graphite utilizing the computed capacitance and previously derived electrochemical windows of the liquids.

Dual analyzer system for surface analysis dedicated for angle-resolved photoelectron spectroscopy at liquid surfaces and interfaces

The investigation of liquid surfaces and interfaces with the powerful toolbox of ultra-high vacuum (UHV)-based surface science techniques generally has to overcome the issue of liquid evaporation within the vacuum system. In the last decade, however, new classes of liquids with negligible vapor pressure at room temperature-in particular, ionic liquids (ILs)-have emerged for surface science studies. It has been demonstrated that particularly angle-resolved X-ray Photoelectron Spectroscopy (ARXPS) allows for investigating phenomena that occur at gas-liquid and liquid-solid interfaces on the molecular level. The results are not only relevant for IL systems but also for liquids in general. In all of these previous ARXPS studies, the sample holder had to be tilted in order to change the polar detection angle of emitted photoelectrons, which restricted the liquid systems to very thin viscous IL films coating a flat solid support. We now report on the concept and realization of a new and unique laboratory "Dual Analyzer System for Surface Analysis (DASSA)" which enables fast ARXPS, UV photoelectron spectroscopy, imaging XPS, and low-energy ion scattering at the horizontal surface plane of macroscopically thick non-volatile liquid samples. It comprises a UHV chamber equipped with two electron analyzers mounted for simultaneous measurements in 0° and 80° emission relative to the surface normal. The performance of DASSA on a first macroscopic liquid system will be demonstrated.

A new QM/MM method oriented to the study of ionic liquids

The interest on room temperature ionic liquids has grown in the last decades because of their use as all-purpose solvent and their low environmental impact. In the present work, a new theoretical procedure is developed to study pure ionic liquids within the framework of the quantum mechanics/molecular mechanics method. Each type of ion (cation or anion) is considered as an independent entity quantum mechanically described that follows a differentiated path in the liquid. The method permits, through an iterative procedure, the full coupling between the polarized charge distribution of the ions and the liquid structure around them. The procedure has been tested with 1-ethyl-3-methylimidazolium tetrafluoroborate. It was found that, similar to non-polar liquids and as a consequence of the low value of the reaction field, the cation and anion charge distributions are hardly polarized by the rest of molecules in the liquid. Their structure is characterized by an alternance between anion and cation shells as evidenced by the coincidence of the first maximum of the anion-anion and cation-cation radial distribution functions with the first minimum of the anion-cation. Some degree of stacking between the cations is also found.