Experimental matrix isolation study and quantum-mechanics-based normal-coordinate analysis of the anharmonic infrared spectrum of picolinic acid N-oxide

This work is, according to our knowledge, the first experimental matrix isolation study of a molecular system with a very short and strong intramolecular OH...O hydrogen bond. It also includes a satisfying interpretation of its entire infrared spectrum. The interpretation relies on the calculation at the DFT/B3LYP/6-31G(d,p) level of the harmonic spectrum and of the anharmonic relaxed potential energy for the stretching motion of the hydrogen-bonded proton, used with our recently modified quantum-mechanics-based normal-coordinate analysis. An important observation about the anharmonic spectrum obtained from this procedure is that the OH stretch coordinate contributes to several normal modes, mixing extensively with other in-plane internal coordinates, in particular OH-bending and C=O-stretching. The two intense normal modes with the largest contributions from the OH-stretching coordinate to the potential energy distribution and to the intensity are located near 1700 and 1500 cm(-1). A calculated anharmonic spectrum obtained from this procedure agrees with the experimental spectrum (frequencies and intensity distribution), within the limits of the estimated uncertainties for the calculation and experiment, allowing the interpretation of the latter. The agreement for the frequencies is about 1-3%. The anharmonic spectrum calculated using the anharmonic keyword in Gaussian 03w is not in satisfactory agreement with experiment insofar as the OH-stretching mode is concerned.     

Properties of strong hydrogen-bonded systems. II. Ab initioSCF-MO study of the hydrogen bond between nitric acid and ammonia

The hydrogen-bonded complex between nitric acid and ammonia molecules has been studied by the ab initio molecular orbital method using the 4-31G basis set. The calculated interaction energy for the complex (ΔE = ?91.4 kJ mole^?1) indicates that one is dealing with the strongest “nonionic” H-bonded complex considered hitherto by theoretical methods. Other properties of the hydrogen-bonded complex such as geometrical parameters, dipole moment, amount of charge transfer, and stretching force constants of the O? H and (OH)… N bonds are calculated and discussed.     

Infrared spectroscopy of the intramolecular hydrogen bond in acethylacetone: a computational approach

The intramolecular hydrogen bond in the enol-acethylacetone (ACAC) is investigated by performing reduced-dimensional quantum calculations. To analyze the shared proton vibrations, two sets of coordinates were employed: normal mode coordinates describing the motion in the vicinity of the most stable configuration, and internal coordinates accounting for the double minimum proton motion. It is proved that the extreme broadness of the OH-stretch band in ACAC is a consequence of the coexistence of two enol-ACAC structures: the global minimum and the transition state for rotation of the distal methyl group. Further, a ground-state tunneling splitting of 116 cm(-1) is found, and it is shown that the inclusion of the kinematic coupling is mandatory when treating large-amplitude proton motion. In the OH-stretch direction a splitting of 853 cm(-1) was predicted.     

Car-Parrinello and path integral molecular dynamics study of the hydrogen bond in the chloroacetic acid dimer system

We have studied the double proton transfer (DPT) reaction in the cyclic dimer of chloroacetic acid using both classical and path integral Car-Parrinello molecular dynamics. We also attempt to quantify the errors in the potential energy surface that arise from the use of a pure density functional. In the classical dynamics a clear reaction mechanism can be identified, where asynchronized DPT arises due to coupling between the O-H stretching oscillator and several low energy intermolecular vibrational modes. This mechanism is considerably altered when quantum tunneling is permitted in the simulation. The introduction of path integrals leads to considerable changes in the thermally averaged molecular geometry, leading to shorter and more centered hydrogen bond linkages.     

Experimental and theoretical study of the polarized infrared spectra of the hydrogen bond in 3‐thiophenic acid crystal

This article presents the results of experimental and theoretical studies of the v_{O H} and v_{O D} band shapes in the polarized infrared spectra of 3‐thiophenic acid crystals measured at room temperature and at 77 K. The line shapes are studied theoretically within the framework of the anharmonic coupling theory, Davydov coupling, Fermi resonance, direct and indirect damping, as well as the selection rule breaking mechanism for forbidden transitions. The adiabatic approximation allowing to separate the high‐frequency motion from the slow one of the H‐bond bridge is performed for each separate H‐bond bridge of the dimer and a strong nonadiabatic correction is introduced via the resonant exchange between the fast‐mode excited states of the two moieties. The spectral density is obtained within the linear response theory by Fourier transform of the damped autocorrelation functions. The approach correctly fits the experimental line shape of the hydrogenated compound and predicts satisfactorily the evolution in the line shapes with temperature and the change in the line shape with isotopic substitution. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Experimental and computational study of crystalline formic acid composed of the higher-energy conformer

Crystalline formic acid (FA) is studied experimentally and by first-principles simulations in order to identify a bulk solid structure composed of the higher-energy (cis) conformer. In the experiments, deuterated FA (HCOOD) was deposited in a Ne matrix and transformed to the cis conformer by vibrational excitation of the ground state (trans) form. Evaporation of the Ne host above 13 K prepared FA in a bulk solid state mainly composed of cis-FA. Infrared absorption spectroscopy at 4.3 K shows that the obtained solid differs from that composed of trans-FA molecules and that the state persists up to the annealing temperature of at least 110 K. The first-principles simulations reveal various energetically stable periodic chain structures containing cis-FA conformers. These chain structures contain either purely cis or both cis and trans forms. The vibrational frequencies of the calculated structures were compared to the experiment and a tentative assignment is given for a novel solid composed of cis-FA.     

Theoretical study of the polarized infrared spectra of the hydrogen bond in 2‐furoic acid crystal dimer

This work presents a theoretical simulation of ν_{O? H} and ν_{O? D} band shapes in the polarized infrared spectra of 2‐furoic acid dimer crystals measured at liquid‐nitrogen temperature. The line shapes are studied theoretically within the framework of the anharmonic couplings between low‐frequency hydrogen‐bond vibrations and degenerate excited states of high‐frequency hydrogen vibrations in hydrogen‐bonded dimers and the anharmonic coupling between the first excited state of the fast mode and the harmonics or band combinations of some low‐frequency bending modes, which lead to Fermi resonances.This approach takes into account the adiabatic approximation, the intrinsic anharmonicity of the low‐frequency mode through a Morse potential, Davydov coupling triggered by resonance exchange between the excited states of the fast modes of the two hydrogen bonds involved in the cyclic dimer, and the direct and indirect damping of the fast‐stretching modes of the hydrogen bonds and of the bending modes. The infrared spectral density was calculated within the linear response theory by Fourier transform of the autocorrelation function of the transition dipole moment operator of the fast mode. Numerical results show that mixing of all these effects allows satisfactory reproduction of the main features of the experimental IR line shapes of crystalline H‐ and D‐bonded 2‐furoic acid at liquid‐nitrogen temperature and for different polarizations. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

On the gas phase hydrogen bond complexes between formic acid and hydroperoxyl radical. A theoretical study

We present a systematic study on the gas-phase hydrogen-bonded complexes formed between formic acid and hydroperoxyl radical, which has been carried out by using B3LYP and CCSD(T) theoretical approaches in connection with the 6-311+G(2df,2p) basis set. For all complexes we have employed the AIM theory by Bader and the NBO partition scheme by Weinhold to analyze the bonding features. We have found 17 stationary points, and 11 of them present a cyclic structure. Their computed stabilities vary from 0.3 to 11.3 kcal/mol, depending on several factors, such as involvement in the hydrogen bond interaction, the geometrical constraints, and the possible concurrence of further effects such as resonance-assisted hydrogen bonds or inductive effects. In addition, three stationary points correspond to transition structures involving a double proton-transfer process whose features are also analyzed.     

Effects of strong and weak hydrogen bond formation on VCD spectra: a case study of 2-chloropropionic acid

The vibrational circular dichroism (VCD) spectrum of S-(-) and R-(+)-2-chloropropionic acid is thoroughly analyzed. Besides the VCD spectrum of the monomer, the dimers (stabilized by strong hydrogen bonds) and the 2-chloropropionic acid-CHCl(3) complexes (stabilized by a weak hydrogen bond) are studied both experimentally (in solution and in low-temperature Ar matrix) and by quantum chemical computations. It is shown that dimer formation drastically changes, and even weak complex formation can also substantially affect the overall shape of the VCD spectrum. The present and previous results can be generalized for the practice of absolute configuration determination of carboxylic acids by VCD spectroscopy. For these measurements, if bulky groups do not block dimer formation, comparison of the computed spectra of the dimers with the experimental spectra recorded in relatively concentrated (～0.1 mol dm(-3)) solutions is suggested. Our study also shows that due to the stabilization of monomers and/or the formation of weak complexes, the VCD spectrum recorded in CHCl(3) is more complex and, like in the present case, can have a lower intensity than that of the spectrum recorded in CCl(4). Therefore, if solubility allows, CCl(4) is a much preferred solvent over CHCl(3).     

Temperature-dependent single-crystal neutron diffraction study of the strong OHN hydrogen bond in pyridinium 2,4-dinitrobenzoate

The structure of pyridinium 2,4-dinitrobenzoate was studied by neutron diffraction at 300, 270, 240, 210, 180, 150, 120, 90, 60, and 30 K. With temperature change, the O...H bond length changes from 1.403(10) A at 300 K to 1.424(4) A at 30 K. The proton shifts in the hydrogen bridge toward the acceptor nitrogen atom. Temperature-dependent changes in the strong OHN hydrogen bond were analyzed by using both the neutron structure and the atom-in-molecule approach. The results are compared with those for other strong OHN hydrogen bonds.     

Infrared and raman spectra of orthorhombic,monoclinic and cubic metaboric acid and their relation to the “strength” of the hydrogen bond present

Infrared and Raman spectra of the three forms of metaboric acid (orthorhombic, monoclinic and cubic) are presented and discussed. The behaviour of the asymmetric stretching vibration v_oh and of the in-plane and out-of-plane deformations δ_OH and γ_OH shows that hydrogen bonds of different “strengths” are present in the three forms of metaboric acid, in accordance with X-ray diffraction measurements. This phenomenon is explained in a qualitative way by a valence bond treatment of the structural centres present in the three forms of metaboric acid. The remaining vibrational modes are related to the analogous modes present in borates and polyborates of known structure.     

Phase transition in hydrazinium hydrogen oxalate. A calorimetric study of the short hydrogen bond system

The heat capacity of hydrazinium hydrogen oxalate crystal has been measured between 14 and 300 K. A first-order phase transition was found at 217.6 K. Total transition enthalpy and entropy are 1.09 kJ mol⁻¹ and 4.18 J K⁻¹ mol⁻¹, respectively. The high temperature phase supercooled to 165−170 K. The heat capacity of the supercooled high temperature phase was significantly larger than that of the low temperature phase. By fitting a Schottky heat capacity function to the heat capacity difference, a Schottky energy level was found at (135 ± 4) cm⁻¹ above the ground state and related to the tunneling splitting of the energy levels of the short acidic hydrogen bond. A thermodynamic model of the first-order transition is proposed in which the Schottky anomaly plays the main role. Far-infrared spectra of hydrazinium hydrogen oxalate are reported for the frequency range 400–30 cm⁻¹ at temperatures of 295 and 90 K.     

Theoretical study of hydrogen and deuterium bond in glutaric acid crystal dimer

In the spirit of the work of Blaise et al. [J Chem Phys, 2005, 122, 64306], we have extended their quantum theoretical approach by accounting for the intrinsic anharmonicity of the slow frequency mode, which is described by a Morse potential to reproduce the polarized infrared line shapes of glutaric acid dimer and its deuterium derivative at different temperatures. In this approach, the adiabatic approximation is performed for each separate H‐bond bridge of the dimer, and a strong nonadiabatic correction is introduced into the model via the resonant exchange between the fast mode excited states of the two moieties. Working within the strong anharmonic coupling theory, according to which the high‐frequency mode is anharmonically coupled to the H‐bond bridge, this approach incorporated the Davydov coupling between the excited states of the two moieties, the quantum direct and indirect dampings and the intrinsic anharmonicity of the H‐bond bridge. The spectral density was obtained within the linear response theory by Fourier transform of the damped autocorrelation functions. The numerical results show that the theoretical line shapes of the glutaric acid dimer are in fairly good agreement with the experimental ones. Using a minimum number of independent parameters, this theoretical approach fits correctly the experimental line shapes of the glutaric acid dimer. The effects of deuteration and temperature have been successfully reproduced by our calculations. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Complex dynamics of proton and deuteron transfer in double hydrogen bond of benzoic acid isotopes

This paper deals with spin-lattice relaxation due to classical jumps and incoherent tunnelling of protons and deuterons in hydrogen bonds in solids. An analysis of experimental spin-lattice relaxation data for carboxylic acids suggests that tunnelling does not contribute to spin-lattice relaxation above the temperature at which the thermal energy of molecules and the potential barrier height are equal. It is also shown that contributions to the spin-lattice relaxation rate due to classical motion and incoherent tunnelling in excited vibrational states are negligible for fast proton transfer. However, for deuterons this contribution to spin-lattice relaxation is significant.     

Spectroscopic studies of ionic solvation. XXI. A raman,infrared, and NMR study of sodium perchlorate solutions in nonaqueous solvents

Sodium perchlorate solutions in several nonaqueous solvents were examined by²³Na,³⁵Cl-NMR, infrared, and Raman spectroscopic techniques. The formation of contact ion pairs lowers the symmetry of ClO ₄ ^– ion from T_d to C_3v or C_2v provided that the interaction is fairly strong. This was manifested for NaClO₄ solutions in acetonitrile, tetrahydrofuran, and pyridine. Combination of the vibrational measurements with NMR shows that in the above three cases the anion-cation interactions are quite strong. Sodium-23 NMR studies confirm the above results and, being a more sensitive technique, also indicates weak cation-anion interaction in propylene carbonate, formic acid, acetone, methanol, ethanol, dimethyl sulfoxide, and water. In all solvents the²³Na resonance shifts upfield with increasing concentration of NaClO₄, indicating that the replacement of solvent by ClO ₄ ^– ion decreases electron density around the cation.     

Effect of temperature and pressure on the structure, dynamics, and hydrogen bond properties of liquid N-methylacetamide: a molecular dynamics study

The structure and dynamical properties of liquid N-methylacetamides (NMA) are calculated at five different temperatures and at four different pressures using classical molecular dynamics simulations. Our results are analyzed in terms of pressure-induced changes in structural properties by investigating the radial distribution functions of different atoms in NMA molecule. It is found that the first peak and also the second peak of C-O and N-H are well defined even at higher temperature and pressure. It is also observed that the number of hydrogen bonds increase with application of pressure at a given temperature. On the other hand, the calculated hydrogen bond energy (E(HB)) shows that the stability of hydrogen bond decreases with increasing of pressure and temperature. Various dynamical properties associated with translational and rotational motion of neat NMA are calculated and the self-diffusion coefficient of NMA is found to be in excellent agreement with the experiment and the behavior is non-Arrhenius at low temperatures with application of pressures. The single particle orientational relaxation time for dipole vector and N-C vector are also calculated and it is found that the orientational relaxation time follows Arrhenius behavior with a variation of temperature and pressure.     

A computational study of thiolate and selenolate oxidation by hydrogen peroxide.

Ab initio molecular orbital calculations have been used to study the effects of the molecular environment on the oxidation of thiolate and selenolate by hydrogen peroxide. The reaction was first examined in vacuo at the QCISD(T)/6-311+G(2df,2pd)//MP2/6-311+G(d,p) level of theory. It was found for both thiolate and selenolate that a reactant aggregate is formed, which has a dissociation rate constant comparable to the activation rate constant (about 10(-3) s(-1) for thiolate and 10(-1) s(-1) for selenolate). Using the polarizable continuum model (PCM) it was then found that the dissociation barrier energy decreases dramatically in water giving a dissociation rate constant of the order of 10(9) s(-1). In this case, the predicted overall rate constant of the thiolate reaction was about 10.2 mol(-1) dm3 s(-1), which is in good agreement with the experimental rate constant of cysteine oxidation in aqueous solution. The calculated rate constant for the selenolate reaction was somewhat higher (about 35.4 mol(-1) dm3 s(-1)). However, this value is several orders of magnitude smaller than the experimental value reported for the oxidation of selenocysteine in glutathione peroxidase. By considering the effect of the PCM dielectric constant on the reaction rate constant it was concluded that the high reactivity of the selenocysteine in glutathione peroxidase, as compared with cysteine, could be mainly due to the molecular environment of the selenocysteine residue.     

A theoretical study of the hydrogen bond donor capability and co-operative effects in the hydrogen bond complexes of the diaza-aromatic betacarbolines

In this work, we present a quantum mechanical investigation on the hydrogen bond interactions of N(9)-methyl-9H-pyrido[3,4-b]indole, MBC, and N(2)-methyl-9H-pyrido[3,4-b]indole, BCA, with different hydrogen bond donors. Thus, it has been analysed the influence that the hydrogen bond donor strength and the co-operative effect of the increasing number of donor molecules have on the shape of the potential energy surfaces versus the N···H distances, r(N–H). To rationalize the nature of the interactions, the Bader theory has been applied and the characteristics of the bond critical points analysed. The results show that two different hydrogen bond complexes can be formed depending on the donor capabilities or the number of donor molecules included in the calculations. The topological parameters from the Bader theory are used to justify the statement that the analysed interactions can be classified as weak or partially covalent hydrogen bond interactions, respectively. As experimentally observed, weak hydrogen bond donors form weak hydrogen bond complexes, called HBC. Upon the increase of the donor strength the N···H proton is shifted nearest to the nitrogen atom giving rise to the observation of a stronger hydrogen bond complex, the proton transfer complex, PTC. The most outstanding result of these studies is the fact that the formation of the PTC can also be managed just by changing the number of donor molecules, that is, by a co-operative effect of the hydrogen bonds.     

Strong intramolecular hydrogen bonds. Part I. Vibrational frequencies of the OH group in some picolinic acid N-oxides predicted from DFT calculated potentials and located in the infrared spectra

Hydrogen bonding in picolinic acid N-oxide (I), its 4-nitro (III), 4-methoxy (IV), 4-amino (V) derivatives and in quinaldic acid N-oxide (II) was characterized by calculations (B3LYP/6-31G(d)) of metric parameters, H-bond energies and one-dimensional proton potential functions with vibrational energy levels. Solvent effects were estimated by the SCRF PCM method of Tomasi and coworkers (J. Tomasi, M. Persico, Chem. Rev. 94 (1994) 2027). The potential functions are strongly asymmetric with the energy minimum placed near the carboxylic oxygen. The inflection near the NO oxygen develops into a second, shallower minimum under the SCRF.

Empirical assignments of the OH stretching and bending modes were made for (I)–(IV). The stretchings of (I, II) and (IV) in various solvents are observed in the region 1600–1300 cm⁻¹, but near 2600 cm⁻¹ for (III). The calculated and observed frequencies are in fairly good agreement with theoretical predictions reflecting the electronic effects of the substituents upon the H-bond strength. The observed trends in the solvent effects upon various parameters characterizing the H-bonding also correspond to predictions.