Cyclotrimerization Reactions Catalyzed by Rhodium(I) Half-Sandwich Complexes: A Mechanistic Density Functional Study

Designing Rh(I)-Half-Sandwich Catalysts for Alkyne [2+2+2] Cycloadditions

Metal-mediated [2+2+2] cycloadditions of unsaturated molecules to cyclic and polycyclic organic compounds are a versatile synthetic route affording good yields and selectivity under mild conditions. In the last two decades, in silico investigations have unveiled important details about the mechanism and the energetics of the whole catalytic cycle. Particularly, a number of computational studies address the topic of half-sandwich catalysts which, due to their structural fluxionality, have been widely employed, since the 1980s. In these organometallic species, the metal is coordinated to an aromatic ring, typically the ubiquitous cyclopentadienyl anion, C5H5 –(Cp) or to the Cp moiety of a larger polycyclic aromatic ligand (Cp′). During the catalytic process, the metal continuously ‘slips’ on the ring, changing its hapticity. This phenomenon of metal slippage and its implications for the catalyst’s performance are discussed in this work, referring to the most important computational mechanistic studies reported in literature for Rh(I) half-metallocenes, with the purpose of providing hints for a rational design of this class of compounds.

1 Introduction

2 Mechanism of Metal-Catalyzed Acetylene [2+2+2] Cycloaddition to Benzene and the Problem of the Indenyl Effect

2.1 Acetylene-Acetonitrile [2+2+2] Co-cycloaddition to 2-Methylpyridine: Evidence of the Indenyl Effect

2.2 Heteroaromatic Catalysts and the Evidence of a Reverse Indenyl Effect

2.3 Booth’s Mechanistic Hypothesis and the Evidence of the Indenyl Effect

3 Structure–Reactivity Correlation: The Slippage-Span Model

4 Conclusions and Perspectives

Computational prediction on the catalytic activity of heterobimetallic complex featuring MM' triple bond in acetylene cyclotrimerization: Mechanistic study

A detailed reaction mechanism of acetylene cyclotrimerization catalyzed by V(ⁱ PrNPMe₂ )₃ Fe-PMe₃ (denote as CAT), a heterobimetallic complex featuring V-Fe triple bond, was computationally investigated using density functional theory. The calculated results show that the first acetylene firstly attaches to the V atom of CAT to get a four-membered ring structure through [2 + 2] cycloaddition reaction. For the second acetylene addition, there are two cyclotrimerization mechanisms, outer sphere mechanism and inner mechanism. The inner sphere reaction pathway is the main reaction pathway. By replacing the V with Nb and Ta, Fe with Ru and Os, a series of new catalysts are screened computationally. The calculated results show that, all of the nine heterobimetallic complexes show high activity at mild condition. The energy barrier of the rate determining step is related to the natural population analysis (NPA) charge of M' and the Wiberg bond index (WBI) of M-M' bond. The more negative NPA charge of M' and the smaller WBI of M-M' bond, the lower energy barrier is.

Highly Stable Silver(I) Complexes with Cyclen-Based Ligands Bearing Sulfide Arms: A Step Toward Silver-111 Labeled Radiopharmaceuticals

With a half-life of 7.45 days, silver-111 (β_max 1.04 MeV, E_γ 245.4 keV [I_γ 1.24%], E_γ 342.1 keV [I_γ 6.7%]) is a promising candidate for targeted cancer therapy with β^- emitters as well as for associated SPECT imaging. For its clinical use, the development of suitable ligands that form sufficiently stable Ag⁺-complexes in vivo is required. In this work, the following sulfur-containing derivatives of tetraazacyclododecane (cyclen) have been considered as potential chelators for silver-111: 1,4,7,10-tetrakis(2-(methylsulfanyl)ethyl)-1,4,7,10-tetraazacyclododecane (DO4S), (2S,5S,8S,11S)-2,5,8,11-tetramethyl-1,4,7,10-tetrakis(2-(methylsulfanyl)ethyl)-1,4,7,10-tetraazacyclododecane (DO4S4Me), 1,4,7-tris(2-(methylsulfanyl)ethyl)-1,4,7,10-tetraazacyclododecane (DO3S), 1,4,7-tris(2-(methylsulfanyl)ethyl)-10-acetamido-1,4,7,10-tetraazacyclododecane (DO3SAm), and 1,7-bis(2-(methylsulfanyl)ethyl)-4,10,diacetic acid-1,4,7,10-tetraazacyclododecane (DO2A2S). Natural Ag⁺ was used in pH/Ag-potentiometric and UV-vis spectrophotometric studies to determine the metal speciation existing in aqueous NaNO₃ 0.15 M at 25 °C and the equilibrium constants of the complexes, whereas NMR and DFT calculations gave structural insights. Overall results indicated that sulfide pendant arms coordinate Ag⁺ allowing the formation of very stable complexes, both at acidic and physiological pH. Furthermore, radiolabeling, stability in saline phosphate buffer, and metal-competition experiments using the two ligands forming the strongest complexes, DO4S and DO4S4Me, were carried out with [¹¹¹Ag]Ag⁺ and promising results were obtained.

Discerning the thermal cyclotrimerizations of fluoro- and chloroacetylenes through ELF, NBO descriptors and QTAIM analysis: pseudodiradical character

In this study the thermal cyclotrimerization reactions of fluoro- and chloroacetylenes involving regioselectively stepwise {2 + 2} and stepwise {4 + 2} cycloadditions were studied using the topological analysis of the electron localization function (ELF), the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses. These methodologies have shown that the electronic reorganization in the regioselectively stepwise {2 + 2} and stepwise {4 + 2} cycloadditions may be considered as {2n+2n} and {2π+2n} pseudodiradical process, respectively. Finally, the last phase of this thermal reaction can be understood as an electronic migration process under the pseudodiradical character in the thermal ring-opening reaction, with the subsequent formation of reaction products. In this sense, new insights are reported on the electronic behavior in the bond formation in the thermal cyclotrimerization of fluoroacetylene.

Diverse ring-opening reactions of rhodium η4-azaborete complexes

Sequential treatment of [Rh(COE)₂Cl]₂ (COE = cyclooctene) with PⁱPr₃, alkyne derivatives and ^tBuN Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> BMes (Mes = 2,4,6-trimethylphenyl) provided functionalized rhodium η⁴-1,2-azaborete complexes of the form (η⁴-azaborete)RhCl(PⁱPr₃). The scope of this reaction was expanded to encompass alkynes with hydrogen, alkyl, aryl, ferrocenyl, alkynyl, azaborinyl and boronate ester substituents. Treatment of these complexes with PMe₃ led to insertion of the rhodium atom into the B–C bond of the BNC₂ ring, forming 1-rhoda-3,2-azaboroles. Addition of N-heterocyclic carbenes to azaborete complexes led to highly unusual rearrangements to rhodium η²,κ¹-allenylborylamino complexes via deprotonation and C–N bond cleavage. Heating and photolysis of an azaborete complex also led to rupture of the C–N bond followed by subsequent rearrangements, yielding an η⁴-aminoborylallene complex and two isomeric η⁴-butadiene complexes.

Rhodium η⁴-azaborete complexes can be transformed into a variety of species with ring-opened, BN-containing ligands by treatment with Lewis bases.

In Silico Acetylene [2+2+2] Cycloadditions Catalyzed by Rh/Cr Indenyl Fragments

Metal-catalyzed alkyne [2+2+2] cycloadditions provide a variety of substantial aromatic compounds of interest in the chemical and pharmaceutical industries. Herein, the mechanistic aspects of the acetylene [2+2+2] cycloaddition mediated by bimetallic half-sandwich catalysts [Cr(CO)3IndRh] (Ind = (C9H7)−, indenyl anion) are investigated. A detailed exploration of the potential energy surfaces (PESs) was carried out to identify the intermediates and transition states, using a relativistic density functional theory (DFT) approach. For comparison, monometallic parent systems, i.e., CpRh (Cp = (C5H5)−, cyclopentadienyl anion) and IndRh, were included in the analysis. The active center is the rhodium nucleus, where the [2+2+2] cycloaddition occurs. The coordination of the Cr(CO)3 group, which may be in syn or anti conformation, affects the energetics of the catalytic cycle as well as the mechanism. The reaction and activation energies and the turnover frequency (TOF) of the catalytic cycles are rationalized, and, in agreement with the experimental findings, our computational analysis reveals that the presence of the second metal favors the catalysis.

Enhanced Catalytic Activity of Nickel Complexes of an Adaptive Diphosphine-Benzophenone Ligand in Alkyne Cyclotrimerization

Adaptive ligands, which can adapt their coordination mode to the electronic structure of various catalytic intermediates, offer the potential to develop improved homogeneous catalysts in terms of activity and selectivity. 2,2'-Diphosphinobenzophenones have previously been shown to act as adaptive ligands, the central ketone moiety preferentially coordinating reduced metal centers. Herein, the utility of this scaffold in nickel-catalyzed alkyne cyclotrimerization is investigated. The complex [( p-tolL1)Ni(BPI)] ( p-tolL1 = 2,2'-bis(di(para-tolyl)phosphino)-benzophenone; BPI = benzophenone imine) is an active catalyst in the [2 + 2 + 2] cyclotrimerization of terminal alkynes, selectively affording 1,2,4-substituted benzenes from terminal alkynes. In particular, this catalyst outperforms closely related bi- and tridentate phosphine-based Ni catalysts. This suggests a reaction pathway involving a hemilabile interaction of the C=O unit with the nickel center. This is further borne out by a comparative study of the observed resting states and DFT calculations.

Half-Sandwich Metal-Catalyzed Alkyne [2+2+2] Cycloadditions and the Slippage Span Model

Half-sandwich RhI compounds display good catalytic activity toward alkyne [2+2+2] cycloadditions. A peculiar structural feature of these catalysts is the coordination of the metal to an aromatic moiety, typically a cyclopentadienyl anion, and, in particular, the possibility to change the bonding mode easily by the metal slipping over this aromatic moiety. Upon modifying the ancillary ligands, or proceeding along the catalytic cycle, hapticity changes can be observed; it varies from η5, if the five metal-carbon distances are identical, through η3+η2, in the presence of allylic distortion, and η3, in the case of allylic coordination, to η1, if a σ metal-carbon bond forms. In this study, we present the slippage span model, derived with the aim of establishing a relationship between slippage variation during the catalytic cycle, quantified in a novel and rigorous way, and the performance of catalysts in terms of turnover frequency, computed with the energy span model. By collecting and comparing new data and data from the literature, we find that the highest performance is associated with the smallest slippage variation along the cycle.

Group 9 Metallacyclopentadienes as Key Intermediates in [2+2+2] Alkyne Cyclotrimerizations. Insight from Activation Strain Analyses

The intramolecular oxidative coupling converting a bis-acetylene complex of formula CpM (C₂ H₂ )₂ (Cp=C₅ H₅^- ; M=Co, Rh, Ir) into a 16-electron metallacycle is studied in silico. This reaction is paradigmatic in acetylene [2+2+2] cycloaddition to benzene catalyzed by CpM fragments, being the step with the highest activation energy, and thus affecting the whole catalysis. Our activation strain and quantitative molecular orbital (MO) analyses elucidate the mechanistic details and reveal why cobalt performs better than rhodium and iridium catalysts outlining general principles for rational design of catalysts to be used in these processes.

A Computational Study of the Intermolecular [2+2+2] Cycloaddition of Acetylene and C60 Catalyzed by Wilkinson's Catalyst

The functionalization of fullerenes helps to modulate their electronic and physicochemical properties, generating fullerene derivatives with promising features for practical applications. Herein, DFT is used to explore the attachment of a cyclohexadiene ring to C₆₀ through a rhodium-catalyzed intermolecular [2+2+2] cycloaddition of C₆₀ and acetylene. All potential reaction paths are analyzed and it can be concluded that the [2+2+2] cycloaddition of C₆₀ and two acetylene molecules catalyzed by [RhCl(PPh₃ )₃ ], yielding a cyclohexadiene ring fused to a [6,6] bond of C₆₀ , is energetically feasible.

Reaction of [TpRh(C2 H4 )2 ] with Dimethyl Acetylenedicarboxylate: Identification of Intermediates of the [2+2+2] Alkyne and Alkyne-Ethylene Cyclo(co)trimerizations

The reaction between the bis(ethylene) complex [TpRh(C2 H4 )2 ], 1, (Tp=hydrotris(pyrazolyl)borate), and dimethyl acetylenedicarboxylate (DMAD) has been studied under different experimental conditions. A mixture of products was formed, in which TpRh(I) species were prevalent, whereas the presence of trapping agents, like water or acetonitrile, allowed for the stabilization and isolation of octahedral TpRh(III) compounds. An excess of DMAD gave rise to a small amount of the [2+2+2] cyclotrimerization product hexamethyl mellitate (6). Although no catalytic application of 1 was achieved, mechanistic insights shed light on the formation of stable rhodium species representing the resting state of the catalytic cycle of rhodium-mediated [2+2+2] cyclo(co)trimerization reactions. Metallacyclopentene intermediate species, generated from the activation of one alkyne and one ethylene molecule from 1, and metallacyclopentadiene species, formed by oxidative coupling of two alkynes to the rhodium centre, are crucial steps in the pathways leading to the final organometallic and organic products.

Synthesis of Functionalized 1,4-Azaborinines by the Cyclization of Di-tert-butyliminoborane and Alkynes

Di-tert-butyliminoborane is found to be a very useful synthon for the synthesis of a variety of functionalized 1,4-azaborinines by the Rh-mediated cyclization of iminoboranes with alkynes. The reactions proceed via [2 + 2] cycloaddition of iminoboranes and alkynes in the presence of [RhCl(PiPr3)2]2, which gives a rhodium η(4)-1,2-azaborete complex that yields 1,4-azaborinines upon reaction with acetylene. This reaction is compatible with substrates containing more than one alkynyl unit, cleanly affording compounds containing multiple 1,4-azaborinines. The substitution of terminal alkynes for acetylene also led to 1,4-azaborinines, enabling ring substitution at a predetermined location. We report the first general synthesis of this new methodology, which provides highly regioselective access to valuable 1,4-azaborinines in moderate yields. A mechanistic rationale for this reaction is supported by DFT calculations, which show the observed regioselectivity to arise from steric effects in the B-C bond coupling en route to the rhodium η(4)-1,2-azaborete complex and the selective oxidative cleavage of the B-N bond of the 1,2-azaborete ligand in its subsequent reaction with acetylene.

The origin of the ligand-controlled regioselectivity in Rh-catalyzed [(2 + 2) + 2] carbocyclizations: steric vs. stereoelectronic effects

Density functional theory calculations demonstrate that the reversal of regiochemical outcome of the addition for substituted methyl propiolates in the rhodium-catalyzed [(2 + 2) + 2] carbocyclization with PPh₃ and (S)-xyl-binap as ligands is both electronically and sterically controlled.

Density functional theory calculations demonstrate that the reversal of regiochemical outcome of the addition for substituted methyl propiolates in the rhodium-catalyzed [(2 + 2) + 2] carbocyclization with PPh₃ and (S)-xyl-binap as ligands is both electronically and sterically controlled. For example, the ester functionality polarizes the alkyne π* orbital to favor overlap of the methyl-substituted terminus of the alkyne with the p_π-orbital of the alkenyl fragment of the rhodacycle during alkyne insertion with PPh₃ as the ligand. In contrast, the sterically demanding xyl-binap ligand cannot accommodate the analogous alkyne orientation, thereby forcing insertion to occur at the sterically preferred ester terminus, overriding the electronically preferred orientation for alkyne insertion.

The activation strain model and molecular orbital theory

The activation strain model is a powerful tool for understanding reactivity, or inertness, of molecular species. This is done by relating the relative energy of a molecular complex along the reaction energy profile to the structural rigidity of the reactants and the strength of their mutual interactions: ΔE(ζ) = ΔE_strain(ζ) + ΔE_int(ζ). We provide a detailed discussion of the model, and elaborate on its strong connection with molecular orbital theory. Using these approaches, a causal relationship is revealed between the properties of the reactants and their reactivity, e.g., reaction barriers and plausible reaction mechanisms. This methodology may reveal intriguing parallels between completely different types of chemical transformations. Thus, the activation strain model constitutes a unifying framework that furthers the development of cross-disciplinary concepts throughout various fields of chemistry. We illustrate the activation strain model in action with selected examples from literature. These examples demonstrate how the methodology is applied to different research questions, how results are interpreted, and how insights into one chemical phenomenon can lead to an improved understanding of another, seemingly completely different chemical process. WIREs Comput Mol Sci 2015, 5:324-343. doi: 10.1002/wcms.1221.

Indenyl effect due to metal slippage? Computational exploration of rhodium-catalyzed acetylene [2+2+2] cyclotrimerization

The mechanism of CpRh (Cp=cyclopentadienyl) and IndRh (Ind=indenyl)-catalyzed acetylene [2+2+2] cyclotrimerization has been revisited aiming at finding an explanation for the better performance of the latter catalyst found experimentally. The hypothesis that an ancillary ligand of the precatalyst remains bonded to the metal center throughout the whole catalytic cycle, based on the experimental evidence that the nature of this ligand can exert some control in cocyclotrimerization of different alkynes, is considered. Strong hapticity variations occur in both the CpRh- and IndRh-catalyzed processes. As the Ind ligand undergoes a more facile slippage than Cp, the energy profile is far smoother in the IndRh-catalyzed cyclotrimerization. This difference in the energetics of the process translates into an enhanced activity of the IndRh catalyst, in nice agreement with experiment.

In silico design of heteroaromatic half-sandwich RhI catalysts for acetylene [2+2+2] cyclotrimerization: evidence of a reverse indenyl effect

A mechanistic density functional theory study of acetylene [2+2+2] cyclotrimerization to benzene catalyzed by Rh(I) half metallocenes is presented. The catalyst fragment contains a heteroaromatic ligand, that is, the 1,2-azaborolyl (Ab) or the 3a,7a-azaborindenyl (Abi) anions, which are isostructural and isoelectronic to the hydrocarbon cyclopentadienyl (Cp) and indenyl (Ind) anions, respectively, but differ from the last ones on having two adjacent carbon atoms replaced with a boron and a nitrogen atom. The better performance of either the classic hydrocarbon or the heteroaromatic catalysts is found to depend on the different mechanistic paths that can be envisioned for the process. The present analyses uncover and explain general structure-reactivity relationships that may serve as rational design principles. In particular, we provide evidence of a reverse indenyl effect.