Ab initio molecular dynamics simulation of proton hopping in a model polymer membrane

Proton Transport in Perfluorinated Ionomer Simulated by Machine-Learned Interatomic Potential

Polymers are a class of materials that are highly challenging to deal with using first-principles methods. Here, we present an application of machine-learned interatomic potentials to predict structural and dynamical properties of dry and hydrated perfluorinated ionomers. An improved active-learning algorithm using a small number of descriptors allows to efficiently construct an accurate and transferable model for this multielemental amorphous polymer. Molecular dynamics simulations accelerated by the machine-learned potentials accurately reproduce the heterogeneous hydrophilic and hydrophobic domains formed in this material as well as proton and water diffusion coefficients under a variety of humidity conditions. Our results reveal pronounced contributions of Grotthuss chains consisting of two to three water molecules to the high proton mobility under strongly humidified conditions.

Influence of Mesoscale Interactions on Proton, Water, and Electrokinetic Transport in Solvent-Filled Membranes: Theory and Simulation

Transport of protons and water through water-filled, phase-separated cation-exchange membranes occurs through a network of interconnected nanoscale hydrophilic aqueous domains. This paper uses numerical simulations and theory to explore the role of the mesoscale network on water, proton, and electrokinetic transport in perfluorinated sulfonic acid (PFSA) membranes, pertinent to electrochemical energy-conversion devices. Concentrated-solution theory describes microscale transport. Network simulations model mesoscale effects and ascertain macroscopic properties. An experimentally consistent 3D Voronoi-network topology characterizes the interconnected channels in the membrane. Measured water, proton, and electrokinetic transport properties from literature validate calculations of macroscopic properties from network simulations and from effective-medium theory. The results demonstrate that the hydrophilic domain size affects the various microscale, domain-level transport modes dissimilarly, resulting in different distributions of microscale coefficients for each mode of transport. As a result, the network mediates the transport of species nonuniformly with dissimilar calculated tortuosities for water, proton, and electrokinetic transport coefficients (i.e., 4.7, 3.0, and 6.1, respectively, at a water content of 8 H₂O molecules per polymer charge equivalent). The dominant water-transport pathways across the membrane are different than those taken by the proton cation. Finally, the distribution of transport properties across the network induces local electrokinetic flows that couple water and proton transport; specifically, local electrokinetic transport induces water chemical-potential gradients that decrease macroscopic conductivity by up to a factor of 3. Macroscopic proton, water, and electrokinetic transport coefficients depend on the collective microscale transport properties of all modes of transport and their distribution across the hydrophilic domain network.

Controlling Hydronium Diffusivity in Model Proton Exchange Membranes

Fuel-cell-based proton exchange membranes (PEMs) show great potential as cost-effective and clean energy conversion devices. In our recent work, we found that for the low-hydrated model PEMs with a inhomogeneous water distribution and a sulfonate anionic functional end group (SO₃^-), the H₃O⁺ reacts with SO₃^- according to SO₃^- + H₃O⁺ ↔ SO₃H + H₂O, indicating that the anions in PEMs become active participants in the hydronium diffusion. In this work, we use fully atomistic ab initio molecular dynamics simulations to elucidate the optimal conditions that would promote the participation of SO₃^- in the hydronium diffusion mechanism by increasing the H₃O⁺/SO₃^- reactivity, thus increasing the hydronium diffusivity along the cell. The results presented in this work allow us to suggest a set of design rules for creating novel, highly conductive PEMs operating at high temperatures under a nonuniform water distribution using a linker/anion with a relatively high pK_a such as (CH₂)₂SO₃. We expect that the discovery of these key design principles will play an important role in the synthesis of high-performing materials for emerging PEM-based fuel cell technologies.

OH- and H3O+ Diffusion in Model AEMs and PEMs at Low Hydration: Insights from Ab Initio Molecular Dynamics

Fuel cell-based anion-exchange membranes (AEMs) and proton exchange membranes (PEMs) are considered to have great potential as cost-effective, clean energy conversion devices. However, a fundamental atomistic understanding of the hydroxide and hydronium diffusion mechanisms in the AEM and PEM environment is an ongoing challenge. In this work, we aim to identify the fundamental atomistic steps governing hydroxide and hydronium transport phenomena. The motivation of this work lies in the fact that elucidating the key design differences between the hydroxide and hydronium diffusion mechanisms will play an important role in the discovery and determination of key design principles for the synthesis of new membrane materials with high ion conductivity for use in emerging fuel cell technologies. To this end, ab initio molecular dynamics simulations are presented to explore hydroxide and hydronium ion solvation complexes and diffusion mechanisms in the model AEM and PEM systems at low hydration in confined environments. We find that hydroxide diffusion in AEMs is mostly vehicular, while hydronium diffusion in model PEMs is structural. Furthermore, we find that the region between each pair of cations in AEMs creates a bottleneck for hydroxide diffusion, leading to a suppression of diffusivity, while the anions in PEMs become active participants in the hydronium diffusion, suggesting that the presence of the anions in model PEMs could potentially promote hydronium diffusion.

Ion partitioning and permeation in charged low-T* membranes

Understanding ion transport in membrane materials is key to engineering and development of desalination and water purification technologies as well as electro-membrane applications. To date, modeling of ion transport has mainly relied on mean-field approaches, originally intended for weak inter-ionic interactions, i.e., high reduced temperature T*. This condition is violated in many membranes, which could explain disagreement between predicted trends and experiments. The paper highlights observed discrepancies and develops a new approach based on the concept of ion association, more adequate in the low-T^⁎ limit. The new model addresses ion binding and mobility consistently within the same physical picture, applied to different types of single and mixed salts. The resulting relations show a significantly weaker connection between ion partitioning and permeability than the standard ones. Estimates using primitive model (PM) of ions in a homogeneous dielectric suggest that non-PM mechanisms, originating from the molecular structure of the ion-solvating environment, might enhance ion association in membranes. PM analysis also predicts that ion solvation and association must be rigidly related, yet non-PM effects may decouple these phenomena and allow a crossover to non-trivial regimes consistent with experiments and simulations. Despite the crude nature of the presented approach and some questions remaining open, it appears to explain most available experimental data and presents a step towards predictive modeling of ion-selective membrane separations in water-, environment- and energy-related applications.

Molecular Modeling of Structure and Dynamics of Nafion Protonation States

We present the results of the atomistic molecular dynamics modeling of different protonation states of Nafion at varying hydration levels. Previous experiments have shown that the degree of deprotonation (DDP) of the sulfonic acid groups in a Nafion membrane varies significantly upon hydration. Our goal is to provide insights into the effects of variable protonation states and water content on the internal structure and vehicular transport inside the Nafion membrane. The Nafion side chain lengths showed a weak increasing trend with increasing DDP at all hydration levels, exposing more of the sulfonic acid groups to the hydrophilic/water phase. The water-phase characteristic size/diameter decreased with increasing DDP, but, interestingly, the average number of water molecules per cluster increased. The probability of water-hydronium hydrogen bond formation decreased with increasing DDP, despite an increase in the total number of such hydrogen bonds. The water diffusion was largely unaffected by the state of deprotonation. In contrast to that, the hydronium ion diffusion slowed down with increasing DDP in the overall membrane. The hydronium ion residence times around the sulfonic acid group increased with increasing DDP. Our simulations show a strong connection between the morphology of the water domains and protonation states of Nafion. Such a connection can also be expected in polyelectrolyte membranes similar to Nafion.

Modeling and simulation of membrane process

The article presents the different approaches to polymer membrane mathematical modeling. Traditional models based on experimental physicochemical correlations and balance models are presented in the first part. Quantum and molecular mechanics models are presented as they are more popular for polymer membranes in fuel cells. The initial part is enclosed by neural network models which found their use for different types of processes in polymer membranes. The second part is devoted to models of fluid dynamics. The computational fluid dynamics technique can be divided into solving of Navier-Stokes equations and into Boltzmann lattice models. Both approaches are presented focusing on membrane processes.

New Insights into Perfluorinated Sulfonic-Acid Ionomers

In this comprehensive review, recent progress and developments on perfluorinated sulfonic-acid (PFSA) membranes have been summarized on many key topics. Although quite well investigated for decades, PFSA ionomers' complex behavior, along with their key role in many emerging technologies, have presented significant scientific challenges but also helped create a unique cross-disciplinary research field to overcome such challenges. Research and progress on PFSAs, especially when considered with their applications, are at the forefront of bridging electrochemistry and polymer (physics), which have also opened up development of state-of-the-art in situ characterization techniques as well as multiphysics computation models. Topics reviewed stem from correlating the various physical (e.g., mechanical) and transport properties with morphology and structure across time and length scales. In addition, topics of recent interest such as structure/transport correlations and modeling, composite PFSA membranes, degradation phenomena, and PFSA thin films are presented. Throughout, the impact of PFSA chemistry and side-chain is also discussed to present a broader perspective.

Moisture and thermal expansion properties and mechanism of interaction between ions of a Nafion-based membrane electrode assembly

The coefficient of moisture and thermal expansion for each layer of membrane electrode assembly.

Water Diffusion Dependence on Amphiphilic Block Design in (Amphiphilic-Hydrophobic) Diblock Copolymer Membranes

Polyelectrolyte membranes (PEMs) are applied in polyelectrolyte fuel cells (PEFC). The proton conductive pathways within PEMs are provided by nanometer-sized water containing pores. Large-scale application of PEFC requires the production of low-cost membranes with high proton conductivity and therefore good connected pore networks. Pore network formation within four alternative model diblock (hydrophobic_amphiphilic) copolymers in the presence of water is studied by dissipative particle dynamics. Each hydrophobic block contains 50 consecutively connected hydrophobic (A) fragments, and amphiphilic blocks contain 40 hydrophobic A beads and 10 hydrophilic C beads. For one amphiphilic block the C beads are distributed uniformly along the backbone. For the other architectures C beads are located at the end of the side chains attached at regular intervals along the backbone. Water diffusion through the pores is modeled by Monte Carlo tracer diffusion through mapped morphologies. Diffusion is highest for the grafted architectures and increases with increase of length of the side chains. A consistent picture emerges in which diffusion strongly increases with the value of ⟨Nbond⟩ within the amphiphilic block, where ⟨Nbond⟩ is the average number of bonds between hydrophobic A beads and the nearest C bead.

Modelling linear and branched amphiphilic star polymer electrolyte membranes and verification of the bond counting method

Water diffusion through hydrated amphiphilic star polymer membranes depends strongly on hydrophilic position within the linear and Y-shaped arms.

Water diffusion within hydrated model grafted polymeric membranes with bimodal side chain length distributions

The effect of bimodal side chain length distributions on pore morphology and solvent diffusion within hydrated amphiphilic polymeric membranes is predicted. Seven polymeric architectures are constructed from hydrophobic backbones from which at regular intervals side chains branch off that are alternatingly short (composed of p hydrophobic A fragments or beads) and long (q A fragments, q > p). The side chains are end-linked with a hydrophilic C fragment. Pore morphologies at a water volume fraction of 0.16 are calculated by dissipative particle dynamics (DPD). Water diffusion through the water containing pores is calculated by tracer diffusion calculations through 140 selected snapshots and from the water bead motions. Diffusion constants decrease with difference in side chain lengths, q - p. Overall, the distance between pores also decreases with q - p. The results are explained by counting for every architecture the average number of bonds 〈N(bond)〉 between an A and the nearest C fragment. These results are in line with a database that contains more than 60 architectures. Diffusion constants tend to increase linearly with 〈N(bond)〉|C|(-1)|A|, where |C| and |A| are the C and A bead fractions within the architecture. 〈N(bond)〉 is therefore expected to be an interesting design parameter for obtaining low percolation thresholds for solvent and/or proton diffusion.

Persistent Subdiffusive Proton Transport in Perfluorosulfonic Acid Membranes

Proton transport (PT) in solutions of small amphiphiles in water has previously been shown to be subdiffusive for long times. The present study analyzes simulations of hydrated perfluorosulfonic acid (PFSA) membranes in order to determine whether PT is also subdiffusive in these important amphiphilic systems. We show that PT is indeed subdiffusive for several hundred picoseconds for all hydration levels examined, and the subdiffusive behavior is highly dependent on water concentration. We also investigate the caging of the excess proton using a recently developed technique and show that the excess proton exhibits caging effects up to at least 1 ns in PFSA systems. In order to fully characterize the long-time behavior of PT in PFSAs, these results demonstrate that multiple nanosecond trajectories are needed, well beyond the current capabilities of ab initio molecular dynamics.