Simulation Study of the Correlation between Structure and Conductivity in Stretched Nafion

Electrostatic Potential of Functional Cations as a Predictor of Hydroxide Diffusion Pathways in Nanoconfined Environments of Anion Exchange Membranes

Nanoconfined anion exchange membranes (AEMs) play a vital role in emerging electrochemical technologies. The ability to control dominant hydroxide diffusion pathways is an important goal in the design of nanoconfined AEMs. Such control can shorten hydroxide transport pathways between electrodes, reduce transport resistance, and enhance device performance. In this work, we propose an electrostatic potential (ESP) approach to explore the effect of the polymer electrolyte cation spacing on hydroxide diffusion pathways from a molecular perspective. By exploring cation ESP energy surfaces and validating outcomes through prior ab initio molecular dynamics simulations of nanoconfined AEMs, we find that we can achieve control over preferred hydroxide diffusion pathways by adjusting the cation spacing. The results presented in this work provide a unique and straightforward approach to predict preferential hydroxide diffusion pathways, enabling efficient design of highly conductive nanoconfined AEM materials for electrochemical technologies.

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.

Electrospun Ion-Conducting Composite Membrane with Buckling-Induced Anisotropic Through-Plane Conductivity

Fuel cell (FC) is an attractive green alternative for today's fuel combustion systems. In common FCs, a polymer electrolyte membrane selectively conducts protons but blocks the passage of electrons and fuel. Nafion, the current benchmark membrane material, has a superior conductivity owing to unique morphology comprising randomly oriented elongated ionic nanochannels within its Teflon-like matrix. Channel orientation enhances Nafion conductivity, yet there has been no facile method to induce a stable alignment in the desired through-plane (TP) direction. Here, we report an approach based on dual electrospun Nafion-PVDF nanofiber composites that yields a stable TP alignment. It utilizes extreme thinness and strong inherent orientation within electrospun nanofibers, which is readily converted to TP alignment by plunging an electrospun nanofiber mat into a thin slit, resulting in nanofiber buckling and subsequent consolidation. Using TEM and SAXS, we demonstrate a pronounced and sustained TP ion channel orientation in prepared membranes, yielding a highly anisotropic swelling and conductivity exceeding that of bulk Nafion when normalized to Nafion content. The analysis also highlights the importance of PVDF as a stabilizing component, preserving orientation upon annealing, while a similarly prepared pure Nafion membrane loses anisotropy. The approach holds potential to advance the FC technology by overcoming current limitations of ionomeric membranes.

Mechanical behavior of a hydrated perfluorosulfonic acid membrane at meso and nano scales

Perfluorosulfonic acid (PFSA) is widely used as the membrane material for proton-exchange membrane fuel cells, and its mechanical properties directly affect the stability and the life of the internal structure of the proton exchange membrane.

Coarse-grained model of nanoscale segregation, water diffusion, and proton transport in Nafion membranes

We present a coarse-grained model of the acid form of Nafion membrane that explicitly includes proton transport. This model is based on a soft-core bead representation of the polymer implemented into the dissipative particle dynamics (DPD) simulation framework. The proton is introduced as a separate charged bead that forms dissociable Morse bonds with water beads. Morse bond formation and breakup artificially mimics the Grotthuss hopping mechanism of proton transport. The proposed DPD model is parameterized to account for the specifics of the conformations and flexibility of the Nafion backbone and sidechains; it treats electrostatic interactions in the smeared charge approximation. The simulation results qualitatively, and in many respects quantitatively, predict the specifics of nanoscale segregation in the hydrated Nafion membrane into hydrophobic and hydrophilic subphases, water diffusion, and proton mobility. As the hydration level increases, the hydrophilic subphase exhibits a percolation transition from a collection of isolated water clusters to a 3D network of pores filled with water embedded in the hydrophobic matrix. The segregated morphology is characterized in terms of the pore size distribution with the average size growing with hydration from ∼1 to ∼4 nm. Comparison of the predicted water diffusivity with the experimental data taken from different sources shows good agreement at high and moderate hydration and substantial deviation at low hydration, around and below the percolation threshold. This discrepancy is attributed to the dynamic percolation effects of formation and rupture of merging bridges between the water clusters, which become progressively important at low hydration, when the coarse-grained model is unable to mimic the fine structure of water network that includes singe molecule bridges. Selected simulations of water diffusion are performed for the alkali metal substituted membrane which demonstrate the effects of the counter-ions on membrane self-assembly and transport. The hydration dependence of the proton diffusivity reproduces semi-qualitatively the trend of the diverse experimental data, showing a sharp decrease around the percolation threshold. Overall, the proposed model opens up an opportunity to study self-assembly and water and proton transport in polyelectrolytes using computationally efficient DPD simulations, and, with further refinement, it may become a practical tool for theory informed design and optimization of perm-selective and ion-conducting membranes with improved properties.

Water sub-diffusion in membranes for fuel cells

We investigate the dynamics of water confined in soft ionic nano-assemblies, an issue critical for a general understanding of the multi-scale structure-function interplay in advanced materials. We focus in particular on hydrated perfluoro-sulfonic acid compounds employed as electrolytes in fuel cells. These materials form phase-separated morphologies that show outstanding proton-conducting properties, directly related to the state and dynamics of the absorbed water. We have quantified water motion and ion transport by combining Quasi Elastic Neutron Scattering, Pulsed Field Gradient Nuclear Magnetic Resonance, and Molecular Dynamics computer simulation. Effective water and ion diffusion coefficients have been determined together with their variation upon hydration at the relevant atomic, nanoscopic and macroscopic scales, providing a complete picture of transport. We demonstrate that confinement at the nanoscale and direct interaction with the charged interfaces produce anomalous sub-diffusion, due to a heterogeneous space-dependent dynamics within the ionic nanochannels. This is irrespective of the details of the chemistry of the hydrophobic confining matrix, confirming the statistical significance of our conclusions. Our findings turn out to indicate interesting connections and possibilities of cross-fertilization with other domains, including biophysics. They also establish fruitful correspondences with advanced topics in statistical mechanics, resulting in new possibilities for the analysis of Neutron scattering data.

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.

Fabrication and electrolyte characterization of uniaxially-aligned anion conductive polymer nanofibers

We report the anion transport properties of anion conductive polymer nanofibers fabricated using an electrospinning method. The aligned nanofibers were prepared to evaluate the anion conductivity of the nanofibers. The aligned nanofibers had 10-15 times higher conductivity (up to 160 mS cm^-1 at 90 °C and 95%RH) and lower activation energy (23-25 kJ mol^-1) than the corresponding membranes, even though the nanofibers showed lower water uptake than the corresponding membranes. The anion conductivity measurement of nanofibers with different IEC values and anion species revealed that the dependency of anion conductivity on these factors was smaller in the nanofibers than in the corresponding membranes. These results indicate that effective ion transport pathways were formed in the nanofibers due to the phase separation and the polymer chain orientation along the nanofiber axis. These nanofibers are expected to be useful for future applications in alkaline fuel cells, air batteries, and other energy- and environment-related devices regardless of the ion species.

Sub-diffusion and population dynamics of water confined in soft environments

We have studied by using molecular dynamics computer simulations the dynamics of water confined in ionic surfactant phases, ranging from well ordered lamellar structures to micelles at low and high water loading, respectively. We have analysed in depth the main dynamical features in terms of mean-squared displacements and intermediate scattering functions, and found clear evidence of sub-diffusive behaviour. We have identified water molecules lying at the charged interface with the hydrophobic confining matrix as the main factor responsible for this unusual feature, and given a comprehensive picture of dynamics based on a very precise analysis of lifetimes at the interface. We conclude by providing, for the first time to our knowledge, a unique framework for rationalizing the existence of important dynamical heterogeneities in fluids adsorbed in soft confining environments.

Structure and property optimization of perfluorinated short side chain membranes for hydrogen fuel cells using orientational stretching

Stretching of membranes with low molecular weight makes structure rearrangement according to neutron scattering data on D₂O-filled membranes.

Water confined in self-assembled ionic surfactant nano-structures

We present a coarse-grained model for ionic surfactants in explicit aqueous solutions, and study by computer simulation both the impact of water content on the morphology of the system, and the consequent effect of the formed interfaces on the structural features of the absorbed fluid. On increasing the hydration level under ambient conditions, the model exhibits a series of three distinct phases: lamellar, cylindrical and micellar. We characterize the different structures in terms of diffraction patterns and neutron scattering static structure factors. We demonstrate that the rate of variation of the nano-metric sizes of the self-assembled water domains shows peculiar changes in the different phases. We also analyse in depth the structure of the water/confining matrix interfaces, the implications of their tunable degree of curvature, and the properties of water molecules in different restricted environments. Finally, we compare our results with experimental data and their impact on a wide range of important scientific and technological domains, where the behavior of water at the interface with soft materials is crucial.

Atomistic simulations of the structure of highly crosslinked sulfonated poly(styrene-co-divinylbenzene) ion exchange resins

The microscopic structures of highly crosslinked sulfonated poly(styrene-co-divinylbenzene) resins have been modeled by generating atomistic microstructures using stochastic-like algorithms, which are subsequently relaxed using molecular dynamics. Two different generation algorithms have been tested. The relaxation of the microstructures generated by the first algorithm, which is based on a homogeneous construction of the resin, leads to a significant overestimation of the experimental density as well as to an unsatisfactory description of the porosity. In contrast, the generation approach that combines algorithms for the heterogeneous growing and branching of the chains enables the formation of crosslinks with different topologies. In particular, the intrinsic heterogeneity observed in these resins is efficiently reproduced when the topological loops, which are defined by two or more crosslinks closing a cycle, are present in their microscopic description. Thus, the apparent density, porosity and pore volume estimated using microstructures with these topological loops, called super-crosslinks, are in very good agreement with the experimental results. Although the backbone dihedral angle distribution of the generated and relaxed models is not influenced by the topology, the number and type of crosslinks affect the medium- and long-range atomic disposition of the backbone atoms and the distribution of sulfonic groups. An analysis of the distribution of the local density indicates that super-crosslinks are responsible for the heterogeneous homogenization observed during the MD relaxation. Finally the π-π stacking interactions have been analyzed. Results indicate that those in which the two rings adopt a T-shaped disposition are considerably more abundant as compared to those with the co-facially oriented rings, independently of the resin topology.

Ultra-high proton conduction in electrospun sulfonated polyimide nanofibers

Uniaxially-aligned sulfonated polyimide nanofibers fabricated by an electrospinning method showed ultra-high proton conductivities due to molecular orientation in the nanofibers.

Large-scale atomistic and quantum-mechanical simulations of a Nafion membrane: Morphology, proton solvation and charge transport

Atomistic and first-principles molecular dynamics simulations are employed to investigate the structure formation in a hydrated Nafion membrane and the solvation and transport of protons in the water channel of the membrane. For the water/Nafion systems containing more than 4 million atoms, it is found that the observed microphase-segregated morphology can be classified as bicontinuous: both majority (hydrophobic) and minority (hydrophilic) subphases are 3D continuous and organized in an irregular ordered pattern, which is largely similar to that known for a bicontinuous double-diamond structure. The characteristic size of the connected hydrophilic channels is about 25-50 Å, depending on the water content. A thermodynamic decomposition of the potential of mean force and the calculated spectral densities of the hindered translational motions of cations reveal that ion association observed with decreasing temperature is largely an entropic effect related to the loss of low-frequency modes. Based on the results from the atomistic simulation of the morphology of Nafion, we developed a realistic model of ion-conducting hydrophilic channel within the Nafion membrane and studied it with quantum molecular dynamics. The extensive 120 ps-long density functional theory (DFT)-based simulations of charge migration in the 1200-atom model of the nanochannel consisting of Nafion chains and water molecules allowed us to observe the bimodality of the van Hove autocorrelation function, which provides the direct evidence of the Grotthuss bond-exchange (hopping) mechanism as a significant contributor to the proton conductivity.

Simulation study of poled low-water ionomers with different architectures

The role of the ionomer architecture in the formation of ordered structures in poled membranes is investigated by molecular dynamics computer simulations. It is shown that the length of the sidechain L(s) controls both the areal density of cylindrical aggregates N(c) and the diameter of these cylinders in the poled membrane. The backbone segment length L(b) tunes the average diameter D(s) of cylindrical clusters and the average number of sulfonates N(s) in each cluster. A simple empirical formula is noted for the dependence of the number density of induced rod-like aggregates on the sidechain length L(s) within the parameter range considered in this study.

Linear coupling of alignment with transport in a polymer electrolyte membrane

Polymer electrolyte membranes (PEMs) selectively transport ions and polar molecules in a robust yet formable solid support. Tailored PEMs allow for devices such as solid-state batteries,'artificial muscle' actuators and reverse-osmosis water purifiers. Understanding how PEM structure and morphology relate to mobile species transport presents a challenge for designing next-generation materials. Material length scales from subnanometre to 1 μm influence bulk properties such as ion conductivity and water transport. Here we employ multi-axis pulsed-field-gradient NMR to measure diffusion anisotropy, and (2)H NMR spectroscopy and synchrotron small-angle X-ray scattering to probe orientational order as a function of water content and of membrane stretching. Strikingly, transport anisotropy linearly depends on the degree of alignment, signifying that membrane stretching affects neither the nanometre-scale channel dimensions nor the defect structure,causing only domain reorientation. The observed reorientation of anisotropic domains without perturbation of the inherent nematic-like domain character parallels the behaviour of nematic elastomers, promises tailored membrane conduction and potentially allows understanding of tunable shape-memory effects in PEM materials. This quantitative understanding will drive PEM design efforts towards optimal membrane transport, thus enabling more efficient polymeric batteries, fuel cells, mechanical actuators and water purification.

Ab initio molecular dynamics of proton networks in narrow polymer electrolyte pores

It is well established that proton conductivity in fuel cell membrane materials such as Nafion decreases strongly with decreasing water content. Proton transport in almost dry membranes is thought to proceed through narrow channels. In the present work we investigate proton structure and dynamics in two narrow cylindrical pores, which differ by their radius and the spacing of SO(3)H groups inside the channel. Pores are modelled through eight CF(3)CF(3) and four CF(3)SO(3)H entities in a helical arrangement. The water content λ (the ratio between the number of water molecules and the number of sulfonic acid groups) in the pores varies between 2.5 and 4.5. We observe a transition from the undissociated acid at very low λ through more or less localized H(3)O(+) entities to more delocalized H(5)O(2)(+) entities for the investigated range of λ. In the narrower pore, where S-S distances vary in a more favourable range (between 6 and 8.5 Å) than in the wider pore, we find that the molecular mobility is significantly higher, even at a rather high density of water molecules inside the pore.

A simulation study of field-induced proton-conduction pathways in dry ionomers

The morphological changes that can be induced in a dry ionomer by application of a strong electric field have been studied by means of computer simulation. The internal energy of the membrane at first slowly decreases with increasing field, but then rapidly increases after a certain threshold field is reached. This effect is interpreted as the reorganization of interacting head group dipoles in response to the external perturbation. The resulting morphology contains continuous channels of hydrophilic material capable of facilitating proton conduction. Upon removal of the poling field, the system does not return to its original morphology, but retains the anisotropic structure of the poled material. The poled structure appears to be thermodynamically stable, as confirmed by calculations of the Helmholtz energy of the original and poled samples.

Simulation study of field-induced morphological changes in a proton-conducting ionomer

A simulation study was made of the effects of strong electric fields on the morphology of a Nafion-like ionomer at various levels of hydration. The results of united-atom molecular-dynamics computations showed a self-organization of the side chain terminal groups into cylindrical clusters. The walls of these clusters contain the sulfonate dipoles, while the interior holds the majority of the water molecules. These cylindrical structures then align to form an hexatic array aligned along the direction of the applied electric field. The hexatic morphology persists after the removal of the field. A calculation by means of the Kirkwood coupling parameter method shows the Helmholtz free energy of the hexatic morphology of the poled membrane to be lower than that of the initial isotropic material, even in the absence of the applied field.

Intermolecular momentum transfer in poly(perfluorosulfonic acid) membrane hydrated by aqueous solution of methanol: A molecular dynamics simulation study

Intermolecular momentum transfer in methanol-water mixture solvated poly(perfluoro-sulfonic acid) membrane is studied in terms of center of mass velocity cross-correlation functions between molecular mass centers in their first coordination shells based on molecular dynamics simulations. Moreover, the center of mass velocity cross-correlation functions are also decomposed into longitudinal and transversal contributions. The fastest momentum transfer is observed between hydronium cation and water molecule due to the strong hydrogen bond interaction. The center of mass velocity cross-correlation functions reach peak value in about 36 fs, corresponding to a single collision with a neighboring molecule. For the momentum transfer between the water molecule and methanol molecule, the peaking time is 70 fs or about twice of that between hydronium cation and water molecule. Oscillation of the center of mass velocity cross-correlation functions between hydronium cation and water molecule is also observed due to the cage effect in their equilibrium positions.

Simulation study of sulfonate cluster swelling in ionomers

We have performed simulations to study how increasing humidity affects the structure of Nafion-like ionomers under conditions of low sulfonate concentration and low humidity. At the onset of membrane hydration, the clusters split into smaller parts. These subsequently swell, but then maintain constant the number of sulfonates per cluster. We find that the distribution of water in low-sulfonate membranes depends strongly on the sulfonate concentration. For a relatively low sulfonate concentration, nearly all the side-chain terminal groups are within cluster formations, and the average water loading per cluster matches the water content of membrane. However, for a relatively higher sulfonate concentration the water-to-sulfonate ratio becomes nonuniform. The clusters become wetter, while the intercluster bridges become drier. We note the formation of unusual shells of water-rich material that surround the sulfonate clusters.

Predicted electric-field-induced hexatic structure in an ionomer membrane

Coarse-grained molecular-dynamics simulations were used to study the morphological changes induced in a Nafion-like ionomer by the imposition of a strong electric field. We observe the formation of structures aligned along the direction of the applied field. The polar head groups of the ionomer sidechains aggregate into clusters, which then form rodlike formations which assemble into a hexatic array aligned with the direction of the field. Occasionally these lines of sulfonates and protons form a helical structure. Upon removal of the electric field, the hexatic array of rodlike structures persists and has a lower calculated free energy than the original isotropic morphology.