Reinvestigation of the Infrared Spectrum of the Gas-Phase Protonated Water Tetramer

Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods

Vibrational spectroscopy is widely used to gain insights into structural and dynamic properties of chemical, biological, and materials systems. Thus, an efficient and accurate method to simulate vibrational spectra is desired. In this paper, we justify and employ a microcanonical molecular simulation scheme to calculate the vibrational spectra of three challenging water clusters: the neutral water dimer (H4O2), the protonated water trimer (H7O3+), and the protonated water tetramer (H9O4+). We find that with the accurate description of quantum nuclear delocalization effects through the constrained nuclear-electronic orbital framework, including vibrational mode coupling effects through molecular dynamics simulations can additionally improve the vibrational spectrum calculations. In contrast, without the quantum nuclear delocalization picture, conventional ab initio molecular dynamics may even lead to less accurate results than harmonic analysis.

Intra-cluster Charge Migration upon Hydration of Protonated Formic Acid Revealed by Anharmonic Analysis of Cold Ion Vibrational Spectra

The rates of many chemical reactions are accelerated when carried out in micron-sized droplets, but the molecular origin of the rate acceleration remains unclear. One example is the condensation reaction of 1,2-diaminobenzene with formic acid to yield benzimidazole. The observed rate enhancements have been rationalized by invoking enhanced acidity at the surface of methanol solvent droplets with low water content to enable protonation of formic acid to generate a cationic species (protonated formic acid or PFA) formed by attachment of a proton to the neutral acid. Because PFA is the key feature in this reaction mechanism, vibrational spectra of cryogenically cooled, microhydrated PFA·(H2O)n=1-6 were acquired to determine how the extent of charge localization depends on the degree of hydration. Analysis of these highly anharmonic spectra with path integral ab initio molecular dynamics simulations reveals the gradual displacement of the excess proton onto the water network in the microhydration regime at low temperatures with n = 3 as the tipping point for intra-cluster proton transfer.

Proton-transfer spectroscopy beyond the normal-mode scenario

A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single [Formula: see text] cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.

Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions

The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200-300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm-1, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm-1 and agrees in location and width with previous experimental mid-infrared spectra.

Towards complete assignment of the infrared spectrum of the protonated water cluster H+(H2O)21

The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H⁺(H₂O)_n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H⁺(H₂O)₂₁ cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster.

Protonated water species have been the subject of numerous experimental and computational studies. Here the authors provide a nearly complete assignment of the experimental IR spectrum of the H⁺(H₂O)₂₁ water cluster based on high-level wavefunction theory and anharmonic vibrational quasi-degenerate perturbation theory.

Comparative studies of IR spectra of deprotonated serine with classical and thermostated ring polymer molecular dynamics simulations

Here we report the vibrational spectra of deprotonated serine calculated from the classical molecular dynamics (MD) simulations and thermostated ring-polymer molecular dynamics (TRPMD) simulation with third-order density-functional tight-binding. In our earlier study [Inakollu and Yu, "A systematic benchmarking of computational vibrational spectroscopy with DFTB3: Normal mode analysis and fast Fourier transform dipole autocorrelation function," J. Comput. Chem. 39, 2067 (2018)] of deprotonated serine, we observed a significant difference in the vibrational spectra with the classical MD simulations compared to the infrared multiple photon dissociation spectra. It was postulated that this is due to neglecting the nuclear quantum effects (NQEs). In this work, NQEs are considered in spectral calculation using the TRPMD simulations. With the help of potential of mean force calculations, the conformational space of deprotonated serine is analyzed and used to understand the difference in the spectra of classical MD and TRPMD simulations at 298.15 and 100 K. The high-frequency vibrational bands in the spectra are characterized using Fourier transform localized vibrational mode (FT-νN AC) and interatomic distance histograms. At room temperature, the quantum effects are less significant, and the free energy profiles in the classical MD and the TRPMD simulations are very similar. However, the hydrogen bond between the hydroxyl-carboxyl bond is slightly stronger in TRPMD simulations. At 100 K, the quantum effects are more prominent, especially in the 2600-3600 cm-1, and the free energy profile slightly differs between the classical MD and TRPMD simulations. Using the FT-νN AC and the interatomic distance histograms, the high-frequency vibrational bands are discussed in detail.

Molecular Characterization of the Surface Excess Charge Layer in Droplets

The surface excess charge layer (SECL) in droplets has often been associated with distinct chemistry. We examine the effect of the nature of ions in the composition and structure of SECL by using molecular dynamics. We find that in the presence of simple ions the thickness of SECL is invariant not only with respect to droplet size but also with respect to the nature of the ions. In the presence of simple ions, this layer has a thickness of ∼1.5-1.7 nm but in the presence of macroions it may extend to ∼2.0 nm. The proportion of ions contained in SECL depends on the nature of the ions and the droplet size. For the same droplet size, I^- and model H₃O⁺ ions show considerably higher concentration than Na⁺ and Cl^- ions. We identify the maximum ion concentration region, which, in nanodrops, may partially overlap with SECL. As the relative shape fluctuations decrease when microdrop size is approached, the overlap between SECL and maximum ion concentration region increases. We suggest the extension of the bilayer droplet structure assumed in the equilibrium partitioning model of Enke to include the maximum ion concentration region that may not coincide with SECL in nanodrops. We compute the ion concentrations in SECL, which are those that should enter the kinetic equation in the ion-evaporation mechanism, instead of the overall drop ion concentration that has been used.

Direct Synthesis of Dimethyl Ether from Syngas on Bifunctional Hybrid Catalysts Based on Supported H3PW12O40 and Cu-ZnO(Al): Effect of Heteropolyacid Loading on Hybrid Structure and Catalytic Activity

The performance of bifunctional hybrid catalysts based on phosphotungstic acid (H3PW12O40, HPW) supported on TiO2 combined with Cu-ZnO(Al) catalyst in the direct synthesis of dimethyl ether (DME) from syngas has been investigated. We studied the effect of the HPW loading on TiO2 (from 1.4 to 2.7 monolayers) on the dispersion and acid characteristics of the HPW clusters. When the concentration of the heteropoliacid is slightly higher than the monolayer (1.4 monolayers) the acidity of the clusters is perturbed by the surface of titania, while for concentration higher than 1.7 monolayers results in the formation of three-dimensional HPW nanocrystals with acidity similar to the bulk heteropolyacid. Physical hybridization of supported heteropolyacids with the Cu-ZnO(Al) catalyst modifies both the acid characteristics of the supported heteropolyacids and the copper surface area of the Cu-ZnO(Al) catalyst. Hybridization gives rise to a decrease in the copper surface area and the disappearance of the strong acidic sites typical of HPW nanocrystals, showing all hybrids similar acid sites of weak or medium strength. The activity of the hybrids was tested for direct DME synthesis from syngas at 30 bar and 250 °C; only the hybrids with HPW loading higher than 1.4 monolayers showed activity for the direct synthesis of DME, showing that the sample loaded with 2.7 monolayers of heteropolyacid had higher activity than the reference hybrid representative of the most widely applied catalysts based on the combination of Cu-ZnO(Al) with HZSM-5. In spite of the high activity of the hybrids, they show a moderate loss in the DME production with TOS that denotes some kind of deactivation of the acidity function under reaction conditions.

Temperature Dependence of Intramolecular Vibrational Bands in Small Water Clusters

Cyclic water clusters are pivotal for understanding atmospheric reactions as well as liquid water, yet the temperature (T) dependence of their dynamics and spectroscopy is poorly studied. The development of highly accurate water potentials, such as MB-pol, partly rectifies this. It remains to account for the quantum nuclear effects (NQE), because quantum nuclear dynamics become increasingly inaccurate at low temperatures. From a practical point of view, we find that NQE can be accounted for simply by subtracting a constant from the frequencies obtained from the velocity autocorrelation functions (VACF) of classical nuclear dynamics, resulting in unprecedented agreement with experiment, mostly within 5 cm^-1. We have performed classical simulations of (H₂O)_n clusters (n = 2-5) from 20 K and up to their melting temperature, calculating both all-atom and partial VACF, thus generating the temperature dependence of the vibrational frequencies (IR and Raman bands). Focusing on the hydrogen-bonded (HBed) OH stretch and HOH bend, we find opposing T dependencies. The HBed OH modes blue shift linearly with T, attributed to ring expansion rather than any specific conformational change. The lowest-frequency Raman concerted mode is predicted to show the largest such shift. In contrast, the HOH bend undergoes a red-shift, with the highest frequency concerted band undergoing the largest red-shift. These results can be explained by a coupled-oscillator model for n hydrogen atoms on a ring, constrained to move either tangentially (stretch) or perpendicularly (bend) to the ring. With increasing temperature and weakening of HBs, the intrinsic force constant increases (stretch) or remains constant (bend), while the nearest-neighbor coupling constant decreases, and this results in the interesting behavior revealed herein. T-dependent Raman studies are required for testing some of these predictions.

Revisit the landscape of protonated water clusters H+(H2O)n with n = 10-17: An ab initio global search

Using a genetic algorithm incorporated with density functional theory, we explore the ground state structures of protonated water clusters H⁺(H₂O)_n with n = 10-17. Then we re-optimize the isomers at B97-D/aug-cc-pVDZ level of theory. The extra proton connects with a H₂O molecule to form a H₃O⁺ ion in all H⁺(H₂O)_10-17 clusters. The lowest-energy structures adopt a monocage form at n = 10-16 and core-shell structure at n = 17 based on the MP2/aug-cc-pVTZ//B97-D/aug-cc-pVDZ+ZPE single-point-energy calculation. Using second-order vibrational perturbation theory, we further calculate the infrared spectra with anharmonic correction for the ground state structures of H⁺(H₂O)_10-17 clusters at the PBE0/aug-cc-pVDZ level. The anharmonic correction to the spectra is crucial since it reproduces the experimental results quite well. The extra proton weakens the O-H bond strength in the H₃O⁺ ion since the Wiberg bond order of the O-H bond in the H₃O⁺ ion is smaller than that in H₂O molecules, which causes a red shift of the O-H stretching mode in the H₃O⁺ ion.

Deconstructing Prominent Bands in the Terahertz Spectra of H7O3+ and H9O4+: Intermolecular Modes in Eigen Clusters

We report experimental vibrational action spectra (210-2200 cm^-1) and calculated IR spectra, using recent ab initio potential energy and dipole moment surfaces, of H₇O₃⁺ and H₉O₄⁺. We focus on prominent far-IR bands, which postharmonic analyses show, arise from two types of intermolecular motions, i.e., hydrogen bond stretching and terminal water wagging modes, that are similar in both clusters. The good agreement between experiment and theory further validates the accuracy of the potential and dipole moment surfaces, which was used in a recent theoretical study that concluded that infrared photodissociation spectra of the cold clusters correspond to the Eigen isomer. The comparison between theory and experiment adds further confirmation of the need of postharmonic analysis for these clusters.

Disentangling the Complex Vibrational Spectrum of the Protonated Water Trimer, H+(H2O)3, with Two-Color IR-IR Photodissociation of the Bare Ion and Anharmonic VSCF/VCI Theory

Vibrational spectroscopy of the protonated water trimer provides a stringent constraint on the details of the potential energy surface (PES) and vibrational dynamics governing excess proton motion far from equilibrium. Here we report the linear spectrum of the cold, bare H⁺(H₂O)₃ ion using a two-color, IR-IR photofragmentation technique and follow the evolution of the bands with increasing ion trap temperature. The key low-energy features are insensitive to both D₂ tagging and internal energy. The D₂-tagged D⁺(D₂O)₃ spectrum is reported for the first time, and the isotope dependence of the band pattern is surprisingly complex. These spectra are reproduced by large-scale vibrational configuration interaction calculations based on a new full-dimensional PES, which treat the anharmonic effects arising from large amplitude motion. The results indicate such extensive mode mixing in both isotopologues that one should be cautious about assigning even the strongest features to particular motions, especially for the absorptions that occur close to the intramolecular bending mode of the water molecule.

High-Level Quantum Calculations of the IR Spectra of the Eigen, Zundel, and Ring Isomers of H+(H2O)4 Find a Single Match to Experiment

The protonated water tetramer H⁺(H₂O)₄, often written as the Eigen cluster, H₃O⁺(H₂O)₃, plays a central role in studies of the hydrated proton. The cluster has been investigated spectroscopically both experimentally and theoretically with some differences and controversies. The major issue stems from the existence of higher-energy Zundel isomers of this cluster and the role these isomers might play in the IR spectra. Settling this fundamental issue is one goal of this Communication, where high-level quantum calculations of the IR spectra of the Eigen and three isomeric forms of this cluster are presented. These calculations make use of a many-body representation of the potential and dipole moment surfaces and VSCF/VCI calculations of vibrational eigenstates and the IR spectrum. The calculated spectra for the Eigen H₃O⁺(H₂O)₃ and D₃O⁺(D₂O)₃ isomers compare very well with experiment. The calculated spectra for the cis and trans-Zundel and ring isomers show prominent features that do not match with experiment but which can guide future experiments to search for these interesting and important isomers.

Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations

We present our novel methodology for computing Raman optical activity (ROA) spectra of liquid systems from ab initio molecular dynamics (AIMD) simulations. The method is built upon the recent developments to obtain magnetic dipole moments from AIMD and to integrate molecular properties by using radical Voronoi tessellation. These techniques are used to calculate optical activity tensors for large and complex periodic bulk phase systems. Only AIMD simulations are required as input, and no time-consuming perturbation theory is involved. The approach relies only on the total electron density in each time step and can readily be combined with a wide range of electronic structure methods. To the best of our knowledge, these are the first computed ROA spectra for a periodic bulk phase system. As an example, the experimental ROA spectrum of liquid (R)-propylene oxide is reproduced very well.