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

Lvsiyujins A–G, new sesquiterpenoids, from Curcuma phaeocaulis Valeton root tuber and their preliminary pharmacological property assessment based on ADME evaluation, molecular docking and in vitro experiments

Seven new sesquiterpenoids were isolated from the root tuber of C. phaeocaulis. A combination of calculations and experiments was used in structural analysis and biological activity exploration.

Reliable prediction of association (free) energies of supramolecular complexes with heavy main group elements – the HS13L benchmark set

We introduce a set of 13 supramolecular complexes featuring diverse non-covalent interactions with heavy main group elements (Zn, As, Se, Te, Br, I), high charges (−2 up to +4), and large systems with up to 266 atoms (HS13L).

Along the road to Crystal Structure Prediction (CSP) of pharmaceutical-like molecules

Computational methods used for predicting crystal structures of organic compounds are mature enough to be routinely used with many rigid and semi-rigid organic molecules. The usefulness of Crystal Structure Prediction...

Quickstart guide to model structures and interactions of artificial molecular muscles with efficient computational methods

Artificial molecular muscles (AMMs) represent an important group of molecular machines. Their theoretical treatment is challenging due to size, element composition, and complex interaction motifs. Moreover, experimentally determined structures often only yield insights into the covalent connectivity of atoms rather than their 3D structure. Accordingly, a reproducible computational modeling of such structures is complicated. In this work we present a standardized, mostly quantum chemical protocol on how to obtain reliable structures from scratch and to compute contraction free energies ΔG_c for daisy-chain rotaxane AMMs efficiently. In this protocol, the recently developed force-field (GFN-FF) and extended tight-binding methods (GFNn-xTB) are employed. For comparison, dispersion-corrected density functional theory (DFT-D) based reference ΔG_c were computed. In one case for which data are available, excellent agreement between theoretical and experimental ΔG_c values within 1-2 kcal mol^-1 is obtained.

2-Amino-4,5-dihydrothiophene-3-carbonitriles: A New Synthesis, Quantum Chemical Studies, and Mannich-Type Reactions Leading to New Hexahydrothieno[2,3-d]pyrimidines

trans-2-Amino-4-aryl-5-benzoyl-4,5-dihydrothiophene-3-carbonitriles were prepared either by the reaction of 3-aryl-2-cyanothioacrylamides with α-thiocyanatoacetophenone or by the Michael-type addition of cyanothioacetamide to α-bromochalcones followed by intramolecular cyclization. The mechanism of the first reaction was studied using high-level quantum chemical calculations. Density functional theory (DFT) studies were carried out to determine the mechanism of the first reaction. A new approach toward the construction of the thieno[2,3-d]pyrimidine core system was demonstrated by the reaction of the prepared dihydrothiophenes with HCHO and RNH2 under noncatalyzed Mannich conditions.

Ligand Protonation at Carbon, not Nitrogen, during H2 Production with Amine-Rich Iron Electrocatalysts

We present monometallic H₂ production electrocatalysts containing electron-rich triamine-cyclopentadienyl (Cp) ligands coordinated to iron. After selective CO extrusion from the iron tricarbonyl precursors, electrocatalysis is observed via cyclic voltammetry in the presence of an exogenous acid. Contrary to the fact that amines in the secondary coordination sphere are often protonated during electrocatalysis, comprehensive quantum-chemical calculations indicate that the amines likely do not function as proton relays; instead, endo-Cp ring protonation is most favorable after 1e^- reduction. This unusual mechanistic pathway emphasizes the need to consider a broad domain of H⁺/e^- addition products by synergistically combining experimental and theoretical resources.

Zn2+ detection of a benzimidazole 8-aminoquinoline fluorescent sensor by inhibited tautomerization

A new fluorescent chemosensor based on 8-aminoquinoline L1 bearing a benzimidazole moiety was synthesized, which exists as two predominant tautomers L1A and L1B in diluted DMSO-d6 solution. Among various metal ions, L1 showed a highly selective and sensitive turn-on fluorescence response to the presence of Zn2+ ions in methanol. The detection limit for Zn2+ by L1 was calculated to be 1.76 × 10-7 M. The 1 : 1 complexation ratio of the L1-Zn complex was confirmed through Job plot measurements. Complexation studies were performed by FT-IR, NMR and HR-ESI MS measurements and DFT calculations. With the gained insight, it was possible to successfully apply L1 in water sample analysis.

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

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

Critical Role of Thermal Fluctuations for CO Binding on Electrocatalytic Metal Surfaces

This work considers the evaluation of density functional theory (DFT) when comparing against experimental observations of CO binding trends on the strong binding Pt(111) and intermediate binding Cu(111) and for weak binding Ag(111) and Au(111) surfaces important in electrocatalysis. By introducing thermal fluctuations using appropriate statistical mechanical NVT and NPT ensembles, we find that the RPBE and B97M-rV DFT functionals yield qualitatively better metal surface strain trends and CO enthalpies of binding for Cu(111) and Pt(111) than found at 0 K, thereby correcting the overbinding by 0.2 to 0.3 eV to yield better agreement with the enthalpies determined from experiment. The importance of dispersion effects are manifest for the weak CO binding Ag(111) and Au(111) surfaces at finite temperatures in which the RPBE functional does not bind CO at all, while the B97M-rV functional shows that the CO-metal interactions are a mixture of chemisorbed and physisorbed species with binding enthalpies that are within ∼0.05 eV of experiment. Across all M(111) surfaces, we show that the B97M-rV functional consistently predicts the correct atop site preference for all metals due to thermally induced surface distortions that preferentially favor the undercoordinated site. This study demonstrates the need to fully account for finite temperature fluctuations to make contact with the binding enthalpies from surface science experiments and electrocatalysis applications.

Assessing Density Functional Theory for Chemically Relevant Open-Shell Transition Metal Reactions

Due to the principle lack of systematic improvement possibilities of density functional theory, careful assessment of the performance of density functional approximations (DFAs) on well-designed benchmark sets, for example, for reaction energies and barrier heights, is crucial. While main-group chemistry is well covered by several available sets, benchmark data for transition metal chemistry is sparse. This is especially the case for larger, chemically relevant molecules. Addressing this issue, we recently introduced the MOR41 benchmark which covers chemically relevant reactions of closed-shell complexes. In this work, we extend these efforts to single-reference open-shell systems and introduce the "reactions of open-shell single-reference transition metal complexes" (ROST61) benchmark set. ROST61 includes accurate coupled-cluster reference values for 61 reaction energies with a mean reaction energy of -42.8 kcal mol^-1. Complexes with 13-93 atoms covering 20 d-block elements are included, but due to the restriction to single-reference open-shell systems, important elements such as iron or platinum could not be taken into account, or only to a small extent. We assess the performance of 31 DFAs in combination with three London dispersion (LD) correction schemes. Further, DFT-based composite methods, MP2, and a few semiempirical quantum chemical methods are evaluated. Consistent with the results for the MOR41 closed-shell benchmark, we find that the ordering of DFAs according to Jacob's ladder is preserved and that adding an LD correction is crucial, clearly improving almost all tested methods. The recently introduced r²SCAN-3c composite method stands out with a remarkable mean absolute deviation (MAD) of only 2.9 kcal mol^-1, which is surpassed only by hybrid DFAs with low amounts of Fock exchange (e.g., 2.3 kcal mol^-1 for TPSS0-D4/def2-QZVPP) and double-hybrid (DH) DFAs but at a significantly higher computational cost. The lowest MAD of only 1.6 kcal mol^-1 is obtained with the DH DFA PWPB95-D4 in the def2-QZVPP basis set approaching the estimated accuracy of the reference method. Overall, the ROST61 set adds important reference data to a sparsely sampled but practically relevant area of chemistry. At this point, it provides valuable orientation for the application and development of new DFAs and electronic structure methods in general.

Intramolecular amination via acid-catalyzed rearrangement of azides: a potent alternative to intermolecular direct electrophilic route

Although acid-catalyzed intramolecular rearrangement of organic azides is an attractive route to amines, its mechanism and synthetic prospective are still debated. Herein, through computational and experimental studies, we demonstrated that azide intramolecular rearrangement could serve as a potent synthetic route to a sought-after amine functionality including preparation of difficult to access and valuable heterocyclic amines. Using quantum chemical calculations at MP2/aug-cc-pVTZ and B3LYP/aug-cc-pVDZ levels, we discovered that this reaction proceeds via a concerted transition state with nitrogen elimination and alkyl/aryl migration occurring at the same time. Two conformers of protonated azides - syn- and anti- - were shown to precede corresponding transition states. It was shown that the reaction follows Curtin-Hammett scenario as the energy gap required for conformer interconversion was substantially lower than activation barrier of either transition state. Intramolecular amination via azide rearrangement was predicted to be a selective process with migratory aptitude increasing in a row alkyl<Ar-EWG<Ar-EDG (EWG - electron withdrawing group; EDG - electron donating group), which was supported by experimental results. We demonstrated experimentally that organic azides can be generated from stilbenes in situ and selectively undergo rearrangement to corresponding amines in a cascade fashion via amino-dealkenylation reaction.

Benchmark test of a dispersion corrected revised Tao-Mo semilocal functional for thermochemistry, kinetics, and noncovalent interactions of molecules and solids

In the density functional theory, dispersion corrected semilocal approximations are often used to benchmark weekly interacting finite and extended systems. Here, the focus is on providing a broad overview of the performance of D3 dispersion corrected revised Tao-Mo (revTM) semilocal functionals [A. Patra et al., J. Chem. Phys. 153, 084 117 (2020)] for thermochemistry and kinetics of molecules, molecular crystals, ice polymorphs, metal-organic systems, atom/molecular adsorption on solids, water interacting with nano-materials, binding energies of layered materials, and properties of weekly and strongly bonded solids. We show that the most suitable "optimized power" function for the revTM functional needs a modification to make it suitable for properties related to the diverse nature of finite and extended systems. The present work is an extension of the previously proposed revTM+D3 method with the motivation to design and benchmark the dispersion corrected cost-effective method based on this semilocal approximation. We show that the revised revTM+D3 functional provides various general purpose molecular and solid properties with the closest to experimental findings than its predecessor. The present assessment and benchmarking can be practically useful for performing cost-effective method based simulations of various molecular and solid-state properties.

Revisiting conformations of methyl lactate in water and methanol

The recently developed efficient protocols to implicit [Grimme et al., J. Phys. Chem. A 125, 4039-4054 (2021)] and explicit quantum mechanical modeling of non-rigid molecules in solution [Katsyuba et al., J. Phys. Chem. B 124, 6664-6670 (2020)] are applied to methyl lactate (ML). Building upon this work, a new combination scheme is proposed to incorporate solvation effects for the computation of infrared (IR) absorption spectra. Herein, Boltzmann populations calculated for implicitly solvated single conformers are used to weight the IR spectra of explicitly solvated clusters with a size of typically ten solvent molecules, i.e., accounting for the first solvation shell. It is found that in water and methanol, the most abundant conformers of ML are structurally modified relative to the gas phase, where the major form is ML1, in which the syn conformation of the -OH moiety is stabilized by a OH⋯O=C intramolecular hydrogen bond (HB). In solution, this syn conformation transforms to the gauche form because the intramolecular HB is disrupted by explicit water molecules that form intermolecular HBs with the hydroxyl and carbonyl groups. Similar changes induced by the gas-solution transition are observed for the minor conformers, ML2 and/or ML3, characterized by OH⋯OCH₃ intramolecular HB in the gas phase. The relative abundance of ML1 is shown to decrease from ∼96% in gas to ∼51% in water and ∼92% in methanol. The solvent strongly influences frequencies, IR intensities, and normal modes, resulting in qualitatively different spectra compared to the gas phase. Some liquid-state conformational markers in the fingerprint region of IR spectra are revealed.

Robust and Efficient Implicit Solvation Model for Fast Semiempirical Methods

We present a robust and efficient method to implicitly account for solvation effects in modern semiempirical quantum mechanics and force fields. A computationally efficient yet accurate solvation model based on the analytical linearized Poisson-Boltzmann (ALPB) model is parameterized for the extended tight binding (xTB) and density functional tight binding (DFTB) methods as well as for the recently proposed GFN-FF general force field. The proposed methods perform well over a broad range of systems and applications, from conformational energies over transition-metal complexes to large supramolecular association reactions of charged species. For hydration free energies of small molecules, GFN1-xTB(ALPB) is reaching the accuracy of sophisticated explicitly solvated approaches, with a mean absolute deviation of only 1.4 kcal/mol compared to the experiment. Logarithmic octanol-water partition coefficients (log K_ow) are computed with a mean absolute deviation of about 0.65 using GFN2-xTB(ALPB) compared to experimental values indicating a consistent description of differential solvent effects. Overall, more than twenty solvents for each of the six semiempirical methods are parameterized and tested. They are readily available in the xtb and dftb+ programs for diverse computational applications.

Titanocene-Catalyzed [2+2] Cycloaddition of Bisenones and Comparison with Photoredox Catalysis and Established Methods

Cp2 Ti(TFA) is a broadly applicable catalyst for the [2+2] cycloaddition of bisenones by inner-sphere electron transfer. The attractiveness of this mechanism is shown by comparison with outer-sphere ET methods. DFT calculations show that the reaction proceeds through a unique unfavorable 5-exo (the rate-determining step) and a favorable 4-exo cyclization.

Efficient Quantum-Chemical Calculations of Acid Dissociation Constants from Free-Energy Relationships

The calculation of acid dissociation constants (pK_a) is an important task in computational chemistry and chemoinformatics. Theoretically and with minimal empiricism, this is possible from computed acid dissociation free energies via so-called linear free-energy relationships. In this study some modifications are introduced to the latter, providing a straightforward, broadly applicable protocol with an adjustable degree of sophistication for quantum chemistry-based calculations of pK_a in water. It targets a wide pK_a range (∼70 units) and medium-sized, flexible molecules. Herein, a focus is set on the recently published r²SCAN-3c and related efficient composite density functionals and the semiempirical GFN2-xTB method, including a newly introduced energy correction for heterolytic dissociation, both in combination with implicit solvation models. The performance is evaluated in comparison with experimental data, showing mean errors often smaller than a targeted 1 pK_a unit accuracy. Larger deviations are observed only upon inclusion of challenging highly negative (<-5) or positive (>15) pK_a values. Among all those tested, it is found that B97-3c is the best performing functional, although rather independently of the density functional theory (DFT) method used; low root-mean-square errors of 0.8-1.0 pK_a units for typical drugs are obtained. For optimal performance, it is recommended to employ DFT functional specific free-energy relationship parameters. Additionally, a significant conformational dependence of the pK_a values is revealed and quantified for some nonrigid drug molecules.

Comment on "The Nature of Chalcogen-Bonding-Type Tellurium-Nitrogen Interactions": Fixing the Description of Finite-Temperature Effects Restores the Agreement Between Experiment and Theory

Mitzel and co-workers recently presented an intriguing molecule displaying a tellurium-nitrogen interaction. Structural data obtained in the solid and in gas phase indicated a large increase of the Te-N equilibrium distance re from 2.64 to 2.92 Å, respectively. Although some DFT calculations appear to support the large re in gas phase, we argue that the lions share of the increase is due to an incomplete description of finite-temperature effects in the back-corrected experimental data. This hypothesis is based on high-level coupled-cluster (CC) and periodic DFT calculations, which consistently point towards a much smaller re in the isolated molecule. Further support comes through MD simulations with a tuned GFN2-xTB Hamiltonian: Calibrated against a CC reference, these show a six-times larger influence of temperature than with the originally used GFN1-xTB. Taking this into account, the back-corrected re in gas phase becomes 2.67±0.08 Å, in good agreement with high-level CC theory and most DFT methods.

Benchmarking London dispersion corrected density functional theory for noncovalent ion-π interactions

The strongly attractive noncovalent interactions of charged atoms or molecules with π-systems are important binding motifs in many chemical and biological systems. These so-called ion-π interactions play a major role in enzymes, molecular recognition, and for the structure of proteins. In this work, a molecular test set termed IONPI19 is compiled for inter- and intramolecular ion-π interactions, which is well balanced between anionic and cationic systems. The IONPI19 set includes interaction energies of significantly larger molecules (up to 133 atoms) than in other ion-π test sets and covers a broad range of binding motifs. Accurate (local) coupled cluster values are provided as reference. Overall, 19 density functional approximations, including seven (meta-)GGAs, eight hybrid functionals, and four double-hybrid functionals combined with three different London dispersion corrections, are benchmarked for interaction energies. DFT results are further compared to wave function based methods such as MP2 and dispersion corrected Hartree-Fock. Also, the performance of semiempirical QM methods such as the GFNn-xTB and PMx family of methods is tested. It is shown that dispersion-uncorrected DFT underestimates ion-π interactions significantly, even though electrostatic interactions dominate the overall binding. Accordingly, the new charge dependent D4 dispersion model is found to be consistently better than the standard D3 correction. Furthermore, the functional performance trend along Jacob's ladder is generally obeyed and the reduction of the self-interaction error leads to an improvement of (double) hybrid functionals over (meta-)GGAs, even though the effect of the SIE is smaller than expected. Overall, the double-hybrids PWPB95-D4/QZ and revDSD-PBEP86-D4/QZ turned out to be the most reliable among all assessed methods for the description of ion-π interactions, which opens up new perspectives for systems where coupled cluster calculations are no longer computationally feasible.

Efficient Quantum Chemical Calculation of Structure Ensembles and Free Energies for Nonrigid Molecules

The application of quantum chemical, automatic multilevel modeling workflows for the determination of thermodynamic (e.g., conformation equilibria, partition coefficients, pK_a values) and spectroscopic properties of relatively large, nonrigid molecules in solution is described. Key points are the computation of rather complete structure (conformer) ensembles with extremely fast but still reasonable GFN2-xTB or GFN-FF semiempirical methods in the CREST searching approach and subsequent refinement at a recently developed, accurate r²SCAN-3c DFT composite level. Solvation effects are included in all steps by accurate continuum solvation models (ALPB, (D)COSMO-RS). Consistent inclusion of thermostatistical contributions in the framework of the modified rigid-rotor-harmonic-oscillator approximation (mRRHO) based on xTB/FF computed PES is also recommended.

Assessing the Applicability of the Geometric Counterpoise Correction in B2PLYP/Double-ζ Calculations for Thermochemistry, Kinetics, and Noncovalent Interactions

We present a proof-of-concept study of the suitability of Kruse and Grimme’s geometric counterpoise correction (gCP) for basis set superposition errors (BSSEs) in double-hybrid density functional calculations with a double-ζ basis set. The gCP approach only requires geometrical information as an input and no orbital/density information is needed. Therefore, this correction is practically free of any additional cost. gCP is trained against the Boys and Bernardi counterpoise correction across a set of 528 noncovalently bound dimers. We investigate the suitability of the approach for the B2PLYP/def2-SVP level of theory, and reveal error compensation effects—missing London dispersion and the BSSE—associated with B2PLYP/def2-SVP calculations, and present B2PLYP-gCP-D3(BJ)/def2-SVP with the reparametrised DFT-D3(BJ) and gCP corrections as a more balanced alternative. Benchmarking results on the S66x8 benchmark set for noncovalent interactions and the GMTKN55 database for main-group thermochemistry, kinetics, and noncovalent interactions show a statistical improvement of the B2PLYP-gCP-D3(BJ) scheme over plain B2PLYP and B2PLYP-D3(BJ). B2PLYP-D3(BJ) shows significant overestimation of interaction energies, barrier heights with larger deviations from the reference values, and wrong relative stabilities in conformers, all of which can be associated with BSSE. We find that the gCP-corrected method represents a significant improvement over B2PLYP-D3(BJ), particularly for intramolecular noncovalent interactions. These findings encourage future developments of efficient double-hybrid DFT strategies that can be applied when double-hybrid calculations with large basis sets are not feasible due to system size.