A molecular dynamics simulation study of hydrated Nafion at water contents ranging from 5 to 20 wt % was performed to examine the structure and dynamics of the hydrated polyelectrolyte system. The simulations show that the system forms segregated hydrophobic regions consisting primarily of the polymer backbone and hydrophilic regions with an inhomogeneous water distribution. We find that the water clustering strongly depends on the water content. At low water content, only isolated small water clusters are formed. As the water content increases, it becomes increasingly possible that a predominant majority of water molecules form a single cluster, suggesting that the hydrophilic regions become connected. We characterize the atomic structures formed within the system by various atomic pair correlation functions. The water structure factor shows a peak at q values corresponding to an intercluster distance about 2.5 nm and greater. With increasing water content, the distance moves to larger values, consistent with findings from scattering experiments. We find that the degree of solvation of hydronium ions by water molecules is a strong function of water content. At 5 wt %, a majority of the hydronium ions are hydrated by no more than two water molecules, prohibiting structural diffusion. As water content increases, the hydronium ions continue to become increasingly hydrated, resulting in structures capable of forming eigen ions, a necessary step in structural diffusion. Addressing the experimentally observed fact that conductivity in these membranes abruptly drops near 5 wt %, we find that both the local structure of the poorly hydrated hydronium ions and the disconnected nature of the global morphology of the water nanonetwork at low water content should contribute to poor conductivity.

Molecular dynamics simulations were performed to investigate the relationship between the molecular structure of perfluorosulfonic acid (PFSA) ionomers and the nanoscale morphology of the hydrated membranes. Three structural features are examined including (i) the length of the side chain to which the sulfonic acid group is attached, (ii) the equivalent weight (EW) of the electrolyte ionomer, and (iii) the molecular weight (MW) of the polymer electrolyte. Membrane morphologies are studied from the water content λ = 3 (λ represents number of water molecules per sulfonate group) to saturation (λ = 22). We find that with the longer side chain, there is more clustering of the sulfonate groups and more local water−water clustering, but a more poorly connected aqueous domain. When one decreases the equivalent weight (EW) in either the short side chain (SSC) PFSA or Nafion, there is more clustering of the sulfonate groups and more local water−water clustering and a better connected aqueous domain. Because connectivity enhances and confinement reduces water mobility, a decrease in EW, which enhances connectivity and reduces confinement, results in an increase in diffusivity. An increase in side chain length diminishes connectivity but reduces confinement, which together result in little change in the observed water diffusivity. For the short chains studied, we find these results to be independent of MW.

Classical molecular dynamics (MD) simulations were performed to determine the hydrated morphology and hydronium ion diffusion coefficients in two different perfluorosulfonic acid (PFSA) membranes as functions of water content. The structural and transport properties of 1143 equivalent weight (EW) Nafion, with its relatively long perfluoroether side chains, are compared to the short-side-chain (SSC) PFSA ionomer at an EW of 977. The separation of the side chains was kept uniform in both ionomers consisting of -(CF 2) 15- units in the backbone, and the degree of hydration was varied from 5 to 20 weight % water. The MD simulations indicated that the distribution of water clusters is more dispersed in the SSC ionomer, which leads to a more connected water-channel network at the low water contents. This suggests that the SSC ionomer may be more inclined to form sample-spanning aqueous domains through which transport of water and protons may occur. The diffusion coefficients for both hydronium ions and water molecules were calculated at hydration levels of 4.4, 6.4, 9.6, and 12.8 H 2O/SO 3H for each ionomer. When compared to experimental proton diffusion coefficients, this suggests that as the water content is increased the contribution of proton hopping to the overall proton diffusion increases.

Molecular dynamics (MD) simulations were performed to investigate the structural and dynamical behavior of water and hydronium ions at the electrode/electrolyte interface of hydrogen polymer electrolyte membrane (PEM) fuel cells. Specifically, we have studied the hydrated Nafion membrane, humidified for four different water contents, 5, 10, 15, and 20%, at 300 K. We analyzed the three-phase interface where the hydrated PEM is in contact with the vapor phase and with either the catalyst surface (platinum in this paper) or the catalyst−support surface (graphite in this paper). These molecular simulations represent portions of interfaces that exist within the PEM fuel cells. We observed significant wetting of the catalyst surface by a mixture of polymer, water, and hydronium ions but not beyond a monolayer. We observed virtually no wetting of the graphite surface. On the catalyst surface, the degree of wetting of the catalyst surface depends strongly on the level of membrane humidity. The pair correlation functions indicate that the water molecules adsorbed in a monolayer on the catalyst surface form small domains of ordered structures, which are bound by fragments of Nafion on the surface. The diffusion of protons from the catalyst surface into the membrane must proceed across this highly inhomogeneous surface.