Dependence of ion hydration on the sign of the ion s charge
|
|
- Giles Skinner
- 5 years ago
- Views:
Transcription
1 THE JOURNAL OF CHEMICAL PHYSICS 122, Dependence of ion hydration on the sign of the ion s charge Alan Grossfield a) Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, Missouri Received 23 August 2004; accepted 13 October 2004; published online 21 December 2004 The solvation of simple ions in water is studied using molecular dynamics simulations with a polarizable force field. Previous simulations using this potential demonstrated that anions are more favorably solvated in water than cations. The present work is an attempt to explain this result by examining the effects of ions on the surrounding water structure, with particular focus on the first solvation shell and its interactions with the surrounding water. We conclude that while the first solvation shell surrounding cations is frustrated by competition between ion-water and water-water interactions, solvation of anions is compatible with good water-water interactions American Institute of Physics. DOI: / I. INTRODUCTION a Current address: Biomolecular Dynamics and Scalable Modeling, IBM T. J. Watson Research Center, 1101 Kitchawan Road, P.O. Box 218, Yorktown Heights, NY Telephone: Fax: Electronic mail: agrossf@us.ibm.com The solvation of ions is of great importance in biological and physical chemistry. In particular, understanding ion solvation is critical to elucidating the mechanism by which ion channels distinguish between different ionic species. Viewed simply, alkali and halide ions are perfectly symmetric spheres distinguished only by their charge, size, and polarizability. While these differences can be quite small, they underlie the remarkable selectivity of biological ion channels. For example, although K and Na differ in size by only a few tenths of an angstrom, potassium channels such as KcsA transport K times faster. 1,2 Moreover, the fluctuations in the KcsA selectivity filter are significantly larger than the size difference between the ions, so the differentiation is not accomplished by simple steric incompatibility between the ion and the protein. 3 Rather, it is the thermodynamics of solvation which dictates selectivity. The solvation of the sodium ion in the channel is not as favorable as solvation in bulk water, resulting in a free energy barrier which slows its transfer. Thus, in order to understand the molecular mechanism of ion selectivity, we must first understand ionic hydration. Significant efforts have been made to explore the physical properties of hydrated ions using a variety of experimental techniques, including x-ray and neutron diffraction and scattering, nuclear magnetic resonance, and vibrational spectroscopy. 4,5 However, these experiments are complicated by a number of factors, one of the most important of which is the requirement for electroneutrality. Specifically, electroneutrality means that it is not possible to study a given ionic species in solution at equilibrium without including neutralizing counterions. Thus, while the solvation thermodynamics of a neutral salt can be measured directly, it is impossible experimentally to separate them into contributions from the individual anion and cation. Instead, an additional assumption must be introduced to perform this separation. For example, one could choose a reference salt and make the assumption that the solvation free energies of the cation and anion are equal. 6 8 Once the solvation free energy for a single ionic species is known unambiguously, all of the others can be determined using the principle of additivity and the data for the thermodynamics of neutral salts. Alternatively, one could perform the separation by making assumptions about the solvation of H 3 O and OH to establish the solvation free energy of the proton, 9 or by extrapolating bulk solvation free energies from thermodynamic data on small water-ion cluster. 10 Unfortunately, the single ion solvation values obtained by these separations are very sensitive to the assumptions used to make the separation. As a result, there is considerable controversy as to the true solvation thermodynamics of individual ionic species with different tabulations in the literature disagreeing by as much as 15 kcal/mol. 6,9 16 Molecular dynamics computer calculations can potentially resolve many of these issues. In contrast to the experiments, simulations can treat a single ion in solvation without difficulty, and as such can be used to compute single ion thermodynamics unambiguously. Moreover, molecular dynamics simulations contain a wealth of structural information which can be used to characterize the results. However, in order to perform these simulations, one must first choose repulsion-dispersion parameters for the ions. This has most commonly been done by choosing parameters which reproduce single ion solvation free energies Unfortunately, this process is complicated by the inconsistencies in the literature tabulations of single ion thermodynamic values, as discussed above, 6,9,10 which introduces large uncertainties into the parametrization process. As a result, the various ionic parameter sets have very different thermodynamic and structural properties when used in solution. 11,20 A better approach would be to parametrize using data which does not suffer from the ambiguities introduced by the requirement for electroneutrality. For example, ab initio calculations 21,22 and mass spectroscopy experiments 23 involving small water-ion clusters in the gas phase can provide the kind of detailed structural and thermodynamic information required for parametrization. However, standard force /2005/122(2)/024506/10/$ , American Institute of Physics
2 Alan Grossfield J. Chem. Phys. 122, (2005) fields use fixed partial charges and represent electronic polarization implicitly, which means that parameters chosen for the gas phase are not appropriate for liquid simulations. For example, the dipole moment for water in gas phase is 1.85 D, while in liquid it is greater than 2.6 D; 24 using a rigid nonpolarizable water model intended for the bulk to fit gas phase ion-water interactions would produce ion parameters which overestimate the effective size of the ion in liquid. Moreover, previous work demonstrated that ion-water dimers computed using a standard nonpolarizable water model and ion parameters intended to reproduce bulk solvation to overestimate both the effective size of the ion and the water-ion dimer interaction energy in gas phase. 11 By contrast, fitting ion parameters to the ion-water dimer using the polarizable AMOEBA atomic multipole optimized energizers for biomolecular applications force field produces single ion solvation free energies which, when summed, correctly reproduce the experimental free energies for neutral salts to within 1 kcal/mol in multiple solvents. 11 By contrast, simulations using other force fields do not reproduce the salt values, erring by at least 5 kcal/mol. One of the oldest and simplest models for the free energy of solvation for an ion is the Born equation 25 q2 TABLE I. Parameters for ions. R min and are the parameters for the buffered 14-7 repulsion-dispersion potential, in angstrom and kcal/mol, respectively, and is the polarizability in Å 3. The parameters for Na,K,and Cl are identical to those used in the work by Grossfield et al. Ref. 11. Ion R min Na K Cl Br A 8 0 R 1 1 1, where A is the free energy, q is the charge, R is the effective radius of the ion, and is the dielectric constant of the medium. This simple model makes several specific predictions regarding relative ionic solvation free energies: i smaller ions should have more favorable solvation free energies; ii the free energy to transfer an ion to a more polar phase should always be negative; and iii the sign of the ion s charge does not matter. Although many groups have suggested improvements to the Born equation 16,26 29 these basic predictions are intrinsic to continuum theory. These predictions are not consistent with our previous results. 11 We found that Cl, though clearly larger than K, had an aqueous solvation free energy that was 12 kcal/mol more favorable. Moreover, Cl is more favorably solvated in water 80 than in formamide 110. Several recent theoretical works have produced analogous results for aqueous solvation. For example, the favorable solvation of Cl is also predicted by a recent extension to the Born equation which takes into account the contribution of favorable dispersion interactions. 30 Moreover, Hummer et al. performed a series of Monte Carlo simulations, where they found that the free energy to charge a van der Waals vdw sphere was approximately the size of a methane molecule to e, with the anion being more favorably solvated by more than 40 kcal/mol. 31 The authors hypothesized that the overly favorable anion solvation is in part an artifact of the absence of a vdw term between the water hydrogens and the ions, which allows the hydrogens to pack very close to the anions. However, this is not the case for the results of Grossfield et al., 11 since the AMOEBA water model includes repulsiondispersion terms on water hydrogens. 32 The present work will explain the solvation free energies of ions in water by examining their effects on solvent structure, using molecular dynamics simulations with the polarizable AMOEBA force field. By comparing the solvent structure around a series of cations, anions, and neutral particles, we will distinguish the effect of the ions size and charge, and explain the observed deviations from the Born equation at a molecular level. II. METHODS Molecular dynamics simulations were performed on systems consisting of a single ion and 216 water molecules embedded in a periodic cubic box Å on a side. The AMOEBA force field was used for all calculations. 24,33 This potential uses a sophisticated electrostatics model involving permanent partial charges, dipoles, and quadrupoles on each atom, with electronic polarization represented via induced dipoles. Long range electrostatics are accounted for using Ewald summation, 33,34 with the real space cutoff set to 9 Å, the Ewald coefficient set to 0.42 Å 1, and tin-foil boundary conditions. Repulsion-dispersion or vdw interactions are represented using a buffered 14-7 potential, 35 and were smoothly truncated at 12 Å. Previous work has shown that this force field accurately reproduces the properties of liquid water under both ambient conditions 24 and at extremes of temperature and pressure, 36 and is able to accurately predict the solvation free energies of neutral salts in multiple solvents. 11 The parameters for the water model are those given by Ren and Ponder, 24 while the ion parameters are shown in Table I. The equations of motion were integrated using the Beeman variant of the velocity Verlet integrater, using a1fs time step. 37 Coordinates were saved every 0.1 ps. The temperature was set to 300 K using the Berendsen weak coupling method. 38 Induced dipoles were iterated until no dipole changed by more than 0.01 D upon successive iterations. Each simulation was run for 200 ps, with the first 50 ps discarded as equilibration. All simulations were run using version 3.9 of the TINKER simulation package, with local modifications. 39 Analysis was performed using various programs from the TINKER package or locally written programs and scripts. Simulations were performed examining the solvation of four naturally occurring ionic species: K,Na,Cl, and Br. In addition, a second set of simulations was performed for K, Na, Cl, and Br ; for these simulations, the charge on the ion was inverted, while the vdw and polarizability parameters were left unchanged. Similarly, two more sets of simulations were performed for these species with
3 Dependence of ion hydration on the sign of the ion s charge J. Chem. Phys. 122, (2005) TABLE II. Effective radii for the ions, measured as the distance to the first minimum of the ion-oxygen radial distribution function. Positive and negative refer to the ions with 1e and 1e charges. Neutral refers to the species with charge and polarizability set to zero, although the results for polarizable uncharged species were virtually identical. Ion Positive Negative Neutral Na K Cl Br their charges set to zero: one with the polarizability set the same as for the naturally occurring species and one with the polarizability set to zero. The electrostatic potential at the neutral ions was estimated as V U q U 0, 2 q where U(q) is the system s potential energy with a test charge q placed on the ion. The potentials were computed with q e, although care was taken to verify that the same answer was obtained across a range of q values, which means that O(q 2 ) terms due to the Born self-energy, finite size corrections, and polarization due to the test charge could be safely neglected. 31,40 43 When analyzing solvent structure, a water-water hydrogen bond was defined to be an interaction where the O H distance was less than 2.8 Å, and the O H O angle was greater than 120. This specific definition is reasonable but somewhat arbitrary, in that a different choice would slightly alter the number of hydrogen bonds we would find. This is not particularly important, as we are mostly concerned with variations in the hydrogen bonding patterns around different ionic species. The statistical uncertainty in the number of hydrogen bonds was computed using Monte Carlo bootstrap error analysis. 11,44 For a time series with N data points, a bootstrap trial consisted of selecting N points randomly from the time series, allowing duplication, and computing the average for that sample. We performed 1000 bootstrap trials and computed the statistical uncertainty as the standard deviation of the averages from the trials. The lifetime of water structure in the first solvation shell was computed from the autocorrelation function of the instantaneous ionic coordination number. The coordination number was computed as the number of waters with oxygens within the effective ionic radius, which we defined to be the distance to the first minimum of the ion-o radial distribution function see Tables II and III. We then computed the autocorrelation function of this time series as C N c N c N c 0 N c, 3 N c where N c (t) is the coordination number for the ion at time t, (N c ) is the variance, and angle brackets indicate the thermodynamic average. For all of the ions studied, the correlation function was not well described as a single exponential, but did fit well to a sum of two exponentials TABLE III. Coordination number for the ions, measured as the number of water oxygens within the effective radius, as defined in Table II. Positive and negative refer to the ions with 1e and 1e charges. The values marked with asterisks differ slightly from those reported in the work by Grossfield et al. Ref. 11 because they were computed by direct analysis of the trajectories instead of summing the RDF. The difference between the reported Cl values is due to an error in the previous work. Ion Positive Negative Na 5.9* K 6.9* Cl * 1.3 Br C t w 1 e t/ 1 w 2 e t/ 2. 4 At long time intervals, the correlation function became noisy due to poor sampling. To prevent this from influencing the fit, only the first 10 ps of the autocorrelation function was used. Because we only recorded coordinate sets every 0.1 ps, we cannot detect events occurring on a faster time scale. As a result, if a water molecule were to leave the first solvation shell and be replaced within 0.1 ps, the event would not be accounted for. III. RESULTS A. Solvent structure around ions The primary theoretical tool for examining molecular structure in solution is the radial distribution function RDF. Figure 1 shows the radial distribution functions for water FIG. 1. Radial distribution functions for water oxygens around the ions. a The standard RDFs and b the RDFs shifted such that the first peaks overlap.
4 Alan Grossfield J. Chem. Phys. 122, (2005) TABLE IV. The heights of the first peak and minimum for the ion-o radial distribution functions. Positive Negative Neutral Ion Peak Well Peak Well Peak Well Na K Cl Br FIG. 2. Radial distribution functions for water oxygens around the potassium. K,K, and K 0 refer to the cationic, anionic, and neutral forms for potassium. The ions had the same vdw parameters in all three calculations. K was assigned the same polarizability as K, while for K 0 the polarizability was 0 Å 3. oxygens around the naturally occurring ionic species. The results in Fig. 1 a are completely as expected: as the ions get larger, the first peak moves outward and becomes less pronounced. In Fig. 1 b, the direct effects of ion size are removed from the plots by shifting the RDFs such that their first peaks are in the same position. We can see that the location of the first minimum and second peak relative to the first peak are largely unaffected by the ion size with the exception of Na, where both appear to be shifted slightly outward. Moreover, the larger the ion, the less structured the solvent appears to be, as evidenced by the decreased depths of the first minimum and height of the second peak. However, these data alone do not allow us to directly separate the effects of ionic size and charge on solvent structure because both variables are changed simultaneously. We approach this problem by examining simulations of several unphysical ions, such as Na and Cl, which we model by inverting the ionic charge while leaving the vdw and polarizability parameters unchanged. Although these ions do not represent realistic physical entities, the resulting calculations cleanly distinguish the effects of ionic size and charge. With this in mind, Fig. 2 shows the distribution of water oxygens around potassium as a function of the ion s charge. The first peaks for K and K have virtually identical heights, with the latter shifted outward by roughly 0.2 Å. By contrast, the first minimum is in virtually the same place in both curves. The oxygen peak shifts because the water oxygen directly solvates the cation, while the anion is solvated by the water hydrogen. This may also account for the fact that the RDF about the anion appears slightly more structured, with a deeper first minimum and higher second peak. The RDF for water around the nonpolarizable and uncharged K 0 is shown as well. It is qualitatively different from that of the two charged species; the first peak is much lower and broader and is located much farther out. This is to be expected, since eliminating the charge allows unfavorable vdw interactions between the ion and the first solvation shell to relax. Moreover, the vdw interactions the neutral particles can make are not sufficiently favorable to compete with water-water hydrogen bonding interactions. The RDFs produced when K 0 was assigned the same polarizability as K were indistinguishable from the nonpolarizable case data not shown. Although Fig. 2 only shows the RDFs for potassium, the other elements behave similarly. These results make it clear that the ion s charge can in principle have a significant effect on the its apparent size in solution. In estimating the ions size, we choose to define the effective radius as the location of the first minimum of the ion-o RDF. Physically, this corresponds to including the first solvation shell around the ion in our definition of its size. Table II shows the effective radii for all of the ions simulated. For each of the four elements, the charged species have similar effective radii regardless of the charge sign, while the radii of the uncharged species are roughly 1.5 Å larger. This does not imply, however, that the solvent packing around the ions is independent of charge. Rather, the coordination number for the ions, computed as the average number of water oxygens located within the first solvation shell of the ion, varies dramatically with the ions sign. Table III summarizes these results: for ions with identical vdw potentials, the cation is coordinated by roughly one more water molecule than the equivalent anion. The coordination numbers computed here are somewhat larger than those computed via quantum molecular dynamics and quasichemical theory. 45,46 However, the values produced by the two ab initio methods differ somewhat, and it is not entirely clear that the optimum coordination number as defined by quasichemical theory should be equal to that produced by analyzing the radial distribution function. The ion coordination number typically shows fluctuations of more than one water molecule, indicating that ion-coordination involves a number of different first shell conformations. Table IV characterizes the ion-o RDFs for the various species in terms of the heights of the first peak and minimum. As before, when the ions get bigger, the degree of structure in the waters around them diminishes, demonstrated by the lower peak heights and higher minima. Again, the RDFs for each element do not appear to be particularly sensitive to the sign of the charge on the ion. This trend is not obeyed by the solvation of the neutral species, where the peak heights and well depths are constant. Indeed, if the RDFs for the neutral particles are shifted such that the loca-
5 Dependence of ion hydration on the sign of the ion s charge J. Chem. Phys. 122, (2005) FIG. 3. Electrostatic potential at the center of the neutral ions, in kcal/ mol-e, as a function of the effective ionic radii from Table II. Curves are shown for neutral species with and without polarizability. The error bars represent standard deviations. tion of the first peaks are equivalent, they superpose almost perfectly, which would appear to indicate that the basic solvation mechanism is maintained across a range of particle sizes. Another measure of the solvent ordering around a solute is the electrostatic potential. By definition, the electrostatic potential is on average constant everywhere in an isotropic solution. However, the presence of a solute breaks the symmetry, and if the waters surrounding a solute are ordered, a net electrostatic potential can be induced, even when the solute makes no electrostatic interactions of its own. When we computed the electrostatic potential at the center of the neutral particles using Eq. 2, we found that the potential was positive for all species. Figure 3 shows the electrostatic potential plotted as a function of the effective radius of the neutral species. These results are qualitatively consistent with those obtained by Ashbaugh. 47 It has been argued that the electrostatic potential at the neutral particle ought to make a roughly linear contribution to the charging free energy for the ion. 31,47 This would amount to a difference of 10 kcal/mol between the solvation free energies of otherwise equivalent anions and cations, which is similar to the difference in solvation free energies for K and Cl computed previously. 11 Figure 3 also compares the electrostatic potentials computed for ions with and without explicit polarizability. The curves differ significantly, in that the potential is initially smaller for the polarizable species but drops more slowly at the particles get bigger. This is interesting, because the water oxygen RDFs appear to be independent of the presence of polarizability on the ion. The effect is likely due to subtle adjustments in the hydrogen packing, which are not readily visible in the ion-oxygen RDFs. B. Solvent structure around the first solvation shell FIG. 4. Oxygen-oxygen RDFs for waters in the first solvation shell. RDFs were computed for water oxygens, including only those pairs where at one of the waters was in the first solvation shell of the ion. a The distribution functions about cations and b for the anions. In both plots, bulk refers to O O RDF from an equivalent neat water simulation. Because of the strong interactions between the ion and the immediately surrounding water, one can effectively consider the first solvation shell to be a part of the ion and examine the solvation of the whole ion-water complex. This is analogous to the logic behind quasichemical calculations of ion solvation. 45,46,48 51 In this case, we ask the question, how do the water-water interactions made by the waters immediately surrounding an ion depend on the ion s charge and size? We define the first solvation shell geometrically as all the waters with oxygens within the effective radius of the ion, as listed in Table II. It should be noted that other definitions are also reasonable; for example, when using quasichemical theory the coordination number and first solvation shell are defined thermodynamically. 15 In order to directly examine the solvent structure about the waters in the first solvation shell, we computed oxygenoxygen RDFs which included only pairs of oxygens where at least one water molecule was in the first shell. The results are shown in Fig. 4. Figure 4 a contains the RDFs for the first shells around cations, while Fig. 4 b shows the same for the anions. For comparison purposes, the equivalent O O RDF from a neat water solution is included in both plots. In analyzing Fig. 4, we begin with the assumption that the bulk water RDF represents the minimum free energy distribution. It then follows that any uncompensated perturbation to that distribution is thermodynamically unfavorable. In Fig. 4 a, all of the cations diminish the height of the first peak to roughly the same degree, with the larger ions having a slightly lesser effect. Moreover, the smaller ions, especially
6 Alan Grossfield J. Chem. Phys. 122, (2005) TABLE V. Average number of water-water hydrogen bonds formed by waters in the first solvation shell around the ions. Intra describes hydrogen bonds between two waters in the first solvation shell, while inter counts hydrogen bonds to waters outside the first shell. The statistical uncertainty in these values is between and Positive Negative Ion Total Intra Inter Total Intra Inter Na K Cl Br FIG. 5. Probability distribution for the cosine of the angle between the permanent dipole moment of waters in the first solvation shell and the ionoxygen vector. a The distribution functions for the cations and b for the anions. Na, eliminate the first minimum. The reduction in the first peak is due to volume excluded by the ion, while the diminished minimum indicates that while the first shell waters do have other waters packed around them, they are not part of the typical water-hydrogen bonding network. The differences from the bulk water RDF decrease with increasing effective ionic radius, implying that as the ion-water interactions grow weaker, the water-water interactions become more dominant. This tradeoff between optimizing ion-water and water-water interactions indicates a degree of frustration in the solvation of small cations. By contrast, the RDFs for the anions, shown in Fig. 4 b, are very similar to the bulk water curve, aside from the fact that the height of the first peak is diminished. With the exception of the smallest anion Na, the first minimum and long range solvent structure are virtually identical to that for bulk water. This indicates that the presence of the anions does not disrupt the waters hydrogen bonding network the way cations do. This dependence on the charge sign is not reflected in the number of hydrogen bonds made by waters in the first solvation shell. As shown in Table V, the total number of hydrogen bonds per water molecule is almost independent of the ionic charge and size. There is a slight shift from hydrogen bonds to bulk waters to hydrogen bonding within the first solvation shell as the ions get larger. Moreover, the increase in intrashell hydrogen bonding is more pronounced for the anions. These trends, while small, do appear to be statistically significant; using Monte Carlo bootstrap error analysis, we estimate that the statistical uncertainties in Table V range from to In order to characterize the structure of the first solvation shell in greater detail, we examined the angles between the waters and the ion. Specifically, we computed the probability distribution for the cosine of the angle between the oxygenion vector and the permanent dipole moment of the water. The permanent dipole moment of the water is essentially parallel to the bisector of the H O H angle; there are small deviations due to fluctuations in the O H bond lengths. We use the cosine rather than the angle itself because the cosine would have a uniform distribution in the absence of any interactions. The results are shown in Fig. 5. Figure 5 a shows the distributions for the cations. The most probable water orientation around the smaller cations is around cos 1, corresponding to a structure with the oxygen pointing directly at the ion. However, the distributions are broad, and even for Na all orientations are seen. As the ions get bigger, the distribution becomes broader, and the maximum moves away from the parallel orientation; Cl and Br have maxima around cos 0.75 or 40. It is interesting to note that the probability distributions for all four cations intersect at cos 0.45 or about 65. One could speculate that the angle distributions are composed of two distinct basins of attraction, where one is dominated by ion-water interactions and the other by water-water interactions. If the relative probability of orienting a water molecule with an angle of 65 was the same in both basins, one would expect to see the probability of this angle remain constant even as the ions became larger and the relative importance of ion-water and water-water interactions shifted. The distribution functions for the anions, shown in Fig. 5 b, are quite different. For all four elements, the distributions are strongly peaked around cos 0.6 corresponding to 125. The complement of this angle, 55, is almost exactly half of the H O H angle, indicating that water primarily points one hydrogen toward the anion. The peaks are broader for the larger anions, but their location does not shift. The qualitative difference between water structure around cations and anions is indicative of the competition between making strong electrostatic interactions with the ion and maintaining favorable water-water interactions. In the case of the cation, these goals appear to be at odds, with the result that the as the ions grow larger, the distribution shifts
7 Dependence of ion hydration on the sign of the ion s charge J. Chem. Phys. 122, (2005) TABLE VI. Lifetimes for ion coordination. The autocorrelation function for the instantaneous coordination number for the ions was fit to a sum of two exponentials see Eq. 4. The decay times 1 and 2 are in picoseconds. Positive Negative Ion w 1 1 w 2 2 w 1 1 w 2 2 Na K Cl Br away from the conformation which optimizes ion-water electrostatics. By contrast, the anion distribution remains centered around the same basic orientation for all of the elements. The weaker electrostatic interactions made by larger ions result in broadening the distribution, but there is no evidence of a competing alternative conformation. Given the strong electrostatic field surrounding an ion, it is interesting to note that angle probability distribution functions computed using the total dipole moment vector, where the effects of polarization are included, are very similar to those using the permanent dipole data not shown. The only exception is Na, where the peak shifts to cos 0.72 or 44. However, this does not mean that the ionic electric field completely dominates the water s interactions. For the rest of the ions considered, the angle between the induced dipole on the water and the ion-o vector follows roughly the same distribution as the permanent dipole, and even in the case of Na the induced dipole does not typically point directly at the ion. Furthermore, these electric fields are not as a rule larger in magnitude than those seen in bulk water if one uses the magnitude of the total dipole moment per molecule as the metric. With the AMOEBA water model, bulk water has an average dipole moment of 2.77 D under the conditions simulated here. Waters in the first solvation shell around all of the cations had a slightly smaller dipole moment around 2.72 D. Similar behavior was seen for waters around Cl and Br, where the average dipole moment was 2.74 D. Only the smaller anions caused the overall dipole moment to increase, K with 2.81 D and Na with 3.11 D. We also examined the kinetics of waters in the first solvation shell. Specifically, we examined the fluctuations in the coordination numbers for different ionic species. As described in Sec. II, we computed the autocorrelation function for the coordination number and fit it to a sum of two exponentials. The results are shown in Table VI. In general, there is a fast component to the decay with a time constant 2 of 0.1 ps, which accounts for roughly half of the decay, and a slower phase with a time constant 1 around 1 ps, which accounts for the rest. However, there is significant variation which is not easily explained by simple trends in size or charge. For Na and K, the behavior is mostly independent of the ionic charge, while for the larger ions the cations and anions behave differently. Cl and Br have longer 1 values than the equivalent cations roughly 2 ps instead of 1 ps. However, the slow component accounts for less of the overall decay, where the other ions weight the two components roughly equally. Physically, we interpret 1 to be the time scale for rearrangement of the first solvation shell, while 2 most likely represents an artifactual flickering as water molecules repeatedly move outside the ion s effective radius without actually leaving the general vicinity of the ion. Thus, we interpret the drop in 1 from Na to K as indicating that the larger K species bind their first solvation shell less tightly. By contrast, the longer 1 for the larger anions most likely has a different physical origin. We believe it is further indication of the lack of frustration in the first solvation shell s interaction with the bulk. IV. DISCUSSION A. Charge sign dependence of ion solvation In a prior paper, 11 we computed the solvation free energies of several ions and observed an interesting phenomenon: the solvation free energy for Cl in water was more favorable than that for K by roughly 12 kcal/mol and only 5 kcal/mol less favorable than that for Na. This is in contradiction to what one would expect from a naive application of continuum theory, which would predict that the solvation free energy should increase monotonically with ionic radius, regardless of the sign of the ion s charge. Since the expected behavior was seen in analogous simulations in formamide, we believe this result points to something interesting about the properties of water as a solvent. The primary goal of the present work has been to explain these results by analyzing the effects of ions on water structure as a function of the ions size and charge. It has long been common to classify ions into kosmotropes structure forming and chaotropes structure breakers, according to their presumed effects on water structure. 52 However, in the present simulations we see little or no evidence for a global effect of the ions on solvent structure. Rather, the structural properties of waters outside the first solvation shell around the ion are indistinguishable from those seen in simulations of bulk water data not shown. This is in accordance with recent vibrational spectroscopy experiments, which also reject the notion of ionic effects on global water structure. 53,54 There are a variety of hypotheses one could attempt to use to explain the surprisingly favorable solvation free energy for Cl. For example, one could hypothesize that its
8 Alan Grossfield J. Chem. Phys. 122, (2005) effective radius is somehow much smaller that what one would expect based on its crystal radius. However, as shown in Table II, the spacing of waters around chloride follows the expected chemical trend. Moreover, there is no evidence that the ion s charge significantly alters the effective radius. The same is true when the peak heights and well depths are examined in Table IV. Given the importance of hydrogen bonding in driving water structure, it is reasonable to look at the effect of ions on the number of hydrogen bonds formed by waters in their immediate vicinity. A recent set of calculations using a stylized water model simulated in two dimensions concluded that anions greatly reduce the number of hydrogen bonds made within the first solvation shell, and that the number of intrashell hydrogen bonds increases with ionic radius, regardless of the charge sign. 55 As seen in Table V, the present simulations confirm the latter result; going from the smallest Na to largest Br ions increases the number of intrashell hydrogen bonds per water molecule by roughly 0.4. However, anions and cations show very similar behavior, which leads us to believe that this is a size effect driven by packing rather than an effect of the ion s charge. Moreover, the total number of hydrogen bonds made by first shell waters remains virtually constant across all of the ions simulated, which makes the observed variation an unlikely candidate to explain the large differences in solvation free energy. Based on the evidence presented above, we believe that the physical origin of favorable anion solvation is the elimination of frustration in the interactions between the first solvation shell and the bulk water. We base this belief primarily on the evidence from the RDFs computed for these waters see Fig. 4 ; while the water-water structure is disrupted around cations, the water-water RDF for the first shell around anions looks almost identical to that for bulk water. The behavior of water around an uncharged sphere is consistent with this explanation. Because a van der Waals sphere cannot make interactions strong enough to favorably compete with water-water hydrogen bonds, it acts as a probe of what the water wants to do, as opposed to what it is forced to do by electrostatic interactions with the ion. As shown in Fig. 3, the water induces a positive electrostatic potential at the center of a neutral sphere. While it is clear that adding a charge changes the water structure see Fig. 2, this does seem to indicate a greater overall compatibility of water structure for anions over cations. This behavior can also be seen in Fig. 5. Where all of the anions show a single peak in the angle between the water permanent dipole and water around 125, there is a clear trend in the cation simulations. Waters around Na are most likely to orient their dipoles directly at the ion, but this tendency diminishes as the cations get bigger. This is a direct evidence of competition between water-water and watercation interactions, which appears to be absent with anions. It also points out the naivety of assuming that the water-ion interaction is simply due to charge-dipole interactions. Indeed, a model which does so will mispredict the angle distribution, as well as the resulting water-water packing around the ion. 55 By contrast, angle distributions from all-atom simulations appear very similar to those presented here. 56,57 B. Comparison of computational methods In recent years, a large number of groups have performed simulations exploring the solvation of ions in water, using a variety of techniques ranging from classical force field calculations involving nonpolarizable 20,57 60 and polarizable 11,56,61 71 force fields, to mixed quantum-classical simulations QM/MM, 15,48,72 74 fully quantum mechanical calculations, 45,46,72,75 and highly simplified models. 30,55,76 Each approach has its advantages and disadvantages. Very simple models tend to make for very inexpensive calculations, so that obtaining sufficient sampling or system size need not be an issue. Moreover, a model with very few terms can make it easy to directly attribute which properties of the simulation are due to which components of the model. However, almost by definition these models involve drastic approximations, and it is not always clear which results are due to the approximations as opposed to the features of the system modeled. At the other end of the spectrum are fully quantum mechanical calculations. In principle, these methods represent the system with great fidelity. However, they are very expensive computationally, with the result that the systems must be quite small and the sampling very limited; for example, Heuft and Meijer recently used quantum molecular dynamics to simulate a system of 64 waters and one anion for a total of 17 ps. 75 Similarly, Rempe and Pratt simulated Na solvated by 32 water molecules for 12 ps. 45 Such simulations typically use small basis sets and either density functional theory or Hartree-Fock theory, both of which neglect attractive dispersion interactions. Still, it is clear that some phenomena, such as solvation of the proton in water, are primarily quantum mechanical in nature and cannot be captured easily without using ab initio methods. Purely classical calculations and QM/MM simulations fall somewhere in between these two extremes, both in terms of expense and accuracy. If one is to use a classical force field, several groups have demonstrated the importance of including explicit electronic polarization in order to accurately describe ionic solvation, especially anions and ions near interfaces, to obtain quantitative accuracy. 11,56,63,67,77,78 However, other groups have shown that it is possible to obtain qualitatively interesting results while only including polarization implicitly. 58 QM/MM simulations clearly capture electronic polarization at least in the quantum mechanical part of the simulation, and in principle everywhere if a polarizable force field is used for the molecular mechanics part, and while they are significantly more expensive than the equivalent classical calculation, they are efficient enough to allow reasonable system sizes and sampling times. However, like the purely quantum mechanical calculations, they typically use theories which do not contain dispersion interactions. The interface with classical models requires that these interactions be included using empirical potentials. 74,79,80 In principle, this could introduce significant uncertainty into the calculation; in practice, it is not yet clear how sensitive these simulations are to a specific choice of dispersion parameters. 80,81 Quasichemical theory provides another framework for combining ab initio and classical methods. 49 In this ap-
9 Dependence of ion hydration on the sign of the ion s charge J. Chem. Phys. 122, (2005) proach, a superion consisting of the ion and a number of inner-shell water molecules are modeled using quantum mechanical methods in gas phase, then transferred into bulk solution using either continuum electrostatics or molecular modeling techniques. Because only a single ab initio optimization is performed for each superion cluster although several cluster sizes must be considered to determine the optimal coordination number, the remainder of the calculation is performed classically, which reduces their computational cost. While this theory correctly accounts for the statistical mechanics of ion transfer, its usual implementation neglects packing and dispersion interactions between the superion and the bulk, 46,51,82 which diminishes the quantitative accuracy of the method. The method also assumes that the fluctuations of the first solvation shell are accurately represented by harmonic fluctuations about the gas phase minimum energy structure for each inner shell coordination number. This approach is analogous to the application of the inherent structures formalism to the interpretation of the behavior of bulk liquids. 83 Moreover, as with other QM/MM methods it is not always clear how to model the interactions between the classical and quantum mechanical portions of the system. Still, this technique has been used to produce a number of interesting results, including calculations of ion solvation free energies and coordination numbers, 45,46,48,50,51,84 and it does provide a clear framework from which one can examine ion solvation. V. CONCLUSION We have presented molecular dynamics simulations using a polarizable force field which explains the differences in solvation between cations and anions. Specifically, the waters in the first solvation shell around an anion are able to maintain water-water interactions which greatly resemble those formed in the bulk liquid, while the waters around smaller cations cannot. However, the present simulations have only treated monovalent atomic ions, and it would be interesting to see if this behavior is retained for more complicated ions such as Ca 2,SO 4 2, or ClO 4. Similarly, the effects of ions on water structure near macromolecular assemblies, such as nucleic acids, protein, and lipid bilayers, are of both biological and physical interest. ACKNOWLEDGMENTS The author would like to thank Jay Ponder, Pengyu Ren, Michael Schnieders, Jonathan Sachs, Horia Petrache, Henry Ashbaugh, Dilip Asthagiri, and Adrian Parsegian for many helpful discussions. This work was funded in part by NIH NRSA fellowship Grant No. 5F32NS D. A. Doyle, J. M. Cabrat, R. A. Pfuetzner, A. Kuo, J. M. Gulbis, S. L. Cohen, B. T. Chait, and R. MacKinnon, Science 280, B. Hille, Ionic Channels of Excitable Membranes Sinauer, Sunderland, MA, Y. Zhou and R. MacKinnon, Biochemistry 43, H. Ohtaki and T. Radnai, Chem. Rev. Washington, D.C. 93, H. Ohtaki, Monatsch. Chem. 132, H. L. Friedman and C. V. Krishnan, in Water: A Comprehensive Treatise, edited by F. Franks Plenum, New York, 1973, Vol. 3, pp R. A. Marcus and N. Sutin, Biochim. Biophys. Acta 811, C. Kalidas, G. Hefter, and Y. Marcus, Chem. Rev. Washington, D.C. 100, R. Schmid, A. M. Miah, and V. N. Sapunov, Phys. Chem. Chem. Phys. 2, M. D. Tissandier, K. A. Cowen, W. Y. Feng, E. Gundlach, M. H. Cohen, A. D. Earhart, J. V. Coe, and T. R. J. Tuttle, J. Phys. Chem. A 102, A. Grossfield, P. Ren, and J. W. Ponder, J. Am. Chem. Soc. 125, C.-G. Zhang and D. A. Dixon, J. Phys. Chem. A 105, C.-G. Zhang and D. A. Dixon, J. Phys. Chem. A 106, C.-G. Zhang and D. A. Dixon, J. Phys. Chem. A 108, D. Asthagiri, L. R. Pratt, and H. S. Ashbaugh, J. Chem. Phys. 119, W. M. Latimer, K. S. Pitzer, and C. M. Slansky, J. Chem. Phys. 7, J. Åqvist, J. Phys. Chem. 94, A. D. MacKerell, Jr., D. Bashford, M. Bellott et al., J. Phys. Chem. B 102, B. Roux and S. Berneche, Biophys. J. 82, M. Patra and M. Karttunen, J. Comput. Chem. 25, D. Feller, E. D. Glendening, D. E. Woon, and M. W. Feyereisen, J. Chem. Phys. 103, S. S. Xantheas, J. Phys. Chem. 100, K. Hiraoka, S. Mizuse, and S. Yamabe, J. Phys. Chem. 92, P. Ren and J. W. Ponder, J. Phys. Chem. B 107, M. Born, Z. Phys. 1, A. A. Rashin and B. Honig, J. Phys. Chem. 89, T. Abe, J. Phys. Chem. 90, S. L. Chan and C. Lim, J. Phys. Chem. 98, J. K. Hyun, C. S. Babu, and T. Ichiye, J. Phys. Chem. 99, M. Boström and B. W. Ninham, J. Phys. Chem. B 108, G. Hummer, L. R. Pratt, and A. E. Garcia, J. Phys. Chem. 100, P. Ren and J. W. Ponder, J. Phys. Chem. B to be published. 33 P. Ren and J. W. Ponder, J. Comput. Chem. 23, J. W. Perram, H. G. Petersen, and S. W. De Leeuw, Mol. Phys. 65, T. A. Halgren, J. Am. Chem. Soc. 114, P. Ren and J. W. Ponder, J. Phys. Chem. B 108, D. Beeman, J. Comput. Phys. 20, H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. DiNola, and J. R. Haak, J. Chem. Phys. 81, J. W. Ponder, TINKER: Software Tools for Molecular Design 3.9, dasher.wustl.edu/tinker/ G. Hummer, L. R. Pratt, and A. E. Garcia, J. Chem. Phys. 107, F. Figueirido, G. S. Del Buono, and R. M. Levy, J. Chem. Phys. 103, S. Bogusz, T. E. I. Cheatham, and B. R. Brooks, J. Chem. Phys. 108, P. H. Hünenberger and J. A. McCammon, Biophys. Chem. 78, W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C, 2nd ed. Cambridge University Press, Cambridge, S. B. Rempe and L. R. Pratt, Fluid Phase Equilib. 183, S. B. Rempe, D. Asthagiri, and L. R. Pratt, Phys. Chem. Chem. Phys. 6, H. S. Ashbaugh, J. Phys. Chem. B 104, P. Grabowski, D. Riccardi, M. A. Gomez, D. Asthagiri, and L. R. Pratt, J. Phys. Chem. A 106, M. E. Paulaitis and L. R. Pratt, Adv. Protein Chem. 62, D. Asthagiri, L. R. Pratt, J. D. Kress, and M. A. Gomez, Chem. Phys. Lett. 380, D. Asthagiri, L. R. Pratt, M. E. Paulaitis, and S. B. Rempe, J. Am. Chem. Soc. 126, K. D. Collins and M. W. Washabaugh, Q. Rev. Biophys. 18, A. W. Omta, M. F. Kropman, S. Woutersen, and H. J. Bakker, Science 301, M. F. Kropman and H. J. Bakker, J. Am. Chem. Soc. 126, B. Hribar, N. T. Southall, V. Vlachy, and K. A. Dill, J. Am. Chem. Soc. 124,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationStatistical Theory and Learning from Molecular Simulations
Statistical Theory and Learning from Molecular Simulations Lawrence R. Pratt 1 and Susan B. Rempe 2 1 Department of Chemical & Biomolecular Engineering 2 Center for Biological and Material Sciences, Sandia
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationBiophysics II. Hydrophobic Bio-molecules. Key points to be covered. Molecular Interactions in Bio-molecular Structures - van der Waals Interaction
Biophysics II Key points to be covered By A/Prof. Xiang Yang Liu Biophysics & Micro/nanostructures Lab Department of Physics, NUS 1. van der Waals Interaction 2. Hydrogen bond 3. Hydrophilic vs hydrophobic
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationTitle Super- and subcritical hydration of Thermodynamics of hydration Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2000), 8109 Issue Date 2000-05-08 URL http://hdl.handle.net/2433/50350
More informationAqueous solutions. Solubility of different compounds in water
Aqueous solutions Solubility of different compounds in water The dissolution of molecules into water (in any solvent actually) causes a volume change of the solution; the size of this volume change is
More information6 Hydrophobic interactions
The Physics and Chemistry of Water 6 Hydrophobic interactions A non-polar molecule in water disrupts the H- bond structure by forcing some water molecules to give up their hydrogen bonds. As a result,
More informationMichael W. Mahoney Department of Physics, Yale University, New Haven, Connecticut 06520
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 23 15 DECEMBER 2001 Quantum, intramolecular flexibility, and polarizability effects on the reproduction of the density anomaly of liquid water by simple potential
More informationTheory and Modeling of Specific Ion Hydration
Theory and Modeling of Specific Ion Hydration Thomas L. Beck, University of Cincinnati November 21, 2013 Thanks: Lawrence Pratt, Hank Ashbaugh, Dilip Asthagiri, PNNL group, Kevin Leung, John Weeks, Barry
More informationMolecular Dynamics Simulation Study of the Ionic Mobility of OH Using the OSS2 Model
1154 Bull. Korean Chem. Soc. 2006, Vol. 27, No. 8 Song Hi Lee Molecular Dynamics Simulation Study of the Ionic Mobility of OH Using the OSS2 Model Song Hi Lee Department of Chemistry, Kyungsung University,
More informationLecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions Van der Waals Interactions
More informationIntroduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th
Abstract When conductive, gramicidin monomers are linked by six hydrogen bonds. To understand the details of dissociation and how the channel transits from a state with 6H bonds to ones with 4H bonds or
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationAtomic and molecular interaction forces in biology
Atomic and molecular interaction forces in biology 1 Outline Types of interactions relevant to biology Van der Waals interactions H-bond interactions Some properties of water Hydrophobic effect 2 Types
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationarxiv:chem-ph/ v1 24 May 1995
On the free energy of ionic hydration Gerhard Hummer, ab Lawrence R. Pratt, c and Angel E. García a (April 6, 1995) The Journal of Physical Chemistry (in press, 1995) arxiv:chem-ph/9505005v1 24 May 1995
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More informationLecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions: Van der Waals Interactions
More informationWater structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
Bull. Mater. Sci., Vol. 31, No. 3, June 2008, pp. 525 532. Indian Academy of Sciences. Water structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
More informationINTERMOLECULAR AND SURFACE FORCES
INTERMOLECULAR AND SURFACE FORCES SECOND EDITION JACOB N. ISRAELACHVILI Department of Chemical & Nuclear Engineering and Materials Department University of California, Santa Barbara California, USA ACADEMIC
More informationMolecular dynamics simulation of limiting conductances for LiCl, NaBr, and CsBr in supercritical water
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 2 8 JANUARY 2000 Molecular dynamics simulation of limiting conductances for LiCl, NaBr, and CsBr in supercritical water S. H. Lee Department of Chemistry,
More informationSolvation and reorganization energies in polarizable molecular and continuum solvents
Solvation and reorganization energies in polarizable molecular and continuum solvents Joel S. Bader CuraGen Corporation, 322 East Main Street, Branford, Connecticut 06405 Christian M. Cortis Department
More informationThe broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals.
Physical Metallurgy The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals. Crystal Binding In our discussions
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationSome properties of water
Some properties of water Hydrogen bond network Solvation under the microscope 1 Water solutions Oil and water does not mix at equilibrium essentially due to entropy Substances that does not mix with water
More informationElectrostatics of membrane adhesion
Electrostatics of membrane adhesion S. Marcelja Department of Applied Mathematics, Research School of Physical Sciences and Engineering, The Australian National University, Canberra ACT 6, Australia ABSTRACT
More informationWater models in classical simulations
Water models in classical simulations Maria Fyta Institut für Computerphysik, Universität Stuttgart Stuttgart, Germany Water transparent, odorless, tasteless and ubiquitous really simple: two H atoms attached
More informationNew Perspective on structure and bonding in water using XAS and XRS
New Perspective on structure and bonding in water using XAS and XRS Anders Nilsson Stanford Synchrotron Radiation Laboratory (SSRL) and Stockholm University, Sweden R. Ludwig Angew. Chem. 40, 1808 (2001)
More informationTheoretical models for the solvent effect
Theoretical models for the solvent effect Benedetta Mennucci Dipartimento di Chimica e Chimica Industriale Web: http://benedetta.dcci.unipi.it Email: bene@dcci.unipi.it 8. Excited electronic states in
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
More informationTheoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes
J. At. Mol. Sci. doi: 10.4208/jams.121011.011412a Vol. 3, No. 4, pp. 367-374 November 2012 Theoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes Xiu-Ying Liu a,, Li-Ying Zhang
More informationarxiv: v1 [physics.chem-ph] 30 Jul 2012
Is ion channel selectivity mediated by confined water? Diego Prada-Gracia and Francesco Rao School of Soft Matter, Freiburg Institute for Advanced Studies, Freiburg, Germany arxiv:27.6953v [physics.chem-ph]
More informationTwo-dimensional hydration shells of alkali metal ions at a hydrophobic surface
JOURNAL OF CHEMICAL PHYSICS VOLUME 121, NUMBER 24 22 DECEMBER 2004 Two-dimensional hydration shells of alkali metal ions at a hydrophobic surface Sheng Meng, D. V. Chakarov, B. Kasemo, and Shiwu Gao a)
More informationDO PHYSICS ONLINE STRUCTURE OF THE ATOM FROM IDEAS TO IMPLEMENTATION ATOMS TO TRANSISTORS STRUCTURE OF ATOMS AND SOLIDS
DO PHYSIS ONLINE FROM IDEAS TO IMPLEMENTATION 9.4.3 ATOMS TO TRANSISTORS STRUTURE OF ATOMS AND SOLIDS STRUTURE OF THE ATOM In was not until the early 1930 s that scientists had fully developed a model
More informationWhy Is Molecular Interaction Important in Our Life
Why Is Molecular Interaction Important in ur Life QuLiS and Graduate School of Science iroshima University http://www.nabit.hiroshima-u.ac.jp/iwatasue/indexe.htm Suehiro Iwata Sept. 29, 2007 Department
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More informationEnduring Understandings & Essential Knowledge for AP Chemistry
Enduring Understandings & Essential Knowledge for AP Chemistry Big Idea 1: The chemical elements are fundamental building materials of matter, and all matter can be understood in terms of arrangements
More informationMONTE CARLO SIMULATION OF I -, Br -, AND Cl - IN WATER USING AB INITIO PAIR POTENSIAL FUNCTIONS
4 Indo. J. Chem., 7, 7 (), 4-9 MONTE CARLO SIMULATION OF I -, Br -, AND Cl - IN WATER USING AB INITIO PAIR POTENSIAL FUNCTIONS Austrian-Indonesian Centre for Computational Chemistry (AIC), Department of
More information2 Structure. 2.1 Coulomb interactions
2 Structure 2.1 Coulomb interactions While the information needed for reproduction of living systems is chiefly maintained in the sequence of macromolecules, any practical use of this information must
More informationLecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions. Let's get together.
Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions Let's get together. Most gases are NOT ideal except at very low pressures: Z=1 for ideal gases Intermolecular interactions come
More informationChapter 11. Intermolecular Forces and Liquids & Solids
Chapter 11 Intermolecular Forces and Liquids & Solids The Kinetic Molecular Theory of Liquids & Solids Gases vs. Liquids & Solids difference is distance between molecules Liquids Molecules close together;
More informationELECTRONIC STRUCTURE OF MAGNESIUM OXIDE
Int. J. Chem. Sci.: 8(3), 2010, 1749-1756 ELECTRONIC STRUCTURE OF MAGNESIUM OXIDE P. N. PIYUSH and KANCHAN LATA * Department of Chemistry, B. N. M. V. College, Sahugarh, MADHIPUR (Bihar) INDIA ABSTRACT
More informationChapter 4. Glutamic Acid in Solution - Correlations
Chapter 4 Glutamic Acid in Solution - Correlations 4. Introduction Glutamic acid crystallises from aqueous solution, therefore the study of these molecules in an aqueous environment is necessary to understand
More informationSolutions and Non-Covalent Binding Forces
Chapter 3 Solutions and Non-Covalent Binding Forces 3.1 Solvent and solution properties Molecules stick together using the following forces: dipole-dipole, dipole-induced dipole, hydrogen bond, van der
More information1.1 The Fundamental Chemistry of life
1.1 The Fundamental Chemistry of life Matter makes up everything in the universe, including all living organisms. Matter is composed of elements, a pure substance that cannot be broken down into simpler
More informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
More information16 years ago TODAY (9/11) at 8:46, the first tower was hit at 9:03, the second tower was hit. Lecture 2 (9/11/17)
16 years ago TODAY (9/11) at 8:46, the first tower was hit at 9:03, the second tower was hit By Anthony Quintano - https://www.flickr.com/photos/quintanomedia/15071865580, CC BY 2.0, https://commons.wikimedia.org/w/index.php?curid=38538291
More information1. Introduction to Clusters
1. Introduction to Clusters 1.1 The Field of Clusters Atomic clusters are aggregates of atoms containing from few to a few thousand atoms. Due to their small size, the properties of the clusters are, in
More informationInorganic Pharmaceutical Chemistry
Inorganic Pharmaceutical Chemistry Lecture No. 4 Date :25/10 /2012 Dr. Mohammed Hamed --------------------------------------------------------------------------------------------------------------------------------------
More informationDeveloping Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College
Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model John Dood Hope College What are MD simulations? Model and predict the structure and dynamics of large macromolecules.
More informationSupporting Text Z = 2Γ 2+ + Γ + Γ [1]
Supporting Text RNA folding experiments are typically carried out in a solution containing a mixture of monovalent and divalent ions, usually MgCl 2 and NaCl or KCl. All three species of ions, Mg, M +
More informationMirror Charges, Hydrogen Bonds, And the Water-Dielectric Interface
Mirror Charges, ydrogen Bonds, And the Water-Dielectric Interface Kevin Cahill cahill@unm.edu http://dna.phys.unm.edu/ 1 Field Energy The energy E of an electric field E( x, t) is the space integral of
More informationBIBC 100. Structural Biochemistry
BIBC 100 Structural Biochemistry http://classes.biology.ucsd.edu/bibc100.wi14 Papers- Dialogue with Scientists Questions: Why? How? What? So What? Dialogue Structure to explain function Knowledge Food
More informationSpecific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate
Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate A. Arslanargin, A. Powers, S. Rick, T. Pollard, T. Beck Univ Cincinnati Chemistry Support: NSF, OSC TSRC 2016 November 2, 2016 A. Arslanargin,
More informationE12 UNDERSTANDING CRYSTAL STRUCTURES
E1 UNDERSTANDING CRYSTAL STRUCTURES 1 Introduction In this experiment, the structures of many elements and compounds are rationalized using simple packing models. The pre-work revises and extends the material
More informationOne Q partial negative, the other partial negative Ø H- bonding particularly strong. Abby Carroll 2
Chemistry Notes v Polarity Experiment Ø Things involved Polarity Solubility Dispersion Ø Polarity Shaving cream has soap steric acid Water is polar Food coloring is polar/ionic because dissolved Like dissolves
More informationLec.1 Chemistry Of Water
Lec.1 Chemistry Of Water Biochemistry & Medicine Biochemistry can be defined as the science concerned with the chemical basis of life. Biochemistry can be described as the science concerned with the chemical
More informationDownloaded from UvA-DARE, the institutional repository of the University of Amsterdam (UvA)
Downloaded from UvA-DARE, the institutional repository of the University of Amsterdam (UvA) http://dare.uva.nl/document/351205 File ID 351205 Filename 5: Vibrational dynamics of the bending mode of water
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More information3. Solutions W = N!/(N A!N B!) (3.1) Using Stirling s approximation ln(n!) = NlnN N: ΔS mix = k (N A lnn + N B lnn N A lnn A N B lnn B ) (3.
3. Solutions Many biological processes occur between molecules in aqueous solution. In addition, many protein and nucleic acid molecules adopt three-dimensional structure ( fold ) in aqueous solution.
More informationBasic Chemistry for Biology. Honors Biology
Basic Chemistry for Biology Honors Biology 2013-2014 Composition of Matter Matter - Everything in universe is composed of matter Matter is anything that occupies space or has mass Mass quantity of matter
More informationEquation of State of Dense Helium
Iowa State University From the SelectedWorks of Richard Alan Lesar October 31, 1988 Equation of State of Dense Helium Richard Alan Lesar, Los Alamos National Laboratory Available at: https://works.bepress.com/richard_lesar/27/
More informationPrinciples of Physical Biochemistry
Principles of Physical Biochemistry Kensal E. van Hold e W. Curtis Johnso n P. Shing Ho Preface x i PART 1 MACROMOLECULAR STRUCTURE AND DYNAMICS 1 1 Biological Macromolecules 2 1.1 General Principles
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationAtoms, Molecules, and Life
1 Atoms, Molecules, and Life The Nature of Matter: Atoms and Molecules Elements All matter is composed of chemical elements, substances that cannot be broken down to other substances by ordinary chemical
More informationLecture Presentation. Chapter 11. Liquids and Intermolecular Forces. John D. Bookstaver St. Charles Community College Cottleville, MO
Lecture Presentation Chapter 11 Liquids and Intermolecular Forces John D. Bookstaver St. Charles Community College Cottleville, MO Properties of Gases, Liquids, and Solids State Volume Shape of State Density
More informationSolvation Dynamics in Aqueous Reverse Micelles: A Computer Simulation Study
Solvation Dynamics in Aqueous Reverse Micelles: A Computer Simulation Study James Faeder and Branka M. Ladanyi* Department of Chemistry, Colorado State UniVersity, Fort Collins, Colorado 80523 ReceiVed:
More informationOther Cells. Hormones. Viruses. Toxins. Cell. Bacteria
Other Cells Hormones Viruses Toxins Cell Bacteria ΔH < 0 reaction is exothermic, tells us nothing about the spontaneity of the reaction Δ H > 0 reaction is endothermic, tells us nothing about the spontaneity
More informationTwo-stage Rydberg charge exchange in a strong magnetic field
Two-stage Rydberg charge exchange in a strong magnetic field M. L. Wall, C. S. Norton, and F. Robicheaux Department of Physics, Auburn University, Auburn, Alabama 36849-5311, USA Received 21 June 2005;
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationCHEM 103 Quantum Mechanics and Periodic Trends
CHEM 103 Quantum Mechanics and Periodic Trends Lecture Notes April 11, 2006 Prof. Sevian Agenda Predicting electronic configurations using the QM model Group similarities Interpreting measured properties
More information2. WATER : THE SOLVENT FOR BIOCHEMICAL REACTIONS
2. WATER : THE SOLVENT FOR BIOCHEMICAL REACTIONS 2.1 Water and Polarity Both geometry and properties of molecule determine polarity Electronegativity - The tendency of an atom to attract electrons to itself
More informationDave S. Walker and Geraldine L. Richmond*
J. Phys. Chem. C 2007, 111, 8321-8330 8321 Understanding the Effects of Hydrogen Bonding at the Vapor-Water Interface: Vibrational Sum Frequency Spectroscopy of H 2 O/HOD/D 2 O Mixtures Studied Using Molecular
More informationThe dative covalent bond acts like an ordinary covalent bond when thinking about shape so in NH 4. the shape is tetrahedral
1.3 Bonding Definition Ionic bonding is the electrostatic force of attraction between oppositely charged ions formed by electron transfer. Metal atoms lose electrons to form ve ions. Non-metal atoms gain
More informationNIH Public Access Author Manuscript Acta Chim Slov. Author manuscript; available in PMC 2010 March 11.
NIH Public Access Author Manuscript Published in final edited form as: Acta Chim Slov. 2009 March 11; 56(1): 196 202. Modeling Hofmeister Effects Barbara Hribar-Lee 1,*, Vojko Vlachy 1, and Ken A. Dill
More informationBiochemistry 530: Introduction to Structural Biology. Autumn Quarter 2014 BIOC 530
Biochemistry 530: Introduction to Structural Biology Autumn Quarter 2014 Course Information Course Description Graduate-level discussion of the structure, function, and chemistry of proteins and nucleic
More informationIonic Bonding - Electrostatic Interactions and Polarization
Ionic Bonding - Electrostatic Interactions and Polarization Chemistry 754 Solid State Chemistry Dr. Patrick Woodward Lecture #13 Born-Haber Cycle for NaCl It is energetically unfavorable for Na metal and
More information- intermolecular forces forces that exist between molecules
Chapter 11: Intermolecular Forces, Liquids, and Solids - intermolecular forces forces that exist between molecules 11.1 A Molecular Comparison of Liquids and Solids - gases - average kinetic energy of
More informationThe Logic of Biological Phenomena
The Logic of Biological Phenomena molecules are lifeless yet, in appropriate complexity and number, molecules compose living things What makes living things distinct? they can grow they can move they can
More informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
More informationInteractions. Chemical misconceptions 137. Target level. Topics. Rationale. Resources. Feedback for students
hemical misconceptions 137 nteractions Target level Topics Rationale This probe is intended for students studying chemistry on post-16 courses. hemical bonding (including: ionic, covalent, metallic, polar,
More informationPhase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationAli Eftekhari-Bafrooei and Eric Borguet. Department of Chemistry, Temple University, Philadelphia PA Supporting information
The Effect of Electric Fields on the Ultrafast Vibrational Relaxation of Water at a Charged Solid-Liquid Interface as Probed by Vibrational Sum Frequency Generation Ali Eftekhari-Bafrooei and Eric Borguet
More information74 these states cannot be reliably obtained from experiments. In addition, the barriers between the local minima can also not be obtained reliably fro
73 Chapter 5 Development of Adiabatic Force Field for Polyvinyl Chloride (PVC) and Chlorinated PVC (CPVC) 5.1 Introduction Chlorinated polyvinyl chloride has become an important specialty polymer due to
More informationElectrolytes. Chapter Basics = = 131 2[ ]. (c) From both of the above = = 120 8[
Chapter 1 Electrolytes 1.1 Basics Here we consider species that dissociate into positively and negatively charged species in solution. 1. Consider: 1 H (g) + 1 Cl (g) + ()+ () = { } = (+ )+ ( ) = 167[
More informationOn the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study
On the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study Wynand Dednam and André E. Botha Department of Physics, University of South Africa, P.O.
More informationCalculations of the Oxygen Isotope Fractionation between Hydration Water of Cations and Free Water
Calculations of the Oxygen Isotope Fractionation between Hydration Water of Cations and Free Water P. Bopp, K. Heinzinger, and P. C. Vogel Max-Planck-Institut für Chemie (Otto-Hahn-Institut), Mainz, Germany
More informationColloid Chemistry. La chimica moderna e la sua comunicazione Silvia Gross.
Colloid Chemistry La chimica moderna e la sua comunicazione Silvia Gross Istituto Dipartimento di Scienze di e Scienze Tecnologie Chimiche Molecolari ISTM-CNR, Università Università degli Studi degli Studi
More informationBig Idea #5: The laws of thermodynamics describe the essential role of energy and explain and predict the direction of changes in matter.
KUDs for Unit 6: Chemical Bonding Textbook Reading: Chapters 8 & 9 Big Idea #2: Chemical and physical properties of materials can be explained by the structure and the arrangement of atoms, ion, or molecules
More informationThe potassium channel: Structure, selectivity and diffusion
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 18 8 MAY 2000 The potassium channel: Structure, selectivity and diffusion T. W. Allen a) Protein Dynamics Unit, Department of Chemistry, Australian National
More informationEarth Solid Earth Rocks Minerals Atoms. How to make a mineral from the start of atoms?
Earth Solid Earth Rocks Minerals Atoms How to make a mineral from the start of atoms? Formation of ions Ions excess or deficit of electrons relative to protons Anions net negative charge Cations net
More informationT6.2 Molecular Mechanics
T6.2 Molecular Mechanics We have seen that Benson group additivities are capable of giving heats of formation of molecules with accuracies comparable to those of the best ab initio procedures. However,
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationMODIFIED PROXIMITY CRITERIA FOR THE ANALYSIS OF THE SOLVATION OF A POLYFUNCTIONAL SOLUTE.
Molecular Simulation 1988, Vol. 1 pp.327-332 c 1988 Gordon and Breach Science Publishers S.A. MODIFIED PROXIMITY CRITERIA FOR THE ANALYSIS OF THE SOLVATION OF A POLYFUNCTIONAL SOLUTE. MIHALY MEZEI Department
More information