Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

Conspectus The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here, we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F– over Cl– in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

This paper reviews local hydration structures of monoand divalent cations and assesses the concept of hydration mimicry for the rapid transport of specific ions through ion channels. Cluster QCT and surface potentials that provide hydration free energies are also reviewed.

INTRODUCTION
This Account describes recent research at the intersection of the topics of ion−water clusters, 5,6 the theory of solutions, 1,7,8 specific ion effects, 9−11 and the selectivity of membrane ion channels. 4 We focus on anions in water because of their central position in classic specific ion effects, so-called Hofmeister effects. 10 The anions considered here have been challenges for the molecular quasi-chemical theory (QCT), 1,7,8,12 which is the most advanced theory available to address hydration while accounting for chemical-level interactions and because recent theoretical progress on those challenges seems decisive. 2,3,13,14 We include HO − (aq) in this discussion because of its centrality in aqueous solution chemistry and the continued theory and simulation interest in this ion 15−17 and because the new results sharpen our understanding of the hydration of that ion.

Context: Selectivity of Membrane Ion Channels
As active components of nearly half of all proteins, 18,19 ions can bind to proteins and stabilize conformational states required for biological function and can participate in enzyme catalysis. 20 As an example, K + ions bind to membrane channel proteins and stabilize functional conformations, thereby catalyzing the permeation of K + across cellular membranes while also rejecting other ions such as Na + . 21 Selective ion transport plays an important role in numerous physiological functions, including electrical signaling and cell volume control. Loss of ion selectivity, or blocking of ion transport, can have either catastrophic or beneficial effects. For example, the loss of selective conduction of K + over Na + by potassium channels in cardiac muscle interferes with the termination of action potentials, which can lead to lifethreatening heart arrhythmia. 21 In beneficial cases of blocked ion transport, drugs that block specific channels hold promise for treating neurological disorders, autoimmune diseases, and cancers. 22, 23 Peptide toxins from several poisonous animals exemplify detrimental possibilities of blocked ion transport. 24 Indeed, simple divalent metal ions can be potent channel blockers by getting trapped in the channel, and both monovalent and divalent ions permeate selectively. Thus, understanding the mechanisms of specific ion binding and transport in proteins is important for understanding protein function critical to health, disease, and therapeutic development. 20 For anions such as Cl − , the regulation of ion concentration is achieved through membrane transport proteins of the CLC family and others, including channelrhodopsins, which facilitate the passage of Cl − through electrochemical potential gradients. 25−27 While those structurally diverse proteins select for Cl − , other channel proteins discriminate against Cl − . An interesting example is FLUC, a family of fluoride-specific ion channels with dual-topology architecture. 28 These channels display an astonishingly high selectivity of 10 4 for F − over Cl − despite their identical charge and their size similarity. Understanding such mechanisms for selectivity in ion transport has been a target for many modeling and simulation studies of high variety, emphasizing the K + /Na + selectivity of potassium ion channels; for example, see refs 29−38. However, computational studies of anion transport mechanisms demonstrated by chloride-selective channels 39−42 that address comparisons to alternatives such as FLUC are less mature. 3,43 A conceptually natural strategy for addressing ion selectivity with computation would be direct simulations controlling for the contrasting cases or perhaps a clear, quantitative theory that permits controlling for mechanistic features of the ion binding or transport. Both of these requests are difficult. Here, we work toward the second of these alternatives by building a statistical molecular theory of ion binding thermodynamics.

Free Energies of Binding of Hydrated Ions Span a Chemical Scale of Energies
QCT 1,7,47 aims to evaluate interaction free energies of ions in s o l u t i o n 4 8 ( F i g u r e 1 ) a n d p r o t e i n b i n d i n g sites. 1,4,29,30,32,33,35,46,47,49−55 That interaction free energy, or excess chemical potential, is obtained from the full chemical potential μ X less the indicated ideal contribution. Here, T is the temperature, R is the molar gas constant, ρ X is the number density of the ion of interest, and X [ ] is the canonical partition function of a molecule X in volume V. 4,56 These free energies�single-ion activities when X is an ion� are knowable and appropriate targets for computation, but they are not measured solely on the basis of classic thermodynamics. Thus, widely available tabulations�Figure 1 utilizes one such tabulation�adopt extrathermodynamic assumptions.
Observe (Figure 1) that the free energies span a chemical energy scale much larger than RT ≈ 0.6 kcal/mol at room temperature. For example, the hydration free energy of Be 2+ is about 1000RT. Clustered below −400 kcal/mol are values for divalent transition-metal ions. 44,51 Including the aqueous ferric ion, Fe 3+ (aq), for which QCT performs satisfactorily, 57 would require expanding the range of Figure 1 by another factor 2. Thus, though the now-canonical van der Waals perspective on liquids is an appropriate definition of the statistical mechanical problem for treating liquids generically, it immediately emphasizes <1%-magnitude free-energy effects. In contrast, the basic concept of QCT is to treat an ion together with innershell partners as an individual molecular species. 1,7,12,57−59 Interactions of typical ions with inner-shell partners are chemical in nature, molecularly intricate, and intense on a thermal scale. The concept of "inner-shell" partners is central to QCT. It identifies near neighbors of a targeted species and will be discussed later for the present applications.
Molecular quasi-chemical theory (QCT) 1,7,8,59 was developed with the explicit goal of including chemical-level interactions within a molecular statistical thermodynamic theory. Initial applications were simple and remarkably accurate. 58 That success obviates a canonical "molecular forcefield fitting → molecular simulation" workflow in the study of liquids.
1.2.1. Direct QCT. The basic status of QCT may be supported by the fact that QCT itself can be implemented through molecular simulation calculations. 8 That approach is termed "direct" QCT. 60 In the direct approach, we acknowledge and exploit the spatial dependence of solute−solvent interaction strengths. The short-range interactions can be chemically involved, and for a suitable choice of the inner shell, the longrange interactions admit a Gaussian statistical model. This direct QCT approach parses the hydration free energy into physically meaningful and computationally well-defined chemical, packing, and long-range nonspecific contributions, thereby becoming a framework for conceptualizing molecular solutions. 12 Direct QCT works naturally with common simulation packages based on either empirical classical or ab initio force fields. 8 On that simulation basis, QCT provides a compelling molecular theory of liquid water itself. 61−63 The direct QCT approach enabled the first direct calculation of the hydration free energy of a protein. 64 Subsequent studies have highlighted the limitations of additive models of free energies that are de rigueur in biophysical speculation and recently led to transformative insights into decades-old assumptions about hydrophobic hydration in proteins. 65 Nevertheless, the initial motivation was the exploitation of molecular electronic structure calculations within statistical thermodynamic modeling. 57 That approach is called "cluster" QCT, wherein the chemical contribution noted above is related to physical solute−solvent clusters. The connection between direct and cluster approaches has been deliberately discussed elsewhere. 1  for the free energies that we seek. This is exact statistical thermodynamics 1 and provides a foundation for mixedresolution approaches, such as QM/MM. We will discuss these three terms in turn. The first term on the right of eq 2 is the inner-shell contribution. It is obtained by studying n water molecules clustering with the X species of interest 2 but without the exterior solution. This contribution involves the equilibrium constant, K n (0) , discussed in section 2.1 below. By utilizing the solution density of the ligands, H O 2 , this term properly assesses the availability of the water molecule ligands.
The rightmost term of eq 2 is the outer-shell contribution. This term involves the hydration of the identified (H 2 O) n X cluster, addressing interactions of the cluster with the exterior solution environment. Previous applications of QCT have utilized the polarizable continuum model (PCM) 67 for this task. (See refs 68 and 69 as examples.) PCM has been incorporated into standard electronic structure codes and focuses on the interactions with the solution at long range. Nevertheless, it is approximate on a molecular scale, and a symptom of that approximate character is the sensitivity of PCM results to radii parameters that are required. We will discuss that issue further below.
The remaining term in eq 2 involves the probability p X (n) that n water molecules contact a distinguished X during its physical motion in solution. This term describes the polydispersity of the populations of an X inner shell; if only one size, say n, were possible, then ln p X (n) = 0. The determination of p X (n) requires the adoption of a proximity criterion describing how a water molecule contacts an X. 3,13,14,70 That polydispersity contribution is typically the smallest of these three and is conceptually the simplest, and here we utilize the AIMD simulation of the solution of interest to compute this term.

Anion Hydration.
QCT applies to both cation and anion hydration cases. In contrast to cations, however, anion hydration clusters often exhibit H-bond donation to the ion ( Figure 2; see also ref 70.). Anion-hydration clusters can be structurally delicate, specifically involving ligand−ligand hydrogen bonds, and that can make hydrated anions more challenging cases.
Initial QCT applications to hydrated anions worked simply with reasonable accuracy 3,13,14,46 compared with experiments. 71 Nevertheless, specifics of the technical ingredients can be perplexing ( Figure 3). The refinement of those initial applications has led to the further considerations discussed here, specifically, treatment of anharmonic effects on free energies of ion hydration clusters and the status of the polarizable continuum model 67 (PCM) for the hydration free energy of those clusters.

Theory Implemented with Simulation Data
We implement QCT here by bringing together quantities available from several different standard computations. The simulation work here thus contrasts with "let's take a look" direct numerical simulations. Indeed, QCT seeks to define minimal clusters that provide the information necessary for the statistical thermodynamic theory; thus, we explicitly do not attempt to construct observational dynamical simulations, nor do we seek a large cluster-size limit in our simulation calculations. 72 Before returning to discuss the theory, we note the details required for the computational procedures in the following. In addition, AIMD simulations of these systems are indeed readily available, and those observational calculations help to secure details that fill-out our understanding of these systems. Here, we note some of that previous work.
The extended work of Heuft and Meijer 73−75 initiated AIMD calculations on halide anions in water. They noted that residence times of water ligands in a halide inner shell spanned approximately 8, 12, and 17 ps for I − , Cl − , and F − , respectively. Those time scales are readily accessible by current AIMD calculations.
The interesting work of Wiktor et al. 76 focused on a basic thermodynamic quantity, the partial molar volumes of ions in water. Such studies are likely to provide fruitful next steps in the understanding of these systems.
The AIMD of Duignan and co-workers 77,78 on F − (aq) also focused on a basic thermodynamic characteristic of simple ions in water, namely, hydration free energies. They made a case for the application of ultrahigh accuracy electronic structure calculations to these problems. Our discussion below will identify aspects of the present efforts that overlap with that previous work but support a different conclusion: specifically, standard electronic structure calculations, properly integrated into statistical thermodynamic theory, are sufficient for experimental accuracy.
Finally, for this section, we note the extensive AIMD work on HO − (aq) that has been exhaustively reviewed. 16,17 Though that work did not proceed to the evaluation of standard thermodynamic characteristics, the discussions below will elaborate on specific points of comparison.
1.3.1. Procedures for Bulk X(aq) Solutions. The data utilized here for F − (aq) and Cl − (aq) was obtained from previous work 3, 14 that treated a single ion and 64 water molecules using the VASP simulation package. 79 The system was a cubic cell of edge 1.24 nm with periodic boundary conditions. The PW91 generalized gradient approximation described the core−valence interactions using the projector augmented wave (PAW) method. Plane waves with a kinetic energy cutoff of 400 eV and a time step of 0.5 fs were used for the simulation in the NVE ensemble. A temperature of 350 K was targeted for the simulation to avoid glassy behavior that can result at lower T. 3,14 After discarding 50 ps of trajectory as aging, our analysis was based on a 50 ps production trajectory.
For HO − (aq), Br − (aq), and I − (aq), the AIMD calculations are new here and used the CP2K simulation package 80 to treat a single ion and 64 water molecules under periodic boundary conditions. We adopted the PBE functional with Goedecker, Teter, and Hutter 81 (GTH) pseudopotentials in the GPW schemes, 82 as broadly used and consistent with our previous cluster results. 70 Molecularly optimized DZVPMOLOPT-SR-GTH basis sets were obtained from the CP2K website. Plane waves with a kinetic energy cutoff of 400 eV and a time step of 0.5 fs were used for the simulation in the NVT ensemble. The cubic cell with edge 1.27 nm reasonably matches the experimental density of water under our standard conditions. T = 300 K was selected 83 through the Nose−Hoover thermostat. Our analysis was based on 50 ps of production trajectory after 50 ps of aging. Figure 4 provides a standard overview of the bulk solution structures observed.

Procedures for Isolated (H 2 O) n X.
Molecular dynamics trajectories of the isolated (H 2 O) n X clusters for 2 ≤ n ≤ 5 and X = F − , HO − , Cl − , Br − , and I − were obtained using CP2K, just as for the HO − , Br − , and I − bulk solution calculations specified above. The pseudopotentials, functionals, and basis sets were the same. These simulations set the temperature at 300 K with the Nose−Hoover thermostat, using the GPW basis with default settings and a kinetic energy cutoff of 400 eV. Five picoseconds of the production trajectory, with a time step of 1 fs, was analyzed after 5 ps of aging.
For our statistical, or rough landscape, analysis of those trajectories, cluster structures were screened for consistency with our clustering definition ( Figure 5). Each clustered configuration is analyzed, according to the right-hand side of eq 5, and the separate structures are subjected to single-point calculations using the Gaussian 84 electronic structure software with the PBE functional and the DEF2TZVP basis set for F − , HO − , Cl − , Br − , and I − configurations. Using the resulting thermal averaging in eq 5 and K 1 (0) from experiment, the resulting stepwise K n (0) (Figure 6) produce the accurate results used below. 2,70

Inner-Shell Contributions
The study of the associative equilibria   is the traditional equilibrium constant, introduced with eq 3 below.

Accounts of Chemical Research
pubs.acs.org/accounts Article judgment might be required for a natural proximity definition, here we defer the discussion of that definition until after subsequent QCT developments. Our scheme for evaluating K n (0) is anchored in classic statistical thermodynamics 3,7,56 and proceeds incrementally following This formulation introduces the energy differences The brackets in eq 5, ⟨···|n⟩, indicate the thermal average utilizing the canonical simulation stream for the (H 2 O) n X cluster. 3 The symbol β stands for 1/RT. The evaluation of the energy combination ΔU n , eq 6, starts with the sampled configuration of the (H 2 O) n X cluster. Each ligand in turn serves to compose the energy difference suggested by the exchange Geometries of species on the right of eq 7 conform to the sampled (H 2 O) n X structure on the left. To appreciate ΔU n , we use the following accounting. Consider first the contribution  Figure 7).
Charge is balanced in eq 3, and the energy combination of eq 6 is not affected by the electrostatic potential of the phase. 4 For n = 1, eq 5 correctly reduces to the trivial case of K 0 (0) = 1. In evaluating K n (0) for n ≥ 2, the value of K 1 (0) can be supplied from experiment 71 or alternative theory. This term incorporates the interaction strength between X and one H 2 O molecule. Carrying out subsequent steps in this scheme then addresses the issues that make anion hydration more challenging, i.e., competing interactions of neighboring H 2 O molecules in those clusters. Figure 3 shows that a modest vertical shift of those harmonic approximation values substantially improves the agreement between harmonic approximation computations and experimental results throughout. Since the approach from eq 5 takes K 1 (0) as input, the present more detailed development should benefit similarly. This approach does not directly address issues of quantum mechanical zero-point motion, except to the extent that a pragmatic inclusion of an external K 1 (0) incorporates zeropoint motion empirically. The empirical grounding of this procedure might be particularly relevant to the multimodal possibilities of the H atom that is interior to (HO)HOH − , 16   similar to the Zundel cation (H 2 O)H(OH 2 ) + . This inner-shell evaluation can use electronic-structure methods of arbitrary sophistication since the calculations need be carried out only for modest-sized clusters. These results reexamine a previous application of QCT to HO − (aq) that was strikingly accurate for the dissociation thermodynamics. 46 2.1.1. Specific Definition of Clustering. In the discussion above, we did not address the consequences of any specific clustering proximity definition. That delinquency amounts to the assumption that complexes encountered in our cluster simulations, i.e., at moderate temperatures, are typically well clustered. Going further, we consider how a specific treatment might be built from here. A simple clustering criterion is that a clustered water molecule should donate at least one H atom within an assigned radius of the ion ( Figure 5).
We introduce a geometric indicator of a clustered configuration. b X (k) indicates whether molecule k is clustered (value 1) with the ion X or not (value 0). χ n takes a value of 1 if the n molecules are clustered with the ion and zero otherwise. Thus, χ n is indeed an indicator function. We will denote by K̅ n (0) as the equilibrium constant obtained by our scheme above without the definition of a restrictive innershell region. Then where ⟨χ n ⟩ indicates the statistical average evaluated as above, without a specific definition of clustered.
Note that 0 ≤ χ n ≤ 1 and thus ln[K n (0) /K̅ n (0) ] ≤ 0. This result makes physical sense because matching the natural clustering should not increase the free energy.
An interesting physical consideration is that the experimental results implicitly describe some specific physical clustering. Poorly formed physical clusters, not conforming to our mathematically defined cluster, might lose a weakly bound ligand, which would then populate the outer shell, perhaps to be pumped away in the experiment. These considerations raise the question, what clustering definition is most appropriate for the experiments?
The estimated average, ln⟨χ n ⟩, sampled from AIMD dynamics of (H 2 O) n X for 2 ≤ n ≤ 5 clusters (Table 1) shows that, generally, as n increases, ⟨χ n ⟩ decreases as repulsive interactions between water molecules force ligands to the outer shell. This behavior is especially evident for (H 2 O) 5 I − on the one hand, where none of the sampled configurations had all five water molecules within the defined inner shell. On the other hand, ln⟨χ n ⟩ ≈ 0 for (H 2 O) 3 HO − (Table 1), which indicates that the present clustering definition effectively encompasses the region of physical clustering for that case.

Outer Shell and PCM
The outer-shell cluster contribution of eq 2 to the hydration free energy can be treated 3 using the polarizable continuum model 67 (PCM) in the Gaussian suite of electronic structure programs. 84 The geometrical structures sampled from the AIMD trajectory for (H 2 O) n X were subjected to two single-point electronic calculations separately: one for the isolated cluster and a second with the external (dielectric) medium described by the PCM tool. The difference, ε ̅ j , is employed in computing Ä where m is the number of configurations from the simulation stream that satisfies the clustering definition. Equation 10 corresponds to the potential distribution theorem (PDT) approach, 1,56 recognizing that thermal fluctuations implicit in the PCM 67 approach complete the PDT averaging.
PCM is an approximate description of molecular-scale aspects of hydration. 85 PCM does include approximate accounts of packing effects and dispersion interactions, secondary to longrange electrostatic interactions. 67 Our discussion has emphasized that QCT is built by the identification of inner-shell clusters, the separate treatment of those clusters, and the integration of those results into the broader-scale solution environment. Inevitably, approximations are required to describe the outer-shell effects, and here those approximations are bundled in the PCM model.
The approximate character of the PCM model is signaled by the sensitive dependence on radii parameters that locate jumps in dielectric responses used in defining the model. It is reassuring that empirical values of those parameters are of reasonable magnitude. Still, the results are sensitive to those radii parameters, and they are not determined by theory 86 or experiment separate from the model. It is striking and important that QCT moderates that sensitivity by the subtraction of the ligand free energies in the formulation of that outer-shell contribution. 1,58 That insensitivity is achieved operationally by the boundaries of X being somewhat buried by the ligands and the ligand boundaries being unchanged in the subtraction (eq 2). The indicated subtraction reflects the appearance of the n H O 2 factors together with K n (0) of eq 2. For the cases of F − (aq) and Cl − (aq), we know already that this PCM-assisted application of QCT works satisfactorily 3 and Table 1. Estimated ln⟨χ n ⟩ from AIMD Calculations of (H 2 O) n X Defined Above a  46,70,87 For the cases of Br − (aq) and I − (aq), that procedure is somewhat degraded, 70 which we surmise as being due to the asymmetric hydration ( Figure 5) of those ions that leaves the central ion more exposed comparatively. Therefore, we tried the alternative approach, Ä (11) where the M clustered structures are sampled from the trajectory of the X(aq) AIMD simulations. This strategy is the well-known inverse formula of the standard potential distribution theorem, 56 but is implemented with the PCM tool. Our physical argument is that clustered (H 2 O) n X structures obtained with this inverse procedure ought to relieve intracluster hydrogen bonding and bury the X ion from contact with the water comparatively better (see also refs 88 and 69); therefore, the severe PCM approximation might perform better. Indeed, that was found to be the case, and the results discussed below (Figure 8) for the cases of Br − (aq) and I − (aq) were obtained with this inverse procedure.
The treatment of long-range electrostatic interactions in this implementation of QCT is worth emphasizing. Though the sampling of the bulk solution structures uses the preferred and commonplace periodic boundary conditions in treating longrang interactions, i.e., Ewald electrostatics, the energetics that enter into the free-energy computations reported do not use Ewald electrostatics.

Polydispersity and Net Hydration Free Energies
We have noted above that the polydispersity contribution is the smallest of the three contributions to cluster QCT. It is conceptually simplest and utilizes direct AIMD simulation in order to estimate p X (n). Thus, at this stage, we display all three contributions and their combination ( Table 2, then Figures 8  and 9) and then proceed to their physical discussion.
Note again that these free energies span a chemical scale of energies. For the least strongly bound case (I − ), the magnitudes of the net free energies are in excess of 60 kcal/mol, roughly 100RT here. The net quantities (eq 2, then Figure 8 and 9) are independent of n on that chemical energy scale. The agreement with the experimental tabulation (Table 2) is excellent and consistent across the anions treated. The latter point shows that  70 (Blue) Polydispersity contribution obtained from the Gaussian model for p(n) from AIMD data. In the cases of Br − (aq) and I − (aq), these results used the PCM in the inverse procedure of eq 11. Otherwise, results were obtained with direct eq 10. − nμ Hd 2 O (ex) )] which rearranges eq 2 and then acknowledges the normalization of p X (n). The root meansquared difference of these two columns is approximately 1.6 kcal/ mol.

Accounts of Chemical Research
pubs.acs.org/accounts Article the agreement is not affected by an assignment of a free-energy value for H + (aq) in the experimental tabulation.

Some Structural Observations
Though this work has explicitly marshalled simulation calculations toward evaluating QCT free energies, some structural observations are also available. The cluster definition ( Figure 5 and Table 1) based on Hatom donation, either singly or doubly, works satisfactorily: small values estimated for ln⟨χ n ⟩ indicate that the defined clustering volumes encompass the clustering observed from the AIMD for (H 2 O) n X cases, especially for the smaller values of n. For example, we see ln⟨χ 3 ⟩ ≈ 0 with (H 2 O) 3 HO − (Table 1). This result suggests simple H-bond donation, and that is supported by the radial distribution functions in the (H 2 O) 3 HO − case ( Figure 4). This observation provides a simple rationale for the remarkable success of QCT free-energy calculations for HO − (aq). 15,46,87 Thus, QCT for hydrated anions works better if the anion inner shell is characterized by Hbond donation. This conclusion is consistent with previous work 3,13,14 but focuses attention specifically on the XH radial distribution functions ( Figure 4).
Nevertheless, the structures of the clusters are different from the bulk aqueous solution case, as depicted by the XH rdfs ( Figure 4). For the clusters, setting aside F − and HO − , asymmetric hydration structures prevail ( Figure 5). Asymmetric inner-shell structures are moderated in the bulk hydration environment but still evident. 93 For the bulk hydration cases overall, the XH rdfs suggest the filling of inner shells to slightly higher average coordination numbers than for the clusters. Using HO − again as an example, the expected coordination number is about 3.8 (Figure 4), which may be compared to the work of ref 15 that estimated 3.7. The predictions of hydration structure agree well with experimental estimates based on X-ray and neutron diffraction and X-ray absorption fine structure studies (Table 3). Additional information, including traditional XO radial distribution functions, and further discussion can be found in the Ph.D. thesis of Diego T. Gomez. 70

CONCLUSIONS
The final free energies (Figures 8 and 9 and Table 2 Structural information obtained by simulations shows that the distinctive asymmetry of anion clusters is moderated in bulk aqueous solution.
The inner-shell free-energy contribution directly tracks available experimental information on gas-phase cluster hydration equilibria, and the polydispersity contribution is a direct structural observation from the AIMD trajectory. Neither of those contributions is expected to be sensitive to the potential of the phase.
The excellent theory−experiment agreement observed for those inner-shell cluster contributions is a breakthrough that supports the approximate remainder of the theory. The ndependent balance of the PCM-approximated, outer-shell contribution with the remaining numerically exact contributions suggests that the PCM approximation performs satisfactorily in these applications.
The implementation of this PCM-approximate outer-shell contribution requires statistical thermodynamic processing for the AIMD results, involving single-point electronic structure calculations of cluster structures sampled from the AIMD trajectory. Since the required electronic structure calculations treat only inner-shell clusters, this electronic structure effort could employ arbitrarily accurate numerical theory. In contrast to earlier recommendations, further effort in that direction is not warranted here because of the observed excellent theory− experiment agreement for the inner-shell cluster contribution.
In summary, the excellent agreement of anion hydration free energies is due to quantum computations accurately checked with experiment, together with the physical statistical thermodynamic theory that enables these computations. Furthermore, the AIMD simulations reveal differences in anion cluster structures compared with structures found in bulk aqueous solution, with the latter having less asymmetry and higher average coordination numbers. Finally, the ability to predict both an accurate solvation free energy and an accurate solvation structure of anions supports future work using QCT to understand the mechanisms of ion transport and selectivity for the large diversity of anion-selective transport proteins.   The authors declare no competing financial interest.