Theoretical model of infrared spectra of hydrogen bonds in molecular crystals and its application to interpretation of infrared spectra of 1-methylthymine

Spectroscopic study of uracil, 1-methyluracil and 1-methyl-4-thiouracil: Hydrogen bond interactions in crystals and ab-initio molecular dynamics

Hydrogen bond networks in uracil, 1-methyluracil and 1-methyl-4-thiouracil were studied by ab initio molecular dynamics as well as analysis of the orbital interactions. The power spectra calculated by ab initio molecular dynamics for atoms involved in hydrogen bonds were analyzed. We calculated spectra by using anharmonic approximation based on the autocorrelation function of the atom positions obtained from the Born-Oppenheimer simulations. Our results show the differences between hydrogen bond networks in uracil and its methylated derivatives. The studied methylated derivatives, 1-methyluracil as well as 1-methyl-4-thiouracil, form dimeric structures in the crystal phase, while uracil does not form that kind of structures. The presence of sulfur atom instead oxygen atom reflects weakness of the hydrogen bonds that build dimers.

Investigation of hydrogen bonds properties in the terephthalic acid crystal, using molecular dynamics method

The aim of this study was to perform calculations using the method of Car-Parrinello molecular dynamics, leading to the optimized geometry of the molecules of 1,4-benzenedicarboxylic acid (terephthalic acid) in crystals, for the hydrogen form and three variants of substitution of deuterium atoms inside a carboxyl group. Based on the results, trajectories and dipole moments were calculated, what makes possible to simulate vibrations in different systems, and to make calculation of theoretical infrared spectra and atomic power spectra. Theoretical results were compared with the experimental spectra, which verifies the correctness of the method and also was compared with the results obtained by quantum-mechanical calculations using DFT for the isolated dimer. Comparison of the spectra of different forms, allowed for in-depth analysis of the effect of isotopic substitution on the frequency of vibrations and shapes of bands, and confirm the presence of possible coupling effects and intra- and intermolecular interactions. Comparison with the DFT results for the dimer show influence of the crystal structure on the spectra.

Quantum chemical and spectroscopic investigations of 3-methyladenine

FTIR, FT-Raman and UV-Vis spectra of 3-methyladenine have been recorded and investigated using quantum chemical calculations. The molecular geometry and vibrational spectra of 3-methyladenine in the ground state are computed by using HF and DFT methods with 6-311G(d,p) basis set. VSCF, CC-VSCF methods based on 2MR-QFF and PT2 (Barone method) have been utilized for computing anharmonic vibrational frequencies. These methods yield results that are in remarkable agreement with the experimental data. The magnitudes of coupling between pair of modes have been also computed. Vibrational modes are assigned with the help of visual inspection of atomic displacements. The electronic spectra, simulated at TD-B3LYP/6-311++G(d,p) level of theory, are compared to the experiment. The global quantities: electronic chemical potential, electrophilicity index, chemical hardness and softness based on HOMO and LUMO energy eigenvalues are also computed at B3LYP/6-311++G(d,p) level of theory.

Car-Parrinello simulation of the vibrational spectrum of a medium strong hydrogen bond by two-dimensional quantization of the nuclear motion: application to 2-hydroxy-5-nitrobenzamide

The nature of medium strong intra- and intermolecular hydrogen bonding in 2-hydroxy-5-nitrobenzamide in the crystal phase was examined by infrared spectroscopy and Car-Parrinello molecular dynamics simulation. The focal point of our study was the part of the infrared spectra associated with the O-H and N-H stretching modes that are very sensitive to the strength of hydrogen bonding. For spectra calculations we used an isolated dimer and the fully periodic crystal environment. We calculated the spectra by using harmonic approximation, the time course of the dipole moment function as obtained from the Car-Parrinello simulation, and the quantization of the nuclear motion of the proton for an instantaneous snapshot of the structures in one and two dimensions. Although quantitative assessment of the agreement between the computed and experimental band contour is difficult due to the fact that the experimental band is very broad, we feel that the most reasonable qualitative agreement with the experiment is obtained from snapshot structures and two-dimensional quantization of the proton motion. We have also critically examined the methods of constructing the one-dimensional proton potential. Perspectives are given for the treatment of nuclear quantum effects in biocatalysis.

Car-Parrinello simulation of hydrogen bond dynamics in sodium hydrogen bissulfate

We studied proton dynamics of a short hydrogen bond of the crystalline sodium hydrogen bissulfate, a hydrogen-bonded ferroelectric system. Our approach was based on the established Car-Parrinello molecular dynamics (CPMD) methodology, followed by an a posteriori quantization of the OH stretching motion. The latter approach is based on snapshot structures taken from CPMD trajectory, calculation of proton potentials, and solving of the vibrational Schrodinger equation for each of the snapshot potentials. The so obtained contour of the OH stretching band has the center of gravity at about 1540 cm(-1) and a half width of about 700 cm(-1), which is in qualitative agreement with the experimental infrared spectrum. The corresponding values for the deuterated form are 1092 and 600 cm(-1), respectively. The hydrogen probability densities obtained by solving the vibrational Schrodinger equation allow for the evaluation of potential of mean force along the proton transfer coordinate. We demonstrate that for the present system the free energy profile is of the single-well type and features a broad and shallow minimum near the center of the hydrogen bond, allowing for frequent and barrierless proton (or deuteron) jumps. All the calculated time-averaged geometric parameters were in reasonable agreement with the experimental neutron diffraction data. As the present methodology for quantization of proton motion is applicable to a variety of hydrogen-bonded systems, it is promising for potential use in computational enzymology.

Theoretical modeling of the O–H stretching IR bands of hydrogen-bonded dimers of benzoic acid in S0 and S1 electronic states

Theoretical model for vibrational interactions in the hydrogen-bonded dimer of benzoic acid is presented. The model takes into account anharmonic-type couplings between the high-frequency O-H and the low-frequency O[cdots, three dots, centered]O stretching vibrations in two hydrogen bonds, resonance interactions (Davydov coupling) between two hydrogen bonds in the dimer, and Fermi resonance between the O-H stretching fundamental and the first overtone of the O-H in-plane bending vibrations. The vibrational Hamiltonians and selection rules for the C(2h) geometry in the S(0) state and for the C(s) in-plane bent geometry in the S(1) state of the dimer are derived. The model is used for theoretical simulation of the O-H stretching IR absorption bands of benzoic acid dimers in the gas phase in the electronic ground and first excited singlet states. Ab initio CIS and CIS(D)6-311++G(d,p) calculations have been performed to determine geometry, frequencies, and excited state energies of benzoic acid dimer in the S(1) state.