Kohn-Sham Density in a Slater Orbital Basis Set.

Finite, atom-centered Slater basis sets are used to determine approximate Kohn-Sham molecular orbitals. This is achieved by minimizing the kinetic energy plus the sum-squared difference between the Kohn-Sham density and the full configuration interaction density. As a result of the finite basis, a weight factor is introduced to balance the two minimization components. Results herein show that this can be done systematically, without sensitive dependence on the choice of scaling factor. In addition, the algorithm is applied to the LiH diatomic for fractional electron counts, where stretching the bond introduces significant reorganization of the electron density. The analysis will show the correct KS orbital structure and reveal the effects of correlation and electron locality on the KS solutions.

Green Light Promoted Iridium(III)/Copper(I)-Catalyzed Addition of Alkynes to Aziridinoquinoxalines Through the Intermediacy of Azomethine Ylides.

This manuscript describes the development of alkyne addition to the aziridine moiety of aziridinoquinoxalines using dual Ir(III)/Cu(I) catalytic system under green light-emitting diode (LED) photolysis (λmax =525 nm). This mild method features high levels of chemo- and regioselectivity and was used to generate 30 highly functionalized substituted dihydroquinoxalines in 36-98 % yield. This transformation was also carried asymmetrically using phthalazinamine-based chiral ligand to provide 9 chiral addition products in 96 : 4 to 86 : 14 e.r. The experimental and quantum chemical explorations of this reaction suggest a mechanism that involves Ir(III)-catalyzed triplet energy transfer followed by a ring-opening reaction ultimately leading to the formation of azomethine ylide intermediates. These azomethine intermediates undergo sequential protonation/copper(I) acetylide addition to provide the products.

Oxidative Cross Dehydrogenative Coupling of N-Heterocycles with Aldehydes through C(sp3)-H Functionalization.

Existing methodologies for metal-catalyzed cross-couplings typically rely on preinstallation of reactive functional groups on both reaction partners. In contrast, C-H functionalization approaches offer promise in simplification of the requisite substrates; however, challenges from low reactivity and similar reactivity of various C-H bonds introduce considerable complexity. Herein, the oxidative cross dehydrogenative coupling of α-amino C(sp3)-H bonds and aldehydes to produce ketone derivatives is described using an unusual reaction medium that incorporates the simultaneous use of di-tert-butyl peroxide as an oxidant and zinc metal as a reductant. The method proceeds with a broad substrate scope, representing an attractive approach for accessing α-amino ketones through the formal acylation of C-H bonds α to nitrogen in N-heterocycles. A combination of experimental investigation and computational modeling provides evidence for a mechanistic pathway involving cross-selective nickel-mediated cross-coupling of α-amino radicals and acyl radicals.

Controlling Catalyst Behavior in Lewis Acid-Catalyzed Carbonyl-Olefin Metathesis.

Lewis acid-catalyzed carbonyl-olefin metathesis has introduced a new means for revealing the behavior of Lewis acids. In particular, this reaction has led to the observation of new solution behaviors for FeCl3 that may qualitatively change how we think of Lewis acid activation. For example, catalytic metathesis reactions operate in the presence of superstoichiometric amounts of carbonyl, resulting in the formation of highly ligated (octahedral) iron geometries. These structures display reduced activity, decreasing catalyst turnover. As a result, it is necessary to steer the Fe-center away from inhibiting pathways to improve the reaction efficiency and augment yields for recalcitrant substrates. Herein, we examine the impact of the addition of TMSCl to FeCl3-catalyzed carbonyl-olefin metathesis, specifically for substrates that are prone to byproduct inhibition. Through kinetic, spectroscopic, and colligative experiments, significant deviations from the baseline metathesis reactivity are observed, including mitigation of byproduct inhibition as well as an increase in the reaction rate. Quantum chemical simulations are used to explain how TMSCl induces a change in catalyst structure that leads to these kinetic differences. Collectively, these data are consistent with the formation of a silylium catalyst, which induces the reaction through carbonyl binding. The FeCl3 activation of Si-Cl bonds to give the silylium active species is expected to have significant utility in enacting carbonyl-based transformations.

Formal Cross-Coupling of Amines and Carboxylic Acids to Form sp3-sp2 Carbon-Carbon Bonds.

Amines and carboxylic acids are abundant synthetic building blocks that are classically united to form an amide bond. To access new pockets of chemical space, we are interested in the development of amine-acid coupling reactions that complement the amide coupling. In particular, the formation of carbon-carbon bonds by formal deamination and decarboxylation would be an impactful addition to the synthesis toolbox. Here, we report a formal cross-coupling of alkyl amines and aryl carboxylic acids to form C(sp3)-C(sp2) bonds following preactivation of the amine-acid building blocks as a pyridinium salt and N-acyl-glutarimide, respectively. Under nickel-catalyzed reductive cross-coupling conditions, a diversity of simple and complex substrates are united in good to excellent yield, and numerous pharmaceuticals are successfully diversified. High-throughput experimentation was leveraged in the development of the reaction and the discovery of performance-enhancing additives such as phthalimide, RuCl3, and GaCl3. Mechanistic investigations suggest phthalimide may play a role in stabilizing productive Ni complexes rather than being involved in oxidative addition of the N-acyl-imide and that RuCl3 supports the decarbonylation event, thereby improving reaction selectivity.

Machine Learning Strategies for Reaction Development: Toward the Low-Data Limit.

Machine learning models are increasingly being utilized to predict outcomes of organic chemical reactions. A large amount of reaction data is used to train these models, which is in stark contrast to how expert chemists discover and develop new reactions by leveraging information from a small number of relevant transformations. Transfer learning and active learning are two strategies that can operate in low-data situations, which may help fill this gap and promote the use of machine learning for tackling real-world challenges in organic synthesis. This Perspective introduces active and transfer learning and connects these to potential opportunities and directions for further research, especially in the area of prospective development of chemical transformations.

Advances in Parallel Heat Bath Configuration Interaction.

Heat-bath configuration interaction (HCI) is a deterministic method that approaches the full CI limit at greatly reduced computational cost. In this work, computational improvements to the HCI algorithm are introduced targeting speed, parallel efficiency, and memory requirements. The new implementation introduces a hash function to distribute determinants and takes advantage of MPI and OpenMP for parallelism allowing for a (22e,168o) active space to be studied, which explicitly includes 2.39 × 107 variational determinants and 8.95 × 1010 perturbative determinants. Benchmarks show up to 86% parallel efficiency of the perturbative step on 32 nodes (4096 cores) and a total efficiency of 74%. The new HCI implementation is benchmarked for accuracy against prior results and applied to study the triplet-quintet gap in the challenging [FeO(NH3)5]2+ complex.

Exchange correlation potentials from full configuration interaction in a Slater orbital basis.

Ryabinkin-Kohut-Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. Examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

The numerical evaluation of Slater integrals on graphics processing units.

This article presents SlaterGPU, a graphics processing unit (GPU) accelerated library that uses OpenACC to numerically compute Slater-type orbital (STO) integrals. The electron repulsion integrals (ERI) are computed under the RI approximation using the Coulomb potential of the Slater basis function. To fully realize the performance capabilities of modern GPUs, the Slater integrals are evaluated in mixed-precision, resulting in speedups for the ERIs of over 80×. Parallelization on multiple GPUs allows for integral throughput of over 3 million integrals per second. This places STO integral throughput within reach of single-threaded, conventional Gaussian integration schemes. To test the quality of the integrals, the fluorine exchange reaction barrier in fluoromethane was computed using heat-bath configuration interaction (HBCI). In addition, the singlet-triplet gap of cyclobutadiene was examined using HBCI in a triple- ζ , polarized basis set. These benchmarks demonstrate the library's ability to generate the full set of integrals necessary for configuration interaction with up to 6 h functions in the auxiliary basis.

Conformer-RL: A deep reinforcement learning library for conformer generation.

Conformer-RL is an open-source Python package for applying deep reinforcement learning (RL) to the task of generating a diverse set of low-energy conformations for a single molecule. The library features a simple interface to train a deep RL conformer generation model on any covalently bonded molecule or polymer, including most drug-like molecules. Under the hood, it implements state-of-the-art RL algorithms and graph neural network architectures tuned specifically for molecular structures. Conformer-RL is also a platform for researching new algorithms and neural network architectures for conformer generation, as the library contains modular class interfaces for RL environments and agents, allowing users to easily swap components with their own implementations. Additionally, it comes with tools to visualize and save generated conformers for further analysis. Conformer-RL is well-tested and thoroughly documented with tutorials for each of the functionalities mentioned above, and is available on PyPi and Github: https://github.com/ZimmermanGroup/conformer-rl.

Glycosyl Exchange of Unactivated Glycosidic Bonds: Suppressing or Embracing Side Reactivity in Catalytic Glycosylations.

While developing boron-catalyzed glycosylations using glycosyl fluoride donors and trialkylsilyl ether acceptors, competing pathways involving productive glycosylation or glycosyl exchange were observed. Experimental and computational mechanistic studies suggest a novel mode of reactivity where a dioxolenium ion is a key intermediate that promotes both pathways through addition to either a silyl ether or to the acetal of an existing glycosidic linkage. Modifications in catalyst structure enable either pathway to be favored, and with this understanding, improved multicomponent iterative couplings and glycosyl exchange processes were demonstrated.

Interplay Between Applied Force and Radical Attack in the Mechanochemical Chain Scission of Poly(acrylic acid).

Sonication and radical attack are both known to contribute to breaking down polymers. Quantum chemical models show how the two can operate together, where radical attack is shown to reduce the effective tensile strength of the material. Using poly(acrylic acid) (PAA) as a model, hydrogen atom abstraction in PAA was found to improve the thermodynamics and kinetics of bond scission. The force needed for bond rupture was estimated to decrease from 4.7 to 2.5 nN. This occurs because hydrogen atom abstraction drastically alters the potential energy surface of the scissile bond. Bond activation was also found to decrease the magnitude of the changes in bond scission geometries and energetics in response to the applied force. While radical abstraction is overall beneficial for mechanical bond scission, the polymer also becomes less responsive to force than the unactivated polymer. This finding places upper limits on the efficacy of the synergy between radical attack and applied force. In addition, the importance of reaction pathway optimization is also shown, where comparisons to the COGEF method show the latter to be qualitatively incapable of describing chain scission after radical activation.

Iron(III)-catalysed carbonyl-olefin metathesis.

The olefin metathesis reaction of two unsaturated substrates is one of the most powerful carbon-carbon-bond-forming reactions in organic chemistry. Specifically, the catalytic olefin metathesis reaction has led to profound developments in the synthesis of molecules relevant to the petroleum, materials, agricultural and pharmaceutical industries. These reactions are characterized by their use of discrete metal alkylidene catalysts that operate via a well-established mechanism. While the corresponding carbonyl-olefin metathesis reaction can also be used to construct carbon-carbon bonds, currently available methods are scarce and severely hampered by either harsh reaction conditions or the required use of stoichiometric transition metals as reagents. To date, no general protocol for catalytic carbonyl-olefin metathesis has been reported. Here we demonstrate a catalytic carbonyl-olefin ring-closing metathesis reaction that uses iron, an Earth-abundant and environmentally benign transition metal, as a catalyst. This transformation accommodates a variety of substrates and is distinguished by its operational simplicity, mild reaction conditions, high functional-group tolerance, and amenability to gram-scale synthesis. We anticipate that these characteristics, coupled with the efficiency of this reaction, will allow for further advances in areas that have historically been enhanced by olefin metathesis.

State-specific solvation for restricted active space spin-flip (RAS-SF) wave functions based on the polarizable continuum formalism.

The restricted active space spin-flip (RAS-SF) formalism is a particular form of single-reference configuration interaction that can describe some forms of strong correlation at a relatively low cost and which has recently been formulated for the description of charge-transfer excited states. Here, we introduce both equilibrium and nonequilibrium versions of a state-specific solvation correction for vertical transition energies computed using RAS-SF wave functions, based on the framework of a polarizable continuum model (PCM). Ground-state polarization is described using the solvent's static dielectric constant and in the nonequilibrium solvation approach that polarization is modified upon vertical excitation using the solvent's optical dielectric constant. Benchmark calculations are reported for well-studied models of photo-induced charge transfer, including naphthalene dimer, C2H4⋯C2F4, pentacene dimer, and perylene diimide (PDI) dimer, several of which are important in organic photovoltaic applications. For the PDI dimer, we demonstrate that the charge-transfer character of the excited states is enhanced in the presence of a low-dielectric medium (static dielectric constant ɛ0 = 3) as compared to a gas-phase calculation (ɛ0 = 1). This stabilizes mechanistic traps for singlet fission and helps to explain experimental singlet fission rates. We also examine the effects of nonequilibrium solvation on charge-separated states in an intramolecular singlet fission chromophore, where we demonstrate that the energetic ordering of the states changes as a function of solvent polarity. The RAS-SF + PCM methodology that is reported here provides a framework to study charge-separated states in solution and in photovoltaic materials.

Studies of Catalyst-Controlled Regioselective Acetalization and Its Application to Single-Pot Synthesis of Differentially Protected Saccharides.

This article describes studies on the regioselective acetal protection of monosaccharide-based diols using chiral phosphoric acids (CPAs) and their immobilized polymeric variants, (R)-Ad-TRIP-PS and (S)-SPINOL-PS, as the catalysts. These catalyst-controlled regioselective acetalizations were found to proceed with high regioselectivities (up to >25:1 rr) on various d-glucose-, d-galactose-, d-mannose-, and l-fucose-derived 1,2-diols and could be carried out in a regiodivergent fashion depending on the choice of chiral catalyst. The polymeric catalysts were conveniently recycled and reused multiple times for gram-scale functionalizations with catalytic loadings as low as 0.1 mol %, and their performance was often found to be superior to the performance of their monomeric variants. These regioselective CPA-catalyzed acetalizations were successfully combined with common hydroxyl group functionalizations as single-pot telescoped procedures to produce 32 regioisomerically pure differentially protected mono- and disaccharide derivatives. To further demonstrate the utility of the polymeric catalysts, the same batch of (R)-Ad-TRIP-PS catalyst was recycled and reused to accomplish single-pot gram-scale syntheses of 6 differentially protected d-glucose derivatives. The subsequent exploration of the reaction mechanism using NMR studies of deuterated and nondeuterated substrates revealed that low-temperature acetalizations happen via a syn-addition mechanism and that the reaction regioselectivity exhibits strong dependence on the temperature. The computational studies indicate a complex temperature-dependent interplay of two reaction mechanisms, one involving an anomeric phosphate intermediate and another via concerted asynchronous formation of an acetal, that results in syn-addition products. The computational models also explain the steric factors responsible for the observed C2 selectivities and are consistent with experimentally observed selectivity trends.

Brønsted-Acid-Catalyzed Intramolecular Carbonyl-Olefin Reactions: Interrupted Metathesis vs Carbonyl-Ene Reaction.

Lewis acid catalysts have been shown to promote carbonyl-olefin metathesis through a critical four-membered-ring oxetane intermediate. Recently, Brønsted-acid catalysis of related substrates was similarly proposed to result in a transient oxetane, which fragments within a single elementary step via a postulated oxygen-atom transfer mechanism. Herein, careful quantum chemical investigations show that Brønsted acid (triflic acid, TfOH) instead invokes a mechanistic switch to a carbonyl-ene reaction, and oxygen-atom transfer is uncompetitive. TfOH's conjugate base is also found to rearrange H atoms and allow isomerization of the carbocations that appear after the carbonyl-ene reaction. The mechanism explains available experimental information, including the skipped diene species that appear transiently before product formation. The present study clarifies the mechanism for activation of intramolecular carbonyl-olefin substrates by Brønsted acids and provides important insights that will help develop this exciting class of catalysts.

Using JPP to Identify Ni Bidentate Phosphine Complexes In Situ.

Identifying intermediates of Ni-containing reactions can be challenging due to the high reactivity of Ni complexes and their sensitivity toward air and moisture. Many Ni bidentate phosphine complexes are diamagnetic and can be analyzed in situ via 31P NMR spectroscopy, but the oxidation state of Ni is difficult to determine using 31P chemical shift analysis alone. The J-coupling between P atoms, JPP, has been proposed to correlate with oxidation state, but few investigations have looked at how JPP is affected by parameters such as length of the linker or identity of the phosphine or other ligands. The present investigation into the JPP values of Ni bidentate phosphine complexes with two-carbon and three-carbon linkers shows that the JPP values observed in 31P NMR spectra, |JPP|, are competent indicators of the oxidation state at Ni. For complexes with two-carbon linkers, |JPP| > 40 Hz is typical of Ni0 while |JPP| < 30 Hz is typical of NiII; this trend is reversed for complexes with three-carbon linkers. Additionally, the Lewis acidity of the Ni and Lewis basicity of the phosphine ligand affect JPP predictably. For example, increased P-to-Ni donation arising from more-donating phosphines or more-withdrawing ligands trans to the P atoms causes a more negative JPP. These results should enable the oxidation state of Ni and properties of ligands in Ni bidentate phosphine complexes to be determined in situ during reactions containing these species.

Toward Full Configuration Interaction for Transition-Metal Complexes.

An efficacious approximation to full configuration interaction (FCI) is adapted to calculate singlet-triplet gaps for transition-metal complexes. This strategy, incremental FCI (iFCI), uses a many-body expansion to systematically add correlation to a simple reference wave function and therefore achieves greatly reduced computational costs compared to FCI. iFCI through the 3-body expansion is demonstrated on four model transition-metal complexes involving the metals Zn, V, and Cu. Screening techniques to increase the computational efficiency of iFCI are proposed and tested, showing reduction in the number of 3-body terms by more than 90% with controlled errors. The largest complex treated herein by iFCI has 142 valence electrons, all of which are correlated among the full set of 444 active orbitals. Computed spin gaps approach experimental results for the four complexes, though room for improvement remains.

Fully variational incremental CASSCF.

The complete-active-space self-consistent field (CASSCF) method is a canonical electronic structure theory that holds a central place in conceptualizing and practicing first principles simulations. For application to realistic molecules, however, the CASSCF must be approximated to circumvent its exponentially scaling computational costs. Applying the many-body expansion-also known as the method of increments-to CASSCF (iCASSCF) has been shown to produce a polynomially scaling method that retains much of the accuracy of the parent theory and is capable of treating full valence active spaces. Due to an approximation made in the orbital gradient, the orbital parameters of the original iCASSCF formulation could not be variationally optimized, which limited the accuracy of its nuclear gradient. Herein, a variational iCASSCF is introduced and implemented, where all parameters are fully optimized during energy minimization. This method is able to recover electronic correlations from the full valence space in large systems, produce accurate gradients, and optimize stable geometries as well as transition states. Demonstrations on challenging test cases, such as the oxoMn(salen)Cl complex with 84 electrons in 84 orbitals and the automerization of cyclobutadiene, show that the fully variational iCASSCF is a powerful tool for describing challenging molecular chemistries.

Mechanochemical Release of N-Heterocyclic Carbenes from Flex-Activated Mechanophores.

We have discovered a new flex-activated mechanophore that releases an N-heterocyclic carbene (NHC) under mechanical load. The mechanophore design is based upon NHC-carbodiimide (NHC-CDI) adducts and demonstrates an important first step toward flex-activated designs capable of further downstream reactivities. Since the flex-activation is non-destructive to the main polymer chains, the material can be subjected to multiple compression cycles to achieve iterative increases in the activation percentage of mechanophores. Two different NHC structures were demonstrated, signifying the potential modularity of the mechanophore design.