Analysis of nuclear quantum effects on hydrogen bonding

The impact of nuclear quantum effects on hydrogen bonding is investigated for a series of hydrogen fluoride (HF)n clusters and a partially solvated fluoride anion, F-(H2O). The nuclear quantum effects are included using the path integral formalism in conjunction with the Car-Parrinello molecular dynamics (PICPMD) method and using the second-order vibrational perturbation theory (VPT2) approach. For the HF clusters, a directional change in the impact of nuclear quantum effects on the hydrogen-bonding strength is observed as the clusters evolve toward the condensed phase. Specifically, the inclusion of nuclear quantum effects increases the F-F distances for the (HF)n=2-4 clusters and decreases the F-F distances for the (HF)n>4 clusters. This directional change occurs because the enhanced electrostatic interactions between the HF monomers become more dominant than the zero point energy effects of librational modes as the size of the HF clusters increases. For the F-(H2O) system, the inclusion of nuclear quantum effects decreases the F-O distance and strengthens the hydrogen bonding interaction between the fluoride anion and the water molecule because of enhanced electrostatic interactions. The vibrationally averaged 19F shielding constant for F-(H2O) is significantly lower than the value for the equilibrium geometry, indicating that the electronic density on the fluorine decreases as a result of the quantum delocalization of the shared hydrogen. Deuteration of this system leads to an increase in the vibrationally averaged F-O distance and nuclear magnetic shielding constant because of the smaller degree of quantum delocalization for deuterium.     

Theoretical analysis of pyridine protonation in water clusters of increasing size.

The protonation of pyridine in water clusters as a function of the number of water molecules was theoretically analyzed as a prototypical case for the protonation of organic bases. We determined the variation of structural, bonding, and energetic properties on protonation, as well as the stabilization of the ionic species formed. Thus, we used supermolecular models in which pyridine interacts with clusters of up to five water molecules. For each complex, we determined the most stable unprotonated and protonated structures from a simulated annealing at the semi ab initio level. The structures were optimized at the B3LYP/cc-pVDZ level. We found that the hydroxyl group formed on protonation of pyridine abstracts a proton from the ortho-carbon atom of the pyridine ring. The "atoms in molecules" theory showed that this C-H group loses its covalent character. However, starting with clusters of four water molecules, the C-H bond recovers its covalent nature. This effect is associated with the presence of more than one ring between the water molecules and pyridine. These rings stabilize, by delocalization, the negative charge on the hydroxyl oxygen atom. Considering the protonation energy, we find that the protonated forms are increasingly stabilized with increasing size of the water cluster. When zero-point energy is included, the variation follows closely an exponential decrease with increasing number of water molecules. Analysis of the vibrational modes for the strongest bands in the IR spectra of the complexes suggests that the protonation of pyridine occurs by concerted proton transfers among the different water rings in the structure. Symmetric water stretching was found to be responsible for hydrogen transfer from the water molecule to the pyridine nitrogen atom.     

Nuclear quantum effects and hydrogen bonding in liquids

We have employed ab initio path integral molecular dynamics simulations to investigate the role of nuclear quantum effects on the strength of hydrogen bonds in liquid hydrogen fluoride. Nuclear quantum effects are shown to be responsible for a stronger hydrogen bond and an enhanced dipole-dipole interaction, which lead, in turn, to a shortening of the H...F intrachain distance. The simulation results are analyzed in terms of the electronic density shifts with respect to a purely classical treatment of the nuclei. The observed enhanced hydrogen-bond interaction, which arises from a coupling of intra- and intermolecular effects, should be a general phenomenon occurring in all hydrogen-bonded systems.     

CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule

The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed.     

Nuclear quadrupole hyperfine structure in HC14N/H14NC and DC15N/D15NC isomerization: a diagnostic tool for characterizing vibrational localization

Large-amplitude molecular motions which occur during isomerization can cause significant changes in electronic structure. These variations in electronic properties can be used to identify vibrationally-excited eigenstates which are localized along the potential energy surface. This work demonstrates that nuclear quadrupole hyperfine interactions can be used as a diagnostic marker of progress along the isomerization path in both the HC14N/H14NC and DC15N/D15NC chemical systems. Ab initio calculations at the CCSD(T)/cc-pCVQZ level indicate that the hyperfine interaction is extremely sensitive to the chemical bonding of the quadrupolar 14N nucleus and can therefore be used to determine in which potential well the vibrational wavefunction is localized. A natural bonding orbital analysis along the isomerization path further demonstrates that hyperfine interactions arise from the asphericity of the electron density at the quadrupolar nucleus. Using the CCSD(T) potential surface, the quadrupole coupling constants of highly-excited vibrational states are computed from a one-dimensional internal coordinate path Hamiltonian. The excellent agreement between ab initio calculations and recent measurements demonstrates that nuclear quadrupole hyperfine structure can be used as a diagnostic tool for characterizing localized HCN and HNC vibrational states.     

Determination of the absolute configuration of chiral alpha-aryloxypropanoic acids using vibrational circular dichroism studies: 2-(2-chlorophenoxy) propanoic acid and 2-(3-chlorophenoxy) propanoic acid

The enantiomers of 2-(2-chlorophenoxy) propanoic acid and 2-(3-chlorophenoxy) propanoic acid were resolved on a chiral HPLC column and investigated using mid-infrared vibrational circular dichroism (VCD). Experimental infrared vibrational absorption and VCD spectra were measured in CDCl3 solution in the 2000-900 cm-1 region and compared with the ab initio predictions of absorption and VCD spectra. The predicted spectra were obtained with density functional theory using B3LYP/6-31G* basis set for the stable and dominant conformers. But the predicted spectra did not provide unambiguous structural information due to intermolecular hydrogen bonding in solution. To eliminate the hydrogen bonding effects, the acids were converted to the corresponding methyl esters and the experimental absorbance and VCD spectra of methyl esters were measured. B3LYP predicted spectra were also obtained for the stable and dominant conformers of the esters. From a comparison of the experimental VCD spectra of methyl esters with corresponding ab initio predictions, the absolute configurations of esters, and therefore of their parent acids, are unambiguously determined to be (+)-(R).     

Near IR overtone spectral investigations of cyclohexanol using local mode model--evidence for variation of anharmonicity with concentration due to hydrogen bonding

The near infrared vibrational overtone absorption spectrum of liquid phase cyclohexanol in carbon tetrachloride in different concentrations are examined in the region Deltav=2, 3 and 4. The free and bonded OH local mode mechanical frequency values and anharmonicity values obtained from fitting the overtones are analysed. The observation supports the conclusions drawn from earlier experimental studies on anharmonicity variation of OH-stretching vibrations of alcohols due to intermolecular hydrogen bonding. Our observation is also in agreement with the ab initio calculations on water dimer and trimer. Mechanical anharmonicity of bonded OH-stretching bands tends to increase as a consequence of strong hydrogen bonding at higher concentrations.     

NIR-FT Raman and infrared spectra and ab initio computations of glycinium oxalate

The single crystals of glycinium oxalate are grown by slow evaporation technique and vibrational spectral analysis is carried out using NIR-FT Raman and FT-IR spectra. The ab initio quantum computations are also performed at HF/6-31 G(d) level to derive the optimized geometry, atomic charges and vibrational frequencies of the glycinium oxalate molecule. Vibrational analysis indicates the presence of peculiar intermolecular C-H...O hydrogen bonding interaction producing "blue shift" of C-H stretching frequency. The vibrational spectra confirm the existence of NH3(+) in glycinium oxalate. Hydroxyl vibrations with different inter and intra molecular H-bonding are analysed, supported by computed results.     

Revealing Intra- and Intermolecular Interactions Determining Physico-Chemical Features of Selected Quinolone Carboxylic Acid Derivatives

The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car–Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion—the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.     

Ab initio prediction of the vibrational spectrum of the hydrogen bonded complex O2N...HONO2

The changes in the structural parameters and vibrational characteristics (vibrational frequencies, infrared intensities and Raman activities) arising from the hydrogen bonding between NO(2) and HONO(2) have been studied employing ab initio 6-31G(d, p)/UHF and 6-31+G(d, p)/UHF, and B3LYP/6-31G(d, p) calculations. The charge rearrangement upon hydrogen bonding have been, estimated using the Mulliken population analyses. It was established that the complexation between NO(2) and HONO(2) leads to changes in the structural parameters and the vibrational characteristics of the monomers. The most sensitive to the hydrogen bond formation are the vibrational characteristics of the normal modes of the monomer bonds participating in the hydrogen bonding. The predicted shifts in the vibrational frequencies by ab initio and B3LYP/6-31G(d, p) calculations are in very good agreement with the experimentally observed, which is an evidence for the reliance of the studied structure.     

Theoretical study of the intermolecular hydrogen bond interaction for furan-HCl and furan-CHCl_3 complexes

The importance of intermolecular interactions in biology and material science has prompted chemists to explore the nature of the variety of such interactions. The strongest of these interac-tions are the hydrogen bonds, which play an important role in determining the molecular confor-mation, crystal packing, and the structure of biological systems such as nucleic acids. Extensive experimental and theoretical efforts[1—5] have been devoted to the studies of this type of interac-tions, such as …     

Modeling competitive interactions in proteins: vibrational spectroscopy of M+(n-methylacetamide)1(H2O)n=0-3, M=Na and K, in the 3 microm region

To properly understand the preferred structures and biological properties of proteins, it is important to understand how they are influenced by their immediate environment. Competitive intrapeptide, peptide...water, ion...water, and ion...peptide interactions, such as hydrogen bonding, play a key role in determining the structures, properties, and functionality of proteins. The primary types of hydrogen bonding involving proteins are intramolecular amide...amide (N-H...O=C) and intermolecular amide...water (O-H...O=C and H-O...H-N). n-Methylacetamide (NMA) is a convenient model for investigating these competitive interactions. An analysis of the IR photodissociation (IRPD) spectra of M+(n-methylacetamide)1(H2O)n=0-3 (M=Na and K) in the O-H and N-H spectral regions is presented. Ab initio calculations (MP2/cc-pVDZ) are used as a guide in identifying both the type and location of hydrogen bonds present. In larger clusters, where several structural isomers may be present in the molecular beam, ab initio calculations are also used to suggest assignments for the observed spectral features. The results presented offer insight to the nature of ion...NMA interactions in an aqueous environment and reveal how different ion...ligand pairwise interactions direct the extent of water...water and water...NMA hydrogen bonding observed.     

On the hydrogen bonding between water and ammonia molecules

The hydrogen bond N·HO between the water and ammonia molecules has been investigated ab initio using the SCF LCAO MO method. The minimal and extended basis sets of Slater type orbitals were used. It was found that the energy of the hydrogen bond is equal to 6.44 kcal/mole and the equilibrium separation of the oxygen and nitrogen atoms in the dimer is 5.72 au. At this intermolecular distance there is only one minimum in the potential energy curve for the motion of proton.     

Triple hydrogen bonding in a circular arrangement: ab initio, DFT and first-principles MD studies of tris-hydroxyaryl enamines

First-principles Car-Parrinello molecular dynamics, ab initio (MP2) and density functional schemes have been used to explore the tautomeric equilibrium in three tris(amino(R)methylidene)cyclohexane-1,3,5-triones (R?=?hydrogen, methyl or phenyl group). The dynamic nature of the cyclic hydrogen bonding has been studied by the first-principles MD method. The comparison of the results obtained by aforesaid methods has been accomplished on the basis of calculations of structural and spectroscopic characteristics of the compounds. The conformational analysis of the studied compounds has been carried out at the MP2/6-31+G(d,p) and B3LYP/6-31+G(d,p) levels of theory. The influence of steric and electronic effects on the cyclic hydrogen bonding has been analysed. The extent of the proton delocalization has been modified by the substituents according to the sequence: hydrogen?<?phenyl?<?methyl. This fact is verified by the spectroscopic and structural data as well as the energy potential curve. A prevalence of the keto-enamine tautomeric form has been observed in the static ab initio and DFT models, and confirmed by the first-principles MD.     

Spectroscopic and quantum chemical study of the novel compound cyclopropylmethylselenol

An investigation into the properties of the novel compound cyclopropylmethylselenol has been undertaken by use of Stark-modulation microwave spectroscopy and high-level quantum chemical calculations. Ground-state spectra belonging to six isotopomers of a single conformer of the molecule were recorded and assigned. This conformer, predicted to be the lowest in energy by a series of quantum chemical calculations, possesses a synclinal arrangement of the H-C-C-Se atoms. In addition to the assignment of these ground-state spectra, transitions attributable to vibrationally excited states of the 78Se- and 80Se-containing isotopomers were identified. A tentative assignment of these excited-state spectra to specific vibrational modes has been made with the assistance of a density functional theory calculation at the B3LYP/6-311++G(3df,2pd) level of theory. Close agreement was found between experimental ground-state rotational constants and ab initio equilibrium values calculated at the MP2/aug-cc-pVTZ level of theory. Good agreement was also noted between certain r(s) principal axis coordinates of atoms in the molecule and the corresponding ab initio r(e) values. Limited evidence in favor of the formation of a weak intramolecular hydrogen bond between the H atom of the selenol group and electron density associated with the cyclopropyl ring is discussed.     

Infrared-visible and visible-visible double resonance spectroscopy of 1-hydroxy-9,10-anthraquinone-(H2O)n (n=1,2) complexes

The structures of hydrated 1-hydroxyanthraquinone complexes (1-HAQ), 1-HAQ(H2O)n=1,2, with intramolecular and intermolecular hydrogen bonding interactions were studied using laser spectroscopic methods such as laser induced fluorescence, fluorescence-detected infrared, infrared-visible hole burning, and visible-visible hole burning spectroscopy. In the 1:1 complex 1-HAQ(H2O)1, the water binds to the free carbonyl group of 1-HAQ not associated with intramolecular hydrogen bond. The second water in the 1:2 complex, 1-HAQ(H2O)2, binds to the first water of the 1:1 complex rather than other hydrogen bonding sites of 1-HAQ. A pair of two geometric isomers was produced in a supersonic jet for each of the 1:1 and 1:2 complexes. Both isomers of each complex have the same vibrational spectra in the region of the OH stretching vibration of water, but have different energies for the 0-0 band of vibronic transition due to the asymmetry of the two phenyl rings in 1-HAQ. The 0-0 bands for all four species of 1-HAQ(H2O)n=1,2 were unambiguously assigned by comparing with the results of ab initio calculations, which yielded the structures, vibrational frequencies, and relative energies of the frontier molecular orbitals.     

Intra- and intermolecular interactions in the crystals of 3,4-diamino-1,2,4-triazole and its 5-methyl derivative. Experimental and theoretical investigations of charge density distribution

In this paper intra- and intermolecular interactions in crystals of 3,4-diamino-1,2,4-triazole (DAT) and its 5-methyl derivative (DAMT) were investigated in details by experimental (high-resolution X-ray diffraction) and theoretical (ab initio quantum chemistry (MP2/aug-cc-pvdz), AIM, and NBO) methods. Influence of n-π conjugation and n→σ* hyperconjugation on the geometry of DAT and DAMT molecules was shown. All intermolecular interactions in crystals of the DAT and DAMT including weak X-H···π and mixed X-H···N/X-H···π hydrogen bonds were considered. Comparison of BCP characteristics of these interactions from experimental and theoretical charge density distribution demonstrates systematic increase of bonding in isolated dimers compared with dimers in the crystal phase. The ability of amino groups in both crystals serve as proton acceptors in hydrogen bonding was confirmed.     

Hydrogen-bond assisted enormous broadening of infrared spectra of phenol-water cationic cluster: an ab initio mixed quantum-classical study

The infrared spectrum of phenol-water cationic cluster, [PhOH.H2O]+, taken by Sawamura et al. [J. Phys. Chem. 100, 8131 (1996)] is puzzling in that the peak due to the stretching mode of the phenolic OH (3657 cm-1 for a neutral monomer and 3524 cm-1 for PhOH.H2O) seemingly disappears and instead an extremely broad tail extending down to 2900 cm-1 is observed. The present authors theoretically ascribe this anomalous spectrum to an inhomogeneous broadening of the OH stretching peak caused by the hydrogen bond, the strength of which has been greatly enhanced by ionization of the phenyl ring. Indeed they estimate that the peak position is at 2300 cm-1 and the spectral width can become as wide as 1000 cm-1 at the cluster energy of 32 kcal/mol. This surprisingly wide broadening can be generic in hydrogen-bond systems, which in turn is useful to study the nature of the hydrogen-bond assisted dynamics in various systems such as those in DNA and proteins. To study the present system quantitatively, the authors have developed an ab initio mixed quantum-classical method, in which the nuclear motions on an adiabatic ab initio potential surface are treated such that only the OH stretching motion is described quantum mechanically, while all the other remaining modes are treated classically with on-the-fly scheme. This method includes the implementation of many numerical methodologies, which enables it to deal with a relatively large molecular system. With this theoretical method, the authors analyze the present anomalous broadening in a great detail. In particular, they suggest that one can extract direct information about the hydrogen-bond dynamics with respect to the clear correlation between the vibrational excitation energy of the OH stretching and intermolecular distance by means of a time-resolved infrared spectroscopy: Reflecting the slow and wide-range variation of the intermolecular distance of the relevant hydrogen bond, the time-resolved spectrum is predicted to vary (shift) largely covering the wide range of frequency domain. Thus, it is found that the short-time average along a selected trajectory sensitively reflects the change of the intermolecular distance. The authors also study the effect of internal energy on the hydrogen bonding and the OH spectrum.     

Analysis of the nuclear-electronic orbital method for model hydrogen transfer systems

Fundamental issues associated with the application of the nuclear-electronic orbital (NEO) approach to hydrogen transfer systems are addressed. In the NEO approach, specified nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wavefunctions are calculated with molecular orbital methods. The positions of the nuclear basis function centers are optimized variationally. In the application of the NEO approach to hydrogen transfer systems, the hydrogen nuclei and all electrons are treated quantum mechanically. Within the NEO framework, the transferring hydrogen atom can be represented by two basis function centers to allow delocalization of the proton vibrational wavefunction. In this paper, the NEO approach is applied to the [He-H-He]+ and [He-H-He]++ model systems. Analyses of technical issues pertaining to flexibility of the basis set to describe both single and double well proton potential energy surfaces, linear dependency of the hydrogen basis functions, multiple minima in the basis function center optimization, convergence of the number of hydrogen basis function centers, and basis set superposition error are presented. The accuracy of the NEO approach is tested by comparison to grid calculations for these model systems.