“…Yan et al (Yan et al, 2008) have studied electroosmotic drag using the classical hydronium cation model in the MD simulations and predicted a value of K drag about 2.5 times greater than the experimental data at 300 K in completely hydrated Nafion membranes. Although few studies incorporate the Grotthuss mechanism into MD simulations to develop reactive MD using, for example, empirical valence bond (EVB) approaches (Feng et al, 2012, Feng and Voth, 2011, Jorn et al, 2012, Petersen and Voth, 2006, Petersen et al, 2005, Tse et al, 2013, Seeliger et al, 2005, Spohr et al, 2002, Savage et al, 2014, they have focused only on the proton diffusion in the equilibrium polymer systems, and therefore the electroosmotic properties have not been fully understood.…”