An efficient and near linear scaling pair natural orbital based local coupled cluster method

Pyrolytic conversion of glucose into hydroxymethylfurfural and furfural: Benchmark quantum-chemical calculations

Quantum chemical methods have been intensively applied to study the pyrolytic conversion of glucose into hydroxymethylfurfural (HMF) and furfural (FF). Herein, we collect the most relevant mechanistic proposals from the recent literature and organize them into a single reaction network. All the transition structures (TSs) and intermediates are characterized using highly accurate ab initio methods and the possible reaction pathways are assessed in terms of the Gibbs energies of the TSs and intermediates with respect to β-glucopyranose, selecting a 2D ideal-gas standard state at 773 K to represent the pyrolysis conditions. Several pathways can lead to the formation of both HMF and FF passing through rate-determining TSs that have ΔG‡ values of ~49-50 kcal/mol. Both water-assisted mechanisms and nonspecific environmental effects have a minor impact on the Gibbs energy profiles. We find that the HMF → FF + CH2O fragmentation has a small ΔrxnG value and an accessible ΔG‡ barrier. Our computational results, which are in consonance with the kinetic parameters derived from lumped models, the results of isotopic labeling experiments and the reported HMF/FF molecular ratios, could be useful for modeling studies including on nonequilibrium kinetic effects that may render more information about product yields and the relevance of the various pathways.

Generation and interception of bicyclo[3.2.1]oct-2-yne: an experimental and theoretical mechanistic study

Bicyclo[3.2.1]oct-2-yne was generated from the Fritsch-Buttenberg-Wiechell rearrangement of 2-norbornylidene carbene. The rearrangement preferentially involves migration of a tertiary carbon over a secondary carbon, a trend that contrasts with rearrangements of acyclic carbenes and which may be attributable to hyperconjugative effects promoted by the bridged structure of the carbene.

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation

This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.

Enzymatic synthesis of azide by a promiscuous N-nitrosylase

Azides are energy-rich compounds with diverse representation in a broad range of scientific disciplines, including material science, synthetic chemistry, pharmaceutical science and chemical biology. Despite ubiquitous usage of the azido group, the underlying biosynthetic pathways for its formation remain largely unknown. Here we report the characterization of an enzymatic route for de novo azide construction. We demonstrate that Tri17, a promiscuous ATP- and nitrite-dependent enzyme, catalyses organic azide synthesis through sequential N-nitrosation and dehydration of aryl hydrazines. Through biochemical, structural and computational analyses, we further propose a plausible molecular mechanism for azide synthesis that sets the stage for future biocatalytic applications and biosynthetic pathway engineering.

A refreshing approach to understanding the action on DNA of vanadium (IV) and (V) complexes derived from the anticancer VCp2Cl2

A computational study based on derivatives of the anticancer VCp2Cl2 compound and their interaction with representative models of deoxyribonucleic acid (DNA) is presented. The derivatives were obtained by substituting the cyclopentadienes of VCp2Cl2 with H2O, NH3, OH-, Cl-, O2- and C2O42- ligands. The oxidation states IV and V of vanadium were considered, so a total of 20 derivative complexes are included. The complexes interactions with DNA were studied using two different models, the first model considers the interactions of the complexes with the pair Guanine-Cytosine (G-C) and the second involves the interaction of the complexes with adjacent pairs, that is, d(GG). This study compares methodologies based on density functional theory with coupled cluster like calculations (DLPNO-CCSD(T)), the gold standard of electronic structure methods. Furthermore, the change in the electron density of the hydrogen bonds that keep bonded the G-C pair and d(GG) pairs, due to the presence of vanadium (IV) and (V) complexes is rationalize. To this aim, quantities obtained from the topology of the electron densities are inspected, particularly the value of the electron density at the hydrogen bond critical points. The approach allowed to identify vanadium complexes that lead to significant changes in the hydrogen bonds indicated above, a key aspect in the understanding, development, and proposal of mechanisms of action between metal complexes and DNA.

Solvent-Directed Social Chiral Self-Sorting in Pd2L4 Coordination Cages

A family of Pd2L4 cages prepared from ligands based on an axially chiral diamino-[1,1'-biazulene] motif (serving as a unique azulene-based surrogate of the ubiquitous BINOL moiety) is reported. We show that preparing a cage starting from the racemate of a shorter bis-monodentate ligand derivative, equipped with pyridine donor groups, leads to integrative ("social") chiral self-sorting, exclusively yielding the meso-trans product, but only in a selection of solvents. This phenomenon is driven by individual solvent molecules acting as hydrogen bonding tethers between the amino groups of neighboring ligands, thereby locking the final coordination cage in a single isomeric form. The experimental (solvent-dependent NMR, single-crystal X-ray diffraction) observations of this cooperative interaction could be explained by computational analyses only when explicit solvation was considered. Furthermore, we prepared a larger chiral ligand with isoquinoline donors, which, unlike the first one, does not undergo social self-sorting from its racemic mixture, further highlighting the importance of solvents bridging short distances between the amino groups. Homochiral cages formed from this larger ligand, however, furnish a cavity that can bind anionic and neutral metal complexes such as [Pt(CN)6]2- and Cr(CO)6 and discriminate between the two enantiomers of chiral guest camphor sulfonate.

Finite-field Cholesky decomposed coupled-cluster techniques (ff-CD-CC): theory and application to pressure broadening of Mg by a He atmosphere and a strong magnetic field

For the interpretation of spectra of magnetic stellar objects such as magnetic white dwarfs (WDs), highly accurate quantum chemical predictions for atoms and molecules in finite magnetic field are required. Especially the accurate description of electronically excited states and their properties requires established methods such as those from coupled-cluster (CC) theory. However, respective calculations are computationally challenging even for medium-sized systems. Cholesky decomposition (CD) techniques may be used to alleviate memory bottlenecks. In finite magnetic field computations, the latter are increased due to the reduction of permutational symmetry within the electron-repulsion-integrals (ERIs) as well as the need for complex-valued data types. CD enables a memory-efficient, approximate description of the ERIs with rigorous error control and thus the treatment of larger systems at the CC level becomes feasible. In order to treat molecules in a finite magnetic field, we present in this work the working equations of the left and right-hand side equations for finite field (ff)-EOM-CD-CCSD for various EOM flavours as well as for the approximate ff-EOM-CD-CC2 method. The methods are applied to the study of the modification of the spectral lines of a magnesium transition by a helium atmosphere that can be found on magnetic WD stars.

Tricyanomethane or Dicyanoketenimine-Silylation makes the Difference

Pseudohalides such as tricyanomethanide, [C(CN)3]-, are well known in chemistry, biochemistry and industrial chemistry. The protonated species HC(CN)3, a classic hydrogen pseudohalide Brønsted acid, is a very strong acid with a pKa value of -5. However, HC(CN)3 is difficult to handle as it tends to decompose rapidly or, more precisely, to oligo- and polymerize. Therefore, silylated pseudohalide compounds with the [Me3Si]+ as the "big organometallic proton" have become interesting, exhibiting similar chemical properties but better kinetic protection. Here, the stepwise silylation of the pseudohalide anion [C(CN)3]- is reported, forming the heavier homologue of HC(CN)3, namely [Me3Si][C(CN)3], and in presence of two additional [Me3Si]+ cations even the dicationic species [(Me3Si-NC)3C]2+ as stable [B(C6F5)4]- salt. Surprisingly, in contrast to the protonated species HC(CN)3, in which the proton is bound to the central carbon atom of [C(CN)3]-, silylation of the [C(CN)3]- anion occurs at one of the three terminal nitrogen atoms, thus forming the long-sought dicyanoketenimine [Me3Si-NC-C(CN)2]. All further silylation steps take place exclusively on the terminal N atoms of the three CN groups and not on the central carbon atom, until the intriguing, highly symmetrical dication, [(Me3Si-NC)3C]2+, is finally generated. The experimental data are supported by quantum chemical calculations in terms of thermodynamics and chemical bonding.

Calculated ionization energies, orbital eigenvalues (HOMO), and related QSAR descriptors of organic molecules: a set of 61 experimental values enables elimination of systematic errors and provides realistic error estimates

Ionization energies (IEs) of organic compounds come in different forms-adiabatic, vertical, as electrode potentials, or as orbital eigenvalues via Koopmans' theorem. They have been linked to the reactivity towards electrophiles and have been used to quantitatively describe electron transfer processes. The de novo prediction of IEs is only meaningful when an estimate of the prediction's uncertainty is included. Pulsed-field ionization (PFI) experiments have quantified adiabatic IEs with unprecedented precision. In this work, a new set of PFI-derived IEs is compiled from the literature as a benchmark for prediction methods. This set includes many common functional groups, a size range from diatomics to two aromatic rings, and IEs between 7 and 14 eV. The first-principles CCSD(T)/CBS protocol presently used reproduces these values within 0.05 eV. For adiabatic IEs and vertical IEs/orbital eigenvalues predicted using approximate density functional theory (DFT), linear regression models are proposed, so that IEs calculated using different methods can be directly compared on a physical scale. This elimination of systematic errors improves the error statistics and allows the performance of predicted IEs to be evaluated if used in quantitative structure-property or -activity relationships, as the latter implicitly correct a descriptor's bias. Owing to the structural scope of the test set, the minimum and maximum deviations from experiment should correspond to those expected for common organic molecules. Deviations from reference values found for orbital eigenvalues but also for IEs calculated explicitly with HF or semi-empirical MO methods were as large as 0.5 eV to 2.0 eV. Such large errors could also propagate into quantitative structure-property models, as shown in illustrative examples of oxidation rate constants in solution.

Spectroscopic and Quantum Chemical Evidence of Amine-CO2 and Alcohol-CO2 Interactions: Confirming an Intriguing Affinity of CO2 to Monoethanolamine (MEA)

A recent broadband rotational spectroscopic investigation of the cross-association mechanisms of CO2 with monoethanolamine (MEA) in molecular beams [F. Xie et al., Angew. Chem. Int. Ed., 2023, 62, e202218539] revealed an intriguing affinity of CO2 to the hydroxy group. These findings have triggered the present systematic vibrational spectroscopic exploration of weakly bound amine··CO2 and alcohol··CO2 van der Waals cluster molecules embedded in inert "quantum" matrices of neon at 4.2 K complemented by high-level quantum chemical conformational analyses. The non-covalent interactions formed between the amino and hydroxy groups and the electron-deficient carbon atom of CO2 are demonstrated to lift the degeneracy of the doubly degenerate intramolecular CO2-bending fundamental significantly with characteristic observed spectral splittings for the amine··CO2 (≈35-45 cm-1) and alcohol··CO2 (≈20-25 cm-1) interactions, respectively, despite the almost identically predicted total association energies (≈12-14 kJ·mol-1) for these van der Waals contacts, as revealed by benchmark Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory. These high-level theoretical predictions reveal significantly higher "geometry preparation energies" for the amine··CO2 systems leading to a more severe distortion of the CO2 linearity upon complexation in agreement with the infrared spectroscopic findings. The systematic combined spectroscopic and quantum chemical evidences for cross-association between CO2 and amines/alcohols in the present work unambiguously confirm an intriguing binding preference of CO2 to the hydroxy group of the important carbon capture agent MEA, with an accurate vibrational zero-point energy corrected association energy (D0) of 13.5 kJ·mol-1 at the benchmark DLPNO-CCSD(T)/aug-cc-pV5Z level of theory.

Rationalizing polymorphism with local correlation-based methods: a case study of pnictogen molecular crystals

A computational workflow is proposed to quantify and rationalize the relative stability of different structures of molecular crystals using cluster models and quantum chemical methods. The Hartree-Fock plus London Dispersion (HFLD) scheme is used to estimate the lattice energy of molecular crystals in various structural arrangements. The fragment-pairwise Local Energy Decomposition (fp-LED) scheme is then employed to quantify the key intermolecular interactions responsible for the relative stability of different crystal structures. The fp-LED scheme provides also in-depth chemical insights by decomposing each interaction into energy components such as dispersion, electrostatics, and exchange. Notably, this analysis requires only a single interaction energy computation per structure on a suitable cluster model. As a case study, two polymorphs of each of the following are considered: naphthyl-substituted dipnictanes (with As, Sb, and Bi as the pnictogen atom) and tris(thiophen-2-yl)bismuthane. The approach outlined offers high accuracy as well as valuable insights for developing design principles to engineer crystal structures with tailored properties, opening up new avenues in the study of molecular aggregates, potentially impacting diverse fields in materials science and beyond.

Cluster perturbation theory. X. A parallel implementation of Lagrangian perturbation series for the coupled cluster singles and doubles ground-state energy through fifth order

We describe an efficient implementation of cluster perturbation and Møller-Plesset Lagrangian energy series through the fifth order that targets the coupled cluster singles and doubles energy utilizing the resolution of the identity approximation. We illustrate the computational performance of the implementation by performing ground state energy calculations on systems with up to 1200 basis functions using a single node and by comparison to conventional coupled cluster singles and doubles calculations. We further show that our hybrid message passing interface/open multiprocessing parallel implementation that also utilizes graphical processing units can be used to obtain fifth order energies on systems with almost 1200 basis functions with a 90 min "time to solution" running on Frontier at Oak Ridge National Laboratory.

"Best" Iterative Coupled-Cluster Triples Model? More Evidence for 3CC

To follow up on the unexpectedly good performance of several coupled-cluster models with approximate inclusion of 3-body clusters [Rishi, V.; Valeev, E. F. J. Chem. Phys. 2019, 151, 064102.] we performed a more complete assessment of the 3CC method [Feller, D. . J. Chem. Phys. 2008, 129, 204105.] for accurate computational thermochemistry in the standard HEAT framework. New spin-integrated implementation of the 3CC method applicable to closed- and open-shell systems utilizes a new automated toolchain for derivation, optimization, and evaluation of operator algebra in many-body electronic structure. We found that with a double-ζ basis set the 3CC correlation energies and their atomization energy contributions are almost always more accurate (with respect to the CCSDTQ reference) than the CCSDT model as well as the standard CCSD(T) model. The mean absolute errors in cc-pVDZ {3CC, CCSDT, and CCSD(T)} electronic (per valence electron) and atomization energies relative to the CCSDTQ reference for the HEAT data set [Tajti, A. . J. Chem. Phys. 2004, 121, 11599-11613.], were {24, 70, 122} μEh/e and {0.46, 2.00, 2.58} kJ/mol, respectively. The mean absolute errors in the complete-basis-set limit {3CC, CCSDT, and CCSD(T)} atomization energies relative to the HEAT model reference, were {0.52, 2.00, and 1.07} kJ/mol, The significant and systematic reduction of the error by the 3CC method and its lower cost than CCSDT suggests it as a viable candidate for post-CCSD(T) thermochemistry applications, as well as the preferred alternative to CCSDT in general.

Formation of atmospheric molecular clusters containing nitric acid with ammonia, methylamine, and dimethylamine

This study investigates the formation of atmospheric molecular clusters containing ammonia (NH3, A), methylamine (CH3NH2, MA), or dimethylamine (CH3NHCH3, DMA) with nitric acid (HNO3, NA) using quantum mechanics. The Atmospheric Cluster Dynamic Code (ACDC) was employed to simulate the total evaporation rate, formation rate, and growth pathways of three types of clusters under dry and hydrated conditions. This study evaluates the enhancing potential of A/MA/DMA for NA-based new particle formation (NPF) at parts per trillion (ppt) levels. The results indicate that A/MA/DMA can enhance NA-based NPF at high nitric acid concentrations and low temperatures in the atmosphere. The enhancing potential of MA is weaker than that of DMA but stronger than that of A. Cluster growth predominantly follows the lowest free energy pathways on the acid-base grid, with the formation of initial acid-base dimers (NA)(A), (NA)(MA), and (NA)(DMA) being crucial. Hydration influences the evaporation rate and formation rate of clusters, especially for initial clusters. When the humidity is at 100%, the formation rate for NA-A, NA-MA, and NA-DMA clusters can increase by approximately 109, 107, and 104-fold compared to the corresponding unhydrated clusters, respectively. These results highlight the significance of nitric acid nucleation in NPF events in low-temperature, high-humidity atmospheres, particularly in regions like China with significant automobile exhaust pollution.

Automatic Orbital Pair Selection for Multilevel Local Coupled-Cluster Based on Orbital Maps

We present an automatic, orbital-map based orbital-pair selection scheme for multilevel local coupled-cluster approaches that exploits the locality of chemical reactions by focusing on the part of the molecule directly involved in the reaction. The previously introduced pair-selected multilevel extension to domain-based local pair natural orbital coupled-cluster with singles, doubles, and semicanonical perturbative triples [DLPNO-CCSD(T0)] partitions the orbital pairs according to relative changes in pair correlation energies [Bensberg, M.; Neugebauer, J. Orbital pair selection for relative energies in the domain-based local pair natural orbital coupled-cluster method. J. Chem. Phys. 2022, 157, 064102. 10.1063/5.0100010]. To this end, maps between localized orbitals are required which in turn require maps between the atoms of structures along reaction paths. So far, these atom maps have been manually determined, which can be a (human) time-consuming procedure. Here, we present an automatic atom mapping algorithm based on the principle of minimum chemical distance that incorporates orientation dependence through dihedral angles. A similar strategy is then introduced to obtain orbital maps, which proves advantageous over the previously used direct orbital selection. Along with a modified orbital pair prescreening, this results in an improved variant of the pair-selected multilevel DLPNO-CCSD(T0) method. The performance of this approach is demonstrated for various reaction types showing a significant efficiency gain and accurate results due to beneficial, systematic error cancellation. The presented method operates in a black-box manner due to its fully automatized algorithms with only the need to specify a single target-accuracy parameter. Additionally, we demonstrate that basis set extrapolation techniques can be applied. In this context, the approach shows deficiencies for the use of large basis sets, especially with diffuse basis functions, which can be traced back to the semicanonical triples correction.

Progress in Multiscale Modeling of Silk Materials

As a result of their hierarchical structure and biological processing, silk fibers rank among nature's most remarkable materials. The biocompatibility of silk-based materials and the exceptional mechanical properties of certain fibers has inspired the use of silk in numerous technical and medical applications. In recent years, computational modeling has clarified the relationship between the molecular architecture and emergent properties of silk fibers and has demonstrated predictive power in studies on novel biomaterials. Here, we review advances in modeling the structure and properties of natural and synthetic silk-based materials, from early structural studies of silkworm cocoon fibers to cutting-edge atomistic simulations of spider silk nanofibrils and the recent use of machine learning models. We explore applications of modeling across length scales: from quantum mechanical studies on model peptides, to atomistic and coarse-grained molecular dynamics simulations of silk proteins, to finite element analysis of spider webs. As computational power and algorithmic efficiency continue to advance, we expect multiscale modeling to become an indispensable tool for understanding nature's most impressive fibers and developing bioinspired functional materials.

Reduced Scaling Correlated Natural Transition Orbitals for Multilevel Coupled Cluster Calculations

Multilevel coupled cluster theory offers reduced scaling computation of intensive properties in systems that are too large for standard coupled cluster calculations. A significant benefit of the multilevel coupled cluster framework is the possibility of calculating intensive properties that are not tightly localized if an appropriate set of active orbitals is used. Correlated natural transition orbitals (CNTOs) are tailored to describe excitation processes. For multilevel coupled cluster singles and doubles (MLCCSD) and singles and perturbative doubles (MLCC2) calculations, the construction of CNTOs generally becomes the computational bottleneck. Here, we demonstrate how CNTOs can be obtained with O ( N 3 ) operations, eliminating the O ( N 5 ) -scaling steps involved in the original approach. This reduction in scaling moves the bottleneck of MLCC2 and MLCCSD calculations from the active orbital space preparation to the MLCC2 and MLCCSD equations with O ( N 4 ) -scaling.

Comparing coupled cluster and composite quantum chemical methods for computing activation energies and reaction enthalpies of radical propagation reactions

Accurate determination of activation energies and reaction enthalpies is essential for understanding the propagation step in free radical polymerization, as it significantly affects polymer chain length and structure. In this study, we compare DLPNO-CCSD(T) to canonical CCSD(T) for 17 radical addition activation energies and 18 reaction enthalpies from Radom and Fischer's test set. Additionally, we compare the computationally efficient composite methods G3(MP2)-RAD and CBS-RAD against CCSD(T)/aug-cc-pVTZ and DLPNO-CCSD(T)/CBS methods. Compared to the CCSD(T)/aug-cc-pVTZ reference, our results indicate that DLPNO-CCSD(T)/CBS with unrestricted Hartree-Fock (UHF) or UB3LYP reference orbitals and NormalPNO parameters consistently achieves chemical accuracy, with mean absolute deviations of 3.5 kJ mol-1 for activation energies and 1.5 kJ mol-1 for reaction enthalpies. Comparing the two composite methods shows that CBS-RAD agrees most closely with coupled cluster reaction enthalpies, while G3(MP2)-RAD tracks the activation energies most closely.

Design of Anti-Hund Organic Emitters Based on Heptazine

Hund's rule, which is powerful in governing the first excited states of closed-shell organic materials, can hardly be violated to get inverted singlet-triplet gap (INVEST) molecules with negative singlet-triplet energy gaps (ΔEST), although INVEST materials have shown extraordinary photophysical properties and promising device performance especially in light-emitting diodes. Here, we propose a facile strategy to construct emissive INVEST molecules by introducing different types of substituents to heptazine in various modes, which can effectively tune the ΔEST to be negative with the enlarged oscillator strength (f) for the high fluorescence rate of the heptazine derivatives. Systematic computational studies show that the double substitution of electron-donating units with another nonconjugated substituent in hybrid substitution mode is the most favorable way in achieving slightly negative ΔEST and large f values; the conjugated substituent will compete with heptazine to make the molecule deviate from the INVEST feature. Especially, a series of high-performance heptazine-based INVEST emitters were constructed, exhibiting ΔEST low to -0.362 eV, f up to 0.0436, as well as a wide range emission color from 339 to 716 nm. Also, the designed molecules were predicted to have fluorescence radiative rates up to 106 s-1, along with efficient reverse intersystem crossing rates reaching 108 s-1. Importantly, the figure of merit (FM) was first proposed as a parameter to wholly evaluate the performance of INVEST emitters, and the highest FM of 0.198 was found in the triazine and double nonconjugated amine-substituted heptazine. These results highlight the great potential of the heptazine chromophore in constructing INVEST emitters, revealing fundamental structure-property understandings for the material design of efficient anti-Hund organic molecules with improved emission properties.

Intermolecular interaction energies with AROFRAG-A systematic approach for fragmentation of aromatic molecules

Intermolecular interactions with polycyclic aromatic hydrocarbons (PAHs) represent an important area of physisorption studies. These investigations are often hampered by a size of interacting PAHs, which makes the calculation prohibitively expensive. Therefore, methods designed to deal with large molecules could be helpful to reduce the computational costs of such studies. Recently we have introduced a new systematic approach for the molecular fragmentation of PAHs, denoted as AROFRAG, which decomposes a large PAH molecule into a set of predefined small PAHs with a benzene ring being the smallest unbreakable unit, and which in conjunction with the Molecules-in-Molecules (MIM) approach provides an accurate description of total molecular energies. In this contribution we propose an extension of the AROFRAG, which provides a description of intermolecular interactions for complexes composed of PAH molecules. The examination of interaction energy partitioning for various test cases shows that the AROFRAG3 model connected with the MIM approach accurately reproduces all important components of the interaction energy. An additional important finding in our study is that the computationally expensive long-range electron-correlation part of the interaction energy, that is, the dispersion component, is well described at lower AROFRAG levels even without MIM, which makes the latter models interesting alternatives to existing methods for an accurate description of the electron-correlated part of the interaction energy.

A static quantum embedding scheme based on coupled cluster theory

We develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller-Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C-C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

Assessing the Noncovalent Interaction of Deucravacitinib and Ethanol with Special Reference to an Independent Gradient Model Based on Hirshfeld Partition

The current study begins by optimizing the deucravacitinib molecule in the gas phase at the ωB97XD/cc-pVDZ level of theory using density functional theory and proceeds to study its intramolecular interactions. Further, a molecule of EtOH was introduced at different locations on the deucravacitinib molecule, and the noncovalent interactions arising from them were also investigated using several computational tools. In this way, eight deucravacitinib-EtOH systems (1-8) were identified and their electronic environment was studied after evaluating their binding energy. Using natural bond orbital analysis, the localization of charges between the donor and acceptor fragments in these interacting systems was examined. The nature of interactions was analyzed using the reduced gradient approach (NCI analysis), and few hydrogen bonding interactions (intermolecular and intramolecular) were found in each system. The strength of these hydrogen bonding interactions was further investigated by using theoretical tools such as atoms in molecules analysis and independent gradient model based on Hirshfeld partition analysis. The binding energy of deucravacitinib with EtOH was decomposed into energy components based on the domain-based local pair natural orbital coupled cluster technique using LED analysis. The results from the hydrogen bonding interaction analysis using different computational tools were found to be consistent with the calculated order of binding energy of systems 1-8 and they also pointed toward the higher stability of system 3.

Reductive Elimination From Tetra-Alkyl Cuprates [MenCu(CF3)4-n]- (n=0-4): Beyond Simple Oxidation States

In recent years, the electronic structures of organocuprates in general and the complex [Cu(CF3)4]- in particular have attracted significant interest. A possible key indicator in this context is the reactivity of these species. Nonetheless, this aspect has received only limited attention. Here, we systematically study the series of tetra-alkyl cuprates [MenCu(CF3)4-n]- and their unimolecular reactivity in the gas phase, which includes concerted formal reductive eliminations as well as radical losses. Through computational studies, we characterize the electronic structures of the complexes and show how these are connected to their reactivity. We find that all [MenCu(CF3)4-n]- ions feature inverted ligand fields and that the distinct reactivity patterns of the individual complexes arise from the interplay of different effects.

Prospects for rank-reduced CCSD(T) in the context of high-accuracy thermochemistry

Obtaining sub-chemical accuracy (1 kJ mol-1) for reaction energies of medium-sized gas-phase molecules is a longstanding challenge in the field of thermochemical modeling. The perturbative triples correction to coupled-cluster single double triple [CCSD(T)] constitutes an important component of all high-accuracy composite model chemistries that obtain this accuracy but can be a roadblock in the calculation of medium to large systems due to its O(N7) scaling, particularly in HEAT-like model chemistries that eschew separation of core and valence correlation. This study extends the work of Lesiuk [J. Chem. Phys. 156, 064103 (2022)] with new approximate methods and assesses the accuracy of five different approximations of (T) in the context of a subset of molecules selected from the W4-17 dataset. It is demonstrated that all of these approximate methods can achieve sub-0.1 kJ mol-1 accuracy with respect to canonical, density-fitted (T) contributions with a modest number of projectors. The approximation labeled Z̃T appears to offer the best trade-off between cost and accuracy and shows significant promise in an order-of-magnitude reduction in the computational cost of the CCSD(T) component of high-accuracy model chemistries.

Improving Predictions of Spin-Crossover Complex Properties through DFT Calculations with a Local Hybrid Functional

We conducted a study on the performance of the local hybrid exchange-correlation functional PBE0r for a set of 95 experimentally characterized iron spin-crossover (SCO) complexes [Vennelakanti, V.; J. Chem. Phys. 2023, 159, 024120]. The PBE0r functional is a variant of PBE0 where the exchange correction is restricted to on-site terms formulated on the basis of local orbitals. We determine the free parameters of the PBE0r functional against the experimental data and other hybrid functionals. With a Hartree-Fock (HF) exchange factor of 4%, the PBE0r functional accurately reproduces the electronic and free-energy trends predicted in prior DFT studies for these 95 complexes by using the B3LYP functional. Larger values of HF exchange stabilize high-spin states. The PBE0r-predicted bond lengths tend to exceed the experimental bond lengths, although bond lengths are less sensitive to HF exchange than in global hybrids. The predicted SCO transition temperatures T1/2 from PBE0r correlate moderately with the experimental transition temperatures, showing a slight improvement compared to the previous modB3LYP-predicted T1/2. This study suggests that the PBE0r functional is computationally cost-effective and offers the possibility of simulating larger complexes with accuracy comparable to global hybrid functionals, provided the HF-exchange parameter is carefully optimized.

Photoelectron Spectroscopy and Theoretical Studies on the Structures and Properties of Ge4C- and Ge4CH- Clusters

We investigate the structures and properties of Ge4C-/0 and Ge4CH-/0 clusters using anion photoelectron spectroscopy and theoretical calculations. Our calculations show that the first two low-lying isomers coexist in the experiments of Ge4C- and Ge4CH-. The first two low-lying isomers of Ge4C- have trigonal bipyramidal structures with the C atom on the equatorial plane and the top vertex, respectively. It is found that the first two low-lying isomers of Ge4CH- can be obtained by adding an H atom to the top and equatorial C atoms of Ge4C-, respectively. The AdNDP analyses reveal that the C atom in Ge4C forms one 4c-2e σ bond, two 4c-2e π bonds, and one 5c-2e σ bond with Ge atoms. The C atom in Ge4C interacts with an H- forming a C-H σ bond in Ge4CH-. AIMD simulation results indicate high dynamic stabilities of Ge4C and Ge4CH- at 300 and 500 K. Our results show that the structures and chemical bonding of Ge4B- and Ge4N+ are similar to those of Ge4C, while those of Ge4BH2- and Ge4NH resemble those of Ge4CH-.

A comprehensive understanding of the mechanism of the biomimetic total synthesis of brevianamide A

Recently, several studies on the chemical synthesis of brevianamide A (BA) were reported. In particular, a highly efficient and remarkably selective synthetic strategy was reported by Lawrence's group. However, a unified mechanistic understanding of these results is still lacking. We have carried out a DFT study and proposed a unified mechanism to understand these experimental results. Starting from intermediate 2, the most favorable reaction sequence is a fast tautomerization, followed by a σ-migration of the base moiety, and a final inverse-electron demanding Diels-Alder reaction, resulting in the formation of the BA product stereoselectively. This reaction mechanism can also be applied to understand the biosynthesis of BA that involves enzymatic catalysis.

From unbound to bound states: Ab initio molecular dynamics of ammonia clusters with an excess electron

Ab initio molecular dynamics simulations of negatively charged clusters of 2-48 ammonia molecules were performed to elucidate the electronic stability of the excess electron as a function of cluster size. We show that while the electronic stability of finite temperature clusters increases with cluster size, as few as 5-7 ammonia molecules can bind an excess electron, reaching a vertical binding energy slightly less than half of the bulk value for the largest system studied. These results, which are in agreement with previous studies wherever available, allowed us to analyze the excess electron binding patterns in terms of its radius of gyration and shape anisotropy and provide a qualitative interpretation based on a particle-in-a-spherical-well model.

Assessment of the similarity-transformed equation of motion (STEOM) for open-shell organic and transition metal molecules

This study benchmarks the newly re-implemented single-reference excited-state methods, IP-EOM-CCSD, EA-EOM-CCSD, and STEOM-CCSD, in ORCA6.0, with a focus on open-shell systems. We compare STEOM against EOM-CCSD, CC3, and CCSDT across a range of systems, including small organic radicals, hydrated transition metal (TM) ions, and TM diatomic systems with both closed and open-shell configurations. For organic radicals, STEOM and EOM-CCSD show comparable performance, aligning closely with CC3 and CCSDT results. In the case of hydrated TM ions, IP-EOM closely matches DLPNO-CCSD results, while deviations from DLPNO-CCSD(T) are consistent. For open-shell TM systems, IP-EOM exhibits a blueshift relative to both the DLPNO-CCSD methods, while EA-EOM-CCSD shows better agreement. When comparing STEOM and CC3 to CCSDT, STEOM shows slightly larger deviations in closed-shell systems but shows excellent agreement in open-shell systems. Computational efficiency is also assessed, revealing a significant speedup in ORCA 6.0 compared to ORCA 5.0, with optimizations improving computation times. This study provides valuable insights into the performance and efficiency of STEOM in various chemical environments, highlighting its strengths and limitations.

Catalytic prenyl conjugate additions for synthesis of enantiomerically enriched PPAPs

Polycyclic polyprenylated acylphloroglucinols (PPAPs) are a class of >400 natural products with a broad spectrum of bioactivity, ranging from antidepressant and antimicrobial to anti-obesity and anticancer activity. Here, we present a scalable, regio-, site-, and enantioselective catalytic method for synthesis of cyclic β-prenyl ketones, compounds that can be used for efficient syntheses of many PPAPs in high enantiomeric purity. The transformation is prenyl conjugate addition to cyclic β-ketoesters promoted by a readily accessible chiral copper catalyst and involving an easy-to-prepare and isolable organoborate reagent. Reactions reach completion in just a few minutes at room temperature. The importance of this advance is highlighted by the enantioselective preparation of intermediates previously used to generate racemic PPAPs. We also present the enantioselective synthesis of nemorosonol (14 steps, 20% yield) and its one-step conversion to another PPAP, garcibracteatone (52% yield).

Exploring the potential of natural orbital functionals

In recent years, Natural Orbital Functional (NOF) theory has gained increasing significance in quantum chemistry, successfully addressing one of the field's most challenging problems: providing an accurate and balanced description of systems with strong electronic correlation. The quest for NOFs that strike the delicate balance between computational tractability and predictive accuracy represents a holy grail for researchers. Today, NOFs provide an alternative formalism to both density functional and wavefunction-based methods, with their appeal rooted in a wonderfully simple conceptual framework. This perspective outlines the basic concepts, strengths and weaknesses, and current status of NOFs, while offering suggestions for their future development.

Probing the structure and dynamics of the heterocyclic PAH xanthene and its water complexes with infrared and microwave spectroscopy

To assess the presence of oxygen-containing polycyclic aromatic hydrocarbons (OPAHs) in the interstellar medium and understand how water aggregates on an OPAH surface, we present a comprehensive gas-phase spectroscopy investigation of the OPAH xanthene (C13H10O) and its complexes with water using IR-UV ion dip spectroscopy and chirped-pulse Fourier transform microwave spectroscopy. The infrared spectrum of xanthene shows weak features at 3.42, 3.43, and 3.47 μm, which have been suggested to partly originate from vibrational modes of PAHs containing sp3 hybridized carbon atoms, in agreement with the molecular structure of xanthene. The high resolution of rotational spectroscopy reveals a tunneling splitting of the rotational transitions, which can be explained with an out-of-plane bending motion of the two lateral benzene rings of xanthene. The nature of the tunnelling motion is elucidated by observing a similar splitting pattern in the rotational transitions of the singly-substituted 13C isotopologues. The rotational spectroscopy investigation is extended to hydrates of xanthene with up to four water molecules. Different xanthene-water binding motifs are observed based on the degree of hydration, with O-H⋯π interactions becoming preferred over O-H⋯Oxanthene interactions as the degree of hydration increases. A structural comparison with water complexes of related molecular systems highlights the impact of the substrate's shape and chemical composition on the arrangement of the surrounding water molecules.

On the Potential Energy Surface of the Pyrene Dimer

Knowledge of reliable geometries and associated intermolecular interaction energy (ΔE) values at key fragments of the potential energy surface (PES) in the gas phase is indispensable for the modeling of various properties of the pyrene dimer (PYD) and other important aggregate systems of a comparatively large size (ca. 50 atoms). The performance of the domain-based local pair natural orbital (DLPNO) variant of the coupled-cluster theory with singles, doubles and perturbative triples in the complete basis set limit [CCSD(T)/CBS] method for highly accurate predictions of the ΔE at a variety of regions of the PES was established for a representative set of pi-stacked dimers, which also includes the PYD. For geometries with the distance between stacked monomers close to a value of such a distance in the ΔE minimum structure, an excellent agreement between the canonical CCSD(T)/CBS results and their DLPNO counterparts was found. This finding enabled us to accurately characterize the lowest-lying configurations of the PYD, and the physical origin of their stabilization was thoroughly analyzed. The proposed DLPNO-CCSD(T)/CBS procedure should be applied with the aim of safely locating a global minimum of the PES and firmly establishing the pertaining ΔE of even larger dimers in studies of packing motifs of organic electronic devices and other novel materials.

Overcoming Challenges in O-Nitration: Selective Alcohol Nitration Deploying N,6-Dinitrosaccharin and Lewis Acid Catalysis

Nitrate esters hold pivotal roles in pharmaceuticals, energetic materials, and atmospheric processes, motivating the development of efficient synthesis routes. Here, we present a novel catalytic method for the synthesis of nitrates via the direct O-nitration of alcohols, addressing limitations of current traditional methods. Leveraging bench-stable and recoverable N,6-dinitrosaccharin reagent, our catalytic strategy employs magnesium triflate to achieve mild and selective O-nitration of alcohols, offering broad substrate scope and unprecedentedly large functional group tolerance (e.g. alkenes, alkynes, carbonyls). DFT mechanistic studies reveal a dual role of the magnesium catalyst in the activation of both the nitrating reagent and the alcohol substrate. They also unveil a barrierless proton transfer upon formation of a widely-accepted - yet elusive in solution - nitrooxonium ion intermediate. Overall, our work contributes to the development of mild, selective, and sustainable approaches to nitrates synthesis, with potential applications in drug discovery, materials science, and environmental chemistry.

Insights into the Structural Modification of Selenium-Doped Derivatives with Narrowband Emissions: A Theory Study

The research on boron/nitrogen (B/N)-based multiresonance thermally activated delayed fluorescence (MR-TADF) emitters has been a prominent topic due to their narrowband emission and high luminous efficiency. However, devices derived from the common types of narrowband TADF materials often experience an efficiency roll-off, which could be ascribed to their relatively slow triplet-singlet exciton interconversion. Since inserting the heavy Se atom into the B/N scheme has been a proven strategy to address the abovementioned issues, herein, extensive density functional theory (DFT) and time-dependent DFT (TD-DFT) simulations have been employed to explore the effects of the structural modification on a series of structurally modified selenium-doped derivatives. Furthermore, the two-layered ONIOM (QM/MM) model has been employed to study the pressure effects on the crystal structure and photophysical properties of the pristine CzBSe. The theoretical results found that the introduced tert-butyl units in Cz-BSeN could result in a shorter charge transfer distance and smaller reorganization energy than the parent CzBSe. In contrast to directly incorporating the o-carborane (Cb) unit to CzBSe, incorporating the bridged phenyl units is important in order to achieve narrowband emissions and high luminous efficiency. The lowest three triplet excited states of CzBSe, Cz-BSeN and PhCb-BSeN all contribute to their triplet-singlet exciton conversions, resulting in a high utilization of triplet excitons. The pressure has an evident influence on the photophysical properties of the aggregated CzBSe and is favored for obtaining narrowband emissions. Our work is promised to provide a feasible strategy for designing selenium-doped derivatives with narrowband emissions and rapid triplet-singlet exciton interconversions.

Probing the alkylidene carbene-strained alkyne equilibrium in polycyclic systems via the Fritsch-Buttenberg-Wiechell rearrangement

Strained cycloalkynes are valuable building blocks in synthetic chemistry due to their high degree of reactivity and ability to form structurally complex scaffolds, common features of many pharmaceuticals and natural products. Alkylidene carbenes provide a pathway to the formation of strained cycloalkynes through Fritsch-Buttenberg-Wiechell rearrangements, but this strategy, like other methods of alkyne generation, is believed to depend upon a thermodynamic equilibrium that favors the alkyne over the carbene. Herein three highly strained, polycyclic alkynes, previously thought to be thermodynamically inaccessible, are generated under mild conditions and intercepted through Diels-Alder cycloaddition with a diene trapping agent. The use of a different trapping agent also allows for the interception of the alkylidene carbene, providing the first instance in which both an exocyclic alkylidene carbene and its cycloalkyne Fritsch-Buttenberg-Wiechell rearrangement product have been trapped.

Electric Field Dependence of EPR Hyperfine Coupling Constants

In this article, we introduce a simple but reliable method to calculate the electric field dependence of the isotropic hyperfine coupling tensor Aiso for free radicals. This dependence, also referred to as the Bloembergen effect, can be of interest in analyzing EPR experiments for solid-state materials but is rarely studied for isolated radicals in the gaseous phase. The proposed method uses the numerical differentiation of the field-perturbed A tensor and, consequently, as a purely numerical method, does not depend on a quantum chemical method used to determine the hyperfine tensor A. To test the performance of the proposed method, we used a set of 28 systems, including seven organic radicals, the majority of which were taken from a recent benchmark study by Bartlett's group. We employed the well-tested and robust density functional theory (DFT) functional, namely CAM-B3LYP, and the variant of the CCSD method based on local pair natural orbitals, namely DLPNO-CCSD to calculate the Aiso tensor itself, and all first derivatives or Bloembergen effect constants ai(1).

Anions Formed by Perchlorofluorobenzenes C6ClnF6-n (n = 0-6)

Anions formed by the perhalobenzene series C6ClnF6-n (n = 0-6) are studied computationally. All members of the series form both stable valence and stable nonvalence anions. At the geometry of the neutral parents, only nonvalence anions are bound, and the respective vertical electron affinities show values in the 20 to 60 meV range. Valence anions show distorted nonplanar structures, and one can distinguish two types of conformers. A-type conformers show puckered-ring structures and excess electrons delocalized over several C-Cl bonds [in the case of C6F6-, C-F bonds], while B-type conformers possess excess electrons essentially localized in a single C-Cl bond, which is accordingly strongly stretched and bent out of plane. For a specific anion, all conformers are close in energy (relative energies of less than 10 kJ/mol) and are connected by low-lying transition states. Accordingly, A-type and B-type conformers possess similar adiabatic electron affinities; however, their vertical detachment energies exhibit drastically different values, which should ease conformer distinction in photoelectron spectroscopy.

The Azide-Wittig Reaction

The introduction of heteroatoms into conjugated organic molecules is an important strategy to tune their reactivity and physical properties. In this realm triazabutadienes (TBDs) of the general from R2C=N-N=NR' are an interesting class of compounds, however, general synthetic protocols for their generation are limited. Based on the serendipitous finding that the sterically encumbered azide Mes*N3 (Mes*=2,4,6-tBu3C6H2) reacted with PMe3 in the presence of an aromatic aldehyde to form a TBD, we now report on the "Azide-Wittig" reaction. This azide-Wittig reaction is shown to be a versatile tool for the synthesis of a variety of TBDs, tolerating a wide range of aldehydes and organic azides as coupling partners. The preference for azide-Wittig, rather than aza-Wittig reactivity was rationalized using computational methods. This study shows how kinetic control can significantly alter the reaction pathway, thereby switching from an aza-Wittig to an azide-Wittig regime.

The Quest for Stable Borozene Core in Main-Group Capped Inverse Sandwich Complexes, [(HE)2B6H6]2- (E=B, Al, Ga, In, and Tl)

The ubiquitous chemistry of benzene led us to explore ways to stabilise analogous borozene, by capping them with appropriate groups. The mismatch in overlap of ring-cap fragment molecular orbitals in [(HB)2B6H6]2- is overcome by replacing the two BH caps with higher congeners of boron. We calculated the relative energies of all the polyhedral structural candidates for [(HE)2B6H6]2- (E=Al-Tl) and found hexagonal bipyramid (HBP) to be more stable with Al-H caps. A global minimum search also gives HBP as the most stable structure for [Al2B6H8]2-. The capped B6H6 ring in [(HAl)2B6H6]2- has aromaticity comparable to that of benzene.

Spontaneous Molecular Bromine Production in Sea-Salt Aerosols

Bromine chemistry is responsible for the catalytic ozone destruction in the atmosphere. The heterogeneous reactions of sea-salt aerosols are the main abiotic sources of reactive bromine in the atmosphere. Here, we present a novel mechanism for the activation of bromide ions (Br-) by O2 and H2O in the absence of additional oxidants. The laboratory and theoretical calculation results demonstrated that under dark conditions, Br-, O2 and H3O+ could spontaneously generate Br and HO2 radicals through a proton-electron transfer process at the air-water interface and in the liquid phase. Our results also showed that light and acidity could significantly promote the activation of Br- and the production of Br2. The estimated gaseous Br2 production rate was up to 1.55×1010 molecules cm-2 ⋅ s-1 under light and acidic conditions; these results showed a significant contribution to the atmospheric reactive bromine budget. The reactive oxygen species (ROS) generated during Br- activation could promote the multiphase oxidation of SO2 to produce sulfuric acid, while the increase in acidity had a positive feedback effect on Br- activation. Our findings highlight the crucial role of the proton-electron transfer process in Br2 production; here, H3O+ facilitates the activation of Br- by O2, serves as a significant source of atmospheric reactive bromine and exerts a profound impact on the atmospheric oxidation capacity.

A switchable green emitting dye and its phenomenal properties: implications for the photoluminescence features of carbon dots

New molecular fluorophores are constantly being discovered in post-synthetic mixtures of carbon dots (CDs), prompting researchers to elucidate their role in the optical properties of these nanomaterials. It has been reported that the green-emitting fluorophore that forms during the synthesis of popular citric acid/urea CDs is HPPT (4-hydroxy-1H-pyrrolo[3,4-c]pyridine-1,3,6(2H,5H)-trione). However, due to the low concentration of HPPT-like molecules within the structure of CDs, their actual binding and contribution to the optical properties of CDs has not so far been convincingly confirmed. In this joint experimental and quantum chemical study, we show that HPPT is a strong acid and only its anionic form, HPPT-, present in solutions of pH 0-10, is emissive. Next, its fluorescence can be switched off rapidly in strongly alkaline environments as a result of HPPT- hydrolysis, leading to the opening of its pyrrole ring and formation of CDPC (3-carbamoyl-2,6-dihydroxypyridine-4-carboxylic acid), existing as the CDPC2- dianion under these conditions. Eventually, we found that the ring opening hydrolysis is reversible and the green emission may be restored in acidic environments. The kinetics and mechanism of this hydrolysis were also revealed. The optical features of citric acid (CA)-urea CDs under various conditions were compared with a simpler CD system prepared by treating the CDs obtained from CA solely with HPPT- (HPPT@CDs). Our results indicate the feasibility of the post-synthetic modification of HPPT- present in the structures of CA-urea CDs and HPPT@CDs. Without HPPT- they emit blue fluorescence only. Thus, this makes the nanosystem switch the PL emission colour reversibly from green to blue owing to the opening and closing of the pyrrole ring in HPPT-like molecules. More importantly, the latter process may be considered a first step toward genuine fine tuning of the PL emission colour from CDs. These findings are of general importance to the further development of citric acid-based CDs with tailored properties.

The Counterintuitive Relationship between Orbital Energy, Orbital Overlap, and Bond Covalency in CeF62- and CeCl62

The 4f orbitals of Ce(IV) have shown appreciably enhanced covalent mixing with ligand orbitals relative to those of Ce(III). Here, X-ray spectroscopy, magnetic susceptibility measurements, and theoretical methods are used to investigate 4f covalency in CeF62- and CeCl62-. These techniques show covalent mixing between Ce 4f and F 2p orbitals to be about 25% less than mixing between Ce 4f and Cl 3p orbitals, placing CeF62- among the most ionic Ce(IV) compounds to-date. However, ligand field analysis using the experimental data shows significantly higher 4f orbital overlap with the F 2p orbitals compared to the Cl 3p. This result is counterintuitive since the Ce-F bonds display less 4f covalency despite their higher orbital overlap, and greater overlap is traditionally associated with enhanced bond covalency. The weaker covalency is attributed to the large energy gap between Ce 4f and F 2p orbitals strongly counteracting the higher orbital overlap. These results highlight that only a concerted consideration of both atomic orbital overlap and energy matching in f-element systems leads to an accurate picture of their bonding.

Bidentate Substrate Binding Mode in Oxalate Decarboxylase

Oxalate decarboxylase is an Mn- and O2-dependent enzyme in the bicupin superfamily that catalyzes the redox-neutral disproportionation of the oxalate monoanion to form carbon dioxide and formate. Its best-studied isozyme is from Bacillus subtilis where it is stress-induced under low pH conditions. Current mechanistic schemes assume a monodentate binding mode of the substrate to the N-terminal active site Mn ion to make space for a presumed O2 molecule, despite the fact that oxalate generally prefers to bind bidentate to Mn. We report on X-band 13C-electron nuclear double resonance (ENDOR) experiments on 13C-labeled oxalate bound to the active-site Mn(II) in wild-type oxalate decarboxylase at high pH, the catalytically impaired W96F mutant enzyme at low pH, and Mn(II) in aqueous solution. The ENDOR spectra of these samples are practically identical, which shows that the substrate binds bidentate (κO, κO') to the active site Mn(II) ion. Domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) calculations of the expected 13C hyperfine coupling constants for bidentate bound oxalate predict ENDOR spectra in good agreement with the experiment, supporting bidentate bound substrate. Geometry optimization of a substrate-bound minimal active site model by density functional theory shows two possible substrate coordination geometries, bidentate and monodentate. The bidentate structure is energetically preferred by ~4.7 kcal/mol. Our results revise a long-standing hypothesis regarding substrate binding in the enzyme and suggest that dioxygen does not bind to the active site Mn ion after substrate binds. The results are in agreement with our recent mechanistic hypothesis of substrate activation via a long-range electron transfer process involving the C-terminal Mn ion.

In silico investigation on the separation of disulfide bonds by N-heterocyclic carbene

Herein, the separation of a disulfide bond using different nucleophilic agents like tri-methyl phosphine (TMP), tris (2-carboxyethyl) phosphine (TCEP), and N-heterocyclic carbene (NHC) has been investigated. Both TMP and TCEP have demonstrated their ability to break disulfide bonds through the SN2 mechanism. However, it is worth noting that these reactions are endothermic. While searching for a suitable nucleophile, it was observed that the NHC-mediated reaction was exothermic. The natural bond orbital (NBO), principal interacting orbital (PIO) and extended transition state-natural orbitals for chemical valence (ETS-NOCV) studies help understand the electron transfer process between interacting orbitals during the chemical reactions.

Accessing the Cloke-Wilson Rearrangement via Conjugate Addition of Phosphoranes to Michael Acceptors: A Route to Cyclopropanes and 5-Membered Ring Heterocycles Investigated by Density Functional and Ab Initio Theory

Conjugate addition of unstabilized Wittig-type phosphonium ylides to 1,1-diacceptor- and 1-acceptor-substituted alkenes is investigated by density functional theory and high-level ab initio (DLPNO-CCSD(T)) calculations. The results indicate that the initial conjugate addition step should be facile with barriers predicted to be between 0 and 21 kcal mol-1. Potential intramolecular follow-up reactions include the formation of acceptor-substituted cyclopropanes as well as the formation of dihydrofuran derivatives via intramolecular SN2-type transition state structures. The barriers calculated for these potentially valuable cyclization reactions are substantial with Gibbs free energies of activation between 19 and 40 kcal mol-1. Competing reaction channels include Wittig olefination (for ketones and aldehydes), as well as Claisen condensation reactions. The reaction offers an alternative entry point to the nucleophile-catalyzed Cloke-Wilson rearrangement.

Range-separated density functional theory using multiresolution analysis and quantum computing

Quantum computers are expected to outperform classical computers for specific problems in quantum chemistry. Such calculations remain expensive, but costs can be lowered through the partition of the molecular system. In the present study, partition was achieved with range-separated density functional theory (RS-DFT). The use of RS-DFT reduces both the basis set size and the active space size dependence of the ground state energy in comparison with the use of wave function theory (WFT) alone. The utilization of pair natural orbitals (PNOs) in place of canonical molecular orbitals (MOs) results in more compact qubit Hamiltonians. To test this strategy, a basis-set independent framework, known as multiresolution analysis (MRA), was employed to generate PNOs. Tests were conducted with the variational quantum eigensolver for a number of molecules. The results show that the proposed approach reduces the number of qubits needed to reach a target energy accuracy.

Prediction of Covalent Metal-Metal Bonding in Cp-M-M'-Nacnac Complexes of Group 2 and 12 Metals (Be, Mg, Ca, Zn, Cd, Hg)

Bimetallic CpMM'Nacnac molecules with group 2 and 12 metals (M=Be, Mg, Ca, Zn, Cd, Hg) that contain novel metal-metal bonding have been investigated in a theoretical study of their molecular and electronic structure, thermodynamic stability, and metal-metal bonding. In all cases the metal-metal bonds are characterized as electron-sharing covalent single bonds from natural bond orbital (NBO) and energy-decomposition analysis with natural orbitals of chemical valence (EDA-NOCV) analysis. The sum of [MM'] charges is relatively constant, with all complexes exhibiting a [MM']2+ core. Quantum theory of atoms in molecules (QTAIM) analysis indicates the presence of non-nuclear attractors (NNA) in the metal-metal bonds of the BeBe, MgMg, and CaCa complexes. There is substantial electron density (0.75-1.33 e) associated with the NNAs, which indicates that these metal-metal bonds, while classified as covalent electron-sharing bonds, retain significant metallic character that can be associated with reducing reactivity of the complex. The predicted stability of these complexes, combined with their novel covalent metal-metal bonding and potential as reducing agents, make them appealing targets for the synthesis of new metal-metal bonds.

Rapid aqueous-phase dark reaction of phenols with nitrosonium ions: Novel mechanism for atmospheric nitrosation and nitration at low pH

Dark aqueous-phase reactions involving the nitrosation and nitration of aromatic organic compounds play a significant role in the production of light-absorbing organic carbon in the atmosphere. This process constitutes a crucial aspect of tropospheric chemistry and has attracted growing research interest, particularly in understanding the mechanisms governing nighttime reactions between phenols and nitrogen oxides. In this study, we present new findings concerning the rapid dark reactions between phenols containing electron-donating groups and inorganic nitrite in acidic aqueous solutions with pH levels <3.5. This reaction generates a substantial amount of nitroso- and nitro-substituted phenolic compounds, known for their light-absorbing properties and toxicity. In experiments utilizing various substituted phenols, we demonstrate that their reaction rates with nitrite depend on the electron cloud density of the benzene ring, indicative of an electrophilic substitution reaction mechanism. Control experiments and theoretical calculations indicate that the nitrosonium ion (NO+) is the reactive nitrogen species responsible for undergoing electrophilic reactions with phenolate anions, leading to the formation of nitroso-substituted phenolic compounds. These compounds then undergo partial oxidation to form nitro-substituted phenols through reactions with nitrous acid (HONO) or other oxidants like oxygen. Our findings unveil a novel mechanism for swift atmospheric nitrosation and nitration reactions that occur within acidic cloud droplets or aerosol water, providing valuable insights into the rapid nocturnal formation of nitrogen-containing organic compounds with significant implications for climate dynamics and human health.

Tensor hypercontraction for fully self-consistent imaginary-time GF2 and GWSOX methods: Theory, implementation, and role of the Green's function second-order exchange for intermolecular interactions

We present an efficient MPI-parallel algorithm and its implementation for evaluating the self-consistent correlated second-order exchange term (SOX), which is employed as a correction to the fully self-consistent GW scheme called scGWSOX (GW plus the SOX term iterated to achieve full Green's function self-consistency). Due to the application of the tensor hypercontraction (THC) in our computational procedure, the scaling of the evaluation of scGWSOX is reduced from O(nτnAO5) to O(nτN2nAO2). This fully MPI-parallel and THC-adapted approach enabled us to conduct the largest fully self-consistent scGWSOX calculations with over 1100 atomic orbitals with only negligible errors attributed to THC fitting. Utilizing our THC implementation for scGW, scGF2, and scGWSOX, we evaluated energies of intermolecular interactions. This approach allowed us to circumvent issues related to reference dependence and ambiguity in energy evaluation, which are common challenges in non-self-consistent calculations. We demonstrate that scGW exhibits a slight overbinding tendency for large systems, contrary to the underbinding observed with non-self-consistent RPA. Conversely, scGWSOX exhibits a slight underbinding tendency for such systems. This behavior is both physical and systematic and is caused by exclusion-principle violating diagrams or corresponding corrections. Our analysis elucidates the role played by these different diagrams, which is crucial for the construction of rigorous, accurate, and systematic methods. Finally, we explicitly show that all perturbative fully self-consistent Green's function methods are size-extensive and size-consistent.