Hydrated excess proton at water-hydrophobic interfaces

The behavior of the hydrated excess proton is investigated at the water-vapor, water-hydrophobic wall, and water-carbon tetrachloride interfaces through molecular dynamics simulations using the third-generation multistate empirical valence bond model (MS-EVB3). The MS-EVB3 simulations show a surface preference of the excess proton at the water-vapor interface, consistent with the discovery of this effect using an earlier version of the MS-EVB model (Petersen et al. J. Phys. Chem. B 2004, 108, 14804) and with the experimental results. The preference of the hydrated excess proton for other water-hydrophobic interfaces is also analyzed for the first time. The hydrated proton structures and charge defect delocalization effects at these interfaces are discussed in detail. By decomposing the free energy profiles into the internal energy and entropic contributions, the thermodynamic (free energy) driving forces for the surface preference of the excess proton are also elaborated. These results indicate that the "rigid" hydrated proton structures at all the interfaces are energetically (as opposed to entropically) stabilized due to the "amphiphilic" nature of the hydrated excess proton, resulting in its overall interfacial concentration enhancement. The effects of acid concentration and nuclear quantization are also explored.     

A molecular dynamics study of water nucleation using the TIP4P/2005 model

Extensive molecular dynamics simulations were conducted using the TIP4P/2005 water model of Abascal and Vega [J. Chem. Phys. 123, 234505 (2005)] to investigate its condensation from supersaturated vapor to liquid at 330 K. The mean first passage time method [J. Wedekind, R. Strey, and D. Reguera, J. Chem. Phys. 126, 134103 (2007); L. S. Bartell and D. T. Wu, 125, 194503 (2006)] was used to analyze the influence of finite size effects, thermostats, and charged species on the nucleation dynamics. We find that the Nose?-Hoover thermostat and the one proposed by Bussi et al. [J. Chem. Phys. 126, 014101 (2007)] give essentially the same averages. We identify the maximum thermostat coupling time to guarantee proper thermostating for these simulations. The presence of charged species has a dramatic impact on the dynamics, inducing a marked change towards a pure growth regime, which highlights the importance of ions in the formation of liquid droplets in the atmosphere. It was found a small but noticeable sign preference at intermediate cluster sizes (between 5 and 30 water molecules) corresponding mostly to the formation of the second solvation shell around the ion. The TIP4P/2005 water model predicts that anions induce faster formation of water clusters than cations of the same magnitude of charge.     

Quantum effects in liquid water and ice: model dependence

This paper explores the influence of choice of potential model on the quantum effects observed in liquid water and ice. This study utilizes standard rigid models and a more formal context for the rigid-body centroid molecular dynamics methodology used to perform the quantum simulations is provided. Quantum and classical molecular dynamics simulations are carried out for liquid water and ice Ih at 298 and 220 K, respectively, with the simple point charge/extended and TIP4P-Ew water models. The results obtained for equilibrium and dynamical properties are compared with those recently reported on TIP4P [L. Hernandez de la Pena and P. G. Kusalik, J. Chem. Phys. 121, 5992 (2004); L. Hernandez de la Pena et al., J. Chem. Phys 123, 144506 (2005)]. For the liquid, an energy shift of about 8% and an average molecular uncertainty of about 11 degrees were found independently of the water model. The self-diffusion coefficient consistently increases by more than 50% when going from the classical to the quantum system and quantum dynamics are found to reproduce the experimental isotopic shifts with the models examined. The ice results compare remarkably well with those previously reported for the TIP4P water model; they confirm that quantum effects are considerable and that the quantum mechanical uncertainty and the energy shifts due to quantization are smaller in ice than in liquid water. The relevance of these findings in the context of the construction of water models is briefly discussed.     

Development of transferable interaction potentials for water. V. Extension of the flexible, polarizable, Thole-type model potential (TTM3-F, v. 3.0) to describe the vibrational spectra of water clusters and liquid water

We present a new parametrization of the flexible, polarizable Thole-type model for water [J. Chem. Phys. 116, 5115 (2002); J. Phys. Chem. A 110, 4100 (2006)], with emphasis in describing the vibrational spectra of both water clusters and liquid water. The new model is able to produce results of similar quality with the previous versions for the structures and energetics of water clusters as well as structural and thermodynamic properties of liquid water evaluated with classical and converged quantum statistical mechanical atomistic simulations. At the same time it yields accurate redshifts for the OH vibrational stretches of both water clusters and liquid water.     

Thermodynamical and structural properties of imidazolium based ionic liquids from molecular simulation

We have performed molecular dynamics simulations to determine the densities and heat of vaporization as well as structural information for the 1-alkyl-3-methyl-imidazolium based ionic liquids [amim][Cl] and [amim][BF(4)] in the temperature range from 298 to 363 K. In this simulation study, we used an united atom model of Liu et al. [Phys. Chem. Chem. Phys. 8, 1096 (2006)] for the [emim(+)] and [bmim(+)] cations, which we have extended for simulation in [hmim]-ILs and combined with parameters of Canongia Lopes et al. [J. Phys. Chem. B 108, 2038 (2004)] for the [Cl(-)] anion. Our simulation results prove that both the original united atoms approach by Liu et al. and our extension yield reasonable predictions for the ionic liquid with a considerably reduced computational expense than that required for all atoms models. Radial distribution functions and spatial distribution functions where employed to analyze the local structure of this ionic liquid, and in which way it is influenced by the type of the anion, the size of the cation, and the temperature. Our simulations give evidence for the occurrence of tail aggregations in these ionic liquids with increasing length of the side chain and also increasing temperature.     

Mechanistic aspects of proton chain transfer in the green fluorescent protein. Part II. A comparison of minimal quantum chemical models

In this paper we report the results of extensive quantum chemical reaction pathway calculations for the electronic ground state of several different cluster models that mimic the proton chain transfer path within the green fluorescent protein (GFP). Our principal objective is to establish the robustness with respect to variations in the model of our recent mechanistic inferences for the ground state proton chain transfer [S. Wang and S. C. Smith, J. Phys. Chem. B, 2006, 110, 5084]. Additionally, comparison of our ground state results with the excited state proton transfer (ESPT) study by Vendrell et al. [O. Vendrell, R. Gelabert, M. Moreno and J. M. Lluch, J. Am. Chem. Soc., 2006, 128, 3564] leads to the conclusion that the mechanism of proton chain transfer may be expected to be analogous in ground and excited states, principally because in both cases the loss of the chromophore's phenolic proton contributes strongly to the reaction coordinate only late in the reaction path.     

Aqueous solvation free energies of ions and ion-water clusters based on an accurate value for the absolute aqueous solvation free energy of the proton

Thermochemical cycles that involve pKa, gas-phase acidities, aqueous solvation free energies of neutral species, and gas-phase clustering free energies have been used with the cluster pair approximation to determine the absolute aqueous solvation free energy of the proton. The best value obtained in this work is in good agreement with the value reported by Tissandier et al. (Tissandier, M. D.; Cowen, K. A.; Feng, W. Y.; Gundlach, E.; Cohen, M. J.; Earhart, A. D.; Coe, J. V. J. Phys. Chem. A 1998, 102, 7787), who applied the cluster pair approximation to a less diverse and smaller data set of ions. We agree with previous workers who advocated the value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton. Considering the uncertainties associated with the experimental gas-phase free energies of ions that are required to use the cluster pair approximation as well as analyses of various subsets of data, we estimate an uncertainty for the absolute aqueous solvation free energy of the proton of no less than 2 kcal/mol. Using a value of -265.9 kcal/mol for the absolute aqueous solvation free energy of the proton, we expand and update our previous compilation of absolute aqueous solvation free energies; this new data set contains conventional and absolute aqueous solvation free energies for 121 unclustered ions (not including the proton) and 147 conventional and absolute aqueous solvation free energies for 51 clustered ions containing from 1 to 6 water molecules. When tested against the same set of ions that was recently used to develop the SM6 continuum solvation model, SM6 retains its previously determined high accuracy; indeed, in most cases the mean unsigned error improves when it is tested against the more accurate reference data.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Nature of the water molecules in the palisade layer of a triton X-100 micelle in the presence of added salts: A solvation dynamics study

The effect of added electrolytes on the nature of water molecules in the palisade layer of a Triton X-100 (TX-100) micelle has been investigated using solvation dynamics studies of C153 dye in the presence of different concentrations of NaCl, KCl, and CsCl salts. In all of the cases, the solvation dynamics is found to be biexponential in nature. It is seen that in the presence of added salts the solvation dynamics becomes slower. As previously reported (Charlton et al. J. Phys. Chem. B 2000, 104, 8327; Molina-Bolivar et al. J. Phys. Chem. B 2002, 106, 870), the presence of salt increases micellar hydration (and also size) for TX-100, mainly due to enhancement in the mechanically trapped water content in the palisade layer. Under normal circumstances, increased micellar hydration was expected to cause faster solvation dynamics (Kumbhakar et al. J. Phys. Chem. B 2004, 108, 19246), though in the present work, a reverse trend is in fact observed with the added salts. In accordance with solvation dynamics results, fluorescence anisotropy studies also indicate an increase in microviscosity for the palisade layer of the TX-100 micelle with the added salts. The present results have been rationalized assuming that the ions reside in the palisade layer, and due to the hydration of the ions, especially the cations, the water molecules in the palisade layer undergo a kind of clustering, causing the microviscosity to in fact increase rather than decrease as expected due to increased micellar hydration. A partial collapse of the surfactant chains due to their dehydration as caused by the hydration of the ions in the palisade layer may also add to the increase in microviscosity and the consequent retardation in relaxation dynamics in the presence of salts.     

Role of hydration in determining the structure and vibrational spectra of L-alanine and N-acetyl L-alanine N′-methylamide in aqueous solution: a combined theoretical and experimental approach

In this work we have utilized recent density functional theory Born-Oppenheimer molecular dynamics simulations to determine the first principles locations of the water molecules in the first solvation shell which are responsible for stabilizing the zwitterionic structure of L-alanine. Previous works have used chemical intuition or classical molecular dynamics simulations to position the water molecules. In addition, a complete shell of water molecules was not previously used, only the water molecules which were thought to be strongly interacting (H-bonded) with the zwitterionic species. In a previous work by Tajkhorshid et al. (J Phys Chem B 102:5899) the L-alanine zwitterion was stabilized by 4 water molecules, and a subsequent work by Frimand et al. (Chem Phys 255:165) the number was increased to 9 water molecules. Here we found that 20 water molecules are necessary to fully encapsulate the zwitterionic species when the molecule is embedded within a droplet of water, while 11 water molecules are necessary to encapsulate the polar region with the methyl group exposed to the surface, where it migrates during the MD simulation. Here we present our vibrational absorption, vibrational circular dichroism and Raman and Raman optical activity simulations, which we compare to the previous simulations and experimental results. In addition, we report new VA, VCD, Raman and ROA measurements for L-alanine in aqueous solution with the latest commercially available FTIR VA/VCD instrument (Biotools, Jupiter, FL, USA) and Raman/ROA instrument (Biotools). The signal to noise of the spectra of L-alanine measured with these new instruments is significantly better than the previously reported spectra. Finally we reinvestigate the causes for the stability of the P_π structure of the alanine dipeptide, also called N-acetyl-L-alanine N′-methylamide, in aqueous solution. Previously we utilized the B3LYP/6-31G* + Onsager continuum level of theory to investigate the stability of the NALANMA4WC Han et al. (J Phys Chem B 102:2587) Here we use the B3PW91 and B3LYP hybrid exchange correlation functionals, the aug-cc-pVDZ basis set and the PCM and CPCM (COSMO) continuum solvent models, in addition to the Onsager and no continuum solvent model. Here by the comparison of the VA, VCD, Raman and ROA spectra we can confirm the stability of the NALANMA4WC due to the strong hydrogen bonding between the four water molecules and the peptide polar groups. Hence we advocate the use of explicit water molecules and continuum solvent treatment for all future spectral simulations of amino acids, peptides and proteins in aqueous solution, as even the structure (conformer) present cannot always be found without this level of theory. Festschift in Honor of Philip J. Stephens’ 65th Birthday. During the proof stage of this article a very relevant article has been published by M. Losada and Y. Xu titled “Chirality transfer through hydrogen-bonding: Experimental and ab initio analyses of vibrational circular dichroism spectra of methyl lactate in water” in Phys Chem Chem Phys 2007, 9: 3127–3135. In that work they confirm that the effects of water are seen in the VCD spectra and hence it is fundamental to include explicit water molecules in modeling studies of the vibrational spectra of biomolecules in aqueous solution.     

Computation of methodology-independent single-ion solvation properties from molecular simulations. IV. Optimized Lennard-Jones interaction parameter sets for the alkali and halide ions in water

The raw single-ion solvation free energies computed from atomistic (explicit-solvent) simulations are extremely sensitive to the boundary conditions and treatment of electrostatic interactions used during these simulations. However, as shown recently [M. A. Kastenholz and P. H. Hu?nenberger, J. Chem. Phys. 124, 224501 (2006); M. M. Reif and P. H. Hu?nenberger, J. Chem. Phys. 134, 144103 (2010)], the application of appropriate correction terms permits to obtain methodology-independent results. The corrected values are then exclusively characteristic of the underlying molecular model including in particular the ion-solvent van der Waals interaction parameters, determining the effective ion size and the magnitude of its dispersion interactions. In the present study, the comparison of calculated (corrected) hydration free energies with experimental data (along with the consideration of ionic polarizabilities) is used to calibrate new sets of ion-solvent van der Waals (Lennard-Jones) interaction parameters for the alkali (Li(+), Na(+), K(+), Rb(+), Cs(+)) and halide (F(-), Cl(-), Br(-), I(-)) ions along with either the SPC or the SPC/E water models. The experimental dataset is defined by conventional single-ion hydration free energies [Tissandier et al., J. Phys. Chem. A 102, 7787 (1998); Fawcett, J. Phys. Chem. B 103, 11181] along with three plausible choices for the (experimentally elusive) value of the absolute (intrinsic) hydration free energy of the proton, namely, ΔG(hyd)(?)[H(+)] = -1100, -1075 or -1050 kJ mol(-1), resulting in three sets L, M, and H for the SPC water model and three sets L(E), M(E), and H(E) for the SPC/E water model (alternative sets can easily be interpolated to intermediate ΔG(hyd)(?)[H(+)] values). The residual sensitivity of the calculated (corrected) hydration free energies on the volume-pressure boundary conditions and on the effective ionic radius entering into the calculation of the correction terms is also evaluated and found to be very limited. Ultimately, it is expected that comparison with other experimental ionic properties (e.g., derivative single-ion solvation properties, as well as data concerning ionic crystals, melts, solutions at finite concentrations, or nonaqueous solutions) will permit to validate one specific set and thus, the associated ΔG(hyd)(?)[H(+)] value (atomistic consistency assumption). Preliminary results (first-peak positions in the ion-water radial distribution functions, partial molar volumes of ionic salts in water, and structural properties of ionic crystals) support a value of ΔG(hyd)(?)[H(+)] close to -1100 kJ·mol(-1).     

Overlapping resonances in the control of intramolecular vibrational redistribution

Coherent control of bound state processes via the interfering overlapping resonance scenario [Christopher et al., J. Chem. Phys. 123, 064313 (2006)] is developed to control intramolecular vibrational redistribution (IVR). The approach is applied to the flow of population between bonds in a model of chaotic OCS vibrational dynamics, showing the ability to significantly alter the extent and rate of IVR by varying quantum interference contributions.     

Adding salt to an aqueous solution of t-butanol: is hydrophobic association enhanced or reduced?

Recent neutron scattering experiments on aqueous salt solutions of amphiphilic t-butanol by Bowron and Finney [Phys. Rev. Lett. 89, 215508 (2002); J. Chem. Phys. 118, 8357 (2003)] suggest the formation of t-butanol pairs, bridged by a chloride ion via O-H...Cl- hydrogen bonds, leading to a reduced number of intermolecular hydrophobic butanol-butanol contacts. Here we present a joint experimental/theoretical study on the same system, using a combination of molecular dynamics (MD) simulations and nuclear magnetic relaxation measurements. Both MD simulation and experiment clearly support the more classical scenario of an enhanced number of hydrophobic contacts in the presence of salt, as it would be expected for purely hydrophobic solutes. [T. Ghosh et al., J. Phys. Chem. B 107, 612 (2003)]. Although our conclusions arrive at a structurally completely distinct scenario, the molecular dynamics simulation results are within the experimental error bars of the Bowron and Finney data.     

Recovering the Crooks equation for dynamical systems in the isothermal-isobaric ensemble: a strategy based on the equations of motion

The Crooks equation [Eq. (10) in J. Stat. Phys. 90, 1481 (1998)] relates the work done on a system during a nonequilibrium transformation to the free energy difference between the final and the initial state of the transformation. Recently, the authors have derived the Crooks equation for systems in the canonical ensemble thermostatted by the Nose-Hoover or Nose-Hoover chain method [P. Procacci et al., J. Chem. Phys. 125, 164101 (2006)]. That proof is essentially based on the fluctuation theorem by Evans and Searles [Adv. Phys. 51, 1529 (2002)] and on the equations of motion. Following an analogous approach, the authors derive here the Crooks equation in the context of molecular dynamics simulations of systems in the isothermal-isobaric (NPT) ensemble, whose dynamics is regulated by the Martyna-Tobias-Klein algorithm [J. Chem. Phys. 101, 4177 (1994)]. Their present derivation of the Crooks equation correlates to the demonstration of the Jarzynski identity for NPT systems recently proposed by Cuendet [J. Chem. Phys. 125, 144109 (2006)].     

A quantum mechanical/molecular dynamics/mean field study of acrolein in aqueous solution: analysis of H bonding and bulk effects on spectroscopic properties

A novel molecular dynamics methodology recently proposed by our group [Rega et al., Chem. Phys. Lett. 422, 367 (2006)], which is based on an integrated hybrid potential rooted in high level quantum mechanical methods using localized basis functions and nonperiodic boundary conditions, has been applied to study acrolein in aqueous solution. The solute structural rearrangement and its hydrogen-bonding pattern due to the interactions with water have been analyzed in some detail. Moreover, the solvent effects on the UV n-->pi* vertical transition and on the NMR 13C and 17O shielding constants of acrolein have been investigated theoretically by performing a posteriori quantum mechanical calculations on a statistically significant number of snapshots extracted from both gas-phase and aqueous solution simulations. Results show that such effective computational strategy can be successfully used to improve our understanding, at atomic level, of important spectroscopic observables.     

New evaluation of reconstructed spatial distribution function from radial distribution functions

Although three dimensional (3D) solvation structure is much more informative than one dimensional structure, its evaluation is difficult experimentally and theoretically. In our previous Communication [Yokogawa et al., J. Chem. Phys. 123, 211102 (2005)], we proposed a new method to present reconstructed spatial distribution function (RC-SDF) from a set of radial distribution functions (RDFs). In this article, we successfully extended the method more accurately with new basis sets. This new method was applied to two liquid solvation structures, methanol and dimethyl sulfoxide, as examples. Their RC-SDFs evaluated here clearly show that the former solvation structure is well defined while the latter one is broad, which agrees well with the SDFs calculated directly from molecular dynamics simulations. These results indicate that the method can reproduce well these 3D solvation structures in reasonable computational cost.     

Simulation of vibrational dephasing of I(2) in solid Kr using the semiclassical Liouville method

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.     

Soft sticky dipole-quadrupole-octupole potential energy function for liquid water: an approximate moment expansion

A new, efficient potential energy function for liquid water is presented here. The new model, which is referred here as the soft sticky dipole-quadrupole-octupole (SSDQO) model, describes a water molecule as a Lennard-Jones sphere with point dipole, quadrupole, and octupole moments. It is a single-point model and resembles the hard-sphere sticky dipole potential model for water by Bratko et al. [J. Chem. Phys. 83, 6367 (1985)] and the soft sticky dipole model by Ichiye and Liu [J. Phys. Chem. 100, 2723 (1996)] except now the sticky potential consists of an approximate moment expansion for the dimer interaction potential, which is much faster than the true moment expansion. The object here is to demonstrate that the SSDQO potential energy function can accurately mimic the potential energy function of a multipoint model using the moments of that model. First, the SSDQO potential energy function using the dipole, quadruple, and octupole moments from SPC/E, TIP3P, or TIP5P is shown to reproduce the dimer potential energy functions of the respective multipoint model. In addition, in Monte Carlo simulations of the pure liquid at room temperature, SSDQO reproduces radial distribution functions of the respective model. However, the Monte Carlo simulations using the SSDQO model are about three times faster than those using the three-point models and the long-range interactions decay faster for SSDQO (1/r(3) and faster) than for multipoint models (1/r). Moreover, the contribution of each moment to the energetics and other properties can be determined. Overall, the simplicity, efficiency, and accuracy of the SSDQO potential energy function make it potentially very useful for studies of aqueous solvation by computer simulations.     

Design of surface active soluble peptide molecules at the air/water interface

Surface active molecules collect at interfaces and have the potential to be used for water evaporation reduction. The objective of this work is to design surface active soluble peptides that collect at the air/water interface using molecular simulations. Rotational isomeric state Monte Carlo (RISMC) sampling together with a solvation model that we recently invented, the AAD solvation model [Gu, C.; Lustig, S.; Trout, B. J. Phys. Chem. B 2006, 110 (3), 1476-1484] was applied to calculate the adsorption free energy of the peptide molecule at the air/water interface. The results were validated by both molecular dynamics simulations with an explicit solvent model and surface tension measurements on synthesized peptides. It was demonstrated that this approach is able to give a reasonable prediction of surface activity with an approximately 50% hit rate in terms of designed surface active molecules actually being surface active. The relationship between the chemical composition and the surface morphology is also discussed.     

Excess proton solvation and delocalization in a hydrophilic pocket of the proton conducting polymer membrane nafion

Solvation properties of the hydrated excess proton are studied in a hydrophilic pocket of Nafion 117 through a series of molecular dynamics simulations. The multistate empirical valence bond (MS-EVB) methodology, which enables the delocalization of the excess proton through the Grotthuss hopping mechanism, was employed for one of the excess protons in the simulation cell. Simulations were performed such that "classical" nondissociable hydronium cations and a single excess proton treated with the MS-EVB methodology were at a concentration ratio of 39:1. Two degrees of hydration of the Nafion polymer electrolyte membrane were simulated, each displaying the same marked difference between the solvation structures of the classical versus MS-EVB treated (Grotthuss shuttling) excess proton species. These differences are attributed to the solvent dynamics needed to transfer the cation between the solvent separated and contact pair positions about the sulfonic acid counterion. The results demonstrate that it is generally impossible to describe the low pH conditions in the hydrophilic domains of Nafion without the explicit treatment of Grotthuss delocalization in the underlying molecular dynamics model for the excess protons.