r2SCAN-3c: A "Swiss army knife" composite electronic-structure method

Toward Reliable Conformational Energies of Amino Acids and Dipeptides─The DipCONFS Benchmark and DipCONL Datasets

Simulating peptides and proteins is becoming increasingly important, leading to a growing need for efficient computational methods. These are typically semiempirical quantum mechanical (SQM) methods, force fields (FFs), or machine-learned interatomic potentials (MLIPs), all of which require a large amount of accurate data for robust training and evaluation. To assess potential reference methods and complement the available data, we introduce two sets, DipCONFL and DipCONFS, which cover large parts of the conformational space of 17 amino acids and their 289 possible dipeptides in aqueous solution. The conformers were selected from the exhaustive PeptideCS dataset by Andris et al. [ J. Phys. Chem. B 2022, 126, 5949-5958]. The structures, originally generated with GFN2-xTB, were reoptimized using the accurate r2SCAN-3c density functional theory (DFT) composite method including the implicit CPCM water solvation model. The DipCONFS benchmark set contains 918 conformers and is one of the largest sets with highly accurate coupled cluster conformational energies so far. It is employed to evaluate various DFT and wave function theory (WFT) methods, especially regarding whether they are accurate enough to be used as reliable reference methods for larger datasets intended for training and testing more approximated SQM, FF, and MLIP methods. The results reveal that the originally provided BP86-D3(BJ)/DGauss-DZVP conformational energies are not sufficiently accurate. Among the DFT methods tested as an alternative reference level, the revDSD-PBEP86-D4 double hybrid performs best with a mean absolute error (MAD) of 0.2 kcal mol-1 compared with the PNO-LCCSD(T)-F12b reference. The very efficient r2SCAN-3c composite method also shows excellent results, with an MAD of 0.3 kcal mol-1, similar to the best-tested hybrid ωB97M-D4. With these findings, we compiled the large DipCONFL set, which includes over 29,000 realistic conformers in solution with reasonably accurate r2SCAN-3c reference conformational energies, gradients, and further properties potentially relevant for training MLIP methods. This set, also in comparison to DipCONFS, is used to assess the performance of various SQM, FF, and MLIP methods robustly and can complement training sets for those.

Molecular interactions of photosystem I and ZIF-8 in bio-nanohybrid materials

Bio-nanohybrid devices featuring natural photocatalysts bound to a nanostructure hold great promise in the search for sustainable energy conversion. One of the major challenges of integrating biological systems is protecting them against harsh environmental conditions while retaining, or ideally enhancing their photophysical properties. In this mainly computational work we investigate an assembly of cyanobacterial photosystem I (PS I) embedded in a metal-organic framework (MOF), namely the zeolitic imidazolate framework ZIF-8. This complex has been reported experimentally [Bennett et al., Nanoscale Adv., 2019, 1, 94] but so far the molecular interactions between PS I and the MOF remained elusive. We show via absorption spectroscopy that PS I remains intact throughout the encapsulation-release cycle. Molecular dynamics (MD) simulations further confirm that the encapsulation has no noticeable structural impact on the photosystem. However, the MOF building blocks frequently coordinate to the Mg2+ ions of chlorophylls in the periphery of the antenna complex. High-level quantum mechanical calculations reveal charge-transfer interactions, which affect the excitonic network and thereby may reversibly change the fluorescence properties of PS I. Nevertheless, our results highlight the stability of PS I in the MOF, as the reaction center remains unimpeded by the heterogeneous environment, paving the way for applications in the foreseeable future.

Emerging Conformational-Analysis Protocols from the RTCONF55-16K Reaction Thermochemistry Conformational Benchmark Set

RTCONF55-16K is a new, reactive conformational data set based on cost-efficient methods to assess different conformational analysis protocols. Our reference calculations underpinned the accuracy of the CENSO (Grimme et al. J. Phys. Chem. A, 2021, 125, 4039) procedure and resulted in alternative recipes with different cost-accuracy compromises. Our general-purpose and economical protocols (CENSO-light and zero, respectively) were found to be 10-30 times faster than the original algorithm, adding only 0.4-0.7 kcal/mol absolute error to the relative free energy estimates.

Isolation of a NHC-stabilized heavier nitrile and its conversion into an isonitrile analogue

Nitriles (R-C≡N) have been investigated since the late eighteenth century and are ubiquitous encounters in organic and inorganic syntheses. In contrast, heavier nitriles, which contain the heavier analogues of carbon and nitrogen, are sparsely investigated species. Here we report the synthesis and isolation of a phosphino-silylene featuring an N-heterocyclic carbene-phosphinidene and a highly sterically demanding silyl group as substituents. Due to its unique structural motif, it can be regarded as a Lewis base-stabilized heavier nitrile. The Si-P bond displays multiple bond character and a bent R-Si-P geometry, the latter indicating fundamental differences between heavier and classical nitriles. In solution, a quantitative unusual rearrangement to a phosphasilenylidene occurs. This rearrangement is consistent with theoretical predictions of rearrangements from heavier nitriles to heavier isonitriles. Our preliminary reactivity studies revealed that both isomers exhibit highly nucleophilic silicon centres capable of oxidative addition and coordination to iron tetracarbonyl.

Flow-Chemistry Based Green Synthesis of Graphene Oxide at Minutes Timescale

Graphene oxide (GO) is broadly investigated in the electrochemical field. However, for industrial applications, it still suffer from high pollution, low efficiency, poor production quality, and safety concerns associated with traditional synthesis methods. Herein, guided by theoretical analyses, a new oxygen-atom-transfer (OAT) mechanism for periodate oxidizing graphite is revealed, exhibiting controllable reaction activity, strong orbital interaction, and abundant electron transfer. Moreover, a flow chemistry strategy with high mass/heat transfer rates is designed to enhance interlayer diffusion and reaction dynamics between oxidants and graphite, ensuring the efficient synthesis of GO within several minutes. As a result, both low oxygen-content GO with large size, and high oxygen-content GO with adequate active sites can be precisely and safely synthesized. Benefitting from the controllability of oxygen content and lateral size, the as-prepared GO sheets can be facilely assembled into fiber/film electrodes that present high mechanical flexibility, large electrical conductivity, and outstanding electrochemical performance.

Exploring Novel GSK-3β Inhibitors for Anti-Neuroinflammatory and Neuroprotective Effects: Synthesis, Crystallography, Computational Analysis, and Biological Evaluation

In the pathogenesis of Alzheimer's disease, the overexpression of glycogen synthase kinase-3β (GSK-3β) stands out due to its multifaced nature, as it contributes to the promotion of amyloid β and tau protein accumulation, as well as neuroinflammatory processes. Therefore, in the present study, we have designed, synthesized, and evaluated a new series of GSK-3β inhibitors based on the N-(pyridin-2-yl)cyclopropanecarboxamide scaffold. We identified compound 36, demonstrating an IC50 of 70 nM against GSK-3β. Subsequently, through crystallography studies and quantum mechanical analysis, we elucidated its binding mode and identified the structural features crucial for interactions with the active site of GSK-3β, thereby understanding its inhibitory potency. Compound 36 was effective in the cellular model of hyperphosphorylated tau-induced neurodegeneration, where it restored cell viability after okadaic acid treatment and showed anti-inflammatory activity in the LPS model, significantly reducing NO, IL-6, and TNF-α release. In ADME-tox in vitro studies, we confirmed the beneficial profile of 36, including high permeability in PAMPA (Pe equals 9.4) and high metabolic stability in HLMs as well as lack of significant interactions with isoforms of the CYP enzymes and lack of considerable cytotoxicity on selected cell lines (IC50 > 100 μM on HT-22 cells and 89.3 μM on BV-2 cells). Based on promising pharmacological activities and favorable ADME-tox properties, compound 36 may be considered a promising candidate for in vivo research as well as constitute a reliable starting point for further studies.

State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms

In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.

Photoelectron Spectroscopy and Computational Study on Microsolvated [B10H10]2- Clusters and Comparisons to Their [B12H12]2- Analogues

Microhydrated closo-boranes have attracted great interest due to their superchaotropic activity related to the well-known Hofmeister effect and important applications in biomedical and battery fields. In this work, we report a combined negative ion photoelectron spectroscopy and quantum chemical investigation on hydrated closo-decaborate clusters [B10H10]2-·nH2O (n = 1-7) with a direct comparison to their analogues [B12H12]2-·nH2O and free water clusters. A single H2O molecule is found to be sufficient to stabilize the intrinsically unstable [B10H10]2- dianion. The first two water molecules strongly interact with the solute forming B-H···H-O dihydrogen bonds while additional water molecules show substantially reduced binding energies. Unlike [B12H12]2-·nH2O possessing a highly structured water network with the attached H2O molecules arranged in a unified pattern by maximizing B-H···H-O dihydrogen bonding, distinct structural arrangements of the water clusters within [B10H10]2-·nH2O are achieved with the water cluster networks from trimer to heptamer resembling free water clusters. Such a distinct difference arises from the variations in size, symmetry, and charge distributions between these two dianions. The present finding again confirms the structural diversity of hydrogen-bonding networks in microhydrated closo-boranes and enriches our understanding of aqueous borate chemistry.

Machine learning demonstrates the impact of proton transfer and solvent dynamics on CO2 capture in liquid ammonia

Direct air capture of CO2 using supported amines provides a promising means to achieve the net-zero greenhouse gas emissions goal; however, many mechanistic details regarding the CO2 adsorption process in condensed phase amines remain poorly understood. This work combines machine learning potentials, enhanced sampling and grand-canonical Monte Carlo simulations to directly compute experimentally relevant quantities to elucidate the mechanism of CO2 chemisorption in liquid ammonia as a model system. Our simulations suggest that CO2 capture in the liquid occurs in a sequential fashion, with the formation of a metastable zwitterion intermediate. Furthermore, we identified the importance of solvent-mediated proton transfer and solvent dynamics, not only in the reaction pathway but also in the efficiency of CO2 chemisorption. Beyond liquid ammonia, the methodology presented here can be readily extended to simulate amines with more complex chemical structures under experimental conditions, paving the way to elucidate the structure-performance of amines for CO2 capture.

Small basis set density functional theory method for cost-efficient, large-scale condensed matter simulations

We present an efficient first-principles based method geared toward reliably predicting the structures of solid materials across the Periodic Table. To this end, we use a density functional theory baseline with a compact, near-minimal min+s basis set, yielding low computational costs and memory demands. Since the use of such a small basis set leads to systematic errors in chemical bond lengths, we develop a linear pairwise correction, available for elements Z = 1-86 (excluding the lanthanide series), parameterized for use with the Perdew-Burke-Ernzerhof exchange-correlation functional. We demonstrate the reliability of this corrected approach for equilibrium volumes across the Periodic Table and the transferability to differently coordinated environments and multi-elemental crystals. We examine relative energies, forces, and stresses in geometry optimizations and molecular dynamics simulations.

Exploring ortho-dianthrylbenzenes for molecular solar thermal energy storage

Molecular solar thermal systems, which absorb light, store it, and release it as heat, have been extensively researched, yet many potential candidates remain unexplored. To expand this range, five specifically designed ortho-dianthrylbenzenes were investigated. Anthracene dimers have been underexplored due to issues like photooxidation and varying photodimerization efficiency. The presented systems address these challenges by aryl-linking two anthracene moieties, achieving photodimerization quantum yields ranging from 11.5% to 16% in mesitylene. The impact of donor or acceptor groups on energy storage time (9-37 years), energy storage density (0.14-0.2 MJ kg-1), and solar energy storage efficiency (0.38-0.66%) was evaluated. The experimental results, supported by density functional theory-based modeling, highlight the potential of anthracene-based photoswitches for molecular solar thermal applications and encourage further exploration of similar systems.

dxtb-An efficient and fully differentiable framework for extended tight-binding

Automatic differentiation (AD) emerged as an integral part of machine learning, accelerating model development by enabling gradient-based optimization without explicit analytical derivatives. Recently, the benefits of AD and computing arbitrary-order derivatives with respect to any variable were also recognized in the field of quantum chemistry. In this work, we present dxtb-an open-source, fully differentiable framework for semiempirical extended tight-binding (xTB) methods. Developed entirely in Python and leveraging PyTorch for array operations, dxtb facilitates extensibility and rapid prototyping while maintaining computational efficiency. Through comprehensive code vectorization and optimization, we essentially reach the speed of compiled xTB programs for high-throughput calculations of small molecules. The excellent performance also scales to large systems, and batch operability yields additional benefits for execution on parallel hardware. In particular, energy evaluations are on par with existing programs, whereas the speed of automatically differentiated nuclear derivatives is only 2 to 5 times slower compared to their analytical counterparts. We showcase the utility of AD in dxtb by calculating various molecular and spectroscopic properties, highlighting its capacity to enhance and simplify such evaluations. Furthermore, the framework streamlines optimization tasks and offers seamless integration of semiempirical quantum chemistry in machine learning, paving the way for physics-inspired end-to-end differentiable models. Ultimately, dxtb aims to further advance the capabilities of semiempirical methods, providing an extensible foundation for future developments and hybrid machine learning applications. The framework is accessible at https://github.com/grimme-lab/dxtb.

Extension of the D3 and D4 London dispersion corrections to the full actinides series

Efficient dispersion corrections are an indispensable component of modern density functional theory, semi-empirical quantum mechanical, and even force field methods. In this work, we extend the well established D3 and D4 London dispersion corrections to the full actinides series, francium, and radium. To keep consistency with the existing versions, the original parameterization strategy of the D4 model was only slightly modified. This includes improved reference Hirshfeld atomic partial charges at the ωB97M-V/ma-def-TZVP level to fit the required electronegativity equilibration charge (EEQ) model. In this context, we developed a new actinide data set called AcQM, which covers the most common molecular actinide compound space. Furthermore, the efficient calculation of dynamic polarizabilities that are needed to construct CAB6 dispersion coefficients was implemented into the ORCA program package. The extended models are assessed for the computation of dissociation curves of actinide atoms and ions, geometry optimizations of crystal structure cutouts, gas-phase structures of small uranium compounds, and an example extracted from a small actinide complex protein assembly. We found that the novel parameterizations perform on par with the computationally more demanding density-dependent VV10 dispersion correction. With the presented extension, the excellent cost-accuracy ratio of the D3 and D4 models can now be utilized in various fields of computational actinide chemistry and, e.g., in efficient composite DFT methods such as r2SCAN-3c. They are implemented in our freely available standalone codes (dftd4, s-dftd3) and the D4 version will be also available in the upcoming ORCA 6.0 program package.

Modulating Hydrogen Shuttling in Ammonia by Neutral and Cationic Boron-Containing Frustrated Lewis Pairs (FLPs)

Xanthene-backbone FLPs featuring secondary borane functions -B(ArX)H (where ArX=C6F5 (ArF) or C6Cl5 (ArCl)) have been targeted through reactions of the dihydroboranes Me2S ⋅ BArXH2 with [4,5-xanth(PR2)Li]2 (R=Ph, iPr), and investigated in the synthesis of related cationic systems via hydride abstraction. The reactivity of these systems (both cationic and charge neutral) with ammonia have been probed, with a view to probing the potential for proton shuttling via N-H bond 'activation.' We find that in the case of four-coordinate boron systems (cationic or change neutral), the N-H linkage remains intact, supported by a NH⋅⋅⋅P hydrogen bond which is worth up to 17 kcal mol-1 thermodynamically, and enabled by planarization of the flexible xanthene scaffold. For cationic three coordinate systems, N-to-P proton transfer is viable, driven by the ability of the boron centre to stabilise the [NH2]- conjugate base through N-to-B π bonding. This proton transfer can be shown to be reversible in the presence of excess ammonia, depending on the nature of the B-bound ArX group. It is viable in the case of C6F5 substituents, but is prevented by the more sterically encumbering and secondary donor-stabilising capabilities of the C6Cl5 substituent.

Chemical accuracy for ligand-receptor binding Gibbs energies through multi-level SQM/QM calculations

Calculating the Gibbs energies of binding of ligand-receptor systems with a thermochemical accuracy of ± 1 kcal mol-1 is a challenge to computational approaches. After exploration of the conformational space of the host, ligand and their resulting complexes upon coordination by semi-empirical GFN2 MD and meta-MD simulations, the systematic refinement through a multi-level improvement of binding modes in terms of electronic energies and solvation is able to give Gibbs energies of binding of drug molecules to CB[8] and β-CD macrocyclic receptors with such an accuracy. The accurate treatment of a small number of structures outperforms system-specific force-matching and alchemical transfer model approaches without an extensive sampling and integration.

Transformation of Metal-Organic Framework from Kinetic to Thermodynamic Product for Controlled Delivery of Vitamin C

A biocompatible metal-organic framework (MOF), named HSTC-4, constructed using the flexible 4,4'-oxybis(benzoic acid) (OBA), was developed to enable efficient loading and controlled release of vitamin C (VC) through a combination of strategies involving ligand length, structure design, and metal selection. The kinetic product HSTC-4 demonstrates a propensity for transforming into the thermodynamically stable HSTC-5 under external stimuli, such as photoillumination and vacuum heating, as witnessed by single-crystal to single-crystal transformation. Density functional theory (DFT) calculations reveal that the VC guest molecules exhibit stronger binding affinity with HSTC-5 due to its narrower pores compared to HSTC-4, resulting in a slower release of VC from VC@HSTC-5. Furthermore, precise control over VC release can be achieved by introducing surface modifications involving SiO2 onto the structure of VC@HSCT-5, while simultaneously adjusting environmental factors such as pH and temperature conditions. Preliminary cell culture experiments and cytotoxicity assays highlight the biocompatibility of HSTC-5, suggesting that it is a promising platform for sustained drug delivery and diverse biomedical applications.

A nucleophilic beryllyl complex via metathesis at [Be-Be]2

Owing to its high toxicity, the chemistry of element number four, beryllium, is poorly understood. However, as the lightest elements provide the basis for fundamental models of chemical bonding, there is a need for greater insight into the properties of beryllium. In this context, the chemistry of the homo-elemental Be-Be bond is of fundamental interest. Here the ligand metathesis chemistry of diberyllocene (1; CpBeBeCp)-a stable complex with a Be-Be bond-has been investigated. These studies yield two complexes with Be-Be bonds: Cp*BeBeCp (2) and [K{(HCDippN)2BO}2]BeBeCp (3; Dipp = 2,6-diisopropylphenyl). Quantum chemical calculations indicate that the Be-Be bond in 3 is polarized to such an extent that the complex could be formulated as a mixed-oxidation state Be0/BeII complex. Correspondingly, it is demonstrated that 3 can transfer the 'beryllyl' anion, [BeCp]-, to an organic substrate, by analogy with the reactivity of sp2-sp3 diboranes. Indeed, this work reveals striking similarities between the homo-elemental bonding linkages of beryllium and boron, despite the respective metallic and non-metallic natures of these elements.

A geminal antimony(iii)/phosphorus(iii) frustrated Lewis pair

The geminal Lewis pair (F5C2)2SbCH2P(tBu)2 (1) was prepared by reacting (F5C2)2SbCl with LiCH2P(tBu)2. Despite its extremely electronegative pentafluoroethyl substituents, the neutral 1 exhibits a relatively soft acidic antimony function according to the HSAB concept (hard-soft acid-base). These properties lead to a reversibility in the binding of CS2 to 1, as observed by VT-NMR spectroscopy, while no reaction with CO2 is observed. The reaction behaviour towards heterocumulenes and the specific interaction situation in the CS2 adduct were analysed by quantum chemical calculations. The FLP-type reactivity of 1 has also been demonstrated by reaction with a variety of small molecules (SO2, PhNCO, PhNCS, (MePh2P)AuCl). The reactions of 1 with PhNCO and PhNCS led to different types of cyclic addition products: PhNCO adds with its N[double bond, length as m-dash]C bond and PhNCS adds preferentially with its C[double bond, length as m-dash]S bond. The reaction of 1 with (MePh2P)AuCl gave an adduct {[(F5C2)2SbCH2(tBu)2P]2Au}+ with a clamp-like structure binding a chloride anion by its two antimony atoms in chelate mode. Compound 1 and its adducts have been characterised by X-ray diffraction experiments, multinuclear NMR spectroscopy, elemental analyses and computational calculations (DFT, QTAIM, IQA).

All-visible-light-driven salicylidene schiff-base-functionalized artificial molecular motors

Light-driven rotary molecular motors are among the most promising classes of responsive molecular machines and take advantage of their intrinsic chirality which governs unidirectional rotation. As a consequence of their dynamic function, they receive considerable interest in the areas of supramolecular chemistry, asymmetric catalysis and responsive materials. Among the emerging classes of responsive photochromic molecules, multistate first-generation molecular motors driven by benign visible light remain unexplored, which limits the exploitation of the full potential of these mechanical light-powered systems. Herein, we describe a series of all-visible-light-driven first-generation molecular motors based on the salicylidene Schiff base functionality. Remarkable redshifts up to 100 nm in absorption are achieved compared to conventional first-generation motor structures. Taking advantage of all-visible-light-driven multistate motor scaffolds, adaptive behaviour is found as well, and potential application in multistate photoluminescence is demonstrated. These functional visible-light-responsive motors will likely stimulate the design and synthesis of more sophisticated nanomachinery with a myriad of future applications in powering dynamic systems.

Kinetics of sulfur-transfer from titanocene (poly)sulfides to sulfenyl chlorides: rapid metal-assisted concerted substitution

The kinetics of sulfur transfer from titanocene (poly)sulfides (RCp2TiS5, Cp2TiS4CMe2, Cp2Ti(SAr)2, Cp2TiCl(SAr)) to sulfenyl chlorides (S2Cl2, RSCl) have been investigated by a combination of stopped-flow UV-Vis/NMR reaction monitoring, titration assays, numerical kinetic modelling and KS-DFT calculations. The reactions are rapid, proceeding to completion over timescales of milliseconds to minutes, via a sequence of two S-S bond-forming steps (k 1, k 2). The archetypical polysulfides Cp2TiS5 (1a) and Cp2TiS4C(Me2) (2a) react with disulfur dichloride (S2Cl2) through rate-limiting intermolecular S-S bond formation (k 1) followed by a rapid intramolecular cyclization (k 2, with k 2 ≫ k 1 [RSCl]). The monofunctional sulfenyl chlorides (RSCl) studied herein react in two intermolecular S-S bond forming steps proceeding at similar rates (k 1 ≈ k 2). Reactions of titanocene bisthiophenolates, Cp2Ti(SAr)2 (5), with both mono- and di-functional sulfenyl chlorides result in rapid accumulation of the monothiophenolate, Cp2TiCl(SAr) (6) (k 1 > k 2). Across the range of reactants studied, the rates are relatively insensitive to changes in temperature and in the electronics of the sulfenyl chloride, moderately sensitive to the electronics of the titanocene (poly)sulfide (ρ (Ti-(SAr)) ≈ -2.0), and highly sensitive to the solvent polarity, with non-polar solvents (CS2, CCl4) leading to the slowest rates. The combined sensitivities are the result of a concerted, polarized and late transition state for the rate-limiting S-S bond forming step, accompanied by a large entropic penalty. Each substitution step {[Ti]-SR' + Cl-SR → [Ti]-Cl + RS-SR'} proceeds via titanium-assisted Cl-S cleavage to generate a transient pentacoordinate complex, Cl-[Cp2TiX]-S(R')-SR, which then undergoes rapid Ti-S dissociation.

Unravelling the kinetics of electro- and photochemical S → O linkage isomerization in Ru(II)-NHC-DMSO complexes utilised for photoinduced substitution reactions

Our recently reported Ru(III)-NHC complexes 1a and 1b were utilized as suitable precursors to prepare new Ru(II)-NHC-(DMSO)2 complexes 2a and 2b. Complexes 2a and 2b reacted with 2,2'-bipyridine to give complexes 3a and 3b, respectively, with substitution of only one of the DMSO ligands. All new complexes were characterized using various spectroscopic techniques and the molecular structures of 2a and 3a were determined using single-crystal X-ray diffraction technique. Complexes 2a, 2b, 3a, and 3b showed the S → O linkage isomerization of the DMSO ligand upon oxidation of the Ru centre from +II to +III, as confirmed by the thermodynamic and kinetic data obtained from cyclic voltammetry experiments. It was observed that in the bisdimethylsulfoxide complexes 2a and 2b, only one DMSO ligand isomerized, which was further corroborated by the computational studies performed to optimize the geometry of the possible linkage isomers of complexes 2a and 3a in +2 and +3 oxidation states, whereas complexes 3a and 3b showed a high preference for the O-bound isomer in the Ru(III) redox state. The role of NHC in stabilizing the mixed isomer in complexes 2a and 2b and preventing the isomerization of both DMSO ligands coordinated to the Ru centre was studied; moreover, NHC provided good solvent compatibility for photochemical S → O isomerization in all the complexes. Taking advantages of the photoinduced linkage isomerization in 2a and 2b, the synthesis of 3a and 3b was revisited and performed using 2a and 2b, respectively, following a photoinduced substitution reaction in the presence of 2,2'-bipyridine. The kinetics of the reversion from the O-bound to S-bound isomer was found to follow the DMSO-assisted intermolecular S → O isomerization pathway.

Assessment of the applicability of DFT methods to [Cp*Rh]-catalyzed hydrogen evolution processes

The present computational study provides a benchmark of density functional theory (DFT) methods in describing hydrogen evolution processes catalyzed by [Cp*Rh]-containing organometallic complexes. A test set was composed of 26 elementary reactions featuring chemical transformations and bonding situations essential for the field, including the emerging concept of non-innocent Cp* behavior. Reference values were obtained from a highly accurate 3/4 complete basis set and 6/7 complete PNO space extrapolated DLPNO-CCSD(T) energies. The performance of lower-level extrapolation procedures was also assessed. We considered 84 density functionals (DF) (including 13 generalized gradient approximations (GGA), nine meta-GGAs, 33 hybrids, and 29 double-hybrids) and three composite methods (HF-3c, PBEh-3c, and r2SCAN-3c), combined with different types of dispersion corrections (D3(0), D3BJ, D4, and VV10). The most accurate approach is the PBE0-DH-D3BJ (MAD of 1.36 kcal mol-1) followed by TPSS0-D3BJ (MAD of 1.60 kcal mol-1). Low-cost r2SCAN-3c composite provides a less accurate but much faster alternative (MAD of 2.39 kcal mol-1). The widely used Minnesota-family M06-L, M06, and M06-2X DFs should be avoided (MADs of 3.70, 3.94, and 4.01 kcal mol-1, respectively).

Oxidation States: Intrinsically Ambiguous?

The oxidation state ( OS ) formalism is a much-appreciated good in chemistry, receiving wide application. However, like all formalisms, limitations are inescapable, some of which have been recently explored. Providing a broader context, we discuss the OS and its interpretation from a computational perspective for transition metal (TM) complexes. We define a broadly applicable and easy-to-use procedure to derive OS s based on quantum chemical calculations, via the use of localized orbitals, dubbed the Intrinsic OS . Applying this approach to a cobalt complex in five OS s, isolated by Hunter and co-workers (Inorg. Chem.2021, 60, 17445), we find that the calculated Intrinsic OS matches the formal OS , consistent with the experimental characterization. Through analysis of the delocalized orbitals, the ligand field of the Co(III) complex is found to be "inverted", despite every cobalt-ligand bond being classically dative from the localized perspective-a bonding scenario very similar to that of [Cu(CF3)4]-. This is not atypical but rather a natural consequence of these metals bonding in the high-valent region, and we propose a more restrictive definition of (locally) inverted bonding. Additionally, two bonding descriptors within the Intrinsic Bonding Orbital (IBO) framework (σ-gain and π-loss) are introduced, which enable facile quantification of electron-sharing covalency across a broad range of TM complexes.

Tips and Tricks in the Modeling of Supramolecular Peptide Assemblies

Supramolecular peptide assemblies (SPAs) hold promise as materials for nanotechnology and biomedicine. Although their investigation often entails adapting experimental techniques from their protein counterparts, SPAs are fundamentally distinct from proteins, posing unique challenges for their study. Computational methods have emerged as indispensable tools for gaining deeper insights into SPA structures at the molecular level, surpassing the limitations of experimental techniques, and as screening tools to reduce the experimental search space. However, computational studies have grappled with issues stemming from the absence of standardized procedures and relevant crystal structures. Fundamental disparities between SPAs and protein simulations, such as the absence of experimentally validated initial structures and the importance of the simulation size, number of molecules, and concentration, have compounded these challenges. Understanding the roles of various parameters and the capabilities of different models and simulation setups remains an ongoing endeavor. In this review, we aim to provide readers with guidance on the parameters to consider when conducting SPA simulations, elucidating their potential impact on outcomes and validity.

Water training initiates spatially regulated microstructures with competitive mechanics in hydroadaptive polymers

The strategy using water as a medium for dynamic modulation of competitive plasticity and viscoelasticity provides a unique perspective to attain adaptive materials. We reveal sustainable polymers, herein cellulose phenoxyacetate as a typical example, with unusual water-responsive dual-mechanic functionalities addressed via a chronological water training strategy. The temporal significance of such water-responsive mechanical behaviors becomes apparent considering that a mere 3-minute exposure or a prolonged 3-hour exposure to water induced different types of mechano-responsiveness. This endows the materials with multiple recoverable shape-changes during water and air training, and consequently even underlines the switchability between the pre-loaded stable water shapes (> 20 months) and the sequentially fixed air shapes. Our discovery exploits the competitive mechanics initiated by water training, enabling polymers with spatially regulated microstructures via their inherently distinct mechanical properties. Insights into the molecular changes represents a considerable fundamental innovation, can be broadly applicable to a diverse array of hydroadaptive polymers.

Antiaromaticity of Cycloheptatrienyl Anions: Structure, Acidity, and Magnetic Properties

Investigations of the nature and degree of antiaromaticity of cycloheptatrienyl anion derivatives using both experimental and computational tools are presented. The ground state of cycloheptatrienyl anion in the gas phase is triplet, planar and Baird-aromatic. In DMSO, it assumes a singlet distorted allylic form with a paratropic ring current. The other derivatives in both phases assume either allylic or diallylic conformations depending on the substituent pattern. A combination of experimental and computational methods was used to determine the pKa values of 16 derivatives in DMSO, which ranged from 36 to -10.7. We revealed that the stronger stabilization of the anionic system, which correlates with acidity, does not necessarily imply a lower degree of antiaromaticity in terms of magnetic properties. Conversely, the substitution pattern first affects the geometry of the ring through the bulkiness of the substituents and their better conjugation with a more distorted system. Consequently, the distortion reduces the cyclic conjugation in the π-system and thereby decreases the paratropic current in a magnetic field, which manifests itself as a decrease in the NICS. The triplet-state geometries and magnetic properties are nearly independent on the substitution pattern, which is typical for simple aromatic systems.

Fresh Impetus in the Chemistry of Calcium Peroxides

Redox-inactive metal ions are essential in modulating the reactivity of various oxygen-containing metal complexes and metalloenzymes, including photosystem II (PSII). The heart of this unique membrane-protein complex comprises the Mn4CaO5 cluster, in which the Ca2+ ion acts as a critical cofactor in the splitting of water in PSII. However, there is still a lack of studies involving Ca-based reactive oxygen species (ROS) systems, and the exact nature of the interaction between the Ca2+ center and ROS in PSII still generates intense debate. Here, harnessing a novel Ca-TEMPO complex supported by the β-diketiminate ligand to control the activation of O2, we report the isolation and structural characterization of hitherto elusive Ca peroxides, a homometallic Ca hydroperoxide and a heterometallic Ca/K peroxide. Our studies indicate that the presence of K+ cations is a key factor controlling the outcome of the oxygenation reaction of the model Ca-TEMPO complex. Combining experimental observations with computational investigations, we also propose a mechanistic rationalization for the reaction outcomes. The designed approach demonstrates metal-TEMPO complexes as a versatile platform for O2 activation and advances the understanding of Ca/ROS systems.

Pushing the Limits of Organometallic Redox Chemistry with an Isolable Mn(-I) Dianion

We report an incredibly reducing and redox-active Mn-I dianion, [Mn(CO)3(Ph2B(tBuNHC)2)]2- (NHC = N-heterocyclic carbene), furnished via 2e- reduction of the parent 16e- MnI complex with Na0 or K0. Cyclic voltammograms show a Mn0/-I redox couple at -3.13 V vs Fc+/0 in tetrahydrofuran (THF), -3.06 V in 1,2-dimethoxyethane, and -2.85 V in acetonitrile. The diamagnetic Mn-I dianion is stable in solution and solid-state at room temperature, tolerating a wide range of countercations ([M(2.2.2)crypt]+, [M(18-crown-6)]+, [nBu4N]+; M = Na, K). Countercation identity does not significantly alter 13C NMR spectral signatures with [nBu4N]+ and Na+, suggesting minimal ion pairing in solution. IR spectroscopy reveals a significant decrease in CO stretching frequencies from MnI to Mn-I (ca. 240 cm-1), consistent with a drastic increase in electron density at Mn. State-of-the-art DFT calculations are in excellent agreement with the observed IR spectral data. Moreover, the Mn-I dianion behaves as a chemical reductant, smoothly releasing 1e- or 2e- to regenerate the oxidized Mn0 or MnI species in solution. The reducing potential of [Mn(CO)3(Ph2B(tBuNHC)2)]2- surpasses the naphthalenide anion in THF (-3.09 V) and represents one of the strongest isolable chemical redox agents.

Redox-Switchable Aromaticity in a Helically Extended Indeno[2,1-c]fluorene

Molecular switches have received major attention to enable the reversible modulation of various molecular properties and have been extensively used as trigger elements in diverse fields, including molecular machines, responsive materials, and photopharmacology. Antiaromaticity is a fascinating property that has attracted not only significant fundamental interest but is also increasingly relevant in different applications, in particular organic (opto)electronics. However, designing systems in which (anti)aromaticity can be judiciously and reversibly switched ON and OFF remains challenging. Herein, we report a helicene featuring an indenofluorene-bridged bisthioxanthylidene as a novel switch wherein a simultaneous two-electron (electro)chemical redox process allows highly reversible modulation of its (anti)aromatic character. Specifically, the two thioxanthylidene rotors, attached to the initially aromatic indenofluorene scaffold via overcrowded alkenes, adopt an anti-folded structure, which upon oxidation convert to singly bonded, twisted conformations. This is not only associated with significant (chir)optical changes but importantly also results in formation of the fully conjugated, formally antiaromatic as-indacene motif in the helical core of the switch. This process proceeds without the buildup of radical cation intermediates and thus enables highly reversible switching of molecular geometry, aromaticity, absorbance, and chiral expression under ambient conditions, as evidenced by NMR, UV-vis, CD, and (spectro)electrochemical analyses, supported by DFT calculations. We expect this concept to be extendable to a wide range of robust antiaromatic-aromatic switches and to provide a basis for modulation of the structure and properties of these fascinating inherently chiral polycyclic π-scaffolds.

pKalculator: A pK a predictor for C-H bonds

Determining the pK a values of various C-H sites in organic molecules offers valuable insights for synthetic chemists in predicting reaction sites. As molecular complexity increases, this task becomes more challenging. This paper introduces pKalculator, a quantum chemistry (QM)-based workflow for automatic computations of C-H pK a values, which is used to generate a training dataset for a machine learning (ML) model. The QM workflow is benchmarked against 695 experimentally determined C-H pK a values in DMSO. The ML model is trained on a diverse dataset of 775 molecules with 3910 C-H sites. Our ML model predicts C-H pK a values with a mean absolute error (MAE) and a root mean squared error (RMSE) of 1.24 and 2.15 pK a units, respectively. Furthermore, we employ our model on 1043 pK a-dependent reactions (aldol, Claisen, and Michael) and successfully indicate the reaction sites with a Matthew's correlation coefficient (MCC) of 0.82.

Organocatalytic asymmetric synthesis of Si-stereogenic silacycles

A strong and confined Brønsted acid catalyzed enantioselective cyclization of bis(methallyl)silanes provides enantioenriched Si-stereogenic silacycles. High enantioselectivities of up to 96.5:3.5 er were obtained for a range of bis(methallyl)silanes. NMR and ESI-MS studies reveal that the formation of a covalent adduct irreversibly inhibits turnover. Remarkably, we found that acetic acid as an additive promotes the collapse of this adduct, enabling full turnover. Experimental investigation and density functional theory (DFT) calculations were conducted to elucidate the origin of this phenomenon and the observed enantioselectivity.

Remote Functionalization by Pd-Catalyzed Isomerization of Alkynyl Alcohols

In recent years, progress has been made in the development of catalytic methods that allow remote functionalizations based on alkene isomerization. In contrast, protocols based on alkyne isomerization are comparatively rare. Herein, we report a general Pd-catalyzed long-range isomerization of alkynyl alcohols. Starting from aryl-, heteroaryl-, or alkyl-substituted precursors, the optimized system provides access preferentially to the thermodynamically more stable α,β-unsaturated aldehydes and is compatible with potentially sensitive functional groups. We showed that the migration of both π-components of the carbon-carbon triple bond can be sustained over several methylene units. Computational investigations served to shed light on the key elementary steps responsible for the reactivity and selectivity. These include an unorthodox phosphine-assisted deprotonation rather than a more conventional β-hydride elimination in the final tautomerization event.

Bartender: Martini 3 Bonded Terms via Quantum Mechanics-Based Molecular Dynamics

Coarse-grained (CG) molecular dynamics (MD) simulations have grown in applicability over the years. The recently released version of the Martini CG force field (Martini 3) has been successfully applied to simulate many processes, including protein-ligand binding. However, the current ligand parametrization scheme is manual and requires an a priori reference all-atom (AA) simulation for benchmarking. For systems with suboptimal AA parameters, which are often unknown, this translates into a CG model that does not reproduce the true dynamical behavior of the underlying molecule. Here, we present Bartender, a quantum mechanics (QM)/MD-based parametrization tool written in Go. Bartender harnesses the power of QM simulations and produces reasonable bonded terms for Martini 3 CG models of small molecules in an efficient and user-friendly manner. For small, ring-like molecules, Bartender generates models whose properties are indistinguishable from the human-made models. For more complex, drug-like ligands, it is able to fit functional forms beyond simple harmonic dihedrals and thus better captures their dynamical behavior. Bartender has the power to both increase the efficiency and the accuracy of Martini 3-based high-throughput applications by producing numerically stable and physically realistic CG models.

Cobalt-Mediated Photochemical C-H Arylation of Pyrroles

Precious metal complexes remain ubiquitous in photoredox catalysis (PRC) despite concerted efforts to find more earth-abundant catalysts and replacements based on 3d metals in particular. Most otherwise plausible 3d metal complexes are assumed to be unsuitable due to short-lived excited states, which has led researchers to prioritize the pursuit of longer excited-state lifetimes through careful molecular design. However, we report herein that the C-H arylation of pyrroles and related substrates (which are benchmark reactions for assessing the efficacy of photoredox catalysts) can be achieved using a simple and readily accessible octahedral bis(diiminopyridine) cobalt complex, [1-Co](PF6)2. Notably, [1-Co]2+ efficiently functionalizes both chloro- and bromoarene substrates despite the short excited-state lifetime of the key photoexcited intermediate *[1-Co]2+ (8 ps). We present herein the scope of this C-H arylation protocol and provide mechanistic insights derived from detailed spectroscopic and computational studies. These indicate that, despite its transient existence, reduction of *[1-Co]2+ is facilitated via pre-assembly with the NEt3 reductant, highlighting an alternative strategy for the future development of 3d metal-catalyzed PRC.

Anion-Coordination Foldamer-Based Polymer Network: from Molecular Spring to Elastomer

Foldamer is a scaled-down version of coil spring, which can absorb and release energy by conformational change. Here, polymer networks with high density of molecular springs were developed by employing anion-coordination-based foldamers as the monomer. The coiling of the foldamer is controlled by oligo(urea) ligands coordinating to chloride ions; subsequently, the folding and unfolding of foldamer conformations endow the polymer network with excellent energy dissipation and toughness. The mechanical performance of the corresponding polymer networks shows a dramatic increase from P-L2UCl (non-folding), to P-L4UCl (a full turn), and then to P-L6UCl (1.5 turns), in terms of strength (2.62 MPa; 14.26 MPa; 22.93 MPa), elongation at break (70 %; 325 %; 352 %), Young's modulus (2.69 MPa; 63.61 MPa; 141.50 MPa), and toughness (1.12 MJ/m3; 21.39 MJ/m3; 49.62 MJ/m3), respectively, which is also better than those without anion centers and the non-foldamer based counterparts. Moreover, P-L6UCl shows enhanced strength and toughness than most of the molecular-spring based polymer networks. Thus, an effective strategy for designing high-performance anion-coordination-based materials is presented.

The synthesis of N,1,4-tri(alkoxy-hydroxybenzyl)-1,4-diazepane-amines: investigations on reaction characteristics and mechanism

1,4-Diazepane-6-amine (DAZA) can be alkylated with three 2-hydroxybenzyl pendant arms, resulting in hexadentate chelators suitable for coordination of radiometals like 68Ga. These chelators, N,1,4-tri(alkoxy-2-hydroxybenzyl)-DAZA, can be produced via a one-pot synthesis, with the first step being a carbonyl amine condensation of DAZA with two respective 4-alkoxy-2-hydroxybenzaldehydes, followed by reductive amination with sodium borohydride. While the first step of this reaction is predictable, the subsequent reductive amination can result in either mono-, di- or tri(alkoxy-hydroxybenzyl)-DAZA compounds. Seeking to identify dependencies that might allow a specific reaction control towards the formation of either of the three possible products, and particularly towards the favoured trialkylated DAZA compounds, a variety of synthesis trials were performed. Additionally, computational methods were employed to evaluate the underlying reaction mechanism. Synthesis trials verified that the trialkylated DAZA compounds are formed via direct reductive amination of the dialkylated DAZA compounds. Subsequently, a synthetic method was established, leading to an increase in the percentage of the trialkylated DAZA compounds, which allowed the successful isolation of those hexadentate chelators. Additionally, an alternative pathway proceeding via aminal C-N bond insertion of an attacking third carbonyl moiety was evaluated by means of quantum chemical calculations but so far remains entirely hypothetical.

Formal Single Atom Editing of the Glycosylated Natural Product Fidaxomicin Improves Acid Stability and Retains Antibiotic Activity

Fidaxomicin (Fdx) constitutes a glycosylated natural product with excellent antibacterial activity against various Gram-positive bacteria but is approved only for Clostridioides difficile infections. Poor water solubility and acid lability preclude its use for other infections. Herein, we describe our strategy to overcome the acid lability by introducing acid-stable S-linked glycosides. We describe the direct, diastereoselective modification of unprotected Fdx without the need to avoid air or moisture. Using our newly established approach, Fdx was converted to the single atom exchanged analogue S-Fdx, in which the acid labile O-glycosidic bond to the noviose sugar was replaced by the acid stable S-glycosidic bond. Studies of the antibacterial activity of a structurally diverse set of thioglycoside derivatives revealed high potency of acyl derivatives of S-Fdx against Clostridioides difficile (MIC range: 0.12-4 μg/mL) and excellent potency against Clostridium perfringens (MIC range: 0.06-0.5 μg/mL).

Utilization of 13C NMR Carbon Shifts for the Attribution of Diastereomers in Methyl-Substituted Cyclohexanes

In this report, we address the challenge of assigning diastereomers for methyl cyclohexanes, particularly those with quaternary centers, which remains nontrivial despite modern NMR techniques. By utilizing a HSQC NMR experiment to identify methyl-carbons coupled with a simple conformational analysis, we identified an effective and quite general method for assigning stereochemistry, even in cases where diastereomeric mixtures are inseparable.

Insight into the nature of carbon-metal bonding for N-heterocyclic carbenes in gold/silver complexes and nanoparticles using DFT-correlated Raman spectroscopy: strong evidence for π-backbonding

N-Heterocyclic carbenes (NHCs) have emerged as promising ligands for stabilizing metallic complexes, nanoclusters, nanoparticles (NPs) and surfaces. The carbon-metal bond between NHCs and metal atoms plays a crucial role in determining the resulting material's stability, reactivity, function, and electronic properties. Using Raman spectroscopy coupled with density functional theory calculations, we investigate the nature of carbon-metal bonding in NHC-silver and NHC-gold complexes as well as their corresponding NPs. While low wavenumbers are inaccessible to standard infrared spectroscopy, Raman detection reveals previously unreported NHC-Au/Ag bond-stretching vibrations between 154-196 cm-1. The computationally efficient r2SCAN-3c method allows an excellent correlation between experimental and predicted Raman spectra which helps calibrate an accurate description of NHC-metal bonding. While π-backbonding should stabilize the NHC-metal bond, conflicting reports for the presence and absence of π-backbonding are seen in the literature. This debate led us to further investigate experimental and theoretical results to ultimately confirm and quantify the presence of π-backbonding in these systems. Experimentally, an observed decrease in the NHC's CN stretching due to the population of the π* orbital is a good indication for the presence of π-backbonding. Using energy decomposition analysis - natural orbitals for chemical valence (EDA-NOCV), our calculations concur and quantify π-backbonding in these NHC-bound complexes and NPs. Surprisingly, we observe that NPs are less stabilized by π-backbonding compared to their respective complexes-a result that partially explains the weaker NHC-NP bond. The protocol described herein will help optimize metal-carbon bonding in NHC-stabilized metal complexes, nanoparticles and surfaces.

Combining low-cost electronic structure theory and low-cost parallel computing architecture

The computational efficiency of low-cost electronic structure methods can be further improved by leveraging heterogenous computing architectures. The software package TeraChem has been developed since 2008 to make use of graphical processing units (GPUs), particularly their strong single-precision performance, for the acceleration of quantum chemical calculations. Here, we present the implementation of three low-cost methods, namely HF-3c, PBEh-3c, and the recently introduced ωB97X-3c. We show that these can benefit in terms of performance when combined with "consumer grade" GPUs by leveraging the mixed precision integral handling in TeraChem. The current limitation of the latter's GPU integral library is that Gaussian integrals only for functions with angular momentum l < 3 can be computed, which generally restricts the achievable accuracy in terms of the one-particle basis set. Particularly, the implementation of the ωB97X-3c method now enables higher accuracy with this setting which, in turn, provides the most efficient implementation accessible with consumer-grade hardware. We furthermore show that the implemented 3c methods can be combined with the hh-TDA formalism. This gives new and efficient low-cost multi-configurational excited states methods, which are benchmarked for the description of lowest vertical excitation energies in this work. All in all, the combination of these efficient electronic structure theory methods with affordable highly parallelized computing hardware provides an optimal computational and monetary cost to accuracy ratio.

Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

Enhancing Molecular Energy Predictions with Physically Constrained Modifications to the Neural Network Potential

Exclusively prioritizing the precision of energy prediction frequently proves inadequate in satisfying multifaceted requirements. A heightened focus is warranted on assessing the rationality of potential energy curves predicted by machine learning-based force fields (MLFFs), alongside evaluating the pragmatic utility of these MLFFs. This study introduces SWANI, an optimized neural network potential stemming from the ANI framework. Through the incorporation of supplementary physical constraints, SWANI aligns more cohesively with chemical expectations, yielding rational potential energy profiles. It also exhibits superior predictive precision compared with that of the ANI model. Additionally, a comprehensive comparison is conducted between SWANI and a prominent graph neural network-based model. The findings indicate that SWANI outperforms the latter, particularly for molecules exceeding the dimensions of the training set. This outcome underscores SWANI's exceptional capacity for generalization and its proficiency in handling larger molecular systems.

Radical and diradical states of bis(molybdenocene dithiolene) complexes

The synthesis and characterization of two bis(dithiolene) proligands involving heteroatomic linkers such as 1,4-dithiine and dihydro-1,4-disiline between the two protected dithiolene moieties are described. Two bimetallic complexes involving these heteroatomic bridges between two redox active bis(cyclopentadienyl)molybdenum dithiolene moieties have been synthesized and characterized by electrochemistry, spectroelectrochemistry, and their properties rationalized with (TD-)DFT. Cyclic voltammetry experiments show sequential oxidation of the two redox centers with ΔE values between successive one-electron transfers varying according to the nature of the bridge. Depending on the nature of the heteroatomic bridge, the bis-oxidized complexes exhibit either a diradical character with both radicals essentially localized on the metallacycles, or a closed-shell dicationic state.

Catalytic Upgrading of Acetaldehyde to Acetoin Using a Supported N-Heterocyclic Carbene Catalyst

We report the catalytic synthesis of 3-hydroxy-2-butanon (acetoin) from acetaldehyde as a key step in the synthesis of C4-molecules from ethanol. Facile C-C bond formation at the α-carbon of the C2 building block is achieved using an N-heterocyclic carbene (NHC) catalyst. The immobilization of the catalyst on a Merrifield's peptide resin and its spectroscopic characterisation using solid-state Nuclear Magnetic Resonance (NMR) is described herein. The immobilization of the NHC catalyst allows for process intensification steps and the reported catalytic system was subjected to batch recycling as well as continuous flow experiments. The robustness of the catalytic system was shown over a maximum of 10 h time-on-stream. Overall, high selectivity S>90 % was observed. The observed deactivation of the catalyst with increasing time-on-stream is explained by ex-situ 1H solution-state, as well as 13C and 15N solid-state NMR spectra allowing us to develop a deeper understanding of the underlying decomposition mechanism of the catalyst.

Synthesis and Characterization of a Copper Dinitrogen Complex Supported by a Weakly Coordinating Anion

We report the synthesis and full characterization of the copper dinitrogen complex [(η1-N2)Cu{Al(ORF)4}] 2 (RF=C(CF3)3) prepared by a cascade metathesis reaction of Ag[Al(ORF)4] with CuI-excess in iso-perfluorohexane (i-pfh) under N2 atmosphere. Title compound 2 features an extraordinarily high N2 stretching frequency at 2313/2314 cm-1 (IR/Raman) and was characterized by single-crystal and powder X-ray diffractometry. Quantum chemical charge displacement analysis based on natural orbitals of chemical valence (CD-NOCV) indicates that the copper-dinitrogen interaction is still governed by weak π-backdonation, but is significantly reduced compared to all literature-known transition metal dinitrogen complexes.

Useful synthetic artifacts? The impact of ubiquitous linker-adjacent groups on the self-assembly of discotic dimers

Although discotic dimers commonly feature bulky ether substituents adjacent to the linking group, the impact of these chains on self-assembly remains unclear. A series of dibenzo[a,c]phenazine dimers with alkoxy groups ortho to the linker were prepared and their solution conformational dynamics and liquid crystalline properties examined. The presence of a methoxy substitutent adjacent to the bridging group increased the phase stability, whereas longer chains dramatically decreased clearing temperatures. NMR solution studies indicated that adjacent groups increased the preference of dimers to adopt unfolded conformers. DFT models indicated that the unfolded structures were nonplanar and hence less compatible with columnar ordering, leading to a destabilization of the mesophases.

Quantitative predictive theories through integrating quantum, statistical, equilibrium, and nonequilibrium thermodynamics

Today's thermodynamics is largely based on the combined law for equilibrium systems and statistical mechanics derived by Gibbs in 1873 and 1901, respectively, while irreversible thermodynamics for nonequilibrium systems resides essentially on the Onsager Theorem as a separate branch of thermodynamics developed in 1930s. Between them, quantum mechanics was invented and quantitatively solved in terms of density functional theory (DFT) in 1960s. These three scientific domains operate based on different principles and are very much separated from each other. In analogy to the parable of the blind men and the elephant articulated by Perdew, they individually represent different portions of a complex system and thus are incomplete by themselves alone, resulting in the lack of quantitative agreement between their predictions and experimental observations. Over the last two decades, the author's group has developed a multiscale entropy approach (recently termed as zentropy theory) that integrates DFT-based quantum mechanics and Gibbs statistical mechanics and is capable of accurately predicting entropy and free energy of complex systems. Furthermore, in combination with the combined law for nonequilibrium systems presented by Hillert, the author developed the theory of cross phenomena beyond the phenomenological Onsager Theorem. The zentropy theory and theory of cross phenomena jointly provide quantitative predictive theories for systems from electronic to any observable scales as reviewed in the present work.

Limiting factors in the accuracy of DFT calculation for redox potentials

In the present study, we have investigated factors affecting the accuracy of computational chemistry calculation of redox potentials, namely the gas-phase ionization energy (IE) and electron affinity (EA), and the continuum solvation effect. In general, double-hybrid density functional theory methods yield IEs and EAs that are on average within ~0.1 eV of our high-level W3X-L benchmark, with the best performing method being DSD-BLYP/ma-def2-QZVPP. For lower-cost methods, the average errors are ~0.2-0.3 eV, with ωB97X-3c being the most accurate (~0.15 eV). For the solvation component, essentially all methods have an average error of ~0.3 eV, which shows the limitation of the continuum solvation model. Curiously, the directly calculated redox potentials show errors of ~0.3 eV for all methods. These errors are notably smaller than what can be expected from error propagation with the two components (IE and EA, and solvation effect). Such a discrepancy can be attributed to the cancellation of errors, with the lowest-cost GFN2-xTB method benefiting the most, and the most accurate ωB97X-3c method benefiting the least. For organometallic species, the redox potentials show large deviations exceeding ~0.5 eV even for DSD-BLYP. The large errors are attributed to those for the gas-phase IEs and EAs, which represents a major barrier to the accurate calculation of redox potentials for such systems.

In Search of the Best Low-Cost Methods for Efficient Screening of Conformers

Locating the lowest energy conformer is crucial for the accurate computation of equilibrium properties of molecular systems. This paper examines the performance of efficient low-cost methods in terms of the alignment and relative energies of their energy minima against the benchmark revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ potential energy surface. The low-cost methods considered include GFN-FF, GFN2-xTB, DFTB3, HF-3c, B97-3c, PBEh-3c, and r2SCAN-3c composite methods against a diverse test set of 20 compounds including alkanes, perfluoroalkyl molecules, peptides, open-shell radicals, and Zn(II) complexes of varying sizes. The "3c" composite methods are generally more accurate, but are at least 2-3 orders of magnitude more expensive than tight-binding methods which have energy minima that align well with the benchmark potential energy surface. The findings of this paper were further exploited to introduce a simple strategy involving Grimme's CENSO energy-sorting algorithm that resulted in up to an order of magnitude reduction in computational time for locating the lowest energy conformer on the revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ surface.

Three-component cascade reaction of 3-ketonitriles, 2-unsubstituted imidazole N-oxides, and aldehydes

A three-component condensation of 2-unsubstituted imidazole N-oxides, 3-ketonitriles, and aldehydes is described. The reaction proceeds via sequential Knoevenagel condensation/Michael addition under mild, catalyst-free conditions with various substrates. Furthermore, the corresponding 2-functionalized imidazole N-oxides can be further dehydrated to (Z)-2-aroyl-3-(1H-imidazol-2-yl)-acrylonitriles, which may also be directly prepared by changing the reaction conditions as a cascade of Knoevenagel condensation/Michael addition/dehydration.