Unraveling the role of entropy in tuning unimolecular vs. bimolecular reaction rates: The case of olefin polymerization catalyzed by transition metals

Synthesis of pyrrolo[2,3-d]pyridazines and pyrrolo[2,3-c]pyridines (6-azaindoles). Experimental and theoretical study

A new synthetic route towards 6-azaindoles (pyrrolo[2,3-c]pyridines) and pyrrolo[2,3-d]pyridazines starting from 4-aroyl pyrroles is described. This overall protocol involves: (i) the Vilsmeier-Haack reaction to obtain pyrrolo-2,3-dicarbonyles and (ii) condensation with hydrazines or glycine methyl ester. The reaction mechanism between pyrrolo-2,3-dicarbonyl with phenyl hydrazine and glycine methyl ester has been modelled using DFT calculations to prove the formation of one from two possible isomers of condensation.

Predicting the propene stereoselectivity on transition metal catalysts: A daunting task for density functional theory

Thanks to recent developments in hardware and software, quantum chemical methods are increasingly used for interpreting the complex mechanisms underlying polymerization reaction by homogeneous catalysis. Unfortunately, the dimensions of even the smallest realistic models are too large to permit the use of state-of-the-art composite wave function methods. Under these circumstances, density functional theory still offers the best compromise between cost and accuracy. However, comprehensive benchmarks of different functionals are not yet available for this important research field. The main aim of the present paper is to fill this gap by performing an unbiased comparison of several density functionals and continuum solvent models for the stereo-control in the propylene polymerization on prototypical catalysts inducing different reaction mechanisms. While it was not possible to define a unique computational protocol providing the best results in all the situations, the B3PW91 functional in conjunction with D3 empirical dispersions and the solvent model density solvent model performs remarkably well for three out of the four investigated catalysts. Under such circumstances, it is recommended to compare the results delivered by different models when approaching additional classes of catalysts.

Sheet Models for Methylaluminoxane (MAO) Activators? A Theoretical Case Study involving rac-Me2Si(η5-C9H6)2Zr (SBIZr) Complexes

Activation of SBIZrMe2 or SBIZrMeCl and a sheet model for an active component of hydrolytic MAO, (MeAlO)16(Me3Al)6, (16,6) has been studied by DFT. Contact ion-pair formation occurs through the intermediacy of SBIZrMe(Cl) or SBIZrMe2 reacting with sheet 16,6 to furnish SBIZrMe-μ-X(MeAlO)16(Me3Al)6 (2, X=Me, Cl). Contact ion-pairs 2 would be in equilibrium with heterodinuclear catalyst precursors [SBIZrMe2AlMe2][(MeAlO)16(Me3Al)6X] (3 (X=Me, Cl) through reversible binding of Me3Al at higher Al : Zr ratios. Calculations show that formation of ion-pairs 3 from contact ion-pairs 2 is more favourable for the SBIZr compared with the parent Cp2Zr complexes. TD-DFT calculations were conducted on relevant SBIZr complexes to relate the results to earlier spectroscopic studies of catalyst activation using UV-Vis spectroscopy. Finally, propene insertion into ion-pairs 2, SBIZrMe-μ-MeB(C6F5)3 (6) and [SBIZrMe][B(C6F5)4] (7) was studied at M06-2X/TZVP level of theory. These studies suggest that contact ion-pairs 2 are significantly less reactive towards insertion than 6 or 7, in disagreement with experiment.

Revisiting Stereoselective Propene Polymerization Mechanisms: Insights through the Activation Strain Model

The stereoelectronic factors responsible for stereoselectivity in propene polymerization with several metallocene and post-metallocene transition metal catalysts have been revisited using a combined approach of DFT calculations, the Activation Strain Model, Natural Energy Decomposition Analysis and a molecular descriptor (%VBur). There are in most cases two different paths leading to the formation of stereoerrors (SE), and the classical model does not suffice to fully understand stereoregulation. Improving stereoselectivity requires raising the energies of both SE insertion transition states. Our analyses show that the degrees of deformation of the active site (catalyst+chain) and the prochiral monomer differ for these two paths, and between different catalyst classes. Based on such analyses we discuss: a) the subtle differences in SE formation between stereoselective catalysts with different ligand frameworks; b) the reason for exceptional stereoselectivity reported for a special ansa-metallocene catalyst; c) the (double) stereocontrol origin for isoselective catalysts; d) the electronic contribution for isoselective catalysts generating SE by a modification of the ligand wrapping mode during the polymerization. Although this study will not immediately suggest new catalyst structures, we believe that understanding stereoregulation in great detail will increase our chances of success.

Theoretical Insight into the Palladium-Catalyzed Prenylation and Geranylation of Oxindoles with Isoprene

This work presents a comprehensive mechanistic study of the ligand-controlled palladium-catalyzed prenylation (with C5 added) and geranylation (with C10 added) reactions of oxindole with isoprene. The calculated results indicate that the prenylation with the bis-phosphine ligand and geranylation with the monophosphine ligand fundamentally share a common mechanism. This mechanism involves the formation of two crucial species: a η3-allyl-Pd(II) cation and an oxindole carbon anion. Furthermore, the reactions necessitate the assistance of a second oxindole molecule, which serves as a Brønsted acid, providing a proton to generate the oxindole nitrogen anion. The oxindole nitrogen anion then acts as a Brønsted base, abstracting a C-H proton from another oxindole molecule to form an oxindole carbon anion. These mechanistic details differ significantly from those proposed in the experimental work. The present calculations do not support the presence of the Pd-H species and the η3, η3-diallyl-Pd(II) intermediate, which were previously suggested in experiments. The theoretical results rationalize the experimental finding that the bis-phosphine ligand favors the prenylation of oxindole, while the monophosphine ligand enables the geranylation of oxindole.

Mechanism of Alternating Poly(lactic-co-glycolic acid) Formation by Polymerization of (S)- and (R)-3-Methyl Glycolide Using an Enantiopure Aluminum Complex

The mechanism(s) of alternating PLGA synthesis by ring-opening polymerization of (S)- and (R)-3-methyl glycolide promoted by enantiopure aluminum complexes have been rationalized by density functional theory (DFT) calculations. The high regioselectivity of the (S)-MeG polymerization is obtained by repetitive ring opening at the glycolyl site by the (R)-catalyst whereas a lower regioselectivity is predicted by the ROP of (R)-MeG. The behavior of the two monomers is rationalized by unveiling the active site fluxionality of the enantiopure catalyst, identifying the rate-limiting steps that encode a preference at the glycolyl site versus the lactyl site, and revealing selection of the opposite monomer enantioface. The microstructure of the PLGA copolymers is predicted by considering the influence of the configuration of the last inserted unit. The identification of the preferred mechanistic paths may allow for a targeted catalyst design to enhance control of the polymer microstructures.

Self-organization of metallo-supramolecular polymer networks by free formation of pyridine-phenanthroline heteroleptic complexes

Nature employs spontaneous self-organization of supramolecular bonds to create complex matter capable of adaptation and self-healing. Accordingly, the self-sorting of unlike ligands towards a cooperative heteroleptic complex or narcistic homoleptic association in a mixed ligand system is frequently employed to form interchangeable stimuli-responsive complex geometries with a wide range of applications. This notion is however just rarely employed in the organization of polymer networks. In this paper, we report the free-formation of heteroleptic complexes between tetra-am poly(ethylene glycol) (tetraPEG) precursors functionalized either with pyridine (tetraPy) or phenanthroline (tetraEPhen). Among a wide range of studied metal ions, tetraPy could form a network only in combination with Pd2+, presumably with a square-planar geometry, highlighting the importance of complex strength and stability in forming gels with monodentate ligands. Also, mixed networks with tetraEPhen form only in combination with Pd2+ and Fe2+, with strengths surpassing those of individual components and stabilities incomparable to those of parent networks, indicative of heteroleptic complexation. Extensive rheological, UV-vis, and DFT simulation studies revealed the coexistence of different coordination geometries, with an octahedral arrangement prevailing in the presence of Fe2+ and a square-planar geometry in the presence of Pd2+. Therefore, this study offers new opportunities for the development of stimuli-responsive topology-switching polymer networks.

Ring-opening polymerization of cyclic oligosiloxanes without producing cyclic oligomers

A long-standing problem associated with silicone synthesis is contamination of the polymer products with 10 to 15% cyclic oligosiloxanes that results from backbiting reactions at the polymer chain ends. This process, in competition with chain propagation through ring-opening polymerization (ROP) of cyclic monomers, was thought to be unavoidable and routinely leads to a thermodynamically controlled reaction mixture (polymer/cyclic oligosiloxanes = 85/15). Here, we report that simple alcohol coordination to the anionic chain ends prevents the backbiting process and that a well-designed phosphonium cation acts as a self-quenching system in response to loss of coordinating alcohols to stop the reaction before the backbiting process begins. The combination of both effects allows a thermodynamically controlled ROP of the eight-membered siloxane ring D4 without producing undesirable cyclic oligosiloxanes.

H-Bonding leading to latent initiators for olefin metathesis polymerization

Ruthenium-NHC based catalysts, with a chelated iminium ligand trans to the N-heterocyclic carbene (NHC) ligand, that polymerize dicyclopentadiene (DCPD) at different temperatures are monitored using Density Functional Theory calculations to unveil the reaction mechanism, and subsequently how important are the geometrical and electronic features vs. the non-covalent interactions in between. The balance is very fragile and H-bonds are fundamental to explain the different behaviour of latent catalysts. This computational study aims to facilitate future studies of new generations of latent initiators for olefin metathesis polymerization, with the 3D and mainly the 2D Non-Covalent Interaction plots the characterization tool for H-bonds.

A Computational Evaluation of the Steric and Electronic Contributions in Stereoselective Olefin Polymerization with Pyridylamido-Type Catalysts

A density functional theory (DFT) study combined with the steric maps of buried volume (%V_Bur) as molecular descriptors and an energy decomposition analysis through the ASM (activation strain model)-NEDA (natural energy decomposition analysis) approach were applied to investigate the origins of stereoselectivity for propene polymerization promoted by pyridylamido-type nonmetallocene systems. The relationships between the fine tuning of the ligand and the propene stereoregularity were rationalized (e.g., the metallacycle size, chemical nature of the bridge, and substituents at the ortho-position on the aniline moieties). The DFT calculations and %V_Bur steric maps reproduced the experimental trend: substituents on the bridge and on the ortho-positions of aniline fragments enhance the stereoselectivity. The ASM-NEDA analysis enabled the separation of the steric and electronic effects and revealed how subtle ligand modification may affect the stereoselectivity of the process.

DFT calculations in solution systems: solvation energy, dispersion energy and entropy

DFT calculations of reaction mechanisms in solution have always been a hot topic, especially for transition-metal-catalyzed reactions. The calculation of solvation energy is performed using either the polarizable continuum model (PCM) or the universal solvation model SMD. The PCM calculation is very sensitive to the choice of atomic radii to form a cavity, where the self-consistent isodensity PCM (SCI-PCM) has been recognized as the best choice and our IDSCRF radii can provide a similar cavity. Moving from a gas-phase case to a solution case, dispersion energy and entropy should be carefully treated. The solvent-solute dispersion is also important in solution systems, and it should be calculated together with the solute dispersion. Only half of the solvent-solute dispersion energy from the PCM calculation belongs to the solute molecules to maintain a thermal equilibrium between a solute molecule and its cavity, similar to the treatment of electrostatic energy. Relative solute dispersion energy should also be shared equally with the newly formed cavity. The entropy change from a gas phase to a liquid phase is quite large, but the modern quantum chemistry programs can only calculate the gas-phase translational entropy based on the idea-gas equation. In this review, we will provide an operable method to calculate the solution translational entropy, which has been coded in our THERMO program.

Unconventional Stereoerror Formation Mechanisms in Nonmetallocene Propene Polymerization Systems Revealed by DFT Calculations

An unconventional mechanism for the stereoerror formation in propene polymerization catalyzed by C₁-symmetric salalen-M systems (M = Zr, Hf) is suggested by DFT calculations. While propagation happens with the ligand in its fac-mer conformation, a change of ligand wrapping mode from fac-mer to fac-fac is the main source of the lower stereoselectivities obtained with Zr and Hf. This is different for the Ti analogues, where the ligand fac-mer wrapping mode does not play a role. Activation strain analysis indicates that the preference for a chain stationary mechanism of the Zr/Hf species is due to the energy required to distort the reactants (ΔE_Strain) rather than to their mutual interaction (ΔE_Int)_.

Hydroaminoalkylation of alkenes using transition metals complexes grafted on silica SBA15 as catalysts

The gas-phase hydroaminoalkylation reaction of propene catalyzed by group 4 (M = Ti, Zr and Hf) metal amido complexes [(≡Si-O-)(M(-NMe₂)₃] was investigated by using PBE0-D3/SVP//TZVP level of theory. The geometrical analysis traced the formation of the metallaaziridines and the azametallacyclopentanes as key intermediates in these reactions. The metallaaziridines were simulated through the activation of α-C-H bonds of the amido groups; while the azametallacyclopentanes were configured by slotting the propene double bond onto the M - C bonds of the metallaaziridines. The latter reaction was considered the rate-determining step. Thermochemical calculations showed that the order of catalytic activity is: Ti ≥ Zr > Hf; while the preference of the azametallacyclopentanes is: Hf > Zr ≥ Ti.

Thermodynamics of metallocene catalyst activation: alignment of theory and experiment

Three equilibria involved in metallocene catalyst activation, including dissociation of R₆Al₂ (R = Me, Et or i-Bu) and related species such as [L₂ZrMe₂AlMe₂][B(C₆F₅)₄] (L₂ = Cp₂, 1,2-ethylenebis(η⁵-indenyl), Me₂C(η⁵-C₅H₄)₂) or [(L₂ZrMe)₂μ-Me][MePBB] (L₂ = (h⁵-1,2-Me₂C₅H₃)₂, [MePBB]^- = [MeB(Ar_F)₃]^- with Ar_F = o-C₆F₅-C₆F₄) are studied by DFT using various approaches to account for the enthalpy and entropy changes in gas and condensed phases. These studies reveal that both low energy vibrations and translational entropy conspire to cause significant deviations between theory and experiment when it comes to the free energy change in condensed or even gas phase. Alignment of theory with experiment requires in addition, consideration of specific solvation of reactants and products.

Ring-Opening Polymerization of L-Lactide Catalyzed by Potassium-Based Complexes: Mechanistic Studies

Two non-toxic potassium compounds, 1 and 2, with a commercial oximate ligand have been prepared and fully spectroscopically characterized. Their activity as catalysts for the ring-opening polymerization (ROP) process of LLA has been studied, showing that they are extremely active and able to polymerize the monomer in a few minutes. For derivative 2, the presence of a crown ether in the potassium coordination sphere affects the nuclearity of the compound and consequently its solubility, with both aspects having an influence in the polymerization process. Detailed studies of the polymerization mechanism have been performed, and an unusual anionic mechanism was observed in absence of a co-initiator. Indeed, the monomer deprotonation generates a lactide enolate, which initiates the polymerization propagation. On the contrary, when a 1:1 ratio of cat:BnOH is used, a mixture of mechanisms is observed, the anionic mechanism and the activated monomer one, while from a cat:BnOH ratio of 1:2 and over, only the activated monomer mechanism is observed.

Non-competitive interactions between hydroxychloroquine and azithromycin: Systematic density functional, molecular dynamics, and docking calculations

In this study, density functional theory (DFT) and docking calculations were systematically performed to study the non-competitive interaction between Hydroxychloroquine (HCQ) and azithromycin (AZTH). The calculated changes in Gibbs free energy and enthalpy (at 310 K) were positive, indicating the non-spontaneous formation of HCQ-AZTH specifically in water media. Docking calculation confirmed the obtained DFT result as evident from the different binding sites of both drugs to the SARS-CoV-2 main protease and human angiotensin-converting enzyme 2 (ACE2) proteins. The HCQ-AZTH structure revealed enhanced electrochemical properties, suggesting the synergy between HCQ and AZTH without affecting their therapeutic efficacy against SARS-CoV-2.

Molecular and supported Ti(iii)-alkyls: efficient ethylene polymerization driven by the π-character of metal-carbon bonds and back donation from a singly occupied molecular orbital

While Ti(iii) alkyl species are the proposed active sites in Ziegler–Natta ethylene polymerization catalysts, the corresponding well-defined homogeneous catalysts are not known. We report that well-defined neutral β-diiminato Ti(iii) alkyl species, namely [Ti(nacnac)(CH₂^tBu)₂] and its alumina-grafted derivative [(Al_sO)Ti(nacnac)(CH₂^tBu)], are active towards ethylene polymerization at moderate pressures and temperatures and possess an electron configuration well-adapted to insertion of ethylene. Advanced EPR spectroscopy showed that ethylene insertion into a Ti(iii)–C bond takes place during polymerization from Ti(nacnac)(CH₂^tBu)₂. A combination of pulsed EPR spectroscopy and DFT calculations, based on a crystal structure of [Ti(nacnac)(CH₂^tBu)₂], enabled us to reveal details about the structure and electronic configurations of both molecular and surface-grafted species. For both compounds, the α-agostic C–H interaction, which involves the singly occupied molecular orbital, indicates a π character of the metal–carbon bond; this π character is enhanced upon ethylene coordination, leading to a nearly barrier-less C₂H₄ insertion into Ti(iii)–C bonds after this first step. During coordination, back donation from the SOMO to the π*(C₂H₄) occurs, leading to stabilization of π-ethylene complexes and to a significant lowering of the overall energy of the C₂H₄ insertion transition state. In d¹ alkyl complexes, ethylene insertion follows an original “augmented” Cossee–Arlman mechanism that involves the delocalization of unpaired electrons between the SOMO, π*(C₂H₄) and σ*(Ti–C) in the transition state, which further favors ethylene insertion. All these factors facilitate ethylene polymerization on Ti(iii) neutral alkyl species and make d¹ alkyl complexes potentially more effective polymerization catalysts than their d⁰ analogues.

Ti(iii) alkyl species polymerize ethylene via an original mechanism, which involves back donation to the π*(C₂H₄) and a delocalization of the unpaired electron in the transition state of C₂H₄ insertion into the partially alkylidenic Ti(iii)–C bond.

N-Heterocyclic Carbene Gold(I) Complexes: Mechanism of the Ligand Scrambling Reaction and Their Oxidation to Gold(III) in Aqueous Solutions

N-Heterocyclic carbene (NHC) gold(I) complexes offer great prospects in medicinal chemistry as antiproliferative, anticancer, and antibacterial agents. However, further development requires a thorough understanding of their reaction behavior in aqueous media. Herein, we report the conversion of the bromido[3-ethyl-4-(4-methoxyphenyl)-5-(2-methoxypyridin-5-yl)-1-propylimidazol-2-ylidene]gold(I) ((NHC)Au^IBr, 1) complex in acetonitrile/water mixtures to the bis[3-ethyl-4-(4-methoxyphenyl)-5-(2-methoxypyridin-5-yl)-1-propylimidazol-2-ylidene]gold(I) ([(NHC)₂Au^I]⁺, 7), which is subsequently oxidized to the dibromidobis[3-ethyl-4-(4-methoxyphenyl)-5-(2-methoxypyridin-5-yl)-1-propylimidazol-2-ylidene]gold(III) ([(NHC)₂Au^IIIBr₂]⁺, 9). By combining experimental data from HPLC, NMR, and (LC-)/HR-MS with computational results from DFT calculations, we outline a detailed ligand scrambling reaction mechanism. The key step is the formation of the stacked ((NHC)Au^IBr)₂ dimer (2) that rearranges to the T-shaped intermediate Br(NHC)₂Au^I–Au^IBr (3). The dissociation of Br^– from 3 and recombination lead to (NHC)₂Au^I–Au^IBr₂ (5) followed by the separation into [(NHC)₂Au^I]⁺ (7) and [Au^IBr₂]⁻ (8). [Au^IBr₂]⁻ is not stable in an aqueous environment and degrades in an internal redox reaction to Au⁰ and Br₂. The latter in turn oxidizes 7 to the gold(III) species 9. The reported ligand rearrangement of the (NHC)Au^IBr complex differs from that found for related silver(I) analogous. A detailed understanding of this scrambling mechanism is of utmost importance for the interpretation of their biological activity and will help to further optimize them for biomedical and other applications.

By means of experimental data from HPLC and (LC-)MS in combination with DFT calculations, we present a detailed mechanism for the ligand scrambling reaction of (NHC)Au^IBr to the corresponding [(NHC)₂Au^I]⁺ complex and the oxidation to the [(NHC)₂Au^IIIBr₂]⁺ species in aqueous solutions.