On the Moisture Absorption Capability of Ionic Liquids

Due to their many attractive physicochemical properties, ionic liquids (ILs) have received extensive attention with numerous applications proposed in various fields of science and technology. Despite this, the molecular origins of many of their properties, such as the moisture absorption capability, are still not well understood. For insight into this, we systematically synthesized 24 types of ILs by the combination of the dimethyl phosphate anion with various types of alkyl group-substituted cyclic cations─imidazolium, pyrazolium, 1,2,3-triazolium, and 1,2,4-triazolium cations─and performed a detailed analysis of the dehumidification properties of these ILs and their aqueous solutions. It was found that these IL systems have a high dehumidification capability (DC). Among the monocationic ILs, the best performance was obtained with 1-cyclohexylmethyl-4-methyl-1,2,4-triazolium dimethyl phosphate, whose DC (per mol) value is 14 times higher than that of popular solid desiccants like CaCl2 and silica gel. Dicationic ILs, such as 1,1'-(propane-1,3-diyl)bis(4-methyl-1,2,4-triazolium) bis(dimethyl phosphate), showed an even better moisture absorption, with a DC (per mol) value about 20 times higher than that of CaCl2. Small- and wide-angle X-ray scattering measurements of eight types of 1,2,4-triazolium dimethyl phosphate ILs were performed and revealed that the majority of these ILs form nanostructures. Such nanostructures, which vary with the identity of the IL and the water content, fall into three main categories: bicontinuous microemulsions, hexagonal cylinders, and micelle-like structures. Water in the solutions exists primarily in polar regions in the nanostructures; these spaces function as water pockets at relatively low water concentrations. Since the structure and stability of the aggregated forms of the ILs are mainly governed by the interactions of nonpolar groups, the alkyl side chains of the cations play an important role in the DC and temperature-dependent equilibrium water vapor pressure of the IL solutions. Our experimental findings and molecular dynamics simulation results shed light on the moisture absorption mechanism of the IL aqueous solutions from a molecular perspective.

Molecular-scale understanding of diluent effects on ligand assembly for metal ion separations

Use of metal-selective ligands in solvent extraction is instrumental in extraction of critical materials and recycling, yet, diluent effects on extraction performance are not well understood. Experimental and empirical solvent parameters have been proposed to correlate with extraction performance, but are often inadequate predictors. We follow the hypothesis that the diluents' primary influence on extraction efficiency is whether or not it hinders assembly of the bulky extracting ligands into a geometry necessary for metal complexation. This behavior is readily accessible with molecular dynamics (MD), where the atomistic description of molecules can be applied to arbitrary extractant-solvent molecules and their mixtures. Several simulated quantities are considered, from both pairwise and graph theoretical analyses, and compared to experimental distribution ratio data for americium extraction by TODGA in a series of inert, non-interacting diluents. These simple properties, especially the formation of closed triplets corresponding to the 3 : 1 ligand : metal stoichiometric solvate, suggest a potential predictive power of this approach. This methodology provides a path forward to comprehensively understand and predict diluent effects in more complex systems involving different extracting ligands and multi-component diluent mixtures.

Solvation Structure of Methanol-in-Salt Electrolyte Revealed by Small-Angle X-ray Scattering and Simulations

The solvation structure of water-in-salt electrolytes was thoroughly studied, and two competing structures─anion solvated structure and anion network─were well-defined in recent publications. To further reveal the solvation structure in those highly concentrated electrolytes, particularly the influence of solvent, methanol was chosen as the solvent for this proposed study. In this work, small-angle X-ray scattering, small-angle neutron scattering, Fourier-transform infrared spectroscopy, and Raman spectroscopy were utilized to obtain the global and local structural information. With the concentration increment, the anion network formed by TFSI- became the dominant structure. Meanwhile, the hydrogen bonds among methanol were interrupted by the TFSI- anion and formed a new connection with them. Molecular dynamic simulations with two different force fields (GAFF and OPLS-AA) are tested, and GAFF agreed with synchrotron small-angle X-ray scattering/wide-angle X-ray scattering (SAXS/WAXS) results well and provided insightful information about molecular/ion scale solvation structure. This article not only deepens the understanding of the solvation structure in highly concentrated solutions, but more importantly, it provides additional strong evidence for utilizing SAXS/WAXS to validate molecular dynamics simulations.

Understanding the emergence of the boson peak in molecular glasses

A common feature of glasses is the "boson peak", observed as an excess in the heat capacity over the crystal or as an additional peak in the terahertz vibrational spectrum. The microscopic origins of this peak are not well understood; the emergence of locally ordered structures has been put forward as a possible candidate. Here, we show that depolarised Raman scattering in liquids consisting of highly symmetric molecules can be used to isolate the boson peak, allowing its detailed observation from the liquid into the glass. The boson peak in the vibrational spectrum matches the excess heat capacity. As the boson peak intensifies on cooling, wide-angle x-ray scattering shows the simultaneous appearance of a pre-peak due to molecular clusters consisting of circa 20 molecules. Atomistic molecular dynamics simulations indicate that these are caused by over-coordinated molecules. These findings represent an essential step toward our understanding of the physics of vitrification.

An Overview of Structure and Dynamics Associated with Hydrophobic Deep Eutectic Solvents and Their Applications in Extraction Processes

Recent development of novel water-immiscible green solvents known as hydrophobic deep eutectic solvents (HDESs) has opened the gates for applications requiring media where presence of water is undesirable. Ever since they were prepared, researchers have used HDESs in diverse fields such as extraction processes, CO 2 sequestration, membrane formation, and catalysis. The microstructure and dynamics associated with the species comprising HDESs guide their suitability for specific applications. For example, varying the alkyl tail length of HDES components significantly affects the dynamics of the components and thus helps in tuning the efficiency of extraction processes. The development of HDESs is still in infancy and very few theoretical studies are available in the literature that help in understanding the structure and dynamics of HDESs. This review highlights the recent work focused on the microscopic structure and dynamics of HDESs and their potential applications, particularly in extraction processes. We have also provided a glimpse of how the integration of experiments and computational techniques can help understand the mechanism of extraction processes.

Extracting free energies of counterion binding to polyelectrolytes by molecular dynamics simulations

We use all-atom molecular dynamics simulations to extract ΔG^eff, the free energy of binding of potassium ions K⁺ to the partially charged polyelectrolyte poly(acrylic acid), or PAA, in dilute regimes. Upon increasing the charge fraction of PAA, the chains adopt more extended conformations, and simultaneously, potassium ions bind more strongly (i.e., with more negative ΔG^eff) to the highly charged chains to relieve electrostatic repulsions between charged monomers along the chains. We compare the simulation results with the predictions of a model that describes potassium binding to PAA chains as a reversible reaction whose binding free energy (ΔG^eff) is adjusted from its intrinsic value (ΔG) by electrostatic correlations, captured by a random phase approximation. The bare or intrinsic binding free energy ΔG, which is an input in the model, depends on the binding species and is obtained from the radial distribution function of K⁺ around the charged monomer of a singly charged, short PAA chain in dilute solutions. We find that the model yields semi-quantitative predictions for ΔG^eff and the degree of potassium binding to PAA chains, α, as a function of PAA charge fraction without using fitting parameters.

Development of a PointNet for Detecting Morphologies of Self-Assembled Block Oligomers in Atomistic Simulations

Molecular simulations with atomistic or coarse-grained force fields are a powerful approach for understanding and predicting the self-assembly phase behavior of complex molecules. Amphiphiles, block oligomers, and block polymers can form mesophases with different ordered morphologies describing the spatial distribution of the blocks, but entirely amorphous nature for local packing and chain conformation. Screening block oligomer chemistry and architecture through molecular simulations to find promising candidates for functional materials is aided by effective and straightforward morphology identification techniques. Capturing 3-dimensional periodic structures, such as ordered network morphologies, is hampered by the requirement that the number of molecules in the simulated system and the shape of the periodic simulation box need to be commensurate with those of the resulting network phase. Common strategies for structure identification include structure factors and order parameters, but these fail to identify imperfect structures in simulations with incorrect system sizes. Building upon pioneering work by DeFever et al. [Chem. Sci. 2019, 10, 7503-7515] who implemented a PointNet (i.e., a neural network designed for computer vision applications using point clouds) to detect local structure in simulations of single-bead particles and water molecules, we present a PointNet for detection of nonlocal ordered morphologies of complex block oligomers. Our PointNet was trained using atomic coordinates from molecular dynamics simulation trajectories and synthetic point clouds for ordered network morphologies that were absent from previous simulations. In contrast to prior work on simple molecules, we observe that large point clouds with 1000 or more points are needed for the more complex block oligomers. The trained PointNet model achieves an accuracy as high as 0.99 for globally ordered morphologies formed by linear diblock, linear triblock, and 3-arm and 4-arm star-block oligomers, and it also allows for the discovery of emerging ordered patterns from nonequilibrium systems.

Heterogeneity in hydrophobic deep eutectic solvents: SAXS prepeak and local environments

The observation of the prepeak in the simulated total X-ray scattering structure function (S(q)) reveals the presence of intermediate-range structural heterogeneity in hydrophobic deep eutectic solvents.

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.

Molecular design of self-coacervation phenomena in block polyampholytes

Coacervation is a common phenomenon in natural polymers and has been applied to synthetic materials systems for coatings, adhesives, and encapsulants. Single-component coacervates are formed when block polyampholytes exhibit self-coacervation, phase separating into a dense liquid coacervate phase rich in the polyampholyte coexisting with a dilute supernatant phase, a process implicated in the liquid-liquid phase separation of intrinsically disordered proteins. Using fully fluctuating field-theoretic simulations using complex Langevin sampling and complementary molecular-dynamics simulations, we develop molecular design principles to connect the sequenced charge pattern of a polyampholyte with its self-coacervation behavior in solution. In particular, the lengthscale of charged blocks and number of connections between oppositely charged blocks are shown to have a dramatic effect on the tendency to phase separate and on the accessible chain conformations. The field and particle-based simulation results are compared with analytical predictions from the random phase approximation (RPA) and postulated scaling relationships. The qualitative trends are mostly captured by the RPA, but the approximation fails catastrophically at low concentration.

The effect of polymorphism on the structural, dynamic and dielectric properties of plastic crystal water: A molecular dynamics simulation perspective

We have employed molecular dynamics simulations based on the TIP4P/2005 water model to investigate the local structural, dynamical, and dielectric properties of the two recently reported body-centered-cubic and face-centered-cubic plastic crystal phases of water. Our results reveal significant differences in the local orientational structure and rotational dynamics of water molecules for the two polymorphs. The probability distributions of trigonal and tetrahedral order parameters exhibit a multi-modal structure, implying the existence of significant local orientational heterogeneities, particularly in the face-centered-cubic phase. The calculated hydrogen bond statistics and dynamics provide further indications of the existence of a strongly heterogeneous and rapidly interconverting local orientational structural network in both polymorphs. We have observed a hindered molecular rotation, much more pronounced in the body-centered-cubic phase, which is reflected by the decay of the fourth-order Legendre reorientational correlation functions and angular Van Hove functions. Molecular rotation, however, is additionally hindered in the high-pressure liquid compared to the plastic crystal phase. The results obtained also reveal significant differences in the dielectric properties of the polymorphs due to the different dipolar orientational correlation characterizing each phase.

The evolution of acidic and ionic aggregates in ionomers during microsecond simulations

We performed microsecond-long, atomistic molecular dynamics simulations on a series of precise poly(ethylene-co-acrylic acid) ionomers neutralized with lithium, with three different spacer lengths between acid groups on the ionomers and at two temperatures. Ionic aggregates form in these systems with a variety of shapes ranging from isolated aggregates to percolated aggregates. At the lower temperature of 423 K, the ionic aggregate morphologies do not reach a steady-state distribution over the course of the simulations. At the higher temperature of 600 K, the aggregates are sufficiently mobile that they rearrange and reach steady state after hundreds of nanoseconds. For systems that are 100% neutralized with lithium, the ions form percolated aggregates that span the simulation box in three directions, for all three spacer lengths (9, 15, and 21). In the partially neutralized systems, the morphology includes lithium ion aggregates that may also include some unneutralized acid groups, along with a coexisting population of acid group aggregates that form through hydrogen bonding. In the lithium ion aggregates, unneutralized acid groups tend to be found on the ends or sides of the aggregates.

Structure and dynamics of liquid CS2: Going from ambient to elevated pressure conditions

Molecular dynamics simulation studies were performed to investigate the structural and dynamic properties of liquid carbon disulfide (CS₂) from ambient to elevated pressure conditions. The results obtained have revealed structural changes at high pressures, which are related to the more dense packing of the molecules inside the first solvation shell. The calculated neutron and X-ray structure factors have been compared with available experimental diffraction data, also revealing the pressure effects on the short-range structure of the liquid. The pressure effects on the translational, reorientational, and residence dynamics are very strong, revealing a significant slowing down when going from ambient pressure to 1.2 GPa. The translational dynamics of the linear CS₂ molecules have been found to be more anisotropic at elevated pressures, where cage effects and librational motions are reflected on the shape of the calculated time correlation functions and their corresponding spectral densities.

Direct Comparison of Atomistic Molecular Dynamics Simulations and X-ray Scattering of Polymerized Ionic Liquids

The design of solid-state electrolytes for electrochemical applications that utilize polymerized ionic liquids (polyILs) would greatly benefit from a molecular-level understanding of structure-function relationships. We herein use atomistic molecular dynamics simulations to investigate the structural properties of a homologous series of poly(n-alkyl-vinylimidzolium bistrifluoromethylsulfonylimide) poly(nVim Tf2N) and present the first direct comparison of the structure factors obtained from X-ray scattering and simulations. Excellent agreement is found in terms of peak position and shape. The backbone-to-backbone correlation length increases at a rate of 1 Å/CH₂. The longer alkyl chains lead to the longer backbone-to-backbone separation and the larger nonpolar nanodomains. This quantitative comparison of atomistic simulations to X-ray scattering will lead to a fundamental understanding in structure and morphology of polyILs and pave a path forward toward the rational design of future polyILs for electrochemical devices.