Intrinsic Atomic Orbitals: An Unbiased Bridge between Quantum Theory and Chemical Concepts

Facile synthesis and bonding of 4-ferrocenyl-1,2,4-triazol-5-ylidene complexes

Ferrocene-substituted carbenes have emerged as attractive, redox-active ligands. However, among the compounds studied to date, ferrocenylated 1,2,4-triazol-5-ylidenes, which are closely related to the archetypal imidazol-2-ylidenes, are still unknown. Here, we demonstrate that the triazolium salt [CHN(Me)NCHN(Fc)]I (2; Fc = ferrocenyl), obtained by alkylation of 4-ferrocenyl-4H-1,2,4-triazole (1) with MeI, reacts selectively with metal alkoxide/hydroxide precursors [(cod)Rh(OMe)]2 and [(IPr)Au(OH)] (cod = cycloocta-1,5-diene, IPr = 1,3-bis(2,6-diisopropylphenyl)imidazol-2-ylidene) to produce the ferrocene-substituted 1,2,4-triazol-5-ylidene complexes [(cod)RhI{CN(Me)NCHN(Fc)}] and [(IPr)Au{CN(Me)NCHN(Fc)}]I in good yields. The complexes were characterised by NMR and IR spectroscopy, mass spectrometry, cyclic voltammetry, and single-crystal X-ray diffraction analysis. Density function theory (DFT) calculations were used to rationalise the electrochemical behaviour of the carbene complexes and to elucidate the bonding situation in these compounds. An analysis using intrinsic bond orbitals (IBOs) revealed that the 1,2,4-triazol-5-ylidene ligand exerted a strong trans influence and showed a synergistic stabilisation by the negative inductive and positive π-donor effects of the nitrogen atoms adjacent to the carbene carbon atom; these effects were enhanced by conjugation with the CHN bond at the exterior, similar to that in imidazol-2-ylidenes.

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Ultrasound-assisted diastereoselective green synthesis of spiro-fused-γ-lactams functionalized with an amide bond heterocyclic bioisostere via the Ugi azide/domino process coupled strategy

A series of linked-type 1,5-disubstituted tetrazoles (1,5-DS-Ts) were synthesised via an isocyanide-based multicomponent reaction (IMCR) and used as synthetic platforms to access bound-type polyheterocycles containing an epoxyisoindol-1(6H)-one scaffold under green conditions. This rapid sonochemical synthetic strategy includes a double domino process using an orthogonal heterocyclic input in the Ugi-azide (UA) reaction. DFT calculations and NBO analysis were performed to understand the pseudopericyclic reaction involved in the 1,5-electrocyclization of the UA mechanism.

Efficiency and Mechanism of Catalytic Siloxane Exchange in Vitrimer Polymers: Modeling and Density Functional Theory Investigations

Of late, siloxane-containing vitrimers have gained significant interest due to their fast dynamic characteristics over a reasonable temperature range (180-220 °C), making them well-suited for diverse applications. The exchange reaction pathway in the siloxane vitrimers is accountable for the covalent adaptive network, with the reaction's effectiveness being regulated by either organic or organometallic catalysts. However, directly studying the exchange reaction pathway in the bulk phase using experimental approaches is challenging because of the intricate and interconnected structure of these vitrimers. Here, we perform comprehensive density functional theory (DFT) and experimental investigations to discover the detailed catalytic efficacy of siloxane exchange and provide direction for the reaction process using a 1,5,7-triazabicyclo[4.4.0]dec-5-ene (TBD) catalyst. The calculated transition barrier energy and catalytic efficiency of hexamethyldisiloxane and dihydroxy-dimethylsilane exchange derived from the nudged elastic band with transition-state calculations strongly agree with the experimental findings. In addition, Fukui indices, along with partial charges, are employed to evaluate the nucleophilic and electrophilic behaviors of silanol and siloxane molecules. Our analysis revealed that by utilizing the Fukui indices of both the acid and the base, we can make an approximate estimation of the respective kinetics of the SN2 process in the siloxane exchange reaction mechanism. These findings establish a foundation for comprehending a crucial aspect of the exchange mechanism in siloxane vitrimer systems and could aid in the development of novel catalysts.

P(V)-bis(amidophenolate) ligand cooperation: stoichiometric CO-bond cleavage in aldehydes and ketones

The cooperation between a geometrically constrained, highly electrophilic phosphorus(V) center, and an electronically rich tetradentate bis(amidophenolate) ligand enables the cleavage of the CO bond from typical aldehydes and ketones delivering iminio phosphoramidate species. The amphiphilic nature of these products, which is demonstrated through their reaction with typical Lewis acids and bases, enables their use as a mild source of silylium cations from silanes, allowing the selective reductive coupling of aldehydes to ethers under catalytic conditions.

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

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

Rabbit muscle pyruvate kinase activators: Synthesis, molecular docking and theoretical studies of N-substituted sulfonamide derivatives

Pyruvate kinase (PK) activators have potential therapeutic applications in diseases such as sickle cell anemia. In this study, N-Substituted sulfonamide derivatives of 1,4-dihydropyridines were synthesized and evaluated as PK activators in vitro and using molecular docking studies. The compounds were synthesized by reacting dicarbonyl compounds with ammonium acetate, 5-nitrobenzaldehyde, and alumina sulfuric acid (ASA), followed by reduction and sulfonylation. The structures of the compounds were analyzed using spectroscopic techniques. DFT calculations provided insights into the electronic properties. Molecular docking of the compounds into the active site of PK showed favorable binding interactions. ADME evaluation indicated suitable solubility, BBB permeation, and lack of CYP450 inhibition. Overall, this study demonstrates the potential of new hybrid 1,4-dihydropyridine substituted sulfonamides as PK activators for further development. According to AC50 values, the compound (DTSF7, 0.97μM) is about 100-fold higher affective than the clinically used sulfonamide compound (AC50 = 90μM) for PK.

Desulfurdioxidative N-N Coupling of N-Arylhydroxylamines and N-Sulfinylanilines: Reaction Development and Mechanism

A highly efficient and chemoselective approach for the divergent assembling of unsymmetrical hydrazines through an unprecedented intermolecular desulfurdioxidative N-N coupling is developed. This metal free protocol employs readily accessible N-arylhydroxylamines and N-sulfinylanilines to provide highly valuable hydrazine products with good reaction yields and excellent functional group tolerance under simple conditions. Computational studies suggest that the in situ generated O-sulfenylated arylhydroxylamine intermediate undergoes a retro-[2π+2σ] cycloaddition via a stepwise diradical mechanism to form the N-N bond and release SO2.

Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes

Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me2-1.1.1, and 36Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3)2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6-2.0 eV) and hence these complexes belong to the superalkali category. The ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3)2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.

N-(Sulfonio)Sulfilimine Reagents: Non-Oxidizing Sources of Electrophilic Nitrogen Atom for Skeletal Editing

The one-pot synthesis of λ4-dibenzothiophen-5-imino-N-dibenzothiophenium triflate (1) in multigram scale is reported. This compound reacts with Rh2(esp)2 (esp=α,α,α',α'-tetramethyl-1,3-benzenedipropionic acid) generating a Rh-coordinated sulfonitrene species, which is able to transfer the electrophilic nitrene moiety to olefins. When indenes are used as substrates, isoquinolines are obtained in good yields. We assumed that after formation of the corresponding N-sulfonio aziridine, a ring expansion occurs via selective C-C bond cleavage and concomitant elimination of dibenzothiophene. Unexpectedly, a similar protocol transforms 1-arylcyclobutenes into 1-cyano-1-arylcyclopropanes. Our calculations indicate that aziridination is not favored in this case; instead, sulfilimine-substituted cyclobutyl carbocations are initially formed, and these evolve to the isolated cyclopropanes via ring contraction. Both procedures are operationally simple, tolerate a range of functional groups, including oxidation-sensitive alcohols and aldehydes, and enable the convenient preparation of valuable 15N-labelled products. These results demonstrate the potential of 1 to provide alternative pathways for the selective transfer of N-atoms in organic molecules.

A crystalline stannyne

The synthesis of heteronuclear alkyne analogues incorporating heavier group 14 elements (R1-C≡E-R2, E = Si, Ge, Sn, Pb) has posed a long-standing challenge. Neutral silynes (R1-C≡Si(L)-R2) and germynes (R1-C≡Ge(L)-R2) stabilized by a Lewis base have achieved sufficient stability for structural characterization at low temperatures. Here we show the isolation of a base-free stannyne (R1-C≡Sn-R2) at room temperature, achieved through the strategic use of a bulky cyclic phosphino ligand in combination with a bulky terphenyl substituent. Despite an allenic structure with strong delocalization of π-electrons, this compound exhibits adjacent ambiphilic carbon and tin centres, forming a carbon-tin multiple bond with ionic character. The stannyne demonstrates reactivity similar to carbenes or stannylenes, reacting with 1-adamantyl isocyanide and 2,3-dimethyl-1,3-butadiene. Additionally, its carbon-tin bond can be saturated by Et3N·HCl or cleaved by isopropyl isocyanate.

Silyliumylidene Ion Stabilized by Two σ-Donating Ni(0)- and Pd(0)-Fragments

A silyliumylidene ion 2 stabilized by two σ-donating Ni(0)- and Pd(0)-fragments was successfully synthesized. Due to the σ-donation of M→Si interactions, 2 presents a pyramidalized cationic silicon center with a localized lone pair. The additional coordination of basic Pd(0) fragment to the mono-Ni(0)-stabilized silyliumylidene 1 results in a higher HOMO level and an unchanged HOMO-LUMO gap and thus, 2 remains highly reactive. Interestingly, the coordination mode at the Si center is closely related to the nature of M-ligands. Indeed, the donor/donor-stabilized silyliumylidene ion 2 has been transformed into a donor/acceptor-stabilized ion 13, featuring a trigonal planar Si center with a vacant orbital, just via a ligand exchange reaction from PCy3/NHC toward PMe3.

Electron Dynamics in Alkane C-H Activation Mediated by Transition Metal Complexes

Alkanes, ideal raw materials for industrial chemical production, typically exhibit limited reactivity due to their robust and weakly polarized C-H bonds. The challenge lies in selectively activating these C-H bonds under mild conditions. To address this challenge, various C-H activation mechanisms have been developed. Yet, classifying these mechanisms depends on the overall stoichiometry, which can be ambiguous and sometimes problematic. In this study, we utilized density functional theory calculations combined with intrinsic bond orbital (IBO) analysis to examine electron flow in the four primary alkane C-H activation mechanisms: oxidative addition, σ-bond metathesis, 1,2-addition, and electrophilic activation. Methane was selected as the representative alkane molecule to undergo C-H heterolytic cleavage in these reactions. Across all mechanisms studied, we find that the CH3 moiety in methane consistently uses an electron pair from the cleaved C-H bond to form a σ-bond with the metal. Yet, the electron pair that accepts the proton differs with each mechanism: in oxidative addition, it is derived from the d-orbitals; in σ-bond metathesis, it resulted from the metal-ligand σ-bonds; in 1,2-addition, it arose from the π-orbital of the metal-ligand multiple bonds; and in electrophilic activation, it came from the lone pairs on ligands. This detailed analysis not only provides a clear visual understanding of these reactions but also showcases the ability of the IBO method to differentiate between mechanisms. The electron flow discerned from IBO analysis is further corroborated by results from absolutely localized molecular orbital energy decomposition analysis, which also helps to quantify the two predominant interactions in each process. Our findings offer profound insights into the electron dynamics at play in alkane C-H activation, enhancing our understanding of these critical reactions.

Light Induced Cobalt(III) Carbene Radical Formation from Dimethyl Malonate As Carbene Precursor

Radical-type carbene transfer catalysis is an efficient method for the direct functionalization of C-H and C=C bonds. However, carbene radical complexes are currently formed via high-energy carbene precursors, such as diazo compounds or iodonium ylides. Many of these carbene precursors require additional synthetic steps, have an explosive nature, or generate halogenated waste. Consequently, the utilization of carbene radical catalysis is limited by specific carbene precursors that access the carbene radical intermediate. In this study, we generate a cobalt(III) carbene radical complex from dimethyl malonate, which is commercially available and bench-stable. EPR and NMR spectroscopy were used to identify the intermediates and showed that the cobalt(III) carbene radical complex is formed upon light irradiation. In the presence of styrene, carbene transfer occurred, forming cyclopropane as the product. With this photochemical method, we demonstrate that dimethyl malonate can be used as an alternative carbene precursor in the formation of a cobalt(III) carbene radical complex.

Propagated (fragment) Pipek-Mezey Wannier functions in real-time time-dependent density functional theory

Localization procedures are an important tool for analysis of complex systems in quantum chemistry, since canonical molecular orbitals are delocalized and can, therefore, be difficult to align with chemical intuition and obscure information at the local level of the system. This especially applies to calculations obeying periodic boundary conditions. The most commonly used approach to localization is Foster-Boys Wannier functions, which use a unitary transformation to jointly minimize the second moment of the orbitals. This procedure has proven to be robust and fast but has a side effect of often mixing σ- and π-type orbitals. σ/π-separation is achieved by the Pipek-Mezey Wannier function (PMWF) approach [Lehtola and Jónsson, J. Chem. Theory Comput. 10, 642 (2014) and Jónsson et al., J. Chem. Theory Comput. 13, 460 (2017)], which defines the spread functional in terms of partial charges instead. We have implemented a PMWF algorithm in the CP2K software package using the Cardoso-Souloumiac algorithm to enable their application to real-time time-dependent density functional theory. The method is demonstrated on stacked CO2 molecules, linear acetylenic carbon, boron and nitrogen co-doped graphene, and nitrogen-vacancy doped diamond. Finally, we discuss its computational scaling and recent efforts to improve it with fragment approaches.

Proton-Coupled Electron Transfer and Hydrogen Tunneling in Olive Oil Phenol Reactions

Olive oil phenols are recognized as molecules with numerous positive health effects, many of which rely on their antioxidative activity, i.e., the ability to transfer hydrogen to radicals. Proton-coupled electron transfer reactions and hydrogen tunneling are ubiquitous in biological systems. Reactions of olive oil phenols, hydroxytyrosol, tyrosol, oleuropein, oleacein, oleocanthal, homovanillyl alcohol, vanillin, and a few phenolic acids with a DPPH• (2,2-diphenyl-1-picrylhydrazyl) radical in a 1,4-dioxane:water = 95:5 or 99:1 v/v solvent mixture were studied through an experimental kinetic analysis and computational chemistry calculations. The highest rate constants corresponding to the highest antioxidative activity are obtained for the ortho-diphenols hydroxytyrosol, oleuropein, and oleacein. The experimentally determined kinetic isotope effects (KIEs) for hydroxytyrosol, homovanillyl alcohol, and caffeic acid reactions are 16.0, 15.4, and 16.7, respectively. Based on these KIEs, thermodynamic activation parameters, and an intrinsic bond orbital (IBO) analysis along the IRC path calculations, we propose a proton-coupled electron transfer mechanism. The average local ionization energy and electron donor Fukui function obtained for the phenolic compounds show that the most reactive electron-donating sites are associated with π electrons above and below the aromatic ring, in support of the IBO analysis and proposed PCET reaction mechanism. Large KIEs and isotopic values of Arrhenius pre-exponential factor AH/AD determined for the hydroxytyrosol, homovanillyl alcohol, and caffeic acid reactions of 0.6, 1.3, and 0.3, respectively, reveal the involvement of hydrogen tunneling in the process.

Isomer-specific photofragmentation of C3H3+ at the carbon K-edge

Individual fingerprints of different isomers of C3H3+ cations have been identified by studying photoionization, photoexcitation, and photofragmentation of C3H3+ near the carbon K-edge. The experiment was performed employing the photon-ion merged-beams technique at the photon-ion spectrometer at PETRA III (PIPE). This technique is a variant of near-edge X-ray absorption fine-structure spectroscopy, which is particularly sensitive to the 1s → π* excitation. The C3H3+ primary ions were generated by an electron cyclotron resonance ion source. C3Hn2+ product ions with n = 0, 1, 2, and 3 were observed for photon energies in the range of 279.0 eV to 295.2 eV. The experimental spectra are interpreted with the aid of theoretical calculations within the framework of time-dependent density functional theory. To this end, absorption spectra have been calculated for three different constitutional isomers of C3H3+. We find that our experimental approach offers a new possibility to study at the same time details of the electronic structure and of the geometry of molecular ions such as C3H3+.

Molecular-strain induced phosphinidene reactivity of a phosphanorcaradiene

Phosphanorcaradienes are an appealing class of phosphorus compounds that can serve as synthons of transient phosphinidenes. However, the synthesis of such species is a formidable task owing to their intrinsic high reactivity. Herein we report straightforward synthesis, characterization and reactivity studies of a phosphanorcaradiene, in which one of the benzene rings in the flanking fluorenyl substituents is intramolecularly dearomatized through attachment to the phosphorus atom. It is facilely obtained by the reduction of phosphorus(III) dichloride precursor with potassium graphite. Despite being thermally robust, it acts as a synthetic equivalent of a transient phosphinidene. It reacts with trimethylphosphine and isonitrile to yield phosphanylidene-phosphorane and 1-phospha-3-azaallene, respectively. When it is treated with one and two molar equivalents of azide, iminophosphane and bis(imino)phosphane are isolated, respectively. Moreover, it is capable of activating ethylene and alkyne to afford [1 + 2] cycloaddition products, as well as oxidative cleavage of Si-H and N-H bonds to yield secondary phosphines. All the reactions proceed smoothly at room temperature without the presence of transition metals. The driving force for these reactions is most likely the high ring-constraint of the three-membered PC2 ring and recovery of the aromaticity of the benzene ring.

Assessment of Embedding Schemes in a Hybrid Machine Learning/Classical Potentials (ML/MM) Approach

Machine learning (ML) methods have reached high accuracy levels for the prediction of in vacuo molecular properties. However, the simulation of large systems solely through ML methods (such as those based on neural network potentials) is still a challenge. In this context, one of the most promising frameworks for integrating ML schemes in the simulation of complex molecular systems are the so-called ML/MM methods. These multiscale approaches combine ML methods with classical force fields (MM), in the same spirit as the successful hybrid quantum mechanics-molecular mechanics methods (QM/MM). The key issue for such ML/MM methods is an adequate description of the coupling between the region of the system described by ML and the region described at the MM level. In the context of QM/MM schemes, the main ingredient of the interaction is electrostatic, and the state of the art is the so-called electrostatic-embedding. In this study, we analyze the quality of simpler mechanical embedding-based approaches, specifically focusing on their application within a ML/MM framework utilizing atomic partial charges derived in vacuo. Taking as reference electrostatic embedding calculations performed at a QM(DFT)/MM level, we explore different atomic charges schemes, as well as a polarization correction computed using atomic polarizabilites. Our benchmark data set comprises a set of about 80k small organic structures from the ANI-1x and ANI-2x databases, solvated in water. The results suggest that the minimal basis iterative stockholder (MBIS) atomic charges yield the best agreement with the reference coupling energy. Remarkable enhancements are achieved by including a simple polarization correction.

DFT Struggles to Predict the Energy Landscape for Iron Pyridine Diimine-Catalyzed [2 + 2] Cycloaddition of Alkenes: Insights into the Problem and Alternative Solutions

In recent years, noninnocent pyridine diimine (PDI) complexes featuring first-row transition metals have emerged as prominent catalysts, demonstrating efficacy in a diverse range of vital organometallic transformations. However, the inherent complexity of the fundamental reactivity paradigm in these systems arises from the presence of a noninnocent ligand and the multispin feasibility of 3d metals. While density functional theory (DFT) has been widely used to unravel mechanistic insights, its limitations as a single-reference method can potentially misrepresent spin-state energetics, compromising our understanding of these intricate systems. In this study, we employ extensive high-level ab initio state averaged-complete active space self-consistent field/N-electron valence state perturbation theory (SA-CASSCF/NEVPT2) calculations in combination with DFT to investigate an iron-PDI-catalyzed [2 + 2] cycloaddition reaction of alkenes. The transformation proceeds through two major steps: oxidative cyclization and reductive elimination. Contrary to the predictions of DFT calculations, which suggest two-state reactivity in the reaction and identify reductive elimination as the turnover-limiting step, SA-CASSCF/NEVPT2-corrected results unequivocally establish a single-state reactivity scenario with oxidative cyclization as the turnover-limiting step. SA-CASSCF/NEVPT2-based insights into electronic ground states and electron distribution elucidate the intriguing interactions between the PDI ligand and the iron center, revealing the highly multiconfigurational nature of these species and providing a precise depiction of metal-ligand cooperativity throughout the transformation. A comparative assessment of several widely recognized DFT functionals against SA-CASSCF/NEVPT2-corrected data indicates that single-point energy calculations using the modern density functional MN15 on TPSSh geometries offer the most reliable density functional methodology, in scenarios where SA-CASSCF/NEVPT2 computational cost is a consideration.

Carbene-Stabilized Phosphagermylenylidene: A Heavier Analog of Isonitrile

Phosphagermylenylidenes (R-P═Ge), as heavier analogs of isonitriles, whether in their free state or as complexes with a Lewis base, have not been previously identified as isolable entities. In this study, we report the synthesis of a stable monomeric phosphagermylenylidene within the coordination sphere of a Lewis base under ambient conditions. This species was synthesized by Lewis base-induced dedimerization of a cyclic phosphagermylenylidene dimer or via Me3SiCl elimination from a phosphinochlorogermylene framework. The deliberate integration of a bulky, electropositive N-heterocyclic boryl group at the phosphorus site, combined with coordination stabilization by a cyclic (alkyl)(amino)carbene at the low-valent germanium site, effectively mitigated its natural tendency toward oligomerization. Structural analyses and theoretical calculations have demonstrated that this unprecedented species features a P═Ge double bond, characterized by conventional electron-sharing π and σ bonds, complemented by lone pairs at both the phosphorus and germanium atoms. Preliminary reactivity studies show that this base-stabilized phosphagermylenylidene demonstrates facile release of ligands at the Ge atom, coordination to silver through the lone pair on P, and versatile reactivity including both (cyclo)addition and cleavage of the P═Ge double bond.

Bispericyclic Ambimodal Dimerization of Pentafulvene: The Origin of Asynchronicity and Kinetic Selectivity of the Endo Transition State

The propensity of fulvenes to undergo dimerization has long been known, although the in-depth mechanism and electronic behavior during dimerization are still elusive. Herein, we made an attempt to gain insights into the reactivity of pentafulvene for Diels-Alder (DA) and [6 + 4]-cycloadditions via conventional and ambimodal routes. The result emphasizes that pentafulvene dimerization preferentially proceeds through a unique bifurcation mechanism where two DA pathways merge together to produce two degenerate [4 + 2]-cycloadducts from a single TS. Despite the [6 + 4]-cycloadduct being thermodynamically preferred, [4 + 2]-cycloaddition reactions are kinetically driven. Singlet biradicaloid is involved in through-space 6e- delocalization as a secondary orbital interaction that originates asynchronicity and stabilizes the bispericyclic transition state (TS). The transformation of various actively participating intrinsic bonding orbitals (IBOs) unambiguously forecasts the formation of multiple products from a single TS and rationalizes the mechanism of ambimodal reactions that are rather difficult to probe with other analyses. The changes in active IBOs clearly distinguish the conventional reactions from bifurcation reactions and can be employed to characterize and confirm the ambimodal mechanism. This report gains a crucial theoretical insight into the mechanism of bifurcation, the origin of asynchronicity, and electronic behavior in ambimodal TS, which will certainly be of enormous value for future studies.

On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment

Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.

Projected Atomic Orbitals As Optimal Virtual Space for Excited State Projection-Based Embedding Calculations

The projected atomic orbital (PAO) technique is presented for the construction of virtual orbital spaces in projection-based embedding (PbE) applications. The proposed straightforward procedure produces a set of virtual orbitals that are used in the final, high-level calculation of the embedded active subsystem. The PAO scheme is demonstrated on intermolecular potentials of bimolecular complexes in ground and excited states, including Rydberg excitations. The results show the outstanding performance of the PbE method when used with PAO virtual orbitals compared with those produced using common orbital localization techniques. The good agreement of the resulting PbE potential curves with those from high-level ab initio dimer calculations, also in diffuse basis sets, confirms that the PAO technique can be suggested for future applications using top-down embedding methods.

The Cyanoketenyl Anion [NC3O]. Angew Chem Int Ed Engl 2024;63:e202403766. [PMID: 38470943 DOI: 10.1002/anie.202403766] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Cumulenes and heterocumulenes with three or more cumulative multiple bonds are usually reactive species that serve as valuable building blocks for more complex molecules but tend to isomerize or cyclize and therefore are difficult to isolate. Using a mild ligand exchange reaction at the carbon in α-metalated ylides, we have now succeeded in the synthesis and gram-scale isolation of the elusive cyanoketenyl anion [NC3O]-. Despite its assumed cumulene-like structure and the delocalization of the negative charge across the whole 5-atom molecule, it features a bent geometry with a nucleophilic central carbon atom. Computational studies reveal an ambiguous bonding situation in the anion, which can be illustrated only by a combination of different resonance structures. Nonetheless, the anion features remarkable stability, thus allowing the storage of its potassium-crown ether salt and its application as a highly functional synthetic building block. The cyanoketenyl anion readily reacts with a series of small molecules to form more complex organic compounds, including industrially valuable compounds such as cyanoacetate. This work demonstrated that reactive species can be generated by novel synthesis methods and open up atom-economic pathways to complex compounds from small abundant molecules.

Computational Organic Chemistry: The Frontier for Understanding and Designing Bioorthogonal Cycloadditions

Computational organic chemistry has become a valuable tool in the field of bioorthogonal chemistry, offering insights and aiding in the progression of this branch of chemistry. In this review, I present an overview of computational work in this field, including an exploration of both the primary computational analysis methods used and their application in the main areas of bioorthogonal chemistry: (3 + 2) and [4 + 2] cycloadditions. In the context of (3 + 2) cycloadditions, detailed studies of electronic effects have informed the evolution of cycloalkyne/1,3-dipole cycloadditions. Through computational techniques, researchers have found ways to adjust the electronic structure via hyperconjugation to enhance reactions without compromising stability. For [4 + 2] cycloadditions, methods such as distortion/interaction analysis and energy decomposition analysis have been beneficial, leading to the development of bioorthogonal reactants with improved reactivity and the creation of orthogonal reaction pairs. To conclude, I touch upon the emerging fields of cheminformatics and machine learning, which promise to play a role in future reaction discovery and optimization.

The p-block challenge: assessing quantum chemistry methods for inorganic heterocycle dimerizations

The elements of the p-block of the periodic table are of high interest in various chemical and technical applications like frustrated Lewis-pairs (FLP) or opto-electronics. However, high-quality benchmark data to assess approximate density functional theory (DFT) for their theoretical description are sparse. In this work, we present a benchmark set of 604 dimerization energies of 302 "inorganic benzenes" composed of all non-carbon p-block elements of main groups III to VI up to polonium. This so-called IHD302 test set comprises two classes of structures formed by covalent bonding and by weaker donor-acceptor (WDA) interactions, respectively. Generating reliable reference data with ab initio methods is challenging due to large electron correlation contributions, core-valence correlation effects, and especially the slow basis set convergence. To compute reference values for these dimerization reactions, after thorough testing, we applied a computational protocol using state-of-the-art explicitly correlated local coupled cluster theory termed PNO-LCCSD(T)-F12/cc-VTZ-PP-F12(corr.). It includes a basis set correction at the PNO-LMP2-F12/aug-cc-pwCVTZ level. Based on these reference data, we assess 26 DFT methods in combination with three different dispersion corrections and the def2-QZVPP basis set, five composite DFT approaches, and five semi-empirical quantum mechanical methods. For the covalent dimerizations, the r2SCAN-D4 meta-GGA, the r2SCAN0-D4 and ωB97M-V hybrids, and the revDSD-PBEP86-D4 double-hybrid functional are found to be the best-performing methods among the evaluated functionals of the respective class. However, since def2 basis sets for the 4th period are not associated to relativistic pseudo-potentials, we obtained significant errors in the covalent dimerization energies (up to 6 kcal mol-1) for molecules containing p-block elements of the 4th period. Significant improvements were achieved for systems containing 4th row elements by using ECP10MDF pseudopotentials along with re-contracted aug-cc-pVQZ-PP-KS basis sets introduced in this work with the contraction coefficients taken from atomic DFT (PBE0) calculations. Overall, the IHD302 set represents a challenge to contemporary quantum chemical methods. This is due to a large number of spatially close p-element bonds which are underrepresented in other benchmark sets, and the partial covalent bonding character for the WDA interactions. The IHD302 set may be helpful to develop more robust and transferable approximate quantum chemical methods in the future.

Toward a Full Configurational Accuracy Calculation of an Arbitrary Molecule via Fragment Embedding and a Stochastic Solver

We demonstrate the feasibility of using a stochastic solver, full configuration interaction quantum Monte Carlo with the initiator approximation (i-FCIQMC), to converge fragment embedding calculations, namely bootstrap embedding (BE). We first propose and test a general protocol for converging BE-i-FCIQMC calculations and then suggest how the quality of the calculation compares against that of deterministic BE-FCI using different numbers of walkers. We then demonstrate that BE-i-FCIQMC can perform as well as BE-FCI in the large walker limit and how different factors, including the size of the Hilbert space of the fragments, the number of walkers, and the nature of the chemical system, affect the achievable matching error. We finally perform BE-FCI calculations in realistic systems like benzene and cyclohexane using a double-ζ basis set. This work demonstrates the potential of performing FCI quality calculations in realistic systems using BE.

The Separability Problem in Molecular Quantum Systems: Information-Theoretic Framework for Atoms in Molecules

Even though molecules are fundamentally quantum entities, the concept of a molecule retains certain classical attributes concerning its constituents. This includes the empirical separability of a molecule into its three-dimensional, rigid structure in Euclidean space, a framework often obtained through experimental methods like X-Ray crystallography. In this work, we delve into the mathematical implications of partitioning a molecule into its constituent parts using the widely recognized Atoms-In-Molecules (AIM) schemes, aiming to establish their validity within the framework of Information Theory concepts. We have uncovered information-theoretical justifications for employing some of the most prevalent AIM schemes in the field of Chemistry, including Hirshfeld (stockholder partitioning), Bader's (topological dissection), and the quantum approach (Hilbert's space definition). In the first approach we have applied the generalized principle of minimum relative entropy derived from the Sharma-Mittal two-parameter functional, avoiding the need for an arbitrary selection of reference promolecular atoms. Within the ambit of topological-information partitioning, we have demonstrated that the Fisher information of Bader's atoms conform to a comprehensive theory based on the Principle of Extreme Physical Information avoiding the need of employing the Schwinger's principle, which has been proven to be problematic. For the quantum approach we have presented information-theoretic justifications for conducting Löwdin symmetric transformations on the density matrix to form atomic Hilbert spaces generating orthonormal atomic orbitals with maximum occupancy for a given wavefunction.

Modelling on a Biomimetic [Cu-O-Cu]2+-mediated Methane-to-Methanol Conversion Unveils the Site for Methane Activation

The Cu-O-Cu core exhibits methane-to-methanol conversion, mirroring the reactivity of the copper-containing enzyme pMMO. Herein, we computationally examined the reactivity of a biomimetic Cu-O-Cu core towards methane-to-methanol conversion. The oxygen atom of the Cu-O-Cu core abstracts hydrogen present in the C-H bond of methane. The spin density at the bridging oxygen helps to abstract hydrogen from the C-H bond. We modulated the spin density of the bridging oxygen by substituting only a single copper atom of the Cu-O-Cu core by metals (M) such as Fe, Co, and Ag. These substitutions result in bimetallic [Cu-O-M]2+ models. We observed that the energy barriers for the C-H activation step and the subsequent rebound step vary with the metal M. [Cu-O-Ag]2+ exhibits the highest reactivity for M2M conversion, while [Cu-O-Fe]2+ displays the lowest reactivity. To understand the different reactivity of these models towards M2M conversion, we employed distortion-interaction analysis, orbital analysis, spin density analysis, and quantum theory of atoms in molecules analysis. Orbital analysis reveals that all four adducts follow a hydrogen atom transfer mechanism for C-H activation. Further, spin density analysis reveals that a higher spin density on the bridging oxygen leads to a lower C-H activation barrier.

Ground and excited electronic structures of electride and alkalide units: The cases of Metal-Tren, -Azacryptand, and -TriPip222 complexes

A systematic electronic structure analysis was conducted for M(L)n molecular electrides and their corresponding alkalide units M(L)n @M' (M/M' = Na, K; L = Tren, Azacryptand, TriPip222; n = 1, 2). All complexes belong to the "superalkali" category due to their low ionization potentials. The saturated molecular electrides display M+ (L)n - form with a greatly diffuse quasispherical electron cloud. They were identified as "superatoms" considering the contours of populating atomic-type molecular orbitals. The observed superatomic Aufbau order of M(Tren)2 is 1S, 1P, 1D, 1F, 2S, 2P, and 1G and it is consistent with those of M(Azacryptand) and M(TriPip222) up to the analyzed 1F level. Their excitation energies decrease gradually moving from M(Tren)2 to M(Azacryptand) and to M(TriPip222). The studied alkalide complexes carry [M(L)n ]+ @M'- ionic structure and their dissociation energies vary in the sequence of K(L)n @Na > Na(L)n @Na > K(L)n @K > Na(L)n @K. Similar to molecular electrides, the anions of alkalide units occupy electrons in diffuse Rydberg-like orbitals. In this work, excited states of [M(L)n @M']0/+/- and their trends are also analyzed.

Ligand-Enabled Oxidative Fluorination of Gold(I) and Light-Induced Aryl-F Coupling at Gold(III)

MeDalphos Au(I) complexes featuring aryl, alkynyl, and alkyl groups readily react with electrophilic fluorinating reagents such as N-fluorobenzenesulfonimide and Selectfluor. The ensuing [(MeDalphos)Au(R)F]+ complexes have been isolated and characterized by multinuclear NMR spectroscopy as well as X-ray diffraction. They adopt a square-planar contra-thermodynamic structure, with F trans to N. DFT/IBO calculations show that the N lone pair of MeDalphos assists and directs the transfer of F+ to gold. The [(MeDalphos)Au(Ar)F]+ (Ar = Mes, 2,6-F2Ph) complexes smoothly engage in C-C cross-coupling with PhCCSiMe3 and Me3SiCN, providing direct evidence for the oxidative fluorination/transmetalation/reductive elimination sequence proposed for F+-promoted gold-catalyzed transformations. Moreover, direct reductive elimination to forge a C-F bond at Au(III) was explored and substantiated. Thermal means proved unsuccessful, leading mostly to decomposition, but irradiation with UV-visible light enabled efficient promotion of aryl-F coupling (up to 90% yield). The light-induced reductive elimination proceeds under mild conditions; it works even with the electron-deprived 2,6-difluorophenyl group, and it is not limited to the contra-thermodynamic form of the aryl Au(III) fluoride complexes.

Chemical Aristocracy: He3 Dication and Analogous Noble-Gas-Exclusive Covalent Compounds

Herein, we predict the first set of covalently bonded triatomic molecular compounds composed exclusively of noble gases. Using a combination of double-hybrid DFT, CCSD(T), and MRCI+Q calculations and a range of bonding analyses, we explored a set of 270 doubly charged triatomics, which included various combinations of noble gases and main group elements. This extensive exploration uncovered nine noble-gas-exclusive covalent compounds incorporating helium, neon, argon, or combinations thereof, exemplified by cases such as He32+ and related systems. This work brings to light a previously uncharted domain of noble gas chemistry, demonstrating the potential of noble gases in forming covalent molecular clusters.

Chemical Bonding and the Role of Node-Induced Electron Confinement

The chemical bond is the cornerstone of chemistry, providing a conceptual framework to understand and predict the behavior of molecules in complex systems. However, the fundamental origin of chemical bonding remains controversial and has been responsible for fierce debate over the past century. Here, we present a unified theory of bonding, using a separation of electron delocalization effects from orbital relaxation to identify three mechanisms [node-induced confinement (typically associated with Pauli repulsion, though more general), orbital contraction, and polarization] that each modulate kinetic energy during bond formation. Through analysis of a series of archetypal bonds, we show that an exquisite balance of energy-lowering delocalizing and localizing effects are dictated simply by atomic electron configurations, nodal structure, and electronegativities. The utility of this unified bonding theory is demonstrated by its application to explain observed trends in bond strengths throughout the periodic table, including main group and transition metal elements.

Dynamically or Kinetically Controlled? Computational Study of the Mechanisms of Electrophilic Aminoalkenylation of Heteroaromatics with Keteniminium Ions

Quantum chemical calculations and molecular dynamics simulations were applied to study the electrophilic aminoalkenylation of heteroaromatics with keniminium ions. Post-transition state bifurcation (PTSB) was found in the electrophilic addition step for the aminoalkenylation of pyrroles and indoles, and the selectivity for these reactions was dynamically controlled. However, the aminoalkenylation of furan was kinetically controlled because no apparent PTSB was found in the electrophilic addition step. The substituents on the keteniminium ions can also affect the dynamic results for the aminoalkenylations to pyrroles: the C2-aminoalkenylated product is much more favored over the C3-aminoalkenylated product for keteniminium ions with electron-donating substituents, while the product ratio (C2 product/C3 product) decreased when stronger electron-withdrawing substituents were applied.

Valence Tautomerism Induced Proton Coupled Electron Transfer:X-H Bond Oxidation with a Dinuclear Au(II) Hydroxide Complex

We report the preparation and characterization of the dinuclear AuII hydroxide complex AuII 2(L)2(OH)2 (L=N,N'-bis (2,6-dimethyl) phenylformamidinate) and study its reactivity towards weak X-H bonds. Through the interplay of kinetic analysis and computational studies, we demonstrate that the oxidation of cyclohexadiene follows a concerted proton-coupled electron transfer (cPCET) mechanism, a rare type of reactivity for Au complexes. We find that the Au-Au σ-bond undergoes polarization in the PCET event leading to an adjustment of oxidation levels for both Au centers prior to C(sp3)-H bond cleavage. We thus describe the oxidation event as a valence tautomerism-induced PCET where the basicity of one reduced Au-OH unit provides a proton acceptor and the second more oxidized Au center serves as an electron acceptor. The coordination of these events allows for unprecedented radical-type reactivity by a closed shell AuII complex.

Deciphering the irradiation induced fragmentation-rearrangement mechanisms in valence ionized CF3CH2F

The increasing presence of 1,1,1,2-tetrafluoroethane (CF3CH2F) in the atmosphere has prompted detailed studies into its complex photodissociation behavior. Experiments focusing on CF3CH2F irradiation have unveiled an array of ions, with the persistent observation of the rearrangement product CHF2+ not yet fully understood. In this work, we combine density functional theory, coupled-cluster calculations with a complete basis set formalism, and atom-centered density matrix propagation molecular dynamics to investigate the energetics and dynamics of different potential pathways leading to CHF2+. We found that the two-body dissociation pathway involving an HF rearrangement, which was previously considered complex for CHF2+ formation, is actually straightforward but not likely due to the facile loss of HF. In contrast, our calculations reveal that the H elimination pathway, once thought of as a potential route to CHF2+, is not only comparably disadvantageous from both thermodynamic and kinetic points of view but also does not align with experimental data, particularly the lack of a rebound peak at m/z 101-102. We establish that the formation of CHF2+ is predominantly via the HF elimination channel, a conclusion experimentally corroborated by studies involving the trifluoroethylene cation CF2CHF+, a key intermediate in this process.

Cobalt-Catalyzed Carbon-Heteroatom Transfer Enables Regioselective Tricomponent 1,4-Carboamination

Tricomponent cobalt(salen)-catalyzed carbofunctionalization of unsaturated substrates by radical-polar crossover has the potential to streamline access to broad classes of heteroatom-functionalized synthetic targets, yet the reaction platform has remained elusive, despite the well-developed analogous hydrofunctionalizations mediated by high-valent alkylcobalt intermediates. We report herein the development of a cobalt(salen) catalytic system that enables carbofunctionalization. The reaction entails a tricomponent decarboxylative 1,4-carboamination of dienes and provides a direct route to aromatic allylic amines by obviating preformed allylation reagents and protection of oxidation-sensitive aromatic amines. The catalytic system merges acridine photocatalysis with cobalt(salen)-catalyzed regioselective 1,4-carbofunctionalization that facilitates the crossover of the radical and polar phases of the tricomponent coupling process, revealing critical roles of the reactants, as well as ligand effects and the nature of the formal high-valent alkylcobalt species on the chemo- and regioselectivity.

Triplet Excitation-Energy Transfer Couplings from Subsystem Time-Dependent Density-Functional Theory

We present an implementation of triplet excitation-energy transfer (TEET) couplings based on subsystem-based time-dependent density-functional theory (sTDDFT). TEET couplings are systematically investigated by comparing "exact" and approximate variants of sTDDFT. We demonstrate that, while sTDDFT utilizing explicit approximate non-additive kinetic energy (NAKE) density functionals is well-suited for describing singlet EET processes, it is inadequate for characterizing TEET. However, we show that projection-based embedding (PbE)-based sTDDFT addresses the challenges faced by NAKE-sTDDFT and emerges as a promising method for accurately describing electronic couplings in TEET processes. We also introduce the mixed PbE-/NAKE-embedding procedure to investigate the TEET effects in solvated pairs of chromophores. This approach offers a good balance between accuracy and efficiency, enabling comprehensive studies of TEET processes in complex environments.

Hydrosilylation of BB triple bonds: catalyst- and reductant-free construction of B-Si bonds and B2Si heterocycles

Hydrosilanes undergo mild, uncatalyzed single and double 1,2-addition across the B-B triple bonds of diborynes, leading to an unsymmetrical silyldiborene and compounds with novel non-cluster three-membered B2Si rings. The reactions are a new addition to the very few catalyst- and alkali-metal-free methods available for the construction of B-Si bonds.

Unveiling Three Interconvertible Redox States of Boraphenalene

Neutral 1-boraphenalene displays the isoelectronic structure of the phenalenyl carbocation and is expected to behave as an attractive organoboron multi-redox system. However, the isolation of new redox states have remained elusive even though the preparation of neutral boron(III)-containing phenalene compounds have been extensively studied. Herein, we have adopted an N-heterocyclic carbene ligand stabilization approach to achieve the first isolation of the stable and ambipolar 1-boraphenalenyl radical 1•. The 1-boraphenalenyl cation 1+ and anion 1- have also been electrochemically observed and chemically isolated, representing new redox forms of boraphenalene for the study of non-Kekulé polynuclear benzenoid molecules. Experimental and theoretical investigations suggest that the interconvertible three-redox-state species undergo reversible electronic structure modifications, which primarily take place on the polycyclic framework of the molecules, exhibiting atypical behavior compared to known donor-stabilized organoboron compounds. Initial reactivity studies, aromaticity evaluations, and photophysical studies show redox-state-dependent trends. While 1+ is luminescent in both the solution and solid states, 1• exhibits boron-centered reactivity and 1- undergoes substitution chemistry on the boraphenalenyl skeleton and serves as a single-electron transfer reductant.

Air- and photo-stable luminescent carbodicarbene-azaboraacenium ions

Substitution of a C=C bond by an isoelectronic B-N bond is a well-established strategy to alter the electronic structure and stability of acenes. BN-substituted acenes that possess narrow energy gaps have attractive optoelectronic properties. However, they are susceptible to air and/or light. Here we present the design, synthesis and molecular structures of fully π-conjugated cationic BN-doped acenes stabilized by carbodicarbene ligands. They are luminescent in the solution and solid states and show high air and moisture stability. Compared with their neutral BN-substituted counterparts as well as the parent all-carbon acenes, these species display improved quantum yields and small optical gaps. The electronic structures of the azabora-anthracene and azabora-tetracene cations resemble higher-order acenes while possessing high photo-oxidative resistance. Investigations using density functional theory suggest that the stability and photo-physics of these conjugated systems may be ascribed to their cationic nature and the electronic properties of the carbodicarbene.

[Br4F21]- - a unique molecular tetrahedral interhalogen ion containing a μ4-bridging fluorine atom surrounded by BrF5 molecules

The reaction of [NMe4][BrF6] with an excess of BrF5 leads to the compound [NMe4][Br4F21]·BrF5. It features molecular [(μ4-F)(BrF5)4]- anions of tetrahedron-like shape containing central μ4-bridging F atoms coordinated by four BrF5 molecules. It is the most BrF5-rich fluoridobromate anion by mass. Quantum-chemical calculations showed that the μ4-F-Br bonds within the anion are essentially ionic in nature. The compound is the first example where F atoms bridge μ4-like neither to metal nor to hydrogen atoms. It was characterized by Raman spectroscopy and by single-crystal X-ray diffraction. The latter showed surprisingly that its crystal structure is related to the intermetallic half-Heusler compound and structure type MgAgAs.

Single electron reduction of NHC-CO2-borane compounds

The carbon dioxide radical anion [CO2˙-] is a highly reactive species of fundamental and synthetic interest. However, the direct one-electron reduction of CO2 to generate [CO2˙-] occurs at very negative reduction potentials, which is often a limiting factor for applications. Here, we show that NHC-CO2-BR3 species - generated from the Frustrated Lewis Pair (FLP)-type activation of CO2 by N-heterocyclic carbenes (NHCs) and boranes (BR3) - undergo single electron reduction at a less negative potential than free CO2. A net gain of more than one volt was notably measured with a CAAC-CO2-B(C6F5)3 adduct, which was chemically reduced to afford [CAAC-CO2-B(C6F5)3˙-]. This room temperature stable radical anion was characterized by EPR spectroscopy and by single-crystal X-ray diffraction analysis. Of particular interest, DFT calculations showed that, thanks to the electron withdrawing properties of the Lewis acid, significant unpaired spin density is localised on the carbon atom of the CO2 moiety. Finally, these species were shown to exhibit analogous reactivity to the carbon dioxide radical anion [CO2˙-] toward DMPO. This work demonstrates the advantage provided by FLP systems in the generation and stabilization of [CO2˙-]-like species.

Phosphinoyl-Substituted Ketenyl Anions: Synthesis and Substituent Effects on the Structural Properties

Ketenyl anions are versatile intermediates in synthetic chemistry and have recently become accessible as isolable reagents from metalated ylides by exchange of the phosphine with CO. Herein, we report on a systematic study of substituent effects on the structure and bonding situation in ketenyl anions. A series of phosphinoyl-substituted ketenyl anions {[R2P(X)CCO]- with X = O, NTol, S, Se} were prepared by carbonylation of the corresponding yldiides and isolated as their corresponding potassium salts. NMR and IR spectroscopic analyses together with computational studies demonstrate that the more electron-withdrawing oxo- and iminophosphinoyl substituents increase the s-character in the bond to the ketene moiety and hence the ynolate character of the anion. This trend is particularly seen in solution, whereas the solid-state properties are influenced by packing effects affecting the bonding situation.

A Linear Two-Coordinate Cr(II) Complex: Synthesis, Characterization, and Reactivity

The synthesis and characterization of a linear two-coordinate Cr(II) amido complex, Cr{N(t Bu)Dipp}2 (Dipp=2,6-diisopropylphenyl), from the reaction of 1 molar equivalent (equiv) of CrCl2 and 2 equiv. of LiN(t Bu)Dipp is reported. Single-crystal X-ray diffractometry (SC-XRD) analysis revealed that it has a short Cr-N bond distance of 1.8878(9) Å, which could be attributed to the relatively less bulky nature of the amido ligand compared with reported systems. Furthermore, the oxidation reaction of the two-coordinate Cr(II) complex was explored. The oxidation reaction of Cr{N(t Bu)Dipp}2 with the one-electron oxidants AgOTf and [FeCp2 ][BArF 4 ] (BArF 4 - =[B{C6 H3 -3,5-(CF3 )2 }4 ]- ) afforded the trigonal planar three- and bent two-coordinate Cr(III) complexes Cr{N(t Bu)Dipp}2 (OTf) and [Cr{N(t Bu)Dipp}2 ][BArF 4 ], respectively. The reaction of Cr{N(t Bu)Dipp}2 with 1 equiv. of the organic azides AdN3 (Ad=1-adamantyl) and PhN3 afforded the three-coordinate Cr(IV) imido complexes Cr{N(t Bu)Dipp}2 (NAd) and Cr{N(t Bu)Dipp}2 (NPh), respectively. The reaction of Cr{N(t Bu)Dipp}2 and two equiv. of Me3 NO afforded the Cr(VI) dioxo complex Cr{N(t Bu)Dipp}2 (O)2 . The reaction of Cr{N(t Bu)Dipp}2 with 1 equiv. of CyN=C=NCy resulted in the insertion of the carbodiimide into the Cr-N bond, with the formation of a three-coordinate Cr(II) complex. Finally, density functional theory (DFT) calculations were used to elucidate the electronic structure of these complexes.

Rigorous Screened Interactions for Realistic Correlated Electron Systems

We derive a widely applicable first-principles approach for determining two-body, static effective interactions for low-energy Hamiltonians with quantitative accuracy. The algebraic construction rigorously conserves all instantaneous two-point correlation functions in a chosen model space at the level of the random phase approximation, improving upon the traditional uncontrolled static approximations. Applied to screened interactions within a quantum embedding framework, we demonstrate these faithfully describe the relaxation of local subspaces via downfolding high-energy physics in molecular systems, as well as enabling a systematically improvable description of the long-range plasmonic contributions in extended graphene.

SPAHM(a,b): Encoding the Density Information from Guess Hamiltonian in Quantum Machine Learning Representations

Recently, we introduced a class of molecular representations for kernel-based regression methods─the spectrum of approximated Hamiltonian matrices (SPAHM)─that takes advantage of lightweight one-electron Hamiltonians traditionally used as a self-consistent field initial guess. The original SPAHM variant is built from occupied-orbital energies (i.e., eigenvalues) and naturally contains all of the information about nuclear charges, atomic positions, and symmetry requirements. Its advantages were demonstrated on data sets featuring a wide variation of charge and spin, for which traditional structure-based representations commonly fail. SPAHM(a,b), as introduced here, expand the eigenvalue SPAHM into local and transferable representations. They rely upon one-electron density matrices to build fingerprints from atomic and bond density overlap contributions inspired from preceding state-of-the-art representations. The performance and efficiency of SPAHM(a,b) is assessed on the predictions for data sets of prototypical organic molecules (QM7) of different charges and azoheteroarene dyes in an excited state. Overall, both SPAHM(a) and SPAHM(b) outperform state-of-the-art representations on difficult prediction tasks such as the atomic properties of charged open-shell species and of π-conjugated systems.

Four-Coordinate Co(III) Imide with an Unusually Tilted Terminal Imido Ligand

We report herein the synthesis and characterization of a terminal Co(III) imido complex supported by an intermediate field N,N,C heteroscorpionate. This chemistry is enabled through the development of an additional member of this ligand type featuring Ph2(CH3)C- substituents, one of which weakly binds and stabilizes Co in the corresponding Co(I) precursor. The Co(III) imide is low-spin with no evidence for thermal population of open-shell excited states. Unusually, the imido ligand in this molecule tilts markedly toward the Calkyl donor. DFT calculations suggest this structural feature to be largely a result of strong Co-C covalency, underscoring the importance of M-C bonding in determining the (electronic) structure of metal centers supported by this class of ligand.

Transferability of atomic energies from alchemical decomposition

We study alchemical atomic energy partitioning as a method to estimate atomization energies from atomic contributions, which are defined in physically rigorous and general ways through the use of the uniform electron gas as a joint reference. We analyze quantitatively the relation between atomic energies and their local environment using a dataset of 1325 organic molecules. The atomic energies are transferable across various molecules, enabling the prediction of atomization energies with a mean absolute error of 23 kcal/mol, comparable to simple statistical estimates but potentially more robust given their grounding in the physics-based decomposition scheme. A comparative analysis with other decomposition methods highlights its sensitivity to electrostatic variations, underlining its potential as a representation of the environment as well as in studying processes like diffusion in solids characterized by significant electrostatic shifts.