Software update: The ORCA program system—Version 5.0

Embedding Beyond Electrostatics: The Extended Polarizable Density Embedding Model

The polarizable density embedding (PDE) model is a focused QM/QM fragment-based embedding model designed to model solvation effects on molecular properties. We extend the PDE model to include exchange and nonadditive exchange-correlation (for DFT) in the embedding potential in addition to the existing electrostatic, polarization, and nonelectrostatic effects already present. The resulting model, termed PDE-X, yields localized electronic excitation energies that accurately capture the range dependence of the solvent interaction and gives close agreement with full quantum mechanical (QM) results, even when using minimal QM regions. We show that the PDE-X embedding description consistently improves the accuracy of excitation energies for a diverse set of organic chromophores. The improved embedding description leads to systematic solvent effects that do not average out when applying configurational sampling.

CO2 Aggregation on Monoethanolamine: Observations from Rotational Spectroscopy

The initial stages of the gas-phase nucleation between CO₂ and monoethanolamine were investigated via broadband rotational spectroscopy with the aid of extensive theoretical structure sampling. Sub-nanometer-scale aggregation patterns of monoethanolamine-(CO₂ )_n , n=1-4, were identified. An interesting competition between the monoethanolamine intramolecular hydrogen bond and the intermolecular interactions between monoethanolamine and CO₂ upon cluster growth was discovered, revealing an intriguing CO₂ binding priority to the hydroxyl group over the amine group. These findings are in sharp contrast to the general results for aqueous solutions. In the quinary complex, a cap-like CO₂ tetramer was observed cooperatively surrounding the monoethanolamine. As the cluster approaches the critical size of new particle formation, the contribution of CO₂ self-assembly to the overall stability increases.

Silmitasertib (CX-4945), a Clinically Used CK2-Kinase Inhibitor with Additional Effects on GSK3β and DYRK1A Kinases: A Structural Perspective

A clinical casein kinase 2 inhibitor, CX-4945 (silmitasertib), shows significant affinity toward the DYRK1A and GSK3β kinases, involved in down syndrome phenotypes, Alzheimer's disease, circadian clock regulation, and diabetes. This off-target activity offers an opportunity for studying the effect of the DYRK1A/GSK3β kinase system in disease biology and possible line extension. Motivated by the dual inhibition of these kinases, we solved and analyzed the crystal structures of DYRK1A and GSK3β with CX-4945. We built a quantum-chemistry-based model to rationalize the compound affinity for CK2α, DYRK1A, and GSK3β kinases. Our calculations identified a key element for CK2α's subnanomolar affinity to CX-4945. The methodology is expandable to other kinase selectivity modeling. We show that the inhibitor limits DYRK1A- and GSK3β-mediated cyclin D1 phosphorylation and reduces kinase-mediated NFAT signaling in the cell. Given the CX-4945's clinical and pharmacological profile, this inhibitory activity makes it an interesting candidate with potential for application in additional disease areas.

Stabilized Molecular Diphosphorus Pentoxide, P2O5L2 (L = N-Donor Base), in the Synthesis of Condensed Phosphate-Organic Molecule Conjugates

Commercial phosphorus pentoxide reacts with some N-donor bases to give the adducts P₂O₅L₂ and P₄O₁₀L₃ (L = DABCO, pyridine, 4-tert-butylpyridine). The DABCO adducts were structurally characterized by single-crystal X-ray diffraction. It is proposed that P₂O₅L₂ and P₄O₁₀L₃ undergo interconversion through a "phosphate-walk" mechanism, which was evaluated using DFT calculations. P₂O₅(pyridine)₂ (1) efficiently transfers monomeric diphosphorus pentoxide to phosphorus oxyanion nucleophiles, yielding substituted trimetaphosphates and cyclo-phosphonate-diphosphates (P₃O₈R)^2- (R¹ = nucleosidyl, phosphoryl, alkyl, aryl, vinyl, alkynyl, H, F). Hydrolytic ring-opening of these compounds forms linear derivatives [R¹(PO₃)₂PO₃H]^3-, and nucleophilic ring-opening gives linear disubstituted [R¹(PO₃)₂PO₂R²]^3- compounds.

A Buckycatcher in Solution-A Computational Perspective

In this work, we study the buckycatcher (C₆₀H₂₈) in solution using quantum chemical models. We investigate the conformational equilibria in several media and the effects that molecules of solvent might have in interconversion barriers between the different conformers. These are studied in a hypothetical gas phase, in the dielectric of a solvent, as well as with hybrid solvation. In the latter case, due to a disruption of π-stacking interactions, the transition states are destabilized. We also evaluate the complexation of the buckycatcher with solvent-like molecules. In most cases studied, there should be no adducts formed because the enthalpy driving force cannot overcome entropic penalties.

Molecular modification effects on the electrochromic and photochromic properties of diarylethene with intramolecular isomerization behavior

Diarylethene (DAE) is one of the most widely used functional units for electrochromic or photochromic materials. To better understand the molecular modification effects on the electrochromic and photochromic properties of DAE, two modification strategies, substitution with functional groups or heteroatoms, were investigated theoretically by density functional theory calculations. It is found that red-shifted absorption spectra caused by a decreased highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap and S₀ → S₁ transition energy during the ring-closing reaction become more significant by adding different functional substituents. In addition, for two isomers, the energy gap and S₀ → S₁ transition energy decreased by heteroatom substitution of S atoms with O or NH, while they increased by replacing two S atoms with CH₂. For intramolecular isomerization, one-electron excitation is the most effective way to trigger the closed-ring (O → C) reaction, while the open-ring (C → O) reaction occurs most readily in the presence of one-electron reduction. Moreover, it is confirmed that substitution with strong electron donating groups (-OCH₃/-NH₂) or with one O/two CH₂ heteroatoms leads to a more favorable closed-ring (O → C) reaction. Functionalized with strong electron-withdrawing groups (-NO₂ and -COOH) or one/two NH heteroatom substitutions, the open-ring (C → O) reaction is easier. Our results confirmed that the photochromic and electrochromic properties of DAE can be tuned effectively by molecular modifications, which provides theoretical guidance for the design of new DAE-based photochromic/electrochromic materials.

The Electronic Nature of Cationic Group 10 Ylidyne Complexes

We report a broad theoretical study on [(PMe3)3MER]+ complexes, with M = Ni, Pd, Pt, E = C, Si, Ge, Sn, Pb, and R = ArMes, Tbb, (ArMes = 2,6-dimesitylphenyl; Tbb = C6H2-2,6-[CH(SiMe3)2]2-4-tBu). A few years ago, our group succeeded in obtaining heavier homologues of cationic group 10 carbyne complexes via halide abstraction of the tetrylidene complexes [(PMe3)3M=E(X)R] (X = Cl, Br) using a halide scavenger. The electronic structure and the M-E bonds of the [(PMe3)3MER]+ complexes were analyzed utilizing quantum-chemical tools, such as the Pipek–Mezey orbital localization method, the energy decomposition analysis (EDA), and the extended-transition state method with natural orbitals of chemical valence (ETS-NOCV). The carbyne, silylidyne complexes, and the germylidyne complex [(PMe3)3NiGeArMes]+ are suggested to be tetrylidyne complexes featuring donor–acceptor metal tetrel triple bonds, which are composed of two strong π(M→E) and one weaker σ(E→M) interaction. In comparison, the complexes with M = Pd, Pt; E = Sn, Pb; and R = ArMes are best described as metallotetrylenes and exhibit considerable M−E−C bending, a strong σ(M→E) bond, weakened M−E π-components, and lone pair density at the tetrel atoms. Furthermore, bond cleavage energy (BCE) and bond dissociation energy (BDE) reveal preferred splitting into [M(PMe3)3]+ and [ER] fragments for most complex cations in the range of 293.3–618.3 kJ·mol−1 and 230.4–461.6 kJ·mol−1, respectively. Finally, an extensive study of the potential energy hypersurface varying the M−E−C angle indicates the presence of isomers with M−E−C bond angles of around 95°. Interestingly, these isomers are energetically favored for M = Pd, Pt; E = Sn, Pb; and R = ArMes over the less-bent structures by 13–29 kJ·mol−1.

Electron Attachment to DNA: The Protective Role of Amino Acids

We have studied the effect of amino acids on the electron attachment properties of a DNA nucleobase, with cytosine as a model system. The equation of motion coupled cluster theory with an extended basis set has been used to simulate the electron-attached state of the DNA model system. Arginine, alanine, lysine, and glycine are the four amino acids considered to investigate their role in electron attachment to a DNA nucleobase. The electron attachment to cytosine in all the four cytosine-amino acid gas-phase dimer complexes follows a doorway mechanism, where the electron gets transferred from the initial dipole-bound doorway state to the final nucleobase-bound state through the mixing of electronic and nuclear degrees of freedom. When cytosine is bulk-solvated with glycine, the glycine-bound state acts as the doorway state, where the initial electron density is localized on the bulk amino acid and away from the nucleobase, thus leading to the physical shielding of the nucleobase from the incoming electron. At the same time, the presence of amino acids can increase the stability of the nucleobase-bound anionic state, which can suppress the sugar-phosphate bond rupture caused by dissociative electron attachment to DNA.

Selectivity Rule of Cryptands for Anions: Molecular Rigidity and Bonding Site

Cryptands utilize inside CH or NH groups as hydrogen bond (H-bond) donors to capture anions such as halides. In this work, the nature and selectivity of confined hydrogen bonds inside cryptands were computationally analyzed with the energy decomposition scheme based on the block-localized wavefunction method (BLW-ED), aiming at an elucidation of governing factors in the binding between cryptands and anions. It was revealed that the intrinsic strengths of inward hydrogen bonds are dominated by the electrostatic attraction, while the anion preferences (selectivity) of inner CH and NH hydrogen bonds are governed by the Pauli exchange repulsion and electrostatic interaction, respectively. Typical conformers of cages are classified into two groups, including the C_3(h) -symmetrical conformers, in which all halide anions are located near the centroids of cages, and the "semi-open" conformers, which exhibit shifted bonding sites for different halide anions. Accordingly, the difference in governing factors of selectivity is attributed to either the rigidity of cages or the binding site of anions for these two groups. In details, the C₃ conformers of NH cryptands can be enlarged more remarkably than the C_3(h) -symmetrical conformers of CH cryptands as the size of anion (ionic radius) increases, resulting in the relaxation of the Pauli repulsion and a dramatic reduction in electrostatic attraction, which eventually rules the selectivity of NH cryptands for halide anions. By contrary, the CH cryptands are more rigid and cannot effectively reduce the Pauli repulsion, which subsequently governs the anion preference. Unlike C₃ conformers whose rigidity determines the selectivity, semi-open conformers exhibit different binding sites for different anions. From F^- to I^- , the bonding site shifts toward the outside end of the pocket inside the semi-open NH cryptand, leading to the significant reduction of the electrostatic interaction that dominates the anion preference. Differently, binding sites are much less affected by the size of anion inside the semi-open CH cryptand, in which the Pauli exchange repulsion remains the key factor for the selectivity of inner hydrogen bonds.

Influence of Rare-Earth Ion Radius on Metal-Metal Charge Transfer in Trinuclear Mixed-Valent Complexes

We report the synthesis and characterization of a highly conjugated bisferrocenyl pyrrolediimine ligand, Fc₂PyrDIH (1), and its trinuclear complexes with rare earth ions─(Fc₂PyrDI)M(N(TMS)₂)₂ (2-M, M = Sc, Y, Lu, La). Crystal structures, nuclear magnetic resonance (NMR) spectra, and ultraviolet/visible/near-infrared (UV/vis-NIR) data are presented. The latter are in good agreement with DFT calculations, illuminating the impact of the rare earth ionic radius on NIR charge transfer excitations. For [2-Sc]⁺, the charge transfer is at 11,500 cm^-1, while for [2-Y]⁺, only a d-d transition at 8000 cm^-1 is observed. Lu has an ionic radius in between Sc and Y, and the [2-Lu]⁺ complex exhibits both transitions. From time-dependent density functional theory (TDDFT) analysis, we assign the 11,500 cm^-1 transition as a mixture of metal-to-ligand charge transfer (MLCT) and metal-to-metal charge transfer (MMCT), rather than pure metal-to-metal CT because it has significant ligand character. Typically, the ferrocenes moieties have high rotational freedom in bis-ferrocenyl mixed valent complexes. However, in the present (Fc₂PyrDI)M(N(TMS)₂)₂ complexes, ligand-ligand repulsions lock the rotational freedom so that rare-earth ionic radius-dependent geometric differences increasingly influence orbital overlap as the ionic radius falls. The Marcus-Hush coupling constant H_AB trends as [2-Sc]⁺ > [2-Lu]⁺ > [2-Y]⁺.

Exploring the Accuracy Limits of PNO-Based Local Coupled-Cluster Calculations for Transition-Metal Complexes

While the domain-based local pair natural orbital coupled-cluster method with singles, doubles, and perturbative triples (DLPNO-CCSD(T)) has proven instrumental for computing energies and properties of large and complex systems accurately, calculations on first-row transition metals with a complex electronic structure remain challenging. In this work, we identify and address the two main error sources that influence the DLPNO-CCSD(T) accuracy in this context, namely, (i) correlation effects from the 3s and 3p semicore orbitals and (ii) dynamic correlation-induced orbital relaxation (DCIOR) effects that are not described by the local MP2 guess. We present a computational strategy that allows us to completely eliminate the DLPNO error associated with semicore correlation effects, while increasing, at the same time, the efficiency of the method. As regards the DCIOR effects, we introduce a diagnostic for estimating the deviation between DLPNO-CCSD(T) and canonical CCSD(T) for systems with significant orbital relaxation.

Virtual screening and structure-activity relationship study of novel BTK inhibitors in Traditional Chinese Medicine for the treatment of rheumatoid arthritis

Bruton tyrosine kinase (BTK) is a known drug target for the treatment of autoimmune diseases, including rheumatoid arthritis (RA). In this study, a series of 1-amino-1H-imidazole-5-carboxamide derivatives with good inhibitory activity against BTK were selected to explore the structure-activity relationships of these BTK inhibitors (BTKIs). Furthermore, we concentrated on 182 prescriptions of Traditional Chinese Medicine with therapeutic effects on RA. 54 herbs with a frequency of ≥10 were counted to establish a database containing 4027 ingredients for virtual screening. Five compounds with relatively higher docking scores and better absorption, distribution, metabolism, elimination and toxicity (ADMET) parameters were then selected for higher precision docking. The results demonstrated that the potentially active molecules form hydrogen bond interactions with the hinge region residues Met477, Glu475, glycine-rich P-loop residue Val416, Lys430 and DFG motif Asp539. In particular, they also interact with the key residues Thr474 and Cys481 of BTK. The molecular dynamics (MD) results demonstrated that all five compounds above could bind with BTK stably as its cognate ligand in dynamic conditions. This work identified several potential BTKIs using a computer-aided drug design approach and may provide crucial information for developing novel BTKIs.Communicated by Ramaswamy H. Sarma.

Curious Case of Cobaltocenium Carbaldehyde

Cobaltocenium carbaldehyde (formylcobaltocenium) hexafluoridophosphate, a long sought-after functionalized cobaltocenium salt, is accessible from cobaltocenium carboxylic acid by a three-step synthetic sequence involving (i) chlorination to the acid chloride, (ii) copper-borohydride reduction to the hydroxymethyl derivative, and (iii) Dess-Martin oxidation to the title compound. Due to the strongly electron-withdrawing cationic cobaltocenium moiety, no standard aldehyde reactivity is observed. Instead, nucleophilic addition followed by haloform-type cleavage prevails, thereby ruling out common useful aldehyde derivatization. One-electron reduction of cobaltocenium carbaldehyde hexafluoridophosphate affords the deep-blue, isolable cobaltocene carbaldehyde 19-valence-electron radical whose spin density is located fully at cobalt and not at the formyl carbon atom. ¹H/¹³C NMR, IR, EPR spectroscopy, high-resolution mass spectrometry, cyclic voltammetry, single crystal structure analysis (XRD), and density functional theory are applied to characterize these unusual formyl-cobaltocenium/cobaltocene compounds.

Formation of 3-Oxa- and 3-Thiacyclohexyne from Ring Expansion of Heterocyclic Alkylidene Carbenes: A Mechanistic Study

The rearrangement pathways of two alkylidene carbenes appended to an oxa or thiacyclopentane into the corresponding heterocyclohexynes were elucidated using ¹³C-labeling experiments. Both carbenes exhibited a preference for migration of the allylic carbon bound to the heteroatom. Anomeric interactions involving a heteroatom lone pair and antibonding orbital of the migrating bond and inductive destabilization of the minor migratory pathway are discussed as plausible reasons for the observed trends.

Base-free synthesis of benchtop stable Ru(III)-NHC complexes from RuCl3·3H2O and their use as precursors for Ru(II)-NHC complexes

A series of Ru(III)-NHC complexes, identified as [Ru^III(PyNHC^R)(Cl)₃(H₂O)] (1a-c), have been prepared, starting from RuCl₃·3H₂O following a base-free route. The Lewis acidic Ru(III) centre operates via a halide-assisted, electrophilic C-H activation for carbene generation. The best results were obtained with azolium salts having the I^- anion, while ligand precursors with Cl^-, BF₄^-, and PF₆^- gave no complex formation and those with Br^- gave a product with mixed halides. The structurally simple, air and moisture-stable complexes represent rare examples of paramagnetic Ru(III)-NHC complexes. Furthermore, these benchtop stable Ru(III)-NHC complexes were shown to be excellent metal precursors for the synthesis of new [Ru^II(PyNHC^R)(Cl)₂(PPh₃)₂] (2a-c) and [Ru^II(PyNHC^R)(CNC^Me)I]PF₆ (3a-c) complexes. All the complexes have been characterised using spectroscopic methods, and the structures of 1a, 1b, 2c, and 3a have been determined using the single-crystal X-ray diffraction technique. This work allows easy access to new Ru-NHC complexes for the study of new properties and novel applications.

Degradation of mevinphos and monocrotophos by OH radicals in the environment: A computational investigation on mechanism, kinetic, and ecotoxicity

Known organophosphorus pesticides are used widely in agriculture to improve the production of crops. Based on the literature, the degradation of some organophosphorus pesticides was studied theoretically. However, the mechanisms and variation of toxicity during the degradation of mevinphos and monocrotophos are still unclear in the environment, especially in wastewater. In this study, the reaction mechanisms for the degradation of the two representative organophosphorus pesticides (i.e., mevinphos and monocrotophos) in presence of OH radicals in the atmosphere and water are proposed using quantum chemical methods wB97-XD/6-311 + +G(3df,2pd)//wB97-XD/6-311 + +G(d,p). Result shows that the dominant channel is OH-addition to the C atom in CC bond with energy barriers being 15.6 and 14.7 kJ/mol, in the atmosphere and water, respectively, for mevinphos. As for monocrotophos, H-abstraction from NH group via barriers of 8.2 and 10.6 kJ/mol is more feasible in both the atmosphere and water. Moreover, the subsequent reactions of the major products in the atmosphere with NO and O₂ were also studied to evaluate the atmospheric chemistry of mevinphos and monocrotophos. Kinetically, the total rate constant is 2.68 × 10^-9 and 3.86 × 10^-8 cm³ molecule^-1·s^-1 for mevinphos and monocrotophos in the atmosphere and 4.91 × 10¹⁰ and 7.77 × 10¹¹ M^-1 s^-1 in the water at 298 K, thus the lifetime is estimated to be 36.46-364.60 s (2.53-25.31 s) in the atmosphere, and 1.41 × 10^-2 - 1.41 × 10^-1 s (8.92 ×10^-4 - 8.92 ×10^-3 s) in the advanced oxidation processes (AOPs) system. Furthermore, ecotoxic predictions for rats and three aqueous organisms imply their toxicity are reduced during degradation by using ECOSAR and T.E.S.T program based quantitative structure and activity relationship (QSAR) method.

Are Accelerated and Enhanced Wave Function Methods Accurate to Compute Static Linear and Nonlinear Optical Properties?

Key components of organic-based electro-optic devices are challenging to design or optimize because they exhibit nonlinear optical responses, which are difficult to model or rationalize. Computational chemistry furnishes the tools to investigate extensive collections of molecules in the quest for target compounds. Among the electronic structure methods that provide static nonlinear optical properties (SNLOPs), density functional approximations (DFAs) are often preferred because of their low cost/accuracy ratio. However, the accuracy of the SNLOPs critically depends on the amount of exact exchange and electron correlation included in the DFA, precluding the reliable calculation of many molecular systems. In this scenario, wave function methods such as MP2, CCSD, and CCSD(T) constitute a reliable alternative to compute SNLOPs. Unfortunately, the computational cost of these methods significantly restricts the size of molecules to study, a limitation that hampers the identification of molecules with significant nonlinear optical responses. This paper analyzes various flavors and alternatives to MP2, CCSD, and CCSD(T) methods that either drastically reduce the computational cost or improve their performance but were scarcely and unsystematically employed to compute SNLOPs. In particular, we have tested RI-MP2, RIJK-MP2, RIJCOSX-MP2 (with GridX2 and GridX4 setups), LMP2, SCS-MP2, SOS-MP2, DLPNO-MP2, LNO-CCSD, LNO-CCSD(T), DLPNO-CCSD, DLPNO-CCSD(T0), and DLPNO-CCSD(T1). Our results indicate that all these methods can be safely employed to calculate the dipole moment and the polarizability with average relative errors below 5% with respect to CCSD(T). On the other hand, the calculation of higher-order properties represents a challenge for LNO and DLPNO methods, which present severe numerical instabilities in computing the single-point field-dependent energies. RI-MP2, RIJK-MP2, or RIJCOSX-MP2 are cost-effective methods to compute first and second hyperpolarizabilities with a marginal average error with respect to canonical MP2 (up to 5% for β and up to 11% for γ). More accurate hyperpolarizabilities can be obtained with DLPNO-CCSD(T1); however, this method cannot be employed to obtain reliable second hyperpolarizabilities. These results open the way to obtain accurate nonlinear optical properties at a computational cost that can compete with current DFAs.

Structure of the Green Heme Isolated from Allylbenzene-Modified Chloroperoxidase: A Novel Heme Architecture Implicating the Mechanisms of CPO Inactivation and Epoxidation

The chemical identification of the modified heme (the green heme) during chloroperoxidase catalyzed epoxidation of allylbenzene remains unestablished due to its high instability within the protein matrix, the absence of paramagnetically shifted signals, and the difficulty in obtaining crystals of the modified enzyme. We have established the unambiguous structure of the modified prosthetic heme group, which was extracted from the protein matrix using ²D NMR spectroscopy and LC-MS spectrometry. The modified heme was isolated as a µ-oxo dimer that can be quantitatively converted to the corresponding monomer. The depolymerized green heme displayed characteristic NMR signatures of iron porphyrin complexes, but no Nuclear Overhauser Effect was observable to assist in signal assignment. An alternative strategy was applied by removing the iron center of the green heme, resulting in a stable demetallated green porphyrin species. Complete assignment of all the NMR resonances in the demetallated green heme allowed us to establish the molecular architecture of the modified species as a novel N-alkylated heme. Decisive space correlations between the propyl protons of allylbenzene and the γ meso proton coupled with clear dipolar connectivities between the propyl-2H of the substrate and the β proton in the side chain of the propionic acid at carbon-6 of the porphyrin ring, clearly indicate that allylbenzene was covalently attached to the nitrogen atom of the pyrrole ring III of the prosthetic heme. In this study, the mechanism of green CPO formation and its relation to CPO catalyzed chiral transformations are also discussed. It is concluded that the double-phenyl clamp formed by two phenylalanine residues at the distal heme pocket plays a critical role in fine-tuning substrate orientation that determines the outcome of CPO catalyzed epoxidation of substituted styrenes.

The structure of pathogenic huntingtin exon 1 defines the bases of its aggregation propensity

Huntington's disease is a neurodegenerative disorder caused by a CAG expansion in the first exon of the HTT gene, resulting in an extended polyglutamine (poly-Q) tract in huntingtin (httex1). The structural changes occurring to the poly-Q when increasing its length remain poorly understood due to its intrinsic flexibility and the strong compositional bias. The systematic application of site-specific isotopic labeling has enabled residue-specific NMR investigations of the poly-Q tract of pathogenic httex1 variants with 46 and 66 consecutive glutamines. Integrative data analysis reveals that the poly-Q tract adopts long α-helical conformations propagated and stabilized by glutamine side chain to backbone hydrogen bonds. We show that α-helical stability is a stronger signature in defining aggregation kinetics and the structure of the resulting fibrils than the number of glutamines. Our observations provide a structural perspective of the pathogenicity of expanded httex1 and pave the way to a deeper understanding of poly-Q-related diseases.

IR/UV Double Resonance Study of the 2-Phenylallyl Radical and its Pyrolysis Products

Isolated 2-phenylallyl radicals (2-PA), generated by pyrolysis from a nitrite precursor, have been investigated by IR/UV ion dip spectroscopy using free electron laser radiation. 2-PA is a resonance-stabilized radical that is considered to be involved in the formation of polycyclic aromatic hydrocarbons (PAH) in combustion, but also in interstellar space. The radical is identified based on its gas-phase IR spectrum. Furthermore, a number of bimolecular reaction products are identified, showing that the self-reaction as well as reactions with unimolecular decomposition products of 2-PA form several PAH efficiently. Possible mechanisms are discussed and the chemistry of 2-PA is compared with the one of the related 2-methylallyl and phenylpropargyl radicals.

Outstanding improvement in removing the delocalization error by global natural orbital functional

This work assesses the performance of the recently proposed global natural orbital functional (GNOF) against the charge delocalization error. GNOF provides a good balance between static and dynamic electronic correlations leading to accurate total energies while preserving spin, even for systems with a highly multi-configurational character. Several analyses were applied to the functional, namely, (i) how the charge is distributed in super-systems of two fragments, (ii) the stability of ionization potentials while increasing the system size, and (iii) potential energy curves of a neutral and charged diatomic system. GNOF was found to practically eliminate the charge delocalization error in many of the studied systems or greatly improve the results obtained previously with PNOF7.

Particle size-dependent flexibility in DUT-8(Cu) pillared layer metal-organic framework

The nature of metal in the isomorphous flexible metal-organic frameworks is often reported to influence flexibility and responsivity. A prominent example of such behaviour is the DUT-8(M) family ([M₂(2,6-ndc)₂(dabco)]_n, 2,6-ndc = 2,6-naphthalene dicarboxylate, dabco = 1,4-diazabicyclo-[2.2.2]-octane), where the isostructural compounds with Ni, Zn, Co, and Cu in the paddle wheel cluster are known. The macro-sized crystals of Ni, Co, and Zn based compounds transform to the closed pore (cp) phase under desolvation and show typical gate opening behaviour upon adsorption. The choice of metal, in this case, allows the adjustment of switching kinetics, selectivity in adsorption, and gate-opening pressures. The submicron-sized crystals of of Ni, Co, and Zn based compounds remain in the open pore (op) phase after desolvation. In this contribution, we demonstrate that the presence of Cu in the paddle wheel leads to fundamentally different flexible behaviour. The DUT-8(Cu) desolvation does not lead to the formation of the cp phase, independent of the particle size regime. However, according to in situ powder diffraction analysis, the desolvated, macro-sized crystals of DUT-8(Cu)_op show breathing upon adsorption of CO₂ at 195 K. The submicron-sized particles show rigid, nonresponsive behaviour.

A Promising 1,3,5-Triazine-Based Anion Exchanger for Perrhenate Binding: Crystal Structures of Its Chloride, Nitrate and Perrhenate Salts

The reaction of pyridine with cyanuric chloride was studied under microwave activation as well as in the presence of silver nitrate. The product of hydrolysis containing two pyridinium rings and chloride anion was isolated. The structures of these anion exchanger salts with chloride, nitrate and perrhenate anions are discussed.

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted-core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller-Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol^-1 relative to the experiment. The two-point QZ/5Z and TZ/QZ CBS 1 / l max 3 extrapolation schemes are inferior to their inverse power counterpart ( 1 / l max + 1 / 2 4 ) by 1.3 and 4.3 kJ mol^-1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol^-1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double-ζ/triple-ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities

Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.

Dressed jeff-1/2 objects in mixed-valence lacunar spinel molybdates

The lacunar-spinel chalcogenides exhibit magnetic centers in the form of transition-metal tetrahedra. On the basis of density-functional computations, the electronic ground state of an Mo₄¹³⁺ tetrahedron has been postulated as single-configuration a₁² e⁴ t₂⁵, where a₁, e, and t₂ are symmetry-adapted linear combinations of single-site Mo t_2g atomic orbitals. Here we unveil the many-body tetramer wave-function: we show that sizable correlations yield a weight of only 62% for the a₁² e⁴ t₂⁵ configuration. While spin-orbit coupling within the peculiar valence orbital manifold is still effective, the expectation value of the spin-orbit operator and the g factors deviate from figures describing nominal t⁵ j_eff = 1/2 moments. As such, our data documents the dressing of a spin-orbit j_eff = 1/2 object with intra-tetramer excitations. Our results on the internal degrees of freedom of these magnetic moments provide a solid theoretical starting point in addressing the intriguing phase transitions observed at low temperatures in these materials.

Benchmarking two-body contributions to crystal lattice energies and a range-dependent assessment of approximate methods

Using the many-body expansion to predict crystal lattice energies (CLEs), a pleasantly parallel process, allows for flexibility in the choice of theoretical methods. Benchmark-level two-body contributions to CLEs of 23 molecular crystals have been computed using interaction energies of dimers with minimum inter-monomer separations (i.e., closest contact distances) up to 30 Å. In a search for ways to reduce the computational expense of calculating accurate CLEs, we have computed these two-body contributions with 15 different quantum chemical levels of theory and compared these energies to those computed with coupled-cluster in the complete basis set (CBS) limit. Interaction energies of the more distant dimers are easier to compute accurately and several of the methods tested are suitable as replacements for coupled-cluster through perturbative triples for all but the closest dimers. For our dataset, sub-kJ mol^-1 accuracy can be obtained when calculating two-body interaction energies of dimers with separations shorter than 4 Å with coupled-cluster with single, double, and perturbative triple excitations/CBS and dimers with separations longer than 4 Å with MP2.5/aug-cc-pVDZ, among other schemes, reducing the number of dimers to be computed with coupled-cluster by as much as 98%.

Computational Study of Mechanisms and Tether Length Effects of Rh-Catalyzed [3+2] and [3+2+1] Reactions of Ene/Yne-Vinylcyclopropanes

DFT calculations have been applied to study the mechanisms of [3+2] and [3+2+1] reactions of ene/yne-vinylcyclopropanes (shorted as ene/yne-VCPs). The [3+2] reactions of ene/yne-VCPs start from C-C cleavage of cyclopropane (CP cleavage) to form six-membered rhodacycle, followed by alkene/alkyne insertion and reductive elimination. The [3+2+1] reactions have two competing pathways, one is the [3+2+1] pathway (CP cleavage, ene/yne insertion, CO insertion and reductive elimination) and the other is the [3+1+2] pathway (CP cleavage, CO insertion, ene/yne insertion and reductive elimination). The length of tether in substrates affects the ene/yne insertion steps in these cycloadditions, making some reactions fail or changing the reaction pathways. The reasons for these tether length effects are discussed.

FCclasses3: Vibrationally-resolved spectra simulated at the edge of the harmonic approximation

We introduce FCclasses3, a code to carry out vibronic simulations of electronic spectra and nonradiative rates, based on the harmonic approximation. Key new features are: implementation of the full family of vertical and adiabatic harmonic models, vibrational analysis in curvilinear coordinates, extension to several electronic spectroscopies and implementation of time-dependent approaches. The use of curvilinear valence internal coordinates allows the adoption of quadratic model potential energy surfaces (PES) of the initial and final states expanded at arbitrary configurations. Moreover, the implementation of suitable projectors provides a robust framework for defining reduced-dimensionality models by sorting flexible coordinates out of the harmonic subset, so that they can then be treated at anharmonic level, or with mixed quantum classical approaches. A set of tools to facilitate input preparation and output analysis is also provided. We show the program at work in the simulation of different spectra (one and two-photon absorption, emission and resonance Raman) and internal conversion rate of a typical rigid molecule, anthracene. Then, we focus on absorption and emission spectra of a series of flexible polyphenyl molecules, highlighting the relevance of some of the newly implemented features. The code is freely available at http://www.iccom.cnr.it/en/fcclasses/.

The Structural Basis of African Swine Fever Virus pS273R Protease Binding to E64 through Molecular Dynamics Simulations

Identification of novel drugs for anti-African swine fever (ASF) applications is of utmost urgency, as it negatively affects pig farming and no effective vaccine or treatment is currently available. African swine fever virus (ASFV) encoded pS273R is a cysteine protease that plays an important role in virus replication. E64, acting as an inhibitor of cysteine protease, has been established as exerting an inhibitory effect on pS273R. In order to obtain a better understanding of the interaction between E64 and pS273R, common docking, restriction docking, and covalent docking were employed to analyze the optimal bonding position between pS273R-E64 and its bonding strength. Additionally, three sets of 100 ns molecular dynamics simulations were conducted to examine the conformational dynamics of pS273R and the dynamic interaction of pS273R-E64, based on a variety of analytical methods including root mean square deviation (RMSD), root mean square fluctuation (RMSF), free energy of ligand (FEL), principal component analysis (PCA), and molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) analysis. The results show that E64 and pS273R exhibited close binding degrees at the activity center of ASFV pS273R protease. The data of these simulations indicate that binding of E64 to pS273R results in a reduction in flexibility, particularly in the ARM region, and a change in the conformational space of pS273R. Additionally, the ability of E64 to interact with polar amino acids such as ASN158, SER192, and GLN229, as well as charged amino acids such as LYS167 and HIS168, seems to be an important factor in its inhibitory effect. Finally, Octet biostratigraphy confirmed the binding of E64 and pS273R with a KD value of 903 uM. Overall, these findings could potentially be utilized in the development of novel inhibitors of pS273R to address the challenges posed by ASFV.

The hydration of Li+ and Mg2+ in subnano carbon nanotubes using a multiscale theoretical approach

The separation of brines with high Mg/Li mass ratios is a huge challenge. To provide a theoretical basis for the design of separation materials, the hydration of Li+ and Mg2+ in confinement using carbon nanotubes (CNTs) as the 1-D nanopore model was investigated using a multiscale theoretical approach. According to the analysis of the first coordination layer of cations, we determined that the coordination shells of two cations exist inside CNTs, while the second coordination shells of the cations are unstable. Moreover, the results of the structure analysis indicate that the hydration layer of Li+ is not complete in CNTs with diameters of 0.73, 0.87, and 1.00 nm. However, this does not occur in the 0.60 nm CNT, which is explained by the formation of contact ion pairs (CIP) between Li+ and Cl- that go through a unstable solvent-shared ion pair [Li(H2O)4]+, and this research was further extended by 400 ns in the 0.60 nm CNT to address the aforementioned results. However, the hydration layer of Mg2+ is complete and not sensitive to the diameter of CNTs using molecular dynamics simulation and an ab initio molecular dynamics (AIMD) method. Furthermore, the results of the orientation distribution of Li+ and Mg2+ indicate that the water molecules around Mg2+ are more ordered than water molecules around Li+ in the CNTs and are more analogous to the bulk solution. We conclude that it is energetically unfavorable to confine Li+ inside the 0.60-nm diameter CNT, while it is favorable for confining Li+ inside the other four CNTs and Mg2+ in all CNTs, which is driven by the strong electrostatic interaction between cations and Cl-. In addition, the interaction between cations and water molecules in the five CNTs was also analyzed from the non-covalent interaction (NCI) perspective by AIMD.

Synthesis, enzyme inhibition, and molecular docking studies of a novel chalcone series bearing benzothiazole scaffold

This study reports the facile synthesis of a novel series of benzothiazole-chalcones, in addition to their inhibitory profile on important metabolic enzymes including human carbonic anhydrases (hCA-I, hCA-II) and paraoxonase (PON-1). The inhibition parameters, IC₅₀ (concentration for 50% inhibition) and Ki (dissociation constant) values, toward the title enzymes were determined for the studied compounds. As a result, IC₅₀ values of hydratase activity were in the range 4.15-5.47 and 2.56-4.58 μM for hCA-I and hCA-II, respectively. At the same time, IC₅₀ values of esterase activity were in the range 24.91-104.00 and 35.25-97.00 μM, while Ki values were in the range 14.43-59.66 and 26.65-73.34 μM for hCA-I and hCA-II, respectively. In addition, PON-1 enzyme inhibition results showed interesting inhibitory effects, with IC₅₀ values between 13.28 and 16.68 μM. Finally, a comprehensive approach was established for the synthesized compounds based on theoretical calculations, which have been done using B3LYP, PBE0 theories and SVP, TVZP, TVZPP basis sets, followed by docking studies by which the outputs proved the harmonically flows with the experimental results.

Exploring the Limits of Intramolecular London Dispersion Stabilization with Bulky Dispersion Energy Donors in Alkane Solution

We present an experimental study of a cyclooctatetraene-based molecular balance disubstituted with increasingly bulky tert-butyl (tBu), adamantyl (Ad), and diamantyl (Dia) substituents in the 1,4-/1,6-positions for which we determined the valence-bond shift equilibrium in n-hexane (hex), n-octane (oct), and n-dodecane (dod). Computations including implicit and explicit solvation support our temperature-dependent NMR equilibrium measurements indicating that the more sterically crowded 1,6-isomer is always favored, irrespective of solvent, and that the free energy is quite insensitive to substituent size.

Theoretical exploration of noncovalent interactions in Sc2C2@C2n (n = 40, 41, and 42)⊂[12]CPP, PF[12]CPP

The encapsulation of fullerenes by carbon nanorings has gained increasing attention because of the unique molecular structure and special properties of the formed complexes. The host-guest interactions between the fullerenes and the carbon nanorings can influence the metal ion orientation and the molecular electronic structure. In this study, we hooped a series of carbide cluster metallofullerenes, namely Sc₂C₂@C_2v(5)-C₈₀, Sc₂C₂@C_3v(8)-C₈₂, and Sc₂C₂@D_2d(23)-C₈₄, with molecular carbon nanorings of [12]cycloparaphenylene ([12]CPP) and perfluoro[12]cycloparaphenylene (PF[12]CPP). The formed complexes were computationally studied via dispersion-corrected density functional theory calculations. The results showed that the deformation rate of PF[12]CPP after the formation of the fullerene-containing complexes was significantly smaller than that of [12]CPP. The binding energy and thermodynamic information showed that PF[12]CPP was more suitable for fullerene encapsulation. Moreover, charge population analysis showed that PF[12]CPP transferred more electrons to Sc₂C₂@C_2n (n = 40, 41, and 42) compared with [12]CPP. Energy decomposition and real-space function analyses of host-guest interactions revealed the characteristics and nature of the noncovalent interactions in the supramolecules. These results provide theoretical support for the study of host-guest systems based on metallofullerenes.

Novel Design of Co-Poly(Hydrazide Imide) and Its Complex with Cu(I) for Membrane Separation of Methanol/Dimethyl Carbonate Mixture

Poly(2,2'-biquinoline-6,6'-dicarbohydrazide)-co-(bistrimelliteimide)methylene-bisanthranylide (PHI) and its metal-polymer complex PHI-Cu(I) containing several types of functional groups (hydrazide, carboxyl, amide, and imide fragments) were synthesized to prepare two types of dense nonporous membranes. The study on morphology using scanning electron microscopy (SEM), measurements of mechanical, thermal, and transport properties of the membrane samples was carried out. The main mechanical properties of both membranes do not differ significantly, but the values of ultimate deformation differ palpably as a result of a non-uniform character of the deformation process for the PHI membrane. The thermal analysis based on the curves of thermogravimetric (TGA) and differential thermal (DTA) analyses of the PHI and PHI-Cu(I) membranes revealed peculiarities of the membrane structure. Transport properties were studied in pervaporation (PV) of methanol (MeOH) and dimethyl carbonate (DMC) mixtures including an azeotropic point. Intrinsic properties of the penetrant-membrane system were also determined. It was found that the total flux is higher through the PHI membrane, but the PHI-Cu(I) membrane exhibits a higher separation factor. Calculation of the pervaporation separation index (PSI) allowed to conclude that the PHI-Cu(I) membrane exhibits better transport properties as compared with the PHI membrane.

Effect of void-carbon on blue-shifted luminescence in TADF molecules by theoretical simulations

Thermally activated delayed fluorescence (TADF) molecules have a theoretical 100% photoluminescence quantum yield in comparison with traditional fluorescent materials, leading to broad application in organic light-emitting diode (OLED). However, the application of TADF molecules with conjugated donor-acceptor structures in blue OLED remains a challenge due to their generally narrow energy gap between frontier molecular orbitals. Recently, a strategy has been approved in the improvement of the performance in TADF, in which void-carbon atoms between donor and acceptor fragments (donor-void-acceptor (D-v-A)) could regulate blue light emission. In this study, we first select three reported isomers followed by two proposed D-v-A TADF isomers to verify the feasibility of the void-carbon strategy through evaluation of the electronic structures in the excited state and photophysical properties. We further proposed a series of TADF molecules by replacing different donor and acceptor fragments to assess the applicability of the void-carbon strategy from the aspect of simulations in electronic structures, different properties of donor and acceptor fragments, photophysical properties, and analysis in the molecular conjugation. The results indicate that void-carbon strategy has conditional feasibility and applicability. Donor-acceptor molecular properties could be tuned through void-carbon strategy on aromatic acceptor fragments during the selection of promising candidates of TADF molecules. However, the void-carbon strategy does not work for the molecules with antiaromatic acceptor fragments, where the steric hindrance of the molecules plays a dominant role. Our work provides insightful guidance for the design of the blue-emission TADF molecules.

Assisted Tandem Pd Catalysis Enables Regiodivergent Heck Arylation of Transiently Generated Substituted Enol Ethers

Two complementary regiodivergent Pd-catalyzed assisted tandem [isomerization/Heck arylation] reactions are reported. They provide access to a broad array of acyclic trisubstituted vinyl ethers starting from readily available alkenyl ethers. In both cases, the isomerization is conducted with a [Pd-H] precatalyst supported by tris-tert-butyl phosphine ligands. When the catalyst is modified by the addition of a chelating bisphosphine ligand (dppp), an organic base (Cy2NMe), sodium acetate, and aryl triflates are used as electrophiles, the α-arylation pathway is promoted preferentially. The β-arylation pathway is favored for electron-deficient and electron-neutral aryl halides when the catalyst is simply modified by the addition of an excess of an organic base (Et3N) after completion of the isomerization reaction. Electron-rich aryl halides lead to reduced levels of regiocontrol. The moderate stereoselectivity obtained are proposed to reflect the absence of stereocontrol in the isomerization step. Computational analyses suggest that migratory insertion is selectivity-determining for both the arylations. For the β-selective arylation, an energy decomposition analysis underscored that electronic factors favor α-regioselectivity and steric effects favor β-regioselectivity. Preliminary investigations show that high levels of stereoselectivity can be achieved for the α-selective arylation by ligand control. Complementarily, reaction conditions for postcatalytic stereo-correction have also been identified for each catalytic system.

Diboraanthracene-Doped Polymer Systems for Colour-Tuneable Room-Temperature Organic Afterglow

Organic ultralong room temperature phosphorescence (RTP), or organic afterglow, is a unique phenomenon, gaining widespread attention due to its far-reaching application potential and fundamental interest. Here, two laterally expanded 9,10-dimesityl-dihydro-9,10-diboraanthracene (DBA) derivatives are demonstrated as excellent afterglow materials for red and blue-green light emission, which is traced back to persistent thermally activated delayed fluorescence and RTP. The lateral substitution of polycyclic DBA scaffold, together with weak transversal electron-donating mesityl groups, ensures the optimal molecular properties for (reverse) intersystem crossing and long-lived triplet states in a rigid poly(methyl methacrylate) matrix. The achieved afterglow emission quantum yields of up to 3 % and 15 %, afterglow lifetimes up to 0.8 s and 3.2 s and afterglow durations up to 5 s and 25 s (for red and blue-green emitters, respectively) are attributed to the properties of single molecules.

Probing the Size Limit of Dispersion Energy Donors with a Bifluorenylidene Balance: Magic Cyclohexyl

We report the synthesis of 14 2,2'-disubstituted 9,9'-bifluorenylidenes as molecular balances for the quantification of London dispersion interactions between various dispersion energy donors. For all balances, we measured ΔG_Z/E at 333 K using ¹H NMR in seven organic solvents. For various alkyl and aryl substituents, we generally observe a preference for the "folded" Z-isomer due to attractive London dispersion interactions. The cyclohexyl-substituted system shows the largest Z-preference in this study with ΔG_Z/E = -0.6 ± 0.05 kcal mol^-1 in all solvents, owing to the rotational freedom of cyclohexyl groups paired with their large polarizability that maximizes London dispersion interactions. On the other hand, rigid and sterically more demanding substituents like tert-butyl unexpectedly favor the unfolded E-isomer. This is a result of the close relative position in which the functional groups are positioned in this molecular balance. This close proximity is the reason for the increase of Pauli repulsion in the Z-isomers with large rigid substituents (tert-butyl, adamantyl, and diamantyl) which leads to an equilibrium shift toward the unfolded E-form. While we were able to reproduce most of our experimental trends qualitatively using contemporary computational chemistry methods, quantitative accuracy of the employed methods still needs further improvement.

Uncovering the Mechanism of Thermally Activated Delayed Fluorescence in Coplanar Emitters Using Potential Energy Surface Analysis

Planarized emitters exhibiting thermally activated delayed fluorescence (TADF) have attracted attention due to their narrow emission spectra, improved photostability, and high quantum yields, but with large singlet-triplet energy gaps (ΔE_ST) and no heavy atoms, the origin of their TADF remains a subject of debate. Here we prepare two isomeric, coplanar donor-acceptor compounds, with HMAT-2PYM performing dual TADF and room-temperature phosphorescence but with HMAT-4PYM exhibiting only prompt fluorescence. Although conventional TADF design principles suggest that neither isomer should exhibit TADF, we reveal differences in the excited state potential energy surfaces that enable spin-flip processes in only one isomer. We also find that hydrogen bonding is absent between the planar units of these emitters, despite earlier claims of intramolecular hydrogen bonding in similar compounds. Overall, this work demonstrates that potential energy surface analysis is a practical strategy for designing coplanar TADF materials that might otherwise be overlooked by conventional TADF design metrics.

Experimental-theoretical study of laccase as a detoxifier of aflatoxins

We investigate laccase-mediated detoxification of aflatoxins, fungal carcinogenic food contaminants. Our experimental comparison between two aflatoxins with similar structures (AFB₁ and AFG₂) shows significant differences in laccase-mediated detoxification. A multi-scale modeling approach (Docking, Molecular Dynamics, and Density Functional Theory) identifies the highly substrate-specific changes required to improve laccase detoxifying performance. We employ a large-scale density functional theory-based approach, involving more than 7000 atoms, to identify the amino acid residues that determine the affinity of laccase for aflatoxins. From this study we conclude: (1) AFB₁ is more challenging to degrade, to the point of complete degradation stalling; (2) AFG₂ is easier to degrade by laccase due to its lack of side products and favorable binding dynamics; and (3) ample opportunities to optimize laccase for aflatoxin degradation exist, especially via mutations leading to π-π stacking. This study identifies a way to optimize laccase for aflatoxin bioremediation and, more generally, contributes to the research efforts aimed at rational enzyme optimization.

Building on the strengths of a double-hybrid density functional for excitation energies and inverted singlet-triplet energy gaps

It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet-triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over "safe" excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet-triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature-the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S₁ and T₁ excited states.

The Halogen Bond in Weakly Bonded Complexes and the Consequences for Aromaticity and Spin-Orbit Coupling

The halogen bond complexes CF3X⋯Y and C2F3X⋯Y, with Y = furan, thiophene, selenophene and X = Cl, Br, I, have been studied by using DFT and CCSD(T) in order to understand which factors govern the interaction between the halogen atom X and the aromatic ring. We found that PBE0-dDsC/QZ4P gives an adequate description of the interaction energies in these complexes, compared to CCSD(T) and experimental results. The interaction between the halogen atom X and the π-bonds in perpendicular orientation is stronger than the interaction with the in-plane lone pairs of the heteroatom of the aromatic cycle. The strength of the interaction follows the trend Cl < Br < I; the chalcogenide in the aromatic ring nor the hybridization of the C−X bond play a decisive role. The energy decomposition analysis shows that the interaction energy is dominated by all three contributions, viz., the electrostatic, orbital, and dispersion interactions: not one factor dominates the interaction energy. The aromaticity of the ring is undisturbed upon halogen bond formation: the π-ring current remains equally strong and diatropic in the complex as it is for the free aromatic ring. However, the spin-orbit coupling between the singlet and triplet π→π* states is increased upon halogen bond formation and a faster intersystem crossing between these states is therefore expected.

Triplets with a Twist: Ultrafast Intersystem Crossing in a Series of Electron Acceptor Materials Driven by Conformational Disorder

Control over the populations of singlet and triplet excitons is key to organic semiconductor technologies. In different contexts, triplets can represent an energy loss pathway that must be managed (i.e., solar cells, light-emitting diodes, and lasers) or provide avenues to improve energy conversion (i.e., photon upconversion and multiplication systems). A key consideration in the interplay of singlet and triplet exciton populations in these systems is the rate of intersystem crossing (ISC). In this work, we design, measure, and model a series of new electron acceptor molecules and analyze them using a combination of ultrafast transient absorption and ultrafast broadband photoluminescence spectroscopies. We demonstrate that intramolecular triplet formation occurs within several hundred picoseconds in solution and is accelerated considerably in the solid state. Importantly, ISC occurs with sufficient rapidity to compete with charge formation in modern organic solar cells, implicating triplets in intrinsic exciton loss channels in addition to charge recombination. Density functional theory calculations reveal that ISC occurs in triplet excited states characterized by local deviations from orbital π-symmetry associated with rotationally flexible thiophene rings. In disordered films, structural distortions, therefore, result in significant increases in spin-orbit coupling, enabling rapid ISC. We demonstrate the generality of this proposal in an oligothiophene model system where ISC is symmetry-forbidden and show that conformational disorder introduced by the formation of a solvent glass accelerates ISC, outweighing the lower temperature and increased viscosity. This proposal sheds light on the factors responsible for facile ISC and provides a simple framework for molecular control over spin states.