Vibrational Circular Dichroism Spectroscopy of Chiral Molecular Crystals: Insights from Theory

Calculation of Vibrational Circular Dichroism Spectra Employing Nuclear Velocity Perturbation or Magnetic Field Perturbation Theory Using an Atomic-Orbital-Based Linear Response Approach

We present the implementation of an AO-based solver for vibrational circular dichroism (VCD) spectra calculations, employing nuclear velocity perturbation theory using velocity-gauge included in atomic orbitals or magnetic field perturbation theory using gauge-including atomic orbitals. The implementations are done within the Gaussian and plane waves framework in the CP2K package. The previously implemented approaches in CP2K are based on MO-solvers, which solve the Sternheimer equation for linear response. Our AO-solver implementations were validated against the MO-solver by performing VCD calculations for the R-enantiomer of mirtazapine. Additionally, we extended the AO-based solver implementation to nuclear displacement perturbation theory in order to calculate infrared absorption spectra. The AO-based solver produced spectra that matched exactly with the MO-based results, confirming the accuracy of the implementation. This allows for efficient calculations of vibrational properties, further extending the capabilities of CP2K for large molecular systems.

Exploring Naproxen Cocrystals Through Solid-State Vibrational Circular Dichroism

Vibrational circular dichroism (VCD) spectroscopy appears as a useful method for characterizing optically active substances in the solid state. This is particularly important for active pharmaceutical ingredients. However, measurement and interpretation of the spectra bring about many difficulties. To assess the experimental and computational methodologies, we explore an anti-inflammatory drug, naproxen. Infrared (IR) and VCD spectra of the pure compound and its cocrystals with alanine and proline were recorded, and the data were interpreted by quantum chemical simulations based on a cluster model and density functional theory. Although unpolarized IR spectroscopy can already distinguish pure ingredients from cocrystals or a mixture, the VCD technique is much more sensitive. For example, the naproxen carboxyl group strongly interacts with the zwitterionic alanine in the cocrystal via two strong hydrogen bonds, which results in a rather rigid structure crystallizing in the chiral P212121 Sohncke group and its VCD is relatively strong. In contrast, the d-proline and (S)-naproxen cocrystal (P21 group) involves a single hydrogen bond between the subunits, which together with a limited motion of the proline ring gives a weaker signal. Solid-state VCD spectroscopy thus appears useful for exploring composite crystal structures and interactions within them, including studies of pharmaceutical compounds.

Mid-infrared giant optical chirality induced by multipole degeneracy in a diamond metasurface

Optical chirality is a fundamental property of light that plays a crucial role in chiroptical spectroscopy, analytical chemistry, spin photonics, and astrobiology. However, giant optical chirality at the wavelength beyond 15 μm remains untapped, which hinders its various applications, such as ultrasensitive vibrational circular dichroism, in the mid-infrared region. Here we propose and theoretically present giant optical chirality with a, to our knowledge, record-high value of ∼330 at the wavelength beyond 15 μm in a diamond metasurface. This is enabled by tailoring the multipole degeneracy and leveraging the low optical loss of the diamond metasurface. This work provides a new method for boosting the optical chirality at the wavelength beyond 15 μm.

Harmonic Scale Factors of Fundamental Transitions for Dispersion-corrected Quantum Chemical Methods

This work provides a procedure and database for obtaining the vibrational frequency scale factors that align quantum chemically computed harmonic frequencies with experimental vibrational spectroscopic data. The database comprises 441 molecules of various sizes, from diatomics to the buckminsterfullerene C60. We provide scale factors for 27 dispersion-corrected methods, 24 of which are DF-Dn/B with DF=BLYP, PBE, B3LYP, PBE0, Dn=D3(BJ), D4, and B=6-31G, def2-SVP, def2-TZVP, and three of them are the 3c-family composite methods (HF-3c, PBEh-3c, and r2SCAN-3c). The two scale factors are derived for each method: the absolute scaling, minimizing the absolute deviation of the scaled harmonic frequency from the experimental value, and the relative scaling, which minimizes an analogous relative deviation. The absolute type of scaling is recommended for frequencies above 2000 cm-1, while the relative scaling is optimal for frequencies below 2000 cm-1.

Nanoscale chirality generated in zinc(II) orthophosphate clusters: evidence by vibrational circular dichroism

Layered zinc(II) hydroxides (LZH) intercalating the deprotonated forms of R-(-) or S-(+)-1,1'-binaphthyl-2,2'-diyl hydrogenphosphate (denoted as R- or S-BNDHPH, respectively) were prepared from Zn(NO3)2 at pH 5 and 60 °C by the mixing method. The obtained hybrid compounds (denoted as R- or S-BNDHP-/LZH, respectively) were heated from room temperature up to 800 °C under nitrogen atmosphere. According to the thermal gravimetric/differential thermal analysis measurements, hydroxyl groups were dehydrated at 270-400 °C, followed by the decomposition of organic components at 420-600 °C. X-ray diffraction patterns and scanning electron microscopy images indicated that the final products were a mixture of α-Zn3(PO4)2, ZnO crystals and non-crystalline zinc(II) orthophosphates. Vibrational circular dichroism (VCD) spectra were recorded before and after calcination. Before calcination, R- or S-BNDHP-/LZH exhibited VCD peaks assigned to intercalated R- or S-BNDHP-. The calcined products exhibited several VCD peaks in the range of 900-1200 cm-1, maintaining the mirror-image relationship between R-BNDHP-/LZH and S-BNDHP-/LZH used as starting materials. The observed peaks were assigned to the PO (symmetric), -POO-, and PO (asymmetric) stretching vibrations of the PO43- groups. According to theoretical simulations, the observed VCD activity can be rationalised in terms of vibrational coupling between two PO43- groups in a generated chiral zinc(II) orthophosphate cluster.

Vibrational circular dichroism of adenosine crystals

Adenosine is one of the building blocks of nucleic acids and other biologically important molecules. Spectroscopic methods have been among the most utilized techniques to study adenosine and its derivatives. However, most of them deal with adenosine in solution. Here, we present the first vibrational circular dichroism (VCD) spectroscopic study of adenosine crystals in solid state. Highly regular arrangement of adenosine molecules in a crystal resulted in a strongly enhanced supramolecular VCD signal originating from long-range coupling of vibrations. The data suggested that adenosine crystals, in contrast to guanosine ones, do not imbibe atmospheric water. Relatively large dimensions of the adenosine crystals resulted in scattering and substantial orientational artifacts affecting the spectra. Several strategies for tackling the artifacts have been proposed and tested. Atypical features in IR absorption spectra of crystalline adenosine (e.g., extremely low absorption in mid-IR spectral range) were observed and attributed to refractive properties of adenosine crystals.

Solid-state vibrational circular dichroism for pharmaceutical applications: Polymorphs and cocrystal of sofosbuvir

X-ray diffraction is a commonly used technique in the pharmaceutical industry for the determination of the atomic and molecular structure of crystals. However, it is costly, sometimes time-consuming, and it requires a considerable degree of expertise. Vibrational circular dichroism (VCD) spectroscopy resolves these limitations, while also exhibiting substantial sensitivity to subtle modifications in the conformation and molecular packaging in the solid state. This study showcases VCD's ability to differentiate between various crystal structures of the same molecule (polymorphs, cocrystals). We examined the most effective approach for producing high-quality spectra and unveiled the intricate link between structure and spectrum via quantum-chemical computations. We rigorously assessed, using alanine as a model compound, multiple experimental conditions on the resulting VCD spectra, with the aim of proposing an optimal and efficient procedure. The proposed approach, which yields reliable, reproducible, and artifact-free results with maximal signal-to-noise ratio, was then validated using a set comprising of three amino acids (serine, alanine, tyrosine), one hydroxy acid (tartaric acid), and a monosaccharide (ribose) to mimic active pharmaceutical components. Finally, the optimized approach was applied to distinguish three polymorphs of the antiviral drug sofosbuvir and its cocrystal with piperazine. Our results indicate that solid-state VCD is a prompt, cost-effective, and easy-to-use technique to identify crystal structures, demonstrating potential for application in pharmaceuticals. We also adapted the cluster and transfer approach to calculate the spectral properties of molecules in a periodic crystal environment. Our findings demonstrate that this approach reliably produces solid-state VCD spectra of model compounds. Although for large molecules with many atoms per unit cell, such as sofosbuvir, this approach has to be simplified and provides only a qualitative match, spectral calculations, and energy analysis helped us to decipher the observed differences in the experimental spectra of sofosbuvir.

The Phenomenon of Vibrational Circular Dichroism Enhancement: A Systematic Survey of Literature Data

While the intensity of vibrational circular dichroism (VCD) signals is commonly 104-105 times smaller than that of corresponding IR signals, several kinds of systems display enhanced VCD spectra with g-values (VCD/IR intensity ratio) above 10-3 and even reaching 5 × 10-2 in some exceptional cases. These systems include transition metal and lanthanide complexes, protein and peptide fibrils, short oligopeptide gels, crystalline compounds, gels and solution aggregates of organic compounds. We review the literature on VCD enhancement, focusing on collecting and analyzing data on enhanced g-values. Special attention is given to the mechanisms proposed to produce these effects.

Stereochemistry of natural products from vibrational circular dichroism

Secondary metabolites from land and marine (micro)organisms have been at the focus of the drug discovery process for many years. One of the reasons for this success is nature's incredible ability to create intricate molecular scaffolds. Such structural richness, however, makes the structural elucidation, and the absolute configuration assignment in particular, a challenging process. Vibrational circular dichroism (VCD) has emerged as one of the most reliable and versatile methods to unambiguously assign both the absolute configuration and conformations of chiral molecules in solution. Although VCD is no longer a curiosity in the field of molecular spectroscopy after 50 years since its first report, it is still underutilized by natural product chemists worldwide for varying reasons. Herein, we highlight the evolution of the application of VCD to natural product chemistry, focusing on its strengths as well as points that still need improvement. General guidelines for the correct application of VCD to stereochemical studies are also provided.

After 50 Years of Vibrational Circular Dichroism Spectroscopy: Challenges and Opportunities of Increasingly Accurate and Complex Experiments and Computations

VCD research continues to thrive, driven by ongoing experimental and theoretical advances. Modern studies deal with increasingly complex samples featuring weak intermolecular interactions and shallow potential energy surfaces. Likewise, the combination of VCD measurements with, for instance, cryo-spectroscopic techniques has significantly increased their sensitivity. The extent to which such modern measurements enhance the informative value of VCD depends significantly on the quality of the theoretical models, which must adequately account for anharmonicity, solvation and molecular dynamics. We herein discuss how experimental advancements engage in a stimulating interplay with recent theoretical developments, pursuing either the static or the dynamic computational route. Both paths have their own strengths and limitations, each addressing fundamentally different problems. We give an outlook on future challenges of VCD research, including the possibility to combine static and dynamic approaches to obtain a full picture of the sample.

Optical Phenomena in Molecule-Based Magnetic Materials

Since the last century, we have witnessed the development of molecular magnetism which deals with magnetic materials based on molecular species, i.e., organic radicals and metal complexes. Among them, the broadest attention was devoted to molecule-based ferro-/ferrimagnets, spin transition materials, including those exploring electron transfer, molecular nanomagnets, such as single-molecule magnets (SMMs), molecular qubits, and stimuli-responsive magnetic materials. Their physical properties open the application horizons in sensors, data storage, spintronics, and quantum computation. It was found that various optical phenomena, such as thermochromism, photoswitching of magnetic and optical characteristics, luminescence, nonlinear optical and chiroptical effects, as well as optical responsivity to external stimuli, can be implemented into molecule-based magnetic materials. Moreover, the fruitful interactions of these optical effects with magnetism in molecule-based materials can provide new physical cross-effects and multifunctionality, enriching the applications in optical, electronic, and magnetic devices. This Review aims to show the scope of optical phenomena generated in molecule-based magnetic materials, including the recent advances in such areas as high-temperature photomagnetism, optical thermometry utilizing SMMs, optical addressability of molecular qubits, magneto-chiral dichroism, and opto-magneto-electric multifunctionality. These findings are discussed in the context of the types of optical phenomena accessible for various classes of molecule-based magnetic materials.

Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases

Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.

VCD spectra of chiral naphthalene-1-carboxamides in the solid-state

The VCD spectra of chiral 2,3-dihydro-1H-benzo[de]isoquinolin-1-one (8-substituted naphthalene-1-carboxamide, BIQ) were studied in KBr pellets. The X-ray diffractometry revealed that the Me, Ph, and pClPh BIQs crystalize in the monoclinic P21, while nBu, pMePh, and oMeOPh BIQs in the orthorhombic P212121 space group. Only the Me-BIQ crystal exhibits the presence of cyclic amide dimers, while the others contain chains of the amid group hydrogen bonds. For all BIQs, except pMePh, the most intense IR band in the 1750-1550 cm-1 region is located at ca. 1680 cm-1 and is accompanied by two weak ones at ca. 1618 and 1590 cm-1. For the pMePh derivative, four almost equally intense IR bands at 1662, 1639, 1614, and 1588 cm-1 are observed. This region of the IR spectra of BIQs, but pMePh, is well reproduced by calculations based on BIQ monomers. On the other hand, the complex IR pattern of pMePh is computationally reproduced when larger crystal fragments, like octamers, are considered. Registration of the VCD spectra enabled recognizing the complexity of IR contours at ca. 1680 cm-1 by the corresponding VCD motives. For (i) Me, Ph and pClPh (R)-enantiomers, two (+)(-) bands were distinguished and for (ii) nBu and pMePh ones, one VCD band with right-side asymmetry was found. For (iii) oMeOPh the VCD pattern cannot be unambiguously assigned. Thus, the VCD spectra in the ν(C=O) range diverse the studied compounds. Among the set of molecules, pMePh has exceptional crystal geometry. Therefore, its most intense ν(C=O) band position and shape can be connected with the geometry of the hydrogen bonds, interactions, and crystal packing. Interpretation of the VCD spectra is based on linear and packed BIQ octamers. This cluster model can reproduce the main features of the solid-state VCD of BIQs.