Nafion Ionomer Dispersion in Mixtures of 1‐Propanol and Water Based on the Martini Coarse‐Grained Model

Colloidal Stability of PFSA-Ionomer Dispersions. Part I. Single-Ion Electrostatic Interaction Potential Energies

Charged colloidal particles neutralized by a single counterion are increasingly important for many emerging technologies. Attention here is paid specifically to hydrogen fuel cells and water electrolyzers whose catalyst layers are manufactured from a perfluorinated sulfonic acid polymer (PFSA) suspended in aqueous/alcohol solutions. Partially dissolved PFSA aggregates, known collectively as ionomers, are stabilized by the electrostatic repulsion of overlapping diffuse double layers consisting of only protons dissociated from the suspended polymer. We denote such double layers containing no added electrolyte as "single ion". Size-distribution predictions build upon interparticle interaction potential energies from the Derjaguin-Landau-Verwey-Overbeek (DLVO) formalism. However, when only a single counterion is present in solution, classical DLVO electrostatic potential energies no longer apply. Accordingly, here a new formulation is proposed to describe how single-counterion diffuse double layers interact in colloidal suspensions. Part II (Srivastav, H.; Weber, A. Z.; Radke, C. J. Langmuir 2024 DOI: 10.1021/acs.langmuir.3c03904) of this contribution uses the new single-ion interaction energies to predict aggregated size distributions and the resulting solution pH of PFSA in mixtures of n-propanol and water. A single-counterion diffuse layer cannot reach an electrically neutral concentration far from a charged particle. Consequently, nowhere in the dispersion is the solvent neutral, and the diffuse layer emanating from one particle always experiences the presence of other particles (or walls). Thus, in addition to an intervening interparticle repulsive force, a backside osmotic force is always present. With this new construction, we establish that single-ion repulsive pair interaction energies are much larger than those of classical DLVO electrostatic potentials. The proposed single-ion electrostatic pair potential governs dramatic new dispersion behavior, including dispersions that are stable at a low volume fraction but unstable at a high volume fraction and finite volume-fraction dispersions that are unstable with fine particles but stable with coarse particles. The proposed single-counterion electrostatic pair potential provides a general expression for predicting colloidal behavior for any charged particle dispersion in ionizing solvents with no added electrolyte.

Colloidal Nafion Particles: Are Cylinders Ubiquitous?

Colloidal Nafion morphology plays a critical role in determining the performance of fuel cells and electrolyzers. While small-angle neutron scattering (SANS) studies previously described Nafion in liquid media as dispersed cylinders, the analysis remains nonunique with multiple possible morphological descriptions of the data. Here, using SANS and all-atomistic molecular dynamics, we confirm that Nafion morphology in liquid media differs substantially depending on dispersing agent and dispersion method. H+ Nafion dispersed in N-methyl pyrrolidone forms swollen cluster particles with physically cross-linked ionic groups. Scattering profiles from dispersed Nafion membrane have a large structure factor feature not observed for redispersed Nafion D-521. H+ Nafion dispersed in water has a highly elongated cylindrical morphology (radius = 10 ± 1.5 Å, height = 358 ± 4.7 Å) with fully dissociated and solvated sulfonic acid groups on the particle wall. These results highlight an important discrepancy between the methods of preparing Nafion dispersions and the use of simplified analysis techniques to describe Nafion morphology.

Deep Learning to Reveal the Distribution and Diffusion of Water Molecules in Fuel Cell Catalyst Layers

Water management in the catalyst layers (CLs) of proton-exchange membrane fuel cells is crucial for its commercialization and popularization. However, the high experimental or computational cost in obtaining water distribution and diffusion remains a bottleneck in the existing experimental methods and simulation algorithms, and further mechanistic exploration at the nanoscale is necessary. Herein, we integrate, for the first time, molecular dynamics simulation with our customized analysis framework based on a multiattribute point cloud dataset and an advanced deep learning network. This was achieved through our workflow that generates simulated transport data of water molecules in the CLs as the training and test dataset. Deep learning framework models the multibody solid-liquid system of CLs on a molecular scale and completes the mapping from the Pt/C substrate structure and Nafion aggregates to the density distribution and diffusion coefficient of water molecules. The prediction results are comprehensively analyzed and error evaluated, which reveals the highly anisotropic interaction landscape between 50,000 pairs of interacting nanoparticles and explains the structure and water transport property relationship in the hydrated Nafion film on the molecular scale. Compared to the conventional methods, the proposed deep learning framework shows computational cost efficiency, accuracy, and good visual display. Further, it has a generality potential to model macro- and microscopic mass transport in different components of fuel cells. Our framework is expected to make real-time predictions of the distribution and diffusion of water molecules in CLs as well as establish statistical significance in the structural optimization and design of CLs and other components of fuel cells.

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.

The Martini Model in Materials Science

The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.