Molecular Dynamics Simulations of Substrate Hydrophilicity and Confinement Effects in Capped Nafion Films

Electro-osmotic Drag and Thermodynamic Properties of Water in Hydrated Nafion Membranes from Molecular Dynamics

One of the important parameters in water management of proton exchange membranes is the electro-osmotic drag (EOD) coefficient of water. The value of the EOD coefficient is difficult to justify, and available literature data on this for Nafion membranes show scattering from in experiments and simulations. Here, we use a classical all-atom model to compute the EOD coefficient and thermodynamic properties of water from molecular dynamics simulations for temperatures between 330 and 420 K, and for different water contents between λ = 5 and λ = 20. λ is the ratio between the moles of water molecules to the moles of sulfonic acid sites. This classical model does not capture the Grotthuss mechanism; however, it is shown that it can predict the EOD coefficient within the range of experimental values for λ = 5 where the vehicular mechanism dominates proton transfer. For λ > 5, the Grotthuss mechanism becomes dominant. To obtain the EOD coefficient, average velocities of water and ions are computed by imposing different electric fields to the system. Our results show that the velocities of water and hydronium scale linearly with the electric field, resulting in a constant ratio of ca. 0.4 within the error bars. We find that the EOD coefficient of water linearly increases from 2 at λ = 5 to 8 at λ = 20 and the results are not sensitive to temperature. The EOD coefficient at λ = 5 is within the range of experimental values, confirming that the model can capture the vehicular transport of protons well. At λ = 20, due to the absence of proton hopping in the model, the EOD coefficient is overestimated by a factor of 3 compared to experimental values. To analyze the interactions between water and Nafion, the partial molar enthalpies and partial molar volumes of water are computed from molecular dynamics simulations. At different water uptakes, multiple linear regression is used on raw simulation data within a narrow composition range of water inside the Nafion membrane. The partial molar volumes and partial molar excess enthalpies of water asymptotically approach the molar volumes and molar excess enthalpies of pure water for water uptakes above 5. This confirms the model can capture the bulklike behavior of water in the Nafion at high water uptakes.

Effect of Nafion content and hydration level on the electrochemical area of a Pt nanocatalyst in the triple-phase boundary

Despite the great scientific effort, there are still some aspects of a polymeric membrane-based fuel cell (PEMFC) operation that are difficult to access experimentally. This is the case of the so-called triple-phase boundary (TPB), where the ionomer (commonly Nafion) interacts with the supported nanocatalyst (commonly Pt) and is key to the catalytic activity of the system. In this work, we use molecular dynamics simulations and electrochemical experiments on a Nafion/Pt/C system. We perform a systematic analysis, at an atomistic level, to evaluate the effect of several fundamental factors and their intercorrelation on the electrochemically active area (ECSA) of the catalysts. Our results reveal that at high Nafion contents, the catalyst utilization is affected due to the strong interaction between the sulfonic groups of the ionomer and the surface of the Pt nanoparticles (NPs). On the other hand, when the hydration level of the membrane decreases, the sulfonic groups have a greater occupation on the NP surface, covering the active area with hydrophobic Nafion chains and therefore increasing the inactive area. Voltammograms can corroborate our calculations. Overall, this investigation allows us to rationalize how the catalyst utilization is affected, which is an important step in establishing the relationship between the environment and the effectiveness and durability of the PEMFC system.

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.

Molecular dynamics of ionic polymer-metal composites

Ionic polymer-metal composites (IPMCs) constitute a promising class of soft, active materials with potentially ubiquitous use in science and engineering. Realizing the full potential of IPMCs calls for a deeper understanding of the mechanisms underpinning their most intriguing characteristics: the ability to deform under an electric field and the generation of a voltage upon mechanical deformation. These behaviours are tightly linked to physical phenomena at the level of atoms, including rearrangements of ions and molecules, along with the formation of sub-nanometre thick double layers on the surface of the metal electrodes. Several continuum theories have been developed to describe these phenomena, but their experimental and theoretical validation remains incomplete. IPMC modelling at the atomistic scale could beget valuable support for these efforts, by affording granular analysis of individual atoms. Here, we present a simplified atomistic model of IPMCs based on classical molecular dynamics. The three-dimensional IPMC membrane is constrained by two smooth walls, a simplified analogue of metal electrodes, impermeable only to counterions. The electric field is applied as an additional force acting on all the atoms. We demonstrate the feasibility of simulating counterions' migration and pile-up upon the application of an electric field, similar to experimental observations. By analysing the spatial configuration of atoms and stress distribution, we identify two mechanisms for stress generation. The presented model offers new insight into the physical underpinnings of actuation and sensing in IPMCs. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Solvation-Driven Electrochemical Actuation

We demonstrate a novel principle of contactless actuation for ionic membranes in salt solution based on solvation. Actuation is driven by differential swelling of the sides of the membrane, due to comigrating water in the solvation shells of mobile ions. We validate our theory through a series of experiments, which unravel a strong dependence of membrane deflection on the hydration numbers of mobile ions in the external solution and membrane. Our study suggests a critical role of solvation in the chemoelectromechanics of natural and artificial selectively permeable membranes.

Mechanism of Water Content on the Electrochemical Surface Area of the Catalyst Layer in the Proton Exchange Membrane Fuel Cell

Combined molecular dynamics (MD) simulation and experiment are adopted to gain the mechanism of water content on the electrochemical surface area (ECSA) of the catalyst layer in a proton exchange membrane fuel cell. The morphology of water domains in the catalyst layer has a strong impact on the ECSA via MD simulation. The morphology of the water domains is isolated water clusters at low water content, resulting in the poor ECSA due to the lack of proton transport paths. The transport paths of protons tend to be quickly established with increasing water content during the transition process of the morphology of water domains from isolated water clusters to the water channel network, thereby leading to the rapid increase of the ECSA. However, the slight increase of the ECSA at high water content mainly results from the improved contact area between water domains and Pt particle instead of the formation of new transport paths. In addition, the stronger binding of water molecules and the Pt particle at low temperature results in a higher ECSA.

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.

Modification of Graphene Oxide Membranes by the Incorporation of Nafion Macromolecules and Conductive Scaffolds

Stable, reproducible and low-cost graphene oxide (GO)/Nafion (N) membranes were fabricated using electronically conductive carbon paper (CP) matts as a scaffold. The presence of polar groups in the Nafion molecule facilitates the strong interaction with functional groups in the GO, which increases GO dispersion and aids the retention of the composite into the CP scaffold. Distribution of GO/N was carefully characterized by X-ray diffraction work function measurements, Raman and scanning electron microscopy analyses. The performance of these membranes was tested with 1 M NaCl at standard conditions, finding 85% ion removal in the best membranes by a mixed ion rejection/retention mechanism. The Nafion provided mechanical stability and fixed negative charge to the membranes, and its micellar organization, segregation and confinement favored ion rejection in Nafion-rich areas. The good electronic conductivity of these membranes was also demonstrated, allowing for the application of a small potential bias to enhance membrane performance in future studies.