Nature and Strength of Lewis Acid/Base Interaction in Boron and Nitrogen Trihalides

In-silico-assisted derivatization of triarylboranes for the catalytic reductive functionalization of aniline-derived amino acids and peptides with H2

Cheminformatics-based machine learning (ML) has been employed to determine optimal reaction conditions, including catalyst structures, in the field of synthetic chemistry. However, such ML-focused strategies have remained largely unexplored in the context of catalytic molecular transformations using Lewis-acidic main-group elements, probably due to the absence of a candidate library and effective guidelines (parameters) for the prediction of the activity of main-group elements. Here, the construction of a triarylborane library and its application to an ML-assisted approach for the catalytic reductive alkylation of aniline-derived amino acids and C-terminal-protected peptides with aldehydes and H2 is reported. A combined theoretical and experimental approach identified the optimal borane, i.e., B(2,3,5,6-Cl4-C6H)(2,6-F2-3,5-(CF3)2-C6H)2, which exhibits remarkable functional-group compatibility toward aniline derivatives in the presence of 4-methyltetrahydropyran. The present catalytic system generates H2O as the sole byproduct.

Acidity of Isomorphic Substituted Zeolites with B, Al and Ga Revisited

Isomorphic substitution of zeolites with B, Al and Ga is a widely used approach in catalysis. The experimentally reported trend of their acidities decreases in the order: Al>Ga>B. However, a consistent explanation is still lacking in the literature. To bring more understanding of this trend, density functional theory computations were conducted on several model systems. First, the acidity of small clusters with two (2T) and five (5T) tetrahedral sites was analyzed. These systems were then projected onto three large void structures: H-[A]-BEA (52T), H-[A]-FAU (84T) and H-[A]-MOR (112T) with A=B, Al, Ga. Our electron density and Interacting Quantum Atom analyses show that the acidity of Al zeolites originates from the much stronger O-Al bond, which is dominated by the electrostatic attraction. The bridging hydroxyl therefore donates more charge density to the metal, the proton becomes more positive and consequently more acidic. Ga zeolites are more acidic than B zeolites due to the greater covalent nature on the O-Ga bond. The resulting acidity, as seen by ammonia, depends on both the acidic oxygen and the charge distribution of the surrounding oxygens exerted by the substituents.

Predicting Lewis Acidity: Machine Learning the Fluoride Ion Affinity of p-Block-Atom-Based Molecules

"How strong is this Lewis acid?" is a question researchers often approach by calculating its fluoride ion affinity (FIA) with quantum chemistry. Here, we present FIA49k, an extensive FIA dataset with 48,986 data points calculated at the RI-DSD-BLYP-D3(BJ)/def2-QZVPP//PBEh-3c level of theory, including 13 different p-block atoms as the fluoride accepting site. The FIA49k dataset was used to train FIA-GNN, two message-passing graph neural networks, which predict gas and solution phase FIA values of molecules excluded from training with a mean absolute error of 14 kJ mol-1 (r2=0.93) from the SMILES string of the Lewis acid as the only input. The level of accuracy is notable, given the wide energetic range of 750 kJ mol-1 spanned by FIA49k. The model's value was demonstrated with four case studies, including predictions for molecules extracted from the Cambridge Structural Database and by reproducing results from catalysis research available in the literature. Weaknesses of the model are evaluated and interpreted chemically. FIA-GNN and the FIA49k dataset can be reached via a free web app (www.grebgroup.de/fia-gnn).

Cobalt-Catalyzed Enantio- and Regioselective C(sp3 )-H Alkenylation of Thioamides

An enantioselective cobalt-catalyzed C(sp3 )-H alkenylation of thioamides with but-2-ynoate ester coupling partners employing thioamide directing groups is presented. The method is operationally simple and requires only mild reaction conditions, while providing alkenylated products as single regioisomers in excellent yields (up to 85 %) and high enantiomeric excess [up to 91 : 9 enantiomeric ratio (er), or up to >99 : 1 er after a single recrystallization]. Diverse downstream derivatizations of the products are demonstrated, delivering a range of enantioenriched constructs. Extensive computational studies using density functional theory provide insight into the detailed reaction mechanism, origin of enantiocontrol, and the unusual regioselectivity of the alkenylation reaction.

The Field of Main Group Lewis Acids and Lewis Superacids: Important Basics and Recent Developments

New developments in the field of Lewis acidity are highlighted, with the focus of novel Lewis acids and Lewis superacids of group 2, 13, 14, and 15 elements. Several important basics, illustrated by modern examples (classification of Donor-Acceptor (DA) complexes, amphoteric nature of any compound in terms of DA interactions, reorganization energies of main group Lewis acids and the role of the energies of frontier orbitals) are presented and discussed. It is emphasized that the Lewis acidity phenomena are general and play vital role in different areas of chemistry: from weak "atomophilic" interactions to the complexes of Lewis superacids.

Blueshift in Trifurcated Hydrogen Bonds: A Tradeoff between Tetrel Bonding and Steric Repulsion

We have quantum chemically investigated the origin of the atypical blueshift of the H-C bond stretching frequency in the hydrogen-bonded complex X- •••H3 C-Y (X, Y=F, Cl, Br, I), as compared to the corresponding redshift occurring in Cl- •••H3 N and Cl- •••H3 C-H, using relativistic density functional theory (DFT) at ZORA-BLYP-D3(BJ)/QZ4P. Previously, this blueshift was attributed, among others, to the contraction of the H-C bonds as the H3 C moiety becomes less pyramidal. Herein, we provide quantitative evidence that, instead, the blueshift arises from a direct and strong X- •••C interaction of the HOMO of A- with the backside lobe on carbon of the low-lying C-Y antibonding σ* LUMO of the H3 C-Y fragment. This X- •••C bond, in essence a tetrel bond, pushes the H atoms towards a shorter H-C distance and makes the H3 C moiety more planar. The blueshift may, therefore, serve as a diagnostic for tetrel bonding.

Atom-Economical Synthesis of Lewis Acidic Boron Containing Porous Organic Polymers via Hydroboration Polymerization for Basic Chemical Capture

Atom economy is one of the main concerns for material synthesis. Here, the facile synthesis of Lewis acidic boron-containing porous organic polymers (B-POPs) via hydroboration polymerization reaction of commercially available borane dimethyl sulfide complex (BH3 ∙SMe2 ) with multi-alkynes under mild reaction conditions is presented. This new synthetic method for B-POPs has the advantage of high atom economy. The resulted porous alkenyl borane polymers (PABPs) have unique features such as high boron content, strong Lewis acidity, and high surface areas. Owing to the strong Lewis acid-base interactions, PABPs exhibit excellent adsorptive capacity toward triethylamine (up to 841 mg g-1 ) and pyridine (up to 1396 mg g-1 ) vapor.

Mechanistic Analysis of Alkyne Haloboration: A DFT, MP2, and DLPNO-CCSD(T) Study

Stereocontrol of the alkyne haloboration reaction has received attention in many experimental but few theoretical studies. Here we present a detailed quantum-chemical study of mechanisms leading to Z versus E isomers of haloboration products, considering acetylene and propyne combined with BCl3, BBr3, and BI3. Calculations using B3LYP-D3, MP2, and DLPNO-CCSD(T) methods are used to study polar reactions between the alkyne and BX3 in the absence and presence of an additional halide anion whose content in the reaction mixture can be controlled experimentally. The formation of anti-haloboration products via radical mechanisms is also explored, namely, by adding BX3 to (Z)-halovinyl radical. For the anti-haloboration of propyne, the radical route is prohibited by the regiochemistry of the initiating halopropenyl radical, while the polar route is unlikely due to a competitive allene generation. In contrast, energetically accessible routes exist for both syn- and anti-bromoboration of acetylene; hence, careful control of reaction conditions is necessary to steer the stereochemical outcome. Methodologically, MP2 results correspond better to the DLPNO-CCSD(T) energies than the B3LYP-D3 results in terms of both reaction barrier heights and relative ordering of energetically close stationary points.

An easy-to-perform evaluation of steric properties of Lewis acids

Steric and electronic effects play a very important role in chemistry, as these effects influence the shape and reactivity of molecules. Herein, an easy-to-perform approach to assess and quantify steric properties of Lewis acids with differently substituted Lewis acidic centers is reported. This model applies the concept of the percent buried volume (%V _Bur) to fluoride adducts of Lewis acids, as many fluoride adducts are crystallographically characterized and are frequently calculated to judge fluoride ion affinities (FIAs). Thus, data such as cartesian coordinates are often easily available. A list of 240 Lewis acids together with topographic steric maps and cartesian coordinates of an oriented molecule suitable for the SambVca 2.1 web application is provided, together with different FIA values taken from the literature. Diagrams of %V _Bur as a scale for steric demand vs. FIA as a scale for Lewis acidity provide valuable information about stereo-electronic properties of Lewis acids and an excellent evaluation of steric and electronic features of the Lewis acid under consideration. Furthermore, a novel LAB-Rep model (Lewis acid/base repulsion model) is introduced, which judges steric repulsion in Lewis acid/base pairs and helps to predict if an arbitrary pair of Lewis acid and Lewis base can form an adduct with respect to their steric properties. The reliability of this model was evaluated in four selected case studies, which demonstrate the versatility of this model. For this purpose, a user-friendly Excel spreadsheet was developed and is provided in the ESI, which works with listed buried volumes of Lewis acids %V _{Bur_LA} and of Lewis bases %V _{Bur_LB}, and no results from experimental crystal structures or quantum chemical calculations are necessary to evaluate steric repulsion in these Lewis acid/base pairs.

Lewis Acidity Trend of Boron and Aluminium Trihalides: If Not π Back-Bonding, What Else?

Lewis acidity trend of boron trihalides is a subject that has received a variety of explanations, and still, the simple π back-bonding based one is believed by most, perhaps because of its simplicity, irrespective of opposing findings. Herein we try to give an alternative explanation based on qualitative Molecular Orbital (MO) theory and support that quantitatively by Generalized Kohn-Sham Energy Decomposition Analysis. While the role of orbital overlap on the orbital interaction energy is widely known, the role of electronegativity of the atoms involved is often overlooked. Here we find that the Lewis acidity trend of boron and aluminium halides can be explained by the Wolfsberg-Helmholz (W-H) formula for resonance integral. The MO theory-based predictions are valid only when the orbital interactions are strong enough. In weakly interacting systems, the effect of orbital interactions can be overshadowed by other effects such as Pauli repulsion, dispersion, etc. Thus the Lewis acidity trend of boron and aluminium halides can vary depending on the strength of the interacting Lewis base. We believe that this work would enable one to gain a better understanding not only on the Lewis acidity of boron trihalides and its heavy analogs but also on a variety of related problems such as the stronger π acidity of CS compared to CO and weaker π bonding between heavy atoms.

Computed ammonia affinity for evaluating Lewis acidity of organoboronates and organoboronamides

Lewis acidity of organoboronates [B(pin), B(neop), B(cat), B(eg), B(nad)] and organoboronamides [B(dan), B(aam), B(mdan)] has been found to be unifiedly evaluated by computed ammonia affinity (AA), while other methods [LUMO energies, global electrophilicity index (GEI), fluoride ion affinity (FIA)] were only partially applicable. The relationships between the AA values and such structural characters including the B-X bond lengths, the X-B-X angles, and the changes in the B-X bond lengths in the formation of the ammonia adducts were also described.

Methyl Substitution Destabilizes Alkyl Radicals

We have quantum chemically investigated how methyl substituents affect the stability of alkyl radicals Me_mH_3−mC⋅ and the corresponding Me_mH_3−mC−X bonds (X = H, CH₃, OH; m = 0 – 3) using density functional theory at M06‐2X/TZ2P. The state‐of‐the‐art in physical organic chemistry is that alkyl radicals are stabilized upon an increase in their degree of substitution from methyl<primary<secondary<tertiary, and that this is the underlying cause for the decrease in C−H bond strength along this series. Here, we provide evidence that falsifies this model and show that, on the contrary, the Me_mH_3−mC⋅ radical is destabilized with increasing substitution. The reason that the corresponding C−H bond nevertheless becomes weaker is that substitution destabilizes the sterically more congested Me_mH_3−mC−H molecule even more.

What Distinguishes the Strength and the Effect of a Lewis Acid: Analysis of the Gutmann-Beckett Method

IUPAC defines Lewis acidity as the thermodynamic tendency for Lewis pair formation. This strength property was recently specified as global Lewis acidity (gLA), and is gauged for example by the fluoride ion affinity. Experimentally, Lewis acidity is usually evaluated by the effect on a bound molecule, such as the induced 31 P NMR shift of triethylphosphine oxide in the Gutmann-Beckett (GB) method. This type of scaling was called effective Lewis acidity (eLA). Unfortunately, gLA and eLA often correlate poorly, but a reason for this is unknown. Hence, the strength and the effect of a Lewis acid are two distinct properties, but they are often granted interchangeably. The present work analyzes thermodynamic, NMR specific, and London dispersion effects on GB numbers for 130 Lewis acids by theory and experiment. The deformation energy of a Lewis acid is identified as the prime cause for the critical deviation between gLA and eLA but its correction allows a unification for the first time.

The inversion of tetrahedral p-block element compounds: general trends and the relation to the second-order Jahn-Teller effect

The tetrahedron is the primary structural motif among the p-block elements and determines the architecture of our bio- and geosphere. However, a broad understanding of the configurational inversion of tetrahedral compounds is missing. Here, we report over 250 energies (DLPNO-CCSD(T)) for square planar inversion of third- and fourth-period element species of groups 13, 14, and 15. Surprisingly low inversion barriers are identified for compounds of industrial relevance (e.g., ≈100 kJ mol-1 for Al(OH)4 -). More fundamentally, the second-order Jahn-Teller theorem is disclosed as suitable to rationalize substituent and central element effects. Bond analysis tools give further insights into the preference of eight valence electron systems with four substituents to be tetrahedral. Hence, this study develops a model to understand, memorize, and predict the angular flexibility of tetrahedral species. Perceiving the tetrahedron not as forcingly rigid but as a dynamic structural entity might leverage new approaches and visions for adaptive matter.

Origin of the Enhanced Reactivity in the ortho C-H Borylation of Benzaldehydes with BBr3

The metal-free ortho C-H borylation of benzaldehyde derivatives using a transient imine directing group was recently developed by our group, providing an efficient strategy for the synthesis of organoboron reagents. Herein, we report on an extensive investigation of the reaction mechanism using density functional theory (DFT) calculations. Computations for the reaction pathway with various imine substrates, as well as the effect of an added base were examined, and the experimentally observed reactivity enhancement is proposed to originate from the tunability of the destabilizing strain energies that results in a reversible complexation process with BBr₃.

Boron Tunneling in the "Weak" Bond-Stretch Isomerization of N-B Lewis Adducts

Some nitrile-boron halide adducts exhibit a double-well potential energy surface with two distinct minima: a "long bond" geometry (LB, a van der Waals interaction mostly based on electrostatics, but including a residual charge transfer component) and a "short bond" structure (SB, a covalent dative bond). This behavior can be considered as a "weak" form of bond stretch isomerism. Our computations reveal that complexes RCN-BX₃ (R=CH₃ , FCH₂ , BrCH₂ , and X=Cl, Br) exhibit a fast interconversion from LB to SB geometries even close to the absolute zero thanks to a boron atom tunneling mechanism. The computed half-lives of the meta-stable LB compounds vary between minutes to nanoseconds at cryogenic conditions. Accordingly, we predict that the long bond structures are practically impossible to isolate or characterize, which agrees with previous matrix-isolation experiments.

New insights in chemical reactivity from quantum chemical topology

Based on the quantum chemical topology of the modified electron localization function ELF_x , an efficient and robust mechanistic methodology designed to identify the favorable reaction pathway between two reactants is proposed. We first recall and reshape how the supermolecular interaction energy can be evaluated from only three distinct terms, namely the intermolecular coulomb energy, the intermolecular exchange-correlation energy and the intramolecular energies of reactants. Thereafter, we show that the reactivity between the reactants is driven by the first-order variation in the coulomb intermolecular energy defined in terms of the response to changes in the number of electrons. Illustrative examples with the formation of the dative bond B-N involved in the BH₃ NH₃ molecule and the typical formation of the hydrogen bond in the canonical water dimer are presented. For these selected systems, our approach unveils a noticeable mimicking of E_dual onto the DFT intermolecular interaction energy surface calculated between the both reactants. An automated reaction-path algorithm aimed to determine the most favorable relative orientations when the two molecules approach each other is also outlined.