A Refined Force Field for Molecular Simulation of Imidazolium-Based Ionic Liquids

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li+) and bis(trifluoromethanesulfonyl)imide (TFSI-) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force-torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion-ion and ion-solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1-4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Overall, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

Molecular insights into the nanoconfinement effect on the structure and dynamics of ionic liquids in carbon nanotubes

The properties and applications of ionic liquids (ILs) have been widely investigated when they are confined within nanochannels such as carbon nanotubes (CNTs). The confined ILs exhibit very different properties from their bulk state due to a nanoconfinement effect, which plays an important role in the performances of devices with ILs. In this work, we studied the effect of the charge carried by CNTs on confined ILs inside CNTs using molecular dynamics simulations. In charged CNTs, cations and anions are distributed separately along the radial directions, and the transition of orientations of the cations between parallel and vertical to CNTs occurs by changing the charge state of CNTs. The number of hydrogen bonds (HBs) formed by the confined ILs can be reduced by switching the surface charge of CNTs from positive to negative due to the contact modes between cations and anions as well as the distributions of cations in CNTs. The diffusivities along and vertical to the axial direction of CNTs were found to be non-monotonic owing to the "trade-off" effect from both ion pair interlocking and anchoring ILs on the CNT walls. Additionally, the region-dependent dynamics of ILs were also related to the intermolecular interactions in different regions of CNTs. Furthermore, the vibrational modes of ILs were obviously influenced in highly charged CNTs as determined by calculating the density of vibrational states, which demonstrated the transitions in the structure and interactions. The density distributions changed from single layer to double layers when increasing the pore size of neutral CNTs while the hydrogen bonds exhibited a non-monotonic tendency versus the pore sizes. Our results might help to understand the structure and dynamics of confined ILs as well as aid optimizing the performance of devices with ILs.

Novel barostat implementation for molecular dynamics

We propose a novel implementation of the extended-dynamics equations for isothermal-isobaric ensemble in molecular dynamics, as the Martyna-Tobias-Klein thermostat and barostat. This method is suitable for systems with constraints and the Verlet-family integrators. Instead of iterations or the Trotter-expansion-based methods, both velocities and box sizes (scaling of bond lengths) are predicted. The algorithm begins with force calculation, requiring neither quarter nor half time steps, and necessitating iterations only inside SHAKE. Several tests demonstrate that the quality is comparable to other implementations. It is found that the formula relating the extended barostat mass to the characteristic time of volume fluctuations is inaccurate for condensed systems, which has consequences for the parameter setup. Emphasis is also put on the verification of the precise isothermal-isobaric ensemble and finite-size effects.

Exploring hydrogen-bond structures in cellulose during regeneration with anti-solvent through two-dimensional correlation infrared spectroscopy

Cellulose, renowned for its excellent biocompatibility, finds extensive applications in both industrial and laboratory settings. However, few studies have specifically addressed the mechanistic evolution of hydrogen bond networks in cellulose during the dissolution and regeneration processes. In this research, the regeneration mechanism of cellulose in water and ethanol is investigated through molecular dynamics simulations. The results indicate that the ability of water molecules to disrupt hydrogen bonds between cellulose and ionic liquids is stronger than that of ethanol, which is more conducive to promoting the regeneration of cellulose. Besides, the Fourier transform infrared spectroscopy coupled with two-dimensional correlation infrared spectroscopy techniques are employed to unveil the evolution sequence of hydrogen bonds during dissolution and regeneration: ν(OH) (absorbed water) → ν(O3-H3···O5) (intrachain) → ν(O6-H6···O3') (interchain) → ν(O2-H2···O6) (intrachain) → ν(OH) (free). This study not only enhances our understanding of the intricate hydrogen bond dynamics in cellulose dissolution and regeneration but also provides a foundation for the expanded application of cellulose in diverse fields.

Liquid Madelung energy accounts for the huge potential shift in electrochemical systems

Achievement of carbon neutrality requires the development of electrochemical technologies suitable for practical energy storage and conversion. In any electrochemical system, electrode potential is the central variable that regulates the driving force of redox reactions. However, quantitative understanding of the electrolyte dependence has been limited to the classic Debye-Hückel theory that approximates the Coulombic interactions in the electrolyte under the dilute limit conditions. Therefore, accurate expression of electrode potential for practical electrochemical systems has been a holy grail of electrochemistry research for over a century. Here we show that the 'liquid Madelung potential' based on the conventional explicit treatment of solid-state Coulombic interactions enables quantitatively accurate expression of the electrode potential, with the Madelung shift obtained from molecular dynamics reproducing a hitherto-unexplained huge experimental shift for the lithium metal electrode. Thus, a long-awaited method for the description of the electrode potential in any electrochemical system is now available.

Temperature-dependent differential capacitance of an ionic liquid-graphene-based supercapacitor

One of the critical factors affecting the performance of supercapacitors is thermal management. The design of supercapacitors that operate across a broad temperature range and at high charge/discharge rates necessitates understanding the correlation of the molecular characteristics of the device (such as interfacial structure and inter-ionic and ion-electrode interactions) with its macroscopic properties. In this study, we use molecular dynamics (MD) simulations to investigate the influence of Joule heating on the structure and dynamics of the ionic liquid (IL)/graphite-based supercapacitors. The temperature-dependent electrical double layer (EDL) and differential capacitance-potential (CD-V) curves of two different ([Bmim][BF4] and [Bmim][PF6]) IL-graphene pairs were studied under various thermal gradients. For the [Bmim][BF4] system, the differential capacitance curves transition from 'U' to bell shape under an applied thermal gradient (∇T) in the range from 3.3 K nm-1 to 16.7 K nm-1. Whereas in [Bmim][PF6], we find a positive dependence of differential capacitance with ∇T with a U-shaped CD-V curve. We examine changes in the EDL structure and screening potential (ϕ(z)) as a function of ∇T and correlate them with the trends observed in the CD-V curve. The identified correlation between the interfacial charge density and differential capacitance with thermal gradient would be helpful for the molecular design of the IL-electrode interface in supercapacitors or other chemical engineering applications.

GAFF-Based Polarizable Force Field Development and Validation for Ionic Liquids

Ionic liquids (ILs) have been used in many applications, including gas separations, electrochemistry, lubrication, and catalysis. Understanding how the different properties of ILs are related to their chemical structure and composition is crucial for these applications. Experimental investigations often provide limited insights and can be tedious in exploring a range of state points. Therefore, molecular simulations have emerged as a powerful tool that not only offers a microscopic perspective but also enables rapid screening and prediction of physical properties. The accuracy of these predictions, however, depends on the quality of the intermolecular potentials (force fields) used. The widely used classical fixed charge models, such as GAFF, OPLS, and CL&P, are popular due to their simplicity and computational efficiency. However, it has been shown that the use of integer charges with these classical models leads to sluggish dynamics. The use of scaled charge models can improve the dynamics, but these mean-field approaches are unable to account for polarization effects explicitly. Several different approaches have been proposed to include polarizability in IL force fields. In this work, we follow the protocol of the CL&Pol model to develop a Drude oscillator model based on the GAFF force field (Goloviznina, K., et al. J. Chem. Theory Comput. 2019, 15, 5858). We compare the performance of the model for eight imidazolium- and pyrrolidinium-based ILs against that of other models. We find that the new model provides reasonable estimations of density, self-diffusivity, and structural properties for these ILs and suggests a relatively simple way of extending the general GAFF model to more ILs.

Invariant dynamics in a united-atom model of an ionic liquid

We study a united-atom model of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethyl)sulfonylamide to determine to what extent there exist curves in the phase diagram along which the microscopic dynamics are invariant when expressed in dimensionless, or reduced, form. The initial identification of these curves, termed isodynes, is made by noting that contours of reduced shear viscosity and reduced self-diffusion coefficient coincide to a good approximation. Choosing specifically the contours of reduced viscosity as nominal isodynes, further simulations were carried out for state points on these, and other aspects of dynamics were investigated to study their degree of invariance. These include the mean-squared displacement, shear-stress autocorrelation function, and various rotational correlation functions. These were invariant to a good approximation, with the main exception being rotations of the anion about its long axis. The dynamical features that are invariant have in common that they are aspects that would be relevant for a coarse-grained description of the system; specifically, removing the most microscopic degrees of freedom in principle leads to a simplification of the potential energy landscape, which allows for the existence of isodynes.

Structure and properties variations of regenerated cellulose fibers induced by metal ion impurity

In the ionic liquids (ILs) method for processing regenerated cellulose fiber (RCF), which is a high-performance ecologically benign product, metal ion impurities (such as Fe3+ and Cu2+) of cellulose might inevitably remain in the recycled ILs and coagulation bath. The presence of metal ions might negatively impact the properties of the manufactured RCFs and obstruct their applications, which are urgent to be made clear. For this research, the solvent for dissolving wood pulp cellulose (WPC) was 1-ethyl-3-methylimidazolium diethyl phosphate ([Emim]DEP) with various metal ion concentrations. The effect of metal ions in IL on the dissolution of cellulose was investigated by Molecular Dynamics simulations. Rheological analysis and degree of polymerization (DP) analysis were adopted to evaluate the influence on fiber spinnability of different spinning solution metal ion concentrations and various dissolving times. Further, the morphology and mechanical performances of the RCFs variation regulation were also thoroughly researched. The findings showed that the presence of metal ions in the spinning solution affected the DP, crystallinity, and orientation factor of RCFs, which will influence their stress more sensitively than the strain. These findings can serve as a practical guide for the commercial manufacture of emerging fiber.

Effect of Ion Mass on Dynamic Correlations in Ionic Liquids

Ionic liquids (ILs) are a class of liquid salts with distinct properties such as high ionic conductivity, low volatility, and a broad electrochemical window, making them appealing for use in energy storage applications. The ion-ion correlations are some of the key factors that play a critical role in the ionic conductivity of ILs. In this work, we present the study of the impact of ion mass on ion-ion correlations in ILs, applying a combination of broadband dielectric spectroscopy measurements and molecular dynamics simulations. We examined three ILs with the same cation but different anions to consider three different cases of cation-anion masses: M+ > M-, M+ ≈ M-, and M+ < M-. We applied the momentum conservation approach to estimate the contribution of distinct ion-ion correlations from experimental data and obtained good agreement with direct calculations of distinct ion-ion correlations from molecular dynamics simulations. Our findings reveal that relative ion mass has a strong effect on the distinct ion-ion correlations, leading to swapping of the relative amplitude of distinct cation-cation and anion-anion correlations.

Influence of Anion Species on Liquid-Liquid Phase Separation in [EMIm+][X-]/Benzene Mixtures

When a molar excess of benzene is mixed with an ionic liquid (IL), liquid-liquid phase separation may appear with a pure liquid phase almost composed of only benzene molecules separated from a liquid clathrate phase with benzene molecules dissolved in the IL. Our previous study (J. Phys. Chem. B, 124, 7929, 2020) on long-chain IL/benzene systems has concluded that benzene molecules, as planar nonpolar ones, majorly dissolve in the IL nonpolar domains consisting of cationic alkyl side chains. Nevertheless, the above mechanism is inadequate for explaining the experimental observations that benzene can also dissolve in IL systems with very short alkyl side chains. In this study, by molecular dynamics simulation of the [EMIm+][X-]/benzene mixtures with X- being Cl-, NO3-, PF6-, or Tf2N-, we still observe liquid-liquid phase separation of the pure benzene phase from a liquid clathrate (mixed IL/benzene) phase where benzene molecules are almost equally distributed near imidazolium rings through π-stacking or near alkyl side chains. The anion species strongly influences the solubility of benzene and the ratio of the two liquid phases via the alteration of anionic charge density, which tunes the strength of the electrostatic interaction among ions and thus the probability of benzene molecules interacting with both imidazolium rings and alkyl side chains: a larger anionic charge density corresponds to a lower solubility of benzene.

Highly efficient zinc electrode prepared by electro-deposition in a salt-induced pre-phase separation region solution

Efficient electrode design is crucial for the electrochemical reduction of CO2 to produce valuable chemicals. The solution used for the preparation of electrodes can affect their overall properties, which in turn determine the reaction efficiency. In this work, we report that transition metal salts could induce the change of two-phase ionic liquid/ethanol mixture into miscible one phase. Pre-phase separation region near the phase boundary of the ternary system was observed. Zinc nanoparticles were electro-deposited along the fibres of carbon paper (CP) substrate uniformly in the salt-induced pre-phase separation region solution. The as-prepared Zn(1)/CP electrode exhibits super-wettability to the electrolyte, rendering very high catalytic performance for CO2 electro-reduction, and the Faradaic efficiency towards CO is 97.6% with a current density of 340 mA cm-2, which is the best result to date in an H-type cell.

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid-liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid-liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. A detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

Spectroscopic and Chemical Properties of Ionic Liquids: Computational Study

A brief account is given of highlights of our computational efforts - often in collaboration with experimental groups - to understand spectroscopic and chemical properties of ionic liquids (ILs). Molecular dynamics, including their inhomogeneous character, responsible for key spectral features observed in dielectric absorption, infra-red (IR) and fluorescence correlation spectroscopy (FCS) measurements are elucidated. Mechanisms of chemical processes involving imidazolium-based ILs are illustrated for CO2 capture and related reactions, transesterification of cellulose, and Au nanocluster-catalyzed Suzuki cross-coupling reaction with attention paid to differing roles of IL ions. A comparison with experiments is also made.

Theory and Simulations of Ionic Liquids in Nanoconfinement

Room-temperature ionic liquids (RTILs) have exciting properties such as nonvolatility, large electrochemical windows, and remarkable variety, drawing much interest in energy storage, gating, electrocatalysis, tunable lubrication, and other applications. Confined RTILs appear in various situations, for instance, in pores of nanostructured electrodes of supercapacitors and batteries, as such electrodes increase the contact area with RTILs and enhance the total capacitance and stored energy, between crossed cylinders in surface force balance experiments, between a tip and a sample in atomic force microscopy, and between sliding surfaces in tribology experiments, where RTILs act as lubricants. The properties and functioning of RTILs in confinement, especially nanoconfinement, result in fascinating structural and dynamic phenomena, including layering, overscreening and crowding, nanoscale capillary freezing, quantized and electrotunable friction, and superionic state. This review offers a comprehensive analysis of the fundamental physical phenomena controlling the properties of such systems and the current state-of-the-art theoretical and simulation approaches developed for their description. We discuss these approaches sequentially by increasing atomistic complexity, paying particular attention to new physical phenomena emerging in nanoscale confinement. This review covers theoretical models, most of which are based on mapping the problems on pertinent statistical mechanics models with exact analytical solutions, allowing systematic analysis and new physical insights to develop more easily. We also describe a classical density functional theory, which offers a reliable and computationally inexpensive tool to account for some microscopic details and correlations that simplified models often fail to consider. Molecular simulations play a vital role in studying confined ionic liquids, enabling deep microscopic insights otherwise unavailable to researchers. We describe the basics of various simulation approaches and discuss their challenges and applicability to specific problems, focusing on RTIL structure in cylindrical and slit confinement and how it relates to friction and capacitive and dynamic properties of confined ions.

A force field for bio-polymers in ionic liquids (BILFF) - part 2: cellulose in [EMIm][OAc]/water mixtures

We present the extension of our force field BILFF (Bio-Polymers in Ionic Liquids Force Field) to the bio-polymer cellulose. We already published BILFF parameters for mixtures of ionic liquid 1-ethyl-3-methylimidazolium acetate ([EMIm][OAc]) with water. Our all-atom force field focuses on a quantitative reproduction of the hydrogen bonds in the complex mixture of cellulose, [EMIm]⁺, [OAc]^- and water when compared to reference ab initio molecular dynamics (AIMD) simulations. To enhance the sampling, 50 individual AIMD simulations starting from different initial configurations were performed for cellulose in solvent instead of one long simulation, and the resulting averages were used for force field optimization. All cellulose force field parameters were iteratively adjusted starting from the literature force field of W. Damm et al. We were able to obtain a very good agreement with respect to both the microstructure of the reference AIMD simulations and experimental results such as the system density (even at higher temperatures) and the crystal structure. Our new force field allows performing very long simulations of large systems containing cellulose solvated in (aqueous) [EMIm][OAc] with almost ab initio accuracy.

Structure and Dynamics of Hydrofluorocarbon/Ionic Liquid Mixtures: An Experimental and Molecular Dynamics Study

The physical properties of four ionic liquids (ILs), including 1-n-butyl-3-methylimidazolium tetrafluoroborate ([C₄C₁im][BF₄]), 1-n-butyl-3-methylimidazolium hexafluorophosphate ([C₄C₁im][PF₆]), 1-n-butyl-3-methylimidazolium thiocyanate ([C₄C₁im][SCN]), and 1-n-hexyl-3-methylimidazolium chloride ([C₆C₁im][Cl]), and their mixtures with hydrofluorocarbon (HFC) gases HFC-32 (CH₂F₂), HFC-125 (CHF₂CF₃), and HFC-410A, a 50/50 wt % mixture of HFC-32 and HFC-125, were studied using molecular dynamics (MD) simulation. Experiments were conducted to measure the density, self-diffusivity, and shear viscosity of HFC/[C₄C₁im][BF₄] system. Extensive analyses were carried out to understand the effect of IL structure on various properties of the HFC/IL mixtures. Density, diffusivity, and viscosity of the pure ILs were calculated and compared with experimental values. The good agreement between computed and experimental results suggests that the applied force fields are reliable. The calculated center of mass (COM) radial distribution functions (RDFs), partial RDFs, spatial distribution functions (SDFs), and coordination numbers (CNs) provide a sense of how the distribution of HFC changes in the liquid mixtures with IL structure. Detailed analysis reveals that selectivity toward HFC-32 and HFC-125 depends on both cation and anion. The molecular insight provided in the current work will help the design of optimal ILs for the separation of azeotropic HFC mixtures.

Ultratough Hydrogen-Bond-Bridged Phosphorene Films

The rapid development of flexible electronic devices, especially based on 2D materials, has triggered the demand for high-strength materials. Mono- or few-layer phosphorene with excellent electronic properties has attracted extensive attention. However, phosphorene is affected by its low Young's modulus when applied to flexible electronic devices. Here, a strategy via ion intercalation to significantly improve the mechanical properties of black phosphorus to generate hydrogen-bond-bridged phosphorene films with Young's modulus as high as 316 GPa is reported. This value is several times larger than the theoretical values of 166 GPa in the zigzag direction, 44 GPa in the armchair direction, and the averaged Young's modulus among all directions of 94 GPa. The impact of intercalation on mechanical properties is also explored. Experimental nanoindentation results obtained by atomic force microscopy indicate that the relationship between the ratio of intercalated ions to phosphorus atoms and the corresponding Young's modulus satisfies the formula E = e a e - [ ln ( x ) + b ] 2 c ( 0 < x ≤ 1.80 ) \[E{\bm{ = }}{e^{a{e^{\frac{{{\bm{ - }}{{[\ln (x){\bm{ + }}b]}^2}}}{c}}}}}(0{\bm{ &lt; }}x{\bm{ \le }}1.80)\] . Furthermore, a flexible NO₂ gas sensor device based on this ultratough material presents excellent performance, even after 10 000 bending cycles. The results provide new insight into the potential for practical applications of black phosphorus devices.

Infrared and Terahertz Spectroscopic Investigation of Imidazolium, Pyridinium, and Tetraalkylammonium Tetrafluoroborate Ionic Liquids

We performed terahertz time-domain spectroscopy and infrared spectroscopy of imidazolium-based, pyridinium-based, and tetraalkylammonium-based tetrafluoroborate ionic liquids to study their characteristic intermolecular and intramolecular vibrational modes to clarify interactions between various cations and the tetrafluoroborate anion. It was found that the central frequency of the intermolecular vibrational band for these ionic liquids has a relatively high frequency, ranging from 90 to 100 cm^-1. In the 900-1150 cm^-1 range, the intramolecular vibrational absorption band of the 3-fold degenerate mode of tetrafluoroborate anions in the ionic liquids was observed. Although the tetrafluoroborate anion is attributable to one of the weakly coordinated anions, the spectroscopic splitting behavior of the 3-fold degenerate mode differs depending on the cation species. It was revealed that the degenerate mode is very sensitive to local interactions between the tetrafluoroborate anion and each cation.

Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

We have monitored the impacts of an increment in the alkyl chain length of the imidazolium-based tetrafluoroborate ionic liquids on the local deuteroxyl probe modes of interest. For this study, we have taken 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIm][BF₄], 1-butyl-3-methylimidazolium tetrafluoroborate [BMIm][BF₄], 1-octyl-3-methylimidazolium tetrafluoroborate [OMIm][BF₄], and 1-decyl-3-methylimidazolium tetrafluoroborate [DMIm][BF₄] ionic liquid solutions with 5% HOD in H₂O as the vibrational reporter of the associated ultrafast system dynamics. Classical molecular dynamics (MD) simulations were employed to determine molecular structure and dynamic properties, while the spectral profiles were derived by applying the wavelet analysis of classical trajectories. Spatial distribution functions reveal the heterogeneity within the molecular structures of the ionic liquids (ILs) with varying alkyl chain lengths. The intense position of the spectral peak, the frequency corresponding to the shoulder peak, and the spectral linewidth of the O-D stretch distribution are not influenced by the increment in the cationic chain length. In addition, the ionic liquid (IL) [BMIm][BF₄] exhibits a notable trend; the dynamic timescales are longer than the other studied systems. Therefore, we have performed the Voronoi decomposition analysis of the ionic and the polar-apolar domains, symmetrically increasing the length of alkyl chains on the IL cations. Domain analysis reveals structural microheterogeneity; the anions form discrete domains, and the ionic liquid constituting cations form continuous domains irrespective of the alkyl chain length on the imidazolium cations. Therefore, this computational ultrafast spectroscopy study aids in forming a molecular-level picture of the ionic liquid cations and anions in the liquid phase, providing a detailed interpretation of the spectral properties of the probe stretching vibrations.

Insight into the dual effect of water on lignin dissolution in ionic liquids

The dual regulation of water on lignin in ionic liquids was studied at the molecular level by molecular dynamics simulation. The simulation results show that a small amount of water will destroy the ion association in ionic liquids, that is, it will produce more free anions and cations. The free ions around lignin are conducive to the dissolution of lignin. On the contrary, excess water will seriously solvate anions and cations. By changing the number of lignin clusters, it is more intuitive to observe that the dissolution of lignin in ILs containing a small amount of water is stronger than that in pure IL, however, the dissolution ability of lignin is reduced after adding a large amount of water in ILs. It is concluded that with the increase of water content, water changes from co-solvent to anti-solvent in the dissolution process. This study provides ideas for the design of IL-water system for economic pretreatment of biomass.

Theoretical Mechanism on the Cellulose Regeneration from a Cellulose/EmimOAc Mixture in Anti-Solvents

The experiments on cellulose dissolution/regeneration have made some achievements to some extent, but the mechanism of cellulose regeneration in ionic liquids (ILs) and anti-solvent mixtures remains elusive. In this work, the cellulose regeneration mechanism in different anti-solvents, and at different temperatures and concentrations, has been studied with molecular dynamics (MD) simulations. The IL considered is 1-ethyl-3-methylimidazolium acetate (EmimOAc). In addition, to investigate the microcosmic effects of ILs and anti-solvents, EmimOAc-nH₂O (n = 0–6) clusters have been optimized by Density Functional Theory (DFT) calculations. It can be found that water is beneficial to the regeneration of cellulose due to its strong polarity. The interactions between ILs and cellulose will become strong with the increase in temperature. The H-bonds of cellulose chains would increase with the rising concentrations of anti-solvents. The interaction energies between cellulose and the anions of ILs are stronger than that of cations. Furthermore, the anti-solvents possess a strong affinity for ILs, cation–anion pairs are dissociated to form H-bonds with anti-solvents, and the H-bonds between cellulose and ILs are destroyed to promote cellulose regeneration.

Insights into the Structure and Dynamics of Imidazolium Ionic Liquid and Tetraethylene Glycol Dimethyl Ether Cosolvent Mixtures: A Molecular Dynamics Approach

In this work, the effect of molecular cosolvents tetraethylene glycol dimethyl ether (TEGDME) on the structure and versatile nature of mixtures of these compounds with imidazolium-based ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF₆]) is analyzed and discussed at a molecular level by means of all-atom molecular dynamics (MD) simulations. In the whole concentration range of the binary mixtures, the structures and properties evolution was studied by means of systematic molecular dynamics simulations of the fraction of hydrogen bonds, the radial and spatial distribution functions for the various molecular ions and molecular species in the system, together with the snapshots visualization of equilibrated simulation boxes with a color-coding scheme and the rotational dynamics of coumarin 153 (C153) in the binary mixtures. The goal of the work is to provide a molecular-level understanding of significant improvement of ionic conductivity and self-diffusion with the presence of TEGDME as a cosolvent, which causes an enhancement to the ion translational motion and fluidity in the [bmim][PF₆] ionic liquids (ILs). Under a mixture concentration change, the microstructure changes of [bmim][PF₆] with the TEGDME molar fraction (X_TEG) above 0.50 show a slight difference from that of neat [bmim][PF₆] IL and concentrated [bmim][PF₆]/TEGDME mixture in terms of the radial and spatial distribution functions. The relative diffusivities of solvent molecules to cations as a function of concentration were found to depend on the solvent but not on the anion. A TEGDME increase is found to be advantageous to the dissipation of the polar regions as well as the nonpolar regions in the [bmim][PF₆] ionic liquids. These conclusions are consistent with the experimental results, which verified that the unique, complex, and versatile nature of [bmim][PF₆]/TEGDME mixture can be correctly modeled and discussed at a molecular level using MD simulation data.

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM+][OAc-]

Reactive molecular dynamics simulations enable a detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the empirical valence bond framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NNs). The physics-based terms ensure the correct asymptotic behavior of electrostatic, polarization, and dispersion interactions and are compatible with existing solvent force fields. NNs are utilized for a versatile description of short-range orbital interactions within the transition state region and accurate rendering of vibrational motion of the reacting complex. We demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium cation with acetate to form 1-ethyl-3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ∼1 kJ mol^-1 mean absolute error accuracy within the transition state region for the gas-phase complex. To characterize the solvent modulation of the reaction profile, we compute potentials of mean force for the gas-phase reaction as well as the reaction within a four-ion cluster and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes the formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

Effect of ionic liquids in carbon nanotube bundles on CO2, H2S, and N2 separation from CH4: A computational study

Carbon nanotube (CNT) bundles are being explored as a support structure for four ionic liquids (ILs) in gas separation. Grand canonical Monte Carlo simulations were performed to investigate the CO₂/CH₄, H₂S/CH₄, and N₂/CH₄ separation performance in CNT bundles and CNT-supported ILs (CNT-ILs) as a function of pressure and IL loading. The results show that by adding ILs to the CNT bundles, the gas separation performance can be significantly increased. Increasing the number of IL molecules in the composites increases the separation performance. Such a phenomenon is more evident for the CO₂/CH₄ mixture in comparison to H₂S/CH₄ and N₂/CH₄. Calculations of isosteric heat of adsorption and selectivities in gas mixtures as a function of pressure show promising gas separation performance for CNT-ILs. Due to the excellent mechanical properties of CNTs, it has been shown that this structure may be used as a strong mechanical support for structures containing ILs with excellent CO₂/CH₄ separation performance.

Solar Desalination Using Thermally Responsive Ionic Liquids Regenerated with a Photonic Heater

Growing global water demand has brought desalination technologies to the forefront for freshwater production from nontraditional water sources. Among these, forward osmosis (FO) is a promising two-step desalination process (draw dilution and regeneration), but it is often overlooked due to the energy requirements associated with draw regeneration. To address this limiting factor, we demonstrate FO desalination using thermally responsive ionic liquids (ILs) that are regenerated using a renewable energy input, that is, solar heat. To efficiently harness sunlight, a simple photonic heater converts incoming irradiation into infrared wavelengths that are directly absorbed by IL-water mixtures, thereby inducing phase separation to yield clean water. This approach is markedly different as it uses radiative heating, a noncontact mode of heat transfer that couples to chemical functional groups within the IL for rapid energy transfer without a heat exchanger or secondary fluid. Overall, a solar-thermal separation efficiency of 50% is achieved under unconcentrated sunlight, which can be increased to 69% with the thermal design. Successful desalination of produced water from oil wells in Southern California highlights the potential of solar-powered IL-FO for energy-efficient and low-cost desalination of complex brines for beneficial water reuse.

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.

A quantum chemical model for a series of self-assembled nanocages: the origin of stability behind the coordination-driven formation of transition metal complexes up to [M12L24]24

Herein, we present a systematic computational model to study the electronic states and free energies of a self-assembled multi-metal complex series. By combining the previously developed model Hamiltonian approach for transition-metal complexes and the generalized Born model, the thermodynamics, optimized geometries, and electronic states of the [Pd12L24]24+ nanocage are revealed, together with [PdnLm]2n+ complex series. The effective model Hamiltonian is a theoretical method to obtain the d-electron wavefunction and potential energy including interaction energy between the transition-metal and ligands. In the present improvement, the electronic state on each transition-metal center is focused as a building unit and solved under the whole electronic field of the assembling system. We realize a reliable and systematic treatment of multi-transition-metal complexes with different sizes and charges. Consequently, our model could reproduce binding energies of the [PdnLm]2n+ complex series quantitatively as compared to density functional theory (DFT). Regarding free energy, we revealed that the assembling solute becomes unstable due to the electrostatic interaction, and effects of the solvent and counter anions mainly compensated it. Optimized geometries were also analysed. The local square-planar coordination structures around the palladium centres were characterized in the complex series. The relationships between the entire symmetrical geometries and the local coordination structures are also discussed. Finally, electronic structures of the [Pd12L24]24+ nanocage were well characterized as a single-determinant, where only dx2-y2 is unoccupied due to the ligand-field effect. We also found that the solvent polarized the electronic states of the Pd ions, whereas the counter anion suppressed the polarization. The present method realizes size-independent reliable and rapid computations, and therefore can be expected to further application studies on self-assembly dynamics.

Infrared Spectroscopy in the Middle Frequency Range for Various Imidazolium Ionic Liquids-Common Spectroscopic Characteristics of Vibrational Modes with In-Plane +C(2)-H and +C(4,5)-H Bending Motions and Peak Splitting Behavior Due to Local Symmetry Breaking of Vibrational Modes of the Tetrafluoroborate Anion

Various alkyl-methylimidazolium ionic liquids (ILs) were inspected using infrared spectroscopy in the middle frequency range. In the 1050-1200 cm^-1 range, there is a skeletal vibrational mode accompanied with a large in-plane ⁺C(2)-H bending motion and ⁺C(4)-H and ⁺C(5)-H motions, and in the 1500-1650 cm^-1 range, there are two skeletal vibrational modes with in-plane ⁺C(4,5)-H bending motions. Interestingly, in both ranges, we found that skeletal vibrational modes with a large in-plane ⁺C(2)-H bending motion and in-plane ⁺C(4,5)-H bending motions are insensitive to increases in the basicity of anions or the strengthening of hydrogen bond-type interactions, and the behaviors are completely different from those in the ⁺C-H stretching vibrational modes in the 3000-3200 cm^-1 range and the skeletal vibrational modes with large out-of-plane ⁺C-H motions in the 700-950 cm^-1 range. Furthermore, in alkyl-methylimidazolium tetrafluoroborate [C _n mim⁺][BF₄ ^-] ILs, we found that absorption due to the (threefold) degenerate vibrational mode of [BF₄ ^-] was observed as a broad absorption band with three splitting peaks in the 900-1150 cm^-1 range as a result of local symmetry breaking due to the cation-anion interactions.

A force field for bio-polymers in ionic liquids (BILFF) - part 1: [EMIm][OAc]/water mixtures

We present BILFF, a novel force field for bio-polymers in ionic liquids. In the first part of our study, we introduce optimized force field parameters for mixtures of the ionic liquid (IL) 1-ethyl-3-methylimidazolium acetate ([EMIm][OAc]) with water. This imidazolium-based IL is of particular practical importance as it can dissolve significant amounts of cellulose even at room temperature. An understanding of this dissolution process via molecular dynamics simulations requires a quantitative description of the microscopic structure and the strong hydrogen bonds with a method able of simulating at least several dozen nanoseconds, which is the main aim of our novel force field. To reach this goal, we optimize the force field parameters to reproduce radial, spatial, and combined distribution functions, hydrogen bond lifetimes, diffusion coefficients, and several other quantities from reference ab initio molecular dynamics (AIMD) simulations. Non-trivial effects such as dispersion interactions between the side chains and π-π stacking of the cations are reproduced very well. We further validate the force field by comparison to experimental data such as thermal expansion coefficients, bulk modulus, and density at different temperatures, which yields good agreement and correct trends. No other force field with optimized parameters for mixtures of [EMIm][OAc] and water has been presented in the literature yet. Optimized force field parameters for cellulose and other ILs will be published in upcoming articles.

Molecular simulation of osmometry in aqueous solutions of the BMIMCl ionic liquid: a potential route to force field parameterization of liquid mixtures

Despite widespread development and use of ionic liquids (ILs) in both academic and industrial research, computational force fields (FFs) for most of those are not available for a precise description of inter-species interactions in aqueous environments. In the scope of this study, by means of molecular simulations, the osmotic coefficient of an aqueous solution of an IL is calculated and used as a basis to reparameterize popular IL-FFs existing in the literature. We first calculate the osmotic coefficients (at 298.15 K and 1 atm pressure) of aqueous solutions of 1-butyl-3-methylimidazolium chloride (BMIMCl), a generic IL, popularly used in biomass processing and the subsequent conversion to value-added intermediates. The performance of two popular atomic, nonpolarizable FFs developed for BMIMCl, one by Lopes, Pádua, and coworkers (FF-LP) and the other by Sambasivarao, Acevedo, and coworkers (FF-SA), when mixed with the SPC/E water model, is tested with respect to their ability to reproduce the experimental osmotic coefficient data. Interestingly, the osmotic coefficient is found to be increasing with a gradual increase in IL molality within the concentration range of our investigation, which is contrary to the experimental trend reported in the literature for the same IL-water mixture. Henceforth, necessary corrections to the nonbonded ion-ion and ion-water interactions are made to match the experimental osmotic coefficient. To further assess the reliability of the new FF, we extensively explore the thermodynamic (density, isothermal compressibility, and thermal expansion coefficient), dynamic (diffusivity and viscosity), and association/dissociation properties (rationalized with the help of radial distribution functions) with both the original and reparameterized FF for a wider range of concentrations up to a molality of 18.50 mol kg-1. The calculated quantities are compared against experimental data wherever available. The modified FF parameters exhibit significant improvements in terms of its ability to match experimental solution properties, such as density, viscosity, association/dissociation, etc. We report that excessive dissociation of BMIMCl in water is responsible for the shortcomings observed in the original FFs and improved prediction of physicochemical properties could be achieved using the modified FFs.

Computational NMR Study of Ion Pairing of 1-Decyl-3-methyl-imidazolium Chloride in Molecular Solvents

The 1H NMR spectra of 10-5 mole fraction solutions of 1-decyl-3-methyl-imidazolium chloride ionic liquid in water, acetonitrile, and dichloromethane have been measured. The chemical shift of the proton at position 2 in the imidazolium ring of 1-decyl-3-methyl-imidazolium (H2) is rather different for all three samples, reflecting the shifting equilibrium between the contact pairs and free fully solvated ions. Classical molecular dynamics simulations of the 1-decyl-3-methyl-imidazolium chloride contact ion pair as well as of free ions in water, acetonitrile, and dichloromethane have been conducted, and the quantum mechanics/molecular mechanics methods have been applied to predict NMR chemical shifts for the H2 proton. The chemical shift of the H2 proton was found to be primarily modulated by hydrogen bonding with the chloride anion, while the influence of the solvents-though differing in polarity and capabilities for hydrogen bonding-is less important. By comparing experimental and computational results, we deduce that complete disruption of the ionic liquid into free ions takes place in an aqueous solution. Around 23% of contact ion pairs were found to persist in acetonitrile. Ion-pair breaking into free ions was predicted not to occur in dichloromethane.

Pressure in Molecular Simulations with Scaled Charges. 1. Ionic Systems

Charge scaling, rationalized as MDEC (molecular dynamics in electronic continuum) or ECC (electronic continuum correction), has become a widely used simple approach to how to avoid self-consistent induced dipoles yet approximately take into account the effects of electronic polarizability. It has been assumed that the continuum permittivity does not depend on density; in turn, pressure is calculated by standard formulas. In this work, we elaborate a complementary approximation of density-independent molecular polarizability and derive formulas for pressure corrections within the MDEC framework; real behavior lies between these two extremes. The pressure corrections for test ionic systems are huge and negative, leading to sizable densities in constant-pressure MDEC simulations. A comparison of MDEC results with equivalent polarizable systems gives a good pressure match for a crystal but very low MDEC pressures for ionic liquids. These results witness about the importance of a correct density dependence not only of continuum permittivity in MDEC simulations but also of polarizability in polarizable simulations.