Mechanisms of Ion Transport in Block Copolymeric Polymerized Ionic Liquids

GCMe: Efficient Implementation of the Gaussian Core Model with Smeared Electrostatic Interactions for Molecular Dynamics Simulations of Soft Matter Systems

In recent years, molecular dynamics (MD) simulations have emerged as an essential tool for understanding the structure, dynamics, and phase behavior of charged soft matter systems. To explore phenomena across greater length and time scales in MD simulations, molecules are often coarse-grained for better computational performance. However, commonly used force fields represent particles as hard-core interaction centers with point charges, which often overemphasizes the packing effect and short-range electrostatics, especially in systems with bulky deformable organic molecules and systems with strong coarse-graining. This underscores the need for an efficient soft-core model to physically capture the effective interactions between coarse-grained particles. To this end, we implement a soft-core model uniting the Gaussian core model with smeared electrostatic interactions that is phenomenologically equivalent to recent theoretical models. We first parametrize it generically using water as the model solvent. Then, we benchmark its performance in the OpenMM toolkit for different boundary conditions to highlight a computational speedup of up to 34 × compared to commonly used force fields and existing implementations. Finally, we demonstrate its utility by investigating how boundary polarizability affects the adsorption behavior of a polyelectrolyte solution on perfectly conducting and nonmetal boundaries.

Polymerized and Colloidal Ionic Liquids─Syntheses and Applications

The breadth and importance of polymerized ionic liquids (PILs) are steadily expanding, and this review updates advances and trends in syntheses, properties, and applications over the past five to six years. We begin with an historical overview of the genesis and growth of the PIL field as a subset of materials science. The genesis of ionic liquids (ILs) over nano to meso length-scales exhibiting 0D, 1D, 2D, and 3D topologies defines colloidal ionic liquids, CILs, which compose a subclass of PILs and provide a synthetic bridge between IL monomers (ILMs) and micro to macro-scale PIL materials. The second focus of this review addresses design and syntheses of ILMs and their polymerization reactions to yield PILs and PIL-based materials. A burgeoning diversity of ILMs reflects increasing use of nonimidazolium nuclei and an expanding use of step-growth chemistries in synthesizing PIL materials. Radical chain polymerization remains a primary method of making PILs and reflects an increasing use of controlled polymerization methods. Step-growth chemistries used in creating some CILs utilize extensive cross-linking. This cross-linking is enabled by incorporating reactive functionalities in CILs and PILs, and some of these CILs and PILs may be viewed as exotic cross-linking agents. The third part of this update focuses upon some advances in key properties, including molecular weight, thermal properties, rheology, ion transport, self-healing, and stimuli-responsiveness. Glass transitions, critical solution temperatures, and liquidity are key thermal properties that tie to PIL rheology and viscoelasticity. These properties in turn modulate mechanical properties and ion transport, which are foundational in increasing applications of PILs. Cross-linking in gelation and ionogels and reversible step-growth chemistries are essential for self-healing PILs. Stimuli-responsiveness distinguishes PILs from many other classes of polymers, and it emphasizes the importance of segmentally controlling and tuning solvation in CILs and PILs. The fourth part of this review addresses development of applications, and the diverse scope of such applications supports the increasing importance of PILs in materials science. Adhesion applications are supported by ionogel properties, especially cross-linking and solvation tunable interactions with adjacent phases. Antimicrobial and antifouling applications are consequences of the cationic nature of PILs. Similarly, emulsion and dispersion applications rely on tunable solvation of functional groups and on how such groups interact with continuous phases and substrates. Catalysis is another significant application, and this is an historical tie between ILs and PILs. This component also provides a connection to diverse and porous carbon phases templated by PILs that are catalysts or serve as supports for catalysts. Devices, including sensors and actuators, also rely on solvation tuning and stimuli-responsiveness that include photo and electrochemical stimuli. We conclude our view of applications with 3D printing. The largest components of these applications are energy related and include developments for supercapacitors, batteries, fuel cells, and solar cells. We conclude with our vision of how PIL development will evolve over the next decade.

Nano sized gallium oxide surface features for enhanced antimicrobial and osteo-integrative responses

Gallium oxide has known beneficial osteo-integrative properties. This may have importance for improving the osteointegration of orthopedic implants. At high concentrations gallium is cytotoxic. Therefore, integration of gallium into implant devices must be carefully controlled to limit its concentration and release. A strategy based on surface doping of gallium although challenging seems an appropriate approach to limit dose amounts to minimize cytotoxicity and maximize osteointegration benefits. In this work we develop a novel form of patterned surface doping via a block copolymer-based surface chemistry that enables very low gallium content but enhanced osteointegration as proven by comprehensive bioassays. Polystyrene-b-poly 4vinyl pyridine (PS-b-P4VP) BCP (block copolymer) films were produced on surfaces. Selective infiltration of the BCP pattern with a gallium salt precursor solution and subsequent UV-ozone treatment produced a surface pattern of gallium oxide nanodots as evidenced by atomic force and scanning electron microscopy. A comprehensive study of the bioactivity was carried out, including antimicrobial and sterility testing, gallium ion release kinetics and the interaction with human marrow mesenchymal stomal cells and mononuclear cells. Comparing the data from osteogenesis media assay tests with osteoclastogenesis tests demonstrated the potential for the gallium oxide nanodot doping to improve osteointegration properties of a surface.

Mesoscale Simulation Based on the Dynamic Mean-Field Density Functional Method on Block-Copolymeric Ionomers for Polymer Electrolyte Membranes

Block copolymers generally have peculiar morphological characteristics, such as strong phase separation. They have been actively applied to polymer electrolyte membranes for proton exchange membrane fuel cells (PEMFCs) to obtain well-defined hydrophilic regions and water channels as a proton pathway. Although molecular simulation tools are advantageous to investigate the mechanism of water channel formation based on the chemical structure and property relationships, classical molecular dynamics simulation has limitations regarding the model size and time scale, and these issues need to be addressed. In this study, we investigated the morphology of sulfonated block copolymers synthesized for PEM applications using a mesoscale simulation based on the dynamic mean-field density functional method, widely applied to investigate macroscopic systems such as polymer blends, micelles, and multi-block/grafting copolymers. Despite the similar solubility parameters of the monomers in our block-copolymer models, very different morphologies in our 3D mesoscale models were obtained. The model with sulfonated monomers, in which the number of sulfonic acid groups is twice that of the other model, showed better phase separation and water channel formation, despite the short length of its hydrophilic block. In conclusion, this unexpected behavior indicates that the role of water molecules is important in making PEM mesoscale models well-equilibrated in the mesoscale simulation, which results in the strong phase separation between hydrophilic and hydrophobic regions and the ensuing well-defined water channel. PEM synthesis supports the conclusion that using the sulfonated monomers with a high sulfonation degree (32.5 mS/cm) will be more effective than using the long hydrophilic block with a low sulfonation degree (25.2 mS/cm).

Predictive Molecular Models for Charged Materials Systems: From Energy Materials to Biomacromolecules

Electrostatic interactions play a dominant role in charged materials systems. Understanding the complex correlation between macroscopic properties with microscopic structures is of critical importance to develop rational design strategies for advanced materials. But the complexity of this challenging task is augmented by interfaces present in the charged materials systems, such as electrode-electrolyte interfaces or biological membranes. Over the last decades, predictive molecular simulations that are founded in fundamental physics and optimized for charged interfacial systems have proven their value in providing molecular understanding of physicochemical properties and functional mechanisms for diverse materials. Novel design strategies utilizing predictive models have been suggested as promising route for the rational design of materials with tailored properties. Here, an overview of recent advances in the understanding of charged interfacial systems aided by predictive molecular simulations is presented. Focusing on three types of charged interfaces found in energy materials and biomacromolecules, how the molecular models characterize ion structure, charge transport, morphology relation to the environment, and the thermodynamics/kinetics of molecular binding at the interfaces is discussed. The critical analysis brings two prominent field of energy materials and biological science under common perspective, to stimulate crossover in both research field that have been largely separated.

The Transport Properties of Semi-Crystalline Polyetherimide BPDA-P3 in Amorphous and Ordered States: Computer Simulations

The effect of polymer chain ordering on the transport properties of the polymer membrane was examined for the semi-crystalline heterocyclic polyetherimide (PEI) BPDA-P3 based on 3,3',4,4'-biphenyltetracarboxylic dianhydride (BPDA) and diamine 1,4-bis [4-(4-aminophenoxy)phenoxy]benzene (P3). All-atom Molecular Dynamics (MD) simulations were used to investigate the gas diffusion process carried through the pores of a free volume several nanometers in size. The long-term (~30 μs) MD simulations of BPDA-P3 were performed at T = 600 K, close to the experimental value of the melting temperature (Tm ≈ 577 K). It was found during the simulations that the transition of the PEI from an amorphous state to an ordered one occurred. We determined a decrease in solubility for both the gases examined (CO2 and CH4), caused by the redistribution of free volume elements occurring during the structural ordering of the polymer chains in glassy state (Tg ≈ 487 K). By analyzing the diffusion coefficients in the ordered state, the presence of gas diffusion anisotropy was found. However, the averaged values of the diffusion coefficients did not differ from each other in the amorphous and ordered states. Thus, permeability in the observed system is primarily determined by gas solubility, rather than by gas diffusion.

Perspective: Morphology and ion transport in ion-containing polymers from multiscale modeling and simulations

Ion-containing polymers are soft materials composed of polymeric chains and mobile ions. Over the past several decades they have been the focus of considerable research and development for their use as the electrolyte in energy conversion and storage devices. Recent and significant results obtained from multiscale simulations and modeling for proton exchange membranes (PEMs), anion exchange membranes (AEMs), and polymerized ionic liquids (polyILs) are reviewed. The interplay of morphology and ion transport is emphasized. We discuss the influences of polymer architecture, tethered ionic groups, rigidity of the backbone, solvents, and additives on both morphology and ion transport in terms of specific interactions. Novel design strategies are highlighted including precisely controlling molecular conformations to design highly ordered morphologies; tuning the solvation structure of hydronium or hydroxide ions in hydrated ion exchange membranes; turning negative ion-ion correlations to positive correlations to improve ionic conductivity in polyILs; and balancing the strength of noncovalent interactions. The design of single-ion conductors, well-defined supramolecular architectures with enhanced one-dimensional ion transport, and the understanding of the hierarchy of the specific interactions continue as challenges but promising goals for future research.

Coarse-Grained Modeling of Ion-Containing Polymers

Ion-containing polymers have continued to be an important research focus for several decades due to their use as an electrolyte in energy storage and conversion devices. Elucidation of connections between the mesoscopic structure and multiscale dynamics of the ions and solvent remains incompletely understood. Coarse-grained modeling provides an efficient approach for exploring the structural and dynamical properties of these soft materials. The unique physicochemical properties of such polymers are of broad interest. In this review, we summarize the current development and understanding of the structure-property relationship of ion-containing polymers and provide insights into the design of such materials determined from coarse-grained modeling and simulations accompanying significant advances in experimental strategies. We specifically concentrate on three types of ion-containing polymers: proton exchange membranes (PEMs), anion exchange membranes (AEMs), and polymerized ionic liquids (polyILs). We posit that insight into the similarities and differences in these materials will lead to guidance in the rational design of high-performance novel materials with improved properties for various power source technologies.

Water Dissolved in a Variety of Polymers Studied by Molecular Dynamics Simulation and a Theory of Solutions

The performance of a polymer medium as a separation membrane is determined by the dissolution free energy ΔG and diffusion coefficient D of the permeant. In this work, ΔG and D of water are investigated with all-atom molecular dynamics simulation in a wide variety of polymer species in the amorphous state. The computed ΔG is shown to agree well with the experimental value for linear homopolymers, and the degrees of polymerization of the homopolymers do not affect ΔG when they are beyond ∼10. The copolymers of ethylene-vinylidene difluoride, ethylene-vinyl acetate, and ethylene-acrylamide are then examined by changing the repeating patterns of the constituent monomers in both the periodic and graft forms. It is found that ΔG is determined primarily by the overall compositions of the monomers and is not affected by the copolymerization topology (periodic or graft). The hydrophobicity of the copolymer is enhanced, furthermore, when the hydrophobicity and hydrophilicity of the ethylene and non-ethylene parts are well contrasted and those parts are fragmented along the polymer chain. According to the computed D, the diffusivity of water tends to be larger when the (co)polymer is more hydrophobic and ΔG is more positive. D is actually seen to vary by orders of magnitude with the polymer structures, while the effect of the polymer species on the water permeation is stronger for ΔG than for D.

Quantitative Evidence of Mobile Ion Hopping in Polymerized Ionic Liquids

Atomistic molecular dynamics simulations were performed, and an extensive set of analyses were undertaken to understand the ion transport mechanism in the polymerized ionic liquid poly(C₂VIm)Tf₂N. The ion hopping events were investigated at different time scales. Ion hopping was examined by monitoring the instantaneous cation-anion association and dissociation. Ion diffusion was subsequently evaluated with correlation functions and the calculation of relaxation times at different time scales. Dynamical heterogeneity in the mobility of the ions was observed with only a small portion of the anions classified as fast mobile ions. The mobile ions were characterized as the ones traveling farther than a certain distance during a characteristic period, which was much longer than the time scale of the instant ion pair dissociation. Effective hopping of the mobile ions contributed to the diffusivity which was dominated by interchain hopping and generally facilitated with five associating cations from two different polymer chains. Mobile anions had relatively fewer associating cations from more associating chains than immobile anions. The stringlike cooperative motion was observed in the mobile anions. The string length was determined to decrease with increasing temperature. These findings provided an in-depth understanding of the ion transport in polymerized ionic liquids and important information for the rational design of novel materials.

Modulating Microphase Separation of Lamellae-Forming Diblock Copolymers via Ionic Junctions

We present a molecular dynamics simulation study investigating the phase behavior of lamellae-forming diblock copolymers with a single ionic junction on the backbone. Our results show qualitative agreement with experimental findings regarding enhanced microphase separation with the introduction of an ionic junction at the conjunction point, while further revealing nonmonotonic changes in domain spacing and order-disorder transition as a function of the electrostatic interaction strength. This highlights the dominant roles of entropic and binding effects of counterions under weak and strong ionic correlations, respectively. The location of the ionic junction is found to effectively modulate the charge distribution and chain conformation in the ordered domains; its presence in the middle of a block promotes folding of the block, leading to a smaller domain size. These findings demonstrate the interplay of ionic coupling with steric hindrance and chain end effects, which enhances our understanding of the delicate control over the microphase domain features.

Ion Mobilities, Transference Numbers, and Inverse Haven Ratios of Polymeric Ionic Liquids

We probe the ion mobilities, transference numbers, and inverse Haven ratio of ionic liquids and polymerized ionic liquids as a function of their molecular weight using a combination of atomistic equilibrium and nonequilibrium molecular dynamics simulations. In contrast to expectations, we demonstrate that the inverse Haven ratio increases with increasing degree of polymerization (N) and then decreases at larger N. For a fixed center of mass reference frame, we demonstrate that such results arise as a consequence of the strong cation-cation correlated motions, which exceed (in magnitude) the self-diffusivity of cations. Together, our findings challenge the premise underlying the pursuit of pure polymeric ionic liquids as high transference number, single-ion conducting electrolytes.