A theoretical and computational study of the anion, neutral, and cation Cu(H(2)O) complexes

An ab initio investigation of the potential energy surfaces and vibrational energies and wave functions of the anion, neutral, and cation Cu(H(2)O) complexes is presented. The equilibrium geometries and harmonic frequencies of the three charge states of Cu(H(2)O) are calculated at the MP2 level of theory. CCSD(T) calculations predict a vertical electron detachment energy for the anion complex of 1.65 eV and a vertical ionization potential for the neutral complex of 6.27 eV. Potential energy surfaces are calculated for the three charge states of the copper-water complexes. These potential energy surfaces are used in variational calculations of the vibrational wave functions and energies and from these, the dissociation energies D(0) of the anion, neutral, and cation charge states of Cu(H(2)O) are predicted to be 0.39, 0.16, and 1.74 eV, respectively. In addition, the vertical excitation energies, that correspond to the 4 (2)P<--4 (2)S transition of the copper atom, and ionization potentials of the neutral Cu(H(2)O) are calculated over a range of Cu(H(2)O) configurations. In hydrogen-bonded, Cu-HOH configurations, the vertical excitation and ionization energies are blueshifted with respect to the corresponding values for atomic copper, and in Cu-OH(2) configurations where the copper atom is located near the oxygen end of water, both quantities are redshifted.     

HCN(HNC)与NH3, H2O和HF分子间相互作用的理论研究

在MP2/aug-cc-pVTZ水平上, 对HCN(HNC)与NH3, H2O和HF分子间可能存在的氢键型复合物进行了全自由度能量梯度优化, 通过在相同水平上的频率验证分析发现了稳定的分子间相互作用形式是HCN(HNC)作为质子供体或作为质子受体形成的复合物. 基组重叠误差对总相互作用能的影响均小于3.34 kJ/mol. 通过自然键轨道(NBO)分析, 研究了单体和复合物中的原子电荷和电荷转移对分子间相互作用的影响. 对称性匹配微扰理论(SAPT, Symmetry Adapted Perturbation Theory)能量分解结果表明, 在分子间相互作用中, 静电作用与诱导作用占主导地位, 而诱导作用与复合物的电荷转移之间具有良好的正相关性.     

Structures and rearrangement reactions of 4-aminophenol(H2O)1+ and 3-aminophenol(H2O)1+ clusters

In this paper the structures of 4-aminophenol(H2O)1+ and 3-aminophenol(H2O)1+ clusters are investigated in molecular beam experiments by different IR/UV-double resonance techniques as well as the mass analyzed threshold ionization spectroscopy yielding both inter- and intramolecular vibrations of the ionic and neutral species. Possible structures are extensively calculated at the level of density functional theory (DFT) or at the ab initio level of theory. From the experimental and theoretical investigations it can be concluded that in the case of 4-aminophenol(H2O)1 one O-H...O hydrogen-bonded structure exists in the neutral cluster but two structures containing either an O-H...O or a N-H...O hydrogen-bonded arrangement are observed in the spectra of the ionic species. This observation is a result of an intramolecular rearrangement reaction within the ion which can only take place if high excess energies are used. A reaction path via the CH bonds is calculated and explains the experimental observations. In the case of 3-aminophenol(H2O)1+ only one O-H...O bound structure is observed both in the neutral and ionic species. Ab initio and DFT calculations show that due to geometrical and energetical reasons a rearrangement cannot be observed in the 3-aminophenol(H2O)1+ cluster ion.     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.     

Infrared spectrum of NH4+(H2O): evidence for mode specific fragmentation

The gas phase infrared spectrum (3250-3810 cm-1) of the singly hydrated ammonium ion, NH4+(H2O), has been recorded by action spectroscopy of mass selected and isolated ions. The four bands obtained are assigned to N-H stretching modes and to O-H stretching modes. The N-H stretching modes observed are blueshifted with respect to the corresponding modes of the free NH4+ ion, whereas a redshift is observed with respect to the modes of the free NH3 molecule. The O-H stretching modes observed are redshifted when compared to the free H2O molecule. The asymmetric stretching modes give rise to rotationally resolved perpendicular transitions. The K-type equidistant rotational spacings of 11.1(2) cm-1 (NH4+) and 29(3) cm-1 (H2O) deviate systematically from the corresponding values of the free molecules, a fact which is rationalized in terms of a symmetric top analysis. The relative band intensities recorded compare favorably with predictions of high level ab initio calculations, except on the nu3(H2O) band for which the observed value is about 20 times weaker than the calculated one. The nu3(H2O)/nu1(H2O) intensity ratios from other published action spectra in other cationic complexes vary such that the nu3(H2O) intensities become smaller the stronger the complexes are bound. The recorded ratios vary, in particular, among the data collected from action spectra that were recorded with and without rare gas tagging. The calculated anharmonic coupling constants in NH4+(H2O) further suggest that the coupling of the nu3(H2O) and nu1(H2O) modes to other cluster modes indeed varies by orders of magnitude. These findings together render a picture of a mode specific fragmentation dynamic that modulates band intensities in action spectra with respect to absorption spectra. Additional high level electronic structure calculations at the coupled-cluster singles and doubles with a perturbative treatment of triple excitations [CCSD(T)] level of theory with large basis sets allow for the determination of an accurate binding energy and enthalpy of the NH4+(H2O) cluster. The authors' extrapolated values at the CCSD(T) complete basis set limit are De [NH4+-(H2O)]=-85.40(+/-0.24) kJ/mol and DeltaH(298 K) [NH4+-(H2O)]=-78.3(+/-0.3) kJ/mol (CC2), in which double standard deviations are indicated in parentheses.     

PH~3...H~2O体系分子间相互作用的从头计算MO研究

在MP2/6-311++G(3d,3p)水平,对PH~3...H~2O体系可能存在的氢键复合物进行了全自由度能量梯度优化,发现PH~3...H~2O体系存在三个能量极小结构A,B和C。其中结构A和B以H~2O作为质子授体,结构C以PH~3作为质子授体,结构A较结构B和C分别稳定6.52kJ/mol和8.18kJ/mol。结构A具有C~s对称性,其结构中P原子和O原子间距离为354.78nm,键角OHP为171.35ⅲ,属于接近于直线的传统型氢键结构。进一步在高级电子相关校正的MP4SDTQ下,用6-311++G(3df,3pd)基加上键函数{3s3p2d1f},通过BSSE校正,精确计算了结构A的结合能为-10.84kJ/mol。     

Computational study on the characteristics of the interaction in naphthalene...(H2X)n=1,2 (X = O,S) clusters

The characteristics of the interaction between the pi cloud of naphthalene and up to two H2O or H2S molecules were studied. Calculations show that clusters formed by naphthalene and one H2O or H2S molecule have similar geometric features, and also present similar interaction energies. Our best estimates for the interaction energy amount to -2.95 and -2.92 kcal/mol for H2O and H2S, respectively, as obtained with the CCSD(T) method. Calculations at the MP2 level employing large basis sets should be avoided because they produce highly overestimated interaction energies, especially for hydrogen sulfide complexes. The MPWB1K functional, however, provides values pretty similar to those obtained with the CCSD(T) method. Although the magnitude of the interaction is similar with both H2X molecules, its nature is somewhat different: the H2O complex presents electrostatic and dispersion contributions of similar magnitude, whereas for H2S the interaction is dominated by dispersion. In clusters containing two H2X molecules several minima were characterized. In water clusters, the total interaction energy is dominated by the presence of a O-H...O hydrogen bond and, as a consequence, structures where this contact is present are the most stable. However, clusters containing H2S show structures with no interaction between H2S moieties which are as stable as the hydrogen bonded ones, because they allow an optimal H2S...naphthalene interaction, which is stronger than the S-H...S contact. Therefore it is possible that in larger polycycles hydrogen sulfide molecules will be spread onto the surface maximizing S-H...pi interactions rather than aggregated, forming H2S clusters.     

Ab initio potential energy and dipole moment surfaces for H5O2 +

Full-dimensional ab initio potential energy surface (PES) and dipole moment surface (DMS) are reported for H(5)O(2) (+). Tens of thousands of coupled-cluster [CCSD(T)] and second-order Moller-Plesset (MP2) calculations of electronic energies, using aug-cc-pVTZ basis, were done. The energies were fit very precisely in terms of all the internuclear distances, using standard least-square procedures, however, with a fitting basis that satisfies permutational symmetry with respect to like atoms. The H(5)O(2) (+) PES is a fit to 48 189 CCSD(T) energies, containing 7962 polynomial coefficients. The PES has a rms fitting error of 34.9 cm(-1) for the entire data set up to 110 000 cm(-1). This surface can describe various internal floppy motions, including the H atom exchanges, monomer inversions, and monomer torsions. First- and higher-order saddle points have been located on the surface and compared with available previous theoretical work. In addition, the PES dissociates correctly (and symmetrically) to H(2)O+H(3)O(+), with D(e)=11 923.8 cm(-1). Geometrical and vibrational properties of the monomer fragments are presented. The corresponding global DMS fit (MP2 based) involves 3844 polynomial coefficients and also dissociates correctly.     

Structure, stability and vibrational spectrum of the hydrogen-bonded complex between HNO3 and H2O. Ab initio and DFT studies

The structure, stability and vibrational spectrum of the binary complex between HONO2 and H2O have been investigated using ab initio calculations at SCF and MP2 levels with different basis sets and B3LYP/6-31G(d,p) calculations. Full geometry optimization was made for the complex studied. It was established that the hydrogen-bonded H2O...HONO2 complex has a planar structure. The corrected values of the dissociation energy at the SCF and MP2 levels and B3LYP calculations are indicative of relatively strong OH...O hydrogen-bonded interaction. The changes in the vibrational characteristics (vibrational frequencies and infrared intensities) arising from the hydrogen bonding between HONO2 and H2O have been estimated by using the ab initio calculations at SCF and MP2 levels and B3LYP/6-31G(d,p) calculations. It was established that the most sensitive to the complexation is the stretching O-H vibration from HONO2. In agreement with the experiment, its vibrational frequency in the complex is shifted to lower wavenumbers. The predicted frequency shift with the B3LYP/6-31G(d,p) calculations (-439 cm(-1)) is in the best agreement with the experimentally measured (-498 cm(-1)). The intensity of this vibration increases dramatically upon hydrogen bonding. The ab initio calculations at the SCF level predict an increase up to five times; at the MP2 level up to 10 times and the B3LYP/6-31G(d,p) predicted increase is up to 17 times. The good agreement between the predicted values of the frequency shifts and those experimentally observed show that the structure of the hydrogen-bonded complex H2O...HONO2 is reliable.     

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

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

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

Structural investigation of asymmetrical dimer radical cation system (H2O-H2S)+: proton-transferred or hemi-bonded?

Ab initio molecular orbital and hybrid density functional methods have been employed to characterize the structure and bonding of (H2O-H2S)+, an asymmetrical dimer radical cation system. A comparison has been made between the two-center three-electron (2c-3e) hemi-bonded system and the proton-transferred hydrogen-bonded systems of (H2O-H2S)+. Geometry optimization of these systems was carried out using unrestricted Hartree Fock (HF), density functional theory with different functionals, and second-order M?ller-Plesset perturbation (MP2) methods with 6-311++G(d,p) basis set. Hessian calculations have been done at the same level to check the nature of the equilibrium geometry. Energy data were further improved by calculating basis set superposition error for the structures optimized through MP2/6-311++G(d,p) calculations. The calculated results show that the dimer radical cation structure with H2O as proton acceptor is more stable than those structures in which H2O acts as a proton donor or the 2c-3e hemi-bonded (H2O thereforeSH2)+ system. This stability trend has been further confirmed by more accurate G3, G3B3, and CCSD(T) methods. On the basis of the present calculated results, the structure of H4OS+ can best be described as a hydrogen-bonded complex of H3O+ and SH with H2O as a proton acceptor. It is in contrast to the structure of neutral (H2O...H2S) dimer where H2O acts as a proton donor. The present work has been able to resolve the ambiguity in the nature of bonding between H2O and H2S in (H2O-H2S)+ asymmetrical dimer radical cation.     

Ab initio and density functional theory reinvestigation of gas-phase sulfuric acid monohydrate and ammonium hydrogen sulfate

We have calculated the thermochemical parameters for the reactions H(2)SO(4) + H(2)O <--> H(2)SO(4).H(2)O and H(2)SO(4) + NH(3) <--> H(2)SO(4).NH(3) using the B3LYP and PW91 functionals, MP2 perturbation theory and four different basis sets. Different methods and basis sets yield very different results with respect to, for example, the reaction free energies. A large part, but not all, of these differences are caused by basis set superposition error (BSSE), which is on the order of 1-3 kcal mol(-1) for most method/basis set combinations used in previous studies. Complete basis set extrapolation (CBS) calculations using the cc-pV(X+d)Z and aug-cc-pV(X+d)Z basis sets (with X = D, T, Q) at the B3LYP level indicate that if BSSE errors of less than 0.2 kcal mol(-1) are desired in uncorrected calculations, basis sets of at least aug-cc-pV(T+d)Z quality should be used. The use of additional augmented basis functions is also shown to be important, as the BSSE error is significant for the nonaugmented basis sets even at the quadruple-zeta level. The effect of anharmonic corrections to the zero-point energies and thermal contributions to the free energy are shown to be around 0.4 kcal mol(-1) for the H(2)SO(4).H(2)O cluster at 298 K. Single-point CCSD(T) calculations for the H(2)SO(4).H(2)O cluster also indicate that B3LYP and MP2 calculations reproduce the CCSD(T) energies well, whereas the PW91 results are significantly overbinding. However, basis-set limit extrapolations at the CCSD(T) level indicate that the B3LYP binding energies are too low by ca. 1-2 kcal/mol. This probably explains the difference of about 2 kcal mol(-1) for the free energy of the H(2)SO(4) + H(2)O <--> H(2)SO(4).H(2)O reaction between the counterpoise-corrected B3LYP calculations with large basis sets and the diffusion-based experimental values of S. M. Ball, D. R. Hanson, F. L Eisele and P. H. McMurry (J. Phys. Chem. A. 2000, 104, 1715). Topological analysis of the electronic charge density based on the quantum theory of atoms in molecules (QTAIM) shows that different method/basis set combinations lead to qualitatively different bonding patterns for the H(2)SO(4).NH(3) cluster. Using QTAIM analysis, we have also defined a proton transfer degree parameter which may be useful in further studies.     

Mechanism of OH radical reactions with HCN and CH3CN: OH regeneration in the presence of O2

A theoretical study on the mechanism of the OH reactions with HCN and CH(3)CN, in the presence of O2, is presented. Optimum geometries and frequencies have been computed at BHandHLYP/6-311++G(2d,2p) level of theory for all stationary points. Energy values have been improved by single-point calculations at the above geometries using CCSD(T)/6-311++G(2d,2p). The initial attack of OH to HCN was found to lead only to the formation of the HC(OH)N adduct, while for CH(3)CN similar proportions of CH(2)CN and CH(3)C(OH)N are expected. A four-step mechanism has been proposed to explain the OH regeneration, experimentally observed for OH + CH(3)CN reaction, when carried out in the presence of O2. The mechanism steps are as follows: (1) OH addition to the C atom in the CN group, (2) O2 addition to the N atom, (3) an intramolecular H migration from OH to OO, and (4) OH elimination. This mechanism is in line with the one independently proposed by Wine et al. for HCN. The results obtained here suggest that for the OH + HCN reaction, the OH regeneration might occur even in larger extension than for OH + CH(3)CN reaction. The agreement between the calculated data and the available experimental evidence on the studied reactions seems to validate the mechanism proposed here.     

Theoretical study of photoacidity of HCN: the effect of complexation with water

The character of the hydrogen bonding and the excited state proton transfer (ESPT) in the model system HCN...H(2)O is investigated. The PES of the two lowest excited states of the H(2)O...HCN complex was calculated using the CASPT2 method. The nonadiabatic coupling of the two states of the (pi-->pi*) and (pi-->sigma*) character is responsible for the excited state proton/hydrogen transfer. Compared to the ground state, the barrier for this process is significantly smaller. An increased number of water molecules in the complex with cyclic hydrogen-bonded network causes a large blue shift of the state of the (pi-->sigma*) character. The question of the dissociation of the complex in its excited state is also addressed.     

Model systems for probing metal cation hydration: the V+(H2O) and ArV+(H2O) complexes

In support of mass-selected infrared photodissociation (IRPD) spectroscopy experiments, coupled-cluster methods including all single and double excitations (CCSD) and a perturbative contribution from connected triple excitations [CCSD(T)] have been used to study the V+(H2O) and ArV+(H2O) complexes. Equilibrium geometries, harmonic vibrational frequencies, and dissociation energies were computed for the four lowest-lying quintet states (5A1, 5A2, 5B1, and 5B2), all of which appear within a 6 kcal mol(-1) energy range. Moreover, anharmonic vibrational analyses with complete quartic force fields were executed for the 5A1 states of V+(H2O) and ArV+(H2O). Two different basis sets were used: a Wachters+f V[8s6p4d1f] basis with triple-zeta plus polarization (TZP) for O, H, and Ar; and an Ahlrichs QZVPP V[11s6p5d3f2g] and Ar[9s6p4d2f1g] basis with aug-cc-pVQZ for O and H. The ground state is predicted to be 5A1 for V+(H2O), but argon tagging changes the lowest-lying state to 5B1 for ArV+(H2O). Our computations show an opening of 2 degrees -3 degrees in the equilibrium bond angle of H2O due to its interaction with the metal ion. Zero-point vibrational averaging increases the effective bond angle further by 2.0 degrees -2.5 degrees, mostly because of off-axis motion of the heavy vanadium atom rather than changes in the water bending potential. The total theoretical shift in the bond angle of about +4 degrees is significantly less than the widening near 9 degrees deduced from IRPD experiments. The binding energies (D0) for the successive addition of H2O and Ar to the vanadium cation are 36.2 and 9.4 kcal mol(-1), respectively.     

低温基质隔离红外光谱和从头算研究CF~2O...IF

在10KAr和Xe基质中测得CF~2O...IF的两种分子间配合物A和B的红外光谱,观察到A和B之间的异构化反应。采用量子化学从头计算,在MP2/LANL2DZ水平上计算得到A和B的平衡构型的结构参数和振动频率。IF与CF~2O的相互作用引起CF~2O的C=O双键伸长、C-F键缩短,导致v~C~=~O的红移和v~C~F~2的蓝移。由势能曲线可看出,A的势能比B高23kJ/mol,表明B比A稳定,升温能够发生A向B的转化。     

Ionized state of hydroperoxy radical-water hydrogen-bonded complex: (HO2-H2O)+

Ab initio molecular orbital calculations have been employed to characterize the structure and bonding of the (HO2-H2O)+ radical cation system. Geometry optimization of this system was carried out using unrestricted density functional theory in conjunction with the BHHLYP functional and 6-311++G(2df,2p) as well as 6-311++G(3df,3p) basis sets, the second-order M?ller-Plesset perturbation (MP2) method with the 6-311++G(3df,3p) basis set, and the couple cluster (CCSD) method with the aug-cc-pVTZ basis set. The effect of spin multiplicity on the stability of the (HO2-H2O)+ system has been studied and also compared with that of oxygen. The calculated results suggest a proton-transferred hydrogen bond between HO2 and H2O in H3O3+ wherein a proton is partially transferred to H2O producing the O2...H3O+ structure. The basis set superposition error and zero-point energy corrected results indicate that the H3O3+ system is energetically more stable in the triplet state; however, the singlet state of H3O3+ is more stable with respect to its dissociation into H3O+ and singlet O2. Since the resulting proton-transferred hydrogen-bonded complex (O2...H3O+) consists of weakly bound molecular oxygen, it might have important implications in various chemical processes and aquatic life systems.     

Ab initio and DFT studies of the vibrational spectra of hydrogen-bonded PhOH...(H2O)4 complexes

The vibrational characteristics (vibrational frequencies and infrared intensities) for the hydrogen-bonded complex of phenol with four water molecules PhOH...(H2O)4 (structure 4A) have been predicted using ab initio and DFT (B3LYP) calculations with 6-31G(d,p) basis set. The changes in the vibrational characteristics from free monomers to a complex have been calculated. The ab initio and B3LYP calculations show that the observed four intense bands at 3299, 3341, 3386 and 3430 cm(-1) can be assigned to the hydrogen-bonded OH stretching vibrations in the complex PhOH...(H2O)4 (4A). The complexation leads to very large red shifts of these vibrations and very strong increase in their IR intensity. The predicted red shifts for these vibrations with B3LYP/6-31G(d,p) calculations are in very good agreement with the experimentally observed. It was established that the phenolic OH stretching vibration is the most sensitive to the hydrogen bonding. The predicted red-shift with the B3LYP/6-31G(d,p) calculations for the most stable ring structure 4A (-590 cm(-1)) is in better agreement with the experimentally observed than the red-shift, predicted with SCF/6-31G(d,p) calculations. The magnitude of the wavenumber shift is indicative of relatively strong OH...H hydrogen-bonded interaction. The complexation between phenol and four water molecules leads to strong increase of the IR intensity of the phenolic OH stretching vibration (up to 38 times).