Proton-Transfer-Driven Water Exchange ... - ACS Publications

Apr 4, 2017 - first solvation (hydration) shell around a sodium ion is only very weakly acidic, ... the simplest dynamical events in the solvation she...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF NEW ENGLAND

Article +

Proton-Transfer-Driven Water Exchange Mechanism in the Na Solvation Shell Matti Hellström, and Jörg Behler J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.7b01490 • Publication Date (Web): 04 Apr 2017 Downloaded from http://pubs.acs.org on April 5, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Proton-Transfer-Driven Water Exchange Mechanism in the Na+ Solvation Shell Matti Hellstr¨om† and J¨org Behler∗,†,‡ Lehrstuhl f¨ ur Theoretische Chemie, Ruhr-Universit¨ at Bochum, 44780 Bochum, Germany, and Institut f¨ ur Physikalische Chemie, Georg-August-Universit¨ at G¨ottingen, Tammannstraße 6, 37077 G¨ ottingen, Germany E-mail: [email protected]



To whom correspondence should be addressed Lehrstuhl f¨ ur Theoretische Chemie, Ruhr-Universit¨at Bochum, 44780 Bochum, Germany ‡ Institut f¨ ur Physikalische Chemie, Georg-August-Universit¨at G¨ ottingen, Tammannstraße 6, 37077 G¨ ottingen, Germany †

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Ligand exchange plays an important role for organic and inorganic chemical reactions. We demonstrate the existence of a novel water exchange mechanism, the “proton transfer pathway” (PTP), around Na+ (aq) in basic (high pH) solution, using reactive molecular dynamics simulations employing a high-dimensional neural network potential. An aqua ligand in the first solvation (hydration) shell around a sodium ion is only very weakly acidic, but if a hydroxide ion is present in the second solvation shell, thermal fluctuations can cause the aqua ligand to transfer a proton to the neighboring OH– , resulting in a transient direct-contact ion pair, Na+ −OH– , which is only weakly bound and easily dissociates. The extent to which water exchange events follow the PTP is pH-dependent: in dilute NaOH(aq) solutions, only very few exchanges do, whereas in saturated NaOH(aq) solutions up to a third of water self-exchange events are induced by proton transfer. The principles and results outlined here are expected to be relevant for chemical synthesis involving bases and alkali metal cations.

Introduction The solvation structure and solvation dynamics around ions in water and in other solvents play important roles in biochemistry, 1–3 electrochemistry, 4 and geochemistry, 5–7 especially during chemical reactions where the solvation shell may experience strong distortions or even break and reform. One of the simplest dynamical events in the solvation shell is water selfexchange around a metal cation in aqueous solution. 8,9 Depending on the chemical element and oxidation state of the cation, the rate of water exchange can span tens of orders of magnitude. Water exchange around alkali metal cations like Na+ , which is the focus of the present work, is quite fast; in fact, it is so fast that the standard nuclear magnetic resonance experimental techniques used to measure exchange rates give quite unreliable results. 8 However, estimating the exchange rate is possible using molecular dynamics (MD) simulations, which have given estimates of the mean residence time for an aqua ligand around Na+ to be between 10-40 ps. 10–14 Exchange mechanisms are frequently classified as associative, dissociative, or interchange, depending on whether the incoming ligand enters the solvation shell before, after, or at the same time as, the leaving ligand leaves the solvation shell. 15 MD simulations have suggested that the exchange around Na+ (aq) is mostly associative, 12,16,17 but there are also other ways to classify ligand exchange mechanisms. For example, the pentaamminechloridocobalt(III) coordination complex, [CoCl(NH3 )5 ]2+ , can exchange its chlorido ligand for a hydroxido ligand via the “Sn1CB” mechanism, where the Cl– leaves as a result of an outside hydroxide ion deprotonating one of the ammine ligands. 18 After a few more reaction steps, the end result is a [Co(NH3 )5 OH]2+ complex, where OH– has substituted the Cl– . In such a mechanism, the deprotonation of one ligand promotes the exchange of another ligand. A similar phenomenon can be seen for many aqua complexes; for example, the rate of exchange of aqua ligands in [Al(OH)(OH2 )5 ]2+ is greater by a factor of 104 as compared to [Al(OH2 )6 ]3+ , 19 which makes the deprotonation of one of the aqua ligands the rate-limiting step for water exchange in [Al(OH2 )6 ]3+ complexes. 20 OH– ligands can thus have a labilizing effect on the OH2 ligands around the same cation. 2

ACS Paragon Plus Environment

Page 2 of 17

Page 3 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In the present work, we will show that the deprotonation of an OH2 ligand around Na (aq), resulting in an OH– ligand, influences the water exchange rate, but in a very different way as compared to the previous AlIII example. We show that around Na+ , the OH– ligand itself is more labile than the OH2 ligands. In dilute sodium hydroxide solutions, direct-contact Na+ −OH– ion pairs are very rare, 21–24 although at very high concentrations, for example near the room-temperature solubility limit of 19 mol L−1 , MD simulations have indicated that Na+ on average coordinates to between 1 and 2 OH– ions in the first solvation shell (FSS). 23,25 At low concentrations, Na+ does not manage to “trap” the OH– in the FSS. However, MD simulations have suggested that OH– ions in dilute NaOH(aq) solutions are frequently found in the second solvation shell (SSS) around the Na+ . 26–28 Hydroxide ions can participate in proton-transfer (PT) events with neighboring water molecules, i.e., HOα H · · · Oβ H− −−→ HOα− · · · HOβ H, (1) +

where the superscripts α and β are used to distinguish between different O atoms. This feature of OH– ions is responsible for their high diffusion coefficients in water. 29 At low concentration, such PT events have been shown to be coupled to specific hydrogen-bonding motifs around the OH– ion; 30 OH– normally accepts four hydrogen bonds from surrounding water molecules, but thermal fluctuations occasionally cause one of the hydrogen bonds to break, so that only three hydrogen bonds remain. This less common hydrogen-bonding environment is associated with a higher rate of PT. Recently, we showed that at higher concentrations, hydrogen-bond fluctuations also around the H2 O molecules significantly influence the PT rate. 28 If the OH– is in the SSS around Na+ , it might happen that the proton-donating original H2 O molecule is an aqua ligand around Na+ . Such a PT reaction would result in a directcontact Na+ −OH– ion pair: + β α − −− Na+ −Oα HH · · · Oβ H− ↽ −⇀ − Na −O H · · · HO H.

(2)

Such direct-contact ion pairs rarely form at low concentrations of NaOH(aq), because the PT barrier is greater for OH2 molecules that are coordinated to Na+ than for those that are “free” in the solution. 28 However, as we mentioned above, at very high concentrations, it is not uncommon to find OH– coordinated to Na+ . The direct-contact ion pair can decompose in one of two ways: (i) proton-transfer, where the hydroxido ligand becomes an aqua ligand by accepting a proton from a neighboring water molecule, i.e., the backward equation 2, or (ii) diffusion, where the OH– leaves the FSS without participating in PT, giving room for some OH2 to take its place in the FSS. We will show that Na+ −OH– bonds are shorter-lived than Na+ −OH2 bonds, i.e., that OH– more readily leaves the FSS than OH2 does (OH– is more labile than OH2 ). This turns out to be the key ingredient for a novel water self exchange mechanism, the “proton transfer pathway”, that we present in this Article for Na+ in basic aqueous solution. With increasing concentration of NaOH(aq), diffusion of hydroxide away from the Na+ turns out to play an increasingly important role for water exchange. That is, a significant fraction of water exchange events take place by the formation of a transient/intermediate direct-

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 4 of 17

Page 5 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

+

Na −OH2 (FSS) aqua direct kleave

aqua kdeprot , −H+ hydrox kprot , + H+

hydrox proton transfer pathway kleave

Na+ · · ·OH– (SSS)

Na+ · · ·OH2 (SSS) a The

Na+ −OH– (FSS)

“direct” and “proton transfer pathway” labels refer only to original OH2 ligands (bold) entering the SSS. For original OH– ligands, the two labels would change places.

Scheme 1: Possible reaction pathways for an aqua ligand in the first solvation shell [FSS; simple dash (−)] around Na+ to enter the second solvation shell [SSS; three dots (· · ·)], together with the corresponding effective rate constants.a depend on the NaOH concentration. The formed Na+ −OH– can then reform Na+ −OH2 hydrox hydrox , or dissociate with a rate constant kleave . In the latter case, the with a rate constant kprot original OH2 ligand has left the FSS via the proton transfer pathway. hydrox hydrox aqua aqua , for different , kprot , and kleave , kdeprot The challenge is to find the rate constants kleave concentrations of NaOH(aq). To this end, we employ a high-dimensional neural network potential 31–34 that we specifically developed to run molecular dynamics (MD) simulations of NaOH(aq) solutions with dispersion-corrected density functional theory quality over the entire NaOH room-temperature solubility range. 25,28 The potential is computationally inexpensive, and allows for proton-transfer events between H2 O and OH– . We have previously validated it 25,28 against x-ray and neutron diffraction data for NaOH solutions, 21–23 as well as results from other molecular dynamics simulations, 23,26,27,35,36 and used it to show, for example, that there is great variation in the structure and coordination numbers for Na+ ions in NaOH solutions of different concentrations. 25 In MD studies on ligand exchange reactions around metal cations, it is often necessary to impose some kind of constraint on the exchanging ligands, for example their distances to the central ion, so as to sample the entire free-energy landscape. 37 However, because both ligand-exchange events around Na+ and PT events involving OH– are quite frequent in NaOH(aq) solutions, we can run completely unbiased molecular dynamics simulations and let the system dynamics evolve naturally. From N V T simulations at 300 K (with about 1100–1500 atoms depending on the concentration) we extracted 20 different starting points per concentration separated by at least 100 ps, and continued from those starting points (positions and velocities) for 1 ns in the N V E ensemble, thus sampling the average dynamics for solutions at constant volume and temperature (300 K). The volumes (densities) had been pre-equilibrated in the N P T ensemble (at 1 bar, 300 K); see also our previous works. 25,28 We use mole fractions xNaOH to identify the concentrations of solutions, i.e., xNaOH =

NNaOH , NNaOH + NH2 O

(3)

where NNaOH and NH2 O are the numbers of NaOH and H2 O formula units in the periodically repeated simulation box. The room-temperature solubility of NaOH in water is about xNaOH = 0.33, i.e., there are two formula units H2 O for every formula unit NaOH. Here, we characterize the reactions in Scheme 1 for xNaOH = 0.016, 0.049, 0.085, 0.143, 0.185, 0.255, 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 6 of 17

Page 7 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the FSS, intermediate, and SSS regions marked with different colors. We consider OH2 and OH– to be in the FSS if the distance Na−O < 2.7 ˚ A, and in the SSS if the distance is between ˚ ˚ 3.6 A and 5.0 A (the upper bound for the SSS has no impact on the results presented in this paper). The two smaller distance criterions were chosen because they encompass the first minimum in the Na−O radial distribution function with a margin of at least 0.2 ˚ A for all the studied concentrations of NaOH(aq). As shown in the snapshot in Figure 2, at xNaOH = 0.016, a typical Na+ is surrounded by 5 OH2 in the FSS, 1 OH2 in the TS region, and 13 H2 O + 1 OH– in the SSS. Similarly, for PT reactions, we also calculate SSP time correlation functions, where “stable” water molecules have two H within 1.1 ˚ A of the central O, and “stable” hydroxide ions ˚ have only one H within 1.4 A of the O. We calculate four different “basic” time correlation functions for each xNaOH : the ligandexchange and proton-transfer functions when the ligand is a stable OH2 in the FSS at t0 , and the same functions when the ligand is a stable OH– in the FSS at t0 . We only apply the conditions at the time origin; thus for the ligand-exchange function, the ligand can participate in an arbitrary number of PT events before entering the SSS, and for the protontransfer function, the ligand may leave and reenter the FSS an arbitrary number of times before participating in PT. This introduces some errors in the calculated reaction rates, that we attempt to correct for (see below). We limit the discussion to molecules that at t0 are in the FSS of only a single Na+ ion. In highly concentrated solutions, some fraction of OH2 molecules and OH– ions bridge two or more Na+ ions, 25 but we do not consider the kinetics of such bridging ligands here, as there may be other, more complicated, mechanisms for these structures.

Results and Discussion Figure 3 shows the calculated SSP time correlation functions for the PT reactions and ligandexchange reactions at xNaOH = 0.085. In Figure 3a, the ligand is a “stable” OH2 at t0 . At this concentration, clearly, the time it takes to deprotonate OH2 is much longer than the time it takes to leave the Na+ . Therefore, many of the OH2 will have already left the Na+ before the PT event takes place. To calculate the rate of deprotonation while the ligand aqua ), we therefore consider a third time correlation function, that is still bound to Na+ (kdeprot decays for OH2 in the FSS that either deprotonates or enters the SSS (thick blue line). This allows us to extract a better estimate of the time scale of the slower of the two possible −1 −1 −1 . In Figure 3a, the “raw” lifetime calculated for the − Rτfast ≈ τeither reactions, kslow = τslow ∞ slower reaction (i.e., PT) is τraw = 0 C(t) dt, while the “modified” (better) one is calculated aqua −1 −1 −1 −1 . ) . The effective rate constant kdeprot = τmod − τfast as τmod = (τeither Using this approach, the time scale associated with PT from H2 O to OH– in the FSS increases from the “raw” estimate of 40 ps to 54 ps. The modified τ is greater, because the barrier for PT is greater when the H2 O is in the FSS of Na+ than when it is free in the solution, 28 and the raw estimate included some H2 O that had left the SSS and then participated in the PT with a lower “free-in-solution”-PT-barrier. For OH– (Figure 3b), both the PT reaction and the ligand exchange event are faster than the corresponding reactions for H2 O (Figure 3a). However, the relative timescales involved 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 8 of 17

Page 9 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 1: Calculated effective rate constants k (in ps−1 ) at 300 K and at different mole fractions xNaOH of NaOH(aq) for aqua and hydroxido ligands in the first FSS solvation shell (FSS) around a single Na+ ion, together with the fraction ξOH − of – a OH ligands in the FSS. FSS xNaOH ξOH − 0.016 0.0018(1) 0.049 0.0066(2) 0.085 0.0155(5) 0.143 0.046(1) 0.185 0.086(1) 0.255 0.197(2) 0.333 0.368(4) a The uncertainties are calculated over

aqua aqua kleave kdeprot 0.100(3) 0.006(1) 0.092(4) 0.013(2) 0.081(3) 0.019(2) 0.066(3) 0.024(1) 0.051(3) 0.026(1) 0.034(2) 0.031(1) 0.019(2) 0.032(2) 20 independent N V E

hydrox hydrox kprot kleave 1.34(10) 0.29(6) 0.94(4) 0.24(2) 0.64(2) 0.22(1) 0.35(2) 0.157(6) 0.22(1) 0.111(7) 0.123(7) 0.066(4) 0.070(6) 0.034(3) trajectories as 2.09×(standard error).

OH– in the solution that can accept the proton. Here, we reiterate that the rate constant aqua kdeprot includes the effect of OH– concentration, and thus describes the overall rate, per OH2 ligand in the FSS, of the reaction. Scheme 1 gives rise to a set of coupled kinetic equations that can be solved analytically for different initial conditions (see the Supporting Information for the full derivation of these solutions). For example, by setting the initial population of Na+ −OH2 to 1, and the population of the other species to 0, the populations of Na+ ···OH2 and Na+ ···OH– as t → ∞ then correspond to the fractions of original OH2 in the FSS that have left via the direct and proton-transfer pathways, respectively. Figure 4a shows the fate of an initial OH2 in the FSS as a function of time at xNaOH aqua = 0.081 ps−1 ≫ = 0.085, based on the calculated rate constants in Table 1. Because kleave aqua = 0.019 ps−1 , most aqua ligands will simply enter the SSS (pink dashed line) withkdeprot out participating in PT. Even for the few that do participate in PT events, most of the formed hydroxido ligands will quickly participate in another PT event and reform the origihydrox hydrox = 0.22 ps−1 , At this nal aqua ligand rather than leaving, since kprot = 0.64 ps−1 ≫ kleave concentration, only five percent (orange dashed line) of original aqua ligands leave as OH– . We also consider the fate of hydroxido ligands, and we denote the fraction of such ligands FSS FSS FSS FSS FSS FSS FSS refer to in the FSS as ξOH = NOH − , i.e., ξOH− − /[NOH− + NOH ]. Here, NOH− and NOH 2 2 – the number of “stable” OH and OH2 (as defined for the SSP time correlation functions) FSS in the FSS. ξOH − is very small at low concentrations (Table 1), which makes exchange FSS events via the proton transfer pathway unlikely, but ξOH − increases rapidly with increasing concentration and the proton transfer pathway may therefore play a more important role FSS at higher concentrations. At xNaOH = 0.085, ξOH − = 0.016, which means that only 1.6% of ligands that can be described as either stable OH– or stable OH2 are OH– , and that the rest, 98.4%, are OH2 . As shown in Figure 4b, the OH– ligand will likely participate in a PT event, which can be seen as a signficant peak in the Na+ −OH2 population at about t = 4 ps (purple line), and the OH2 can then either leave as OH2 or as OH– as discussed above. Because PT is so fast for OH– in the FSS, only 29% of ligands that are OH– leave as hydrox OH– , despite the fact that breaking the Na+ −OH– bond is much easier (kleave = 0.22 ps−1 ) aqua than breaking the Na+ −OH2 bond (kleave = 0.081 ps−1 ) that results from PT. Figure 4c 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 10 of 17

Page 11 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

participate in PT (the population of Na+ −OH– approaches 0.2 at about t = 15 ps in Figure 4d), and, consequently, that more of the initial OH2 ligands leave as OH– . In fact, 35%, or about one third, of OH2 ligands follow the proton transfer pathway and leave as OH– . For ligands that are initially OH– , most (56%) will also leave as OH– , unlike the case at xNaOH = 0.085 where 71% of OH– ligands left as OH2 because of the fast PT reaction. The weighted averages in Figure 4f reveal that in a saturated solution, about 57% of ligands leave as OH2 , and the rest as OH– . Figure 4d–e also show the total population of ligands leaving (dashed black line, which is the sum of the Na+ · · · OH2 and Na+ · · · OH– populations), together with the population of ligands that would have left had only the “direct”, and no PT, pathway been available hydrox aqua ·t) ·t) for original OH2 in Figure 4d, and 1−exp(−kleave (dashed green line; i.e., 1−exp(−kleave for original OH– in Figure 4e). As a result of the PT pathway, the total rate of water exchange increases (because the black dashed line is above the green dashed line in Figure 4d), whereas the total rate of hydroxido exchange decreases, where the rate of “exchange” refers only to breaking the Na+ −O bond through diffusion, i.e., the rate at which an O atom in the FSS enters the SSS. In Figure 5, we divide all exchange events into four categories: H2 O leaving as H2 O (pink solid bar), OH– leaving as H2 O (pink striped bar), H2 O leaving as OH– (orange striped bar), and OH– leaving as OH– (orange solid bar). Of course, there can be an arbitrary number of PT events before a ligand exchange. Therefore, for a ligand leaving as, for example, H2 O, we use the relative time periods that the ligand spent as OH– and as H2 O to weight the contribution to either of the “OH– leaving as H2 O” and “H2 O leaving as H2 O” categories. Interestingly, for all xNaOH , the number of ligands leaving as H2 O (sum of pink solid and FSS striped bars) is smaller than the average fraction of H2 O ligands (1 − ξOH − , black line and points). This means that the PT mechanism gives a net increase in the rate of H2 O ligands leaving (as already shown for xNaOH = 0.333 in Figure 4d); in contrast, it gives the corresponding net decrease for OH– ligands leaving. FSS At all concentrations, most ligands are OH2 (ξOH − < 0.5, Table 1). At low concentration, most of the OH2 ligands also leave as OH2 (solid pink in Figure 5), because the probability aqua is small, Table 1). However, with increasing concentration, more for PT is very low (kdeprot and more of the OH2 ligands leave as OH– (brown striped bar); indeed, at the highest xNaOH = 0.333, the number of OH2 that leave as OH– is about half of those that leave as OH2 , i.e., about one third of water exchange events take place via the proton-transfer pathway. The OH– ligands, on the other hand, are, independently of the concentration, roughly equally likely to leave as OH– (orange solid bar) as they are to leave as OH2 (pink dashed bar). In some ways, one can consider the OH– ions to “catalyze” water self-exchange, since they provide an alternative reaction mechanism, the proton-transfer pathway, which increases the total rate of water exchange. However, even in the saturated solution, the majority of OH2 ligands still follow the “direct” pathway, possibly because there are not enough OH– ions in the solution (as compared to the number of OH2 ligands).

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 12 of 17

Page 13 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The two kinds of pathways (proton-transfer and direct) in the present article differ from the more common distinction between associative, dissociative, and interchange mechanisms. To what relative extent the direct and proton-transfer pathways proceed via associative, dissociative, and interchange mechanisms, remains an open question, as does whether there are substantial differences for ligands that bridge multiple Na+ ions as compared to those that only coordinate to a single Na+ ion. The present work illustrates how the presence of hydroxide ions, i.e., a high pH value, influences the rate of water exchange around one of the most common ions in water, Na+ (aq).

Acknowledgement This work was supported by the Cluster of Excellence RESOLV (EXC 1069) funded by the Deutsche Forschungsgemeinschaft, the DFG Heisenberg fellowship Be3264/6-1, and the DFG Heisenberg professorship Be3264/11-1.

Supporting Information Available Derivation of the solutions to the kinetic equations arising from Scheme 1, tabulated data for Figure 5, and snapshots (cartesian coordinates) for all xNaOH . This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Sigel, A., Sigel, H., Sigel, R., Eds. The Alkali Metal Ions: Their Role for Life; Metal Ions in Life Sciences; Springer, 2016. (2) Lambry, J.-C.; Stranava, M.; Lobato, L.; Martinkova, M.; Shimizu, T.; Liebl, U.; Vos, M. H. Ultrafast Spectroscopy Evidence for Picosecond Ligand Exchange at the Binding Site of a Heme Protein: Heme-Based Sensor YddV. J. Phys. Chem. Lett. 2016, 7, 69–74. (3) Stachura, M.; Chakraborty, S.; Gottberg, A.; Ruckthong, L.; Pecoraro, V. L.; Hemmingsen, L. Direct Observation of Nanosecond Water Exchange Dynamics at a Protein Metal Site. J. Am. Chem. Soc. 2017, 139, 79–82. (4) Schmickler, W.; Santos, E. Interfacial Electrochemistry, 2nd ed.; Springer: Berlin, 2010. (5) Stack, A. G.; Raiteri, P.; Gale, J. D. Accurate Rates of the Complex Mechanisms for Growth and Dissolution of Minerals Using a Combination of Rare-Event Theories. J. Am. Chem. Soc. 2012, 134, 11–14. (6) Hofmann, A. E.; Bourg, I. C.; DePaolo, D. J. Ion desolvation as a mechanism for kinetic isotope fractionation in aqueous systems. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 18689–18694.

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(7) Loganathan, N.; Yazaydin, A. O.; Bowers, G. M.; Kalinichev, A. G.; Kirkpatrick, R. J. Cation and Water Structure, Dynamics, and Energetics in Smectite Clays: A Molecular Dynamics Study of Ca–Hectorite. J. Phys. Chem. C 2016, 120, 12429–12439. (8) Helm, L.; Merbach, A. E. Inorganic and Bioinorganic Solvent Exchange Mechanisms. Chem. Rev. 2005, 105, 1923–1960. (9) Richens, D. T. Ligand Substitution Reactions at Inorganic Centers. Chem. Rev. 2005, 105, 1961–2002. (10) Rey, R.; Hynes, J. T. Hydration Shell Exchange Kinetics: An MD Study for Na+ (aq). J. Phys. Chem. 1996, 100, 5611–5615. (11) Koneshan, S.; Rasaiah, J. C.; Lynden-Bell, R. M.; Lee, S. H. Solvent Structure, Dynamics, and Ion Mobility in Aqueous Solutions at 25 ◦ C. J. Phys. Chem. B 1998, 102, 4193–4204. (12) Kerisit, S.; Rosso, K. M. Transition path sampling of water exchange rates and mechanisms around aqueous ions. J. Chem. Phys. 2009, 131, 114512. (13) Bankura, A.; Carnevale, V.; Klein, M. L. Hydration structure of salt solutions from ab initio molecular dynamics. J. Chem. Phys. 2013, 138, 014501. (14) Kelley, M.; Donley, A.; Clark, S.; Clark, A. Structure and Dynamics of NaCl Ion Pairing in Solutions of Water and Methanol. J. Phys. Chem. B 2015, 119, 15652–15661. (15) Langford, C. H.; Gray, H. B. Ligand Substitution Processes; W. A. Benjamin, Inc.: New York, 1966. (16) Sp˚ angberg, D.; Wojcik, M.; Hermansson, K. Pressure dependence and activation volume for the water exchange mechanism in NaCl(aq) from MD simulations. Chem. Phys. Lett. 1997, 276, 114 – 121. (17) Hermansson, K.; Wojcik, M. Water Exchange around Li+ and Na+ in LiCl(aq) and NaCl(aq) from MD Simulations. J. Phys. Chem. B 1998, 102, 6089–6097. (18) Pearson, R. G.; Basolo, F. Mechanism of Substitution Reactions of Complex Ions. X. πBonding in Dissociation Reactions of Octahedral Complexes. J. Am. Chem. Soc. 1956, 78, 4878–4883. (19) Nordin, J. P.; Sullivan, D. J.; Phillips, B. L.; Casey, W. H. An 17 O-NMR Study of the Exchange of Water on AlOH(H2 O)2+ 5 (aq). Inorg. Chem. 1998, 37, 4760–4763. (20) Swaddle, T. W.; Rosenqvist, J.; Yu, P.; Bylaska, E.; Phillips, B. L.; Casey, W. H. Kinetic Evidence for Five-Coordination in AlOH(aq)2+ Ion. Science 2005, 308, 1450–1453. (21) Botti, A.; Bruni, F.; Imberti, S.; Ricci, M. A.; Soper, A. K. Ions in water: The microscopic structure of concentrated NaOH solutions. J. Chem. Phys. 2004, 120, 10154– 10162. 14

ACS Paragon Plus Environment

Page 14 of 17

Page 15 of 17

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(22) Imberti, S.; Botti, A.; Bruni, F.; Cappa, G.; Ricci, M. A.; Soper, A. K. Ions in water: The microscopic structure of concentrated hydroxide solutions. J. Chem. Phys. 2005, 122, 194509. (23) Megyes, T.; B´alint, S.; Gr´osz, T.; Radnai, T.; Bak´o, I.; Sipos, P. The structure of aqueous sodium hydroxide solutions: A combined solution x-ray diffraction and simulation study. J. Chem. Phys. 2008, 128, 044501. (24) Aziz, E. F.; Eisebitt, S.; Eberhardt, W.; Cwiklik, L.; Jungwirth, P. Existence of Oriented Ion–Hydroxide Clusters in Concentrated Aqueous NaCl Solution at pH 13. J. Phys. Chem. B 2008, 112, 1262–1266. (25) Hellstr¨om, M.; Behler, J. Structure of aqueous NaOH solutions: insights from neuralnetwork-based molecular dynamics simulations. Phys. Chem. Chem. Phys. 2017, 19, 82–96. (26) Chen, B.; Park, J. M.; Ivanov, I.; Tabacchi, G.; Klein, M. L.; Parrinello, M. FirstPrinciples Study of Aqueous Hydroxide Solutions. J. Am. Chem. Soc. 2002, 124, 8534– 8535. (27) Chen, B.; Ivanov, I.; Park, J. M.; Parrinello, M.; Klein, M. L. Solvation Structure and Mobility Mechanism of OH− : A Car-Parrinello Molecular Dynamics Investigation of Alkaline Solutions. J. Phys. Chem. B 2002, 106, 12006–12016. (28) Hellstr¨om, M.; Behler, J. Concentration-Dependent Proton Transfer Mechanisms in Aqueous NaOH Solutions: From Acceptor-Driven to Donor-Driven and Back. J. Phys. Chem. Lett. 2016, 7, 3302–3306. (29) Agmon, N.; Bakker, H. J.; Campen, R. K.; Henchman, R. H.; Pohl, P.; Roke, S.; Th¨amer, M.; Hassanali, A. Protons and Hydroxide Ions in Aqueous Systems. Chem. Rev. 2016, 116, 7642–7672. (30) Tuckerman, M. E.; Marx, D.; Parrinello, M. The nature and transport mechanism of hydrated hydroxide ions in aqueous solution. Nature 2002, 417, 925–929. (31) Behler, J.; Parrinello, M. Generalized Neural-Network Representation of HighDimensional Potential-Energy Surfaces. Phys. Rev. Lett. 2007, 98, 146401. (32) Behler, J. Atom-centered symmetry functions for constructing high-dimensional neural network potentials. J. Chem. Phys. 2011, 134, 074106. (33) Behler, J. Neural network potential-energy surfaces in chemistry: a tool for large-scale simulations. Phys. Chem. Chem. Phys. 2011, 13, 17930–17955. (34) Behler, J. Representing potential energy surfaces by high-dimensional neural network potentials. J. Phys.: Condens. Matter 2014, 26, 183001. (35) White, J. A.; Schwegler, E.; Galli, G.; Gygi, F. The solvation of Na+ in water: Firstprinciples simulations. J. Chem. Phys. 2000, 113, 4668–4673. 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(36) Bankura, A.; Carnevale, V.; Klein, M. L. Hydration structure of Na+ and K+ from ab initio molecular dynamics based on modern density functional theory. Mol. Phys. 2014, 112, 1448–1456. (37) Annapureddy, H. V. R.; Dang, L. X. Understanding the Rates and Molecular Mechanism of Water-Exchange around Aqueous Ions Using Molecular Simulations. J. Phys. Chem. B 2014, 118, 8917–8927. (38) Laage, D.; Hynes, J. T. On the Residence Time for Water in a Solute Hydration Shell: Application to Aqueous Halide Solutions. J. Phys. Chem. B 2008, 112, 7697–7701. (39) Kumar, P.; Kulkarni, A. D.; Yashonath, S. Influence of a Counterion on the Ion Atmosphere of an Anion: A Molecular Dynamics Study of LiX and CsX (X = F− , Cl− , I− ) in Methanol. J. Phys. Chem. B 2015, 119, 10921–10933. (40) Saha, D.; Supekar, S.; Mukherjee, A. Distribution of Residence Time of Water around DNA Base Pairs: Governing Factors and the Origin of Heterogeneity. J. Phys. Chem. B 2015, 119, 11371–11381.

16

ACS Paragon Plus Environment

Page 16 of 17

Page 17 of 17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment