In Situ Cocrystallization via Spray Drying with Polymer as a Strategy to Prevent Cocrystal Dissociation.

The aim of the present study was to investigate how different polymers affect the dissociation of cocrystals prepared by co-spray-drying active pharmaceutical ingredient (API), coformer, and polymer. Diclofenac acid-l-proline cocrystal (DPCC) was selected in this study as a model cocrystal due to its previously reported poor physical stability in a high-humidity environment. Polymers investigated include polyvinylpyrrolidone (PVP), poly(1-vinylpyrrolidone-co-vinyl acetate) (PVPVA), hydroxypropyl methyl cellulose, hydroxypropylmethylcellulose acetate succinate, ethyl cellulose, and Eudragit L-100. Terahertz Raman spectroscopy (THz Raman) and powder X-ray diffraction (PXRD) were used to monitor the cocrystal dissociation rate in a high-humidity environment. A Raman probe was used in situ to monitor the extent of the dissociation of DPCC and DPCC in crystalline solid dispersions (CSDs) with polymer when exposed to pH 6.8 phosphate buffer and water. The solubility of DPCC and solid dispersions of DPCC in pH 6.8 phosphate buffer and water was also measured. The dissociation of DPCC was water-mediated, and more than 60% of DPCC dissociated in 18 h at 40 °C and 95% RH. Interestingly, the physical stability of the cocrystal was effectively improved by producing CSDs with polymers. The inclusion of just 1 wt % polymer in a CSD with DPCC protected the cocrystal from dissociation over 18 h under the same conditions. Furthermore, the CSD with PVPVA was still partially stable, and the CSD with PVP was stable (undissociated) after 7 days. The superior stability of DPCC in CSDs with PVP and PVPVA was also demonstrated when systems were exposed to water or pH 6.8 phosphate buffer and resulted in higher dynamic solubility of the CSDs compared to DPCC alone. The improvement in physical stability of the cocrystal in CSDs was thought to be due to an efficient mixing between polymer and cocrystal at the molecular level provided by spray drying and in situ gelling of polymer. It is hypothesized that polymer chains could undergo gelling in situ and form a physical barrier, preventing cocrystal interaction with water, which contributes to slowing down the water-mediated dissociation.

Low-Frequency Vibrational Spectroscopy and Quantum Mechanical Simulations of the Crystalline Polymorphs of the Antiviral Drug Ribavirin.

Crystal polymorphism is a common phenomenon in pharmaceutical solids and a critical issue when considering the formulation of therapeutics since multiple polymorphs may form during drug manufacturing. Low-frequency vibrational spectroscopy is sensitive to polymorphic content, and in this work, terahertz time-domain spectroscopy and low-frequency Raman spectroscopy were utilized in the study of crystalline ribavirin, a widely applicable antiviral. Characteristic spectra with numerous peaks in the sub-200 cm-1 region were obtained of the more common polymorph of ribavirin (Form II). Solid-state density functional theory (ss-DFT) simulations were then used to optimize the crystal structure of this polymorph and calculate the frequencies and spectral intensities of the lattice vibrations in the low-frequency region. The near-harmonic thermal behavior of the sample with cooling enabled excellent agreement between experiment and theory to be achieved, emphasizing the quality of the applied model, and the observed spectral peaks could be assigned to specific atomic motions in the solid. Form I and Form II polymorphs of ribavirin were both investigated with ss-DFT to understand the different aspects governing the relative stabilities of these solids. The ss-DFT simulations of the polymorph energies revealed that Form II is more stable at all temperatures due to a stronger cohesive energy than Form I; however, ribavirin in Form I has a significantly lower conformational energy. The finding of monotropism appears to conflict with the reported enantiotropism of the ribavirin polymorphs but ultimately confirms that crystal defects in the real samples greatly affect the thermodynamic relationship of the crystals.

New Insights into the Crystallographic Disorder in the Polymorphic Forms of Aspirin from Low-Frequency Vibrational Analysis.

Terahertz time-domain spectroscopy (THz-TDS) is applied to two polymorphs of acetylsalicylic acid (aspirin), and the experimental spectra are compared to lattice dynamical calculations using high accuracy density functional theory. The calculations confirm that forms I and II have very close energetic and thermodynamic properties and also that they show similar spectral features in the far-infrared region, reflecting the high degree of similarity in their crystal structures. Unique vibrational modes are identified for each polymorph which allow them to be distinguished using THz-TDS measurements. The observation of spectral features attributable to both polymorphic forms in a single sample, however, provides further evidence to support the hypothesis that crystalline aspirin typically comprises intergrown domains of forms I and II. Differences observed in the baseline of the measured THz-TDS spectra indicate a greater degree of structural disorder in the samples of form II. Calculated Gibbs free-energy curves show a turning point at 75 K, inferring that form II is expected to be more stable than form I above this temperature as a result of its greater vibrational entropy. The calculations do not account for any differences in configurational entropy that may arise from expected structural defects. Further computational work on these structures, such as ab initio molecular dynamics, would be very useful to further explore this perspective. Here, aspirin is a model system to show how the additional insight from the low-frequency vibrational information complements the structural data and allows for quantitative thermodynamic information of pharmaceutical polymorphs to be extracted. The methodology is directly applicable to other polymorphic systems.

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.

Quantitative Analysis of Minium and Vermilion Mixtures Using Low-Frequency Vibrational Spectroscopy.

Low-frequency vibrational spectroscopy offers a compelling solution for the nondestructive and noninvasive study of pigments in historical artifacts by revealing the characteristic sub-200 cm-1 spectral features of component materials. The techniques of terahertz time-domain spectroscopy (THz-TDS) and low-frequency Raman spectroscopy (LFRS) are complementary approaches to accessing this spectral region and are valuable tools for artifact identification, conservation, and restoration. In this investigation of historical pigments, pure and mixed samples of minium (Pb3O4) and vermilion (HgS) were studied using a combination of THz-TDS and LFRS experiments to determine the limits of detection (LOD) and quantitation (LOQ) for each compound with both methods. The measurements were also supported using solid-state density functional theory simulations of the pigment structures and vibrations, enabling spectral peaks to be assigned to specific atomic motions in these solids. The THz-TDS LOD was found to be similar for both minium and vermilion at 6% by mass on average. In comparison, LFRS was found to be more sensitive to both pigments, particularly to the presence of vermilion with an LFRS LOD of 0.2%. These results demonstrate that low-frequency vibrational spectroscopy can be used for successful quantitative analysis of pigment mixtures and provide reliable new data for use in heritage science.

Identifying and explaining vibrational modes of quinacridones via temperature-resolved terahertz spectroscopy: absorption experiments and solid-state density functional theory simulations.

Quinacridone and its substituted analogs are pigments widely used in art and industry. The temperature dependence of the crystal structures of two quinacridone polymorphs (ß and γ), along with the common variant 2,9-dimethylquinacridone, were investigated using powder X-ray diffraction and terahertz spectroscopy. These were then compared with solid-state density functional theory simulations of both structures and vibrations. X-ray patterns were collected at eight temperatures in the range 13-298 K and terahertz spectra at fifteen temperatures in the range 20-300 K. Simulations were at absolute zero and at appropriate expansions to model room temperature. It was found that some of the powder X-ray diffraction features in only ß-quinacridone (15.7°, 19.7° and 31.2° at 13 K) underwent anomalous shifting with temperature change. We attribute this to the unique coplanar hydrogen bonding pattern of ß-quinacridone compared to the other solids, with the unusual diffraction peaks originating from crystallographic planes perpendicular to the a axis intermolecular hydrogen bonds. This observation coincides with a contraction of the a axis with heating and results from its relatively weak N-HO hydrogen bonds and significant C-HH-C repulsions. Associated with this anomalous contraction, for ß-quinacridone only spectral peaks are seen to increase in energy with temperature.

Simulation and Assignment of the Terahertz Vibrational Spectra of Enalapril Maleate Cocrystal Polymorphs.

The identification of crystalline drug polymorphs using terahertz vibrational spectroscopy is a powerful approach for the nondestructive and noninvasive characterization of solid-state pharmaceuticals. However, a complete understanding of the terahertz spectra of molecular solids is challenging to obtain because of the complex nature of the low-frequency vibrational motions found in the sub-3 THz (sub-100 cm-1) range. Unambiguous assignments of the observed spectral features can be achieved through quantum mechanical solid-state simulations of crystal structures and lattice vibrations utilizing the periodic boundary condition approach. The terahertz spectra of two polymorphs of enalapril maleate are presented here to demonstrate that even large pharmaceuticals can be successfully modeled using solid-state density functional theory, including cocrystalline solids comprised of multiple distinct species. These simulations enable spectral assignments to be made, but also provide insights into the conformational and cohesion energies that contribute to the polymorph stabilities. The results reveal that the Form II polymorph of enalapril maleate is the more stable of the two under ambient conditions, and that this stability is driven by a greater intermolecular cohesion energy as compared to Form I.

Anomalous Temperature Dependence of the Lowest-Frequency Lattice Vibration in Crystalline γ-Aminobutyric Acid.

Crystalline γ-aminobutyric acid (GABA) exhibits unusual thermal behavior in a low-frequency lattice vibration that occurs at 37.2 cm-1 at 290 K but decreases dramatically by 34.0% when the sample is cooled to 78 K. Lattice vibrations in molecular crystals are indicators of intermolecular force characteristics, and the extraordinary temperature sensitivity of this vibration offers new insight into the local environment within the solid. Solid-state density functional theory simulations of the GABA crystal have found this anomalous frequency shift is based in unexpected differences in the strengths of the intermolecular hydrogen bonds that are cursorily the same. This was accomplished through mapping of the potential energy surfaces governing the terahertz-frequency motions of the GABA solid and use of the quasi-harmonic approximation to model the response of all the lattice vibrations to temperature-induced unit cell volume changes brought about through the anharmonic character of the intermolecular interactions. The analysis reveals that the vibration in question is rotational in nature and involves the significant distortion of a specific weak intermolecular N-H···O hydrogen bond in the crystal that results in its unique thermal response.

Terahertz Spectroscopy and Quantum Mechanical Simulations of Crystalline Copper-Containing Historical Pigments.

Terahertz spectroscopy, a noninvasive and nondestructive analytical technique used in art conservation and restoration, can provide compelling data concerning the composition and condition of culturally valuable and historical objects. Terahertz spectral databases of modern and ancient artists' pigments exist but lack explanations for the origins of the unique spectral features. Solid-state density functional theory simulations can provide insight into the molecular and intermolecular forces that dominate the observed absorption features as well as reveal deviations from simple harmonic vibrational behavior that can complicate these spectra. The characteristic terahertz spectra of solid azurite, malachite, and verdigris are presented here, along with simulations of their crystalline structures and sub-3.0 THz lattice vibrations. The powerful combination of theory and experiment enables unambiguous spectral assignment of these complex materials and highlights the challenges that anharmonic peak broadening in organic-containing materials may present in the construction of reference pigment databases.

Optical Properties of Meloxicam in the Far-Infrared Spectral Region.

One of the most commonly used nonsteroidal anti-inflammatory active pharmaceutical ingredient called Meloxicam has been characterized spectroscopically both by Terahertz (THz) time domain spectroscopy (THz-TDS) and by Fourier Transform Infrared (FTIR) spectroscopy in far-IR regions of electromagnetic spectrum; 0.2 THz to 20 THz. While many relatively sharp features are observed in the far-IR range between 2 THz to 20 THz as expected for being an organic substance, very distinct and relatively strong absorption bands are also observed at 1.00, 1.66, 2.07 and 2.57 THz in the THz range. These well separated, defined, and fairly strong spectral features can be used for discrimination and quantification of Meloxicam in drug analysis. Frequency dependent refractive index of the drug was determined in a range of 0.2 THz and 2.7 THz, where an almost constant index was observed with an average index of 1.75. Powder XRD, and solid-state Density Functional Theory (SS-DFT) calculations were utilized to determine the crystalline form of the Meloxicam sample in its enolic crystalline form. Single molecule DFT calculations were also performed in all four possible structures of Meloxicam. In addition, the capability of THz waves transmission through common packaging materials is demonstrated for possibility of future on-site analysis. The results suggest that drug analysis will be possible to perform not only at every stage of manufacturing without destruction but also directly at the shelf of a market after development of portable THz technologies.

Concomitant polymorphism and the martensitic-like transformation of an organic crystal.

Crystalline polymorphism is a phenomenon that occurs in many molecular solids, resulting in a diverse range of possible bulk structures. Temperature and pressure can often be used to thermodynamically control which crystal form is preferred, and the associated transitions between polymorphic phases are often discontinuous and complete. N-Methyl-4-carboxypyridinium chloride is a solid that undergoes an apparent continuous temperature-dependent phase transition from an orthorhombic to a monoclinic polymorph. However, a hybrid characterization approach using single-crystal X-ray diffraction, terahertz time-domain spectroscopy, and solid-state density functional theory reveals the transformation to be actually a slowly changing ratio of the two discrete polymorphic forms. The potential energy surface of this process can be directly accessed using terahertz radiation, and the data show that a very low barrier (43.3 J mol-1) exists along the polymorph transformation coordinate.

Quantification of cation-anion interactions in crystalline monopotassium and monosodium glutamate salts.

Crystalline salt compounds composed of metal cations and organic anions are becoming increasingly popular in a number of fields, including the pharmaceutical and food industries, where such formulations can lead to increased product solubility. The origins of these effects are often in the interactions between the individual components in the crystals, and understanding these forces is paramount for the design and utilisation of such materials. Monosodium glutamate monohydrate and monopotassium glutamate monohydrate are two solids that form significantly different structures with correspondingly dissimilar dynamics, while their chemistry only differs in cation identity. Crystals of each were characterised experimentally with single-crystal X-ray diffraction and terahertz time-domain spectroscopy and theoretically using solid-state density functional theory simulations, in order to explain the observed differences in their bulk properties. Specifically, crystal orbital overlap and Hamiltonian population analyses were performed to examine the role that the individual interactions between the cation and anion played in the solid-state structures and the overall energetic profiles of these materials.

The significance of the amorphous potential energy landscape for dictating glassy dynamics and driving solid-state crystallisation.

The fundamental origins surrounding the dynamics of disordered solids near their characteristic glass transitions continue to be fiercely debated, even though a vast number of materials can form amorphous solids, including small-molecule organic, inorganic, covalent, metallic, and even large biological systems. The glass-transition temperature, Tg, can be readily detected by a diverse set of techniques, but given that these measurement modalities probe vastly different processes, there has been significant debate regarding the question of why Tg can be detected across all of them. Here we show clear experimental and computational evidence in support of a theory that proposes that the shape and structure of the potential-energy surface (PES) is the fundamental factor underlying the glass-transition processes, regardless of the frequency that experimental methods probe. Whilst this has been proposed previously, we demonstrate, using ab initio molecular-dynamics (AIMD) simulations, that it is of critical importance to carefully consider the complete PES - both the intra-molecular and inter-molecular features - in order to fully understand the entire range of atomic-dynamical processes in disordered solids. Finally, we show that it is possible to utilise this dependence to directly manipulate and harness amorphous dynamics in order to control the behaviour of such solids by using high-powered terahertz pulses to induce crystallisation and preferential crystal-polymorph growth in glasses. Combined, these findings provide compelling evidence that the PES landscape, and the corresponding energy barriers, are the ultimate controlling feature behind the atomic and molecular dynamics of disordered solids, regardless of the frequency at which they occur.

Comparison of Intermolecular Forces in Anhydrous Sorbitol and Solvent Cocrystals.

The hygroscopicity of solid sorbitol is important for its utilization as a sweetener in the pharmaceutical and food industries. The molecular foundations of sorbitol hydration characteristics are explored here using two solvated cocrystals, sorbitol-water and sorbitol-pyridine. In this work, solid-state density functional theory and terahertz time-domain spectroscopy were used to evaluate the relative stabilities of these cocrystals as compared to anhydrous sorbitol in terms of conformational and cohesive energies. The modification of the hydrogen-bonding network in crystalline sorbitol by solvent molecules gives new insight into the origins of the notable stability of sorbitol-water as compared to similar solids such as mannitol-water. In particular, the energy analysis reveals that the relative instability of the mannitol hydrate is based primarily in the lack of water-water interactions which provide considerable stabilization in the sorbitol-water crystal.

Distinguishing Quinacridone Pigments via Terahertz Spectroscopy: Absorption Experiments and Solid-State Density Functional Theory Simulations.

Through a combined experimental and theoretical investigation we determine that the fundamental modes of three quinacridones fall in the terahertz spectral range (1-10 THz, ∼30-300 cm-1). In each spectrum the terahertz resonances correspond to wagging, rocking, or twisting of the quinacridone rings, with the most intense absorption being an in-plane rocking vibration of the carbonyl oxygens. In spite of these spectral similarities, we demonstrate that terahertz measurements readily differentiate ß-quinacridone, γ-quinacridone, and 2,9-dimethylquinacridone. The spectrum of ß-quinacridone has a group of closely spaced modes at ∼4 THz, whereas in contrast the spectrum of γ-quinacridone displays a widely spaced series of modes spread over the range ∼1-5 THz. Both of these have the strongest mode at ∼9 THz, whereas in contrast 2,9-dimethylquinacridone exhibits the strongest mode at ∼7 THz. Because quinacridones are the basis of widely used synthetic pigments of relatively recent origin, our findings offer promising applications in the identification and dating of modern art.

The crucial role of water in shaping low-barrier hydrogen bonds.

Low-barrier hydrogen bonds (LBHBs) are key components in a range of chemical processes, often appearing in metal-mediated catalytic applications. The formation of LBHBs has generally been attributed to the particular geometric arrangement of the surrounding atoms and molecules, yet few specific atomic-level details have emerged. In this study, several metal maleate tetrahydrates have been studied using a combination of solid-state density functional theory, terahertz spectroscopy, and X-ray diffraction to evaluate the significance of both water molecules and metal cations in guiding LBHB formation and function. The findings reveal the assumption that metal identity is of paramount importance to be incomplete, and that the metal cation does not directly influence the LBHB in the maleate ligand. Rather, the characteristics are regulated by water molecule positioning, asserting the critical role of water in governing LBHBs and providing new insight into their formation mechanisms.

Origins of the Relative Stabilities of Anhydrous and Hydrated d-Mannitol Crystals.

d-Mannitol, a common bulking agent used in the pharmaceutical and food industries, is known to crystallize in three anhydrous polymorphs (α, ß, and δ) and a hemihydrate form. In this work, solid-state density functional theory simulations and experimental terahertz spectroscopy were used to evaluate the stabilities, thermodynamic relationships, and spectral signatures of these various solids. The thermodynamic analyses of the anhydrous polymorphs identified enantiotropic transitions for ß-δ at approximately room temperature and α-δ near the α melting point, but showed α and ß to be monotropically related. Computational refinement of the d-mannitol hemihydrate crystal structure and comparison of its Gibbs free energy to the anhydrous forms showed the hemihydrate to be stable only at low temperatures (<253 K), in agreement with experiment, and reveals its conversion to the δ polymorph upon dehydration. Quantification of the intramolecular and intermolecular energies within these solids highlights the competition that exists between molecular conformation and cohesion energies and how a balance of these forces dictates the observed behaviors.

Uncovering the Terahertz Spectrum of Copper Sulfate Pentahydrate.

Terahertz vibrational spectroscopy has evolved into a powerful tool for the detection and characterization of transition metal sulfate compounds, specifically for its ability to differentiate between various hydrated forms with high specificity. Copper(II) sulfate is one such system where multiple crystalline hydrates have had their terahertz spectra fully assigned, and the unique spectral fingerprints of the forms allows for characterization of multicomponent systems with relative ease. Yet the most commonly occurring form, copper(II) sulfate pentahydrate (CuSO4·5H2O), has proven elusive due to the presence of a broad absorption across much of the terahertz region, making the unambiguous identification of its spectral signature difficult. Here, it is shown that the sub-100 cm(-1) spectrum of CuSO4·5H2O is obscured by absorption from adsorbed water and that controlled drying reveals sharp underlying features. The crystalline composition of the samples was monitored in parallel by X-ray diffraction as a function of drying time, supporting the spectroscopic results. Finally, the terahertz spectrum of CuSO4·5H2O was fully assigned using solid-state density functional theory simulations, helping attribute the additional absorptions that appear after excessive drying to formation of CuSO4·3H2O.

Evaluation of Range-Corrected Density Functionals for the Simulation of Pyridinium-Containing Molecular Crystals.

The problem of nonlocal interactions in density functional theory calculations has in part been mitigated by the introduction of range-corrected functional methods. While promising solutions, the continued evaluation of range corrections in the structural simulations of complex molecular crystals is required to judge their efficacy in challenging chemical environments. Here, three pyridinium-based crystals, exhibiting a wide range of intramolecular and intermolecular interactions, are used as benchmark systems for gauging the accuracy of several range-corrected density functional techniques. The computational results are compared to low-temperature experimental single-crystal X-ray diffraction and terahertz spectroscopic measurements, enabling the direct assessment of range correction in the accurate simulation of the potential energy surface minima and curvatures. Ultimately, the simultaneous treatment of both short- and long-range effects by the ωB97-X functional was found to be central to its rank as the top performer in reproducing the complex array of forces that occur in the studied pyridinium solids. These results demonstrate that while long-range corrections are the most commonly implemented range-dependent improvements to density functionals, short-range corrections are vital for the accurate reproduction of forces that rapidly diminish with distance, such as quadrupole-quadrupole interactions.

Examination of l-Glutamic Acid Polymorphs by Solid-State Density Functional Theory and Terahertz Spectroscopy.

The ability of l-glutamic acid to crystallize in two different forms has long been the subject of study due to its commercial importance. While a solvent-mediated phase transformation between the α and ß polymorphs is the prevailing theory, recent reports indicate a thermal solid-solid transformation between the two may be possible. However, determining accurate thermodynamic stabilities of these crystals has been challenging. Here new low-temperature single-crystal X-ray diffraction data coupled to solid-state density functional theory simulations have enabled a detailed description to be achieved for the energetic parameters governing the stabilization of the two l-glutamic acid solids. The temperature-dependent Gibbs free-energy curves show that α-glutamic acid is the preferred form at low temperatures (<222 K) and the ß form is most stable at ambient temperatures. Terahertz time-domain spectroscopy was utilized to evaluate the quality of the intermolecular force modeling as well as to provide characteristic low-frequency spectral data that can be used for quantification of polymorph mixtures or crystal growth monitoring.