Car–Parrinello Molecular Dynamics Simulations of Infrared Spectra of Crystalline Vitamin C with Analysis of Double Minimum Proton Potentials for Medium-Strong Hydrogen Bonds

Renewed spectroscopic and theoretical research of hydrogen bonding in ascorbic acid

The studies of two isomers of ascorbic acid and their deuteroanalogues, presented in the paper, have been accomplished by vibrational spectroscopy methods and quantum-chemical simulations. The spectroscopic research of L-ascorbic and D-isoascorbic acids have been carried out by the infrared (IR) and Raman (R) techniques. On the basis of the obtained results the spectral interpretation of the hydrogen bonded groups of ascorbic acids has been performed. Car-Parrinello Molecular Dynamics (CPMD) and Density Functional Theory (DFT) have been employed to support spectroscopic experimental findings and shed light onto the bridged proton dynamics in the L- and D- isomers of ascorbic acids. The accurate assignments of the hydrogen bond modes have been accomplished with the application of deuterosubstitution, CPMD-solid state simulations and Potential Energy Distribution (PED) analysis. The spectral and structural results have shown that dependency ν(OH) = f(γ(OH)) is the most common for the OHO hydrogen bond, whereas dependency d(OO) = f(γ(OH)) differs as for the ionic and resonance assisted hydrogen bonds.

Very Strong Hydrogen Bond in Nitrophthalic Cocrystals

This work presents the studies of a very strong hydrogen bond (VSHB) in biologically active phthalic acids. Research on VSHB comes topical due to its participation in many biological processes. The studies cover the modelling of intermolecular interactions and phthalic acids with 2,4,6-collidine and N,N-dimethyl-4-pyridinamine complexes with aim to obtain a VSHB. The four synthesized complexes were studied by experimental X-ray, IR, and Raman methods, as well as theoretical Car-Parrinello Molecular Dynamics (CP-MD) and Density Functional Theory (DFT) simulations. By variation of the steric repulsion and basicity of the complex' components, a very short intramolecular hydrogen bond was achieved. The potential energy curves calculated by the DFT method were characterized by a low barrier (0.7 and 0.9 kcal/mol) on proton transfer in the OHN intermolecular hydrogen bond for 3-nitrophthalic acid with either 2,4,6-collidine or N,N-dimethyl-4-pyridinamine cocrystals. Moreover, the CP-MD simulations exposed very strong bridging proton dynamics in the intermolecular hydrogen bonds. The accomplished crystallographic and spectroscopic studies indicate that the OHO intramolecular hydrogen bond in 4-nitrophthalic cocrystals is VSHB. The influence of a strong steric effect on the geometry of the studied cocrystals and the stretching vibration bands of the carboxyl and carboxylate groups was elaborated.

Anharmonicity and Spectra-Structure Correlations in MIR and NIR Spectra of Crystalline Menadione (Vitamin K3)

Mid-infrared (MIR) and near-infrared (NIR) spectra of crystalline menadione (vitamin K3) were measured and analyzed with aid of quantum chemical calculations. The calculations were carried out using the harmonic approach for the periodic model of crystal lattice and the anharmonic DVPT2 calculations applied for the single molecule model. The theoretical spectra accurately reconstructed the experimental ones permitting for reliable assignment of the MIR and NIR bands. For the first time, a detailed analysis of the NIR spectrum of a molecular system based on a naphthoquinone moiety was performed to elucidate the relationship between the chemical structure of menadione and the origin of the overtones and combination bands. In addition, the importance of these bands during interpretation of the MIR spectrum was demonstrated. The overtones and combination bands contribute to 46.4% of the total intensity of menadione in the range of 3600-2600 cm-1. Evidently, these bands play a key role in shaping of the C-H stretching region of MIR spectrum. We have shown also that the spectral regions without fundamentals may provide valuable structural information. For example, the theoretical calculations reliably reconstructed numerous overtones and combination bands in the 4000-3600 and 2800-1800 cm-1 ranges. These results, provide a comprehensive origin of the fundamentals, overtones and combination bands in the NIR and MIR spectra of menadione, and the relationship of these spectral features with the molecular structure.

Understanding the Structure of the Hydrogen Bond Network and Its Influence on Vibrational Spectra in a Prototypical Aprotic Ionic Liquid

Analysis of the hydrogen bond network in aprotic ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) has been performed based on structures obtained from ab initio or classical molecular dynamics simulations. Statistics of different donor and acceptor atoms and the amount of chelating or bifurcated bonds has been presented. Most of the hydrogen bonds in EMIM-TFSI are formed with oxygen atoms as hydrogen acceptors; and the most probable bifurcated bonds are those with a mixed pair of oxygen and nitrogen acceptors. Spectral graph analysis has shown that the cations may form hydrogen bonds with up to five different anions and the connectivity of the whole hydrogen bond network is supported mainly by H-O bonds. In the structures of the liquid simulated via force field-based dynamics, the number of hydrogen bonds is smaller and fluorine atoms are the most favored hydrogen acceptors. One-dimensional potential energy profiles for hydrogen atom displacements and corresponding vibrational frequencies have been calculated for selected C-H bonds. Individual C-H stretching frequencies vary by 200-300 cm^-1, indicating differences in local environment of hydrogen atoms forming C-H···O hydrogen bonds.

Physics behind the Barrier to Internal Rotation of an Acetyl Chloride Molecule: A Combined Approach from Density Functional Theory, Car-Parrinello Molecular Dynamics, and Time-Resolved Wavelet Transform Theory

The physics behind the barriers to internal rotation of acetyl chloride (AC) molecule has been reported. The AC molecule closely resembles the molecular structure of acetaldehyde; the only subtle difference is the presence of a heavy chlorine atom in place of the hydrogen atom of the aldehyde group for the latter. This paper aims to study the effect of substitution of the heavy chlorine atom on the barrier energetics of the AC molecule. The reason behind the barrier for the AC molecule has been estimated for the first time from the unified approach using barrier energetics, natural bond orbital, nuclear virial, and relaxation analyses using density functional theory, Car-Parrinello molecular dynamics, and wavelet transform theory. Complete analyses reveal the concomitant relaxations of both the in-plane C_methyl-C₁ and C_methyl-H₄ bonds toward understanding the origin of the barrier due to internal rotation for the AC molecule. The large negative value of "V ₆" further suggests that both the abovementioned degrees of freedom are coupled with the -CH₃ torsional vibration of the molecule. The coupling matrix (H ₁₂) element has also been estimated. Time-resolved band stretching frequencies of C_methyl-C₁ and C₁-Cl₃ bonds of the AC molecule, as obtained from wavelet transformation analysis, primarily preclude the possibility of coupling between the C₁-Cl₃ bond and the torsional motion associated with the methyl group of the molecule.

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.

Specific features of supramolecular organisation and hydrogen bonding in proline cocrystals: a case study of fenamates and diclofenac

New zwitterionic cocrystals of fenamate drugs and diclofenac with the naturally occurring amino acid l-proline have been obtained and thoroughly characterised by a variety of experimental and theoretical techniques.

Born-Oppenheimer Molecular Dynamics Study on Proton Dynamics of Strong Hydrogen Bonds in Aspirin Crystals, with Emphasis on Differences between Two Crystal Forms

In this study, the proton dynamics of hydrogen bonds for two forms of crystalline aspirin was investigated by the Born-Oppenheimer molecular dynamics (BOMD) method. Analysis of the geometrical parameters of hydrogen bonds using BOMD reveals significant differences in hydrogen bonding between the two crystalline forms of aspirin, Form I and Form II. Analysis of the trajectory for Form I shows spontaneous proton transfer in cyclic dimers, which is absent in Form II. Quantization of the O-H stretching modes allows a detailed discussion on the strength of hydrogen-bonding interactions. The focal point of our study is examination of the hydrogen bond characteristics in the crystal structure and clarification of the influence of hydrogen bonding on the presence of the two crystalline forms of aspirin. In the BOMD method, thermal motions were taken into account. Solving the Schrödinger equation for the snapshots of 2D proton potentials, extracted from MD, gives the best agreement with IR spectra. The character of medium-strong hydrogen bonds in Form I of aspirin was compared with that of weaker hydrogen bonds in aspirin Form II. Two proton minima are present in the potential function for the hydrogen bonds in Form I. The band contours, calculated by using one- and two-dimensional O-H quantization, reflect the differences in the hydrogen bond strengths between the two crystalline forms of aspirin, as well as the strong hydrogen bonding in the cyclic dimers of Form I and the medium-strong hydrogen bonding in Form II.