Effect of Waters of Crystallization on Terahertz Spectra: Anhydrous Oxalic Acid and Its Dihydrate

Strong Hydrogen Bonds in Acetylenedicarboxylic Acid Dihydrate

Acetylenedicarboxylic acid dihydrate (ADAD) represents a complex with strong hydrogen bonding between the carboxylic OH and the water molecule. An X-ray re-examination of the ADAD crystal structure confirms the O^…O distance of the short hydrogen bonds, and clearly shows different bond lengths between the two oxygen atoms with respect to the carbon atom in the carboxyl group, indicating a neutral structure for the complex. The neutral structure was also confirmed by vibrational spectroscopy, as no proton transfer was observed. The diffraction studies also revealed two polymorph modifications: room temperature (α) and low temperature (β), with a phase transition at approximately 4.9 °C. The calculated vibrational spectra are in satisfactory agreement with the experimental spectra. A comparison of the structure and the vibrational spectra between the ADAD and the oxalic acid dihydrate reveals some interesting details. The crystal structures of both crystal hydrates are almost identical; only the O^…O distances of the strongest hydrogen bonds differ by 0.08 Å. Although it was expected that a larger O^…O spacing in the ADAD crystal may significantly change the infrared and Raman spectra, especially for the frequency and the shape of the acidic OH stretching vibration, both the shape and frequency are almost identical, with all subpeaks topped on the broad OH stretching vibration. The O^…O distance dependent are only in- and out-of-plane OH deformations modes. The presence of polarons due to the ionized defects was not observed in the vibrational spectra of ADAD. Therefore, the origin of the broad OH band shape was explained in a similar way to the acid dimers. The anharmonicity of a potential enhances the coupling of the OH stretching with the low-frequency hydrogen bond stretching, which, in addition to the Fermi resonance, structures the band shape of the OH stretching. The fine structure found as a superposition of a broad OH stretching is attributed to Davydov coupling.

Study on Hydration and Dehydration of Ezetimibe by Terahertz Spectroscopy with Humidity-Controlled Measurements and Theoretical Analysis

Understanding the solid-state transitions of active pharmaceutical ingredients (APIs) is essential for quality control since differences in their forms affect the bioavailability of APIs. Terahertz (THz) frequency-domain spectroscopy is suitable for such an application since it can sensitively probe the lattice phonon modes originating in the crystal structures. THz absorption spectra were obtained for ezetimibe (EZT) and ezetimibe monohydrate (EZT-MH), which have similar crystalline structures and belong to the same space group. The observed absorption spectrum of EZT matched well with the solid-state density functional theory (ss-DFT)-simulated spectrum for the structures at 0 K and room temperature (modeled using constrained unit cell volumes). However for EZT-MH, the ss-DFT spectrum of the room-temperature structure showed better correlation with the experimental THz spectrum than that of the simulated spectrum of the 0 K structures, suggesting that the EZT-MH crystal has greater anharmonic character. Gibbs free-energy curves were calculated, and EZT-MH was found to be more stable than pure EZT and water in a broad temperature range. The hydrate stability may be influenced by the existence of more hydrogen bonds in EZT-MH. The hydration and dehydration of EZT in a pure API tablet and formulation tablets were monitored using a THz spectrometer with a humidity-controlled sample chamber. The effect of the excipient in the formulation tablet on hydration was successfully confirmed by showing that the solid-state transition of the API with excipients is significantly slower than that without it. Under a relative humidity of 60%, hydration of EZT in a pure EZT tablet occurred in 200 min, while the hydration of EZT in a formulation tablet was 50 times slower.

Terahertz spectroscopy for interpreting the formation and hierarchical structures of silk fibroin oligopeptides

Determining the configuration and conformation of peptides is crucial for interpreting their structure-property relationships. In this work, we present nondestructive terahertz time-domain spectroscopy combined with density functional theory (DFT) and potential energy distribution (PED) analysis to identify the hierarchical structures of oligopeptides. The characteristic THz spectra of silk fibroin oligopeptides have been measured. Supported by DFT and PED analysis, the intrinsic differences among the dipeptides were identified by the collective vibrational modes of "R" groups and terminal groups linked by molecular chains of amido bonds or benzene rings. For tetrapeptides and hexapeptides, a few weak resonances and intensity differences were distinguished by the vibration mode of the molecular collective network formed by the interaction of amide planes and intramolecular hydrogen bond interactions. According to the THz absorption analyses of amide planes and intramolecular interactions within the molecular chains of silk fibroin oligopeptide isomer pairs, the formation and hierarchical structures were successfully interpreted using THz spectroscopy. This investigation develops a better understanding of the peptide formation mechanism, which further provides guidance in interpreting the formation of silk.

Examination of proline, hydroxyproline and pyroglutamic acid with different polar groups by terahertz spectroscopy

Hydroxyproline (HYP) and pyroglutamic acid (PGA), as amino acid derivatives, are highly similar in structure to proline (Pro). However, their low-frequency vibrations show significant differences in the range of 0.25-2.6 THz. Therefore, this study investigated the reasons for the differences combined with terahertz time domain spectroscopy (THz-TDS) and density functional theory (DFT). The results show that HYP and PGA have stronger absorption of terahertz waves due to the existence of polar substituents. Furthermore, the absorption peaks of HYP and PGA are significant red shifted and blue shifted, respectively. We believe that this is caused by the change in the strength of intermolecular hydrogen bonds. Our findings demonstrate that dipole and hydrogen bond effects play a significant role in low-frequency vibrations.

Low-frequency vibrational spectroscopy: a new tool for revealing crystalline magnetic structures in iron phosphate crystals

In this report, the strong-dependence of low-frequency (terahertz) vibrational dynamics on weak and long-range forces in crystals is leveraged to determine the bulk magnetic configuration of iron phosphate - a promising material for cathodes in lithium ion batteries. We demonstrate that terahertz time-domain spectroscopy - coupled with quantum mechanical simulations - can discern between various spin configurations in FePO₄. Furthermore, the results of this work unambiguously show that the well-accepted space group symmetry for FePO₄ is incorrect, and the low-frequency spectroscopic measurements provide a clearer picture of the correct structure over the gold-standard of X-ray diffraction. This work opens the door for characterizing, predicting, and interpreting crystalline magnetic ordering using low-frequency vibrational spectroscopy.

Discrimination of Pharmaceutical Tablets Based on the Analysis of Solid-State Structures of Ingredients Using Terahertz Transmission Spectroscopy with the Injection-Seeded Parametric Generation Technique

A frequency-domain terahertz (THz) spectrometer that uses a tunable source, called an injection-seeded THz parametric generator, was applied to the analysis of solid-state structures of ingredients in pharmaceutical tablets, and its performance on discriminating pharmaceutical products was evaluated. The spectrometer has a dynamic range of 70 dB at 2 THz and is suitable for analyzing materials such as pharmaceutical ingredients that often have characteristic absorption peaks between 0.5 and 2.5 THz. Nine ofloxacin (racemate) and four levofloxacin (levorotatory enantiomer) tablet products commercially available in Japan were used as samples. They contain 8-12 additives in addition to the API. The sample tablets were filed down to a thickness of 1.2 mm (ofloxacin tablets) and 1.6 mm (levofloxacin tablets) to obtain transmission spectra over the wide spectral range of 0.8-2.1 THz. The absorption spectra obtained from the spectrometer were preprocessed by the second derivative; then, principal component analysis (PCA) was conducted on the results. Next, quadratic discriminant analysis (DA) was conducted on the scores of the three PCA components. The accuracy of the DA for all 13 products was 96.1%. In addition to the difference in crystal forms of the active ingredient, the small differences in the formulation were clearly discriminated using the THz absorption spectra. The spectrometer combined with data analysis shows potential for applications such as identifying pharmaceutical tablets, monitoring the stability of production processes, evaluating the stability of formulations during storage, and detecting counterfeit drugs on the market.

The necessity of periodic boundary conditions for the accurate calculation of crystalline terahertz spectra

Terahertz vibrational spectroscopy has emerged as a powerful spectroscopic technique, providing valuable information regarding long-range interactions - and associated collective dynamics - occurring in solids. However, the terahertz sciences are relatively nascent, and there have been significant advances over the last several decades that have profoundly influenced the interpretation and assignment of experimental terahertz spectra. Specifically, because there do not exist any functional group or material-specific terahertz transitions, it is not possible to interpret experimental spectra without additional analysis, specifically, computational simulations. Over the years simulations utilizing periodic boundary conditions have proven to be most successful for reproducing experimental terahertz dynamics, due to the ability of the calculations to accurately take long-range forces into account. On the other hand, there are numerous reports in the literature that utilize gas phase cluster geometries, to varying levels of apparent success. This perspective will provide a concise introduction into the terahertz sciences, specifically terahertz spectroscopy, followed by an evaluation of gas phase and periodic simulations for the assignment of crystalline terahertz spectra, highlighting potential pitfalls and good practice for future endeavors.

Infrared spectra of hydrogen bond network in lamellar perfluorocarboxylic acid monohydrates

The infrared spectra of the long-chain perfluorocarboxylic acid monohydrates differ markedly from those of the anhydrous dimers. Consequently, the structure of the solid perfluorocarboxylic acid monohydrates must differ from any known dimer-containing carboxylic acid crystals. Consideration of the significant features of the infrared spectra of the long-chain perfluorocarboxylic acid monohydrates, supplemented by their Raman spectra, and comparison with the spectra of auxiliary substances have led us to conclude that the rather strong neutral carboxyl-hydroxyl to water bonding can best explain the observations. The infrared spectra indicate the presence of fairly short hydrogen bonds connecting the water molecules to the carbonyl groups. In the construction of the hydrogen bonding pattern of the perfluorocarboxylic acid monohydrates, the oxalic acid dihydrate plays the key role. The striking similarity between the infrared spectra of the oxalic acid dihydrates and the perfluorocarboxylic acid monohydrates in the regions characteristic of water and OH⋯O vibration suggests that the structure of the hydrated carboxyl groups is the same in both crystals. These regions are characterized by the sharp doublet at 3539 cm^-1 and 3464 cm^-1, which is due to the H₂O ν₁ and ν₃ stretching vibrations, respectively, and the broad absorption between 3000 cm^-1 and 1500 cm^-1 with the intense band at 1970 cm^-1, both associated with the vibration of the OH⋯O group. The later peak consists of two band components at near 1980 cm^-1 and 2020 cm^-1. These band components show different behaviour when the temperature, polarization or deuteration is changed. In general, the infrared spectra of long-chain perfluorocarboxylic acids represent the system with very short hydrogen bonds connecting the water molecules to the carboxylates. This hydrogen bond pattern should be very similar to that found in the crystals of α-oxalic acid dihydrate.

Transformation of terahertz vibrational modes of cytosine under hydration

Cytosine and cytosine monohydrate are representative biomolecules for investigating the effect of hydrogen bonds in deoxyribonucleic acid. To better understand intermolecular interactions, such as hydrogen bonds, between nucleobases it is necessary to identify the low-frequency vibrational modes associated with intermolecular interactions and crystalline structures. In this study, we investigated the characteristic low-frequency vibrational modes of cytosine and cytosine monohydrate using terahertz time-domain spectroscopy (THz-TDS). The crystal geometry was obtained by the powder X-ray diffraction technique. The optimized atomic positions and the normal modes in the terahertz region were calculated using density functional theory (DFT), which agreed well with the experimental results. We found that overall terahertz absorption peaks of cytosine and cytosine monohydrate consist of collective vibrations mixed with intermolecular and intramolecular vibrations in mode character analysis, and that the most intense peaks of both samples involve remarkable intermolecular translational vibration. These results indicate that THz-TDS combined with DFT calculations including mode character analysis can be an effective method for understanding how water molecules contribute to the characteristics of the low-frequency vibrational modes by intermolecular vibrations with hydrogen bonding in biological and biomedical applications.

Analysis of the polarized IR reflectance spectra of the monoclinic α-oxalic acid dihydrate

In contrast to the extensive molecular and crystal structure investigations on oxalic acid dihydrate (C₂H₂O₄·2H₂O, α-POX) and its deuterated analogues (α-DOX), stands the absence of a complete vibrational spectra analysis. Such analysis is desirable in view of the proton dynamics in α-POX crystals. In this communication we report the room temperature polarized IR reflectance spectra of a single crystal of α-POX recorded from the ac crystal plane, and from the plane containing the b-crystallographic axis, with polarization along the axis. Dispersion analysis of the reflectance spectra of both B_u and A_u symmetry type modes, using model dielectric and reflectance function valid for the monoclinic case, have been performed and the results are discussed. Some aspects of the obtained fit results for some of the spectral regions and the peculiar change of the reflectance function with polarization angle are also discussed in this work. A correlation between the crystal structure and measured spectra, together with the results of the performed dispersion analysis, gave answers to some of the problems concerned with the orientation of the transition dipole moments of the IR active modes. The assignment of modes is assisted by DFT calculations. Another aspect covered in this work is the model reflectance functions using different averaging theories that have been applied in obtaining the reflectance spectrum of a polycrystalline sample. The results of the comparison between these spectra and the recorded reflectance spectra from a polycrystalline sample were further discussed.

Characteristic fingerprint spectrum of neurotransmitter norepinephrine with broadband terahertz time-domain spectroscopy

A broadband THz absorption spectrum of NE corresponds to specific vibrational modes.

Experimental and theoretical study of fingerprint spectra of 2-(4-fluorophenyl)benzimidazole and 2-(4-chlorophenyl)benzimidazole in terahertz range

Due to wide variety of biological and pharmacological activities of benzimidazole derivatives, the terahertz fingerprint spectra of 2-(4-fluorophenyl)benzimidazole and 2-(4-chlorophenyl)benzimidazole are researched by employing terahertz time-domain spectroscopy and density functional theory systematically. Although the substituent at the para position on the benzene moiety are the same family elements (F, Cl) in the periodic table, both experimental and theoretical results show that there are substantial differences in their fingerprint spectra in the range of 0.2-2.5 THz, such as the amount, amplitude and frequency position of absorption peaks. The validity of these results was confirmed by isolated-molecule and solid-state calculations based on density functional theory. The possible reasons of these differences are different intramolecular hyperconjugative interactions, different elongation of the corresponding bond lengths, different HOMO-LUMO energy gaps, as well as different atomic motions within in the unit cell owing to the electron-withdrawing effects of different halogen atoms at the para position on the benzene moiety. These results indicate the importance of this spectral range as a conformational fingerprint region where even minor changes in the molecular configuration lead to major differences in its THz absorption.

Evidence of Polaron Excitations in Low Temperature Raman Spectra of Oxalic Acid Dihydrate

Low temperature Raman spectra of oxalic acid dihydrate (8-300 K) for both the polycrystalline and single crystal phase show strong variation with temperature in the interval from 1200 to 2000 cm(-1). Previous low temperature diffraction studies all confirmed the stability of the crystal P21/n phase with no indications of any phase transition, reporting the existence of a strong hydrogen bond between the oxalic acid and a water molecule. A new group of Raman bands in the 1200-1300 cm(-1) interval below 90 K is observed, caused by possible loss of the center of inversion. This in turn could originate either due to disorder in hydroxyl proton positions or due to proton transfer from carboxylic group to water molecule. The hypothesis of proton transfer is further supported by the emergence of new bands centered at 1600 and 1813 cm(-1), which can be explained with vibrations of H3O(+) ions. The broad band at 1600 cm(-1) looses intensity, while the band at 1813 cm(-1) gains intensity on cooling. The agreement between quantum calculations of vibrational spectra and experimentally observed Raman bands of hydronium ions in oxalic acid sesquihydrate crystal corroborates this hypothesis.

Vibrational spectroscopic study of polymorphism and polymorphic transformation of the anti-viral drug lamivudine

Vibrational spectra of hydrated and anhydrous lamivudines, and also the dynamic process of polymorphic transformation have been characterized by Fourier transform infrared (FT-IR) and Raman spectroscopic techniques. The vibrational modes of both polymorphic lamivudines are assigned. FT-IR and Raman spectral results show that the interaction between crystalline water and lamivudine molecular has an important effect on the molecular vibration motions of polymorphic lamivudines. The two characteristic Raman peaks at 783 and 798 cm(-1) represent hydrated and anhydrous lamivudine respectively. The relationship between changes of two characteristic peak normalized areas and heating time could be fitted with single exponential functions, and the dynamic information of polymorphic transformation of lamivudine drug is obtained. The decay rate of characteristic peak for hydrated lamivudine and the growth rate of that for anhydrous lamivudine are consistent during dehydration transformation process. The reported results provide us important benchmark for qualitatively monitoring different polymorphic drugs and also establishing the corresponding model for the polymorphic transformation of drugs in related pharmaceutical research fields.

The applicability of the dimeric heterosynthon concept to molecules with equivalent binding sites. A DFT study of crystalline urea–H2O2

The limited applicability of the dimeric heterosynthon concept to a two-component urea–H₂O₂crystal is reported.

Terahertz time-domain and low-frequency Raman spectroscopy of organic materials

With the ongoing proliferation of terahertz time-domain instrumentation from semiconductor physics into applied spectroscopy over the past decade, measurements at terahertz frequencies (1 THz ≡ 10(12) Hz ≡ 33 cm(-1)) have attracted a sustained growing interest, in particular the investigation of hydrogen-bonding interactions in organic materials. More recently, the availability of Raman spectrometers that are readily able to measure in the equivalent spectral region very close to the elastic scattering background has also grown significantly. This development has led to renewed efforts in performing spectroscopy at the interface between dielectric relaxation phenomena and vibrational spectroscopy. In this review, we briefly outline the underlying technology, the physical phenomena governing the light-matter interaction at terahertz frequencies, recent examples of spectroscopic studies, and the current state of the art in assigning spectral features to vibrational modes based on computational techniques.

Cl···Cl interactions in molecular crystals: insights from the theoretical charge density analysis

The structure, IR harmonic frequencies and intensities of normal vibrations of 20 molecular crystals with the X-Cl···Cl-X contacts of different types, where X = C, Cl, and F and the Cl···Cl distance varying from ~3.0 to ~4.0 Å, are computed using the solid-state DFT method. The obtained crystalline wave functions have been further used to define and describe quantitatively the Cl···Cl interactions via the electron-density features at the Cl···Cl bond critical points. We found that the electron-density at the bond critical point is almost independent of the particular type of the contact or hybridization of the ipso carbon atom. The energy of Cl···Cl interactions, E(int), is evaluated from the linking E(int) and local electronic kinetic energy density at the Cl···Cl bond critical points. E(int) varies from 2 to 12 kJ/mol. The applicability of the geometrical criterion for the detection of the Cl···Cl interactions in crystals with two or more intermolecular Cl···Cl contacts for the unique chlorine atom is not straightforward. The detection of these interactions in such crystals may be done by the quantum-topological analysis of the periodic electron density.

How strong is hydrogen bonding in ionic liquids? Combined X-ray crystallographic, infrared/Raman spectroscopic, and density functional theory study

Hydrogen bonding in ionic liquids based on the 1-(2'-hydroxylethyl)-3-methylimidazolium cation ([C₂OHmim](+)) and various anions ([A](-)) of differing H-bond acceptor strength, viz. hexafluorophosphate [PF6](-), tetrafluoroborate [BF₄](-), bis(trifluoromethanesulfonimide) [Tf₂N](-), trifluoromethylsulfonate [OTf](-), and trifluoroacetate [TFA](-), was studied by a range of spectroscopic and computational techniques and, in the case of [C₂OHmim][PF6], by single crystal X-ray diffraction. The first quantitative estimates of the energy (E(HB)) and the enthalpy (-ΔH(HB)) of H-bonds in bulk ILs were obtained from a theoretical analysis of the solid-state electron-density map of crystalline [C₂OHmim][PF6] and an analysis of the IR spectra in crystal and liquid samples. E(HB) for OH···[PF6](-) H-bonds amounts to ~3.4-3.8 kcal·mol(-1), whereas weaker H-bonds (2.8-3.1 kcal·mol(-1)) are formed between aromatic C2H group of imidazolium ring and the [PF6](-) anion. The enthalpy of the OH···[A](-) H-bonds follows the order: [PF6] (2.4 kcal·mol(-1)) < [BF₄] (3.3 kcal·mol(-1)) < [Tf₂N] (3.4 kcal·mol(-1)) < [OTf] (4.7 kcal·mol(-1)l) < [TFA] (6.2 kcal·mol(-1)). The formation of aggregates of self-associated [C₂OHmim](+) cations is present in liquid [C₂OHmim][PF6], [C₂OHmim][BF₄], and [C₂OHmim][Tf₂N], with the energy of the OH···OH H-bonds amounting to ~6 kcal·mol(-1). Multiple secondary interactions in the bulk ILs influence their structure, vibrational spectra, and H-bond strength. In particular, these interactions can blue-shift the stretching frequencies of the CH groups of the imidazolium ring in spite of red-shifting CH···[A](-) H-bonds. They also weaken the H-bonding in the IL relative to the isolated ion pairs, with these anticooperative effects amounting to ca. 50% of the E(HB) value.

Using terahertz time-domain spectroscopical technique to monitor cocrystal formation between piracetam and 2,5-dihydroxybenzoic acid

Far-infrared vibrational absorption of cocrystal formation between 2,5-dihydroxybenzoic acid (2,5-DHBA) and piracetam compounds under solvent evaporation and grinding methods have been investigated using terahertz time-domain spectroscopy (THz-TDS) at room temperature. The experimental results show large difference among absorption spectra of the formed cocrystals and the involved individual parent molecules in 0.20-1.50 THz region, which probably originated from the intra-molecular and inter-molecular hydrogen bonds due to the presence of two hydroxyl groups in 2,5-DHBA and amide moieties in piracetam compound. The THz absorption spectra of two formed cocrystals with different methods are almost identical. With grinding method, the reaction process can be monitored directly from both time-domain and frequency-domain spectra using THz-TDS technique. The results indicate that THz-TDS technology can absolutely offer us a high potential method to identify and characterize the formed cocrystals, and also provide the rich information about their reaction dynamic process involving two or more molecular crystals in situ to better know the corresponding reaction mechanism in pharmaceutical fields.