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 simulation of a dichain surfactant + water + carbon dioxide (solvent) system is performed to study the structural properties of reversed micelle-like surfactant aggregates formed in the system. The simulations use a detailed and realistic molecular model for the surfactant molecule and explicit representation of the water and solvent molecules to enable quantitative comparisons with a prior experimental (small-angle neutron scattering) study. The results of the simulation are found to be in reasonable agreement with experimental values. The simulations show that the size and shape of the surfactant aggregates depends on their water-to-surfactant ratio. A higher water-to-surfactant ratio results in larger and more spherical aggregates. The two distinct tails of the surfactant molecule exhibit different conformations in carbon dioxide indicating contrasting CO2-philic behavior. The perfluoroalkane tails assume more extended conformation than the alkane tails. The microstructure of the aqueous core reveals that the water molecules in the interfacial region are strongly oriented in response to the electric fields of the anionic headgroups and sodium counterions, while water near the center of the core approaches bulklike properties with the presence of a hydrogen-bonded network.