PNNL

Abstract

We report the results of ab initio molecular dynamics simulations of a model NafionTM polymer membrane initially equilibrated using classical molecular dynamics simulations. We studied three hydration levels (?) of 3, 9, and 15 H2O/SO3- corresponding to dry, hydrated and saturated fuel cell membrane, respectively. The barrier for proton transfer from the SO3-–H3O+ contact ion pair to a solvent-separated ion pair decreased from 2.3 kcal/mol for ? = 3 to 0.8 kcal/mol for ? = 15. The barrier for proton transfer between two water molecules was in the range from 0.7 to 0.8 kcal/mol for the ? values studied. The number of proton shuttling events between a pair of water molecules is an order of magnitude more than the number of proton hops across three distinct water molecules. The proton diffusion coefficient at ? = 15 is about 0.9x10-5 cm2/s, which is in good agreement with experiment and our previous quantum hopping molecular dynamics simulations.