Cost-effective composite methods for large-scale solid-state calculations

Unlocking Diabetic Acetone Vapor Detection by A Portable Metal-Organic Framework-Based Turn-On Optical Sensor Device

Despite exhaled human breath having enabled noninvasive diabetes diagnosis, selective acetone vapor detection by fluorescence approach in the diabetic range (1.8-3.5 ppm) remains a long-standing challenge. A set of water-resistant luminescent metal-organic framework (MOF)-based composites have been reported for detecting acetone vapor in the diabetic range with a limit of detection of 200 ppb. The luminescent materials possess the ability to selectively detect acetone vapor from a mixture comprising nitrogen, oxygen, carbon dioxide, water vapor, and alcohol vapor, which are prevalent in exhaled breath. It is noteworthy that this is the first luminescent MOF material capable of selectively detecting acetone vapor in the diabetic range via a turn-on mechanism. The material can be reused within a matter of minutes under ambient conditions. Industrially pertinent electrospun luminescent fibers are likewise fabricated alongside various luminescent films for selective detection of ultratrace quantities of acetone vapor present in the air. Ab initio theoretical calculations combined with in situ synchrotron-based dosing studies uncovered the material's remarkable hypersensitivity toward acetone vapor. Finally, a freshly designed prototype fluorescence-based portable optical sensor is utilized as a proof-of-concept for the rapid detection of acetone vapor within the diabetic range.

Frontiers of molecular crystal structure prediction for pharmaceuticals and functional organic materials

The reliability of organic molecular crystal structure prediction has improved tremendously in recent years. Crystal structure predictions for small, mostly rigid molecules are quickly becoming routine. Structure predictions for larger, highly flexible molecules are more challenging, but their crystal structures can also now be predicted with increasing rates of success. These advances are ushering in a new era where crystal structure prediction drives the experimental discovery of new solid forms. After briefly discussing the computational methods that enable successful crystal structure prediction, this perspective presents case studies from the literature that demonstrate how state-of-the-art crystal structure prediction can transform how scientists approach problems involving the organic solid state. Applications to pharmaceuticals, porous organic materials, photomechanical crystals, organic semi-conductors, and nuclear magnetic resonance crystallography are included. Finally, efforts to improve our understanding of which predicted crystal structures can actually be produced experimentally and other outstanding challenges are discussed.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

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.

Performance of Screened-Exchange Functionals for Band Gaps and Lattice Constants of Crystals

Kohn-Sham density functional theory is the most widely used method for electronic structure calculations of solid-state systems. The screened-exchange functionals developed following the influential work of Scuseria and co-workers in 2003-2006 have significantly improved the accuracy of the predictions of solid-state properties. This work assesses six screened-exchange density functionals for the prediction of 60 band gaps (database BG60) and 68 lattice constants (database LC68). The band gaps are calculated with both consistently calculated lattice constants and experimental lattice constants. Results for the nonlocal screened-exchange functionals are compared with those for six widely used or recently developed local functionals. The results show that all the screened-exchange functionals have smaller mean absolute errors (MAEs) than any of the local functionals. All the functionals except HLE17 overestimate (on average) the lattice constants, and M06-SX gives the best performance among the compared functionals, with a MAE of 0.051 Å. All the functionals underestimate (on average) the band gaps, and M06-SX outperforms all other functionals, with a MAE of 0.47 eV. M06-SX also has the lowest root-mean-squared error for both LC68 and BG60. For the subdatabases of BG60, M06-SX shows better performance for ionic crystals and systems with large band gaps, while HSE12s gives better results for semiconductors and systems with small band gaps. Overall, M06-SX shows the best performance for solid-state systems, followed by N12-SX and HSE12s. The best-performing local functionals are M06-L, revM06-L, and HLE17 for band gaps and M06-L and revM06-L for lattice constants. We found that M06-SX, revM06-L, and N12-SX not only are well optimized for a broad array of chemical properties but also have very good performance for the databases in this paper, making them well-suited for applications involving heterogeneous chemistry.

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.

In Situ Observation of the Structure of Crystallizing Magnesium Sulfate Heptahydrate Solutions with Terahertz Transmission Spectroscopy

Terahertz time-domain spectroscopy in a transmission geometry combined with visual analysis was used to investigate the crystallization process of MgSO₄ solution. Careful spectral analysis of both a feature at 1.6 THz and the overall magnitude of absorption allowed the extraction of information about the liquid phase before and during crystallization, aiding the investigation of solvation dynamics and the behavior of molecular species at phase boundaries. The method was reproducibly applied to a number of measurements on a series of solutions of three chosen concentrations at different temperatures. When increasing temperature at the end of the measurement, the dissolution of crystals was observed as well. The temperature-dependent absorption data of the semicrystalline systems were converted to the solvent concentrations using a recently developed method. Solutions of a series of concentrations were also investigated in the temperature range of 4-25 °C. The results were compared to the theoretical calculated values, and the consistent differences proved the existence of a hydration shell around the salt ions whose behavior is different from bulk water. Future work will focus on triggering nucleation at specific positions in order to study the very beginning of the crystallization process. MgSO₄ heptahydrate is used as a model system in this study, while the concept and the setup can be applied to other systems.

Defect Engineering in Metal-Organic Framework Nanocrystals: Implications for Mechanical Properties and Performance

The growth process of metal-organic framework (MOF) nanocrystals defines their properties and functions. However, defects may be prevalent during the crystallization of even seemingly perfect MOFs, such as zeolitic imidazolate framework-8 (ZIF-8), and yet direct probing of such structural defects has been challenging because of the lack of nanoscale techniques to locally examine individual nanocrystals. Here, we directly study local defects, such as missing linkers or metal vacancies, in ZIF-8 nano- and microcrystals with near-field IR nanospectroscopy combined with density functional theory calculations. We track the chemical changes during crystallization and show that structural defects like zinc cations that are bound to molecules of the reactant gradually disappear with ripening of the crystals, while dangling and missing linker defects prevail. The resulting defect-terminating groups or open-metal sites produce mechanical anisotropy and reduce the Young's modulus, as measured via tip force microscopy with nanoscale resolution and supported by theoretical modeling. However, these structural defects also open the door for defect engineering to tune the performance of ZIF-8 by offering additional adsorption sites for targeted catalytic reactions, chemical sensing, or gas capture.

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.

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.