Vibrational spectroscopy of (SO4(2-)).(H2O)n clusters, n=1-5: harmonic and anharmonic calculations and experiment

The vibrational spectroscopy of (SO4(2-)).(H2O)n is studied by theoretical calculations for n=1-5, and the results are compared with experiments for n=3-5. The calculations use both ab initio MP2 and DFT/B3LYP potential energy surfaces. Both harmonic and anharmonic calculations are reported, the latter with the CC-VSCF method. The main findings are the following: (1) With one exception (H2O bending mode), the anharmonicity of the observed transitions, all in the experimental window of 540-1850 cm(-1), is negligible. The computed anharmonic coupling suggests that intramolecular vibrational redistribution does not play any role for the observed linewidths. (2) Comparison with experiment at the harmonic level of computed fundamental frequencies indicates that MP2 is significantly more accurate than DFT/B3LYP for these systems. (3) Strong anharmonic effects are, however, calculated for numerous transitions of these systems, which are outside the present observation window. These include fundamentals as well as combination modes. (4) Combination modes for the n=1 and n=2 clusters are computed. Several relatively strong combination transitions are predicted. These show strong anharmonic effects. (5) An interesting effect of the zero point energy (ZPE) on structure is found for (SO4(2-)).(H2O)(5): The global minimum of the potential energy corresponds to a C(s) structure, but with incorporation of ZPE the lowest energy structure is C2v, in accordance with experiment. (6) No stable structures were found for (OH-).(HSO4-).(H2O)n, for n相似文献   

Zero point energy of polyhedral water clusters

Polyhedral water clusters (PWCs) are cage-like (H2O)n clusters where every O participates in exactly three H bonds. For a database of 83 PWCs, 8 < or = n < or = 20, geometry was optimized and zero point energy (ZPE) was calculated at the B3LYP/6-311++G** level. ZPE correlates negatively with electronic energy (E0): each increase of 1 kcal/mol in E0 corresponds to a decrease of about 0.11 kcal/mol in ZPE. For each n, a set of four connectivity parameters accounts for 98% or more of the variance in ZPE. Linear regression of ZPE against n and this set gives an RMS error of 0.13 kcal/mol. The contributions to ZPE from stretch modes only (ZPE(S)) and from torsional modes only (ZPE(T)) also correlate strongly with E0 and with each other.     

Ab initio theoretical calculations of the electronic excitation energies of small water clusters

A direct ab initio molecular dynamics method has been applied to a water monomer and water clusters (H(2)O)(n) (n = 1-3) to elucidate the effects of zero-point energy (ZPE) vibration on the absorption spectra of water clusters. Static ab initio calculations without ZPE showed that the first electronic transitions of (H(2)O)(n), (1)B(1)←(1)A(1), are blue-shifted as a function of cluster size (n): 7.38 eV (n = 1), 7.58 eV (n = 2) and 8.01 eV (n = 3). The inclusion of the ZPE vibration strongly affects the excitation energies of a water dimer, and a long red-tail appears in the range of 6.42-6.90 eV due to the structural flexibility of a water dimer. The ultraviolet photodissociation of water clusters and water ice surfaces is relevant to these results.     

Structures, energetics, vibrational spectra of NH4+ (H2O)(n=4,6) clusters: Ab initio calculations and first principles molecular dynamics simulations

Important structural isomers of NH(4) (+)(H(2)O)(n=4,6) have been studied by using density functional theory, Moller-Plesset second order perturbation theory, and coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The zero-point energy (ZPE) correction to the complete basis set limit of the CCSD(T) binding energies and free energies is necessary to identify the low energy structures for NH(4) (+)(H(2)O)(n=4,6) because otherwise wrong structures could be assigned for the most probable structures. For NH(4) (+)(H(2)O)(6), the cage-type structure, which is more stable than the previously reported open structure before the ZPE correction, turns out to be less stable after the ZPE correction. In first principles Car-Parrinello molecular dynamics simulations around 100 K, the combined power spectrum of three lowest energy isomers of NH(4) (+)(H(2)O)(4) and two lowest energy isomers of NH(4) (+)(H(2)O)(6) explains each experimental IR spectrum.     

Theoretical study of pyridine and 4,4'-bipyridine adsorption on the lewis acid sites of alumina surfaces based on ab initio and density functional cluster calculations

The interactions of pyridine and 4,4'-bipyridine with the Lewis acid sites of alumina surfaces are investigated using ab initio and density functional calculations. Four cluster models of different sizes and shapes are chosen to represent the Lewis acid sites: three hydrogenated clusters Al(OH)(3), Al(4)O(9)H(6), and Al(10)O(21)H(12) and one non-hydrogenated cluster Al(4)O(6). The Hartree-Fock (HF) and B3LYP approaches with two basis sets 6-31G and 6-31+G are used to calculate the geometries, the electronic structures, the vibrational frequencies, and the adsorption energies of the complexes formed upon interaction of pyridine or 4,4'-bipyridine ligands on the cluster surfaces. Electronic structures are determined by the electrostatic potential (ESP) analysis of charges. Adsorption energies are calculated with corrections made for zero-point energies (ZPE) and basis set superposition error (BSSE). The ESP analysis of atomic charges reveals that the charge-transfer effects are more important in Lewis complexes formed with Al(4)O(6) cluster than in those formed with hydrogenated clusters Al(OH)(3), Al(4)O(9)H(6), and Al(10)O(21)H(12). The significantly larger charge transferred from pyridine or 4,4'-bipyridine ligand to Al(4)O(6) cluster should increase the adsorption energy of these complexes. Consequently, at all levels of calculation, the adsorption energies of pyridine and 4,4'-bipyridine complexed to Al(4)O(6) cluster ( approximately 46 kcal/mol), which compare very well to experiment, are strongly larger than those obtained for both pyridine and 4,4'-bipyridine ligands complexed to Al(OH)(3) (32 kcal/mol), Al(4)O(9)H(6) (24 kcal/mol) and Al(10)O(21)H(12) (25 kcal/mol) clusters. The corrected adsorption energy is found to be insensitive to basis set and electron correlation effects. It essentially depends on the ionic character of the cluster model rather than on its size. For 4,4'-bipyridine complexes, similar results to those obtained for pyridine are found, and the geometry and the amount of charge of the unbound pyridyl ring are unchanged upon complexation. The calculated vibrational frequencies and frequency shifts are little sensitive to the size and shape of the cluster model. The two ring stretching modes 8a and 19b of pyridine and 4,4'-bipyridine observed in the 1400-1600 cm(-1) region are the most affected modes upon adsorption, in good agreement with the available infrared and Raman data.     

Stepwise hydration of ionized aromatics. Energies, structures of the hydrated benzene cation, and the mechanism of deprotonation reactions

The stepwise binding energies (DeltaHdegree(n-1,n)) of 1-8 water molecules to benzene(.+) [Bz(.+)(H2O)n] were determined by equilibrium measurements using an ion mobility cell. The stepwise hydration energies, DeltaHdegree(n-1,n), are nearly constant at 8.5 +/- 1 kcal mol-1 from n = 1-6. Calculations show that in the n = 1-4 clusters, the benzene(.+) ion retains over 90% of the charge, and it is extremely solvated, that is, hydrogen bonded to an (H2O)n cluster. The binding energies and entropies are larger in the n = 7 and 8 clusters, suggesting cyclic or cage-like water structures. The concentration of the n = 3 cluster is always small, suggesting that deprotonation depletes this ion, consistent with the thermochemistry since associative deprotonation Bz(.+)(H2O)(n-1) + H2O-->C6H5. + (H2O)nH+ is thermoneutral or exothermic for n > or = 4. Associative intracluster proton transfer Bz(.+)(H2O)(n+1) + H2O-->C6H5.(H2O)nH+ would also be exothermic for n > or = 4, but lack of H/D exchange with D2O shows that the proton remains on C6H6(.+) in the observed Bz(.+)(H2O)n clusters. This suggests a barrier to intracluster proton transfer, and as a result, the [Bz(.+)(H2O)n]* activated complexes either undergo dissociative proton transfer, resulting in deprotonation and generation of (H2O)nH+, or become stabilized. The rate constant for the deprotonation reaction shows a uniquely large negative temperature coefficient of K = cT(-67+/-4) (or activation energy of -34+/- 1 kcal mol-1), caused by a multibody mechanism in which five or more components need to be assembled for the reaction.     

Exploring the rich energy landscape of sulfate-water clusters SO4(2-) (H2O)(n=3-7): an electronic structure approach

We present a reinvestigation of sulfate-water clusters SO4(2-) (H2O)(n=3-7), which involves several new aspects. Using a joint molecular mechanics/first principles approach, we perform exhaustive searches for stable cluster geometries, showing that the sulfate-water landscape is much richer than anticipated previously. We check the compatibility of the new structures with experiment by comparing vertical detachment energies (VDEs) calculated at the B3LYP/6-311++G** level of theory and determine the energetic ordering of the isomers at the RI-MP2/aug-cc-pVTZ level. Our results are bench-marked carefully against reference energies of estimated CCSD(T)/aug-cc-VTZ quality and VDEs of CCSD(T)/aug-cc-pVDZ quality. Furthermore, we calculate anharmonic vibrational corrections for up to the n = 6 clusters, which are shown to be significant for isomer energy ordering. We use energy decomposition analysis (EDA) based on the absolutely localized fragment (ALMO) expansion to gain chemical insight into the binding motifs.     

Study of the structure, energetics, and vibrational properties of small ammonia clusters (NH3)n (n = 2-5) using correlated ab initio methods

Equilibrium geometries, interaction energies, and harmonic frequencies of (NH3)n isomers (n = 2-5) have been computed using correlated calculations (MP2) in conjunction with Dunning's aug-cc-pVXZ (X = D, T, Q) basis sets and the Counterpoise procedure. Whenever available, literature values for the binding energy and geometry of dimers and trimers agree well with our data. Low lying isomers for (NH3)4 and (NH3)5 have been found to have similar binding energies (roughly 16 and 20 kcal/mol for the tetramer and pentamer, respectively), perhaps suggesting the presence of a very smooth energy landscape. Using BSSE corrected forces or freezing the monomer structure to its gas phase geometry have been found to have only a weak impact on the energetic and structural properties of the clusters. The effect of zero-point energy (ZPE) on the relative stability of the clusters has been estimated using harmonic frequencies. The latter also highlighted the presence of vibrational fingerprints for the presence of double acceptor ammonia molecules. Many-body effects for (NH3)n isomers (n = 2-4) have been investigated to explore the possibility of building a pairwise interaction model for ammonia. In the frame of the work presented, we have found the 3-body effect to account for 10-15% of the total interaction energy, whereas the 4-body effects may be neglected as first approximation.     

Structure and spectral features of H+(H2O)7: Eigen versus Zundel forms

The two dimensional (2D) to three dimensional (3D) transition for the protonated water cluster has been controversial, in particular, for H(+)(H(2)O)(7). For H(+)(H(2)O)(7) the 3D structure is predicted to be lower in energy than the 2D structure at most levels of theory without zero-point energy (ZPE) correction. On the other hand, with ZPE correction it is predicted to be either 2D or 3D depending on the calculational levels. Although the ZPE correction favors the 3D structure at the level of coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] using the aug-cc-pVDZ basis set, the result based on the anharmonic zero-point vibrational energy correction favors the 2D structure. Therefore, the authors investigated the energies based on the complete basis set limit scheme (which we devised in an unbiased way) at the resolution of the identity approximation Moller-Plesset second order perturbation theory and CCSD(T) levels, and found that the 2D structure has the lowest energy for H(+)(H(2)O)(7) [though nearly isoenergetic to the 3D structure for D(+)(D(2)O)(7)]. This structure has the Zundel-type configuration, but it shows the quantum probabilistic distribution including some of the Eigen-type configuration. The vibrational spectra of MP2/aug-cc-pVDZ calculations and Car-Parrinello molecular dynamics simulations, taking into account the thermal and dynamic effects, show that the 2D Zundel-type form is in good agreement with experiments.     

Exploring the correlation between network structure and electron binding energy in the (H(2)O)(7)(-) cluster through isomer-photoselected vibrational predissociation spectroscopy and ab initio calculations: addressing complexity beyond types I-III

We report a combined photoelectron and vibrational spectroscopy study of the (H(2)O)(7)(-) cluster anions in order to correlate structural changes with the observed differences in electron binding energies of the various isomers. Photoelectron spectra of the (H(2)O)(7)(-) . Ar(m) clusters are obtained over the range of m=0-10. These spectra reveal the formation of a new isomer (I') for m>5, the electron binding energy of which is about 0.15 eV higher than that of the type I form previously reported to be the highest binding energy species [Coe et al., J. Chem. Phys. 92, 3980 (1990)]. Isomer-selective vibrational predissociation spectra are obtained using both the Ar dependence of the isomer distribution and photochemical depopulation of the more weakly (electron) binding isomers. The likely structures of the isomers at play are identified with the aid of electronic structure calculations, and the electron binding energies, as well as harmonic vibrational spectra, are calculated for 28 low-lying forms for comparison with the experimental results. The HOH bending spectrum of the low binding type II form is dominated by a band that is moderately redshifted relative to the bending origin of the bare water molecule. Calculations trace this feature primarily to the bending vibration localized on a water molecule in which a dangling H atom points toward the electron cloud. Both higher binding forms (I and I') display the characteristic patterns in the bending and OH stretching regions signaling electron attachment primarily to a water molecule in an AA binding site, a persistent motif found in non-isomer-selective spectra of the clusters up to (H(2)O)(50)(-).     

Infrared spectroscopy of water cluster anions, (H2O)n=3-24-)in the HOH bending region: persistence of the double H-bond acceptor (AA) water molecule in the excess electron binding site of the class I isomers

We report vibrational predissociation spectra of water cluster anions, (H(2)O)(n=)()(3)(-)(24)(-) in the HOH bending region to explore whether the characteristic red-shifted feature associated with electron binding onto a double H-bond acceptor (AA) water molecule survives into the intermediate cluster size regime. The spectra of the "tagged" (H(2)O)(n)()(-).Ar clusters indeed exhibit the signature AA band, but assignment of this motif to a particular isomer is complicated by the fact that argon attachment produces significant population of three isomeric forms (as evidenced by their photoelectron spectra). We therefore also investigated the bare clusters since they can be prepared exclusively in the high binding (isomer class I) form. Because the energy required to dissociate a water molecule from the bare complexes is much larger than the transition energies in the bending region, the resulting (linear) action spectroscopy selectively explores the properties of clusters with most internal energy content. The (H(2)O)(15)(-) predissociation spectrum obtained under these conditions displays a more intense AA feature than was found in the spectra of the Ar tagged species. This observation implies that not only is the AA motif present in the class I isomer, but also that it persists when the clusters contain considerable internal energy.     

Isomer-selective detection of microsolvated oxonium and carbenium ions of protonated phenol: infrared spectra of C6H7O+-Ln clusters (L=Ar/N2, n</=6)

Infrared photodissociation (IRPD) spectra of clusters composed of protonated phenol (C(6)H(7)O(+)) and several ligands L are recorded in the O-H and C-H stretch ranges using a tandem mass spectrometer coupled to a cluster ion source. The C(6)H(7)O(+)-L(n) complexes (L=Ar/N(2), n=1-6) are generated by chemical ionization of a supersonic expansion. The IRPD spectra of mass selected C(6)H(7)O(+)-L(n) clusters obtained in various C(6)H(7)O(+)-L(m) fragment channels (m相似文献   

A comparative ab initio study of Br2*- and Br2 water clusters

The work presents ab initio results on structure and electronic properties of Br2*-.nH2O(n=1-10) and Br2.nH2O(n=1-8) hydrated clusters to study the effects of an excess electron on the microhydration of the halide dimer. A nonlocal density functional, namely, Becke's half-and-half hybrid exchange-correlation functional is found to perform well on the present systems with a split valence 6-31++G(d,p) basis function. Geometry optimizations for all the clusters are carried out with several initial guess structures and without imposing any symmetry restriction. Br2*-.nH2O clusters prefer to have symmetrical double hydrogen-bonding structures. Results on Br2.nH2O(n>or=2) cluster show that the O atom of one H2O is oriented towards one Br atom and the H atom of another H2O is directed to other Br atom making Br2 to exist as Br+-Br- entity in the cluster. The binding and solvation energies are calculated for the Br2*-.nH2O and Br2.nH2O clusters. Calculations of the vibrational frequencies show that the formation of Br2*- and Br2 water clusters induces significant shifts from the normal modes of isolated water. Excited-state calculations are carried out on Br2*-.nH2O clusters following configuration interaction with single electron excitation procedure and UV-VIS absorption profiles are simulated. There is an excellent agreement between the present theoretical UV-VIS spectra of Br2*-.10H2O cluster and the reported transient optical spectra for Br2*- in aqueous solution.     

Benchmark structures and binding energies of small water clusters with anharmonicity corrections

For (H(2)O)(n) where n = 1-10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.     

Theoretical study of the geometries and dissociation energies of molecular water on neutral aluminum clusters Al(n) (n = 2-25)

Geometries and dissociation energies of water molecules on Al(n) (n = 2-25) clusters were investigated using density functional theory with all electron relativistic spin-polarized calculations under the generalized gradient approximation. An extensive structure search was performed to identify the low-energy conformations of Al(n)H(2)O complexes for each size. Optimal adsorption sites were assigned for low-energy isomers of the clusters. Size and site specific dependences were studied for the Al(n)H(2)O complexes in stabilities, geometries, adsorption energies, dissociation energies, Al-O bond lengths, and other characteristic quantities. The stabilities and geometries revealed that H atom in H(2)O is not inclined to bond with Al atoms. The most stable Al(n)H(2)O configurations for each size tend to correspond to the most stable bare Al(n) cluster except of Al(6) and Al(24) clusters. The HO bond lengths increase generally 0.01 ? with respect to the isolated H(2)O in all of the adsorption complexes. The dissociation energy of an isolated H(2)O into HO and H was 5.39 eV, which decreased about two-thirds to the energy range of 0.83-2.12 eV with the help of Al(n) clusters. In spite of the fluctuations, the dissociation energies of Al(n)H(2)O complexes rise with the size increasing as a whole. In addition, we also found that the bare Al(n) clusters with high vertical ionization potentials usually have high dissociation energies of H(2)O in the corresponding adsorption models. The energetically preferred spin-multiplicity of all the odd-n Al(n)H(2)O complexes is doublet, and it is singlet for all the even-n complexes with exception of Al(2)H(2)O which is triplet.     

The structures and electronic states of zinc-water clusters Zn(n)(H2O)(m) (n = 1-32 and m = 1-3)

Ab initio and Density Functional Theory (DFT) calculations have been carried out for zinc-water clusters Zn(n)-(H2O)(m) (n = 1-32 and m = 1-3, where n and m are the numbers of zinc atoms and water molecules, respectively) to elucidate the structure and electronic states of the clusters and the interaction of zinc cluster with water molecules. The binding energies of H2O to zinc clusters were small at n = 2-3 (2.3-4.2 kcal mol(-1)), whereas the energy increased significantly in n = 4 (9.0 kcal mol(-1)). Also, the binding nature of H2O was changed at n = 4. The cluster size dependency of the binding energy of H2O accorded well with that of the natural population of electrons in the 4p orbital of the zinc atom. In the larger clusters (n > 20), it was found that the zinc atoms in surface regions of the zinc cluster have a positive charge, whereas those in the interior region have a negative charge with the large electron population in the 4p orbital. The interaction of H2O with the zinc clusters were discussed on the basis of the theoretical results.     

Fragmentation cross sections of protonated water clusters

We have measured fragmentation cross sections of protonated water cluster cations (H(2)O)(n=30-50)H(+) by collision with water molecules. The clusters have well-defined sizes and internal energies. The collision energy has been varied from 0.5 to 300 eV. We also performed the same measurements on deuterated water clusters (D(2)O)(n=5-45)D(+) colliding with deuterated water molecules. The main fragmentation channel is shown to be a sequential thermal evaporation of single molecules following an initial transfer of relative kinetic energy into internal energy of the cluster. Unexpectedly, that initial transfer is very low on average, of the order of 1% of collision energy. We evaluate that for direct collisions (i.e., within the hard sphere radius), the probability for observing no fragmentation at all is more than 35%, independently of cluster size and collision energy, over our range of study. Such an effect is well known at higher energies, where it is attributed to electronic effects, but has been reported only in a theoretical study of the collision of helium atoms with sodium clusters in that energy range, where only vibrational excitation occurs.     

Charge penetration and the origin of large O-H vibrational red-shifts in hydrated-electron clusters, (H2O)n-

The origin of O-H vibrational red-shifts observed experimentally in (H2O)n(-) clusters is analyzed using electronic structure calculations, including natural bond orbital analysis. The red-shifts are shown to arise from significant charge transfer and strong donor-acceptor stabilization between the unpaired electron and O-H sigma* orbitals on a nearby water molecule in a double hydrogen-bond-acceptor ("AA") configuration. The extent of e(-) --> sigma* charge transfer is comparable to the n --> sigma* charge transfer in the most strongly hydrogen-bonded X(-)(H2O) complexes (e.g., X = F, O, OH), even though the latter systems exhibit much larger vibrational red-shifts. In X(-)(H2O), the proton affinity of X(-) induces a low-energy XH...(-)OH diabatic state that becomes accessible in v = 1 of the shared-proton stretch, leading to substantial anharmonicity in this mode. In contrast, the H + (-)OH(H2O)(n-1) diabat of (H2O)n(-) is not energetically accessible; thus, the O-H stretching modes of the AA water are reasonably harmonic, and their red-shifts are less dramatic. Only a small amount of charge penetrates beyond the AA water molecule, even upon vibrational excitation of these AA modes. Implications for modeling of the aqueous electron are discussed.     

Proton and deuteron position preferences in water clusters: an ab initio study

In order to explore the effect of H-to-D substitution on the zero-point energy (ZPE) of water clusters, Hessians were computed for a database of 53 optimized (H2O)n clusters, 5 < or = n < or = 21, at the B3LYP6-311 + + G** level. The 53 clusters contained 1524 protons, which were sorted into 18 categories according to the type of their donor O and (if not free) acceptor O. Letting deltaZPE[H]* denote the change in ZPE when the proton H* is replaced by D, mean values for deltaZPE[H*] for the H-bonded categories ranged from -2172 cal mol(-1) for H* in a DDAA-DDAA bond to -2118 for H* in a DAA-DDA bond. Mean value for H* free on DAA (respectively, DA) was -2018 (respectively, -1969). For DAA-DDA bonds, and for short H bonds in general, there was a strong inverse correlation between /deltaZPE[H*]/ and the O-H* distance. deltaZPE for multiple H-to-D substitutions was additive, except for a cooperativity effect of -13.7 to -19.7 cal mol(-1) when two substituted protons were in the same H2O unit and a much smaller cooperativity when one proton's donor was the other's acceptor. Implications of these data include a relative preference for D to occupy H bonded rather than free positions in finite water clusters, a value of 3.82 for the disproportionation equilibrium constant of mixed ice at 150 K, increased occupation by H at surface positions of mixed ice, and a larger average coordination number for liquid D2O than for liquid H2O.     

Hydrogen bonding interactions of pyridine*+ with water: stepwise solvation of distonic cations

The solvation energies of the pyridine*+ radical cation by 1-4 H2O molecules were determined by equilibrium measurements in a drift cell. The binding energies of the pyridine*+(H2O)n clusters are similar to the binding energies of protonated pyridine-water clusters, (C5H5NH+)(H2O)n, which involve NH+..OH2 bonds and different from those of the solvated benzene radical cation-water clusters, C6H6*+(H2O)n, which involve CHdelta+..OH2 bonds. These relations indicate that the observed pyridine*+ ions have the distonic *C5H4NH+ structures that can form NH+..OH2 bonds. The observed thermochemistry and ab initio calculations show that these bonds are not affected significantly by an unpaired electron at another site of the ion. Similar observations also identify the 2-fluoropyridine*+ distonic ion. The distonic structure is also consistent with the reactivity of pyridine*+ in H atom transfer, intra-cluster proton transfer and deprotonation reactions. The results present the first measured stepwise solvation energies of distonic ions, and demonstrate that cluster thermochemistry can identify distonic structures.