Stereocontrol in Diphenylprolinol Silyl Ether Catalyzed Michael Additions: Steric Shielding or Curtin–Hammett Scenario?

Development of analytic gradients for the Huzinaga quantum embedding method and its applications to large-scale hybrid and double hybrid DFT forces

The theory of analytic gradients is presented for the projector-based density functional theory (DFT) embedding approach utilizing the Huzinaga-equation. The advantages of the Huzinaga-equation-based formulation are demonstrated. In particular, it is shown that the projector employed does not appear in the Lagrangian, and the potential risk of numerical problems is avoided at the evaluation of the gradients. The efficient implementation of the analytic gradient theory is presented for approaches where hybrid DFT, second-order Møller-Plesset perturbation theory, or double hybrid DFT are embedded in lower-level DFT environments. To demonstrate the applicability of the method and to gain insight into its accuracy, it is applied to equilibrium geometry optimizations, transition state searches, and potential energy surface scans. Our results show that bond lengths and angles converge rapidly with the size of the embedded system. While providing structural parameters close to high-level quality for the embedded atoms, the embedding approach has the potential to relax the coordinates of the environment as well. Our demonstrations on a 171-atom zeolite and a 570-atom protein system show that the Huzinaga-equation-based embedding can accelerate (double) hybrid gradient computations by an order of magnitude with sufficient active regions and enables affordable force evaluations or geometry optimizations for molecules of hundreds of atoms.

Linear-Scaling Local Natural Orbital CCSD(T) Approach for Open-Shell Systems: Algorithms, Benchmarks, and Large-Scale Applications

The extension of the highly optimized local natural orbital (LNO) coupled cluster (CC) with single-, double-, and perturbative triple excitations [LNO-CCSD(T)] method is presented for high-spin open-shell molecules based on restricted open-shell references. The techniques enabling the outstanding efficiency of the closed-shell LNO-CCSD(T) variant are adopted, including the iteration- and redundancy-free second-order Møller-Plesset and (T) formulations as well as the integral-direct, memory- and disk use-economic, and OpenMP-parallel algorithms. For large molecules, the efficiency of our open-shell LNO-CCSD(T) method approaches that of its closed-shell parent method due to the application of restricted orbital sets for demanding integral transformations and a novel approximation for higher-order long-range spin-polarization effects. The accuracy of open-shell LNO-CCSD(T) is extensively tested for radicals and reactions thereof, ionization processes, as well as spin-state splittings, and transition-metal compounds. At the size range where the canonical CCSD(T) reference is accessible (up to 20-30 atoms), the average open-shell LNO-CCSD(T) correlation energies are found to be 99.9 to 99.95% accurate, which translates into average absolute deviations of a few tenths of kcal/mol in the investigated energy differences already with the default settings. For more extensive molecules, the local errors may grow, but they can be estimated and decreased via affordable systematic convergence studies. This enables the accurate modeling of large systems with complex electronic structures, as illustrated on open-shell organic radicals and transition-metal complexes of up to 179 atoms as well as on challenging biochemical systems, including up to 601 atoms and 11,000 basis functions. While the protein models involve difficulties for local approximations, such as the spin states of a bounded iron ion or an extremely delocalized singly occupied orbital, the corresponding single-node LNO-CCSD(T) computations were feasible in a matter of days with 10s to 100 GB of memory use. Therefore, the new LNO-CCSD(T) implementation enables highly accurate computations for open-shell systems of unprecedented size and complexity with widely accessible hardware.

Highly efficient morpholine-based organocatalysts for the 1,4-addition reaction between aldehydes and nitroolefins: an unexploited class of catalysts

Many studies have demonstrated how the pyrrolidine nucleus is more efficient than the corresponding piperidine or morpholine as organocatalysts in the condensation of aldehydes with electrophiles via enamine. Focussing on morpholine-enamines, their low reactivity is ascribed to the presence of oxygen on the ring and to the pronounced pyramidalisation of nitrogen, decreasing the nucleophilicity of the enamine. Thus, the selection of efficient morpholine organocatalysts appears to be a difficult challenge. Herein, we reported on the synthesis of new organocatalysts belonging to the class of ß-morpholine amino acids that were tested in a model reaction, i.e., the 1,4-addition reaction of aldehydes to nitroolefins. Starting from commercially available amino acids and epichlorohydrin, we designed an efficient synthesis for the aforementioned catalysts, controlling the configuration and the substitution pattern. Computational studies indeed disclosed the transition state of the reaction, explaining why, despite all the limitations of the morpholine ring for enamine catalysis, our best catalyst works efficiently, affording condensation products with excellent yields, diastereoselection and good-to-exquisite enantioselectivity.

Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach

Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).

Origins of Stereospecificity and Divergent Reactivity of Pd-Catalyzed Cross Coupling with α,α-Disubstituted Alkenyl Hydrazones

This article presents an exploration of stereospecificity and divergent reactivity of Pd-catalyzed α,α-disubstituted alkenyl hydrazones to synthesize 1,4-dienes in the Z configuration and vinylcyclopropane. We calculated the energy profiles of four α,α-disubstituted alkenyl hydrazones. The results show that the energy profiles of the whole catalytic cycle are basically the same before the syn-carbopalladation step. Subsequent syn-β-C elimination yields skipping dienes, or direct β-H elimination yields vinylcyclopropane. Current theoretical calculations reveal that the stereospecificity and the divergent reactivity of reactions result from the competition between syn-β-C elimination and β-H elimination. The C-C bond rotation and subsequent syn-β-C elimination step control the stereospecificity of the reaction by changing the olefin stereostructure from E to Z configuration. The steric factor of α-substituted groups mediates the transformation between syn-β-C elimination and β-H elimination. The results are of great significance for the scientific design of substrates to achieve accurate synthesis of target products.

Providing direction for mechanistic inferences in radical cascade cyclization using Transformer model

Even in modern organic chemistry, predicting or proposing a reaction mechanism and speculating on reaction intermediates remains challenging. For example, it is challenging to predict the regioselectivity of radical attraction...

Insight into the mechanism of the arylation of arenes via norbornene relay palladation through meta- to para-selectivity

A theoretical insight was shown into the origin of site-selectivity in the arylation of arenes by a norbornene relay palladation through meta- to para-selectivity.

Accurate Reduced-Cost CCSD(T) Energies: Parallel Implementation, Benchmarks, and Large-Scale Applications

The accurate and systematically improvable frozen natural orbital (FNO) and natural auxiliary function (NAF) cost-reducing approaches are combined with our recent coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] implementations. Both of the closed- and open-shell FNO-CCSD(T) codes benefit from OpenMP parallelism, completely or partially integral-direct density-fitting algorithms, checkpointing, and hand-optimized, memory- and operation count effective implementations exploiting all permutational symmetries. The closed-shell CCSD(T) code requires negligible disk I/O and network bandwidth, is MPI/OpenMP parallel, and exhibits outstanding peak performance utilization of 50-70% up to hundreds of cores. Conservative FNO and NAF truncation thresholds benchmarked for challenging reaction, atomization, and ionization energies of both closed- and open-shell species are shown to maintain 1 kJ/mol accuracy against canonical CCSD(T) for systems of 31-43 atoms even with large basis sets. The cost reduction of up to an order of magnitude achieved extends the reach of FNO-CCSD(T) to systems of 50-75 atoms (up to 2124 atomic orbitals) with triple- and quadruple-ζ basis sets, which is unprecedented without local approximations. Consequently, a considerably larger portion of the chemical compound space can now be covered by the practically "gold standard" quality FNO-CCSD(T) method using affordable resources and about a week of wall time. Large-scale applications are presented for organocatalytic and transition-metal reactions as well as noncovalent interactions. Possible applications for benchmarking local CCSD(T) methods, as well as for the accuracy assessment or parametrization of less complete models, for example, density functional approximations or machine learning potentials, are also outlined.

Amino-Catalyzed Reactions of Aldehydes with Chiral Nitroalkenes

Chiral nitroalkenes are used for the first time in Michael additions of aldehydes, catalyzed by pyrrolidine derivatives. They yield the same major stereoisomer with either (S)-proline or (R)-proline, but this asymmetric induction does not overcome the effect of sterically more congested catalysts. Nitrocyclobutane intermediates are often formed, which are more stable than those from (E)-1-nitro-2-phenylethene. The cyclobutanes and final products were characterized by 2D NMR and chemical correlations.

Deactivation of Secondary Amine Catalysts via Aldol Reaction-Amine Catalysis under Solvent-Free Conditions

Despite intense interest in amine-catalyzed stereoselective reactions, high catalyst loadings of ≥10 mol % are still common and either due to low reactivity or catalyst deactivation. Yet, few deactivation pathways are well understood. Here, we unraveled the deactivation of secondary amines by undesired aldol reaction. Mechanistic studies with peptide and prolinol silyl ether catalysts showed the generality of this so-far underappreciated catalyst deactivation pathway. The insights enabled conjugate addition reactions between aldehydes and nitroolefins on a multigram scale in the absence of solvent-conditions that are attractive as environmentally benign processes-with excellent product yields and stereoselectivities in the presence of as little as 0.1 mol % of a chemoselective peptidic catalyst.

Effect of the enamine pyramidalization direction on the reactivity of secondary amine organocatalysts

Chiral secondary amines are valuable catalysts for reactions that proceed through an enamine intermediate. Here, we explored the importance of the pyramidalization direction of the enamine-N on the reactivity of chiral enamines with a combination of computational, NMR spectroscopic, and kinetic experiments. Studies with peptidic catalysts that bear cyclic amines with different ring sizes revealed that endo-pyramidalized enamines are significantly more reactive compared to exo-pyramidalized analogs. The results show that the pyramidalization direction can have a greater effect than n→π* orbital overlap on the reactivity of chiral enamines. The data enabled the development of a catalyst with higher reactivity compared to the parent catalyst.

Endo-pyramidalisation at nitrogen bestows enamines derived from α-substituted amines with higher reactivity compared to exo-pyramidalisation.

The origin of regioselectivity in Cu-catalyzed hydrocarbonylative coupling of alkynes with alkyl halides

Steric interactions mediate a switch between a ketone and allylic alcohol in Cu-catalyzed hydrocarbonylative coupling of alkynes with alkyl halides.

NMR and Computational Studies on the Reactions of Enamines with Nitroalkenes That May Pass through Cyclobutanes

The addition of aldehyde enamines to nitroalkenes affords cyclobutanes in all solvents, with all of the pyrrolidine and proline derivatives tested by us and with all of the substrates we have examined. Depending on the temperature, concentration of water, solvent polarity, and other factors, the opening and hydrolysis of such a four-membered ring may take place rapidly or last for several days, producing the final Michael-like adducts (4-nitrobutanals). Thirteen new cyclobutanes have now been characterized by NMR spectroscopy. As could be expected, s-trans-enamine conformers give rise to all-trans-(4S)-4-nitrocyclobutylpyrrolidines, while s-cis-enamine conformers afford all-trans-(4R)-4-nitrocyclobutylpyrrolidines. These four-membered rings can isomerize to adduct enamines, which should be hydrolyzed via their iminium ions. MP2 and M06-2X calculations predict that one iminium ion is more stable than the other iminium species, so that protonation of the adduct enamines can be quite stereoselective; in the presence of water, the so-called syn adducts (e.g., OCH-*CHR-*CHPh-CH₂NO₂, with R and Ph syn) eventually become the major products. Why one syn adduct is obtained with aldehydes, whereas cyclic ketones (the predicted ring-fused cyclobutanes of which isomerize to their enamines more easily) produce the other syn adduct, is also explained by means of molecular orbital calculations. Nitro-Michael reactions of aldehyde enamines that "stop" at the nitrocyclobutane stage and final enamine stage do not work catalytically, as known, but those of cyclic ketone enamines that do not work stop at the final enamine stage (if their hydrolysis to the corresponding nitroethylketones is less favorable than expected). These and other facts are accounted for, and the proposals of the groups led by Seebach and Hayashi, Blackmond, and Pihko and Papai are reconciled.

A Computational Mechanistic Study of Pd(II)-Catalyzed Enantioselective C(sp3)-H Borylation: Roles of APAO Ligands

A computational mechanistic study has been performed on Pd(II)-catalyzed enantioselective reactions involving acetyl-protected aminomethyl oxazolines (APAO) ligands that significantly improved reactivity and selectivity in C(sp³)-H borylation. The results support a mechanism including initiation of C(sp³)-H bond activation generating a five-membered palladacycle and ligand exchange, followed by HPO₄^2--promoted transmetalation. These resulting Pd(II) complexes further undergo sequential reductive elimination by coordination of APAO ligands and protonation to afford the enantiomeric products and deliver Pd(0) complexes, which will then proceed by oxidation and deprotonation to regenerate the catalyst. The C(sp³)-H activation is found to be the rate- and enantioselectivity-determining step, in which the APAO ligand acts as the proton acceptor to form the two enantioselectivity models. The results demonstrate that the diverse APAO ligands control the enantioselectivity by differentiating the distortion and interaction between the major and minor pathways.

Organocatalytic Asymmetric Addition of Aldehyde to Nitroolefin by H-d-Pro-Pro-Glu-NH2: A Mechanistic Study

The mechanism of the asymmetric addition of aldehyde (butanal) to nitroolefin (β-nitrostyrene) catalyzed by H-d-Pro-Pro-Glu-NH₂ (dPPE-NH₂; 1) was explored using density functional theory methods in chloroform. By conformational search, it was confirmed that catalyst 1 and its enamine intermediate adopted a dominant conformation with a βI structure stabilized by a C₁₀ H-bond between the C=O of d-Pro1 and C-terminal NH₂ proton and by an additional H-bond between the side chain and the backbone of Glu3. This βI turn structure was conserved all along the catalytic cycle. Consistently with the kinetic studies, the C-C bond formation between the enamine and electrophile was also confirmed as the rate-determining step. The stereoselectivity results from a re → re prochiral approach of enamine and β-nitrostyrene with a gauche^- orientation of the double bonds. Although it was suggested as the possible formation of dihydrooxazine oxide species, this process was confirmed to be kinetically less accessible than the formation of acyclic nitronate. In particular, our calculated results supported that the carboxylic acid group of Glu3 in 1 played a central role by acting as general acid/base all along the catalytic cycle and orienting the asymmetric C-C bond formation.

Mechanistic Studies on the Organocatalytic α-Chlorination of Aldehydes: The Role and Nature of Off-Cycle Intermediates

Herein we report the isolation and characterization of aminal intermediates in the organocatalytic α-chlorination of aldehydes. These species are stable covalent ternary adducts of the substrate, the catalyst and the chlorinating reagent. NMR-assisted kinetic studies and isotopic labeling experiments with the isolated intermediate did not support its involvement in downstream stereoselective processes as proposed by Blackmond. By tuning the reactivity of the chlorinating reagent, we were able to suppress the accumulation of rate-limiting off-cycle intermediates. As a result, an efficient and highly enantioselective catalytic system with a broad functional group tolerance was developed.

Asymmetric Conjugate Addition of α,α-Disubstituted Aldehydes to Nitroalkenes Organocatalyzed by Chiral Monosalicylamides from trans-Cyclohexane-1,2-Diamines

Primary amine-salicylamides derived from chiral trans-cyclohexane-1,2-diamines are used as organocatalysts for the enantioselective conjugate addition of α,α-disubstituted aldehydes to arylated and heteroarylated nitroalkenes. The reaction is performed in the presence of 4-dimethylaminopyridine as an additive in dichloromethane as a solvent at room temperature. The corresponding enantioenriched γ-nitroaldehydes are obtained with enantioselectivities up to 95%. Theoretical calculations are used to justify the reasons of the stereoinduction.