Optical Rotation Calculated with Time-Dependent Density Functional Theory: The OR45 Benchmark

First Principles Simulations of Optical Rotation of Chiral Molecular Crystals

In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT-PBC). Calculations are compared with experimental measurements and show semi-quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X-ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.

Helicene Appended Benzothiadiazoles as Chiral Emitters

A homologous series of 4,7-bis(aryl) substituted benzothiadiazole (BTD) compounds, containing the helicenic derivatives bis([4]helicene), bis([5]helicene) and bis([6]helicene), have been prepared upon a double Suzuki coupling between 3,6-bis(pinacolyl-borane)-BTD and the corresponding bromo-aryl precursors. The single crystal X-ray structure of the bis([4]helicene) compound shows the existence of both helicities (M) and (P) on the same molecule. All the compounds of the series are highly emissive in solution, with quantum yields of the emission ranging from 50 to 91 %. The enantiopure compounds (M,M) and (P,P) for the BTD-bis([6]helicene) have been prepared from the corresponding enantiopure 2-bromo-[6]helicene precursors. Their chiroptical properties have been investigated in correlation with density functional theory (DFT) calculations, which allowed to confidently assign the absolute configuration of the helicene arms and to characterize the different electronic transitions, including the low energy charge transfer excitation from helicenes to BTD. The enantiomerically pure fluorophores (M,M)- and (P,P)-BTD-bis([6]helicene), which exist in solution as two main conformers, according to the DFT calculations, show CPL activity in solution, with glum factors of ≈1.7×10-3 at λem=525 nm, and also in the solid state, with glum factors of ≈1.2×10-3 in spite of the strong decrease of the quantum efficiency.

A Perspective on the Simulation of Electronic Circular Dichroism and Circularly Polarized Luminescence Spectra in Chiral Solid Materials

Chiral materials have shown tremendous potential for many technological applications, such as optoelectronics, sensing, magnetism, information technology, and imaging. Characterization of these materials is mostly based on chiroptical spectroscopies, such as electronic circular dichroism (ECD) and circularly polarized luminescence (CPL). These experimental measurements would greatly benefit from theoretical simulations for interpretation of the spectra as well as predictions on new materials. While ECD and CPL simulations are well established for molecular systems, they are not for materials. In this Perspective, we describe the theoretical quantities necessary to simulate ECD and CPL spectra in oriented systems. Then, we discuss the approximate strategies currently used to perform these calculations, what computational machinery is already available to develop more general approaches, and some of the open challenges for the simulation of ECD and CPL spectra in solid materials. When methods that are as reliable and computationally efficient as those for molecules are developed, these simulations will provide invaluable insight and guidance for the rational design of optically active materials.

Acid/Base-Triggered Photophysical and Chiroptical Switching in a Series of Helicenoid Compounds

A series of molecules that possess two quinolines, benzoquinolines, or phenanthrolines connected in a chiral fashion by a biaryl junction along with their water-soluble derivatives was developed and characterized. The influence of the structure on the basicity of the nitrogen atoms in two heterocycles was examined and the photophysical and chiroptical switching activity of the compounds upon protonation was studied both experimentally and computationally. The results demonstrated that changes in the electronic structure of the protonated vs. neutral species, promoting a bathochromic shift of dominant electronic transitions and alternation of their character from π-to-π* to charge-transfer-type, when additionally accompanied by the high structural flexibility of a system, leading to changes in conformational preferences upon proton binding, produce particularly pronounced modifications of the spectral properties in acidic medium. The latter combined with reversibility of the read-out make some of the molecules in this series very promising multifunctional pH probes.

On the choice of coordinate origin in length gauge optical rotation calculations

In this work, we explore the issue of origin dependence in optical rotation (OR) calculations in the length dipole gauge (LG) using standard approximate methods belonging to density functional theory (DFT) and coupled cluster (CC) theory. We use the origin-invariant LG approach, LG(OI), that we recently proposed as reference for the calculations, and we study whether a proper choice of coordinate origin and molecular orientation can be made such that diagonal elements of the LG-OR tensor match those of the LG(OI) tensor. Using a numerical search algorithm, we show that multiple spatial orientations can be found where the LG and LG(OI) results match. However, a simple analytical procedure provides a spatial orientation where the origin of the coordinate system is close to the center of mass of the molecule. At the same time, we also show that putting the origin at the center of mass is not an ideal choice for every molecule (relative errors in the OR up to 70% can be obtained in out test set). Finally, we show that the choice of coordinate origin based on the analytical procedure is transferable across different methods and it is superior to putting the origin in the center of mass or center of nuclear charge. This is important because the LG(OI) approach is trivial to implement for DFT, but not necessarily for nonvariational methods in the CC family. Therefore, one can determine an optimal coordinate origin at DFT level and use it for standard LG-CC response calculations.

X-ray and Optical Circular Dichroism as Local and Global Ultrafast Chiral Probes of [12]Helicene Racemization

Chirality is a fundamental molecular property that plays a crucial role in biophysics and drug design. Optical circular dichroism (OCD) is a well-established chiral spectroscopic probe in the UV-visible regime. Chirality is most commonly associated with a localized chiral center. However, some compounds such as helicenes (Figure 1) are chiral due to their screwlike global structure. In these highly conjugated systems, some electric and magnetic allowed transitions are distributed across the entire molecule, and OCD thus probes the global molecular chirality. Recent advances in X-ray sources, in particular the control of their polarization and spatial profiles, have enabled X-ray circular dichroism (XCD), which, in contrast to OCD, can exploit the localized and element-specific nature of X-ray electronic transitions. XCD therefore is more sensitive to local structures, and the chirality probed with it can be referred to as local. During the racemization of helicene, between opposite helical structures, the screw handedness can flip locally, making the molecule globally achiral while retaining a local handedness. Here, we use the racemization mechanism of [12]helicene as a model to demonstrate the capabilities of OCD and XCD as time-dependent probes for global and local chiralities, respectively. Our simulations demonstrate that XCD provides an excellent spectroscopic probe for the time-dependent local chirality of molecules.

Enantiopure cycloplatinated pentahelicenic N-heterocyclic carbenic complexes that display long-lived circularly polarized phosphorescence

The preparation of the first enantiopure cycloplatinated complexes bearing a bidentate, helicenic N-heterocyclic carbene and a diketonate ancillary ligand is presented, along with their structural and spectroscopic characterization based on both experimental and computational studies. The systems exhibit long-lived circularly polarized phosphorescence in solution and in doped films at room temperature, and also in a frozen glass at 77 K, with dissymmetry factor g_lum values ≥10^-3 in the former and around 10^-2 in the latter.

Polychlorinated cyclopentenes from a marine derived Periconia sp. (strain G1144)

Studies on an organic extract of a marine fungus, Periconia sp. (strain G1144), led to the isolation of three halogenated cyclopentenes along with the known and recently reported rhytidhyester D; a series of spectrometric and spectroscopic techniques were used to elucidate these structures. Interestingly, two of these compounds represent tri-halogenated cyclopentene derivatives, which have been observed only rarely from Nature. The relative and absolute configurations of the compounds were established via mass spectrometry (MS), nuclear magnetic resonance (NMR) spectroscopy, Mosher's esters method, optical rotation and GIAO NMR calculations, including correlation coefficient calculations and the use of both DP4+ and dJ DP4 analyses. Several of the isolated compounds were tested for activity in anti-parasitic, antimicrobial, quorum sensing inhibition, and cytotoxicity assays and were shown to be inactive.

Diterpenoids and p-methoxycinnamic acid diol esters from Kaempferia saraburiensis Picheans. (Zingiberaceae): Structural assignment of saraburol and their biological activities

Five undescribed compounds, including three diterpenoids namely, saraburol, saraburanes A and B, and two p-methoxycinnamic acid monoterpene diol esters, named E/Z-saraburinic esters, together with ten known oxygenated isopimarane diterpenoids, were isolated from the whole plant of Kaempferia saraburiensis Picheans. Among these compounds, saraburol possesses an unusual 6/9/6 tricyclic ring system bearing a 1,3-dioxonane-4-one scaffold, which is rarely found in natural products. The structure of isolated compounds was elucidated by spectroscopic methods, including HRESIMS, FTIR, 1D and 2D-NMR, and by comparison with published data, and their absolute configurations were determined by comparison of experimental with calculated ECD spectra and hydrolysis reaction. Using gauge-independent atomic orbital (GIAO) NMR shift calculations coupled with DP4+ probability analyses, biogenetic considerations, and optical rotation allowed for the complete characterization of saraburol. A plausible biosynthetic pathway for saraburol and saraburane A was proposed. The cytotoxicity result indicated that E-saraburinic ester exhibited the most potent activity with an IC₅₀ value of 12.0 μM against MOLT-3 cells with a selectivity index of 12.5. Saraburane B exhibited the most potent activity against Gram-positive bacteria strain Staphylococcus epidermidis with MIC (MBC) value of 25 (50) μg/mL.

Finding reliable methodology for optical rotation and correct predictions of (s)-methyloxirane and (1R,5R)-β-pinene

Theoretical calculations of optical rotation (OR), although important to predict absolute configurations (ACs) and corroborate experiments, require efficient methodology able to reproduce enantiomer specificity and real OR values. Also, troublesome molecules are recurring in the literature, such as (S)-methyloxirane and (1R,5R)-β-pinene. This study evaluates DFT functionals B3LYP, CAM-B3LYP, ωB97X-D, M06-2X, and PBE0 considering basis sets aug-cc-pVDZ, aug-cc-pVTZ, 6-311++G(2d,p), and 6-311++G(3df,2p) in OR prediction of 42 rigid organic molecules assessing cases with wrong enantiomeric determination comparing to available experimental data at wavelengths 355, 589, and 633 nm. Functionals CAM-B3LYP and ωB97X-D with aug-cc-pVTZ are indicated here to reproduce experimental values more accurately considering fewer number of wrong AC predictions, normalized RMSD values below 0.70, and a good approximation to experimental values in hierarchical cluster analysis. Methyloxirane AC was reproduced in CAM-B3LYP and PBE0, with [ α ] 355 = 6 . 94 $$ {\left[\alpha \right]}_{355}=6.94 $$ for CAM-B3LYP/aug-cc-pVTZ close to experimental value [ α ] 355 = 7 . 49 ± 0 . 30 $$ {\left[\alpha \right]}_{355}=7.49\pm 0.30 $$ . Good results were found for AC of β-pinene in M06-2X, CAM-B3LYP, and ωB97X-D while the latter in 6-311++G(3df,2p) obtained OR values of [ α ] 589 = 3 . 44 $$ {\left[\alpha \right]}_{589}=3.44 $$ and [ α ] 689 = 4 . 20 $$ {\left[\alpha \right]}_{689}=4.20 $$ close to experimental values [ α ] 589 = 2 . 8 $$ {\left[\alpha \right]}_{589}=2.8 $$ and [ α ] 689 = 4 . 66 ± 0 . 60 $$ {\left[\alpha \right]}_{689}=4.66\pm 0.60 $$ . The two molecules aforementioned are, for the first time, reported to give valid theoretical OR values in such simple methodologies. OR calculations were all performed after geometry optimization at the same level of theory, and analysis of different functional combinations for each step in β-pinene showed it can interfere with AC prediction even in rigid molecules.

Comparative Study of Vibrational Raman Optical Activity with Different Time-Dependent Density Functional Approximations: The VROA36 Database

A new database, VROA36, is introduced to investigate the performance of computational approaches for vibrational Raman optical activity (VROA) calculations. The database is composed of 36 molecules with known experimental VROA spectra. It includes 93 conformers. Normal modes calculated with B3LYP-D3(BJ)/def2-TZVP are used to compute the VROA spectra with four functionals, B3LYP-D3(BJ), ωB97X-D, M11, and optimally tuned LC-PBE, as well as several basis sets. SimROA indices and frequency scaling factors are used to compare calculated spectra with each other and with experimental data. The four functionals perform equally well independently of the basis set and usually achieve good agreement with the experimental data. For molecules in near- or at-resonance conditions, the inclusion of a complex (damped) linear response approach is important to obtain physically meaningful VROA intensities. The use of any of the tested functional approximations with the def2-SVPD Gaussian-type basis set, or a basis of similar flexibility, can be recommended for efficient and reliable theoretical VROA studies.

New Cytotoxic Secondary Metabolites from Two Deep-Sea-Derived Fungi and the Co-Culture Impact on the Secondary Metabolic Patterns

In this study, chemical profiles for two co-existing deep-sea-derived Penicillium fungal strains were thoroughly investigated. Two new compounds and 11 known compounds were identified from Penicillium sp. LXY140-R, while one new compound and 12 known compounds were isolated from Penicillium sp. LXY140-3. Their structures were elucidated by extensive 1D and 2D NMR experiments, which were supported by HR-ESI-MS data. The antiproliferative activities of all isolates against HCT-116, A549 and Bel-7402 cell lines were also evaluated. Compounds 2, 5, 6, 10 and 13 showed potent antiproliferative activity. To reveal the metabolic relationship of the two strains, we conducted co-culture experiments to discover cross-talk molecules by a device that allows only small molecule to communicate. Extensive HPLC/MS² experiments were applied to identify the disturbance of the chemical profiles within the synthetic Penicillium-Penicillium community. The fungal strain LXY140-R was found to accumulate mono or multiple-acetylation derivatives of deoxynivalenol (DON) sesquiterpenes as responsible molecules by the disturbance of the metabolites produced by the LXY140-3 strain.

Basis Set Dependence of Optical Rotation Calculations with Different Choices of Gauge

In this work, the basis set dependence of optical rotation (OR) calculations is examined for various choices of gauge/level of theory. The OR is calculated for a set of 50 molecules using B3LYP and CAM-B3LYP and 17 molecules using coupled cluster with single and double excitations (CCSD). The calculations employ the correlation-consistent basis sets, aug-cc-pVζZ with ζ = D, T, Q. An inverse-power extrapolation formula is then utilized to obtain OR values at the complete basis set (CBS) limit. We investigate the basis set convergence for these methods and three choices of gauge: length gauge (with gauge-including atomic orbitals, LG(GIAOs), for DFT), the origin-invariant length gauge [LG(OI)], and the modified velocity gauge (MVG). The results show that all methods converge smoothly to the CBS limit and that the LG(OI) approach has a slightly faster convergence rate than the other choices of gauge. While the DFT methods reach gauge invariance at the CBS limit, CCSD does not. The significant difference between the MVG and LG(OI) results at the CBS limit, 26%, indicates that CCSD is not quite at convergence in the description of electron correlation for this property. On the other hand, gauge invariance at the CBS limit for DFT does not lead to the same OR values for the two density functionals, which is also due to electron correlation incompleteness. A limited comparison to gas-phase experimental OR values for the DFT methods shows that CAM-B3LYP seems more accurate than B3LYP. Overall, this study shows that the LG(OI) approach with the aug-cc-pVTZ basis set for DFT, and with the CBS(DT) extrapolation for CCSD, provides a good cost/accuracy balance.

Automated Quantum Chemistry-Based Calculation of Optical Rotation for Large Flexible Molecules

The calculation of optical rotation (OR, [α]_D) for nonrigid molecules was limited to small systems due to the challenging problem of generating reliable conformer ensembles, calculating accurate Boltzmann populations and the extreme sensitivity of the OR to the molecules' three-dimensional structure. Herein, we describe and release the crenso workflow for the automated computation of conformer ensembles in solution and corresponding [α]_D values for flexible molecules. A comprehensive set of 28 organic drug molecules (28-144 atoms) with experimentally determined values is used in our assessment. In all cases, the correct OR sign is obtained with an overall mean relative deviation of 72% (mean absolute deviation of 82 °[dm(g/cm³)]^-1 for experimental values in the range -160 to 287 °[dm(g/cm³)]^-1). We show that routine [α]_D computations for very flexible, biologically active molecules are both feasible and reproducible in about a day of computation time on a standard workstation computer. Furthermore, we observed that the effect of energetically higher-lying structures in the ensemble on the OR is often averaged out and that in 23 out of 28 cases, the correct OR sign is obtained by just considering only the lowest free energy conformer. In four example cases, we show that the approach can also describe the OR of pairs of flexible diastereomers properly. In summary, even very sensitive, multifactorial physicochemical properties appear reliably predictable with minimal user input from efficiently automated quantum chemical methods.

Opportunities and Limitations for Assigning Relative Configurations of Antibacterial Bislactones using GIAO NMR Shift Calculations

Four new bislactones, dihydroacremonol (1), clonostachyone (2), acremodiol B (3), and acremodiol C (4), along with one known compound, hymeglusin (5), were isolated from cultures of two fungal strains (MSX59876 and MSX59260). Both strains were identified based on phylogenetic analysis of molecular data as Clonostachys spp.; yet, they biosynthesized a suite of related, but different, secondary metabolites. Given the challenges associated with elucidating the structures and configurations of bislactones, GIAO NMR calculations were tested as a complement to traditional NMR and HRESIMS experiments. Fortuitously, the enantiomer of the new natural product (4) was known as a synthetic compound, and the predicted configuration from GIAO NMR calculations (i.e., for the relative configuration) and optical rotation calculations (i.e., for the absolute configuration) matched those of the synthesis product. These results engendered confidence in using similar procedures, particularly the mixture of GIAO NMR shift calculations coupled with an orthogonal technique, to predict the configuration of 1-3; however, there were important limitations, which are discussed for each of these. The metabolites displayed antimicrobial activities, with compounds 1 and 4 being the most potent against Staphylococcus aureus with MICs of 1 and 4 μg/mL, respectively.

A Density Functional Theory Study of Optical Rotation in Some Aziridine and Oxirane Derivatives

We present time-dependent density functional theory (TDDFT) calculations of the electronic optical rotation (ORP) for seven oxirane and two aziridine derivatives in the gas phase and in solution and compare the results with the available experimental values. For seven of the studied molecules it is the first time that their optical rotation was studied theoretically and we have therefore investigated the influence of several settings in the TDDFT calculations on the results. This includes the choice of the one-electron basis set, the exchange-correlation functional or the particular polarizable continuum model (PCM). We can confirm that polarized quadruple zeta basis sets augmented with diffuse functions are necessary for converged results and find that the aug-pc-3 basis set is a viable alternative to the frequently employed aug-cc-pVQZ basis set. Based on our study, we cannot recommend the generalized gradient functional KT3 for calculations of the ORP in these compounds, whereas the hybrid functional PBE0 gives results quite similar to the long-range correct CAM-B3LYP functional. Finally, we observe large differences in the solvent effects predicted by the integral equation formalism of PCM and the SMD variant of PCM. For the majority of solute/solvent combinations in this study, we find that the SMD model in combination with the PBE0 functional and the aug-pc-3 basis set gives the best agreement with the experimental values.

Full optical rotation tensor at coupled cluster with single and double excitations level in the modified velocity gauge

This work presents the first simulations of the full optical rotation (OR) tensor at coupled cluster with single and double excitations (CCSD) level in the modified velocity gauge (MVG) formalism. The CCSD-MVG OR tensor is origin independent, and each tensor element can in principle be related directly to experimental measurements on oriented systems. We compare the CCSD results with those from two density functionals, B3LYP and CAM-B3LYP, on a test set of 22 chiral molecules. The results show that the functionals consistently overestimate the CCSD results for the individual tensor components and for the trace (which is related to the isotropic OR), by 10%-20% with CAM-B3LYP and 20%-30% with B3LYP. The data show that the contribution of the electric dipole-magnetic dipole polarizability tensor to the OR tensor is on average twice as large as that of the electric dipole-electric quadrupole polarizability tensor. The difficult case of (1S,4S)-(-)-norbornenone also reveals that the evaluation of the former polarizability tensor is more sensitive than the latter. We attribute the better agreement of CAM-B3LYP with CCSD to the ability of this functional to better reproduce electron delocalization compared with B3LYP, consistent with previous reports on isotropic OR. The CCSD-MVG approach allows the computation of reference data of the full OR tensor, which may be used to test more computationally efficient approximate methods that can be employed to study realistic models of optically active materials.

Monoterpenoid indole alkaloids from the fruits of Gelsemium elegans and their anti-inflammatory activities

Two novel monoterpenoid indole alkaloids (MIAs), gelsechizines A-B (1-2), along with four known ones (3-6) were isolated from the fruits of Gelsemium elegans. Compound 1 features a new carbon skeleton with two additional carbon atoms forming a 4-methylpyridine unit. Their structures with absolute configurations were elucidated by NMR, MS, X-ray diffraction and electronic circular dichroism (ECD) calculations. Compounds 1-3 showed significant anti-inflammatory effects in vivo and in vitro, which may be related to the inhibition of the trecruitment of neutrophils and macrophages as well as the secretion of TNF-α and IL-6. Preliminary structure-activity relationship analysis revealed that the β-N-acrylate moiety plays an important role in the anti-inflammatory effect.

Dinuclear Rhenium Complexes with a Bridging Helicene-bis-bipyridine Ligand: Synthesis, Structure, and Photophysical and Chiroptical Properties

By attaching pyridine groups to a diaza[6]helicene, a helical, bis-ditopic, bis-N N-coordinating ligand can be accessed. Dinuclear rhenium complexes featuring this bridging ligand, of the form [{Re(CO)3 Cl}2 (N N-N N)], have been prepared and resolved to give enantiopure complexes. These complexes are phosphorescent in solution at room temperature under one- and two-photon excitation. Their experimental chiroptical properties (optical rotation, electronic circular dichroism and circularly polarized emission) have been measured. They show, for instance, emission dissymmetry factors of c.a. ±3x10-3 . Quantum-chemical calculations indicate the importance of stereochemistry on the optical activity, pointing towards further design improvements in such types of complexes.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

The ONIOM/PMM Model for Effective Yet Accurate Simulation of Optical and Chiroptical Spectra in Solution: Camphorquinone in Methanol as a Case Study

This paper deals with the development and first validation of a composite approach for the simulation of chiroptical spectra in solution aimed to strongly reduce the number of full QM computations without any significant accuracy loss. The approach starts from the quantum mechanical computation of reference spectra including vibrational averaging effects and taking average solvent effects into account by means of the polarizable continuum model. Next, the snapshots of classical molecular dynamics computations are clusterized and one reference configuration from each cluster is used to compute a reference spectrum. Local fluctuation effects within each cluster are then taken into account by means of the perturbed matrix model. The performance of the proposed approach is tested on the challenging case of the optical and chiroptical spectra of camphorquinone in methanol solution. Although further validations are surely needed, the results of this first study are quite promising also taking into account that agreement with experimental data is reached by just a couple of full quantum mechanical geometry optimizations and frequency computations.

Amplification of the chiroptical response of UV-transparent amines and alcohols by N-phthalimide derivatization enabling absolute configuration determination through ECD computational analysis

The stereoselective transformation of chiral UV-transparent amines and alcohols to phthalimides has proved to be a simple and efficient method to enhance the chiroptical response of these substrates allowing their reliable absolute configuration determination by computational analysis of ECD spectra. Such a transformation also leads to a significant reduction in the molecular conformational flexibility thus simplifying the conformational analysis required by the computational treatment. The method described herein thus allows the absolute configuration assignment to these challenging substrates to be much easier and reliable.

Computational prediction of chiroptical properties in structure elucidation of natural products

Covering: up to 2019This review covers the current status of the quantum mechanical prediction of chiroptical properties, such as electronic CD and optical rotation, as needed for stereochemical assignments in new natural products. The reliability of the prediction of chiroptical properties is steadily increasing, with a parallel decrease in the required computational resources. Now, quantum mechanical calculations for a medium-sized natural product can be reliably performed by natural product chemists on a mainstream PC. This review is aimed to guide natural product chemists through the numerous steps involved in such calculations. Through a concise, but comprehensive, discussion of the current computational practice, enriched by a few illustrative examples, this review provides readers with the theoretical background and practical knowledge needed to select the most appropriate parameters for performing the calculations, to anticipate possible problems, and to critically evaluate the reliability of their computational results. Common reasons for mistakes are also discussed; in particular, the importance of the correct evaluation of conformational ensembles of flexible molecules (an aspect often overlooked in current research) is stressed.

Isolation of Peribysins O, P, and Q from Periconia macrospinosa KT3863 and Configurational Reinvestigation of Peribysin E Diacetate from Periconia byssoides OUPS-N133

Peribysins O (1), P (3), and Q (4) were isolated from Periconia macrospinosa KT3863. The relative configuration of the 6,7-epoxide of 1 was elucidated by performing quantitative NOE experiments. The structure of 2, which is a tautomer of 1 present in CDCl₃ solutions in 5% abundance, was also fully characterized by NMR analysis. Their absolute configurations were independently determined by the modified Mosher's method (for 1 and 3), the electronic circular dichroism (ECD) exciton coupling theory after conversion into dibenzoate 9 (for 3), and theoretical ECD calculations (for 1, 3, and 4). The obtained relative structures 1, 3, and 4 were verified by calculating their ¹³C chemical shifts using density functional theory (DFT). Although the established (4 S)-enantiomer for 1-4 is in accordance with that of other peribysins isolated from the related fungus Periconia byssoides OUPS-N133, Danishefsky's total synthesis of peribysin E (5) led to the subsequent revision of the (2 R,4 S,5 R,6 S,7 S,8 R,10 S)-enantiomer to the (2 S,4 R,5 S,6 R,7 R,8 S,10 R)-enantiomer. This discordance led us to reinvestigate the configuration using time-dependent DFT-based ECD spectral calculations, which supported the original (4 S)-enantiomer.

Phosphahelicenes: From Chiroptical and Photophysical Properties to OLED Applications

Herein, the experimental physicochemical and chiroptical properties of a series of phosphahelicenes are reported, focusing on their UV/Vis absorption, luminescence, electronic circular dichroism, optical rotations, and circularly polarized luminescence. Furthermore, detailed analysis of absorption and ECD spectra performed with the help of quantum-chemical calculations allowed us to highlight general features of these helicenic phosphines. Finally, due to well-suited electrochemical properties and thermal stability, the systems were successfully used as emitters in organic light-emitting diodes.

Chiroptical Protocol for the Absolute Configurational Assignment of Alkyl-Substituted Epoxides Using Bis(zinc porphyrin) as a CD-Sensitive Bidentate Host

The absolute configurations of simple alkyl-substituted chiral epoxides not bearing other ligating groups are readily determined via the exciton-coupled circular dichroism (ECCD) protocol using bidentate bis(zinc porphyrin) host system BP1 as a CD-sensitive chirality probe. In this situation, chiral epoxides can successfully be incorporated into the cleft of V-shaped host BP1 by double coordination of both oxygen lone pairs of the guest to the two central zinc ions of the host. We also propose a working model based on an MM2 optimized structure of the substrates that enables nonempirical prediction of the chirality of the bound epoxide.

Resolution-of-identity accelerated relativistic two- and four-component electron dynamics approach to chiroptical spectroscopies

We present an implementation and application of electron dynamics based on real-time time-dependent density functional theory (RT-TDDFT) and relativistic 2-component X2C and 4-component Dirac-Coulomb (4c) Hamiltonians to the calculation of electron circular dichroism and optical rotatory dispersion spectra. In addition, the resolution-of-identity approximation for the Coulomb term (RI-J) is introduced into RT-TDDFT and formulated entirely in terms of complex quaternion algebra. The proposed methodology was assessed on the dimethylchalcogenirane series, C₄H₈X (X = O, S, Se, Te, Po, Lv), and the spectra obtained by non-relativistic and relativistic methods start to disagree for Se and Te, while dramatic differences are observed for Po and Lv. The X2C approach, even in its simplest one-particle form, reproduces the reference 4c results surprisingly well across the entire series while offering an 8-fold speed-up of the simulations. An overall acceleration of RT-TDDFT by means of X2C and RI-J increases with system size and approaches a factor of almost 25 when compared to the full 4c treatment, without compromising the accuracy of the final spectra. These results suggest that one-particle X2C electron dynamics with RI-J acceleration is an attractive method for the calculation of chiroptical spectra in the valence region.

New Basis Set for the Evaluation of Specific Rotation in Flexible Biological Molecules in Solution

A detailed theoretical investigation of specific rotation is carried out in solution for nine flexible molecules of biological importance. Systematic search for the main conformers is followed by time-dependent density functional theory (TD-DFT) calculations of specific rotation employing a wide range of basis sets. Due to conformational flexibility of the compounds under study, the possibility of basis set size reduction without deterioration of the results is investigated. The increasing size (d-)aug-cc-pVXZ (X = D, T, Q) bases of Dunning et al., and the ORP basis set, recently developed to efficiently provide molecular specific rotation, are used for this purpose. The polarizable continuum model is employed at all steps of the investigation. Comparison of the present results with the available data obtained in a vacuum reveals considerable differences, the values in solution being much closer to the experimental specific rotation data available. The ORP basis set proves to be competitive with the d-aug-cc-pVDZ set of Dunning in specific rotation calculations carried out in solution. While having the same number of functions, the former yields, in general, results considerably closer to the reference triple-ζ values. We can thus recommend the ORP basis set to study the optical rotation in conformationally flexible molecules in solution.

Comparison of measured and predicted specific optical rotation in gas and solution phases: A test for the polarizable continuum model of solvation

A comparative theoretical and experimental study of dispersive optical activity is presented for a set of small, rigid organic molecules in gas and solution phases. Target species were chosen to facilitate wavelength-resolved measurements of specific rotation in rarefied vapors and in organic solvents having different polarities, while avoiding complications due to conformational flexibility. Calculations were performed with two density functionals (B3LYP and CAM-B3LYP) and with the coupled-cluster singles and doubles (CCSD) ansatz, and solvent effects were included through use of the polarizable continuum model (PCM). Across the various theoretical methods surveyed, CCSD with the modified velocity gauge provided the best overall performance for both isolated and solvated conditions. Zero-point vibrational corrections to equilibrium calculations of chiroptical response tended to improve agreement with gas-phase experiments, but the quality of performance realized for solutions varied markedly. Direct comparison of measured and predicted specific-rotation suggests that PCM, in general, is not able to reproduce attendant solvent shifts (neither between gas and solution phases nor among solvents) and fares better in estimating actual medium-dependent values of this property (although the error is rather system dependent). Thus, more elaborate solvation models seem necessary for a proper theoretical description of solvation in dispersive optical activity.

Attosecond Charge Migration with TDDFT: Accurate Dynamics from a Well-Defined Initial State

We investigate the ability of time-dependent density functional theory (TDDFT) to capture attosecond valence electron dynamics resulting from sudden X-ray ionization of a core electron. In this special case the initial state can be constructed unambiguously, allowing for a simple test of the accuracy of the dynamics. The response following nitrogen K-edge ionization in nitrosobenzene shows excellent agreement with fourth-order algebraic diagrammatic construction (ADC(4)) results, suggesting that a properly chosen initial state allows TDDFT to adequately capture attosecond charge migration. Visualizing hole motion using an electron localization picture (ELF), we provide an intuitive chemical interpretation of the charge migration as a superposition of Lewis dot resonance structures.

Absolute configuration determination through the unique intramolecular excitonic coupling in the circular dichroisms of o,p′-DDT and o,p′-DDD. A combined experimental and theoretical study

Circular dichroisms (CDs) of theo,p′-isomers of 1,1,1-trichloro- and 1,1-dichloro-2,2-bis(chlorophenyl)ethanes (DDTandDDD) were investigated experimentally and theoretically.