Conformational and Dynamic Properties of Poly(ethylene oxide) in an Ionic Liquid: Development and Implementation of a First-Principles Force Field

Molecular dynamics simulation of salt diffusion in constituting phosphazene-based polymer electrolyte

A growing demand to visualize polymer models in liquid poses a computational challenge in molecular dynamics (MD) simulation, as this requires emerging models under suitable force fields (FFs) to capture the underlying molecular behaviour accurately. In our present study, we have employed TIP3P potential on water and all atomistic optimized potentials for liquid simulations FFs to study the liquid electrolyte behavior of phosphazene-based polymer by considering its potential use in lithium-ion polymer batteries. We have explored the polymer's local structure, chain packing, wettability, and hydrophobic tendencies against the silicon surface using a combination of a pseudocontinuum model in MD simulation, and surface-sensitive sum frequency generation (SFG) vibrational spectroscopy. The finding yields invaluable insights into the molecular architecture of phosphazene. This approach identifies the importance of hydrophobic interactions with air and hydrophilic units with water molecules in understanding the behavior and properties of phosphazene-based polymers at interfaces, contributing to its advancements in materials science. The MD study uniquely captures traces of the polymer-ion linkage, which is observed to become more pronounced with the increase in polymer weight fraction. The theoretical observation of this linkage's influence on lithium-ion diffusion motion offers valuable insights into the fundamental physics governing the behavior of atoms and molecules within phosphazene-based polymer electrolytes in aqueous environments. Further these predictions are corroborated in the molecular-level depiction at the air-aqueous interface, as evidenced from the OH-oscillator strength variation measured by the SFG spectroscopy.The fundamental findings from this study open new avenues for utilizing MD simulation as a versatile methodology to gain profound insights into intermolecular interactions of polymer. It could be useful in the application of biomedical and energy-related research, such as polymer lithium-ion batteries, fuel cells, and organic solar cells.

Topology Automated Force-Field Interactions (TAFFI): A Framework for Developing Transferable Force Fields

Force-field development has undergone a revolution in the past decade with the proliferation of quantum chemistry based parametrizations and the introduction of machine learning approximations of the atomistic potential energy surface. Nevertheless, transferable force fields with broad coverage of organic chemical space remain necessary for applications in materials and chemical discovery where throughput, consistency, and computational cost are paramount. Here, we introduce a force-field development framework called Topology Automated Force-Field Interactions (TAFFI) for developing transferable force fields of varying complexity against an extensible database of quantum chemistry calculations. TAFFI formalizes the concept of atom typing and makes it the basis for generating systematic training data that maintains a one-to-one correspondence with force-field terms. This feature makes TAFFI arbitrarily extensible to new chemistries while maintaining internal consistency and transferability. As a demonstration of TAFFI, we have developed a fixed-charge force-field, TAFFI-gen, from scratch that includes coverage for common organic functional groups that is comparable to established transferable force fields. The performance of TAFFI-gen was benchmarked against OPLS and GAFF for reproducing several experimental properties of 87 organic liquids. The consistent performance of these force fields, despite their distinct origins, validates the TAFFI framework while also providing evidence of the representability limitations of fixed-charge force fields.

Deep Eutectic Solvents: Molecular Simulations with a First-Principles Polarizable Force Field

The unique properties of deep eutectic solvents make them useful in a variety of applications. In this work we develop a first-principles force field for reline, which is composed of choline chloride and urea in the molar ratio 1:2. We start with the symmetry adapted perturbation theory (SAPT) protocol and then make adjustments to better reproduce the structure and dynamics of the liquid when compared to first-principles molecular dynamics (FPMD) simulations. The resulting force field is in good agreement with experiments in addition to being consistent with the FPMD simulations. The simulations show that primitive molecular clusters are preferentially formed with choline-chloride ionic pairs bound with a hydrogen bond in the hydroxyl group and that urea molecules coordinate the chloride mainly via the trans-H chelating hydrogen bonds. Incorporating polarizability qualitatively influences the radial distributions and lifetimes of hydrogen bonds and affects long-range structural order and dynamics. The polarizable force field predicts a diffusion constant about an order of magnitude larger than the nonpolarizable force field and is therefore less computationally intensive. We hope this study paves the way for studying complex hydrogen-bonding liquids from a first-principles approach.

Phase Behavior of Poly(ethylene oxide) in Room Temperature Ionic Liquids: A Molecular Simulation and Deep Neural Network Study

The phase behavior of polymers in room temperature ionic liquids is a topic of considerable interest. In this work we study the phase diagram of poly(ethylene oxide) in four imidazolium ionic liquids (ILs) using molecular simulation. We develop united atom models for 1-butyl-2,3-dimethylimidazolium ([BMMIM]), 1-ethyl-2,3-dimethylimidazolium ([EMMIM]), and 1-ethyl-3-methylimidazolium ([EMIM]) in an analogous fashion to previously developed models for 1-butyl-3-methylimidazolium ([BMIM]) and tetrafluoroborate ([BF₄]) using symmetry-adapted perturbation theory. At high temperatures we obtain the coexistence concentrations using an interface method where the polymer and IL are simulated in a large elongated box, and an interface between coexisting phases is formed. At lower temperatures we use a deep neural network (DNN) method. The input descriptors for the DNN are the cohesive energy of mixing, the volume change of mixing, and the coordination numbers between cation and polymer, all of which are obtained from simulations of mixed systems at a series of temperatures. The DNN is trained by using the phase-separated systems at high temperatures and a mixed phase at low temperatures. The method predicts a lower critical solution temperature which decreases as the alkyl chain length on the cation is decreased, consistent with experiment. The simulations show that methylation of the cation has little effect on the phase diagram. This is in contrast to what is seen in experiments but could be because the polymer chains in the simulations are too short. At low temperatures the chains display two conformational motifs, namely a crown ether conformation and a ring conformation, each of which can wrap the chain around a single cation. This provides the entropic penalty for mixing and a reason for demixing as the temperature is raised. Such conformations might not be possible for longer chains. The combination of data-driven techniques and molecular simulation shows promise in the study of the phase behavior and physical properties of complex fluids.

Solvation and diffusion of poly(vinyl alcohol) chains in a hydrated inorganic ionic liquid

While the behavior of polyelectrolyte chains in aqueous salt solutions has been extensively studied, little is known about polar polymer chains in solvents with extremely high concentrations of inorganic ions, such as those found in ionic liquids (ILs). Here, we report on expansion, solvation and diffusion of poly(vinyl alcohol), PVA, chains in dilute solutions of a hydrated inorganic IL phase change material (PCM), lithium nitrate trihydrate (LNH). This solvent has an extremely high concentration of inorganic ions (≈18 M) with a low concentration of water molecules largely forming solvation shells of Li⁺ and NO₃^- ions, as shown using ATR-FTIR spectroscopy. Diffusion and hydrodynamic size of PVA chains of different molecular weights in this unusual solvent were studied using fluorescence correlation spectroscopy (FCS). A higher scaling exponent obtained from the molecular weight dependences of the diffusion coefficients of PVA chains as well as a lower overlap concentration (c*) of PVA in LNH solutions as measured by FCS suggest an expansion of the polymer coils in this solvent. We argue that enhanced solubility of PVA in LNH solutions is likely a result of increased rigidification of polymer chains due to the binding of solvated Li⁺ ions, which is demonstrated using ⁷Li NMR spectroscopy. We believe that an understanding of solvation and ion-binding capability can offer crucial insight into designing polymer-based shape stabilization matrices for inorganic PCMs.

Proper Thermal Equilibration of Simulations with Drude Polarizable Models: Temperature-Grouped Dual-Nosé-Hoover Thermostat

An explicit treatment of electronic polarization is critically important to accurate simulations of highly charged or interfacial systems. Compared to the iterative self-consistent field (SCF) scheme, extended Lagrangian approaches are computationally more efficient for simulations that employ a polarizable force field. However, an appropriate thermostat must be chosen to minimize heat flow and ensure an equipartition of kinetic energy among all unconstrained system degrees of freedom. Here we investigate the effects of different thermostats on the simulation of condensed phase systems with the Drude polarizable force field using several examples that include water, NaCl/water, acetone, and an ionic liquid (IL) BMIM⁺/BF₄^-. We show that conventional dual-temperature thermostat schemes often suffer from violations of equipartitioning and adiabatic electronic state, leading to considerable errors in both static and dynamic properties. Heat flow from the real degrees of freedom to the Drude degrees of freedom leads to a steady temperature gradient and puts the system at an incorrect effective temperature. Systems with high-frequency internal degrees of freedom such as planar improper dihedrals or C-H bond stretches are most vulnerable; this issue has been largely overlooked in the literature because of the primary focus on simulations of rigid water molecules. We present a new temperature-grouped dual-Nosé-Hoover thermostat, where the molecular center of mass translations are assigned to a temperature group separated from the rest degrees of freedom. We demonstrate that this scheme predicts correct static and dynamic properties for all the systems tested here, regardless of the thermostat coupling strength. This new thermostat has been implemented into the GPU-accelerated OpenMM simulation package and maintains a significant speedup relative to the SCF scheme.

Rheology of Concentrated Polymer/Ionic Liquid Solutions: An Anomalous Plasticizing Effect and a Universality in Nonlinear Shear Rheology

An anomalous plasticizing effect was observed in polymer/ionic liquid (IL) solutions by applying broad range of rheological techniques. Poly(ethylene oxide)(PEO)/IL solutions exhibit stronger dynamic temperature dependence than pure PEO, which is in conflict with the knowledge that lower-Tg solvent increases the fractional free volume. For poly(methy methacrylate)(PMMA)/IL solutions, the subtle anomaly was detected from the fact that the effective glass transition temperature Tg,eff of PMMA in IL is higher than the prediction of the self-concentration model, while in conventional polymer solutions, Tg,eff follows the original Fox equation. Observations in both solutions reveal retarded segmental dynamics, consistent with a recent simulation result (Macromolecules, 2018, 51, 5336) that polymer chains wrap the IL cations by hydrogen bonding interactions and the segmental unwrapping delays their relaxation. Start-up shear and nonlinear stress relaxation tests of polymer/IL solutions follow a universal nonlinear rheological behavior as polymer melts and solutions, indicating that the segment-cation interaction is not strong enough to influence the nonlinear chain orientation and stretch. The present work may arouse the further theoretical, experimental, and simulation interests in interpreting the effect of complex polymer-IL interaction on the dynamics of polymer/IL solutions.

Influence of molecular weight on ion-transport properties of polymeric ionic liquids

We report the results of atomistic molecular dynamics simulations on polymerized 1-butyl-3-vinylimidazolium-hexafluorophosphate ionic liquids, studying the influence of the polymer molecular weight on the ion mobilities and the mechanisms underlying ion transport, including ion-association dynamics, ion hopping, and ion-polymer coordinations. With an increase in polymer molecular weight, the diffusivity of the hexafluorophosphate (PF₆^-) counterion decreases and plateaus above seven repeat units. The diffusivity is seen to correlate well with the ion-association structural relaxation time for pure ionic liquids, but becomes more correlated with ion-association lifetimes for larger molecular weight polymers. By analyzing the diffusivity of ions based on coordination structure, we unearth a transport mechanism in which the PF₆^- moves by "climbing the ladder" while associated with four polymeric cations from two different polymers.