Ion Pairing Controls Physical Properties of Ionic Liquid-Solvent Mixtures

Development of Heteroatomic Constant Potential Method with Application to MXene-Based Supercapacitors

We describe a method for modeling constant-potential charges in heteroatomic electrodes, keeping pace with the increasing complexity of electrode composition and nanostructure in electrochemical research. The proposed "heteroatomic constant potential method" (HCPM) uses minimal added parameters to handle differing electronegativities and chemical hardnesses of different elements, which we fit to density functional theory (DFT) partial charge predictions in this paper by using derivative-free optimization. To demonstrate the model, we performed molecular dynamics simulations using both HCPM and conventional constant potential method (CPM) for MXene electrodes with Li-TFSI/AN (lithium bis(trifluoromethane sulfonyl)imide/acetonitrile)-based solvent-in-salt electrolytes. Although the two methods show similar accumulated charge storage on the electrodes, the results indicated that HCPM provides a more reliable depiction of electrode atom charge distribution and charge response compared with CPM, accompanied by increased cationic attraction to the MXene surface. These results highlight the influence of elemental composition on electrode performance, and the flexibility of our HCPM opens up new avenues for studying the performance of diverse heteroatomic electrodes including other types of MXenes, two-dimensional materials, metal-organic frameworks (MOFs), and doped carbonaceous electrodes.

Structural evolution of water-in-propylene carbonate mixtures revealed by polarized Raman spectroscopy and molecular dynamics

The liquid structure of systems wherein water is limited in concentration or through geometry is of great interest in various fields such as biology, materials science, and electrochemistry. Here, we present a combined polarized Raman and molecular dynamics investigation of the structural changes that occur as water is added incrementally to propylene carbonate (PC), a polar, aprotic solvent that is important in lithium-ion batteries. Polarized Raman spectra of PC solutions were collected for water mole fractions 0.003 ≤ χwater ≤ 0.296, which encompasses the solubility range of water in PC. The novel approach taken herein provides additional hydrogen bond and solvation characterization of this system that has not been achievable in previous studies. Analysis of the polarized carbonyl Raman band in conjunction with simulations demonstrated that the bulk structure of the solvent remained unperturbed upon the addition of water. Experimental spectra in the O-H stretching region were decomposed through Gaussian fitting into sub-bands and comparison to studies of dilute HOD in D2O. With the aid of simulations, we identified these different bands as water arrangements having different degrees of hydrogen bonding. The observed water structure within PC indicates that water tends to self-aggregate, forming a hydrogen bond network that is distinctly different from the bulk and dependent on concentration. For example, at moderate concentrations, the most likely aggregate structures are chains of water molecules, each with two hydrogen bonds.

Anion and ether group influence in protic guanidinium ionic liquids

Ionic liquids are attractive liquid materials for many advanced applications. For targeted design, in-depth knowledge about their structure-property-relations is urgently needed. We prepared a set of novel protic ionic liquids (PILs) with a guanidinium cation with either an ether or alkyl side chain and different anions. While being a promising cation class, the available data is insufficient to guide design. We measured thermal and transport properties, nuclear magnetic resonance (NMR) spectra as well as liquid and crystalline structures supported by ab initio computations and were able to obtain a detailed insight into the influence of the anion and the ether substitution on the physical and spectroscopic properties. For the PILs, hydrogen bonding is the main interaction between cation and anion and the H-bond strength is inversely related to the proton affinity of the constituting acid and correlated to the increase of ¹H and ¹⁵N chemical shifts. Using anions from acids with lower proton affinity leads to proton localization on the cation as evident from NMR spectra and self-diffusion coefficients. In contrast, proton exchange was evident in ionic liquids with triflate and trifluoroacetate anions. Using imide-type anions and ether side groups decreases glass transitions as well as fragility, and accelerated dynamics significantly. In case of the ether guanidinium ionic liquids, the conformation of the side chain adopts a curled structure as the result of dispersion interactions, while the alkyl chains prefer a linear arrangement.

Structure-Dynamics Interrelation Governing Charge Transport in Cosolvated Acetonitrile/LiTFSI Solutions

Concentrated ionic solutions present a potential improvement for liquid electrolytes. However, their conductivity is limited by high viscosities, which can be attenuated via cosolvation. This study employs a series of experiments and molecular dynamics simulations to investigate how different cosolvents influence the local structure and charge transport in concentrated lithium bis(trifluoromethane-sulfonyl)imide (LiTFSI)/acetonitrile solutions. Regardless of whether the cosolvent's dielectric constant is low (for toluene and dichloromethane), moderate (acetone), or high (methanol and water), they preserve the structural and dynamical features of the cosolvent-free precursor. However, the dissimilar effects of each case must be individually interpreted. Toluene and dichloromethane reduce the conductivity by narrowing the distribution of Li⁺-TFSI^- interactions and increasing the activation energies for ionic motions. Methanol and water broaden the distributions of Li⁺-TFSI^- interactions, replace acetonitrile in the Li⁺ solvation, and favor short-range Li⁺-Li⁺ interactions. Still, these cosolvents strongly interact with TFSI^-, leading to conductivities lower than that predicted by the Nernst-Einstein relation. Finally, acetone preserves the ion-ion interactions from the cosolvent-free solution but forms large solvation complexes by joining acetonitrile in the Li⁺ solvation. We demonstrate that cosolvation affects conductivity beyond simply changing viscosity and provide fairly unexplored molecular-scale perspectives regarding structure/transport phenomena relation in concentrated ionic solutions.

Blistering failure of elastic coatings with applications to corrosion resistance

A variety of polymeric surfaces, such as anti-corrosion coatings and polymer-modified asphalts, are prone to blistering when exposed to moisture and air. As water and oxygen diffuse through the material, dissolved species are produced, which generate osmotic pressure that deforms and debonds the coating. These mechanisms are experimentally well-supported; however, comprehensive macroscopic models capable of predicting the formation osmotic blisters, without extensive data-fitting, is scant. Here, we develop a general mathematical theory of blistering and apply it to the failure of anti-corrosion coatings on carbon steel. The model is able to predict the irreversible, nonlinear blister growth dynamics, which eventually reaches a stable state, ruptures, or undergoes runaway delamination, depending on the mechanical and adhesion properties of the coating. For runaway delamination, the theory predicts a critical delamination length, beyond which unstable corrosion-driven growth occurs. The model is able to fit multiple sets of blister growth data with no fitting parameters. Corrosion experiments are also performed to observe undercoat rusting on carbon steel, which yielded trends comparable with model predictions. The theory is used to define three dimensionless numbers which can be used for engineering design of elastic coatings capable of resisting visible deformation, rupture, and delamination.

Binary mixtures of ionic liquids: Ideal, non-ideal, or quasi-ideal?

The mixing of ILs provides an opportunity for fine tuning the physiochemical properties of ILs for various applications. However, a suitable mixture having desired properties can only be designed when the physiochemical properties of the mixtures of ILs along with their spectroscopic properties are well understood. With an aim to achieve this objective, three different mixtures with a common anion, namely, [C₂C₁im][C₄C₁im][NTf₂], [C₃C₁pyr][C₄C₁pyr][NTf₂], and [C₃C₁im][C₃C₁pyr][NTf₂], have been investigated in the current study. Investigations have been carried out at the macroscopic level by observing the thermophysical properties, such as molar volume and thermal expansion coefficient, and at the microscopic level with time-resolved fluorescence measurements and the pulse field gradient nuclear magnetic resonance (NMR) technique. The results obtained from the thermophysical study have indicated that excess molar volume for imidazolium-based IL-IL mixtures may be linked to the free volume created by the alkyl chain of the imidazolium cation whereas for the mixture of pyrrolidinium ILs, lowering of density can give rise to free volume. Analysis of time-resolved fluorescence anisotropy data has provided clear evidence in favor of the presence of free volume in the binary mixture of ILs. NMR studies have also supported the fluorescence anisotropy data. The outcome of the present investigation reveals that the mixtures show appreciable deviation from ideal behavior and the deviation from the ideal behavior is caused due to the generation of free volume in the resultant mixture, describing these IL mixtures as quasi-ideal rather than ideal or non-ideal.

Tailoring amphotericin B as an ionic liquid: an upfront strategy to potentiate the biological activity of antifungal drugs

Aspergillus species are the primary cause of invasive aspergillosis, which afflicts hundreds of thousands of patients yearly, with high mortality rates. Amphotericin B is considered the gold standard in antifungal drug therapy, due to its broad-spectrum activity and rarely reported resistance. However, low solubility and permeability, as well as considerable toxicity, challenge its administration. Lipid formulations of amphotericin B have been used to promote its slow release and diminish toxicity, but these are expensive and adverse health effects of their prolonged use have been reported. In the past decades, great interest emerged on converting biologically active molecules into an ionic liquid form to overcome limitations such as low solubility or polymorphisms. In this study, we evaluated the biological activity of novel ionic liquid formulations where the cholinium, cetylpyridinium or trihexyltetradecylphosphonium cations were combined with an anionic form of amphotericin B. We observed that two formulations increased the antifungal activity of the drug, while maintaining its mode of action. Molecular dynamics simulations showed that higher biological activity was due to increased interaction of the ionic liquid with the fungal membrane ergosterol compared with amphotericin B alone. Increased cytotoxicity could also be observed, probably due to greater interaction of the cation with cholesterol, the main sterol in animal cells. Importantly, one formulation also displayed antibacterial activity (dual functionality), likely preserved from the cation. Collectively, the data set ground for the guided development of ionic liquid formulations that could improve the administration, efficacy and safety of antifungal drugs or even the exploitation of their dual functionality.

Pyr1,xTFSI Ionic Liquids (x = 1–8): A Computational Chemistry Study

Pyrrolidinium-based (Pyr) ionic liquids are a very wide family of molecular species. Pyrrolidinium cations are electrochemically stable in a large potential interval and their molecular size hinders their transport properties. The corresponding ionic liquids with trifluoromethyl sulphonyl imide anions are excellent solvents for lithium/sodium salts and have been demonstrated as electrolytes in aprotic batteries with enhanced safety standards. In this study, the analysis of the physicochemical properties of a homologous series of pyrrolidinium-based ionic liquids with general formula Pyr1,xTFSI (x = 1–8) have been tackled by first principles calculations based on the density functional theory. The molecular structures of isolated ions and ion pairs have been predicted by electronic structure calculations at B3LYP level of theory in vacuum or in simulated solvents. Thermodynamic properties have been calculated to evaluate the ion pairs dissociation and oxidation/reduction stability. This is the first systematic computational analysis of this series of molecules with a specific focus on the impact of the length of the alkyl chain on the pyrrolidinium cation on the overall physicochemical properties of the ion pairs.

Diffusivity and Structure of Room Temperature Ionic Liquid in Various Organic Solvents

Room-temperature ionic liquids (RTILs) hold promise for applications in electric double layer capacitors (EDLCs), owing to a much wider potential window, lower vapor pressure, and better thermal and chemical stabilities compared to conventional aqueous and organic electrolytes. However, because the low diffusivity of ions in neat RTILs negates the EDLCs' advantage of high power density, the ionic liquids are often used in mixture with organic solvents. In this study, we measured the diffusivity of cations and anions in RTIL, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl) ([BMIM⁺][TFSI^-]), mixed with 10 organic solvents, by using the pulsed-field gradient NMR method. The ion diffusivity was found to follow that of neat solvents and in most studied solvents showed an excellent agreement with the predicted values reported in the recent molecular dynamics (MD) study [Thompson, M. W.; J. Phys. Chem. B 2019, 123, 1340-1347]. In two solvents consisting of long-chain molecules, however, the MD simulations predictions slightly underestimated the ionic diffusivities. The degree of ion dissociation was also estimated for each solvent by comparing the ionic conductivity with the molar conductivity derived from the diffusion measurements. The degree of ion dissociation and the hydrodynamic radius of ions suggest that the ions are coordinated by ∼1 solvent molecule. The scarcity of solvent-ion interactions explains the fact that the diffusivity of ions in the mixture significantly depends on the viscosity of the solvent.

Towards Molecular Simulations that are Transparent, Reproducible, Usable By Others, and Extensible (TRUE)

Systems composed of soft matter (e.g., liquids, polymers, foams, gels, colloids, and most biological materials) are ubiquitous in science and engineering, but molecular simulations of such systems pose particular computational challenges, requiring time and/or ensemble-averaged data to be collected over long simulation trajectories for property evaluation. Performing a molecular simulation of a soft matter system involves multiple steps, which have traditionally been performed by researchers in a "bespoke" fashion, resulting in many published soft matter simulations not being reproducible based on the information provided in the publications. To address the issue of reproducibility and to provide tools for computational screening, we have been developing the open-source Molecular Simulation and Design Framework (MoSDeF) software suite. In this paper, we propose a set of principles to create Transparent, Reproducible, Usable by others, and Extensible (TRUE) molecular simulations. MoSDeF facilitates the publication and dissemination of TRUE simulations by automating many of the critical steps in molecular simulation, thus enhancing their reproducibility. We provide several examples of TRUE molecular simulations: All of the steps involved in creating, running and extracting properties from the simulations are distributed on open-source platforms (within MoSDeF and on GitHub), thus meeting the definition of TRUE simulations.