Transferable coarse-grained model for perfluorosulfonic acid polymer membranes

Morphology and Transport Study of Acid-Base Blend Proton Exchange Membranes by Molecular Simulations: Case of Chitosan/Nafion

Blending a basic polymer (e.g., chitosan) with Nafion can modify some membrane properties in direct methanol fuel cell applications, e.g., controlling methanol crossover, by regulating the morphology of hydrophilic channels. Unraveling the mechanisms by which the channel morphology is modified is essential to formulate design strategies for acid-base blend membrane development. Thus, we use molecular simulations to analyze the morphological features of a blend membrane (at 75/25 chitosan/Nafion wt %), i.e., (i) water/polymer phase organizations, (ii) number and size of water clusters, and (iii) quantitative morphological measures of hydrophilic channels, and compare them to the pure Nafion in a wide range of water contents. It is found that the affinity of water to different hydrophilic groups in the blend membrane can result in more distorted and dispersed hydrophilic phase and fewer bulk water-like features compared to pure Nafion. Also, the width of the hydrophilic network bottleneck, i.e., pore limiting diameter (PLD), is found to be almost five times smaller for the blend membrane compared to Nafion at their maximum water contents. Moreover, by changing the chitosan/Nafion weight ratio from 75/25 to 0/100, we show that as Nafion content increases, all channel morphological characteristics alter monotonically except PLD. This is mainly due to the strong acid-base interactions between Nafion and chitosan, which hinder the monotonic growth of PLD. Interestingly, water and methanol diffusion coefficients are strongly correlated with PLD, suggesting that PLD can be used as a single parameter for tailoring the blending ratio for achieving the desired diffusion properties of acid-base membranes.

Integrating Machine Learning in the Coarse-Grained Molecular Simulation of Polymers

Machine learning (ML) is having an increasing impact on the physical sciences, engineering, and technology and its integration into molecular simulation frameworks holds great potential to expand their scope of applicability to complex materials and facilitate fundamental knowledge and reliable property predictions, contributing to the development of efficient materials design routes. The application of ML in materials informatics in general, and polymer informatics in particular, has led to interesting results, however great untapped potential lies in the integration of ML techniques into the multiscale molecular simulation methods for the study of macromolecular systems, specifically in the context of Coarse Grained (CG) simulations. In this Perspective, we aim at presenting the pioneering recent research efforts in this direction and discussing how these new ML-based techniques can contribute to critical aspects of the development of multiscale molecular simulation methods for bulk complex chemical systems, especially polymers. Prerequisites for the implementation of such ML-integrated methods and open challenges that need to be met toward the development of general systematic ML-based coarse graining schemes for polymers are discussed.

Combined Separator Based on a Porous Ion-Exchange Membrane for Zinc-Halide Batteries

In this work, we report on a comparative analysis of the bromine permeability for three separator groups under the operating conditions of a non-flow zinc-bromine battery. A new method for the synthesis of porous heterogeneous membranes based on a cation-exchange resin followed by treatment with tetrabutylammonium bromide is proposed. It was shown that the modified membrane significantly reduced the bromine permeability (crossover) with an acceptable increase in the ionic conductivity of the separator group. Leakage currents not exceeding 10-20 µA/cm² were achieved, and the Coulomb efficiency was over 90%. The ionic conductivity (at AC) of a membrane soaked with water was compared for different pretreatment conditions. The frequency dependence of the membrane resistance is shown. The features of the conduction mechanism of the modified membrane are discussed.

Modelling Sorption and Transport of Gases in Polymeric Membranes across Different Scales: A Review

Professor Giulio C. Sarti has provided outstanding contributions to the modelling of fluid sorption and transport in polymeric materials, with a special eye on industrial applications such as membrane separation, due to his Chemical Engineering background. He was the co-creator of innovative theories such as the Non-Equilibrium Theory for Glassy Polymers (NET-GP), a flexible tool to estimate the solubility of pure and mixed fluids in a wide range of polymers, and of the Standard Transport Model (STM) for estimating membrane permeability and selectivity. In this review, inspired by his rigorous and original approach to representing membrane fundamentals, we provide an overview of the most significant and up-to-date modeling tools available to estimate the main properties governing polymeric membranes in fluid separation, namely solubility and diffusivity. The paper is not meant to be comprehensive, but it focuses on those contributions that are most relevant or that show the potential to be relevant in the future. We do not restrict our view to the field of macroscopic modelling, which was the main playground of professor Sarti, but also devote our attention to Molecular and Multiscale Hierarchical Modeling. This work proposes a critical evaluation of the different approaches considered, along with their limitations and potentiality.

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.

Inverse Boltzmann Iterative Multi-Scale Molecular Dynamics Study between Carbon Nanotubes and Amino Acids

Our work uses Iterative Boltzmann Inversion (IBI) to study the coarse-grained interaction between 20 amino acids and the representative carbon nanotube CNT55L3. IBI is a multi-scale simulation method that has attracted the attention of many researchers in recent years. It can effectively modify the coarse-grained model derived from the Potential of Mean Force (PMF). IBI is based on the distribution result obtained by All-Atom molecular dynamics simulation; that is, the target distribution function and the PMF potential energy are extracted, and then, the initial potential energy extracted by the PMF is used to perform simulation iterations using IBI. Our research results have been through more than 100 iterations, and finally, the distribution obtained by coarse-grained molecular simulation (CGMD) can effectively overlap with the results of all-atom molecular dynamics simulation (AAMD). In addition, our work lays the foundation for the study of force fields for the simulation of the coarse-graining of super-large proteins and other important nanoparticles.

Longitudinal strand ordering leads to shear thinning in Nafion

Proton-exchange membrane fuel cells (PEMFC) offer a promising energy generation alternative for a wide range of technologies thanks to their ecological friendliness and unparalleled efficiency. At the heart of these electrochemical cells lies the membrane electrode assembly with its most important energy conversion components, the Proton Exchange Membrane. This component is created through the use of printing techniques and Nafion inks. The physicochemical properties of the ink, such as its viscosity under shear, are critical for the finished product. In this work we present non-equilibrium Molecular Dynamics simulations using a MARTINI based coarse-grained model for Nafion to understand the mechanism governing the shear viscosity of Nafion solutions. By simulating a Couette flow and calculating density maps of the Nafion chains in these simulations we shed light on the process that leads to the experimentally observed shear thinning effects of Nafion solutions under flow. We observe rod-shaped Nafion microstructures, 3 nm in size on average, when shear flow is absent or low. Higher shear rates instead break these structures and align Nafion strands along the direction of the flow, resulting in lower shear viscosities. Our work paves the way for a deeper understanding of the dynamic and mechanical properties of Nafion including studies of more complex CL and PEM inks.

CHARMM-GUI Polymer Builder for Modeling and Simulation of Synthetic Polymers

Molecular modeling and simulations are invaluable tools for polymer science and engineering, which predict physicochemical properties of polymers and provide molecular-level insight into the underlying mechanisms. However, building realistic polymer systems is challenging and requires considerable experience because of great variations in structures as well as length and time scales. This work describes Polymer Builder in CHARMM-GUI (http://www.charmm-gui.org/input/polymer), a web-based infrastructure that provides a generalized and automated process to build a relaxed polymer system. Polymer Builder not only provides versatile modeling methods to build complex polymer structures, but also generates realistic polymer melt and solution systems through the built-in coarse-grained model and all-atom replacement. The coarse-grained model parametrization is generalized and extensively validated with various experimental data and all-atom simulations. In addition, the capability of Polymer Builder for generating relaxed polymer systems is demonstrated by density calculations of 34 homopolymer melt systems, characteristic ratio calculations of 170 homopolymer melt systems, a morphology diagram of poly(styrene-b-methyl methacrylate) block copolymers, and self-assembly behavior of amphiphilic poly(ethylene oxide-b-ethylethane) block copolymers in water. We hope that Polymer Builder is useful to carry out innovative and novel polymer modeling and simulation research to acquire insight into structures, dynamics, and underlying mechanisms of complex polymer-containing systems.

Sensitivity of pair statistics on pair potentials in many-body systems

We study the sensitivity and practicality of Henderson's theorem in classical statistical mechanics, which states that the pair potential v(r) that gives rise to a given pair correlation function g₂(r) [or equivalently, the structure factor S(k)] in a classical many-body system at number density ρ and temperature T is unique up to an additive constant. While widely invoked in inverse-problem studies, the utility of the theorem has not been quantitatively scrutinized to any large degree. We show that Henderson's theorem has practical shortcomings for disordered and ordered phases for certain densities and temperatures. Using proposed sensitivity metrics, we identify illustrative cases in which distinctly different potential functions give very similar pair correlation functions and/or structure factors up to their corresponding correlation lengths. Our results reveal that due to a limited range and precision of pair information in either direct or reciprocal space, there is effective ambiguity of solutions to inverse problems that utilize pair information only, and more caution must be exercised when one claims the uniqueness of any resulting effective pair potential found in practice. We have also identified systems that possess virtually identical pair statistics but have distinctly different higher-order correlations. Such differences should be reflected in their individually distinct dynamics (e.g., glassy behaviors). Finally, we prove a more general version of Henderson's theorem that extends the uniqueness statement to include potentials that involve two- and higher-body interactions.

Martini Coarse-Grained Model of Hyaluronic Acid for the Structural Change of Its Gel in the Presence of Monovalent and Divalent Salts

Hyaluronic acid (HA) has a wide range of biomedical applications including the formation of hydrogels, microspheres, sponges, and films. The modeling of HA to understand its behavior and interaction with other biomolecules at the atomic level is of considerable interest. The atomistic representation of long HA polymers for the study of the macroscopic structural formation and its interactions with other polyelectrolytes is computationally demanding. To overcome this limitation, we developed a coarse grained (CG) model for HA adapting the Martini scheme. A very good agreement was observed between the CG model and all-atom simulations for both local (bonded interactions) and global properties (end-to-end distance, a radius of gyration, RMSD). Our CG model successfully demonstrated the formation of HA gel and its structural changes at high salt concentrations. We found that the main role of CaCl₂ is screening the electrostatic repulsion between chains. HA gel did not collapse even at high CaCl₂ concentrations, and the osmotic pressure decreased, which agrees well with the experimental results. This is a distinct property of HA from other proteins or polynucleic acids which ensures the validity of our CG model. Our HA CG model is compatible with other CG biomolecular models developed under the Martini scheme, which allows for large-scale simulations of various HA-based complex systems.

Molecular Modeling Investigations of Sorption and Diffusion of Small Molecules in Glassy Polymers

With a wide range of applications, from energy and environmental engineering, such as in gas separations and water purification, to biomedical engineering and packaging, glassy polymeric materials remain in the core of novel membrane and state-of the art barrier technologies. This review focuses on molecular simulation methodologies implemented for the study of sorption and diffusion of small molecules in dense glassy polymeric systems. Basic concepts are introduced and systematic methods for the generation of realistic polymer configurations are briefly presented. Challenges related to the long length and time scale phenomena that govern the permeation process in the glassy polymer matrix are described and molecular simulation approaches developed to address the multiscale problem at hand are discussed.