Switching the Structural Force in Ionic Liquid-Solvent Mixtures by Varying Composition

Wave mechanics in an ionic liquid mixture

Experimental measurements of interactions in ionic liquids and concentrated electrolytes over the past decade or so have revealed simultaneous monotonic and oscillatory decay modes. These observations have been hard to interpret using classical theories, which typically allow for just one electrostatic decay mode in electrolytes. Meanwhile, substantial progress in the theoretical description of dielectric response and ion correlations in electrolytes has illuminated the deep connection between density and charge correlations and the multiplicity of decay modes characterising a liquid electrolyte. The challenge in front of us is to build connections between the theoretical expressions for a pair of correlation functions and the directly measured free energy of interaction between macroscopic surfaces in experiments. Towards this aim, we here present measurements and analysis of the interactions between macroscopic bodies across a fluid mixture of two ionic liquids of widely diverging ionic size. The measured oscillatory interaction forces in the liquid mixtures are significantly more complex than for either of the pure ionic liquids, but can be fitted to a superposition of two oscillatory and one monotonic mode with parameters matching those of the pure liquids. We discuss this empirical finding, which hints at a kind of wave mechanics for interactions in liquid matter.

Counterion Distribution in the Stern Layer on Charged Surfaces

Counterion adsorption at the solid-liquid interface affects numerous applications. However, the counterion adsorption density in the Stern layer has remained poorly evaluated. Here we report the direct determination of surface charge density at the shear plane between the Stern layer and the diffuse layer. By the Grahame equation extension and streaming current measurements for different solid surfaces in different aqueous electrolytes, we are able to obtain the counterion adsorption density in the Stern layer, which is mainly related to the surface charge density but is less affected by the bulk ion concentration. The charge inversion concentration is further found to be sensitive to the ion type and ion valence rather than to the charged surface, which is attributed to the ionic competitive adsorption and ion-ion correlations. Our findings offer a framework for understanding ion distribution in many physical and chemical processes where the Stern layer is ubiquitous.

Integrative studies of ionic liquid interface layers: bridging experiments, theoretical models and simulations

Ionic liquids (ILs) are a class of salts existing in the liquid state below 100 °C, possessing low volatility, high thermal stability as well as many highly attractive solvent and electrochemical capabilities, etc., making them highly tunable for a great variety of applications, such as lubricants, electrolytes, and soft functional materials. In many applications, ILs are first either physi- or chemisorbed on a solid surface to successively create more functional materials. The functions of ILs at solid surfaces can differ considerably from those of bulk ILs, mainly due to distinct interfacial layers with tunable structures resulting in new ionic liquid interface layer properties and enhanced performance. Due to an almost infinite number of possible combinations among the cations and anions to form ILs, the diversity of various solid surfaces, as well as different external conditions and stimuli, a detailed molecular-level understanding of their structure-property relationship is of utmost significance for a judicious design of IL-solid interfaces with appropriate properties for task-specific applications. Many experimental techniques, such as atomic force microscopy, surface force apparatus, and so on, have been used for studying the ion structuring of the IL interface layer. Molecular Dynamics simulations have been widely used to investigate the microscopic behavior of the IL interface layer. To interpret and clarify the IL structure and dynamics as well as to predict their properties, it is always beneficial to combine both experiments and simulations as close as possible. In another theoretical model development to bridge the structure and properties of the IL interface layer with performance, thermodynamic prediction & property modeling has been demonstrated as an effective tool to add the properties and function of the studied nanomaterials. Herein, we present recent findings from applying the multiscale triangle "experiment-simulation-thermodynamic modeling" in the studies of ion structuring of ILs in the vicinity of solid surfaces, as well as how it qualitatively and quantitatively correlates to the overall ILs properties, performance, and function. We introduce the most common techniques behind "experiment-simulation-thermodynamic modeling" and how they are applied for studying the IL interface layer structuring, and we highlight the possibilities of the IL interface layer structuring in applications such as lubrication and energy storage.

The surface force balance: direct measurement of interactions in fluids and soft matter

Over the last half-century, direct measurements of surface forces have been instrumental in the exploration of a multitude of phenomena in liquid, soft, and biological matter. Measurements of van der Waals interactions, electrostatic interactions, hydrophobic interactions, structural forces, depletion forces, and many other effects have checked and challenged theoretical predictions and motivated new models and understanding. The gold-standard instrument for these measurements is thesurface force balance(SFB), orsurface forces apparatus, where interferometry is used to detect the interaction force and distance between two atomically smooth planes, with 0.1 nm resolution, over separations from about 1 µm down to contact. The measured interaction forcevs.distance gives access to the free energy of interaction across the fluid film; a fundamental quantity whose general form and subtle features reveal the underlying molecular and surface interactions and their variation. Motivated by new challenges in emerging fields of research, such as energy storage, biomaterials, non-equilibrium and driven systems, innovations to the apparatus are now clearing the way for new discoveries. It is now possible to measure interaction forces (and free energies) with control of electric field, surface potential, surface chemistry; to measure time-dependent effects; and to determine structurein situ. Here, we provide an overview the operating principles and capabilities of the SFB with particular focus on the recent developments and future possibilities of this remarkable technique.

Deep Eutectic Solvents for Efficient Drug Solvation: Optimizing Composition and Ratio for Solubility of β-Cyclodextrin

Deep eutectic solvents (DESs) show promise in pharmaceutical applications, most prominently as excellent solubilizers. Yet, because DES are complex multi-component mixtures, it is challenging to dissect the contribution of each component to solvation. Moreover, deviations from the eutectic concentration lead to phase separation of the DES, making it impractical to vary the ratios of components to potentially improve solvation. Water addition alleviates this limitation as it significantly decreases the melting temperature and stabilizes the DES single-phase region. Here, we follow the solubility of β-cyclodextrin (β-CD) in DES formed by the eutectic 2:1 mole ratio of urea and choline chloride (CC). Upon water addition to DES, we find that at almost all hydration levels, the highest β-CD solubility is achieved at DES compositions that are shifted from the 2:1 ratio. At higher urea to CC ratios, due to the limited solubility of urea, the optimum composition allowing the highest β-CD solubility is reached at the DES solubility limit. For mixtures with higher CC concentration, the composition allowing optimal solvation varies with hydration. For example, β-CD solubility at 40 wt% water is enhanced by a factor of 1.5 for a 1:2 urea to CC mole ratio compared with the 2:1 eutectic ratio. We further develop a methodology allowing us to link the preferential accumulation of urea and CC in the vicinity of β-CD to its increased solubility. The methodology we present here allows a dissection of solute interactions with DES components that is crucial for rationally developing improved drug and excipient formulations.

Anomalous Underscreening in the Restricted Primitive Model

Underscreening is a collective term for charge correlations in electrolytes decaying slower than the Debye length. Anomalous underscreening refers to phenomenology that cannot be attributed alone to steric interactions. Experiments with concentrated electrolytes and ionic fluids report anomalous underscreening, which so far has not been observed in simulation. We present Molecular Dynamics simulation results exhibiting anomalous underscreening that can be connected to cluster formation. A theory that accounts for ion pairing confirms the trend. Our results challenge the classic understanding of dense electrolytes impacting the design of technologies for energy storage and conversion.

Gelation, Clustering and Crowding in the Electrical Double Layer of Ionic Liquids

Understanding the bulk and interfacial properties of super-concentrated electrolytes, such as ionic liquids (ILs), has attracted significant attention lately for their promising applications in supercapacitors and batteries. Recently, McEldrew et al. developed a theory for reversible ion associations in bulk ILs, which accounted for the formation of all possible Cayley tree clusters and a percolating ionic network (gel). Here we adopt and develop this approach to understand the associations of ILs in the electrical double layer at electrified interfaces. With increasing charge of the electrode, the theory predicts a transition from a regime dominated by a gelled or clustered state to a crowding regime dominated by free ions. This transition from gelation to crowding is conceptually similar to the overscreening to crowding transition.

Structural Forces in Ionic Liquids: The Role of Ionic Size Asymmetry

Ionic liquids (ILs) are charged fluids composed of anions and cations of different size and shape. The ordering of charge and density in ILs confined between charged interfaces underlies numerous applications of IL electrolytes. Here, we analyze the screening behavior and the resulting structural forces of a representative IL confined between two charge-varied plates. Using both molecular dynamics simulations and a continuum theory, we contrast the screening features of a more-realistic asymmetric system and a less-realistic symmetric one. The ionic size asymmetry plays a nontrivial role in charge screening, affecting both the ionic density profiles and the disjoining pressure distance dependence. Ionic systems with size asymmetry are stronger coupled systems, and this manifests itself both in their response to the electrode polarization and spontaneous structure formation at the interface. Analytical expressions for decay lengths of the disjoining pressure are obtained in agreement with the pressure profiles computed from molecular dynamics simulations.

Ionic Liquids and Deep Eutectics as a Transformative Platform for the Synthesis of Nanomaterials

Ionic liquids (ILs) are becoming a revolutionary synthesis medium for inorganic nanomaterials, permitting more efficient, safer and environmentally benign preparation of high quality products. A smart combination of ILs and...

Controlling adhesion using AC electric fields across fluid films

We demonstrate reversible and switchable actuation using AC electric fields to bring two surfaces separated by a thin film of ionic fluid in and out of adhesive contact. Using a surface force balance we apply electric fields normal to a crossed-cylinder contact and measure directly the adhesive force and surface separation with sub-molecular resolution. Taking advantage of the oscillatory structural force acting between the surfaces across the fluid, which we show to be unaffected by the AC field, we pick between the distinct (quantized) adhesive states through precise tuning of the field. This proof-of-concept indicates exquisite control of surface interactions using an external field.

Bayesian unsupervised learning reveals hidden structure in concentrated electrolytes

Electrolytes play an important role in a plethora of applications ranging from energy storage to biomaterials. Notwithstanding this, the structure of concentrated electrolytes remains enigmatic. Many theoretical approaches attempt to model the concentrated electrolyte by introducing the idea of ion pairs, with ions either being tightly "paired" with a counter-ion or "free" to screen charge. In this study, we reframe the problem into the language of computational statistics and test the null hypothesis that all ions share the same local environment. Applying the framework to molecular dynamics simulations, we find that this null hypothesis is not supported by data. Our statistical technique suggests the presence of two distinct local ionic environments at intermediate concentrations, whose differences surprisingly originate in like charge correlations rather than unlike charge attraction. Through considering the effect of these "aggregated" and "non-aggregated" states on bulk properties including effective ion concentration and dielectric constant, we identify a scaling relation between the effective screening length and theoretical Debye length, which applies across different dielectric constants and ion concentrations.

Primitive model electrolytes in the near and far field: Decay lengths from DFT and simulations

Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT) and compare these with the decay lengths of the corresponding two-body correlation functions in bulk systems, obtained in previous Integral Equation Theory (IET) studies. Extensive Molecular Dynamics (MD) simulations are employed to complement the DFT and IET predictions. Our DFT calculations incorporate electrostatic interactions between the ions using three different (existing) approaches: one is based on the simplest mean-field treatment of Coulomb interactions (MFC), while the other two employ the Mean Spherical Approximation (MSA). The MSAc invokes only the MSA bulk direct correlation function, whereas the MSAu also incorporates the MSA bulk internal energy. Although MSAu yields profiles that are in excellent agreement with MD simulations in the near field, in the far field, we observe that the decay lengths are consistent between IET, MSAc, and MD simulations, whereas those from MFC and MSAu deviate significantly. Using DFT, we calculated the solvation force, which relates directly to surface force experiments. We find that its decay length is neither qualitatively nor quantitatively close to the large decay lengths measured in experiments and conclude that the latter cannot be accounted for by the primitive model. The anomalously large decay lengths found in surface force measurements require an explanation that lies beyond primitive models.

Structures and non-covalent interaction behaviours of binary systems containing the ionic liquid 1-(2'-hydroxylethyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide and chloroform

Mixtures of ionic liquids (ILs) and molecular solvents can overcome the drawbacks (high viscosity, high polarity, and high cost) of pure ILs and extend their practical use. The structural and interaction properties of ILs form the bases for understanding their properties. In this work, the structural properties of the mixtures of an IL, 1-(2'-hydroxylethyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([C2OHMIM][Tf₂N]), with chloroform, a molecular solvent of weak polarity, in various concentrations were analysed using Fourier transform infrared spectroscopy and density functional theory calculations. Excess spectra were used to analyse the infrared spectra. The IL forms a stable ion cluster-CDCl₃ complex with CDCl₃ in the concentration range investigated. In the ion cluster-CDCl₃ complex, the hydrogen atom of CDCl₃ forms hydrogen-bonds with the fluorine atoms of the anion. In addition, the chlorine atom of CDCl₃ forms a halogen-bond with the oxygen atom of the anion. All the hydrogen and halogen-bonds identified between the [C2OHMIM][Tf₂N] ion cluster and CDCl₃ exhibit low strength, closed shells, and electrostatically dominant interactions.

Water-induced mica/ionic liquid interfacial nanostructure switches revealed by AFM

Glovebox-AFM-based force curve measurements have been employed to investigate the effect of controlled small amounts of water on the interfacial structure of mica/a pyrrolidinium-based ionic liquid. A close examination reveals that with the increase of water content, the long-range monotonic force, which is beyond the region of the short-range oscillatory structure, switches from van der Waals attraction-dominated force to double layer repulsion-dominated force.

Strong Propensity of Ionic Liquids in Their Aqueous Solutions for an Organic-Modified Metal Surface

Understanding ionic structure and electrostatic environments near a surface has both fundamental and practical value. In electrochemistry, especially when room temperature ionic liquids (ILs) are involved, the complex ionic structure near the interface is expected to crucially influence reactions. Here we report evidence that even in dilute aqueous solutions of several ILs, the ions aggregate near the surface in ways that are qualitatively different from simple electrolytes. We have used a vibrational probe molecule, 4-mercaptobenzonitrile (MBN), tethered to a metal surface to monitor the behavior of the ionic layers. The characteristic nitrile vibrational frequency of this molecule has distinct values in the presence of pure water (∼2232 cm^-1) and pure IL (for example, ∼2226 cm^-1 for ethylmethylimidazolium tetrafluoroborate, [EMIM][BF₄]). This difference reflects the local electrostatic field and the hydrogen-bonding variations between these two limiting cases. We tracked this frequency shift as a function of IL concentration in water all the way from pure water to pure IL. We report two important findings. First, only one nitrile peak is observed for the entire concentration range, indicating that at least on the length scale of the probe molecule water and ILs do not phase separate within the interface, and no heterogeneously distinct electrostatic environments are formed. Second, and more importantly, we find that even up to a significant mole fraction of bulk water (x ∼ 0.95), the nitrile frequency does not change from that indicative of a pure IL for [EMIM][BF₄], indicating preferential aggregation of the ions near the surface. Because this behavior is very similar to surfactants, we chose an imidazolium cation with a longer side chain which resulted in behavior expected from a surfactant, with a preferential layer of the ions on the surface even in dilute water solutions (x ∼ 0.995). This observation indicates that even those ILs that are not nominally categorized as surfactants have a strong tendency to aggregate at the surface. Because ILs serve as electrolytes in a range of electrochemical reactions, including those requiring water, our results are likely useful for mechanistic understanding and tuning of such reactions.

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.

Uncovering the Underlying Mechanisms Governing the Solidlike Layering of Ionic Liquids (ILs) on Mica

Significant progress has been made in understanding the IL-solid interface in the past three decades, and a key finding is that ILs can form solidlike layers at the interface. It has been recognized that the electrostatic forces at the solid-IL interface and self-assembly of ILs are key enablers of the IL layering. However, regarding the layering structure of ILs, research from different laboratories is not consistent; i.e., the number of solidlike layers could range from 0 to ∼60, indicating the complexity of the underlying mechanisms and/or the existence of overlooked key parameters. In the current review, we will discuss the underlying mechanisms and key parameters governing the layering of ILs on mica, the most studied model solid. First, we will present the experimental findings from various laboratories, both consistent and contradictory ones, and summarize the current understanding of the governing mechanisms. Then, we will discuss the possible key parameters, including the structure of ILs, surface modification and contamination of mica, and cosolvent impacting the solidlike layering of ILs. Finally, we will discuss future research directions in uncovering the underlying mechanisms.

Correlation Length in Concentrated Electrolytes: Insights from All-Atom Molecular Dynamics Simulations

We study the correlation length of the charge-charge pair correlations in concentrated electrolyte solutions by means of all-atom, explicit-solvent molecular dynamics simulations. We investigate LiCl and NaI in water, which constitute highly soluble, prototypical salts for experiments, as well as two more complex, molecular electrolyte systems of lithium bis(trifluoromethane)sulfonimide (LiTFSI), a salt commonly employed in electrochemical storage systems, in water, and in an organic solvent mixture of dimethoxyethane and dioxolane. Our simulations support the recent experimental observations as well as theoretical predictions of a nonmonotonic behavior of the correlation length with increasing salt concentration. We observe a Debye-Hückel like regime at low concentration, followed by a minimum reached when d/λ_D ≃ 1, where λ_D is the Debye correlation length and d is the effective ionic diameter, and an increasing correlation length with salt concentration in very concentrated electrolytes. As in the experiments, we find that the screening length in the concentrated regime follows a universal scaling law as a function d/λ_D for all studied salts. However, the scaling exponent is significantly lower than the experimentally measured one and lies in the range of the theoretical predictions based on much simpler electrolyte models.

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.

A new methodology for a detailed investigation of quantized friction in ionic liquids

When confined at the nanoscale between smooth surfaces, an ionic liquid forms a structured film responding to shear in a quantized way, i.e., with a friction coefficient indexed by the number of layers in the gap. So far, only a few experiments have been performed to study this phenomenon, because of the delicate nature of the measurements. We propose a new methodology to measure friction with a surface force balance, based on the simultaneous application of normal and lateral motions to the surfaces, allowing for a more precise, comprehensive and rapid determination of the friction response. We report on proof-of-concept experiments with an ionic liquid confined between mica surfaces in dry or wet conditions, showing the phenomenon of quantized friction with an unprecedented resolution. First, we show that the variation of the kinetic friction force with the applied load for a given layer is not linear, but can be quantitatively described by two additive contributions that are respectively proportional to the load and to the contact area. Then, we find that humidity improves the resistance of the layers to be squeezed-out and extends the range of loads in which the liquid behaves as a superlubricant, interpreted by an enhanced dissolution of the potassium ions on the mica leading to a larger surface charge. There, we note a liquid-like friction behavior, and observe in certain conditions a clear variation of the kinetic friction force over two decades of shearing velocities, that does not obey a simple Arrhenius dynamics.

The structure and hydrogen-bond properties of N-alkyl-N-methyl-pyrrolidinium bis(trifluoromethylsulfonyl)imide and DMSO mixtures

Mixing ionic liquids (ILs) with molecular solvents can extend the practical applications of ILs and overcome the drawbacks of neat ILs. Knowledge on the structure and hydrogen-bond interaction properties of IL-molecular solvent mixtures is essential for chemical applications. In this work, the structure and hydrogen-bond features of N-alkyl-N-methyl-pyrrolidinium bis(trifluoromethylsulfonyl)imide ([CnMPyr][Tf2N], n = 3, 4, 6 and 8) and DMSO mixtures were studied using Fourier transform infrared spectroscopy (FTIR) and density functional theory (DFT) calculations. Excess infrared absorption spectroscopy and two-dimensional correlation spectroscopy (2D-COS) were employed to extract structural information on the mixtures from the C-D systematic stretching vibrational (νs(C-D)) region of the methyl groups in DMSO-d6. It was found that the mixing process of [CnMPyr][Tf2N] and DMSO is non-ideal and interaction complexes form between [CnMPyr][Tf2N] and DMSO-d6. They are ion cluster-DMSO-d6 complexes and ion pair-DMSO-d6 complexes. In the mixing processes, the species present in pure DMSO gradually decrease from DMSO dimer to DMSO monomer with an increase in ILs. Besides, the ion cluster-DMSO complexes gradually increase, while the ion pair-DMSO complexes decrease due to the strong electrostatic interaction between the cation and anion. In the ion cluster-DMSO complexes and ion pair-DMSO complexes, the ring hydrogen atoms of the methylene group directly attached to the nitrogen atom are the preferred interaction sites of the [CnMPyr]+ cations. All the hydrogen bonds in the identified complexes are closed-shell, electrostatically dominant and weak.

The structure and hydrogen-bond behaviours of binary systems containing ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate and methanol/ethanol

Mixing ionic liquids (ILs) with a molecular cosolvent can largely reduce the high viscosities, high polarities, and high costs of ILs. The macroscopic properties of IL-cosolvent mixtures have been studied extensively. However, some fundamental questions regarding the microscopic properties of the binary mixtures still remain to be answered. In this work, the structural and hydrogen-bond features of binary systems containing 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF₄]) and methanol/ethanol were studied by using the combination of Fourier-transform infrared spectroscopy (FTIR) and density functional theory (DFT) calculations. Excess infrared absorption spectroscopy with enhanced spectral resolution was used to analyse the original IR spectra. The alcohol tetramer/larger multimers, alcohol trimer, anion-alcohol, ion pair-alcohol, and ion cluster-alcohol complexes were identified in the excess spectra. With the increasing [BMIM][BF₄], the alcohol multimers gradually broke out from the larger multimers into smaller multimers. The hydrogen-bonded complex related with anion [BF₄]^- and alcohol gradually changes from anion-alcohol complex to ion pair-alcohol complex. The ion cluster-alcohol appears when the x(alcohol) is <0.50. The most stable optimized geometries of anion-alcohol, ion pair-alcohol, and ion cluster-alcohol were carefully analysed, and the hydrogen-bonds were identified. All of the hydrogen-bonds in these studied complexes had weak strength, closed shells and electrostatically dominant interactions.

Interfacial Forces across Ionic Liquid Solutions: Effects of Ion Concentration and Water Domains †

Using the surface force apparatus (SFA), the interaction forces between mica surfaces across ionic liquid (IL) solutions are studied. The IL solution, 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide in propylene carbonate solvent, is used at different concentrations to elucidate the ions' conformation at the interface from the analysis of short-range structural forces. A direct correlation between the ion layer thickness at the interface and the IL molar fraction in the solution is observed, suggesting conformational changes relative to the ion packing density. In addition, effects of large microscopic and macroscopic water domains at the interface are investigated. The microscopic water domains induced significant adhesion at contact because of the long-range capillary forces, which are found to depend on solvent concentration. The macroscopic water domains entirely cover the interaction area, ensuring that the long-range interfacial interactions occur entirely across the aqueous electrolyte solution with dissolved IL ions as the electrolyte. These results help elucidate the interfacial interactions in IL-charged solid interfaces with practical importance in green energy storage, catalysis, and lubrication.

Structural Forces in Mixtures of Ionic Liquids with Organic Solvents

Using molecular dynamics simulations, we study the impact of electrode charging and addition of solvent (acetonitrile, ACN) on structural forces of the BMIM PF₆ ionic liquid (IL) confined by surfaces at nanometer separations. We establish relationships between the structural forces and the microscopic structure of the confined liquid. Depending on the structural arrangements of cations and anions across the nanofilm, the load-induced squeeze-out of liquid layers occurs via one-layer or bilayer steps. The cations confined between charged plates orient with their aliphatic chain perpendicular to the surface planes and link two adjacent IL layers. These structures facilitate the squeeze-out of single layers. For both pure IL and IL-ACN mixtures, we observe a strong dependence of nanofilm structure on the surface charge density, which affects the simulated pressure-displacement curves. Addition of solvent to the IL modifies the layering in the confined film. At high electrode charges and high dilution of IL (below 10% molar fraction), the layered structure of the nanofilm is less well defined. We predict a change in the squeeze-out mechanism under pressure, from a discontinuous squeeze-out (for high IL concentrations) to an almost continuous one (for low IL concentrations). Importantly, our simulations show that charged electrodes are coated with ions even at low IL concentrations. These ion-rich layers adjacent to the charged plate surfaces are not squeezed out even under very high normal pressures of ∼5 GPa. Hence, we demonstrate the high performance of IL-solvent mixtures to protect surfaces from wear and to provide lubrication at high loads.

Solvent-Mediated Affinity of Polymer-Wrapped Single-Walled Carbon Nanotubes for Chemically Modified Surfaces

Semiconducting single-walled carbon nanotube (s-CNT) arrays are being explored for next-generation semiconductor electronics. Even with the multitude of alignment and spatially localized s-CNT deposition methods designed to control s-CNT deposition, fundamental understanding of the driving forces for s-CNT deposition is still lacking. The individual roles of the dispersant, solvent, target substrate composition, and the s-CNT itself are not completely understood because it is difficult to decouple deposition parameters. Here, we study poly[(9,9-dioctylfluorenyl-2,7-diyl)-alt-co-(6,6'-[2,2'-{bipyridine}])] (PFO-BPy)-wrapped s-CNT deposition from solution onto a chemically modified substrate. We fabricate various self-assembled monolayers (SAMs) to gain a greater understanding of substrate effects on PFO-BPy-wrapped s-CNT deposition. We observe that s-CNT deposition is dependent on both the target substrate and s-CNT dispersion solvent. To complement the experiments, molecular dynamics simulations of PFO-BPy-wrapped s-CNT deposition on two different SAMs are performed to obtain mechanistic insights into the effect of the substrate and solvent on s-CNT deposition. We find that the global free-energy minimum associated with favorable s-CNT adsorption occurs for a configuration in which the minimum of the solvent density around the s-CNT coincides with the minimum of the solvent density above a SAM-grafted surface, indicating that solvent structure near a SAM-grafted surface determines the adsorption free-energy landscape driving s-CNT deposition. Our results will help guide informative substrate design for s-CNT array fabrication in semiconductor devices.

Large Variations in the Composition of Ionic Liquid-Solvent Mixtures in Nanoscale Confinement

Mixtures of an ionic liquid with an organic solvent are widely used as electrolytes in supercapacitors where they are often confined in porous electrodes with pore widths only slightly larger than the sizes of bare ions or solvent molecules. The composition of the electrolyte inside these pores, which may depend on the pore width and choice of electrolyte, can affect supercapacitor performance but remains poorly understood. Here, we perform all-atom molecular dynamics simulations of solutions of two different ionic liquids in acetonitrile under confinement between graphene sheets forming slit pores of various widths. We observe significant oscillations in the in-pore ionic liquid mole fraction with varying pore widths. Ions are excluded from very narrow pores, while for pore widths that tightly fit a single layer of ions, we observe an in-pore ionic liquid mole fraction over three times greater than that in the bulk. At slightly larger pore widths, we observe for different ionic liquids either a nearly complete exclusion of ions from the pore or a slight depletion of ions, while ion population again increases as pore width further increases. We develop an analytical model that can qualitatively predict the in-pore ionic liquid mole fraction based on the effective molar volumes and the pore wall interaction energies of each species. Our work suggests a new avenue for tuning the ionic liquid mole fraction in nanopores with potentially significant implications for designing systems involving nanoconfined liquid electrolytes such as supercapacitors where in-pore ion population can affect charging dynamics.

The intimate relationship between the dielectric response and the decay of intermolecular correlations and surface forces in electrolytes

A general, exact theory for the decay of interactions between any particles immersed in electrolytes, including surface forces between macroscopic bodies, is derived in a self-contained, physically transparent manner. It is valid for electrolytes at any density, including ionic gases, molten salts, ionic liquids, and electrolyte solutions with molecular solvent at any concentration. The ions, the solvent and any other particles in the system can have any sizes, any shapes and arbitrary internal charge distributions. The spatial propagation of the interactions in electrolytes has several decay modes with different decay lengths that are given by the solutions, κν, ν = 1, 2,…, to a general equation for the screening parameter κ; an equation that describes the dielectric response. There can exist simultaneous decay modes with plain exponential decay and modes with damped oscillatory exponential decay, as observed experimentally and theoretically. In the limit of zero ionic density, the decay length 1/κν of the mode with the longest range approaches the Debye length 1/κD. The coupling between fluctuations in number density and charge density, described by the density-charge correlation function HNQ(r), makes all decay modes of pair correlations and interaction free energies identical to those of the screened electrostatic potential, and hence they have the same values for the screening parameters. The density-density and charge-charge correlation functions, HNN(r) and HQQ(r), also have these decay modes. For the exceptional case of charge-inversion invariant systems, HNQ(r) is identically zero for symmetry reasons and HNN(r) and HQQ(r) have, instead, decay modes with different decay lengths.

Bridging the gap between two different scaling laws for structuring of liquids under geometrical confinement

Structural forces are a phenomena obtained in liquids of one-component (e.g. for organic solvents) and two-components (colloidal dispersions), alike. So far, those two systems were discussed separately, using two different scaling laws. In this review article, an attempt is made to bridge the gap between both scaling laws by defining the scaling limit for two-component systems. Colloidal probe atomic force microscopy (CP-AFM) is used to measure structural forces in suspensions of silica nanoparticles (NPs) of three different sizes. In these two-component systems (solid NPs suspended in water), oscillatory behaviour can be obtained in the force vs. separation profiles. The wavelength λ is larger than the actual particle diameter d and rather depends on the particles' volume fraction ϕ following the inverse cubic root law λ∝ϕ^-13. It is shown that the real particle diameter d can be determined by a gedankenexperiment by extrapolating the fitted wavelength λ from the structural force measurements at a specific particle concentration to a particle volume fraction ϕ of 52% - the packing factor for simple cubic packing - using the well-known inverse cubic root scaling law. This extrapolation can be interpreted as a transition from a two-component system towards a one-component-like problem. In this case, particles are in contact and the wavelength λ is equal to the particle diameter d, λ = d as for one-component systems. The determined diameters d of the different silica nanoparticles agree well with independent measurements using transmission electron microscopy (TEM), validating the used approach. The proposed method can be extended to numerous dispersions of spherical nano-sized objects, for which structural forces can be measured.

Unexpected stability of aqueous dispersions of raspberry-like colloids

Aqueous colloidal suspensions, both man-made and natural, are part of our everyday life. The applicability of colloidal suspensions, however, is limited by the range of conditions over which they are stable. Here we report a novel type of highly monodisperse raspberry-like colloids, which are prepared in a single-step synthesis that relies on simultaneous dispersion and emulsion polymerisation. The resulting raspberry colloids behave almost like hard spheres. In aqueous solutions, such prepared raspberries show unexpected stability against aggregation over large variations of added salt concentrations without addition of stabilisers. We present simple Derjaguin–Landau–Verwey–Overbeek (DLVO) calculations performed on raspberry-like and smooth colloids showing that this stability results from our raspberries’ unique morphology, which extends our understanding of colloidal stability against salting. Further, the raspberries’ stability facilitates the formation of superspheres and thin films in which the raspberry colloids self-assemble into hexagonally close-packed photonic crystals with exquisite reproducibility.

The ability to stabilise colloidal suspensions in solution against salt-induced aggregation is critical to many industrial applications, but it remains challenging at high salt concentration. To overcome this problem, Lan et al. introduce a raspberry-like colloidal particle with controllable morphology.

Screening Lengths in Ionic Fluids

The decay of correlations in ionic fluids is a classical problem in soft matter physics that underpins applications ranging from controlling colloidal self-assembly to batteries and supercapacitors. The conventional wisdom, based on analyzing a solvent-free electrolyte model, suggests that all correlation functions between species decay with a common decay length in the asymptotic far field limit. Nonetheless, a solvent is present in many electrolyte systems. We show using an analytical theory and molecular dynamics simulations that multiple decay lengths can coexist in the asymptotic limit as well as at intermediate distances once a hard sphere solvent is considered. Our analysis provides an explanation for the recently observed discontinuous change in the structural force across a thin film of ionic liquid-solvent mixtures as the composition is varied, as well as reframes recent debates in the literature about the screening length in concentrated electrolytes.

Recent highlights in nanoscale and mesoscale friction

Friction is the oldest branch of non-equilibrium condensed matter physics and, at the same time, the least established at the fundamental level. A full understanding and control of friction is increasingly recognized to involve all relevant size and time scales. We review here some recent advances on the research focusing of nano- and mesoscale tribology phenomena. These advances are currently pursued in a multifaceted approach starting from the fundamental atomic-scale friction and mechanical control of specific single-asperity combinations, e.g., nanoclusters on layered materials, then scaling up to the meso/microscale of extended, occasionally lubricated, interfaces and driven trapped optical systems, and eventually up to the macroscale. Currently, this "hot" research field is leading to new technological advances in the area of engineering and materials science.

Experimental evaluation of additional short ranged repulsion in structural oscillation forces

The paper addresses additional short ranged repulsion in structural oscillation forces between silica surfaces across a suspension of silica nanoparticles. Fit and prediction of the structural oscillation forces usually involve an exponentially decreasing harmonic as introduced by Israelachvili [Israelachvili, Intermolecular & surface forces, Academic Press, San Diego, USA, 1985]. Recently we demonstrated, for aqueous suspensions of silica nanoparticles at various concentrations, that this fit equation is insufficient to describe the structural oscillation forces in its whole range [Schön et al., Beilstein J. Nanotechnol., 2018, 9, 1095-1107]. An additional force acting on short separations leads to the fit parameters scattering widely as well as being dependent on each other and the starting point of the fit. An additional repulsive term was introduced to solve these problems. The additional repulsive force has also been observed by others, in ionic liquids and polyelectrolyte solutions at high ionic strength. It was attributed to the diffusive double layer forces. The rise of the additional repulsion with increasing particle concentration seems to conflict with this interpretation. In this work, colloidal probe atomic force microscopy is used in aqueous suspensions of silica nanoparticles to investigate other contributing factors such as the increasing hydrodynamic drag in the normal direction to the confining surface with increasing particle concentration. A kinetic component to the structural oscillation forces is observed. Furthermore, sodium chloride is used to adjust the ionic strength of two different concentrated silica nanoparticle suspensions. For these systems the additional decay length is compared to the Debye length in the range from low to high ionic strength. A master curve of the additional decay length over Debye length at different ionic strengths, approach speed and particle concentration is produced. It affirms the link between the two and the connection between the additional force and the diffusive double layer forces. The increasing trend for the additional repulsion with increasing particle concentration reveals a synergistic effect of diffusive double layer forces and structural oscillation forces at low to medium ionic strength, which cannot be observed at high ionic strength.