Chemical Physics 514 (2018) 164–175
Contents lists available at ScienceDirect
Chemical Physics
journal homepage: www.elsevier.com/locate/chemphys
The protonated water trimer and its giant Fermi resonances ⇑ Nagaprasad Reddy Samala, Noam Agmon
The Fritz Haber Research Center, Institute of Chemistry, The Hebrew University of Jerusalem, 91904, Israel article info abstract
Article history: The protonated water trimer is a prototype ‘‘proton wire” that can transport two protons nearly concert- Received 8 January 2018 edly. This ‘‘proton transfer mode” (PTM) is an important contributor to the infrared spectrum of the iso- In final form 1 April 2018 lated gas-phase cluster. We have simulated its infrared spectrum for both the hydrated and deuterated Available online 3 April 2018 isotopologues, using vibrational 2nd order perturbation theory (VPT2) and ab initio molecular dynamics (AIMD) trajectories. VPT2 calculations explain quantitatively the experimental spectra at both high and Keywords: low frequencies, provided that high-level quantum chemistry is utilized. In the D2-tagged hydrated Cluster cluster, the PTM undergoes giant Fermi resonances (FR’s) with two combination bands. In the deuterated Fermi resonance analogue, one observes a single FR of ‘‘normal” intensity, manifested as the doublet recently Infrared Proton reported experimentally. We provide band assignment for both isotopologues, with and without the Water D2 tag applied in experiment. We discuss possible manifestations of the giant resonance on proton transfer through water wires. Ó 2018 Elsevier B.V. All rights reserved.
1. Introduction of occurrences involve longer-range concerted multi-proton hops [8,9]. Such events are prominent in trajectories of hydronium- Protons in water diffuse abnormally fast, which is attributed to hydroxide recombination [10], and are manifested experimentally the Grotthuss mechanism [1]. In bulk water, this is depicted as an in the inverse temperature (T) dependence of the rate coefficient þð Þ [11] for PT between photoacid and base in liquid water along ‘‘isomerization” reaction from a solvated hydronium, H3O H2O 3 (the ‘‘Eigen cation”), centered on the donor oxygen, via a shared- chains of hydrogen-bonded water molecules (‘‘water-wires”) þð Þ [12]. The shortest water wire beyond the Zundel cation is the pro- proton intermediate, H H2O 2 (the ‘‘Zundel cation”), to an Eigen cation localized on the acceptor oxygen [2]. While the rate limiting tonated water trimer. step involves hydrogen-bond rearrangement in the 2nd solvation Infrared (IR) studies demonstrated that traces of aqueous acids þ + ð Þ shell of the H3O , the actual proton transfer (PT) step in the Zundel in organic solvents consist predominantly of H H2O 3 cations cation can be as short as a vibrational period, or involve barrier [13], rather than the traditional Eigen or Zundel ions. In narrow recrossings during several vibrational periods. Thus, vibrational carbon nanotubes, which host a single-file of water molecules, PT dynamics plays a role in the kinetics of proton mobility. may involve isomerization between two Zundel cations via a tran- þð Þ Indeed, low temperature vibrational predissociation spectra of sient H H2O 3 structure [14]. In the green fluorescent protein þ (GFP) a water-trimer was found to connect the chromophore to the H5O2 Zundel cation show an intense ‘‘proton transfer mode” (PTM) at 1047 cm 1 [3], arising from the proton shuttling between the bulk [15]. It may efficiently conduct the proton, generated the two oxygen centers, while at 928 cm 1 appears a strong com- upon chromophore photoexcitation, outside the protein barrel bination band, gaining its unusual intensity from the PTM via a structure. Consequently, it is important to understand the dynam- ‘‘Fermi resonance” (FR) [4]. This combination, of OO stretch and a ics of the protonated water trimer and, particularly, its IR spectrum + wag overtone [5], tends to localize the proton as H3O on one side [16–23]. of the Zundel cation, at the expense of the shared proton mode. While it is easy to identify the IR absorption due to the flanking Here we propose that FRs play an important role also in larger pro- water free OH stretch, that of the hydrogen bonded OH (HB–OH) þ stretches is more elusive. Schwarz suggested that some small tonated water clusters, specifically the H ðH2OÞ cluster. 3 1 While the stepwise Eigen-Zundel-Eigen mechanism may absorption in the range 2000–2300 cm might be due to the þð Þ describe the majority of PT events in liquid water [6,7], a minority H H2O 3 cluster [16]. Subsequently, Okumura et al., who first applied vibrational predissociation spectroscopy to this cluster, identified a peak near 2300 cm 1 as originating from the HB–OH ⇑ Corresponding author. stretch [17]. Headrick et al. [18] extended the accessible spectral E-mail address: [email protected] (N. Agmon). https://doi.org/10.1016/j.chemphys.2018.04.003 0301-0104/Ó 2018 Elsevier B.V. All rights reserved. N.R. Samala, N. Agmon / Chemical Physics 514 (2018) 164–175 165 range to below 2000 cm 1, locating a very strong peak at 1880 cm 1. Although they have also observed a broad weak feature near 2300 cm 1, they have assigned the 1880 cm 1 band to the HB–OH asymmetric stretch (as) mode [18] (denoted below HB–OHas). This concerted 2-proton motion is depicted in Fig. 1 (that also defines a numbering scheme for the atoms), and may be viewed as the PTM responsible for concerted PT in the liquid phase. From its high intensity one might expect facile concerted PT whenever a þð Þ H H2O 3 grouping is formed. As we argue below, due to the FRs this may not always be true. While the above assignment prevailed for a decade [19,20], the- ory has recently cast some doubt on this simplistic picture. Notice- ably, Duong et al. [22] have improved upon the analytical multidimensional potential energy surface (PES) and dipole surface of Yu & Bowman [21], for which a vibrational self-consistent field (VSCF) analysis combined with vibrational configuration interac- þð Þ tion (VCI) of the untagged (‘‘bare”) H H2O 3 cluster was per- formed. They observed that the dominant 1880 cm 1 band arises ‘‘from a congested series of transitions to ...mixed levels” [22]. This is lamentable, because a major goal in performing vibrational anal- ysis is to provide the physical insight gleaned from clear-cut assignments. Interestingly, in a recent M.Sc. thesis [24], Shtainmetz suggested a FR of the HB–OHas with a combination band involving the two normal modes (NM’s) depicted in Fig. 2. Even more recently, it became possible to measure this spec- trum in the intermolecular, terahertz regime (below 1000 cm 1) using a free electron laser [23]. The most intense low frequency bands were assigned to the water wag and hydrogen bond stretch. However, it was no longer possible to treat this regime with VSCF/ VCI, so quasiclassical molecular dynamics (QCMD) was used instead for the same analytical PES. The challenge then is to explain the whole spectrum from one theory.
Similarly intriguing is the novel D2-tagged IR predissociation þð Þ spectrum for the deuterated analogue, D2 D3O D2O 2 [22]. Here 1 the 1878 cm PTM band (Fig. 3a) is replaced by a doublet, at 1680 Fig. 2. The hydronium (a) umbrella and (b) rocking modes (harmonic B3LYP-D3/6- and 1452 cm 1 (Fig. 3b). Unfortunately, a simple assignment was 311++G⁄⁄) are (a) symmetric and (b) antisymmetric to reflection. They combine and again not feasible, and the isotope effect was notably smaller than interact with the PTM of Fig. 1 to produce the unusually strong combination bands expected. It was consequently suggested that ‘‘calculations that are in the spectra shown in Fig. 3. simpler than the present VSCF/VCI ones (e.g. ...vibrational second order perturbation theory (VPT2)) could be performed to shed light on the unusual spectral features seen in these clusters” [22]. The problem with the VPT2 approach is that bands involving The present work lives up to this challenge by testing VPT2 the excess proton motion appear to be unstable to perturbations, from Gaussian 09 using 2nd-order Møller-Plesset perturbation the- varying sometimes wildly with method and bs. This has been ory (MP2) [with a frozen core (fc)] and density functional theory clearly demonstrated for the protonated water dimer and tetramer (DFT) with several density functionals and basis sets (bs’s) for [33]. Consequently, it has been postulated that for the trimer þ ‘‘VPT2 as implemented in Gaussian 09 will not even qualitatively the H ðH2OÞ cluster over the whole frequency range, as well as 3 reproduce the essential physics controlling the vibrational band for its deuterated analogue. These calculations have been per- pattern” [22]. It turns out that the recent availability of detailed formed with and without the D2 tag. experimental spectra [22,23] now allows to select, from a plethora of diverse spectra [24], those VPT2 calculations that reproduce the experimental results. These, in turn, provide the lacking under- standing of the essential physics controlling the vibrational band pattern.
The calculations with the D2 tag are the closest to experiment, suggesting that the 1878 cm 1 band is not due to the PTM funda- mental (which remains where Schwarz has predicted), but rather + to a combination of the H3O umbrella and rocking modes (Fig. 2), which enters into an abnormally strong FR with the PTM [24]. As in previous studies of protonated water clusters [25–33],we have also performed ab initio molecular dynamics (AIMD) simula- tions, in which the forces are calculated from DFT (using the B3LYP-D3 functional), while the atoms are propagated classically, using Newton’s equations of motion. This method [34] thus lacks nuclear quantum effects (NQE), which we replace by frequency Fig. 1. The hydronium asymmetric stretch NM (termed herein HB–OHas or PTM) of þð Þ ⁄⁄ scaling, and specialized analysis is required to infer about mode the H H2O 3 cluster (harmonic B3LYP-D3/6-311++G ), with atom labels indicated. Note the antisymmetry to reflection. couplings [27,35]. With this, we obtain a consistent interpretation 166 N.R. Samala, N. Agmon / Chemical Physics 514 (2018) 164–175
(a) Hydrated isomer (with D tag)
2 3726 Fundamental 1883 MP2/aug-cc-pVTZ
Overtones 1878 3580 Combination bands 3640 2084 Expt 2284 1943 3540 2383 2109 1534 3726 1639 1059 2307 3814 2961 2410 2693 1992 2367 3628 2544 2664 2175 1553 1080 1621 1907 992 3039
1000 1500 2000 2500 3000 3500 4000 1680 Intensity (b) Deuterated isomer (with D tag) 1452 2 1740 1770 2631 2770 2660 800 1477 1202 1791 1143 1991 2631 2149 2314 2978 1540 1913 818 2063 2648 2767 1424 3038 1186 731 1145
1000 1500 2000 2500 3000 3500 4000
-1 Wavenumber (cm )
þð Þ þð Þ Fig. 3. The experimental D2-tagged IR predissociation spectrum of (a) the H H2O 3 cluster, and (b) the D D2O 3 cluster at 12 K (blue lines, digitized from Ref. [22]), compared with the corresponding (unscaled) VPT2 spectra calculated at the MP2(fc)/aug-cc-pVTZ level, with the tag at the H3 atom of the hydronium core (see Fig. 11b in Ref. [17]). Black, red and green sticks mark fundamental, overtone and combination bands, with their frequencies in cm 1. Omitted from the figure are all combination bands involving a fundamental with a negative or near-zero frequency (modes 22 and 23 for both clusters). For visibility, the fundamental intensities in panel (a) were scaled by a factor of 2. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)
for both VPT2 and AIMD results that sheds new light on the infra- The harmonic/anharmonic frequencies and intensities, for red spectrum of this important cluster. fundamental and combination bands, are collected in the Excel file of the Supplementary Material (SM). The Excel sheet 2. Theoretical methods names identify the tagging, deuteration, method and bs. For exam- ple, mp2-avtz-D2-H indicates that the IR spectrum of the hydrated 2.1. VPT2 calculations tagged cluster was computed with MP2/aug-cc-pVTZ. Similarly, b3lyp-pople-D indicates that the IR spectrum of the deuterated The calculation of the anharmonic VPT2 spectra was performed bare cluster was computed with B3LYP-D3/6-311++G⁄⁄. On the using the Gaussian 09 (Revision D.01) suite of programs [36]. Two top of each sheet appear the accuracy conditions. methods were used: (a) MP2(fc), a quantum chemistry method The above mentioned calculations were initially performed only that accounts for electron correlations in the valence shell, and on the bare clusters (no tag), because experiment showed notice- (b) DFT with a variety of exchange–correlation functionals. The able effect of the tag only for the O3H3 vibration [22]: It appears functionals tested, B3LYP-D3, M06-2X, mPW2-PLYP and xB97, between the flanking water symmetric stretch (ss) and as modes 1 were chosen from a recent benchmarking of DFT functional perfor- of the bare cluster, moving some 90 cm to the red when a D2 mance in predicting harmonic frequencies of water clusters [37], molecule gets attached at this position. where D3 is Grimme’s third generation dispersion-correction In a second phase, we have performed additional calculations ð Þ [38]. For the H2O 3 cluster, the mean absolute deviation (MAD) with explicit inclusion of the D2 tag on the H3 (D3) atom of the + + with these functionals from the CCSD(T) near complete basis set H3O (D3O ) core [17]. Here we applied what appeared to us the limit was 31, 13, 14 and 20 cm 1, respectively, as compared to most promising of the above methods, MP2(fc) and B3LYP-D3. 13 cm 1 with MP2 [37]. The above methods were combined here These additional calculations were performed because, in-spite of ⁄⁄ with either Pople’s 6-311++G or Dunning’s aug-cc-pVTZ bs’s. All an apparently limited tag effect on the experimental spectrum, þð Þ these calculations were performed for both the H H2O 3 and the high sensitivity of perturbation expansion methods (such as þð Þ VPT2) to the environment might lead to more noticeable changes D D2O 3 clusters. þð Þ in the computed spectrum. The geometry of the H H2O 3 cluster was first optimized with the tight convergence criterion of Gaussian 09, confirming that the For the calculations showing the most promising agreement obtained structures are stable minima on the PES. Deuterated iso- with experiment (B3LYP-D3 and MP2 with the Dunning bs), we þð Þ have performed additional tests to check the numerical stability mers were reoptimized starting from the optimized H H2O 3 coordinates (which did not change). The harmonic NM spectra of the results. One set of calculations used the ‘‘VeryTight” instead 12 were then calculated at each of the optimized geometries with of ‘‘Tight” optimization condition, and 10 instead of the default the corresponding quantum method/bs. Anharmonic analysis was 10 10 integral accuracy. For B3LYP-D3 we have checked the effect subsequently performed with VPT2. of using a tighter, ‘‘Ultrafine” integration grid for the exchange N.R. Samala, N. Agmon / Chemical Physics 514 (2018) 164–175 167 functional. Finally, we have checked the effect of using a 3. Results ‘‘FourPoint” numerical differentiation formula instead of the default two-point differentiation. Results from these additional calcula- þð Þ We begin by presenting the results for the H H2O 3 cluster, tions are also collected in the SM Excel file. The results show that using VPT2 and AIMD, then proceed to its deuterated analogue. the very low frequency modes (below, say, 200 cm 1) are not con- The spectra are reported for both the high and low frequency vibra- sistently reproduced, while the variations in the higher frequency tional modes, without and with explicit inclusion of the D2 tag. This 1 1 modes is modest (under 15 cm with the D2 tag, and 5 cm with- exceeds the scope of other recent computations [22,23], which did out). Moreover the improvement is not systematic. Consequently, not include the tag, and reported VSCF results only at high frequen- we report below the results obtained with the Gaussian defaults. cies, and QCMD data only at low frequencies. All harmonic and anharmonic spectra are collected in the Excel file of the SM.
2.2. AIMD setup 3.1. The hydrated cluster Classical AIMD simulations were performed in the CP2K/Quick- step software package [39], see http://www.cp2k.org/quickstep. 3.1.1. Harmonic normal modes This code describes the electronic structure using DFT with a Fig. 4 shows the harmonic IR spectra of this cluster using MP2 mixed Gaussian and plane wave basis. We used the B3LYP-D3 and DFT/B3LYP-D3 with the Dunning bs. As expected from the functional with Grimme’s third generation dispersion-correction water cluster benchmarking [37], the methods indeed give similar (D3) [38] and the augmented triple-f valence polarization (aug- results. Around 3800 cm 1 we have the free OH stretches, around TZVP) bs, both of which were not utilized in previous AIMD com- 2500 cm 1 appear the HB–OH stretches, between 1600 and 1700 putations for this cluster [25–27]. The plane wave energy cutoff cm 1 the 4 bending modes, and near 1250 cm 1 – the umbrella was set to 350 Ry. Self-interaction correction was applied with mode. Animations of all the fundamental modes are shown in the Martyna-Tuckerman Poisson equation solver [40]. The orbital the Word file of the SM. transformation method was applied for faster convergence, with The normal mode visualization allows several observations. The the convergence criterion of 1 10 6 a.u. at every timestep (0.5 fs). isolated water asymmetric stretch (as) and symmetric stretch (ss) The trajectories generated in CP2K for this study are summa- modes each combine into two nearly degenerate modes, in which rized in Table 1. In the optimization stage, initial coordinates were the two water molecules vibrate either in phase or out of phase. + selected from the MP2 optimized geometries. In the equilibration The HB–OHas mode involves also some rocking of the H3O free stage, these structures were run in the canonical (NVT) ensemble, OH. The 4 bending modes are clearly divided into 2 symmetric with a target temperature of 50 K maintained by the Nosé-Hoover bends (bs) with respect to a perpendicular mirror plane, and 2 chain thermostat. From each NVT trajectory, we branched out three asymmetric bends (ba). The higher energy ba and bs pair is often + NVE production trajectories (at 9, 9.5 and 10 ps). These were con- ascribed to the H3O core whereas the lower energy pair arises tinued at constant energy for 20 ps each. While the measurements from the flanking water molecules [27]. However, these bending in Fig. 3 were performed at 12 K, equilibration at such a low tem- modes are rather thoroughly mixed and their localized origin is perature is prohibitively long, and not much change is expected hard to discern. The umbrella mode (U) involves also in-phase between 12 and 50 K or, more accurately, up to the average temper- rocking of the 2 terminal water molecules (hence this mode is ature, hTi, along the NVE trajectories. Coordinates and velocities symmetric to reflection), whereas the hydronium rocking (R, near from each production run were saved every timestep (0.5 fs) to 1060 cm 1) and rotation (r, near 600 cm 1) modes involve out of ensure that the spectra at high frequencies do not get distorted. phase water rocking (asymmetric to reflection), see Fig. 2.D2 trans- The dipole moment was calculated using Wannier localization, lations and rotations couple into several low frequency modes. and then the IR absorption coefficient was computed from the There are, of course, no overtones and combination bands on the dipole moment (time derivative) autocorrelation function (DACF), harmonic level. by taking its temporal Fourier transform [34]. A frequency-depen- dent nuclear ‘‘harmonic quantum correction” factor was applied to the computed IR intensities (see Eq. (1) in Ref. [33]). The spectrum 3.1.2. VPT2: High frequencies obtained contains anharmonic effects that are generated as the tra- In comparison with the experimental spectrum in Fig. 3a,itis jectory moves on the (anharmonic) multidimensional PES. Noise clear that all the indicated harmonic peaks in Fig. 4 should red- was reduced by averaging together the normalized spectra from shift with inclusion of anharmonicity. The different methods agree the three trajectories that were started from different initial condi- on the extent of the red-shift of the observed fundamental bands, tions. Missing from the picture, of course, are the NQE because of except for the two HB–OH modes. For example, the PTM frequency the use of Newton’s equations. varies within a frequency window of ca. 600 cm 1 (Fig. 5). The In order to assign the vibrational frequencies, partial velocity accuracy of predicting the harmonic frequencies [37] seems to pro- autocorrelation functions (VACF) were calculated from each trajec- vide no guide on selecting the most accurate anharmonic method. tory, by restricting the Fourier transform to the velocity of the Another problem with a perturbation theory like VPT2 is that specified coordinate (when interatomic charge transfer can be for some modes it becomes so unstable that frequencies become neglected, the velocity is proportional to dipole moment negative or intensities become (nearly) infinite. When this occurs derivative). only for a few low frequency modes (e.g., modes 22 and 23 in
Table 1 AIMD (B3LYP-D3/aug-TZVP) trajectories run for this study.
Optimized cluster Equilibration Production type/name target T,K Dt; ps type # a trajectory name branching times, ps Dt,ps hTi,K þð Þ H H2O 3 NVT1 50 10 NVE 3 1, 2, 3 @ 50 K 9, 9.5, 10 20 68, 69, 89 þð Þ D D2O 3 NVT2 50 10 NVE 3 4, 5, 6 @ 50 K 9, 9.5, 10 20 60, 51, 48
a Number of trajectories of the given type. 168 N.R. Samala, N. Agmon / Chemical Physics 514 (2018) 164–175
+ and the VSCF/VCI calculation [22]. Clearly, the VPT2 results are in Harmonic IR spectra (a) H (H2O) (without D tag) 3 2 better agreement with experiment than any previous computation. For intensities, VSCF/VCI is possibly more reliable. MP2/aug-cc-pVTZ 2464 The assignment of most peaks becomes much clearer now. Par- B3LYP-D3/aug-cc-pVTZ 1 ticularly, band a8 at 1878 cm , long believed to be the PTM [18– 21], is actually a combination band, U þ R, of the hydronium umbrella (U) and rocking (R) modes. [For a non-linear H–O–H moi- ety, twist, wag and rock are defined as the arrested rotations around the symmetry axis, the perpendicular in-plane axis, and the perpendicular out-of-plane axis, respectively]. The HB–OH fun- 2620 1 damentals contribute near 2300 cm (band a5), as suggested by 3897 3789 403 254 Schwarz [16], and they are comparatively weak (note that we have 369 1265 3812 1077 1702 scaled the fundamentals in Fig. 3 by 2, for visibility). Band a6 is
again a combination band, r þ bs, of hydronium rotation (r) around
0 500 1000 1500 2000 2500 3000 3500 4000 its main symmetry axis and its symmetric bending (bs) mode. Band Intensity + (b) H (H2O) (with D tag) a7 is the least certain, and could be the U overtone.
3 2 2558 Note that the VPT2 computations only report frequencies for transitions involving up to two quanta. For higher order transitions one may consult the VSCF/VCI results (Table S2 in Ref. [22]). The only 3-quanta transition reported there is the third harmonic of the hydronium arrested rotation mode, and it occurs within the
envelope of band a8. 2693 3708
240 364 3.1.3. VPT2: Low frequencies 1057 1244 3901 1705 Low frequency measurements (below 1000 cm 1) are difficult, 3792 434 585 requiring the intensity of a free electron laser. Such measurements 0 500 1000 1500 2000 2500 3000 3500 4000 were just published [23]. ‘‘Unfortunately, in this range the quantum -1 approach becomes prohibitively difficult” so that VSCF/VCI results Wavenumber (cm ) could not be obtained, and classical dynamic results were reported þð Þ instead [23]. We find that also VPT2 has more problems at low fre- Fig. 4. Harmonic IR spectra for (a) the bare and (b) the tagged H H2O 3 cluster for two different quantum methods (indicated). Note that (a) with MP2 was previously quencies, such as negative anharmonic frequencies (while all the reported in Table 2 of Ref. [27], though with somewhat different assignments for modes have positive real harmonic frequencies, as expected for a the low frequency modes. minimum on the PES). Another problem is diverging intensities. However, not all is lost. From the many different computations, ⁄⁄ þð Þ we found that B3LYP-D3/6-311++G for the H H2O 3 cluster with 2600 aD2 tag has no negative anharmonic frequencies, and only two Dunning bs (without D tag) 2 anharmonic intensities which seem unphysically large (modes 22 Pople bs (without D2 tag) Dunning bs (with D tag) 2400 2 and 23). We replace these by the harmonic intensities and present ) Pople bs (with D tag) -1 2 the results in Fig. 6. The frequencies are listed in Table 3.
2200 With VPT2 one is able to assign these modes. We agree with the + assignment [23] of band a12 to frustrated H3O rotation. This mode is weak experimentally, and even weaker in the calculation. The 2000 exact intensity is difficult to ascertain here. For example, in our Wavenumber (cm MP2(fc)/aug-cc-pVTZ calculation (SM Excel file) the intensity is 1800 50 times larger, but the frequency moves down to 618 cm 1. The
previously unassigned band a13 is best attributed to a water wag 1600 overtone, although its intensity is not quite right.
P The most intense low frequency peak, a14, is assigned to the 97 -D3 2X B LY P MP2 OOO as, denoted herein OOO . The OOO is a weak peak predicted Y w as ss M06- 1 W2-P B3L to occur 67 cm to the red. The frequency of the OOOas band (347 Method mP cm 1) is in excellent agreement with experiment (344 cm 1). It is þ also reproducible between different methods, for example in the Fig. 5. Anharmonic PTM frequency of the H ðH2OÞ cluster for different quantum 3 1 methods, bs’s and tagging. above MP2 calculation it is observed at 338 cm , and then OOOss is at 311 cm 1. Experimentally, there is a weak unmarked feature at nearly the same frequency (309 cm 1), which is thus assigned the MP2 calculation in Fig. 3a), these modes and their combination to OOOss. bands can simply be excluded from the analysis. In agreement with Ref. [23], the other high-intensity low fre-
Considering the different computations for the hydrated cluster quency peak (a15) is due to the water wag. From the MP2 calcula- (Excel file of the SM), we find that a quantitative agreement with tion, we can only cite the harmonic frequency, 224 cm 1 , which is experiment is obtained only with the highest level method (MP2) remarkably close to experiment (234 cm 1). In fact, the B3LYP and largest bs (aug-cc-pVTZ) (Fig. 3a). The DFT PES’s appear to be (251 cm 1) and MP2 results bracket this value, as they do for insufficiently accurate for quantitative IR spectroscopy of this clus- OOOas. For the D2-tagged cluster the wag is coupled to D2 transla- ter, even though they have performed well for neutral water clus- tion, either in- or out-of-phase with the two (in-phase) wagging ters [37]. The MP2/aug-cc-pVTZ frequencies are collected in water molecules. Only one of these is observed in the experimen- Table 2, where they are compared with those from experiment tally accessed spectral range. In contrast to the good agreement N.R. Samala, N. Agmon / Chemical Physics 514 (2018) 164–175 169
Table 2 1 þð Þ Calculated frequencies (in cm ) and assignments of the D2-tagged H H2O 3 IR spectrum.
Ref. [22] This work Label Experiment VSCF/VCIa modeb VPT2c AIMDd Assignment
a1 3726 3748 1/2 3726 3707 H2O as a2 3640 3658 3/4 3628 3617 H2O ss a3 3580 5 3540 O3H3 stretch D2 2961 6 3039 D2 stretch a4 2410 2554 7 2383 2450 HB–OHss a5 2307 2364 8 2284 2305 HB–OHas e a6 2109 2240 15 + 10 2084 FR (8) f a7 1943 2034 2 13 1992 FR (7) a8 1878 1950 14 + 13 1883 FR (8) a9 1639 1701 9 1621 1597 bend (ba) a10 1534 1477 12 1553 1527 bend (ba) + a11 1059 1202 13 1080 1157 H3O umbrella (U) + 14 992 1022 H3O rock (R) + 15 618 579 H3O rotation (r)
20 338 352 OOOas MADg 0 72 15.6 33
a Peak frequencies. This calculation does not include the D2 tag, hence we leave out the O3H3 frequency. b Fundamental normal modes for the D2 tagged cluster, numbered from the highest (harmonic) frequency to the lowest (see Excel file in the SM). c MP2/aug-cc-pVTZ, D2-tagged. d B3LYP-D3/aug-TZVP, not D2-tagged, VACF frequencies scaled by 0.966. e Combination band undergoing Fermi resonance with mode number 8. f Overtone undergoing FR with mode number 7. g 1 ; ; Mean absolute deviation from the experimental frequencies (in cm ), averaged over the fundamentals of the bare cluster and the FR (i.e., excluding the a3 D2 a6 and a7 modes).
+ in Refs. [18–20]. The U þ R combination has now inconspicuous B3LYP-D3/6-311++G** (anharmonic) D tagged H (H O) 2 2 3 intensity, marked by the dashed orange arrow in Fig. 7b. Fundamental 1 a 1875
8 1982 In both cases there is at least one HB–OH band near 2300 cm 344 Overtones 234 a Combination bands [16]. Thus the main change with the added D2 tag is in the reas- a 14 15 Expt 2149 signment of a to a combination band, rather than the PTM (both 2100 8 1942 a6 occurring at similar frequencies). Since the origin of an IR band a7 1533 508 (fundamental vs. combination) is not an experimental observable, 1640 a 1060 a 13 668 10 current experiments could not check this, whereas previous com- a11 a 9 2174 a12 putations did not explicitly include the tag. 2232 Intensity 251
347 3.1.5. Indications for Fermi resonances 2020 1075 171
603 FR’s are omnipresent in vibrational spectroscopy [41], for exam- 63 321 127 388 1601 1057 2047 1537 529 ple in CO2 [4], methanol [42] and the Zundel cation [5]. A FR arises 646 when a combination band (or an overtone) becomes nearly degen- 0 500 1000 1500 2000 erate (i.e. of equal frequency) with a fundamental mode having the -1 Wavenumber (cm ) same molecular symmetry. As a result, the combination gains intensity at the expense of the fundamental, and the two bands þð Þ Fig. 6. The experimental IR predissociation spectrum of the D2-tagged H H2O 3 move away from each other along the frequency axis. cluster at low frequencies [23] compared with the (unscaled) VPT2 spectrum In the present case, U is symmetric while R and the PTM are calculated at the B3LYP-D3/6-311++G⁄⁄ level. Black, red and green sticks mark 1 antisymmetric to reflection in a mirror perpendicular to the fundamental, overtone and combination bands, with their frequencies in cm . þ Experimental data (blue line with mode labels) courtesy of Tim Esser and Knut H3’O3H3” plane (Cs point group). Consequently, U R is also Asmis. (For interpretation of the references to colour in this figure legend, the antisymmetric, like the PTM, allowing a FR to occur between them. reader is referred to the web version of this article.) Consider the B3LYP-D3 and MP2/aug-cc-pVTZ calculations that were found to agree best with experiment. In the bare cluster cal- culation, the U þ R combination must be too remote from the PTM, just described, the main band a8 is not accurately predicted by the so that no such FR occurs and the PTM is near 1800 cm 1 (Fig. 5). B3LYP calculation, which we attribute to the FR discussed below. The small perturbation induced by the D2 tag appears to be suffi- cient to bring the PTM into FR with U þ R, pushing the PTM to 3.1.4. The bare cluster the blue and U þ R to the red. This shift is only partially developed Fig. 7 compares the computed spectra of the bare and tagged in the B3LYP-D3/6-311++G⁄⁄ calculation in Fig. 6. clusters using the same MP2 method for both. The bare cluster fre- An indication that this is indeed the case is obtained by plotting quencies are collected in Table 4. According to Ref. [22],D2 binding the intensity of the U þ R combination band as a function of the affects predominantly the O3H3 bond to which it binds. Upon fundamental frequency difference, PTM U R, for all the compu- removing the D2 tag, the O3H3 stretch blue-shifts by approximately tations performed in this study (Fig. 8). For the bare cluster, it is 90 cm 1 (in excellent agreement with 92 cm 1 in VPT2). Here we seen that when this difference (in absolute value) is small, the
find that, in addition, the two HB–OH bands red shift: HB–OHss intensity of the U þ R combination becomes very strong, dropping 1 1 by 90 cm , while HB–OHas (PTM) jumps nearly 500 cm . Conse- abruptly when this frequency gap widens. The latter occurs for the quently, it identifies in the bare cluster with band a8, as suggested B3LYP and MP2 methods that reproduce experiment, and thus we 170 N.R. Samala, N. Agmon / Chemical Physics 514 (2018) 164–175
Table 3 1 þð Þ Measured and calculated band positions in the low frequency IR spectrum (in cm ) of the D2-tagged H H2O 3 cluster with band assignments.
Ref. [23] This work Label Experimenta QCMDb mode # VPT2c AIMDd Assignment
e + a11 1060 1176 13 1075 1157 H3O umbrella (U) + ? 1015 14 1057 1022 H3O rock (R) + a12 668 15 646 579 H3O rotation (r) a13 508 456 2 23 529 wag overtone a14 344 355 19 347 352 OOOas
309? 21 321 322 OOOss e a15 234 241 23 251 228 water wag 22 171 195 water wag MADf 0471543
a Free electron laser spectrum. b Quasiclassical molecular dynamics. This calculation does not include the D2 tag. c DFT based spectrum using the B3LYP functional with the D3 correction and Pople’s 6-311++G⁄⁄ bs (see Excel file in the SM). d B3LYP-D3/aug-TZVP, not D2-tagged, VACF frequencies scaled by 0.966. e Looks more like noise than a real peak, see Fig. 2b in Ref. [23]. f Mean absolute deviation from the experimental frequencies (in cm 1), averaged over all the available frequencies.
a ^ 1 (a) Hydrated isomer (with D2 tag) Fundamental a MP2/aug-cc-pVTZ Overtones 8 1883 a3 Combination bands a 2084 2 Expt 2284 a7
a a 3540 10 6 2383 3726 a9 a11 a5 a D 4 v 2 2693 1992 2367 3628 2544 2664 2175 1553 1907 1080 1621 992 3039
1000 1500 2000 2500 3000 3500 4000 Intensity 2263 3393 1613 2563 1943
2157 ^ 1539 1018 3632 1985 1113 2454 3623 3721 2293 2373
(b) Hydrated isomer (without D2 tag) 1790 -1 Wavenumber (cm )
þð Þ Fig. 7. Comparison of the computed VPT2 IR spectra for (a) the tagged and (b) the bare H H2O 3 cluster (with an upside-down intensity scale) using MP2(fc)/aug-cc-pVTZ [24]. Black, red and green bars indicate the calculated fundamental, overtone and combination bands, respectively. Orange arrows mark the location of the U þ R combination band. In panel (a) the experimental spectrum is reproduced (blue) with the peak labels according to Ref. [22], and the VPT2 fundamental intensities are scaled by a factor of 2. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) assume that it is the correct result. Curiously enough, for the com- These effects are connected by a theory [43] of FR intensity putations with a D2 tag the U þ R intensity is large (as expected ratios (q), as a function of the (hypothetical) frequency split (d0) q when the FR is operative) but nearly independent of the frequency and intensity ratio ( 0) in the absence of the resonance. Assuming gap. q 3 that 0 0 (e.g., for the bare cluster it is 3 10 ), one obtains Eq. What is unusual here is the large intensity of the combination (3) in Ref. [42], namely: band (Ic), which is larger than the intensity of the fundamental d d0 (If ) with which it is in resonance. Defining the intensity ratio as q ¼ : ð1Þ d þ d0 q Ic=If , we find for the D2-tagged cluster that q ¼ 4:05 (see Excel file in SM), which is rather exceptional [43]. Another curious This can be rearranged to read observation is the crossover of the two bands as a function of the 1 q interaction at the O3H3 position. Denoting the frequency differ- d ¼ d; ð2Þ 0 1 þ q ence by d mc mf , we find that it changes sign between the tagged and bare clusters: For the tagged cluster, d ¼ 402 cm 1, from which one immediately sees that d0=d changes sign when d ¼ 1 whereas for the bare cluster 472 cm . q ¼ 1. For the tagged cluster, d ¼ 402 cm 1 and q ¼ 4:05, so that N.R. Samala, N. Agmon / Chemical Physics 514 (2018) 164–175 171