Simulation of equivalent weight dependence of Nafion morphologies and predicted trends regarding water diffusion

Microphase separation within hydrated Nafion® membranes was simulated using Dissipative Particle Dynamics (DPD). Morphologies were obtained at branching densities corresponding with equivalent weights ranging from ∼800 to ∼1400 (g/mole SO3) and water percentage volume contents ([H2O]) varying between 10% and 30%. All cases showed pronounced microphase separation involving a hydrophobic Teflon phase and a hydrophilic phase in which water is associated with SO3 groups that are located near the phase boundaries. Pore morphologies were found to depend strongly on water content and branching density. The average pore radius (Rpore) and the distance between the pores (Dcl-cl) were found to increase with water content obeying the relations Rpore = 1.3 + α[H2O] (nm), and Dcl-cl = 3.2 + β[H2O] (nm). The values of the expansion coefficients α and β decrease linearly with branching density with α = 5.3 × 10−5 × (EW-450) and β = 1.3 × 10−4 × (EW-450) nm/vol%. For decreasing branching density the pores obtain a more spherical character. The consequence of this on water diffusion is estimated by employing Monte Carlo trajectory calculations in which we assume that water movement is confined within the hydrophilic phase and local water mobility to be equal to that of pure water. The estimated diffusion constants increase linearly with branching density (i.e. linear decrease with equivalent weight). Experimental water diffusion constants obtained from literature for Nafion1100 membrane are in good agreement with our calculations. A counterintuitive picture evolves in which smaller pores lead to enhanced water diffusion.