DFT Studies on the Mechanisms of Carboamination/Diamination of Unactivated Alkenes Mediated by Pd(IV) Intermediates.

Density functional theory (DFT) calculations have been employed to investigate the mechanism of carboamination and diamination of unactivated alkenes mediated by Pd(IV) intermediates. Both reactions share a common Pd(IV) intermediate, serving as the starting point for either the carboamination or the diamination pathway. The formation of this Pd(IV) intermediate encompasses a transition state that substantially impacts the turnover frequency (TOF) of catalytic cycles, with an apparent activation free-energy barrier of 26.1 kcal mol-1. Carboamination of unactivated alkenes proceeds through the coordination of a toluene molecule, C-H activation, inner reductive elimination, and the separation of the carboamination product from this intermediate, while diamination of unactivated alkenes involves the formation of the ion nucleophile, SN2 attack, and the separation of the diamination product. A comparison of the free-energy profiles for carboamination and diamination of unactivated alkenes can elucidate the origin of the chemoselectivity, and Bader's atoms in molecules (AIM) wave function analyses have been performed to analyze the contributions of the outer C-N bonding in the diamination process.

The treatment of dispersion terms for solution systems.

DFT calculations of reaction mechanisms in solution have always been a hot topic, especially for transition-metal-catalyzed reactions, in which the traditional DFT-D3 method has been extensively employed. The overestimation of the dispersion from the traditional DFT-D3 method leads to a quite low activation free-energy barrier, so it is worth finding a proper way to deal with the dispersion for solution systems. The solvent-solute dispersion is also important for solution systems, and thus it should be calculated together with the solute dispersion. The newly generated solute-solute dispersion energy should be shared equally with the newly formed cavity between two interacting species; therefore, only half of the solute-solute and solvent-solute dispersion terms belong to the solute molecule. The detailed treatment of dispersion correction for solution systems has been fully addressed, and this method has been confirmed with the examples of ligand exchange reactions and catalytic reactions.

High-efficient capture and degradation of formaldehyde based on the electric-field-enhanced catalytic effect.

Enhancing the generation of active groups is of great significance for alleviating the catalyst deactivation of formaldehyde (HCHO) by accelerating the decomposition of intermediate products. Herein, an electric-field-enhanced catalytic effect was proposed for the efficient capture and degradation of HCHO base on carbon cloth loaded manganese oxide catalyst (MnOx-CC). Under the action of electric field, MnOx can generate more hydroxyl radicals (•OH) and superoxide radicals (•O2-), thus accelerating the degradation of HCHO and intermediates at room temperature. After the introduction electric field (∼1 ×104 V/m), •O2- and •OH radical on the surface of MnOx-CC catalyst can be increased by 8 times and 23 times, respectively. At weight hourly space velocity of 300,000 mL/(gcat h) for ∼15 ppm HCHO, MnOx-CC-Electric Field catalyst reached the removal efficiency of 99.4%, and the CO2 conversion efficiency of 81.2%, without decrease significantly within 80 h. Theoretical calculation shows that the electric field can increase the electron state density of Mn atom at the Fermi level and reduce the adsorption energy of HCHO, O2 and H2O, thus promoting the generation of active groups and degradation of intermediate products. The electric-field-enhancement catalytic effect provides a new approach for the degradation of Volatile Organic Compounds.

α-Stereoselective 3-Deoxy-d-manno-oct-2-ulosonoic Acid (Kdo) O-Glycosylation with a p-Toluenethioglycoside Donor by the (p-Tol)2SO/Tf2O Preactivation Strategy.

A convenient and efficient approach was developed to synthesize α-Kdo O-glycosides based on the Tf2O/(p-Tol)2SO preactivation strategy using peracetylated Kdo thioglycoside as a donor. Under the optimized reaction conditions, several O-glycoside products, including α-(2 → 1)-, α-(2 → 2)-, α-(2 → 3)-, and α-(2 → 6)-Kdo products, were stereoselectively synthesized in high yields. Remarkably, a series of aromatic α-Kdo O-glycosides were first and successfully constructed in high yields. An SN2-like mechanism was revealed by DFT calculations and experimental results.

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.

DFT Case Study on the Comparison of Ruthenium-Catalyzed C-H Allylation, C-H Alkenylation, and Hydroarylation.

Density functional calculations at the B3LYP-D3+IDSCRF/TZP-DKH(-dfg) level of theory have been performed to understand the mechanism of ruthenium-catalyzed C-H allylation reported in the literature in depth. The plausible pathway consisted of four sequential processes, including C-H activation, migratory insertion, amide extrusion, and recovery of the catalyst, in which C-H activation was identified as the rate-determining step. The amide extrusion step could be promoted kinetically by trifluoroacetic acid since its mediation lowered the free-energy barrier from 32.1 to 12.2 kcal/mol. Additional calculations have been performed to explore other common pathways between arenes and alkenes, such as C-H alkenylation and hydroarylation. A comparison of the amide extrusion and ß-H elimination steps established the following reactivity sequence of the leaving groups: protonated amide group > ß-H group > unprotonated amide group. The suppression of hydroarylation was attributed to the sluggishness of the Ru-C protonation step as compared to the amide extrusion step. This study can unveil factors favoring the C-H allylation reaction.

Tunnelling assisted hydrogen elimination mechanisms of FeCl3/TEMPO.

Metal-TEMPO hybrids are a family of novel and promising catalysts for aerobic oxidation of alcohols, yet the underlying mechanisms have not been understood theoretically. Using density functional theory, we probe the hydrogen abstraction mechanisms of FeCl3/TEMPO on two characteristic substrates, 9,10-dihydroanthracene and benzyl alcohol. We found that the low spin state of FeCl3/TEMPO is more favourable, and that the N atom is the preferred hydrogen acceptor. Moreover, dispersion interactions assist the reaction, as well as nuclear tunnelling, which even at room temperature can speed up the process by almost two orders of magnitude. We also predict that pronounced kinetic isotope effects (KIE) could be observed due to tunnelling. Our findings provide insights into improving the substrate scope and the development of new transformations for the FeCl3/TEMPO system.

Dissociative electron transfer of copper(ii) complexes of glycyl(glycyl/alanyl)tryptophan in vacuo: IRMPD action spectroscopy provides evidence of transition from zwitterionic to non-zwitterionic peptide structures.

We report herein the first detailed study of the mechanism of redox reactions occurring during the gas-phase dissociative electron transfer of prototypical ternary [CuII(dien)M]˙2+ complexes (M, peptide). The two final products are (i) the oxidized non-zwitterionic π-centered [M]˙+ species with both the charge and spin densities delocalized over the indole ring of the tryptophan residue and with a C-terminal COOH group intact, and (ii) the complementary ion [CuI(dien)]+. Infrared multiple photon dissociation (IRMPD) action spectroscopy and low-energy collision-induced dissociation (CID) experiments, in conjunction with density functional theory (DFT) calculations, revealed the structural details of the mass-isolated precursor and product cations. Our experimental and theoretical results indicate that the doubly positively charged precursor [CuII(dien)M]˙2+ features electrostatic coordination through the anionic carboxylate end of the zwitterionic M moiety. An additional interaction exists between the indole ring of the tryptophan residue and one of the primary amino groups of the dien ligand; the DFT calculations provided the structures of the precursor ion, intermediates, and products, and enabled us to keep track of the locations of the charge and unpaired electron. The dissociative one-electron transfer reaction is initiated by a gradual transition of the M tripeptide from the zwitterionic form in [CuII(dien)M]˙2+ to the non-zwitterionic M intermediate, through a cascade of conformational changes and proton transfers. In the next step, the highest energy intermediate is formed; here, the copper center is 5-coordinate with coordination from both the carboxylic acid group and the indole ring. A subsequent switch back to 4-coordination to an intermediate IM1, where attachment to GGW occurs through the indole ring only, creates the structure that ultimately undergoes dissociation.

Synthesis of novel technetium-99m tricarbonyl-HBED-CC complexes and structural prediction in solution by density functional theory calculation.

HBED-CC (N,N'-bis[2-hydroxy-5-(carboxyethyl)benzyl]ethylene diamine-N,N'-diacetic acid, L1 ) is a common bifunctional chelating agent in preparation of 68Ga-radiopharmaceuticals. Due to its high stability constant for the Ga3+ complex (logKGaL = 38.5) and its acyclic structure, it is well known for a rapid and efficient radiolabelling at ambient temperature with Gallium-68 and its high in vivo stability. [99mTc][Tc(CO)3(H2O)3]+ is an excellent precursor for radiolabelling of biomolecules. The aim of this study was to develop a novel preparation method of 99mTc-HBED-CC complexes. In this study, HBED-CC-NI (2,2'-(ethane-1,2-diylbis((2-hydroxy-5-(3-((2-(2-nitro-1H-imidazol-1-yl)ethyl)amino)-3-oxopropyl)benzyl)-azanediyl))-diacetic acid, L2 ), a derivative of HBED-CC, was designed and synthesized. Both L1 and L2 were radiolabelled by [99mTc][Tc(CO)3(H2O)3]+ successfully for the first time. In order to explore the coordination mode of metal and chelates, non-radioactive Re(CO)3 L1 and Re(CO)3 L2 were synthesized and characterized spectroscopically. Tc(CO)3 L1 and Tc(CO)3 L2 in solution were calculated by density functional theory and were analysed with radio-HPLC chromatograms. It showed that [99mTc]Tc(CO)3 L2 forms two stable diastereomers in solution, which is similar to those of [68Ga]Ga-HBED-CC complexes. Natural bond orbital analysis through the natural population charges revealed a charge transfer between [99mTc][Tc(CO)3]+ and L1 or L2 . The experimental results showed that tricarbonyl technetium might form stable complex with HBED-CC derivatives, which is useful for the future application of using HBED-CC as a bifunctional chelating agent in developing new 99mTc-radiopharmaceuticals as diagnostic imaging agents.

Electronic Effect-Guided, Palladium-Catalyzed Regioselective B-H Activation and Multistep Diarylation of o-Carboranes with Aryl Iodides.

Density functional theory calculations at IDSCRF-B3LYP/DZVP computational level were conducted on palladium-catalyzed regioselective B-H activation and diarylation of o-carboranes with aryl iodides in solution. Computational results indicate that this reaction follows a multistep mechanism and needs to get over several transition states before the final B(4,5)-diarylated o-carborane derivatives are formed. B-H activation, oxidation addition, and successive reduction of the Pd(II) catalyst involving a Pd(II)-Pd(IV)-Pd(II) catalytic cycle has been confirmed, in which AgOAc plays a crucial role. Electron-donating group on the cage carbon of o-carboranes is verified to be beneficial for its B-H activation and diarylation, while steric hindrance between the aryl and o-carboranyl groups retards it. Natural population analysis and Gibbs free energetic results predict consistent regioselectivities with experiments and manifest the pivotal role of electronic effect in controlling regioselective B-H activation of o-carboranes. These results are expected to shed some light on further improvement of experimental conditions and better controlling of regioselectivities.

Mechanistic Insights into the Ru(II)-Catalyzed Intramolecular Formal [3 + 2] Cycloaddition of (E)-1,6-Enynes.

Design of a unique reaction pathway in transition-metal-catalyzed 1,6-enynes cyclization to construct valuable synthetic motifs is a significant challenge in organic chemistry. Herein, we report a Ru(II)-catalyzed formal [3 + 2] cycloaddition as an efficient method to prepare unprecedented bicyclo[3.3.0]octenes from readily available (E)-1,6-enynes. Mechanistic studies based on the deuterium labeling experiments and the DFT calculation disclose a reasonable mechanistic pathway, where a ruthenacyclopentene generated by an ene-yne oxidative cyclization undergoes a sequential ß-hydride elimination and intramolecular hydroruthenation to form a ruthenacyclohexene, producing the desirable bicyclo[3.3.0]octenes.

Grinding-Triggered Single Crystal-to-Single Crystal Transformation of a Zinc(II) Complex: Mechanochromic Luminescence and Aggregation-Induced Emission Properties.

We first report single crystal X-ray analysis of ground crystals of mechanochromic luminescence (MCL) that shows single crystal-to-single crystal transformation (SCSCT). Single crystals of [ZnL2] (1-SG, HL = 2-[[[4-(2-benzoxazolyl)phenyl]imino]-methyl]-5-(diethylamino)-phenol) were obtained upon slight grinding of single crystals of [ZnL2]·0.5CH3OH (1), both of which were characterized by single crystal X-ray diffraction. Crystals of 1 showed emission centered at 647 nm (red color), while crystals of 1-SG showed emission band at 624 nm (orange-red color) under UV light, indicating MCL property of the Zn(II) complex. Reversible MCL property with emission color change between red and yellow for 1 was observed upon high grinding and recrystallization. Single crystal X-ray analysis suggested that it is due to the alteration of molecular conformation of ligands in ZnL2 instead of weak intermolecular interaction that 1 exhibits MCL. Investigation of the control Zn(II) complexes (2-4) indicated that flexible substituents and rotated aromatic rings are desirable to generate the MCL-active complexes. In addition, 1 was highly fluorescent in THF solution, but its fluorescence quenched upon addition of water. DFT calculations suggested that this is due to the formation of the excited hydrated ZnL2 species via Zn-O coordination bond, which results in electron-driven proton transfer (EDPT). Aggregates formed as water fraction ( fw) in THF/H2O (v/v) reached 70%, and fluorescence emission was enhanced. This phenomenon continued until fw was 90%, indicating aggregation-induced emission (AIE) property. The mechanism of AIE of ZnL2 in THF/H2O is the restriction of intramolecular rotation (RIR).

Photophysical properties of a D-π-A Schiff base and its applications in the detection of metal ions.

A D-π-A Schiff-base compound, 2-amino-3({[4-(diethylamino)-2-hydroxylphenyl]methylene}amino)-but-2-enedinitrile (H2L), was prepared using diaminomaleonitrile and 4-(N,N-diethylamino)salicylaldehyde. Compared with H2L at 293 K, a low temperature of 100 K makes the parallel aromatics in H2L come closer and fluorescence emission becomes weaker because of π-π interaction-caused quenching. After grinding the crystals of H2L, the colour changed from dark brown-red to yellow under room light and the fluorescence emission enhanced about 9-fold due to the damage of the intermolecular hydrogen bonds, leading to a decrease of non-radiative transition. H2L showed aggregation-induced emission enhancement (AIEE) characteristics in THF/H2O, whose mechanism is attributed to the restriction of intramolecular rotation (RIR). The UV-Vis spectra of H2L with Cu2+ in THF/H2O showed that at first a CuL complex was formed and subsequently a CuL' complex (H2L' = N',N'-bis(4-N',N'-diethylsalicylidene) ethylenediamine) was obtained. The CuL complex turned into the CuL' complex as time prolonged. H2L acted as a dual channel chemosensor for Cu2+ ions in THF/HEPES (v/v: 2 : 8, pH = 8.0) and the CuL complex was stable in this medium. H2L is also a naked-eye probe for Mn2+ ions in CH3CN. The limits of detection are much lower than the allowable level of copper(ii) and manganese(ii) in drinking water set by the World Health Organization (WHO).

The pivotal role of electronics in preferred alkene over alkyne Ni-carboryne insertions and absolute regioselectivities.

The in situ formation mechanisms of active Ni-carboryne species (COM1) and subsequent alkene/alkyne Ni-C bond insertion priorities, as well as relevant cycloaddition regioselectivities and kinetics, were investigated using the IDSCRF-B3LYP density functional theory (DFT) method, and all atoms were equitably treated at the DGDZVP level. The results reveal the o-carborane species to be energetically hedged into a four-step path (barrier heights 5.3, 19.7, 18.4 and 0.3 kcal mol-1, respectively) prior to being transferred into the active Ni-carboryne species (COM1) with the assistance of nBuLi and NiCl2(PPh3)2 at room temperature. In direct agreement with empirical trends, alkene insertion into Ni-C bonds on COM1 is exclusively favoured over the competing alkyne insertion. Electronic structure analyses of the corresponding transition structures showed that the preference of alkenes to alkynes is due to different bonding characteristics during this insertion process, namely, back donation for alkenes but donation for alkyne insertion, as evidenced by molecular graphics and NBO charge distributions. Subsequent alkyne additions (i.e. post alkene insertion) arise as the rate-determining step (RDS) for each of the five different reactions (a-e) explored. The solution free-energy barriers of these RDSs (30.5-38.5 kcal mol-1) were in quantitative agreement with their corresponding experimental yields, evidencing the reliability of the DFT results to reproduce chemical phenomena and energetic trends in real Ni-catalysed carboryne-alkene/alkyne cycloadditions.

DFT studies on the distinct mechanisms of C-H activation and oxidation reactions mediated by mononuclear- and binuclear-palladium.

A series of density functional theory calculations have been carried out to investigate the detailed mechanisms of C-H activation and oxidation reactions, and further to disclose the distinct effects of mononuclear- and binuclear-palladium on these reaction pathways. The results of calculations demonstrated that the C-H activation of 2-phenylpyridine with mononuclear Pd(OAc)2 prefers the inner-shell proton-abstraction mechanism, while that with binuclear Pd2(µ-OAc)4 is biased to the outer-shell proton-abstraction mechanism. The rate-determining free-energy barriers of the two mechanisms were calculated to be 24.2 and 24.8 kcal mol-1, respectively. More importantly, we have simulated the oxidation pathways of PdII → PdIII and PdII → PdIV with strong oxidants including PhI(OAc)2, PhICl2 and NCS, and found that a binuclear PdII-precursor would be oxidized to the corresponding binuclear PdIII-complex while a mononuclear PdII-precursor was deemed to evolve to the corresponding mononuclear PdIV-complex. In addition, the oxidation of PdII with PhI(OAc)2 has been characterized as a radical mechanism, in sharp contrast to the ion-pair mechanism prevalent for the oxidation of PdII with PhICl2. The calculated kinetic and thermodynamic parameters could be qualitatively consistent with the related experimental observations. This molecular modelling can provide valuable insights into the understanding of the distinct effects of the resting state and oxidant on these important transformations.

PdIV Species Mediation in PdII-Catalyzed Direct Alkylation of Arenes with Oxiranes: A DFT Study.

The reaction mechanisms of Pd(OAc)2-catalyzed dehydrogenative alkylation of 2-phenylpyridine with oxirane were investigated using DFT calculations. The most plausible reaction pathway was confirmed as a PdII/IV/II catalytic cycle consisting of four processes: C-H activation, ring-opening oxidative addition of oxirane, reductive elimination, and recovery of the catalyst. According to the B2PLYP/DGDZVP computational data, the oxidative addition of oxirane for converting PdII to PdIV was assigned to be the rate-determining step with a free-energy barrier of 28.1 kcal·mol-1. For comparison, we also studied the alternative PdII-only pathway without a change of oxidation state and found that it was hindered kinetically by a high free-energy barrier of 75.1 kcal·mol-1 occurring for the ring-opening migratory insertion of oxirane. In addition, the small-ring strain of oxirane should be responsible for the feasible C-O bond-cleavage and subsequent PdII → PdIV conversion, because the designed four-, five-, and six-membered-ring reagents did not display such an oxidative addition reactivity. Lastly, an extended reactivity order among oxirane, PhI, PhBr, and PhCl toward oxidative addition onto PdII to form PdIV was proposed in this article based on the computed kinetic parameters.

Energy transfer between rare earths in layered rare-earth hydroxides.

Energy transfer between rare earths in layered rare-earth hydroxides (LRHs) is worth the intensive study because the hydroxyls that act as the bridge connecting the neighbouring rare earths would generate non-radiative transitions. This study focuses on the energy transfer in the intralayer and the adjacent layers of LRHs. A series of LEu x Tb1-x Hs (x = 0, 0.05, 0.2, 0.5, 0.8, and 0.95) was synthesized, the basal spacing (d basal) was adjusted from 8.3 to 46 Å through ion-exchange process, and unilamellar nanosheets were prepared through a delamination process. The luminescence behaviours of the samples demonstrated the following: (1) for the delaminated nanosheets, the quenching effect of both Eu3+ and Tb3+ was hardly observed. This implies that in the intralayer, the efficiency of energy transfer is extremely low, so that highly-concentrated co-doping does not influence the luminescence and by controlling the Eu/Tb molar ratio, white light can be obtained. (2) For small d basal, e.g., 27 Å, the fluorescence quenching of Tb3+ and Eu3+ was remarkable, while for large d basal, e.g., 46 Å, the emission of Tb3+ emerged and the self-quenching between Eu3+ ions weakened. (3) The energy transfer efficiency deceased with an increase in the distance between adjacent layers. In other words, either the energy transfer between Eu3+ and Tb3+ or the energy migration between Eu3+ ions was more efficient when they were located in adjacent layers than in intralayers even when they were the nearest neighbours.

Coordination-Directed Stacking and Aggregation-Induced Emission Enhancement of the Zn(II) Schiff Base Complex.

2-(Trityliminomethyl)-quinolin-8-ol (HL) and its Zn(II) complex were synthesized and characterized by single-crystal X-ray diffraction. HL is an unsymmetrical molecule and coordinated with Zn(II) ion to form ZnL2 in the antiparallel-mode arrangement via Zn-O (hydroxyl group) and Zn-N (quinoline ring) of HL. A high degree of ZnL2 molecules ordering stacking is formed by the coordination bonds and intermolecular π-π interactions, in which head-to-tail arrangement (J-mode stacking) for L- is found. HL is nonfluorescent and ZnL2 is weakly fluorescent in THF. The fluorescence emission of ZnL2 enhances in THF/H2O as H2O% (volume %) is above 60% and aggregates particles with several hundred nanometers are formed, which is confirmed by DLS data and TEM images. The J-aggregates stacking for L- in ZnL2 results in aggregation-induced emission enhancement (AIEE) for ZnL2 in THF/H2O. Theoretical computations based on B3LYP/6-31G(d, p) and TD-B3LYP/6-31G(d, p) methods were carried out. ESIPT is the supposed mechanism for fluorescent silence of HL, and fluorescence emission of ZnL2 is attributed to the restriction of ESIPT process. The oscillator strength of ZnL2 increases from 0.017 for monomer to 0.032 for trimer. It indicates that a high degree of ZnL2 molecules ordering stacking in THF/H2O is of benefit to fluorescence enhancement. HL is an ESIPT-coupled AIEE chemosensor for Zn(II) with high selectivity and sensitivity in aqueous medium. HL can efficiently detect intracellular Zn(II) ions because of ESIPT-coupled AIEE property of ZnL2 in mixed solvent.

A DFT kinetic study on 1,3-dipolar cycloaddition reactions in solution.

Several popular density functional theory (DFT) methods have been employed to characterize a series of 1,3-dipolar cycloaddition reactions, including the exploration of reaction mechanisms and the calculations of kinetic parameters. Both the gas- and solution-phase translational entropy models have been used to calculate the activation entropies, and the results obtained from the latter method are quite close to the experimental measurements. For some of the reactions studied, e.g., a1 + b9, a1 + b10, a5 + b9 and a12 + b5, the explicit + implicit solvation model, namely, a cluster-continuum type model, should be employed to account for the specific solvent-solute interactions. The quasi rigid-rotor-harmonic-oscillator (qRRHO) small frequency vibrational entropy correction, in conjunction with the conformational entropy correction, could further improve the calculated activation entropy data. The comparison between calculation data with experimental measurements, using 23 activation entropies and 160 reaction rate constants as test benchmark, demonstrated that our present strategy could calibrate the root-mean-square-deviation (RMSD) of activation entropies to be 1.8 cal mol-1 K-1 and that of Gibbs free energy barriers to be 1.8 kcal mol-1.

Multi-Pathway Consequent Chemoselectivities of CpRuCl(PPh3 )2 /MeI-Catalysed Norbornadiene Alkyne Cycloadditions.

Chemoselectivities of five experimentally realised CpRuCl(PPh3 )2 /MeI-catalysed couplings of 7-azabenzo-norbornadienes with selected alkynes were successfully resolved from multiple reaction pathway models. Density functional theory calculations showed the following mechanistic succession to be energetically plausible: (1) CpRuI catalyst activation; (2) formation of crucial metallacyclopentene intermediate; (3) cyclobutene product (P2) elimination (ΔGRel(RDS) ≈11.9-17.6 kcal mol-1 ). Alternative formation of dihydrobenzoindole products (P1) by isomerisation to azametalla-cyclohexene followed by subsequent CpRuI release was much less favourable (ΔGRel(RDS) ≈26.5-29.8 kcal mol-1 ). Emergent stereoselectivities were in close agreement with experimental results for reactions a, b, e. Consequent investigations employing dispersion corrections similarly support the empirical findings of P1 dominating in reactions c and d through P2→P1 product transformations as being probable (ΔG≈25.3-30.1 kcal mol-1 ).