Atomistic simulation of water percolation and proton hopping in Nafion fuel cell membrane

Effects of Water Content on the Transport and Thermodynamic Properties of Phosphonium Ionic Liquids

We present a numerical investigation of the influence of water content on the dynamic properties of a family of phosphonium-based room-temperature ionic liquids. The study presents a compelling correlation between structural changes in water-ionic liquid solutions and thermodynamic and transport properties across diverse systems. The results for phosphonium ionic liquids are compared with 1-butyl-3-methylimidazolium hexaphosphate ([bmim]PF6) as a reference. Through this approach, phosphonium cation structure-related characteristics can be identified and placed within the broader context of ionic liquids. These insights are underpinned by observed changes in interaction energy, boiling point, diffusion rate, and viscosity, highlighting the crucial role of water molecules in weakening the strength of interactions between ions within the ionic liquid. The investigation also explains temperature-dependent trends in phosphonium cations, showing that alkyl group length and molecular symmetry are important tuning parameters for the strength of Coulomb interactions. These results contribute to a refined understanding of phosphonium ionic liquid behavior in the presence of water, offering valuable insights for optimizing their use in diverse fields.

How fork-length asymmetry affects solvent connectivity and diffusion in grafted polymeric model membranes

The hydrophilic pore morphology and solvent diffusion within model (amphiphilic) polymer membranes are simulated by dissipative particle dynamics (DPD). The polymers are composed of a backbone of 18 covalently bonded A beads to which at regular intervals side chains are attached. The side chains are composed of linear Ap chains (i.e., -A1-A2…Ap) from which two branches, [AsC] and [ArC], split off (s ≤ r). C beads serve as functionalized hydrophilic pendent sites. The branch lengths (s + 1 and r + 1) are varied. Five repeat unit designs (with general formula A3[Ap[AsC][ArC]]) are considered: A2[A3C][A3C] (symmetric branching), A2[A2C][A4C], A2[AC][A5C], A2[C][A6C] (highly asymmetric branching), and A4[AC][A3C]. The distribution of water (W) and W diffusion through nanophase segregated hydrophilic pores is studied. For similar primary length p, an increase in side chain symmetry favors hydrophilic pore connectivity and long-range water transport. C beads located on the longer [ArC] branches reveal the highest C bead mobility and are more strongly associated with water than the C beads on the shorter [AsC] branches. The connectivity of hydrophilic (W and W + C) phases through mapped replica of selected snapshots obtained from Monte Carlo tracer diffusion simulations is in line with trends found from the W bead diffusivities during DPD simulations. The diffusive pathways for protons (H+) in proton exchange membranes and for hydronium (OH-) in anion exchange membranes are the same as for solvents. Therefore, control of the side chain architecture is an interesting design parameter for optimizing membrane conductivities.

Miniaturized Capsule System Toward Real-Time Electrochemical Detection of H2 S in the Gastrointestinal Tract

Hydrogen sulfide (H2 S) is a gaseous inflammatory mediator and important signaling molecule for maintaining gastrointestinal (GI) homeostasis. Excess intraluminal H2 S in the GI tract has been implicated in inflammatory bowel disease and neurodegenerative disorders; however, the role of H2 S in disease pathogenesis and progression is unclear. Herein, an electrochemical gas-sensing ingestible capsule is developed to enable real-time, wireless amperometric measurement of H2 S in GI conditions. A gold (Au) three-electrode sensor is modified with a Nafion solid-polymer electrolyte (Nafion-Au) to enhance selectivity toward H2 S in humid environments. The Nafion-Au sensor-integrated capsule shows a linear current response in H2 S concentration ranging from 0.21 to 4.5 ppm (R2 = 0.954) with a normalized sensitivity of 12.4% ppm-1 when evaluated in a benchtop setting. The sensor proves highly selective toward H2 S in the presence of known interferent gases, such as hydrogen (H2 ), with a selectivity ratio of H2 S:H2 = 1340, as well as toward methane (CH4 ) and carbon dioxide (CO2 ). The packaged capsule demonstrates reliable wireless communication through abdominal tissue analogues, comparable to GI dielectric properties. Also, an assessment of sensor drift and threshold-based notification is investigated, showing potential for in vivo application. Thus, the developed H2 S capsule platform provides an analytical tool to uncover the complex biology-modulating effects of intraluminal H2 S.

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.

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. The effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

Simulated and Experimental Trends Regarding Water Uptake in Polymeric Electrolyte Membranes

Polymeric membranes in an anion or a proton exchange membrane fuel cell need sufficient hydration in order to provide a high hydroxide ion or proton conductivity. The water uptake for six model ionomer membranes, all of the same ion exchange capacity, is modeled by dissipative particle dynamics. The architectures cover three types of families that are of potential interest in fuel cell membrane research. All architectures consist of connected hydrophobic backbone A beads, to which side chains are grafted. For the type I family, the hydrophilic (functional) C beads are pendent on (amphiphilic) [AxC] side chains. The type II architecture contains both hydrophobic [A4] and short hydrophilic [C] side chains. For type III, the C beads are embedded along various locations within the [AxCAy] side chains (x + y = constant). For similar equilibrium time, the membrane water volume fraction increases with side chain length x for type I, and for type III, it increases with the distance x that C beads are separated from the backbone. Among the architectures (types I and III) for which the number of covalent C-A bonds are the same, the water uptake increases with the average number of A-A and A-C bonds (dpd springs) between A beads and the nearest C bead. A picture emerges in which for similar ion exchange capacity model membranes water uptake increases as a function of ⟨Nbondphob-phyl⟩.

Novel side chain functionalized polystyrene/O-PBI blends with high alkaline stability for anion exchange membrane water electrolysis (AEMWE)

We report the synthesis of a polystyrene-based anion exchange polymer bearing the cationic charge at a C6-spacer. The polymer is prepared by a functionalized monomer strategy. First, a copper halide catalyzed C-C coupling reaction between a styryl Grignard and 1,6-dibromohexane is applied, followed by quaternization with N-methylpiperidine and free radical polymerization. The novel polymer is blended with the polybenzimidazole O-PBI to yield mechanically stable blend membranes representing a new class of anion exchange membranes. In this regard, the ratio of the novel anion exchange polymer to O-PBI is varied to study the influence on water uptake and ionic conductivity. Blend membranes with IECs between 1.58 meq. OH- g-1 and 2.20 meq. OH- g-1 are prepared. The latter shows excellent performance in AEMWE, reaching 2.0 A cm-2 below 1.8 V in 1 M KOH at 70 °C, with a minor degradation rate from the start. The blend membranes show no conductivity loss after immersion in 1 M KOH at 85 °C for six weeks indicating high alkaline stability.

Perspective: Morphology and ion transport in ion-containing polymers from multiscale modeling and simulations

Ion-containing polymers are soft materials composed of polymeric chains and mobile ions. Over the past several decades they have been the focus of considerable research and development for their use as the electrolyte in energy conversion and storage devices. Recent and significant results obtained from multiscale simulations and modeling for proton exchange membranes (PEMs), anion exchange membranes (AEMs), and polymerized ionic liquids (polyILs) are reviewed. The interplay of morphology and ion transport is emphasized. We discuss the influences of polymer architecture, tethered ionic groups, rigidity of the backbone, solvents, and additives on both morphology and ion transport in terms of specific interactions. Novel design strategies are highlighted including precisely controlling molecular conformations to design highly ordered morphologies; tuning the solvation structure of hydronium or hydroxide ions in hydrated ion exchange membranes; turning negative ion-ion correlations to positive correlations to improve ionic conductivity in polyILs; and balancing the strength of noncovalent interactions. The design of single-ion conductors, well-defined supramolecular architectures with enhanced one-dimensional ion transport, and the understanding of the hierarchy of the specific interactions continue as challenges but promising goals for future research.

Translational Jump-Diffusion of Hydroxide Ion in Anion Exchange Membrane: Deciphering the Nature of Vehicular Diffusion

Earlier, ab initio and reactive force-field-based molecular dynamics (MD) simulation studies suggested an overwhelming contribution of the vehicular diffusion in the total diffusion of hydroxide ions rather than structural diffusion. But does the vehicular diffusion occur via small-step displacement? This question is important to have an understanding of the real characteristics of vehicular diffusion. To answer this question, we perform a classical molecular dynamics simulation of a system containing a hydroxide ion exchange membrane polymer and hydroxide ion at different hydration levels and temperatures using the same molecular force field (Dubey, V. Chem. Phys. Lett. 2020, 755, 137802), which successfully captured the microscopic structure and dynamics of the system. We use the translational jump-diffusion approach, used previously in supercooled water for understanding the origin of breakdown of the Stokes-Einstein relation, to calculate the jump-diffusion coefficient of hydroxide ion and water in the anion exchange membrane. We have seen a significant role of hydration level and temperature in the mechanism of vehicular diffusion. In overhydrated membrane, both hydroxide ions and water molecules diffuse via both small- and large-step displacement. With decreasing hydration level and temperature, the diffusion is increasingly governed by the jump-diffusion mechanism. The larger contribution of jump-diffusion comes from the stronger caging of the diffusing species by the solvent at lower hydration levels and temperature. These results, therefore, suggest that the hydration level and temperature of the hydroxide ion exchange membrane determine the detailed mechanism of the vehicular diffusion of hydroxide ion, especially whether the diffusion follows hydrodynamics or not.

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.

A broad-range variable-temperature solid state NMR spectral and relaxation investigation of the water state in Nafion 117

Understanding the water state in Nafion is not only crucial for operating a proton-exchange membrane (PEM)-based fuel cell, but also intimately related to the elucidation of the proton transport mechanism in a PEM. Although many studies have been published on this subject, some controversies and ambiguities remain unresolved. In this work, we design three different types of Nafion samples by substituting protons with lithium or sodium cations. We also pay special attention to the preparation of samples for carrying out broad-range variable temperature solid state NMR experiments so that no membrane dehydration occurs during the long experimental time at low temperatures. With these precautions and improvements, clear and largely straightforward information could be obtained to ensure minimal ambiguity and complexity in the interpretation of the experimental data. Our results show that about 40-60% of water remains unfrozen at -70 °C, depending on the type of the substituting cation. Both the ¹H and ²H spectral and relaxation results indicate that water freezing starts from the center of the nanopores inside Nafion and increases gradually as the temperature decreases. The protons remain dissociated with sulfonate groups even at the lowest temperature we reached (-70 °C), whereas both lithium and sodium are associated with sulfonate groups at most temperatures below 0 °C. The experimental data also suggest that besides frozen and unfrozen water, there is broad distribution of water state and dynamics in Nafion as the temperature is lowered from above zero down to -70 °C. The effect of the size of the substituting cation significantly affects the properties of supercooled water by modifying the cation-water interaction and impeding the rotation of sulfonate groups. These novel results not only help us in establishing a better understanding of the water state in Nafion and its performance as a proton exchange mebrane, but also provide insights into water freezing, antifreeze and supercooling in other nanoscopic environments.

Exclusion Zone Phenomena in Water-A Critical Review of Experimental Findings and Theories

The existence of the exclusion zone (EZ), a layer of water in which plastic microspheres are repelled from hydrophilic surfaces, has now been independently demonstrated by several groups. A better understanding of the mechanisms which generate EZs would help with understanding the possible importance of EZs in biology and in engineering applications such as filtration and microfluidics. Here we review the experimental evidence for EZ phenomena in water and the major theories that have been proposed. We review experimental results from birefringence, neutron radiography, nuclear magnetic resonance, and other studies. Pollack theorizes that water in the EZ exists has a different structure than bulk water, and that this accounts for the EZ. We present several alternative explanations for EZs and argue that Schurr's theory based on diffusiophoresis presents a compelling alternative explanation for the core EZ phenomenon. Among other things, Schurr's theory makes predictions about the growth of the EZ with time which have been confirmed by Florea et al. and others. We also touch on several possible confounding factors that make experimentation on EZs difficult, such as charged surface groups, dissolved solutes, and adsorbed nanobubbles.

Dissipative Particle Dynamics Modeling of Polyelectrolyte Membrane-Water Interfaces

Previous experiments of water vapor penetration into polyelectrolyte membrane (PEM) thin films have indicated the influence of the water concentration gradient and polymer chemistry on the interface evolution, which will eventually affect the efficiency of the fuel cell operation. Moreover, PEMs of different side chains have shown differences in water cluster structure and diffusion. The evolution of the interface between water and polyelectrolyte membranes (PEMs), which are used in fuel cells and flow batteries, of three different side-chain lengths has been studied using dissipative particle dynamics (DPD) simulations. Higher and faster water uptake is usually beneficial in the operation of fuel cells and flow batteries. The simulated water uptake increased with the increasing side chain length. In addition, the water uptake was rapid initially and slowed down afterwards, which is in agreement with the experimental observations. The water cluster formation rate was also found to increase with the increasing side-chain length, whereas the water cluster shapes were unaffected. Water diffusion in the membranes, which affects proton mobility in the PEMs, increased with the side-chain length at all distances from the interface. In conclusion, side-chain length was found to have a strong influence on the interface water structure and water penetration rates, which can be harnessed for the better design of PEMs. Since the PEM can undergo cycles of dehydration and rehydration, faster water uptake increases the efficiency of these devices. We show that the longer side chains with backbone structure similar to Nafion should be more suitable for fuel cell/flow battery usage.

Green Synthetic Fuels: Renewable Routes for the Conversion of Non-Fossil Feedstocks into Gaseous Fuels and Their End Uses

Innovative renewable routes are potentially able to sustain the transition to a decarbonized energy economy. Green synthetic fuels, including hydrogen and natural gas, are considered viable alternatives to fossil fuels. Indeed, they play a fundamental role in those sectors that are difficult to electrify (e.g., road mobility or high-heat industrial processes), are capable of mitigating problems related to flexibility and instantaneous balance of the electric grid, are suitable for large-size and long-term storage and can be transported through the gas network. This article is an overview of the overall supply chain, including production, transport, storage and end uses. Available fuel conversion technologies use renewable energy for the catalytic conversion of non-fossil feedstocks into hydrogen and syngas. We will show how relevant technologies involve thermochemical, electrochemical and photochemical processes. The syngas quality can be improved by catalytic CO and CO2 methanation reactions for the generation of synthetic natural gas. Finally, the produced gaseous fuels could follow several pathways for transport and lead to different final uses. Therefore, storage alternatives and gas interchangeability requirements for the safe injection of green fuels in the natural gas network and fuel cells are outlined. Nevertheless, the effects of gas quality on combustion emissions and safety are considered.

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.

Elucidating the morphological aspects and proton dynamics in a hybrid perfluorosulfonic acid membrane for medium-temperature fuel cell applications

A perfluorosulfonic acid (PFSA) membrane, i.e. Nafion® 117, was doped with the heteropoly salts (HPS) Cs3PW12O40, Rb3PW12O40, and (NH4)3PW12O40. Also, composite membranes with CsxH3-xPW12O40 (x = 1, 2, and 3) as dopants were investigated, which were rendered insoluble by substituting protons with larger cations. Morphological assessment and a detailed analysis of the hopping events via SAXS measurement and analysis of the hydrogen bond networks were performed using classical and quantum hopping molecular dynamics simulation. The phase segregation decreased by increasing the extent proton substitution in HPA. HPS containing cations with a larger ionic radius induced smaller phase segregation in the membrane, as confirmed by the RDF plots. SAXS simulation revealed that the hydrophilic phase domains in the HPS-doped Nafion® membrane were spaced further apart than that in the HPA-doped membrane. Although there was a greater number of isolated clusters for the Cs3PW12O40-doped Nafion®, the average number of cluster decreased with an increase in the substitution cation/proton ratio and ionic radius of the cation. The analysis of the H-bond network stability revealed that the proton hops slower when the membrane contains HPS particles and the mean residence time of a proton on water molecules increases with an increase in the extent of proton substitution in H3PW12O40. Indeed, for the HPA-doped membrane, the diffusion of water molecules is lower than that in the HPS-doped system.

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

Nafion nanocomposites for energy-related applications are being used extensively because of the attractive properties such as enhanced water retention, low unwanted crossover of electrolytes, and high proton conductivity. We present the results of the molecular dynamics modeling of Nafion films confined between two walls (substrates) of different polymer-wall interaction strengths and of different separation distances to model Nafion nanocomposites. Our goal is to provide insights into the effects of varying hydrophilicity and volume fraction of fillers/nanoparticles on the internal structure and water transport inside the Nafion membrane. The sulfur-sulfur radial distribution function first peak distance and the sulfur-oxygen (water) coordination number in the first hydration shell were negligibly affected by the wall (substrate) hydrophilicity or the film thickness. The Nafion side chains were found to bend toward the substrates with high hydrophilicity which is in qualitative agreement with existing experiments. The amount of bending was observed to reduce with increasing film thickness. However, the side-chain length did not show any noticeable variation with wall (substrate) hydrophilicity or film thickness. The water clusters became smaller and more isolated clusters emerged for highly hydrophilic substrates. In addition, the water cluster sizes showed a decreasing trend with decreasing film thickness in the case of hydrophilic substrates, which has also been observed in experiments of supported Nafion films. The in-plane water diffusion was enhanced considerably for hydrophilic substrates, and this mechanism has also been proposed previously in experiments. The in-plane water diffusion was also found to be a strong function of the substrate selectivity toward the hydrophilic phase. Our simulations can help provide more insights to experimentalists for choosing or modifying nanoparticles for Nafion nanocomposites.

Theoretical analyses on water cluster structures in polymer electrolyte membrane by using dissipative particle dynamics simulations with fragment molecular orbital based effective parameters

The mesoscopic structures of polymer electrolyte membrane (PEM) affect the performances of fuel cells. Nafion® with the Teflon® backbone has been the most widely used of all PEMs, but sulfonated poly-ether ether-ketone (SPEEK) having an aromatic backbone has drawn interest as an alternative to Nafion. In the present study, a series of dissipative particle dynamics (DPD) simulations were performed to compare Nafion and SPEEK. These PEM polymers were modeled by connected particles corresponding to the hydrophobic backbone and the hydrophilic moiety of sulfonic acid group. The water particle interacting with Nafion particles was prepared as well. The crucial interaction parameters among DPD particles were evaluated by a series of calculations based on the fragment molecular orbital (FMO) method in a non-empirical way (Okuwaki et al., J. Phys. Chem. B, 2018, 122, 338–347). Through the DPD simulations, the water and hydrophilic particles aggregated, forming cluster networks surrounded by the hydrophobic phase. The structural features of formed water clusters were investigated in detail. Furthermore, the differences in percolation behaviors between Nafion and SPEEK revealed much better connectivity among water clusters by Nafion. The present FMO-DPD simulation results were in good agreement with available experimental data.

The mesoscopic structures of polymer electrolyte membrane (PEM) affect the performances of fuel cells.

Classical Molecular Dynamics with Mobile Protons

An important limitation of standard classical molecular dynamics simulations is the inability to make or break chemical bonds. This restricts severely our ability to study processes that involve even the simplest of chemical reactions, the transfer of a proton. Existing approaches for allowing proton transfer in the context of classical mechanics are rather cumbersome and have not achieved widespread use and routine status. Here we reconsider the combination of molecular dynamics with periodic stochastic proton hops. To ensure computational efficiency, we propose a non-Boltzmann acceptance criterion that is heuristically adjusted to maintain the correct or desirable thermodynamic equilibria between different protonation states and proton transfer rates. Parameters are proposed for hydronium, Asp, Glu, and His. The algorithm is implemented in the program CHARMM and tested on proton diffusion in bulk water and carbon nanotubes and on proton conductance in the gramicidin A channel. Using hopping parameters determined from proton diffusion in bulk water, the model reproduces the enhanced proton diffusivity in carbon nanotubes and gives a reasonable estimate of the proton conductance in gramicidin A.

Molecular Structuring and Percolation Transition in Hydrated Sulfonated Poly(ether ether ketone) Membranes

The extent of phase separation and water percolation in sulfonated membranes are the key to their performance in fuel cells. Toward this, the effect of hydration on the morphology and transport characteristics of sulfonated poly(ether ether ketone), sPEEK, membrane is investigated using atomistic molecular dynamics simulation at various hydration levels(λ: number of water molecules per sulfonate group). The evolution of local morphology is investigated using structural correlations and minimum pair distances. Transport properties are probed using mean squared displacements and diffusion coefficients. The water-sulfonate interaction in sPEEK is found to be stronger than that in Nafion, as observed in experiments. Analyses indicate the presence of narrow connected path of water and hydronium at λ = 4 and large domains, spanning half the simulation box, at λ = 15. The behavior of membrane water remains far from bulk as indicated by its diffusion coefficient. The persistence of small isolated water clusters demonstrates the extent of phase separation in sPEEK to be lesser than that in Nafion. Analyses at molecular and collective levels suggest the occurrence of a percolation transition between λ = 8 and 10, which leads to a connected network of water channels in the membrane, thereby boosting the hydronium mobility.

Proton Transfer in Perfluorosulfonic Acid Fuel Cell Membranes with Differing Pendant Chains and Equivalent Weights

Proton transfer in the nanoscopic water channels of polyelectrolyte fuel cell membranes was studied using a photoacid, 8-hydroxypyrene-1,3,6-trisulfonic acid sodium salt (HPTS), in the channels. The local environment of the probe was determined using 8-methoxypyrene-1,3,6-trisulfonic acid sodium salt (MPTS), which is not a photoacid. Three fully hydrated membranes, Nafion (DuPont) and two 3M membranes, were studied to determine the impact of different pendant chains and equivalent weights on proton transfer. Fluorescence anisotropy and excited state population decay data that characterize the local environment of the fluorescent probes and proton transfer dynamics were measured. The MPTS lifetime and anisotropy results show that most of the fluorescent probes have a bulk-like water environment with a relatively small fraction interacting with the channel wall. Measurements of the HPTS protonated and deprotonated fluorescent bands' population decays provided information on the proton transport dynamics. The decay of the protonated band from ∼0.5 ns to tens of nanoseconds is in part determined by dissociation and recombination with the HPTS, providing information on the ability of protons to move in the channels. The dissociation and recombination is manifested as a power law component in the protonated band fluorescence decay. The results show that equivalent weight differences between two 3M membranes resulted in a small difference in proton transfer. However, differences in pendant chain structure did significantly influence the proton transfer ability, with the 3M membranes displaying more facile transfer than Nafion.

Parameterization of a coarse-grained model with short-ranged interactions for modeling fuel cell membranes with controlled water uptake

The coarse-grained model FF_pvap reproduces the experimental activity coefficient of water in tetramethylammonium chloride solutions over a wide range of concentrations, with a hundred-fold gain in computing efficiency with respect to atomistic models.

Molecular Dynamics Simulations Reveal that Water Diffusion between Graphene Oxide Layers is Slow

Membranes made of stacked layers of graphene oxide (GO) hold the tantalizing promise of revolutionizing desalination and water filtration if selective transport of molecules can be controlled. We present the findings of an integrated study that combines experiment and molecular dynamics simulation of water intercalated between GO layers. We simulated a range of hydration levels from 1 wt.% to 23.3 wt.% water. The interlayer spacing increased upon hydration from 0.8 nm to 1.1 nm. We also synthesized GO membranes that showed an increase in layer spacing from about 0.7 nm to 0.8 nm and an increase in mass of about 15% on hydration. Water diffusion through GO layers is an order of magnitude slower than that in bulk water, because of strong hydrogen bonded interactions. Most of the water molecules are bound to OH groups even at the highest hydration level. We observed large water clusters that could span graphitic regions, oxidized regions and holes that have been experimentally observed in GO. Slow interlayer diffusion can be consistent with experimentally observed water transport in GO if holes lead to a shorter path length than previously assumed and sorption serves as a key rate-limiting step.

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.

Elucidating Water Transport Mechanisms in Nafion Thin Films

Ion-exchange membranes are critical components of hydrogen fuel cells, where these ionomers can be confined to nanoscale thicknesses, altering the physical properties of these films from that of bulk membranes. Therefore, it is important to develop methods capable of measuring and elucidating the transport mechanisms under thin film confinement compared to bulk Nafion. In this study, water sorption and diffusion in a Nafion thin film were measured using time-resolved in situ polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS). Interfacial mass transport limitations were confirmed to be minimal, while restricted water diffusion was observed, where the effective diffusion coefficient of water in the thin Nafion film was many orders of magnitude lower (between 4 and 5 orders of magnitude) than those reported for bulk membranes and was dependent on the initial hydration state of the Nafion. Furthermore, the response of the hydrophobic domains (Teflon backbone) to the swelling of the hydrophilic domains (ionic clusters) was shown to be orders of magnitude slower than that of bulk Nafion.