Gas Phase Vibrational Spectroscopy of the Protonated Water Pentamer : The Role of Isomers and Nuclear Quantum Effects

,


Introduction
How ions are hydrated has intrigued chemists since the birth of the field of Physical Chemistry. 1 Proton hydration, in particular, is of fundamental importance in chemical, biological and atmospheric processes, but the molecular level interpretation of the anomalously high mobility of protons in aqueous solutions remains hotly debated. 2 Size-selected protonated water clusters, H + (H 2 O) n , isolated in the gas phase, play an important role in clarifying the spectroscopic markers that encode the collective mechanics involved in long range proton translocation.
The advantage of isolated clusters is that they allow direct spectroscopic characterization of key local hydration motifs that are metastable in solution, but can occur naturally in the size-and temperature-dependent structures of the gas-phase clusters.Since protonated water clusters can be characterized experimentally with high selectivity and sensitivity and are also amenable to higher level quantum chemical calculations, such measurements deliver stringent benchmarks for evaluating the validity of approximations necessary for treating ions in extended solvation environments.However, the widely accepted structural assignment of small protonated water clusters, in particular that of the protonated water tetramer and pentamer, has recently been drawn into question based on the implications of ab initio molecular dynamics (AIMD) simulations. 3,4n this study, we characterize the gas phase vibrational spectroscopy of the protonated water pentamer over nearly the complete IR spectral range and show experimentally, using an isomer-specific, double resonance technique as well as temperature-dependent measurements (13-250 K), that the observed IR transitions are due to a single spectroscopically unique species.We then address the assignments of the strong bands in the cold spectrum to fundamentals expected for the branched structure originally proposed by Jiang et al. 5 with the aid of anharmonic spectra from electronic structure calculations.

View Article Online View Journal | View Issue
This journal is © the Owner Societies 2016 Finally, we use RRKM theory to explore the kinetic stability of the isomers at low temperatures, highlighting the importance of considering nuclear quantum effects.
The vibrational spectroscopy of the protonated water pentamer has been studied by several groups.Jiang et al. 5 measured the infrared photodissociation (IRPD) spectrum of H + (H 2 O) 5 in the O-H stretching region (2700-3900 cm À1 ) at an estimated temperature of 170 K. Aided by harmonic vibrational frequencies and intensities from electronic structure calculations, they identified an acyclic, branched isomer (labelled B in Fig. 1) as the signal carrier in the IRPD spectrum.Temperature-dependent measurements confirmed that the IRPD spectrum was free of contributions from other isomers.Subsequently, Headrick et al. 6 probed colder clusters using argon-tagging.They measured IRPD spectra down into the fingerprint region (1000-3900 cm À1 ), identifying the evolution of the spectral signature of the hydrated proton with cluster size.These studies confirmed the assignment to the branched isomer B. 6 The effect of messenger-tagging on the IRPD spectra of protonated water clusters was systematically studied by Mizuse et al. also confirming the presence of a single isomer (B). 7agging with H 2 , Ar and Ne leads to sharper vibrational features due to colder clusters but hardly perturbs the IRPD spectrum.Douberly et al. then studied the effect of deuteration on the band pattern in the IRPD spectrum in the O-D stretching region (2300-2900 cm À1 ). 8 They found that the spectra of H + (H 2 O) 5 ÁAr and D + (D 2 O) 5 ÁAr both agree best with the predicted harmonic spectrum of the branched isomer B and that the contribution from other isomers was insignificant.This is particularly remarkable, since their electronic structure calculations predict that B is not the global minimum energy structure. 8lectronic structure calculations 3,5,[8][9][10][11][12][13] (Monte Carlo, DFT and MP2) of bare H + (H 2 O) 5 consistently predicted a cyclic fourmember ring with an Eigen core (labelled R in Fig. 1) as the global minimum-energy structure, closely followed by two noncyclic Eigen isomers: one branched (B: +3.7 kJ mol À1 ) 13 and the other non-branched (C: +17 kJ mol À1 ). 13 When zero-point energy (ZPE) is considered, the energetic ordering changes and B is found to be 0.3 kJ mol À1 more stable than R. 8 Such a change in ordering is not found for the perdeuterated analogue because its ZPE is smaller.In both cases, however, B is favoured by entropy and it represents the global minimum above 50 K in terms of free energy. 8The comparison of the computed harmonic spectra with the experimental IR spectra, in particular in the O-H (O-D) stretching region, supports the assignment of the protonated water pentamer and of its fully deuterated isotopologue to the branched isomer B. [5][6][7][8] The agreement is less satisfactory in the fingerprint region, namely for the Eigen core O-H stretching mode, which involves a displacement towards the single acceptorsingle donor (1-AD) H 2 O.It is predicted by MP2 calculations to occur about 460 cm À1 above the experimentally observed band at 1885 cm À1 . 6However, Vibrational Self-Consistent Field calculations (VSCF) accurately recover the red shift of this band, thus supporting the assignment to isomer B. 6 The complexity of this assignment is apparent, however, by the fact that application of the VPT2 perturbative treatment supports the assignment of a previously unassigned, relatively strong feature at 1470 cm À1 (labelled a 9 in Fig. 3) to the key OH stretching band of the embedded H 3 O + moiety. 14The appearance of such unassigned ''extra features'' (i.e., those not readily assignable to fundamentals) 6,14 has, in fact raised the spectre that distinct isomers may be in play, and indeed the assignment of the vibrational spectrum to the branched isomer B has recently been challenged by Kulig and Agmon, based on insights derived from AIMD simulations. 3These authors suggest that two other isomers contribute to the observed spectrum: the cyclic Eigen isomer R and the higher energy chain isomer C. A major motivation for the present study was, therefore, to address the heterogeneity in the spectrum of the light isotopologue and to follow how the bands evolve with increasing temperature as well as upon isotopic substitution.

Experimental
IRPD experiments are conducted using a 6 K ion-trap triple mass spectrometer described previously. 15,16Clusters are produced in a nanospray ion-source using a 10 mM H 2 SO 4 in CH 3 OH/H 2 O (1 : 1, v : v) solution and D 2 SO 4 in a CD 3 OD/D 2 O (2 : 1, v : v) solution for the all-H and all-D isotopologues, respectively.The ions of interest are thermalized at room temperature in a gasfilled ion guide and then mass-selected using a quadrupole mass filter and focused in a radio-frequency (RF) ring-electrode iontrap.The trap, filled with pure D 2 buffer gas, is cooled with a closed-cycle He cryostat and held at temperatures between 13 K and 250 K.Many collisions of the trapped ions with the buffer gas provide gentle cooling of the internal degrees of freedom close to the ambient temperature and typically avoid the production of kinetically trapped species.At sufficiently low ion-trap temperatures, ion-messenger complexes are formed via threebody collisions. 17,18All ions are extracted either every 100 ms or 200 ms, depending on the lasers used, and focused both temporally and spatially into the centre of the extraction region of an orthogonally mounted reflectron time-of-flight (TOF) tandem photofragmentation mass spectrometer.Here, the ions are irradiated with a counter-propagating IR laser pulse produced either by one of the two OPO/OPA/AgGaSe 2 IR laser systems 19 (700-4500 cm À1 ) or by the FHI-FEL (210-3000 cm À1 ). 20All parent and photofragment ions are then accelerated towards an MCP detector.An IRPD spectrum is recorded by averaging over 50-200 TOF mass spectra per wavelength and scanning the wavelength of the laser pulse.The photodissociation cross section s IRPD is determined as described previously. 15,16he covered spectral range from 220-4000 cm À1 requires the use of three different laser configurations.The individual spectra covering these overlapping spectral regions (200-1000 cm À1 , 900-2300 cm À1 and 2000-4000 cm À1 ) are then scaled in intensity such that the intensity of common peaks match, in order to correct for changes in laser beam/ion cloud overlap.
Alternatively, isomer-specific IR 2 MS 2 double resonance measurements [21][22][23] can be performed using an analogous setup, employing two IR laser pulses. 24,25The first laser pulse generates a set of photofragment ions, which are separated from parent ions in the 1801 reflectron stage of the tandem TOF mass-spectrometer, where they separate in time and space according to their mass-to-charge ratio.They are then refocused between the extraction plates, and the parent ion packet is irradiated with a second, properly timed laser pulse to generate a second set of photofragment ions.All ions located between the acceleration plates are then reaccelerated and mass analysed using the linear TOF stage. 15,16e used two ways to measure IR 2 MS 2 spectra, which differ in the sequence of the probe and scan laser pulses.If the probe pulse (fixed wavelength) interacts with the ions after the scan pulse (tuneable wavelength), the IR 2 MS 2 spectrum manifests itself in the form of dips in the otherwise constant fragment ion signal from the probe pulse (ion-dip technique).If the scanned pulse is applied after the probe pulse, then two TOF spectra are acquired, one with and the other without the probe pulse using an additional mechanical chopper (a population-modulation or hole-burning technique).In this case, the IR 2 MS 2 spectrum is obtained by subtracting these two TOF spectra.Probe laser pulses with wavenumbers larger than 2300 cm À1 are typically generated with one of the OPO/OPA systems while those at smaller wavenumbers (o2300 cm À1 ) are produced by the FHI FEL.The macropulse duration of the IR FEL pulses (B10 ms) requires that the FEL laser beam is aligned counter-propagating to the ion cloud extracted from the trap and prior to acceleration into the tandem TOF mass spectrometer to ensure a sufficient temporal overlap.Hence, the FHI FEL laser pulse always needs to be timed prior to the laser pulse from the OPO/OPA system, which is sufficiently short in duration (B7 ns) that the fast (few keV) ions can be efficiently irradiated.

Quantum chemistry
Electronic structure calculations were performed using the Gaussian 09 rev.D.01 software package. 269][30][31] Additionally, we checked the influence of dispersion interactions by using the D3 version of Grimme's dispersion correction. 32Tight convergence criteria and the ''superfine'' grid option are employed.We identified transition states (TS) of the unimolecular isomerization reactions by locating first-order saddle points along the reaction coordinate.The reaction paths were followed using the intrinsic reaction coordinate (IRC) method. 33,34ibrational harmonic frequencies of minima and TS configurations were computed at the same level of theory.To account for anharmonic effects as well as systematic errors on the harmonic force constants, we determined scaling factors on the basis of a comparison of the IRPD spectra with harmonic spectra of isomer B for both the all-H and all-D isotopologues.Two scaling factors for each method (see Table S4, ESI †), one for the spectral region below 1500 cm À1 (1100 cm À1 for the perdeuterated complex) and one for the higher energy region, were calculated with a least-squares procedure by minimizing the residuals as proposed by Scott and Radom. 35Additionally, we performed anharmonic frequency calculations using the Second-order Vibrational Perturbation Theory (VPT2) method as implemented by Bloino and Barone. 36To obtain more reliable electronic energies and further improve the energetic description of the isomers, we performed CCSD(T)/aug-cc-pVTZ single point calculations at the DFT geometries.Optimized geometries, energies, vibrational spectra, IRC paths and scaling factors are reported in the ESI † (Fig. S1, S2, S4 and S5; Tables S1-S13, ESI †).

RRKM rate constants and kinetics
Microcanonical isomerization rate constants of the protonated water pentamer k(E), where E refers to the internal energy of the system, were calculated using RRKM theory. 37The sequential kinetic scheme consisting of three isomerization equilibria is depicted in Fig. 2. We used CCSD(T) energies of the B3LYP-D3 structures plus the B3LYP-D3 ZPEs, derived from scaled harmonic wavenumbers.To obtain the rate constants, we used the code written by Zhu and Hase, 38 using the standard formula: 37 where s is the reaction path degeneracy, E 0 is the unimolecular threshold, N # (E À E 0 ) is the TS sum of states and r(E) the reactant density of states.N # (E À E 0 ) and r(E) refer only to active degrees of freedom.The reaction path is doubly degenerate for all isomerizations considered.Note, the contribution from rotational excitation was neglected, since the minima and TS configurations have very similar moments of inertia and thus cancelled out.Indeed, test calculations with different J values provided consistent results.Classical RRKM rate constants were calculated disregarding quantization of vibrational modes with classical state counting.In contrast, to obtain the quantum rate constants, ZPEs were taken into account and state counting was performed with the direct count method using the Beyer-Swinehart algorithm. 39rom microcanonical rate constants and densities of states, the canonical rate constants can be expressed as 40 kðTÞ ¼ 1 QðTÞ L½kðEÞrðEÞ where L is the Laplace transform operator and Q(T) is the partition function, which can be expressed as the Laplace transform of density of states, QðTÞ ¼ L½rðEÞ.
This journal is © the Owner Societies 2016 We then solved the four kinetic differential equations associated with the system for different initial conditions to study the kinetic stability of the four isomers: The kinetic constants k 1 , k À1 , k 2 , k À2 , k 3 and k À3 of the isomerization reactions considered are defined in Fig. 2. The system of differential equations was numerically solved using the ODE45 solver as implemented in Matlab 41,42 and yields the time evolution of the relative isomer populations.

Results and discussion
Fig. 3 shows a comparison of the IRPD spectra of the D 2 -tagged protonated water pentamer and its fully perdeuterated isotopologue from 210 cm À1 to 4000 cm À1 .Band positions are listed in Table 1.For better comparability, a scaling factor (1/1.36) is used for the alignment of the two wavenumber scales, which brings the sharp water bending fundamentals of the two isotopologues (a 8 and b 8 ) into alignment.The scaling factor is close to the theoretical value of 1/1.37 that is expected for isolated O-H vs. O-D oscillators.][7][8] Additionally, the D 2 -tagged IRPD spectrum of D + (D 2 O) 5 (1000-3000 cm À1 ) measured at the Yale photofragmentation mass spectrometer 43 confirms the reproducibility of the experimental results presented here (see Fig. S3, ESI †).The striking resemblance of the two spectra shown in Fig. 3 provides unambiguous evidence that both spectra probe signal carriers with similar structures and, in particular, identical connectivity.Note that the absorption near the key a 9 feature is substantially less intense in the spectrum of the perdeuterated species, as are the weaker features that appear higher in energy, labelled by * and †.Such selectivity in the affected bands could arise from the appearance of a second isomer in the all-H isotopologue, or from the suppression of anharmonic bands in the heavy isotopologue.As mentioned above, Kulig and Agmon 3 have recently invoked an assignment scheme for the H + (H 2 O) 5 based on two isomers (R and C in Fig. 2) and not involving isomer B. We address these key issues of isomeric composition and temperature dependence in the next two sections.

Double resonance IR 2 MS 2 measurements: identifying the number of spectroscopically unique species
2][23] More precisely, the method is sensitive to laser-induced population changes that are maintained longer than the delay between the two IR laser pulses (B10 À5 s).Hence, the IR spectra of more quickly interconverting isomers cannot be isolated using this double resonance technique, but their presence may be indicated by broadening of spectral features (see below).
The results of the IR 2 MS 2 experiments are shown in Fig. 4, along with the IRPD spectrum of H + (H 2 O) 5 ÁD 2 (top trace).The IR 2 MS 2 traces are colour-coded to easily identify the probe wavenumber (green: a 4 , 3205 cm À1 ; red: a 5 , 2890 cm À1 ; blue: a 7 , 1879 cm À1 ), which is also indicated by coloured arrows.The choices of the probe wavenumbers were motivated by a previous assignment of these bands to isomers C (a 4 and a 7 ) and R (a 5 ) based on AIMD simulations. 3Depending on the laser used for the probe pulse (either OPO/OPA or FEL) we applied different experimental schemes to obtain the isomer specific spectra (see Experimental section), explaining the differences in the signal-to-noise ratios of the traces in Fig. 4. The measured spectra show no distinct differences, thus providing strong evidence for the presence of a single, spectroscopically unique species and hence a single constitutional isomer under the present experimental conditions.

Empirical band assignments based on isomer B
The persistent bands across both isotopologues are readily assigned in the context of those expected for isomer B.      22 The weak features just above the highest energy free O-H (O-D) stretches (see † band in Fig. 5) are often observed in cationic water clusters.
They originate from excitation of the frustrated rotations of the single acceptor water molecules in combination with the antisymmetric O-H stretch. 44he strongly isotope-dependent feature, the a 9 /b 9 pair near 1400 cm À1 , is much more intense in the all-H isotopologue.This is important because, as mentioned above, the assignment of a 7 and a 9 has been controversial, with VSCF theory supporting the assignment of a 7 to the OH stretch of the embedded Eigen motif 6 in B, while a VPT2 treatment 14 suggested that a 9 was a more likely candidate.In either case, one of these features remains unassigned.The persistence of the a 7 /b 7 pair, however, strongly argues for that feature as the bonded OH stretch, because it appears with the shift and intensity expected for an anharmonic fundamental, while b 9 is much weaker in the spectrum of the heavy isotope.This suggests that the transitions contributing to band a 9 predominantly arise from anharmonicities that are suppressed as the range of the vibrational zero-point wavefunction is reduced upon deuteration.The reduction of anharmonicities is also reflected by the difference in the scaling factors obtained for the two isotopologues (See Table 1 and Table S4, ESI †).When comparing the calculated values at the same level of theory for the two species, the method-related contributions to the scaling factor cancel out.The result is an evaluation of the change in anharmonicity upon deuteration, which amounts to B20-30% smaller anharmonicities in the heavier isotopologue.
a Bands position are given with an error of AE5 cm À1 .b B3LYP-D3/aug-cc-pVTZ harmonic wavenumbers of the all-H (all-D) isotopologue, scaled by 0.962 (0.973) above and 0.908 (0.927) below 1500 cm À1 (1100 cm À1 ).c Anharmonic VPT2/B3LYP/6-311++G(3df,3pd).d Only combination bands with intensities higher than 2 km mol À1 are considered.e Tentative assignment.Temperature-dependent IRMPD spectra: the role of entropy In order to probe how the anharmonicity-induced bands respond to temperature, we measured vibrational spectra at five ion trap temperatures in-between 50 and 250 K. Fig. 5 shows the temperature-dependent infrared multiple photon photodissociation (IRMPD) spectra of H + (H 2 O) 5 in the O-H stretching region (2650-4000 cm À1 ).The IRMPD spectrum measured at 50 K is similar to the IRPD spectrum of D 2 -tagged H + (H 2 O) 5 obtained at 13 K, confirming that the D 2 -tag has a minor influence on the vibrational spectrum as well as the energetic ordering of the lowest-lying isomers.At higher temperatures, band broadening due to thermal rotational and low-energy vibrational excitation is observed leading, for example, to a coalescence of bands a 1 and a 2 above 150 K. Interestingly, the combination bands, labelled with an asterisks and a dagger in Fig. 5, gain substantially in relative intensity with increasing temperature, which is consistent with their oscillator strengths being derived from large amplitude displacements.Here, increasing temperature plays a role analogous to that attributed to isotopic modulation of the displacements available in the vibrational ground state.[7][8] Assignment of the vibrational spectrum of D + (D 2 O) 5 ÁD 2 Much of the confusion regarding the assignment of the vibrational spectra of protonated water clusters arises from the fact that standard quantum computational approaches, due to their approximate nature, do not fully recover the pronounced anharmonic effects in these systems. 45With this in mind, we first evaluate which single isomer yields the best fit vibrational spectrum before discussing the remaining discrepancies, and then assess if this is a reasonable choice regarding its thermodynamic and kinetic stability.
6][7][8] Therefore, we focus on the description of the IRPD spectrum of the fully deuterated isotopologue here.Fig. 6 shows a comparison of the IRPD spectrum of the D 2 -tagged perdeuterated water pentamer (top trace) with the B3LYP-D3/aug-cc-pVTZ harmonic vibrational spectra of the B, R and C isomers (see Fig. 1) in the range from 215 to 3000 cm À1 .The experimental spectrum shows seventeen distinct absorption features labelled with b 1 to b 17 .Band positions and assignments are summarized in Table 1.The calculated harmonic spectrum of B provides the best fit, reproducing nearly all the main features observed in the experiment, in particular the new experimental data in the far IR region (see also Fig. 8).Once again, this discrepancy highlights the fact that the structure associated with these motions is not well described within either  1 for band labels and assignments.are better reproduced by the anharmonic calculations.Summarizing, the anharmonic calculations support our spectral assignments based on the presence of a single constitutional isomer and provide further insight into some of the details of the vibrational spectrum of D + (D 2 O) 5 .However, besides overestimating the exceptional red shift of band b 7 , two low energy fundamental modes, as well as their overtones and combination bands, are predicted to have negative energy.Furthermore, the obtained intensities of two other low energy modes are unreasonably high (on the order of 200 000 km mol À1 ).The VPT2 calculation showed poorer performance when applied to the light isotopologue or employing the aug-cc-pVTZ basis set, resulting in additional negative energy and high intensities.The a 7 mode is most unstable with respect to correction at the VPT2 level, raising the possibility that the levels involving this motion require inclusion of higher order terms in the potential to converge on the experimentally observed band location.A comparison of the harmonic and VPT2 anharmonic frequencies and intensities obtained with the two basis sets is reported in Tables S5 and S6 (ESI †).
A conformational isomer of B, referred to as B2, with an almost identical vibrational spectrum (see Fig. 8), is obtained  by rotating the 1-AD/2-A waters by roughly 180 degrees around an axis parallel to the hydrogen bond between the D 3 O + and the 1-AD water (see Fig. 2).The two branched rotamers are approximately isoenergetic and the low barrier between them can easily be overcome, even at the low temperatures of the experiment, as is supported by the kinetic model (see below).The system is thus expected to be so fluxional that it can explore all the configurations connecting the two structures.This fluxionality is likely the cause of the broadening of bands b 10 and b 11 .The corresponding vibrational modes show a shift in wavenumber of 39 and 37 cm À1 , respectively, upon isomerization via TS B-B2 (see Fig. 8).We thus refine our assignment of the structure of D + (D 2 O) 5 and note the fluxional character of the branched structure.
The assignments of the fundamentals in the all-D isotopologue spectrum are useful in unravelling the origin of the more complex band patterns in the all-H system.Because all of the bands in the latter are demonstrated by IR 2 MS 2 to arise from the same carrier, the ''extra'' bands in the light isotopologue must be due to overtone and combination bands induced by anharmonic coupling.For example, the feature labelled with an asterisk denotes a combination band and, as mentioned above, is present only in the IRPD spectrum of H + (H 2 O) 5 ÁD 2 .Anharmonic VPT2 calculations on the all-H isotopologue predict several active combination bands in this region.The main contribution stems from the combination of the hydrogen-bonded O-H stretch 1-AD-H 2 O (a 4 ) with two low energy modes.One of these is the frustrated rotation of a 1-A-H 2 O combined with the hydrogen bond stretch between the 1-AD and the 2-A water molecules and the second is the frustrated rotation of the 2-A water.In the spectrum of the perdeuterated species only the latter combination band has some IR intensity but falls under the broad feature b 4 .

RRKM calculations and kinetic modelling
The relative B3LYP and CCSD(T) energies of the three lowest energy isomers of the protonated water pentamer B, R and C and the first-order transition states TS R-B and TS B-C , connecting the corresponding potential minima, are listed in Table 2 (see also Table S1, ESI †).For the branched isomer, two nearly isoenergetic structures, labelled B and B2 in Fig. 2 and Table 1, exist, separated by a shallow barrier (TS B-B2 : 1.0 kJ mol À1 ).As we show later, only these conformational isomers B and B2 are predicted to be kinetically stable at low temperature, while the constitutional isomers R and C are not.R represents the global minimum energy structure on both the B3LYP and presumably the CCSD(T) potential energy hypersurfaces (see Table 2).However, B is stabilized relative to R, when (harmonic) ZPE is taken into account reducing the energy difference to less than 1 kJ mol À1 , making them nearly isoenergetic.While the level of theory used here is not sufficiently accurate to predict the energetic ordering of B and R at 0 K, it is apparent that B is entropically favoured over R at higher temperatures (see Gibbs free energies DG 0 , equilibrium constants K 0 and Boltzmann weighted populations P 0 in Table 2 and Tables S2, S3, ESI †).Moreover, the free energy difference between B and C remains nearly constant and larger than 10 kJ mol À1 , indicating that C can only be spectroscopically a ZPE is calculated using scaled harmonic B3LYP-D3/aug-cc-pVTZ wavenumbers.b The thermal energy contributions are calculated at 0.001 atm and with a vibrational scaling factor of 0.962 (0.973) for the all-H (all-D) isotopologue.Note the different scaling below 1500 in Table 1.This journal is © the Owner Societies 2016 probed at low temperature if it is kinetically trapped behind the B5 kJ mol À1 barrier.The RRKM canonical rate constants are shown in Fig. 9: panel (a) depicts the results obtained from the quantum treatment (i.e.quantum density and sum of vibrational states and inclusion of ZPE) while panel (b) displays the ones from the classical approach (i.e.classical density and sum of rovibrational states and no inclusion of ZPE).In particular, rate constants are reported as a function of the temperature for the forward and backward elementary reactions of the three isomerization pathways considered (see Fig. 2).Microcanonical rate constants are reported in Fig. S6 (ESI †).First, we analyse the quantum case.It is helpful to compare the two rate constants for opposite directions of the same isomerization process at the same temperature.For example, the fact that k 1 exceeds k À1 by one to three orders of magnitude indicates a preference for the branched isomer B over the ring isomer R across the entire energy range considered.On the other hand, similar values for k 2 and k À2 suggest that the conformational isomers B and B2 will rapidly interconvert.k 3 and k À3 show the largest difference, in particular at low energies (or temperature), revealing that the chain isomer C will convert nearly completely to the branched conformers B and B2.The main difference between the classical and the quantum treatment is found in the low temperature region, where the classical rate constants are smaller.As expected, the values of the classical and quantum rate constants approach each other with increasing temperature.
The most striking difference between the classical and quantum RRKM calculations manifests itself in the rate constants for the R/B isomerization reaction, which leads to a qualitatively different kinetic stability at low energies (temperatures).In the absence of ZPE R, the global minimum energy structure, is also the kinetically preferred isomer, which is reflected in a larger isomerization rate constant for the back reaction (k À1 4 k 1 ) forming R at low energies (temperatures).The rate constants are equal at 40 kJ mol À1 (111 K) and above this energy (temperature) the forward reaction leading to B is preferred (k 1 4 k À1 ).
The kinetic stability of each isomer at 50 K is evaluated in Fig. 10 for the quantum (upper row) and the classical case (lower row).Independent of the initially prepared constitutional isomer, Fig. 10 The kinetic stability of the isomers R, B, B2 and C is tested with a kinetic model that uses either quantum (top row) or classical (bottom row) rate constants obtained from RRKM calculations at 50 K.Boltzmann weighted populations at 50 K obtained from relative Gibbs free energies (see Table S2 R (left column), B/B2 (center column) or C (right column) the thermodynamic limit is reached within a few tens of nanoseconds with the exception of the quantum case for isomer C, where it is reached in several microseconds.B/B2 are predominantly formed for the quantum case, while interconversion to R is dominant in the classical case.
We also tested the effect of the quality of the potential energy surface on the RRKM calculation by using the energies obtained without dispersion correction or with different basis sets (e.g.6-311++G(3df,3pd)).This leads to small changes in the relative energies and thus the respective equilibrium constants, which in return can lead to slightly larger numbers for the population of R in the quantum case.The equilibration times are also affected, but overall the results remain in qualitative agreement.A more thorough description of the details and the results of the RRKM calculations, also involving higher energy isomerization reactions, lies outside the scope of the present study and will be part of a separate publication.
Summarizing, the results of the RRKM calculations lead to two important insights.First, the kinetic stability of the considered isomers below 100 K changes dramatically, when nuclear quantum effects in the form of ZPE are considered.Second, isomer C is not kinetically stable (on a microsecond timescale) at 50 K.

Summary and conclusions
We have shown that the IRPD spectra of the D 2 -tagged protonated water pentamer and its fully deuterated isotopologue probe the same spectroscopic species that is most consistent with a branched structure, consisting of two conformational isomers that interconvert readily on the time scale of the experiment.
Harmonic DFT calculations of the branched isomer B reproduce the IRPD spectra of the D 2 -tagged H + (H 2 O) 5 and D + (D 2 O) 5 with reasonable agreement, with the exception of bands a 7 /b 7 and a 9 /b 9 .The strong feature a 9 just below the intramolecular HOH bending mode in H + (H 2 O) 5 , previously attributed to an H-bonded OH fundamental of the Eigen core, 14 is dramatically suppressed in the perdeutero isotopologue, indicating that it contains significant contributions arising from anharmonic coupling, while the higher energy feature (a 7 /b 7 ) is the fundamental associated with this oscillator.Anharmonic VPT2 calculations accurately account for most of the features across the D + (D 2 O) 5 spectrum.However, this approach overestimates the pronounced red shift of mode b 7 .It also performs considerably worse in reproducing the H + (H 2 O) 5 spectrum, since anharmonic effects are more pronounced here as a result of the larger displacements in both the ground and excited vibrational states.In summary, higher level anharmonic calculations are necessary to achieve satisfactory convergence between experiment and theory for the IR spectra of the protonated water pentamer, similar to the archetypal case of H 5 O 2 + . 45simple kinetic treatment shows that for a reliable understanding of the equilibrium dynamics between the various isomers of the protonated water pentamer, nuclear quantum effects must be taken into account at low temperature.Hence, conclusions derived from standard AIMD simulations, which make use of quantum mechanics (at the DFT level) for the electronic structure, but treat nuclear motion in a pure Newtonian way, can be unreliable when dealing with (i) effects that result from the quantum nature of protons and/or (ii) kinetically unstable isomers.A step closer to convergence between experiment and theory could certainly be achieved by the application of path integral molecular dynamics simulations (as well as related methods) [46][47][48][49][50] to the present system, thus including the quantum dynamics of the nuclei.

Fig. 1
Fig. 1 B3LYP-D3/aug-cc-pVTZ minimum-energy structures of the three constitutional isomers discussed in the text.Intramolecular hydronium ion O-H and intermolecular hydrogen bond lengths are reported in pm.Water molecules are classified according to their location in the first (1) or second (2) hydration shell of H 3 O + , and their function as hydrogen bond acceptor (A) or donor (D).

Fig. 2
Fig.2B3LYP-D3/aug-cc-pVTZ minimum-energy and TS structures of the protonated water pentamer discussed in the text.The energy diagram reports relative energies DE 0 , determined from CCSD(T) energies plus DFT ZPEs, and relative Gibbs free energies (DG 0 ), determined from CCSD(T) energies plus thermal corrections to Gibbs free energies at 50, 150 and 250 K.All energies are given in kJ mol À1 (see Table2).The bottom panel illustrates the sequence of the equilibrium reactions assumed for the RRKM calculations.
Fig.2B3LYP-D3/aug-cc-pVTZ minimum-energy and TS structures of the protonated water pentamer discussed in the text.The energy diagram reports relative energies DE 0 , determined from CCSD(T) energies plus DFT ZPEs, and relative Gibbs free energies (DG 0 ), determined from CCSD(T) energies plus thermal corrections to Gibbs free energies at 50, 150 and 250 K.All energies are given in kJ mol À1 (see Table2).The bottom panel illustrates the sequence of the equilibrium reactions assumed for the RRKM calculations.
Three sharp peaks with characteristic spacing and relative intensities (a 1 -a 3 /b 1 -b 3 ) are found at highest energy in the free O-H/O-D stretch region.Further to the red, the next four features (a 4 -a 7 /b 4 -b 7 ) are attributed to hydrogen-bonded (H-bonded) O-H/O-D stretches and appear broader, as is typical for modes that contain embedded OH groups.The sharp transitions (a 8 /b 8 ) are associated with the water bending modes.Characteristic triplets (a 10 , a 11 , a 13 /b 10 , b 11 , b 13 ) emerge in the high energy part of the libration region, of which band a 10 /b 10 corresponds to the umbrella mode of the Eigen core (see

Fig. 3
Fig. 3 IRPD spectra of D 2 -tagged isotopologues of the protonated water pentamer: D + (D 2 O) 5 (top) and H + (H 2 O) 5(bottom).For better comparability the spectra are plotted against opposing ordinates and the abscissa of the top spectrum has been scaled by 1.36 (see text).Bands coloured red indicate modes assigned to the hydronium ion.See Table1for band assignments.

Fig. 4
Fig. 4 IRPD (top row) and IR 2 MS 2 (middle and bottom row) spectra of D 2tagged H + (H 2 O) 5 from 1000 to 2100 cm À1 and from 2600 to 3900 cm À1 .Probe energies are indicated by vertical coloured arrows and IR 2 MS 2 traces are coloured accordingly.
Three sharp peaks, b 1 , b 2 and b 3 , are found at the high energy end of the spectrum, in the free O-D stretching region.The band positions and relative intensities of these modes are very sensitive to the environment of a particular water molecule.With this in mind, it is worth noting that the calculated harmonic spectrum of B is the only one that reproduces the characteristic spectral pattern of these bands.Features b 4 , b 5 , b 6 and b 7 correspond to the hydrogen-bonded O-D stretches.b 4 arises from the O-D stretch of the 1-AD water.The symmetric and antisymmetric combination of the two O-D stretches of the D 3 O + core towards the 1-A waters are responsible for peaks b 5 and b 6 .b 7 corresponds to the third O-D stretch of the D 3 O + core (towards 1-AD-D 2 O) and shows the largest deviation (+235 cm À1 ) from the experiment.

Fig. 5
Fig. 5 (b-f) Temperature-dependent IRMPD spectra of H + (H 2 O) 5 in the O-H stretching region recorded by varying the trap temperature in 50 K steps from 50 to 250 K. (a) The D 2 -predissociation spectrum at 13 K is reported for comparison.See Table1for band labels and assignments.

Fig. 6
Fig. 6 IRPD spectra of D 2 -tagged D + (D 2 O) 5 in the region from 215 to 3000 cm À1 compared to the scaled (0.922 below 1100 cm À1 and 0.974 above) B3LYP-D3/aug-cc-pVTZ harmonic spectra for the three isomers B (branched), R (ring) and C (chain) discussed in the text.Traces below 1400 cm À1 and above 2580 cm À1 are shown with a vertical magnification of Â5 for better visibility.

Fig. 7
Fig. 7 IRPD spectrum of D 2 -tagged D + (D 2 O) 5 in the region 215 to 3000 cm À1 (top trace) compared to the anharmonic (middle trace) and the scaled (0.92/0.97) harmonic spectrum (bottom trace) of the branched isomer B calculated at the B3LYP/6-311++G(3df,3pd) level of theory.Combination bands and overtones in the VPT2 spectrum are coloured in red and blue, respectively.Traces below 1085 cm À1 and above 2600 cm À1 are shown with a vertical magnification of Â5 for better visibility.

Fig. 8
Fig. 8 IRPD spectrum of D 2 -tagged D + (D 2 O) 5 in the region 215 to 1000 cm À1 compared to the scaled B3LYP-D3/aug-cc-pVTZ harmonic spectra for the two versions of the branched isomer and the TS between them.

Fig. 9
Fig. 9 Quantum (a) and classical (b) canonical rate constants obtained with RRKM calculations for the forward and backward elementary reactions of the three isomerization pathways considered (see Fig.2).
Fig.10The kinetic stability of the isomers R, B, B2 and C is tested with a kinetic model that uses either quantum (top row) or classical (bottom row) rate constants obtained from RRKM calculations at 50 K.Boltzmann weighted populations at 50 K obtained from relative Gibbs free energies (see TableS2, ESI †) are indicated by colour-coded triangles.

Table 2
Relative B3LYP and CCSD(T) electronic energies (in kJ mol À1 ) without ZPE (DE), CCSD(T) electronic energies with ZPE (D E 0 ), and Gibbs free energies (DG 0 ) of the B3LYP-D3/aug-cc-pVTZ minimum energy and first order transition state structures considered in this work