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.

Conformational analysis and spectroscopic properties of antichagasic nifurtimox

Considering the health relevance of Chagas' disease, recent research efforts have focused on developing more efficient drug delivery systems containing nifurtimox (NFX). This paper comprehensively investigates NFX through conformational analysis and spectroscopic characterization. Using a conformer-rotamer ensemble sampling tool (CREST-xtb), five distinct conformers of NFX were sampled within a 3.0 kcal mol-1 relative energy window. Subsequently, such structures were used as inputs for geometry optimization by density functional theory (DFT) at B3LYP-def2-TZVP level of theory. Notably, harmonic vibrational frequencies were calculated to establish an in-depth comparison with experimental results and existing literature for the NFX or similar molecules and functional groups, thereby achieving a widely reasoned assignment of the mid-infrared band absorptions for the first time. Moreover, UV-VIS spectra of NFX were obtained in several solvents, enabling the determination of the molar absorptivity coefficient for the two electronic transitions observed for NFX. Among the aprotic solvents, a bathochromic effect was observed in the function of the dielectric constants. Furthermore, a hypochromic effect was observed when the drug was dissolved in protic solvents. These findings offer crucial support for new drug delivery systems containing NFX while demonstrating the potential of spectrophotometric studies in establishing quality control assays for NFX drug products.

Carbon-Rich Plasma-Deposited Silicon Oxycarbonitride Films Derived from 4-(Trimethylsilyl)morpholine as a Novel Single-Source Precursor

4-(trimethylsilyl)morpholine O(CH2CH2)2NSi(CH3)3 (TMSM) was investigated as a single-source precursor for SiCNO films synthesis. Optical emission spectroscopy of plasma generated from TMSM/He, TMSM/H2, and TMSM/NH3 gas mixtures revealed the presence of N2, CH, H, CN, and CO species. The last two are suggested to be responsible for the lowering of carbon concentration in the films in comparison with the precursor. The refractive index ranged from 1.5 to 2.0, and bandgap varied from 2.0 to 4.6 eV, which pointed that some of the films can be used as antireflective coatings in silicon photovoltaic cell technologies and dielectric layers in electronic devices.

Migration of Condensed Aromatic Hydrocarbons During Alkyne-Vinylidene Rearrangements

Diarylacetylenes ArC≡CAr featuring condensed aromatic hydrocarbon fragments (Ar) such as naphthalene, anthracene, phenanthrene and pyrene were converted into vinylidene ligands by 1,2-migration reactions within the coordination sphere of half-sandwich complexes [MII(dppe)Cp]+ (MII = RuII, FeII). Comparison of the extent of conversion of the alkyne substrates to the vinylidene complexes [Ru{=C=CAr2}(dppe)Cp]+ with those obtained from acetylenes functionalized by smaller groups (H, CH3, Ph) show that the molecular volume (VM) of the migrating group and relief of steric congestion plays a role during the rearrangement process. Conversely, the H-atoms from the larger condensed ring aryl groups that are in close proximity to the migrating sites also have a significant influence on the efficacy and extent of the reaction by restricting access of the alkyne to the metal center, resulting in a less effective migration reaction. This combination of competing steric factors (acceleration due to relief of steric congestion and restricted access of the alkyne moiety to the reaction site) is exemplified by the facile migration of 1-pyryl entities and the low yields of vinylidene products formed from 1,2-bis(9-anthryl)acetylene.

Sequential nucleophilic aromatic substitutions on cyanuric chloride: synthesis of BODIPY derivatives and mechanistic insights

Herein we report a study on the sequential substitution of different nucleophiles on cyanuric chloride to obtain potential candidates for metal sensors (5a-c). The set of nucleophiles on the 1,3,5-triazine ring includes a phenolic BODIPY, an aminoalkyl pyridine and aminoalkyl phosphoramidates, each one designed to play a specific role in the final fluoroionophore. Three new triazine triads were synthesized in similar yields: 5a (45%), 5b (43%) and 5c (52%) after a methodical sequential combination of the nucleophiles via thermodependent nucleophilic aromatic substitution of the three chlorine atoms of cyanuric chloride. To ratify the synthetic results we simulated the reaction mechanisms for the different nucleophiles, aiming to address the distinctive orthogonality and temperature control inherent in this process, identifying and providing a sound rationale for any preferential sequence of nucleophiles inserted into the triazine core. According to our experimental and computational analysis (thermo- and kinetic preferences), we have identified the following preferential order for the sequential substitution: p-hydroxybenzaldehyde > 2-(pyridin-2-yl)ethanamine > aminoalkyl phosphoramidate, indicating that all steps follow a single-step process (concerted) in two stages, where nucleophilic addition precedes leaving group dissociation. The Meisenheimer σ-complex was identified as a transition state structure, with insufficient stability to exist as an intermediate. We observed a consistent and progressive increase in barrier height: 2-8 kcal mol-1 for the first step, 9-15 kcal mol-1 for the second step, and >15 kcal mol-1 for the third substitution. These findings align with the experimental observation of thermodependency in the sequential substitution.

Tricyclic 1,4-Diphosphinines: Local vs. Global Aromaticity

1D and 2DNICS π . zz SOM ${{{\rm { NICS}}}_{{\rm { {\rm \pi}{}}}{\rm { .zz}}}^{{\rm { SOM}}}}$ (Sigma only model) calculations were performed on recently established tricyclic 1,4-diphosphinines as well as related benzene and pyrazine derivatives. The study was extended to evaluate the effect of the fused rings on the overall aromatic properties with a special focus on functional groups such as carbenes. The effect of non-aromatic heterocycles on the local ring current of the central ring is small, while aromatic heterocycles (e. g. NHC, imidazolium) lead to a global aromaticity. A higher sulfur content of the adjacent five-membered rings reduces the central ring current. The comparison to related tricyclic benzene and pyrazine derivatives showed that the 1,4-diphosphinine systems resemble more closely the situation in the benzene derivatives than the pyrazines. The effect of charged systems was studied using bis(TTF)-fused 1,4-diphosphinines and, according toNICS π . zz SOM ${{{\rm { NICS}}}_{{\rm { {\rm \pi}{}}}{\rm { .zz}}}^{{\rm { SOM}}}}$ values, the neutral form doesn't possess significant aromaticity but the tetracation resembles the global aromatic situation observed for other heterocycles in this study.

Unveiling the Antioxidative Potential of Galangin: Complete and Detailed Mechanistic Insights through Density Functional Theory Studies

A comprehensive quantum mechanical investigation delved into the antioxidative activity of galangin (Glg). Thermochemical and kinetic data were used to assess antiradical, chelating, and renewal potential under physiological conditions. A brief comparison with reference antioxidants and other flavonoids characterized Glg as a moderate antioxidative agent. The substance showed significantly lower performance in lipid compared to aqueous solvent─the reaction rates for scavenging •OOH in both media were established at 3.77 × 103 M-1 s-1 and 6.21 × 104 M-1 s-1, respectively, accounting for the molar fraction of both interacting molecules at the given pH. The impact of pH value on the kinetics was assessed. Although efficient at chelating Cu(II) ions, the formed complexes can still undergo the Fenton reaction. On the other hand, they persistently scavenge •OH in statu nascendi. The flavonoid effectively repairs oxidatively damaged biomolecules except model lipid acids. All Glg radicals are readily restored by physiologically prevailing O2•-. Given this, the polyphenol is expected to participate in antiradical and regenerating activities multiple times, amplifying its antioxidative potential.

On the origin of optical rotation changes during the κ-carrageenan disorder-to-order transition

It is well established that solutions of both polymeric and oligomeric κ-carrageenan exhibit a clear change in optical rotation (OR), in concert with gel-formation for polymeric samples, as the solution is cooled in the presence of certain ions. The canonical interpretation - that this OR change reflects a 'coil-to-helix transition' in single chains - has seemed unambiguous; the solution- or 'disordered'-state structure has ubiquitously been assumed to be a 'random coil', and the helical nature of carrageenan in the solid-state was settled in the 1970s. However, recent work has found that κ-carrageenan contains substantial helical secondary structure elements in the disordered-state, raising doubts over the validity of this interpretation. To investigate the origins of the OR, density-functional theory calculations were conducted using atomic models of κ-carrageenan oligomers. Changes were found to occur in the predicted OR owing purely to dimerization of chains, and - together with the additional effects of slight changes in conformation that occur when separated helical chains form double-helices - the predicted OR changes are qualitatively consistent with experimental results. These findings contribute to a growing body of evidence that the carrageenan 'disorder-to-order' transition is a cooperative process, and have further implications for the interpretation of OR changes demonstrated by macromolecules in general.

Evaluation of DFT methods for predicting geometries and NMR spectra of Bi(III) dithiocarbamate complexes with antitumor properties

CONTEXT

Bismuth complexes with dithiocarbamate ligands have attracted attention because of their biological applications, such as antimicrobial, antileishmanial, and anticancer properties. These complexes have high cytotoxic activity against cancer cells, being more active than the standard drugs cisplatin, doxorubicin, and tamoxifen. In the present study, we investigated the ability of some DFT methods to reproduce the geometries and NMR spectra of the Bi(III) dithiocarbamate complexes, selected based on their proven antitumor activity. Our investigation revealed that the M06-L/def2-TZVP/ECP/CPCM method presented good accuracy in predicting geometries, while the TPSSh/def2-SVP/ECP/CPCM method proved effective in analyzing the 13C NMR spectra of these molecules. In general, all examined methods exhibited comparable performance in predicting 1H NMR signals.

METHODS

Calculations were performed with the Gaussian 09 program using the def2-SVP and def2-TZVP basis sets, employing relativistic effective core potential (ECP) for Bi and using the CPCM solvent model. The exchange-correlation functionals BP86, PBE, OLYP, M06-L, B3LYP, B3LYP-D3, M06-2X, TPSSh, CAM-B3LYP, and ωB97XD were used in the study. Geometry optimizations were started from crystallographic structures available at the Cambridge Structural Database. The theoretical results were compared with experimental data using the mean root-mean-square deviation (RMSD), mean absolute deviations (MAD), and linear correlation coefficient (R2).

Predicting valence tautomerism in diverse cobalt-dioxolene complexes: elucidation of the role of ligands and solvent

The ability of molecular switches to reversibly interconvert between different forms promises potential applications at the scale of single molecules up to bulk materials. One type of molecular switch comprises cobalt-dioxolene compounds that exhibit thermally-induced valence tautomerism (VT) interconversions between low spin Co(iii)-catecholate (LS-CoIII-cat) and high spin Co(ii)-semiquinonate (HS-CoII-sq) forms. Two families of these compounds have been investigated for decades but have generally been considered separately: neutral [Co(diox)(sq)(N2L)] and cationic [Co(diox)(N4L)]+ complexes (diox = generic dioxolene, N2L/N4L = bidentate/tetradentate N-donor ancillary ligand). Computational identification of promising new candidate compounds prior to experimental exploration is beneficial for environmental and cost considerations but requires a thorough understanding of the underlying thermochemical parameters that influence the switching. Herein, we report a robust approach for the analysis of both cobalt-dioxolene families, which involved a quantitative density functional theory-based study benchmarked with reliable quasi-experimental references. The best-performing M06L-D4/def2-TZVPP level of theory has subsequently been verified by the synthesis and experimental investigation of three new complexes, two of which exhibit thermally-induced VT, while the third remains in the LS-CoIII-cat form across all temperatures, in agreement with prediction. Valence tautomerism in solution is markedly solvent-dependent, but the origin of this has not been definitively established. We have extended our computational approach to elucidate the correlation of VT transition temperature with solvent stabilisation energy and change in dipole moment. This new understanding may inform the development of VT compounds for applications in soft materials including films, gels, and polymers.

Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder

The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.

Role of Explicit Solvation in Computational Modeling of Chemical Reactions: Mechanism of Cu(I) Transfer Between Thiolate-Based Chelators in Water

Solvation plays important roles in controlling the thermodynamic and kinetic aspects of chemical reactions. The conventional approaches to treat solvation via electronic structure methods are likely to become inadequate, when the reacting solutes have strong electrostatic and hydrogen bonding interactions with the solvent and undergo significant structural changes during the course of the reaction. In this article, we present evidence of such solvent and structural effects in the computational study of the Cu(I) transfer reaction between thiolate-based chelators dithiobutylamine (DTBA) and dithiotheritol (DTT) in water, inspired from biological copper trafficking phenomena. We propose a general solution to the problem by combining classical molecular dynamics (MD) simulations of the bulk system and static quantum chemistry calculations. The fluctuating solvation shell was estimated from MD, and energetics was assessed by averaging QM energies of a series of molecular clusters constructed from the MD snapshots. Applying this approach, we propose a reaction pathway with estimates of relative intermediate stabilities and barriers, which suggest the overall reaction to be reversible in nature and likely to go through both two and three coordinated intermediates, confirming previous studies of similar protein analogues. An interesting fact that emerged from our study was the strong indication that the rate-determining step is the deprotonation of initial thiol bound Cu(I) complex, without involving any Cu(I)-S bonds. The proposed method will lead to a better treatment of solvations, and these mechanistic insights will aid our understanding of biological copper(I) trafficking.

Benchmark Study on the Calculation of 207Pb NMR Chemical Shifts

A benchmark set for the computation of 207Pb nuclear magnetic resonance (NMR) chemical shifts is presented. The PbS50 set includes conformer ensembles of 50 lead-containing molecular compounds and their experimentally measured 207Pb NMR chemical shifts. Various bonding motifs at the Pb center with up to seven bonding partners are included. Six different solvents were used in the measurements. The respective shifts lie in the range between +10745 and -5030 ppm. Several calculation settings are assessed by evaluating computed 207Pb NMR shifts for the use with different density functional approximations (DFAs), relativistic approaches, treatment of the conformational space, and levels for geometry optimization. Relativistic effects were included explicitly with the zeroth order regular approximation (ZORA), for which only the spin-orbit variant was able to yield reliable results. In total, seven GGAs and three hybrid DFAs were tested. Hybrid DFAs significantly outperform GGAs. The most accurate DFAs are mPW1PW with a mean absolute deviation (MAD) of 429 ppm and PBE0 with an MAD of 446 ppm. Conformational influences are small as most compounds are rigid, but more flexible structures still benefit from Boltzmann averaging. Including explicit relativistic treatments such as SO-ZORA in the geometry optimization does not show any significant improvement over the use of effective core potentials (ECPs).

Ab Initio and Kinetic Modeling of β-d-Xylopyranose under Fast Pyrolysis Conditions

Lignocellulosic biomass is an abundant renewable resource that can be upgraded to chemical and fuel products through a range of thermal conversion processes. Fast pyrolysis is a promising technology that uses high temperatures and fast heating rates to convert lignocellulose into bio-oils in high yields in the absence of oxygen. Hemicellulose is one of the three major components of lignocellulosic biomass and is a highly branched heteropolymer structure made of pentose, hexose sugars, and sugar acids. In this study, β-d-xylopyranose is proposed as a model structural motif for the essential chemical structure of hemicellulose. The gas-phase pyrolytic reactivity of β-d-xylopyranose is thoroughly investigated using computational strategies rooted in quantum chemistry. In particular, its thermal degradation potential energy surfaces are computed employing Minnesota global hybrid functional M06-2X in conjunction with the 6-311++G(d,p) Pople basis set. Electronic energies are further refined by performing DLPNO-CCSD(T)-F12 single-point calculations on top of M06-2X geometries using the cc-pVTZ-F12 basis set. Conformational analysis for minima and transition states is performed with state-of-the-art semiempirical quantum chemical methods coupled with metadynamics simulations. Key thermodynamic quantities (free energies, barrier heights, enthalpies of formation, and heat capacities) are computed. Rate coefficients for the initial steps of thermal decomposition are computed by means of reaction rate theory. For the first time, a detailed elementary reaction kinetic model for β-d-xylopyranose is developed by utilizing the thermodynamic and kinetic information acquired from the aforementioned calculations. This model specifically targets the initial stages of β-d-xylopyranose pyrolysis in the high-pressure limit, aiming to gain a deeper understanding of its reaction kinetics. This approach establishes a systematic strategy for exploring reactive pathways, evaluating competing parallel reactions, and selectively accepting or discarding pathways based on the analysis. The findings suggest that acyclic d-xylose plays a significant role as an intermediary in the production of key pyrolytic compounds during the pyrolysis of xylose. These compounds include furfural, anhydro-d-xylopyranose, glycolaldehyde, and dihydrofuran-3(2H)-one.

Do Optimally Tuned Range-Separated Hybrid Functionals Require a Reparametrization of the Dispersion Correction? It Depends

For ground- and excited-state studies of large molecules, it is the state of the art to combine (time-dependent) DFT with dispersion-corrected range-separated hybrid functionals (RSHs), which ensures an asymptotically correct description of exchange effects and London dispersion. Specifically for studying excited states, it is common practice to tune the range-separation parameter ω (optimal tuning), which can further improve the accuracy. However, since optimal tuning essentially changes the functional, it is unclear if and how much the parameters used for the dispersion correction depend on the chosen ω value. To answer this question, we explore this interdependency by refitting the DFT-D4 dispersion model for six established RSHs over a wide range of ω values (0.05-0.45 a0-1) using a set of noncovalently bound molecular complexes. The results reveal some surprising differences among the investigated functionals: While PBE-based RSHs and ωB97M-D4 generally exhibit a weak interdependency and robust performance over a wide range of ω values, B88-based RSHs, specifically LC-BLYP, are strongly affected. For these, even a minor reduction of ω from the default value manifests in strong systematic overbinding and poor performance in the typical range of optimally tuned ω values. Finally, we discuss strategies to mitigate these issues and reflect the results in the context of the employed D4 parameter optimization algorithm and fit set, outlining strategies for future improvements.

Investigating the accuracy of density functional methods for molecules in electric fields

The use of oriented external electric fields (OEEFs) as a potential tool for catalyzing chemical reactions has gained traction in recent years. Electronic structure calculations using OEEFs are commonly done using methods based on density functional theory (DFT), but until now, the performance of DFT methods for calculating molecules in OEEFs had not been assessed in a more general scope. Looking at the accuracy of molecular geometries, electronic energies, and electric dipole moments compared to accurate coupled-cluster with perturbative triples data, we have investigated a wide variety of density functionals using different basis sets to determine how well the individual functionals perform on various types of chemical bonds. We found that most functionals accurately calculate geometries in OEEFs and that small basis sets are sufficient in many cases. Calculations of electronic energies show a significant error introduced by the OEEF, which the use of a larger basis set helps mitigate. Our findings show that DFT methods can be used for accurate calculations in OEEFs, allowing researchers to make full use of the advantages that they bring.

Feed-Forward Neural Network for Predicting Enantioselectivity of the Asymmetric Negishi Reaction

Density functional theory (DFT) is a powerful tool to model transition state (TS) energies to predict selectivity in chemical synthesis. However, a successful multistep synthesis campaign must navigate energetically narrow differences in pathways that create some limits to rapid and unambiguous application of DFT to these problems. While powerful data science techniques may provide a complementary approach to overcome this problem, doing so with the relatively small data sets that are widespread in organic synthesis presents a significant challenge. Herein, we show that a small data set can be labeled with features from DFT TS calculations to train a feed-forward neural network for predicting enantioselectivity of a Negishi cross-coupling reaction with P-chiral hindered phosphines. This approach to modeling enantioselectivity is compared with conventional approaches, including exclusive use of DFT energies and data science approaches, using features from ligands or ground states with neural network architectures.

Extended Conductor-like Polarizable Continuum Solvation Model (CPCM-X) for Semiempirical Methods

We have developed a new method to accurately account for solvation effects in semiempirical quantum mechanics based on a polarizable continuum model (PCM). The extended conductor-like polarizable continuum model (CPCM-X) incorporates a computationally efficient domain decomposition conductor-like screening model (ddCOSMO) for extended tight binding (xTB) methods and uses a post-processing approach based on established solvation models, like the conductor-like screening model for real solvents (COSMO-RS) and the universal solvent model based on solute electron density (SMD). According to various benchmarks, the approach performs well across a broad range of systems and applications, including hydration free energies, non-aqueous solvation free energies, and large supramolecular association reactions of neutral and charged species. Our method for computing solvation free energies is much more accurate than the current methods in the xtb program package. It improves the accuracy of solvation free energies by up to 40% for larger supramolecular association reactions to match even the accuracy of higher-level DFT-based solvation models like COSMO-RS and SMD while being computationally more than 2 orders of magnitude faster. The proposed method and the underlying ddCOSMO model are readily available for a wide variety of solvents and are accessible in xtb for use in various computational applications.

Two Shared Icosahedral Metallacarboranes through Iron: A Joint Experimental and Theoretical Refinement of Mössbauer Spectrum in [Fe(1,2-C2B9H11)2]Cs

Mössbauer and X-ray photoelectron spectroscopies (XPS) are complemented with high-level quantum-chemical computations in the study of the geometric and electronic structure of the paramagnetic salt of the metallacarborane sandwich complex [Fe(1,2-C₂B₉H₁₁)₂]Cs = FeSanCs. Experimental ⁵⁷Fe isomer shifts and quadrupole splitting parameters are compared with the theoretical prediction, with good agreement. The appearance of two sets of Cs(3d) doublets in the XPS spectrum, separated by 2 eV, indicates that Cs has two different chemical environments due to ease of the Cs⁽⁺⁾ cation moving around the sandwich complex with low-energy barriers, as confirmed by quantum-chemical computations. Several minimum-energy geometries of the FeSanCs structure with the corresponding energies and Mössbauer parameters are discussed, in particular the atomic charges and spin population and the surroundings of the Fe atom in the complex. The Mössbauer spectra were taken at different temperatures showing the presence of a low-spin Fe atom with S = 1/2 and thus confirming a paramagnetic Fe^III species.

Three-Component Synthesis of Dioxaphosphorane-Fused Diphosphacycles Exhibiting Unique Dynamic Fluorescence "On/Off" Properties

Rigidly planar polycyclic phosphacycles featuring an internal dioxaphosphorane are promising photofunctional materials. However, the lack of efficient synthetic methods resulted in limited structural diversities which significantly hampered extensive study. Herein, we report a straightforward three-component synthesis of novel dioxaphosphorane-fused diphosphacycles with distinctive photophysical properties. Control experiments and theory calculations were performed to account for a plausible reaction mechanism. We also systematically investigated the structure-property relationships of these unprecedented platforms by combining experiments (X-ray analysis, optical and redox properties) and theoretical computations. Based on their unique structure and properties, a novel fluorescent switch for pH sensing was revealed by a dynamic ring-opening/ring-closing process.

Angewandte Chemie: A Team Effort

While the number of submissions to the journal continued to increase in 2022, the Editorial Team also implemented a number of changes with support of the GDCh and the journal's advisory bodies. For example, "Team Profiles" are author profiles that celebrate joint efforts as the basis of an inclusive research culture. In looking forward, two new members will join the International Advisory Board in 2023.

A Guide to In Silico Drug Design

The drug discovery process is a rocky path that is full of challenges, with the result that very few candidates progress from hit compound to a commercially available product, often due to factors, such as poor binding affinity, off-target effects, or physicochemical properties, such as solubility or stability. This process is further complicated by high research and development costs and time requirements. It is thus important to optimise every step of the process in order to maximise the chances of success. As a result of the recent advancements in computer power and technology, computer-aided drug design (CADD) has become an integral part of modern drug discovery to guide and accelerate the process. In this review, we present an overview of the important CADD methods and applications, such as in silico structure prediction, refinement, modelling and target validation, that are commonly used in this area.

An Efficient Approach to Large-Scale Ab Initio Conformational Energy Profiles of Small Molecules

Accurate conformational energetics of molecules are of great significance to understand maby chemical properties. They are also fundamental for high-quality parameterization of force fields. Traditionally, accurate conformational profiles are obtained with density functional theory (DFT) methods. However, obtaining a reliable energy profile can be time-consuming when the molecular sizes are relatively large or when there are many molecules of interest. Furthermore, incorporation of data-driven deep learning methods into force field development has great requirements for high-quality geometry and energy data. To this end, we compared several possible alternatives to the traditional DFT methods for conformational scans, including the semi-empirical method GFN2-xTB and the neural network potential ANI-2x. It was found that a sequential protocol of geometry optimization with the semi-empirical method and single-point energy calculation with high-level DFT methods can provide satisfactory conformational energy profiles hundreds of times faster in terms of optimization.

Multiple Facets of Modeling Electronic Absorption Spectra of Systems in Solution

In this Perspective, we outline the essential physicochemical aspects that need to be considered when building a reliable approach to describe absorption properties of solvated systems. In particular, we focus on how to properly model the complexity of the solvation phenomenon, arising from dynamical aspects and specific, strong solute-solvent interactions. To this end, conformational and configurational sampling techniques, such as Molecular Dynamics, have to be coupled to accurate fully atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) methodologies. By exploiting different illustrative applications, we show that an effective reproduction of experimental spectral signals can be achieved by delicately balancing exhaustive sampling, hydrogen bonding, mutual polarization, and nonelectrostatic effects.

"Vermellogens" and the Development of CB[8]-Based Supramolecular Switches Using pH-Responsive and Non-Toxic Viologen Analogues

We present herein the "vermellogens", a new class of pH-responsive viologen analogues, which replace the direct linking between para-substituted pyridinium moieties within those by a hydrazone functional group. A series of such compounds have been efficiently synthesized in aqueous media by hydrazone exchange reactions, displaying a marked pH-responsivity. Furthermore, the parent N,N'-dimethylated "vermellogen": the "red thread", an analogue of the herbicide paraquat and used herein as a representative model of the series, showed anion-recognition abilities, non-reversible electrochemical behavior, and non-toxicity of the modified bis-pyridinium core. The host-guest chemistry for the "red thread" with the CB[7,8] macrocyclic receptors has been extensively studied experimentally and by dispersion corrected density functional theory methods, showing a parallel behavior to that previously described for the herbicide but, crucially, swapping the well-known redox reactive capabilities of the viologen-based inclusion complexes by acid-base supramolecular responsiveness.