B97-3c: A revised low-cost variant of the B97-D density functional method

Fluoride Additive as a Simple Tool to Qualitatively Improve Performance of Nickel-Catalyzed Asymmetric Michael Addition of Malonates to Nitroolefins

Nowadays, design of the new chiral ligands for organometallic catalysts is often based on the step-by-step increase in their complexity to improve efficiency. Herein we describe that simple in situ addition of the fluoride source to the asymmetric organometallic catalyst can improve not only activity but also enantioselectivity. Bromide-nickel diimine complexes were found to catalyze asymmetric Michael addition in low yields and ee, but activation with fluoride leads to a significant improvement in catalyst performance. The developed approach was applied to prepare several enantioenriched GABA analogues.

Molecular recognition between membrane epitopes and nearly free surface silanols explains silica membranolytic activity

Inhaled crystalline silica causes inflammatory lung diseases, but the mechanism for its unique activity compared to other oxides remains unclear, preventing the development of potential therapeutics. Here, the molecular recognition mechanism between membrane epitopes and "nearly free silanols" (NFS), a specific subgroup of surface silanols, is identified and proposed as a novel broad explanation for particle toxicity in general. Silica samples having different bulk and surface properties, specifically different amounts of NFS, are tested with a set of membrane systems of decreasing molecular complexity and different charge. The results demonstrate that NFS content is the primary determinant of membrane disruption causing red blood cell lysis and changes in lipid order in zwitterionic, but not in negatively charged liposomes. NFS-rich silica strongly and irreversibly adsorbs zwitterionic self-assembled phospholipid structures. This selective interaction is corroborated by density functional theory and supports the hypothesis that NFS recognize membrane epitopes that exhibit a positive quaternary amino and negative phosphate group. These new findings define a new paradigm for deciphering particle-biomembrane interactions that will support safer design of materials and what types of treatments might interrupt particle-biomembrane interactions.

On the separation of enantiomers by drift tube ion mobility spectrometry

Racemic mixtures of twelve common α-amino acids and three chiral drugs were tested for the separation of their enantiomers by drift tube ion mobility spectrometry (IMS)-quadrupole mass spectrometry (QMS) by introducing chiral selectors into the buffer gas of the IMS instrument. (R)-α-(Trifluoromethyl)benzyl alcohol, (L)-ethyl lactate, methyl (S)-2-chloropropionate, and the R and S enantiomers of 2-butanol and 1-phenyl ethanol were evaluated as chiral selectors. Experimental conditions were varied during the tests, including buffer gas temperature, concentration and type of chiral selectors, analyte concentration, electrospray (ESI) voltage, ESI solvent pH, and buffer gas flow rate. The individual enantiomers yielded the same drift times and the racemic mixtures could not be separated. Energy calculations of the chiral selector-ion interactions showed that these separations are unlikely using 2-butanol as a chiral selector but they might be theoretically feasible depending on the chiral selector nature and the type of enantiomers. Several plausible explanations for not achieving separations were analyzed. Recommendations for potential enantiomer separations by IMS are proposed.

High-Level Quantum Chemistry Reference Heats of Formation for a Large Set of C, H, N, and O Species in the NIST Chemistry Webbook and the Identification and Validation of Reliable Protocols for Their Rapid Computation

A recent study has examined the accuracy of NIST heats of formation for a set of C, H, and O-containing species with a proposed low-cost quantum chemistry approach. In the present study, we have used high-level methods such as W1X-2 to obtain these data more rigorously, which we have then used to assess the NIST and the previously computed values. We find that many of these NIST data that are as suggested to be unreliable by the previous study are indeed inconsistent with our high-level reference values. However, we also find substantial deviations for the previously computed values from our benchmark. Thus, we have assessed the performance of alternative low-cost methods. In our assessment, we have additionally examined C, H, N, and O-containing species for which heats of formation are available from the NIST database. We find the ωB97M-V/ma-def2-TZVP, DSD-PBEP86/ma-def2-TZVP, and CCSD(T)-F12b/aug'-cc-pVDZ methods to be adequate for obtaining heats of formation with the atomization approach, once their atomic energies are optimized with our benchmark. Notably, the low-cost ωB97M-V method yields values that agree to be within 10 kJ mol^-1 for more than 90% of the (∼1500) species. A higher 20 kJ mol^-1 threshold captures 98% of the data. The outlier species typically contain many electron-withdrawing (nitro) groups. In these cases, the use of isodesmic-type reactions rather than the atomization approach is more reliable. Our assessment has also identified significant outliers from the NIST database, for which experimental re-determination of the heats of formation would be desirable.

Informing geometric deep learning with electronic interactions to accelerate quantum chemistry

Predicting electronic energies, densities, and related chemical properties can facilitate the discovery of novel catalysts, medicines, and battery materials. However, existing machine learning techniques are challenged by the scarcity of training data when exploring unknown chemical spaces. We overcome this barrier by systematically incorporating knowledge of molecular electronic structure into deep learning. By developing a physics-inspired equivariant neural network, we introduce a method to learn molecular representations based on the electronic interactions among atomic orbitals. Our method, OrbNet-Equi, leverages efficient tight-binding simulations and learned mappings to recover high-fidelity physical quantities. OrbNet-Equi accurately models a wide spectrum of target properties while being several orders of magnitude faster than density functional theory. Despite only using training samples collected from readily available small-molecule libraries, OrbNet-Equi outperforms traditional semiempirical and machine learning-based methods on comprehensive downstream benchmarks that encompass diverse main-group chemical processes. Our method also describes interactions in challenging charge-transfer complexes and open-shell systems. We anticipate that the strategy presented here will help to expand opportunities for studies in chemistry and materials science, where the acquisition of experimental or reference training data is costly.

Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C60 Dimer, DNA-Ellipticine, and HIV-Indinavir

In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C₆₀ dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C₆₀. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r²SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.

Fragmentation-Based Decomposition of a Metalloenzyme-Substrate Interaction: A Case Study for a Lytic Polysaccharide Monooxygenase

We present a novel decomposition scheme for electronic interaction energies based on the flexible formulation of fragmentation schemes through fragment combination ranges (FCRs; J. Chem. Phys., 2021, 155, 164105). We devise a clear additive decomposition with contribution of nondisjoint fragments and correction terms for overlapping fragments and apply this scheme to the metalloenzyme-substrate complex of a lytic polysaccharide monooxygenase (LPMO) with an oligosaccharide. By this, we further illustrate the straightforward adaptability of the FCR-based schemes to novel systems. Our calculations suggest that the description of the electronic structure is a larger error source than the fragmentation scheme. In particular, we find a large impact of the basis set size on the interaction energies. Still, the introduction of three-body interaction terms in the fragmentation setup improves the agreement to the supermolecular reference. Yet, the qualitative results for the decomposition scheme with two-body terms only largely agree within the investigated electronic-structure approaches and basis sets, which are B97-3c, DFT (TPSS and B3LYP), and MP2 methods. The overlap contributions are found to be small, allowing analysis of the interaction energy into individual amino acid residues: We find a particularly strong interaction between the substrate and the LPMO copper active site.

Manganese(I) Complex with Monodentate Arylisocyanide Ligands Shows Photodissociation Instead of Luminescence

Recently reported manganese(I) complexes with chelating arylisocyanide ligands exhibit luminescent metal-to-ligand charge-transfer (MLCT) excited states, similar to ruthenium(II) polypyridine complexes with the same d6 valence electron configuration used for many different applications in photophysics and photochemistry. However, chelating arylisocyanide ligands require substantial synthetic effort, and therefore it seemed attractive to explore the possibility of using more readily accessible monodentate arylisocyanides instead. Here, we synthesized the new Mn(I) complex [Mn(CNdippPhOMe2)6]PF6 with the known ligand CNdippPhOMe2 = 4-(3,5-dimethoxyphenyl)-2,6-diisopropylphenylisocyanide. This complex was investigated by NMR spectroscopy, single-crystal structure analysis, high-resolution electrospray ionization mass spectrometry (HR-ESI-MS) measurements, IR spectroscopy supported by density functional theory (DFT) calculations, cyclic voltammetry, and time-resolved as well as steady-state UV-vis absorption spectroscopy. The key finding is that the new Mn(I) complex is nonluminescent and instead undergoes arylisocyanide ligand loss during continuous visible laser irradiation into ligand-centered and charge-transfer absorption bands, presumably owed to the population of dissociative d-d excited states. Thus, it seems that chelating bi- or tridentate binding motifs are essential for obtaining emissive MLCT excited states in manganese(I) arylisocyanides. Our work contributes to understanding the basic properties of photoactive first-row transition metal complexes and could help advance the search for alternatives to precious metal-based luminophores, photocatalysts, and sensors.

Calculation of Lipophilicity of Organophosphate Pesticides Using Density Functional Theory

Higher lipophilicity facilitates the passage of a substance across lipid cell membranes, the blood–brain barrier and protein binding, and may also indicate its toxicity. We proposed eight methods for predicting the lipophilicity of the 22 most commonly used organophosphate pesticides. In this work, to determine the lipophilicity and thermodynamic parameters of the solvation of pesticides, we used methods of density functional theory with various basis sets, as well as modern Grimm methods. The prediction models were evaluated and compared against eight performance statistics, as well as time and RAM used in the calculation. The results show that the PBE-SVP method provided the best of the proposed predictive capabilities. In addition, this method consumes relatively less CPU and RAM resources. These methods make it possible to reliably predict the ability of pesticide molecules to penetrate cell membranes and have a negative effect on cells and the organism as a whole.

Optimization of the r2SCAN-3c Composite Electronic-Structure Method for Use with Slater-Type Orbital Basis Sets

The "Swiss army knife" composite density functional electronic-structure method r²SCAN-3c (J. Chem. Phys. 2021, 154, 064103) is extended and optimized for the use with Slater-type orbital basis sets. The meta generalized-gradient approximation (meta-GGA) functional r²SCAN by Furness et al. is combined with a tailor-made polarized triple-ζ Slater-type atomic orbital (STO) basis set (mTZ2P), the semiclassical London dispersion correction (D4), and a geometrical counterpoise (gCP) correction. Relativistic effects are treated explicitly with the scalar-relativistic zeroth-order regular approximation (SR-ZORA). The performance of the new implementation is assessed on eight geometry and 74 energy benchmark sets, including the extensive GMTKN55 database as well as recent sets such as ROST61 and IONPI19. In geometry optimizations, the STO-based r²SCAN-3c is either on par with or more accurate than the hybrid density functional approximation M06-2X-D3(0)/TZP. In energy calculations, the overall accuracy is similar to the original implementation of r²SCAN-3c with Gaussian-type atomic orbitals (GTO), but basic properties, intermolecular noncovalent interactions, and barrier heights are better described with the STO approach, resulting in a lower weighted mean absolute deviation (WTMAD-2(STO) = 7.15 vs 7.50 kcal mol^-1 with the original method) for the GMTKN55 database. The STO-optimized r²SCAN-3c outperforms many conventional hybrid/QZ approaches in most common applications at a fraction of their cost. The reliable, robust, and accurate r²SCAN-3c implementation with STOs is a promising alternative to the original implementation with GTOs and can be generally used for a broad field of quantum chemical problems.

Conformational Energy Benchmark for Longer n-Alkane Chains

We present the first benchmark set focusing on the conformational energies of highly flexible, long n-alkane chains, termed ACONFL. Unbranched alkanes are ubiquitous building blocks in nature, so the goal is to be able to calculate their properties most accurately to improve the modeling of, e.g., complex (biological) systems. Very accurate DLPNO-CCSD(T1)/CBS reference values are provided, which allow for a statistical meaningful evaluation of even the best available density functional methods. The performance of established and modern (dispersion corrected) density functionals is comprehensively assessed. The recently introduced r²SCAN-V functional shows excellent performance, similar to efficient composite DFT methods like B97-3c and r²SCAN-3c, which provide an even better cost-accuracy ratio, while almost reaching the accuracy of much more computationally demanding hybrid or double hybrid functionals with large QZ AO basis sets. In addition, we investigated the performance of common wave function methods, where MP2/CBS surprisingly performs worse compared to the simple D4 dispersion corrected Hartree-Fock. Furthermore, we investigate the performance of several semiempirical and force field methods, which are commonly used for the generation of conformational ensembles in multilevel workflows or in large scale molecular dynamics studies. Outstanding performance is obtained by the recently introduced general force field, GFN-FF, while other commonly applied methods like the universal force field yield large errors. We recommend the ACONFL as a helpful benchmark set for parametrization of new semiempirical or force field methods and machine learning potentials as well as a meaningful validation set for newly developed DFT or dispersion methods.

Sniffing out camphor: the fine balance between hydrogen bonding and London dispersion in the chirality recognition with α-fenchol

Binary complexes between the chiral monoterpenoids camphor and α-fenchol were explored with vibrational and rotational jet spectroscopy as well as density functional theory in order to explore how chirality can influence the binding preferences in gas-phase complexes. The global minimum structures of the two diastereomers were assigned. It is found that chirality recognition leads to different compromises in the fine balance between intermolecular interactions. While one isomer features a stronger hydrogen bond, the other one is more tightly arranged and stabilized by larger London dispersion interactions. These new spectroscopic results help understand the influence of chirality in molecular aggregation and unveil the kind of interactions involved between a chiral alcohol and a chiral ketone with large dispersion contributions.

Jatrolignans C and D: New Neolignan Epimers from Jatropha curcas

Two new neolignans jatrolignans, C (1) and D (2), a pair of epimers, were isolated from the whole plants of Jatropha curcas L. (Euphorbiaceae). Their structures were determined with HRESIMS, IR, and NMR data analysis, and electronic circular dichroism (ECD) experiments via a comparison of the experimental and the calculated ECD spectra. Their antichlamydial activity was evaluated in Chlamydia abortus. They both showed dose-dependent antichlamydial effects. Significant growth inhibitory effects were observed at a minimum concentration of 40 μM.

Computational Modeling of Supramolecular Metallo-organic Cages-Challenges and Opportunities

Self-assembled metallo-organic cages have emerged as promising biomimetic platforms that can encapsulate whole substrates akin to an enzyme active site. Extensive experimental work has enabled access to a variety of structures, with a few notable examples showing catalytic behavior. However, computational investigations of metallo-organic cages are scarce, not least due to the challenges associated with their modeling and the lack of accurate and efficient protocols to evaluate these systems. In this review, we discuss key molecular principles governing the design of functional metallo-organic cages, from the assembly of building blocks through binding and catalysis. For each of these processes, computational protocols will be reviewed, considering their inherent strengths and weaknesses. We will demonstrate that while each approach may have its own specific pitfalls, they can be a powerful tool for rationalizing experimental observables and to guide synthetic efforts. To illustrate this point, we present several examples where modeling has helped to elucidate fundamental principles behind molecular recognition and reactivity. We highlight the importance of combining computational and experimental efforts to speed up supramolecular catalyst design while reducing time and resources.

Supramolecular Mitigation of the Cyanine Limit Problem

Currently, there is a substantial research effort to develop near-infrared fluorescent polymethine cyanine dyes for biological imaging and sensing. In water, cyanine dyes with extended conjugation are known to cross over the "cyanine limit" and undergo a symmetry breaking Peierls transition that favors an unsymmetric distribution of π-electron density and produces a broad absorption profile and low fluorescence brightness. This study shows how supramolecular encapsulation of a newly designed series of cationic, cyanine dyes by cucurbit[7]uril (CB7) can be used to alter the π-electron distribution within the cyanine chromophore. For two sets of dyes, supramolecular location of the surrounding CB7 over the center of the dye favors a nonpolar ground state, with a symmetric π-electron distribution that produces a sharpened absorption band with enhanced fluorescence brightness. The opposite supramolecular effect (i.e., broadened absorption and partially quenched fluorescence) is observed with a third set of dyes because the surrounding CB7 is located at one end of the encapsulated cyanine chromophore. From the perspective of enhanced near-infrared bioimaging and sensing in water, the results show how that the principles of host/guest chemistry can be employed to mitigate the "cyanine limit" problem.

Guaianolides and unusual 3-oxa-guaianolides from Artemisia macrocephala

Eight undescribed guaianolides (macrocephalolides A-H) and two known analogues (parishin C, artabsinolide E) were isolated from the whole plant of Artemisia macrocephala growing in Xinjiang, China. Their structures were determined on the basis of extensive spectroscopic analysis, with absolute configurations established by comparison of experimental and calculated ECD data, as well as confirmation of single-crystal X-ray diffraction crystallography. Macrocephalolides A-B featured an unusual type of 3-oxa-guaianolide with a cyclopentenone moiety. Macrocephalolides C-E possessed a dihydro-2H-pyran acetal segment, representing an unprecedented 2, 3-secoguaianolide skeleton with 6/7/5 tricyclic ring system in natural sesquiterpenes. The X-ray crystal structures of parishin C and artabsinolide E were reported for the first time.

Searching for a Reliable Density Functional for Molecule-Environment Interactions, Found B97M-V/def2-mTZVP

In the present study, we have examined density functional theory methods for the calculation of the interaction energy between a small molecule and its environment. For simple systems such as a neutral solute in a neutral solvent, good accuracy can be attained using low-cost "3c" methods, in particular r²SCAN-3c. When part(s) of the system is charged, the accurate computation of the interactions is more challenging. In these cases, we find the B97M-V/def2-mTZVP method to agree well with reference values; it also shows good accuracy for the more straightforward neutral systems. Thus, B97M-V/def2-mTZVP provides a means for accurate and low-cost computation of interaction energies, notably the binding between a substrate or a drug molecule and an enzyme, which may facilitate rational drug design.

Investigation of the Nature of Intermolecular Interactions in Tetra(thiocyanato)corrolato-Ag(III) Complexes: Agostic or Hydrogen Bonded?

Tetra(thiocyanato)corrolato-Ag(III) complexes presented here constitute a new class of metallo-corrole complexes. The spectroscopic properties of these complexes are quite unusual and interesting. For example, the absorption spectra of these β-substituted corrolato-Ag(III) complexes are very different from those of the β-unsubstituted corrolato-Ag(III) derivatives. Single-crystal XRD analysis of a representative tetra(thiocyanato)corrolato-Ag(III) derivative reveals C-H···Ag interactions. The C-H···Ag interactions are rarely demonstrated in the crystal lattice of a discrete coordination/organometallic compound. Optimization of the hydrogen positions of the crystal structure discloses the geometrical parameters of the said interaction as a Ag···H distance of 2.597 Å and ∠C-H···Ag of 109.62°. The natural bond orbital analysis provides information about the donor-acceptor orbitals involved in the interactions and their interaction energies. It was observed that the σ_C-H orbital overlaps with the vacant d-orbital of Ag with an interaction energy of 17.93 kJ/mol. The filled d-orbital of Ag overlaps with the σ*_C-H orbital with an interaction energy of 4.79 kJ/mol. The highlights of this work are that the H···Ag distance is outside of the distance range for the typical agostic interaction but fitted with the weak H-bond distance. However, the ∠C-H···Ag angle is within the range of the agostic interaction. Both crystallographic data and electronic structure calculations reveal that these kinds of intermolecular interactions in square-planar d⁸ Ag(III) complexes are intermediate in nature. Thus, they cannot be categorically called either hydrogen bonding or agostic interaction.

Open-Shell Variant of the London Dispersion-Corrected Hartree-Fock Method (HFLD) for the Quantification and Analysis of Noncovalent Interaction Energies

The London dispersion (LD)-corrected Hartree–Fock (HF) method (HFLD) is an ab initio approach for the quantification and analysis of noncovalent interactions (NCIs) in large systems that is based on the domain-based local pair natural orbital coupled-cluster (DLPNO-CC) theory. In the original HFLD paper, we discussed the implementation, accuracy, and efficiency of its closed-shell variant. Herein, an extension of this method to open-shell molecular systems is presented. Its accuracy is tested on challenging benchmark sets for NCIs, using CCSD(T) energies at the estimated complete basis set limit as reference. The HFLD scheme was found to be as accurate as the best-performing dispersion-corrected exchange-correlation functionals, while being nonempirical and equally efficient. In addition, it can be combined with the well-established local energy decomposition (LED) for the analysis of NCIs, thus yielding additional physical insights.

Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials

Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.

Hyper-CEST NMR of metal organic polyhedral cages reveals hidden diastereomers with diverse guest exchange kinetics

Guest capture and release are important properties of self-assembling nanostructures. Over time, a significant fraction of guests might engage in short-lived states with different symmetry and stereoselectivity and transit frequently between multiple environments, thereby escaping common spectroscopy techniques. Here, we investigate the cavity of an iron-based metal organic polyhedron (Fe-MOP) using spin-hyperpolarized 129Xe Chemical Exchange Saturation Transfer (hyper-CEST) NMR. We report strong signals unknown from previous studies that persist under different perturbations. On-the-fly delivery of hyperpolarized gas yields CEST signatures that reflect different Xe exchange kinetics from multiple environments. Dilute pools with ~ 104-fold lower spin numbers than reported for directly detected hyperpolarized nuclei are readily detected due to efficient guest turnover. The system is further probed by instantaneous and medium timescale perturbations. Computational modeling indicates that these signals originate likely from Xe bound to three Fe-MOP diastereomers (T, C3, S4). The symmetry thus induces steric effects with aperture size changes that tunes selective spin manipulation as it is employed in CEST MRI agents and, potentially, impacts other processes occurring on the millisecond time scale.

Benchmarking Semiempirical QM Methods for Calculating the Dipole Moment of Organic Molecules

The dipole moment is a simple descriptor of the charge distribution and polarity and is important for understanding and predicting various molecular properties. Semiempirical (SE) methods offer a cost-effective way to calculate dipole moment that can be used in high-throughput screening applications although the accuracy of the methods is still in question. Therefore, we have evaluated AM1, GFN0-xTB, GFN1-xTB, GFN2-xTB, PM3, PM6, PM7, B97-3c, HF-3c, and PBEh-3c SE methods, which cover a variety of SE approximations, to directly assess the performance of SE methods in predicting molecular dipole moments and their directions using 7211 organic molecules contained in the QM7b database. We find that B97-3c and PBEh-3c perform best against coupled-cluster reference values yielding dipole moments with a mean absolute error (MAE) of 0.10 and 0.11 D, respectively, which is similar to the MAE of density functional theory (DFT) methods (∼0.1 D) reported earlier. Analysis of the atomic composition shows that all SE methods show good performance for hydrocarbons for which the spread of error was within 1 D of the reference data, while the worst performances are for sulfur-containing compounds for which only B97-3c and PBEh-3c show acceptable performance. We also evaluate the effect of SE optimized geometry, instead of the benchmark DFT geometry, that shows a dramatic drop in the performance of AM1 and PM3 for which the range of error tripled. Based on our overall findings, we highlight that there is an optimal compromise between accuracy and computational cost using GFN2-xTB (MAE: 0.25 D) that is 3 orders of magnitude faster than B97-3c and PBEh-3c. Thus, we recommend using GFN2-xTB for cost-efficient calculation of the dipole moment of organic molecules containing C, H, O, and N atoms, whereas, for sulfur-containing organic molecules, we suggest PBEh-3c.

What quantum chemical simulations tell us about the infrared spectra, structure and interionic interactions of a bulk ionic liquid

The recently developed efficient protocol to explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is applied to ionic liquid 1-ethyl-3-methyl-imidazolium tetrafluoroborate [Emim][BF₄], and its C2-deuterated analog [Emim-d][BF₄]. It is shown that the solvation strongly modifies the frequencies and IR intensities of both cationic and anionic components of the ionic liquids. The main features of the bulk spectra are reproduced by the simulations for cluster ([Emim][BF₄])₈, representing an ion pair solvated by the first solvation shell. The geometry of the cluster closely resembles the solid-state structure of the actual ionic liquid and is characterized by short contacts of all CH moieties of the imidazolium ring with [BF₄]^- anions. Both structural and spectroscopic analyses allow the contacts to be interpreted as hydrogen bonds of approximately equal strength. The enthalpies of these liquid-state H-bonds, estimated with the use of empirical correlations, amount to 1.2-1.5 kcal mol^-1, while the analogous estimates obtained for the gas-phase charged species [Emim][BF₄]₂^- and [Emim]₂[BF₄]⁺ increase to 3.6-3.9 kcal mol^-1.

Thermally Stable Array of Discrete C60s on a Two-Dimensional Crystalline Adlayer of Macrocycles both in Vacuo and under Ambient Pressure

A periodic monolayer array of discrete C₆₀s was generated on an atomically flat Au(111) surface with the aid of a template adlayer. The template was a two-dimensional (2D) array of molecular pits prepared on an Au(111) surface through 2D crystallization of shape-persistent macrocycles composed of four carbazole and four salphens/Ni-salphens with a 1 nm hollow. Scanning tunneling microscopy imaging under ultra-high vacuum revealed that the square-shaped macrocycles, with 1.5 nm sides, were arranged with a periodic spacing of approximately 4.0 nm on the Au(111) surface, where the orientation and periodicity of the macrocycles were dependent on their chemical structures. After sublimation of C₆₀s onto the adlayer, a single C₆₀ molecule was entrapped in each pit, and an ordered molecular array of C₆₀s was attained with a pattern similar to that of the macrocycles. The periodic pattern of C₆₀s on the surface was thermally stable up to approximately 200 °C, even under ambient pressure. Scanning tunneling spectroscopy suggested the existence of an electronic interaction between the C₆₀s and the Au(111) surface that was influenced by the macrocycle template on the surface.

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Unlocking the computational design of metal-organic cages

Metal-organic cages are macrocyclic structures that can possess an intrinsic void that can hold molecules for encapsulation, adsorption, sensing, and catalysis applications. As metal-organic cages may be comprised from nearly any combination of organic and metal-containing components, cages can form with diverse shapes and sizes, allowing for tuning toward targeted properties. Therefore, their near-infinite design space is almost impossible to explore through experimentation alone and computational design can play a crucial role in exploring new systems. Although high-throughput computational design and screening workflows have long been known as powerful tools in drug and materials discovery, their application in exploring metal-organic cages is more recent. We show examples of structure prediction and host-guest/catalytic property evaluation of metal-organic cages. These examples are facilitated by advances in methods that handle metal-containing systems with improved accuracy and are the beginning of the development of automated cage design workflows. We finally outline a scope for how high-throughput computational methods can assist and drive experimental decisions as the field pushes toward functional and complex metal-organic cages. In particular, we highlight the importance of considering realistic, flexible systems.

Structural, Electronic, and Charge Transport Properties of New Materials based on 2-(5-Mercapto-1,3,4-Oxadiazol-2-yl) Phenol for Organic Solar Cells and Light Emitting Diodes by DFT and TD-DFT

This work reports on the density functional theory (DFT) and its time-dependent extension (TD-DFT) study of the structural, electronic, and charge transport properties of 2-(5-mercapto-1,3,4-oxadiazol-2-yl) phenol (MODP) and some of its transition M2+ complexes (M = Fe, Co, Cu, Ni, Zn, Pd, Pt). Reorganization energy, integral charge transfer, mobility, open circuit voltage, and electronic properties of these compounds have been calculated by employing the global hybrid functional PBE0 in conjunction with the Karlsruhe basis set def2-TZVP. Results show that MODP and its transition metal complexes are good electron donors for organic solar cells (OSC) owing to their relatively higher HOMO and LUMO energies compared to the prototypical (6, 6)-phenyl-C61-butyric acid methyl ester (PCBM). Energy gaps ranging between 2.502 and 4.455 eV, energy driving forces (∆EL-L) ranging between 2.08 and 2.44 eV, and large open circuit voltages ( $V_{\text{OC}}$ ) ranging from 1.12 to 2.05 eV were obtained. The results also revealed that MODP and its Pd(II) and Pt(II) complexes could serve as ambipolar charge transport materials owing to their very small reorganization energies, integral charge transfers, high rate charge transfers, and mobilities. All studied molecules showed OSC donor and hole/electron transport characteristics required by organic light-emitting diodes (OLEDs). Based on these results, new ways for designing charge transport materials for OLEDs as well as donor materials in OSCs are proposed.

Strategies for Computer-Aided Discovery of Novel Open-Shell Polymers

Organic π-conjugated polymers with a triplet ground state have been the focus of recent research for their interesting and unique electronic properties, arising from the presence of the two unpaired electrons. These compounds are usually built from alternating electron-donating and electron-accepting monomer pairs which lower the HOMO-LUMO gap and yield a triplet state instead of the typical singlet ground state. In this paper, we use density functional theory calculations to explore the design rules that govern the creation of a ground-state triplet conjugated polymer and find that a small HOMO-LUMO gap in the singlet state is the best predictor for the existence of a triplet ground state, compared to previous use of a pro-quinoidal bonding character. This work can accelerate the discovery of new stable triplet materials by reducing the computational resources needed for electronic-state calculations and the number of potential candidates for synthesis.

Modeling the Binding Free Energy of Large Atmospheric Sulfuric Acid-Ammonia Clusters

Sulfuric acid and ammonia are believed to account for a large fraction of new-particle formation in the atmosphere. However, it remains unclear how small clusters grow to larger sizes, eventually ending up as stable aerosol particles. Here we present the largest sulfuric acid-ammonia clusters studied to date using quantum chemical methods by calculating the binding free energies of (SA) _n (A) _n clusters, with n up to 20. Based on benchmark calculations, we apply the B97-3c//GFN1-xTB level of theory to calculate the cluster structures and thermochemical parameters. We find that the cluster structures drastically evolve at larger sizes. We identify that an ammonium ion is fully coordinated in the core of the cluster at n = 7, and at n = 13 we see the emergence of the first fully coordinated bisulfate ion. We identify multiple ammonium and bisulfate ions that are embedded in the core of the cluster structure at n = 19. The binding free energy per acid-base pair levels out around n = 8-10, indicating that at a certain point the thermochemistry of the clusters converges toward a constant value.

Trigonal Bipyramidal Rhodium(I) Methyl and Phenyl Complexes: Precursors of Oxidative Methyl and Phenyl Radical Generation

The new complexes [Rh(Me)(trop3P)] (2) and [Rh(Ph)(trop3P)] (3) (trop = 5H-dibenzo[a,d]cyclohepten-5-yl) were synthesised by addition of organolithium reagents (MeLi and PhLi) to the parent pentacoordinated chloride complex [RhCl(trop3P)]. The compounds have a trigonal bipyramidal structure with olefin-only ligands in the equatorial position and the methyl or phenyl substituent in the axial position. Oxidation of complexes 2 and 3 leads to the liberation of methyl and phenyl radicals, which were indirectly detected by reaction with common spin trapping reagents.

Quantification of Geometric Errors Made Simple: Application to Main-Group Molecular Structures

Nearly all electronic structure simulations begin with obtaining approximate geometries, making a systematic quantification of errors in approximate molecular structures of key importance. Recently, the geometric energy offset (GEO) framework based on a single and natural measure for quantifying and analyzing these errors has been proposed ( J. Phys. Chem. Lett. 2020, 11, 99579964). An accurate and far less costly approximation to GEO is utilized here to readily quantify errors in main-group structures and analyze them in a chemically intuitive way. The use of semiexperimental geometries as a reference further simplifies the analysis. The analysis reveals new insights into the geometric performance of methods, their rankings, as well as patterns across different classes of methods and basis sets that arise from the analysis.

Benchmark Study on the Calculation of 119Sn NMR Chemical Shifts

A new benchmark set termed SnS51 for assessing quantum chemical methods for the computation of ¹¹⁹Sn NMR chemical shifts is presented. It covers 51 unique ¹¹⁹Sn NMR chemical shifts for a selection of 50 tin compounds with diverse bonding motifs and ligands. The experimental reference data are in the spectral range of ±2500 ppm measured in seven different solvents. Fifteen common density functional approximations, two scalar- and one spin-orbit relativistic approach are assessed based on conformer ensembles generated using the CREST/CENSO scheme and state-of-the-art semiempirical (GFN2-xTB), force field (GFN-FF), and composite DFT methods (r²SCAN-3c). Based on the results of this study, the spin-orbit relativistic method combinations of SO-ZORA with PBE0 or revPBE functionals are generally recommended. Both yield mean absolute deviations from experimental data below 100 ppm and excellent linear regression determination coefficients of ≤0.99. If spin-orbit calculations are not affordable, the use of SR-ZORA with B3LYP or X2C with ωB97X or M06 may be considered to obtain qualitative predictions if no severe spin-orbit effects, for example, due to heavy nuclei containing ligands, are expected. An empirical linear scaling correction is demonstrated to be applicable for further improvement, and respective empirical parameters are given. Conformational effects on chemical shifts are studied in detail but are mostly found to be small. However, in specific cases when the ligand sphere differs substantially between conformers, chemical shifts can change by up to several hundred ppm.

Attaching Onto or Inserting Into an Intramolecular Hydrogen Bond: Exploring and Controlling a Chirality-Dependent Dilemma for Alcohols

Prereactive complexes in noncovalent organocatalysis are sensitive to the relative chirality of the binding partners and to hydrogen bond isomerism. Both effects are present when a transiently chiral alcohol docks on a chiral α-hydroxy ester, turning such 1:1 complexes into elementary, non-reactive model systems for chirality induction in the gas phase. With the help of linear infrared and Raman spectroscopy in supersonic jet expansions, conformational preferences are investigated for benzyl alcohol in combination with methyl lactate, also exploring p-chlorination of the alcohol and the achiral homolog methyl glycolate to identify potential London dispersion and chirality effects on the energy sequence. Three of the four combinations prefer barrierless complexation via the hydroxy group of the ester (association). In contrast, the lightest complex predominantly shows insertion into the intramolecular hydrogen bond, such as the analogous lactate and glycolate complexes of methanol. The experimental findings are rationalized with computations, and a uniform helicality induction in the alcohol by the lactate is predicted, independent of insertion into or association with the internal lactate hydrogen bond. p-chlorination of benzyl alcohol has a stabilizing effect on association because the insertion motif prevents a close contact between the chlorine and the hydroxy ester. After simple anharmonicity and substitution corrections, the B3LYP-D3 approach offers a fairly systematic description of the known spectroscopic data on alcohol complexes with α-hydroxy esters.

Computer-aided simulation of infrared spectra of ethanol conformations in gas, liquid and in CCl4 solution

The recently developed efficient protocol combining implicit and explicit, accurate quantum-mechanical modeling of the condensed state (Katsyuba et al., J. Chem. Phys. 155, 024507 [2021]) is used to describe the IR spectra of liquid ethanol and its solutions in CCl₄ . The relative abundance of the anti and gauche conformers of ethanol is shown to increase from ~40:60 in the gas phase to ~55:45 in the liquid phase. In spite of a moderate impact of media effects on the conformational composition of the liquid, the solvent strongly influences vibrational frequencies, IR intensities, and normal modes of each conformer, producing qualitatively different spectra compared to the gas phase and CCl₄ solution. Further, these solvent effects affecting IR frequencies and intensities depend not only on the conformation of the solvated molecule but also on the solvating species. Nevertheless, vibrational frequencies of anti and gauche conformers of liquid ethanol and its several isotopomers practically coincide with each other. Convenient liquid-state conformational markers in the fingerprint region of IR spectra are revealed for the hydroxyl-deuterated species: CH₃ CH₂ OD, CH₃ CHDOD, CH₃ CD₂ OD, and CD₃ CD₂ OD.

Organic Emitters Showing Excited-States Energy Inversion: An Assessment of MC-PDFT and Correlation Energy Functionals Beyond TD-DFT

The lowest-energy singlet (S1) and triplet (T1) excited states of organic conjugated chromophores are known to be accurately calculated by modern wavefunction and Time-Dependent Density Functional Theory (TD-DFT) methods, with the accuracy of the latter heavily relying on the exchange-correlation functional employed. However, there are challenging cases for which this cannot be the case, due to the fact that those excited states are not exclusively formed by single excitations and/or are affected by marked correlation effects, and thus TD-DFT might fall short. We will tackle here a set of molecules belonging to the azaphenalene family, for which research recently documented an inversion of the relative energy of S1 and T1 excited states giving rise to a negative energy difference (ΔEST) between them and, thereby, contrary to most of the systems thus far treated by TD-DFT. Since methods going beyond standard TD-DFT are not extensively applied to excited-state calculations and considering how challenging this case is for the molecules investigated, we will prospectively employ here a set of non-standard methods (Multi-Configurational Pair Density Functional Theory or MC-PDFT) and correlation functionals (i.e., Lie–Clementi and Colle–Salvetti) relying not only on the electronic density but also on some modifications considering the intricate electronic structure of these systems.

Increased Antiaromaticity through Pentalene Connection in [n]Cyclo-1,5-dibenzopentalenes

Conjugated nanohoops incorporating nonalternant hydrocarbons have altered optoelectronic properties compared to [n]cycloparaphenylenes or other purely aromatic hoops. We synthesized [n]cyclo-1,5-dibenzopentalenes (n = 4, 5), in which nonalternant dibenzo[a,e]pentalenes are connected through their pentalene units. This leads to an increase in antiaromatic character and low-lying LUMO energies. Calculations show puckered or entangled conformations of the precursor macrocyclic Pt-complexes. Our study proves dibenzopentalene as a versatile nonalternant building block for conjugated nanohoops with modifiable antiaromaticity and optoelectronic properties.

Synthesis and Computationally Assisted Spectroscopic Study of Tautomerism in 3-(Phenyl(2-arylhydrazineylidene)methyl)quinoxalin-2(1H)-ones

The recently developed efficient protocol combining implicit and explicit, accurate quantum mechanical modeling of the condensed state [Katsyuba et al., J. Chem. Phys. 155, 024507 (2021)] is used to describe...