Investigating the Anharmonicity of Lattice Vibrations in Water-Containing Molecular Crystals through the Terahertz Spectroscopy of l-Serine Monohydrate

THz spectroscopy on the amino acids L-serine and L-cysteine

We present a detailed study on the temperature-dependent THz spectra of the polycrystalline amino acids, L-serine and L-cysteine, for wavenumbers from 20 to 120 cm-1 and temperatures from 4 to 300 K. Even though the structure of these two amino acids is very similar, with a sulfur atom in the side chain of cysteine instead of an oxygen atom in serine, the excitation spectra are drastically different. Obviously, the vibrational dynamics strongly depend on the ability of cysteine to form sulfur-hydrogen bonds. In addition, cysteine undergoes an order-disorder type phase transition close to 80 K, documented by additional specific heat experiments, with accompanying anomalies in the THz results. On increasing temperatures, well-defined vibrational excitations exhibit significant shifts in the eigenfrequencies with concomitant line-broadening yielding partly overlapping modes. Interestingly, several modes completely lose all their dipolar strength and are unobservable under ambient conditions. Comparing the recent results to the published work utilizing THz, Raman, and neutron-scattering techniques, as well as with ab initio simulations, we aim at a consistent analysis of the results ascribing certain eigenfrequencies to distinct collective lattice modes. We document that THz spectra can be used to fine-tune the parameters of model calculations and as fingerprint properties of certain amino acids. In addition, we analyzed the low-temperature heat capacity of both the compounds and detected strong excess contributions compared to the canonical Debye behavior of crystalline solids, indicating soft excitations and a strongly enhanced phonon-density of states at low frequencies.

Anharmonic Coupling of Stretching Vibrations in Ice: A Periodic VSCF and VCI Description

The anharmonicity of O-H stretching vibrations of water ice is characterized by use of a periodic implementation of the vibrational self-consistent field (VSCF) and vibrational configuration interaction (VCI) methods, which take phonon-phonon couplings explicitly into account through numerical evaluation of high-order terms of the nuclear potential. The low-temperature, proton-ordered phase of water ice (namely, ice XI) is investigated. The net effect of a coupled anharmonic treatment of stretching modes is not just a rigid blue-shift of the respective harmonic spectral frequencies but rather a complex change of their relative spectral positions, which cannot be captured by simple scaling strategies based on harmonic calculations. The adopted techniques allow for a hierarchical treatment of anharmonic terms of the nuclear potential, which is key to an effective identification of leading factors. We show that the anharmonic independent-mode approximation─only describing the "intrinsic anharmonicity" of the O-H stretches─is unable to capture the correct physics, and that couplings among O-H stretches must be described. Inspection of harmonic normal coordinates allows identification of specific features of the O-H stretching motions which most likely enable strong mode-mode couplings. Finally, by coupling O-H stretches to all other possible modes of ice XI (THz collective vibrations, molecular librations, bendings), we identify specific types of motion which significantly affect O-H stretching states: in particular, molecular librations are found to affect the stretching states more than molecular bendings.

Understanding the Metastability of Theophylline FIII by Means of Low-Frequency Vibrational Spectroscopy

While theophylline has been extensively studied with multiple polymorphs discovered, there is still currently no conclusive structure for the metastable theophylline form III. In this present work, by combining more widely used techniques such as X-ray diffraction and thermogravimetric analysis with more emerging techniques like low-frequency Raman and terahertz time-domain spectroscopy, to analyze the structure and dynamics of a crystalline system, it was possible to provide further evidence that the form III structure has a theophylline monohydrate structure with the water molecules removed. Solid-state density functional theory simulations were paramount in proving that this proposed structure is correct and explain how vibrational modes within the crystal structures feature and govern polymorphic transitions and the metastable form III. Through the insight provided by both simulated and experimental results, it was possible to decisively conclude the elusive crystal structure of theophylline form III. It was also shown that the correct space group for theophylline monohydrate is not P21/n but, in fact, Pc.

Condensation Effects on Electron Chiral Asymmetries in the Photoionization of Serine: From Free Molecules to Nanoparticles

Structural changes at the molecular level, occurring at the onset of condensation, can be probed by angle-resolved valence photoelectron spectroscopy, which is inherently sensitive to the electronic structure. For larger condensed systems like aerosol particles, the observation of intrinsic angular anisotropies in photoemission (β parameters) is challenging due to the strong reduction of their magnitude by electron transport effects. Here, we use a less common, more sensitive observable in the form of the chiral asymmetry parameter to perform a comparative study of the VUV photoelectron spectroscopy and photoelectron circular dichroism (PECD) between pure gas phase enantiomers of the amino acid serine and their corresponding homochiral nanoparticles. We observe a relatively large (1%) and strongly kinetic energy-dependent asymmetry, discussed in terms of the emergence of local order and conformational changes potentially counterbalancing the loss of angular information due to electron transport scattering. This demonstrates the potential of PECD as a sensitive probe of the condensation effects from the gas phase to bulk-like chiral aerosol particles surpassing the potential of conventional photoemission observables such as β parameters.

Resolving Anharmonic Lattice Dynamics in Molecular Crystals with X-Ray Diffraction and Terahertz Spectroscopy

Molecular crystals are increasingly being used for advanced applications, ranging from pharmaceutics to organic electronics, with their utility dictated by a combination of their three-dimensional structures and molecular dynamics-with anharmonicity in the low-frequency vibrations crucial to numerous bulk phenomena. Through the use of temperature-dependent x-ray diffraction and terahertz time-domain spectroscopy, the structures and dynamics of a pair of isomeric molecular crystals exhibiting nearly free rotation of a CF_{3} functional group at ambient conditions are fully characterized. Using a recently developed solid-state anharmonic vibrational correction, and applying it to a molecular crystal for the first time, the temperature-dependent spatial displacements of atoms along particular terahertz modes are obtained, and are found to be in excellent agreement with the experimental observations, including the assignment of a previously unexplained absorption feature in the low-frequency spectrum of one of the solids.

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.

Effects of Solvent Stabilization on Pharmaceutical Crystallization: Investigating Conformational Polymorphism of Probucol Using Combined Solid-State Density Functional Theory, Molecular Dynamics, and Terahertz Spectroscopy

Solid-state density functional theory (DFT), molecular dynamics (MD), and terahertz (THz) spectroscopy were used to study the formation of enantiotropically related conformational Form I and Form II polymorphs of the pharmaceutical compound, probucol. DFT calculations were performed on the crystal systems to compare relative lattice energies and the solvent stabilization of the metastable Form II structure. The thermodynamics of solvent inclusion in the Form II·MeOH crystal system were determined from MD simulations, as was the favored conformation of molecular probucol in methanol and ethanol solutions. The findings from both solid-state DFT and MD calculations suggest that the preferred molecular orientations of the probucol molecule in solution and the probable inclusion of methanol in the crystal lattice during the crystallization process lead to the solvent selectivity of the probucol polymorph formation. The additional stabilization energy provided by the crystallization solvent facilitates the nucleation and growth of the Form II polymorph under conditions that favor this metastable crystal form over the thermodynamically stable Form I, despite the higher energy molecular and crystalline configurations of probucol Form II. We demonstrate the influence of solvent on the formation of pharmaceutical polymorphs and provide a molecular-level view of complex interactions leading to polymorphism using a combination of computational methods and THz spectral data.

The effect of the flexibility of hydrogen bonding network on low-frequency motions of amino acids. Evidence from Terahertz spectroscopy and DFT calculations

Low-frequency modes of L-Asp and L-Asn were studied in the range from 0.1 to 3.0THz using time-domain Terahertz spectroscopy and density functional theory calculation. The results show that PBE-D2 shows more success than BLYP-D2 in prediction of THz absorption spectra. To compare their low-frequency modes, we adopted "vibrational character ID strips" proposed by Schmuttenmaer and coworkers [Journal of Physical Chemistry B, 117, 10444(2013)]. We found that the most intense THz absorption peaks of two compounds both involve severe distortion of their hydrogen bonding networks. Due to less rigid hydrogen bonding network in L-Asp, the side chain (carboxyl group) of L-Asp exhibits larger motions than that (carboxamide group) of L-Asn in low-frequency modes.

Intermolecular anharmonicity in molecular crystals: interplay between experimental low-frequency dynamics and quantum quasi-harmonic simulations of solid purine

A combined experimental spectroscopic and theoretical quasi-harmonic investigation of intermolecular anharmonicity in crystalline purine.

Resolving the Origins of Crystalline Anharmonicity Using Terahertz Time-Domain Spectroscopy and ab Initio Simulations

Anharmonicity has been shown to be an important piece of the fundamental framework that dictates numerous observable phenomena. In particular, anharmonicity is the driving force of vibrational relaxation processes, mechanisms that are integral to the proper function of numerous chemical processes. However, elucidating its origins has proven difficult due to experimental and theoretical challenges, specifically related to separating the anharmonic contributions from other unrelated effects. While no one technique is particularly suited for providing a complete picture of anharmonicity, by combining multiple complementary methods such a characterization can be made. In this study the role of individual atomic interactions on the anharmonic properties of crystalline purine, the building block of many DNA and RNA nucleobases, is studied by experimental terahertz time-domain spectroscopy and first-principles density functional theory (DFT) and ab initio molecular dynamics simulations (AIMD). In particular, the detailed vibrational information provided by the DFT calculations is used to interpret the atomic origins of anharmonic-related effects as determined by the AIMD calculations, which are in good agreement with the experimental data. The results highlight that anharmonicity is especially pronounced in the intermolecular interactions, particularly along the amine hydrogen bond coordinate, and yields valuable insight into what is similarly observed complex biosystems and crystalline solids.

Intramolecular vibrations in low-frequency normal modes of amino acids: L-alanine in the neat solid state

This paper presents a theoretical analysis of the low-frequency phonons of L-alanine by using the solid-state density functional theory at the Γ point. We are particularly interested in the intramolecular vibrations accessing low-frequency phonons via harmonic coupling with intermolecular vibrations. A new mode-analysis method is introduced to quantify the vibrational characteristics of such intramolecular vibrations. We find that the torsional motions of COO(-) are involved in low-frequency phonons, although COO(-) is conventionally assumed to undergo localized torsion. We also find the broad distributions of intramolecular vibrations relevant to important functional groups of amino acids, e.g., the COO(-) and NH3(+) torsions, in the low-frequency phonons. The latter finding is illustrated by the concept of frequency distribution of vibrations. These findings may lead to immediate implications in other amino acid systems.

Intermolecular hydrogen bond energies in crystals evaluated using electron density properties: DFT computations with periodic boundary conditions

The hydrogen bond (H-bond) energies are evaluated for 18 molecular crystals with 28 moderate and strong O-H···O bonds using the approaches based on the electron density properties, which are derived from the B3LYP/6-311G** calculations with periodic boundary conditions. The approaches considered explore linear relationships between the local electronic kinetic G(b) and potential V(b) densities at the H···O bond critical point and the H-bond energy E(HB). Comparison of the computed E(HB) values with the experimental data and enthalpies evaluated using the empirical correlation of spectral and thermodynamic parameters (Iogansen, Spectrochim. Acta Part A 1999, 55, 1585) enables to estimate the accuracy and applicability limits of the approaches used. The V(b)-E(HB) approach overestimates the energy of moderate H-bonds (E(HB) < 60 kJ/mol) by ~20% and gives unreliably high energies for crystals with strong H-bonds. On the other hand, the G(b)-E(HB) approach affords reliable results for the crystals under consideration. The linear relationship between G(b) and E(HB) is basis set superposition error (BSSE) free and allows to estimate the H-bond energy without computing it by means of the supramolecular approach. Therefore, for the evaluation of H-bond energies in molecular crystals, the G(b) value can be recommended to be obtained from both density functional theory (DFT) computations with periodic boundary conditions and precise X-ray diffraction experiments.

Noncovalent interactions between modified cytosine and guanine DNA base pair mimics investigated by terahertz spectroscopy and solid-state density functional theory

Modified cytosine and guanine nucleobases cocrystallize in a hydrogen bonding configuration similar to that observed in native DNA. The noncovalent interactions binding these base pairs in the crystalline solid were investigated using terahertz (THz) spectroscopy and solid-state density functional theory (DFT). While stronger hydrogen bonding interactions are responsible for the general molecular orientations in the crystalline state, it is the weaker dipole-dipole and dispersion forces that determine the overall packing arrangement. The inclusion of dispersion interactions in the DFT calculations was found to be necessary to accurately simulate the unit cell structure and THz vibrational spectrum. Using properly modeled intermolecular potentials, the lattice vibrational motions of the cytosine and guanine derivatives were calculated. The vibrational characters of the modes exhibited by the DNA base pair mimic in the THz region were primarily rotational motions and are indicative of the energies and the nature of vibrations that would likely be observed between similar base pairs in DNA molecules.

H-bond network in amino acid cocrystals with H2O or H2O2. The DFT study of serine-H2O and serine-H2O2

The structure, IR spectrum, and H-bond network in the serine-H(2)O and serine-H(2)O(2) crystals were studied using DFT computations with periodic boundary conditions. Two different basis sets were used: the all-electron Gaussian-type orbital basis set and the plane wave basis set. Computed frequencies of the IR-active vibrations of the titled crystals are quite different in the range of 10-100 cm(-1). Harmonic approximation fails to reproduce IR active bands in the 2500-2800 frequency region of serine-H(2)O and serine-H(2)O(2). The bands around 2500 and 2700 cm(-1) do exist in the anharmonic IR spectra and are caused by the first overtone of the OH bending vibrations of H(2)O and a combination vibration of the symmetric and asymmetric bendings of H(2)O(2). The quantum-topological analysis of the crystalline electron density enables us to describe quantitatively the H-bond network. It is much more complex in the title crystals than in a serine crystal. Appearance of water leads to an increase of the energy of the amino acid-amino acid interactions, up to ~50 kJ/mol. The energy of the amino acid-water H-bonds is ~30 kJ/mol. The H(2)O/H(2)O(2) substitution does not change the H-bond network; however, the energy of the amino acid-H(2)O(2) contacts increases up to 60 kJ/mol. This is caused by the fact that H(2)O(2) is a much better proton donor than H(2)O in the title crystals.