Building on the previously developed multistate empirical valence bond model [U. W. Schmitt and G. A. Voth, J. Chem. Phys 111, 9361 (1999)] for the dynamics and energetics of an excess proton in bulk phase water, a second generation model is described. This model is shown to produce similar dynamic and structural properties to the previous model, while allowing for the use of the full hydronium charge. This characteristic of the model is required for its implementation in a host of realistic applications beyond bulk water. An improved state selection algorithm is also presented, resulting in a significantly reduced energy drift during microcanonical molecular dynamics simulations. The unusually high self diffusion constant of an excess proton in water due to the proton hopping (Grotthuss) process is observed in the simulation data and is found to be quantitatively in the same range as the experimental value if a quantum correction is taken into consideration. Importantly, a more complete analysis of proton transport process is also presented.