Activation of CH4 by gas-phase Ni+ and the thermochemistry of Ni-ligand complexes

Infrared Spectroscopic Characterization of Early 4d Transition Metal Carbene Cations, ZrCH2+ and NbCH2. J Phys Chem A 2024;128:6658-6667. [PMID: 39083656 DOI: 10.1021/acs.jpca.4c03553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

IR multiple-photon dissociation (IRMPD) action spectroscopy is combined with quantum chemical calculations to examine the [M,C,2H]+ species for the early 4d metals, M = Zr and Nb. These ions were formed by reacting laser ablated M+ ions with cyclopropane (c-C3H6) in a molecular beam apparatus. Both IRMPD spectra exhibit one major band near 700 cm-1 and a second weaker band at about twice that wavenumber, more evident when irradiated in focus. The [Nb,C,2H]+ species also has a sharp band at 800 cm-1. Comparison with B3LYP calculations allow assignment of the [M,C,2H]+ structures to agostic carbenes, which is similar to the structures found for the 5d analogues, WCH2+ and TaCH2+. A molecular orbital analysis traces the reasons for the agostic deformation from a classic C2v symmetric carbene.

IR spectroscopic characterization of [M,C,2H]+ (M = Ru and Rh) products formed by reacting 4d transition metal cations with oxirane: Spectroscopic evidence for multireference character in RhCH2. Phys Chem Chem Phys 2024;26:11445-11458. [PMID: 38572552 PMCID: PMC11022548 DOI: 10.1039/d4cp00012a]

A combination of infrared multiple-photon dissociation (IRMPD) action spectroscopy and quantum chemical calculations was employed to investigate the [M,C,2H]+ (M = Ru and Rh) species. These ions were formed by reacting laser ablated M+ ions with oxirane (ethylene oxide, c-C2H4O) in a room-temperature ion trap. IRMPD spectra for the Ru species exhibit one major band and two side bands, whereas spectra for the Rh species contain more distinct bands. Comparison with density functional theory (DFT), coupled-cluster (CCSD), and equation-of-motion spin-flip CCSD (EOM-SF-CCSD) calculations allows assignment of the [M,C,2H]+ structures. For the spectrum of [Ru,C,2H]+, a combination of HRuCH+ and RuCH2+ structures reproduces the observed spectrum at all levels of theory. The well-resolved spectrum of [Rh,C,2H]+ could not be assigned unambiguously to any calculated structure using DFT approaches. The EOM-SF-CCSD calculations showed that the ground-state surface has multireference electronic character, and symmetric carbenes in both the 1A1 and 3A2 states are needed to reproduce the observed spectrum.

IR spectroscopic characterization of 3d transition metal carbene cations, FeCH2+ and CoCH2+: periodic trends and a challenge for DFT approaches

A combination of IR multiple-photon dissociation (IRMPD) action spectroscopy and quantum chemical calculations was employed to investigate the [M,C,2H]+ (M = Fe and Co) species. These were formed by reacting laser ablated M+ ions with oxirane (ethylene oxide, c-C2H4O) in a room temperature ion trap. IRMPD spectra for the Fe and Co species are very similar and exhibit one major band. Comparison with density functional theory (DFT) and coupled cluster with single and double excitations (CCSD) calculations allows assignment of the spectra to MCH2+ carbene structures. For these 3d transition metal systems, experimental IRMPD spectra compare relatively poorly with DFT calculated IR spectra, but CCSD calculated spectra are a much better match primarily because the M-C stretch gains significant intensity. The origins of this behavior are explored in some detail. The present results are also compared to previous results for the 4d and 5d congeners and the periodic trends in these structures are evaluated.

On the origin of reactivity variation upon sequential ligation: the [Re(Cl)x]+/CH4 (x = 1-3) couples

The potential of [ReCl_x]⁺ (x = 1-3) in activating methane has been explored by using a combination of gas-phase experiments and high-level quantum calculations. When the number of Cl ligands increases, the reactivity towards methane activation varies accordingly. While [ReCl_x]⁺ (x = 1-2) are able to dehydrogenate methane by a three-state reactivity scenario, [ReCl₃]⁺ shows inertness towards methane at ambient conditions. Furthermore, the product ion [ClRe(H)CH]⁺ of the [ReCl]⁺/CH₄ couple could continue to activate methane and liberate molecular dihydrogen but another product ion [Cl₂ReCH₂]⁺ is unreactive with methane. Obviously, the nature and the number of ligands make a difference to the reactivity towards methane activation. The associated reaction mechanism and the electron origins for the rather different reactivities are discussed in detail. Finally and more importantly, instructive information concerning the rational design of Re-catalysts for methane conversion is obtained.

High-Level ab Initio Predictions for the Ionization Energies, Bond Dissociation Energies, and Heats of Formation of Vanadium Methylidene, Vanadium Methyl Species, and Their Cations (VCH2/VCH2+, VCH3/VCH3+)

The ionization energies of VCH₂ and VCH₃, the various 0 K bond dissociation energies (D₀s) in their neutrals and cations, and their respective heats of formation at 0 and 298 K are computed by the single-reference, wave function-based CCSDTQ/CBS procedure. The core of the composite method is the approximation to the complete basis set (CBS) limit at the coupled cluster (CC) level which includes up to full quadruple excitations. The zero-point vibrational energy, core-valence correlation, spin-orbit coupling, and scalar relativistic effects have their contributions incorporated in an additive manner. For the species in the current study, this protocol requires geometry optimizations and harmonic frequency calculations practically no higher than the CCSD(T)/aug-cc-pwCVTZ and CCSD(T)/aug-cc-pVTZ levels, respectively. The present calculations successfully predict D₀(V⁺-CH₃) = 2.126 eV and D₀(V⁺-CH₂) = 3.298 eV in remarkable agreement with the data recently measured by a spin-orbit state selected V⁺ + CH₄ collision experiment (Phys. Chem. Chem. Phys. 2021, 23, 273-286). The good accord encourages the use of CCSDTQ/CBS protocol in thermochemical predictions of various feasible product channels identified in methane activation by transition metal species.

Infrared Spectroscopy of Gold Carbene Cation (AuCH2+): Covalent or Dative Bonding?

The present work explores the structure of the gold carbene cation, AuCH₂⁺, using infrared multiple photon dissociation action spectroscopy and density functional theory (DFT). Unlike several other 5d transition-metal cations (M⁺ = Ta⁺, W⁺, Os⁺, Ir⁺, and Pt⁺) that react with methane by dehydrogenation to form MCH₂⁺ species, gold cations are unreactive with methane at thermal energies. Instead, the metal carbene is formed by reacting atomic gold cations formed in a laser ablation source with ethylene oxide (cC₂H₄O) pulsed into a reaction channel downstream. The resulting [Au,C,2H]⁺ product photofragmented by loss of H₂ as induced by radiation provided by the free-electron laser for intracavity experiments in the 300-1800 cm^-1 range. Comparison of the experimental spectrum, obtained by monitoring the appearance of AuC⁺, and DFT calculated spectra leads to the identification of the ground-state carbene, AuCH₂⁺ (¹A₁), as the species formed, as previously postulated theoretically. Unlike the covalent double bonds formed by the lighter, open-shell 5d transition metals, the closed-shell Au⁺ (¹S, 5d¹⁰) atom binds to methylene by donation of a pair of electrons from CH₂(¹A₁) into the empty 6s orbital of gold coupled with π back-bonding, i.e., dative bonding, as explored computationally. Contributions to the AuC⁺ appearance spectrum from larger complexes are also considered, and H₂CAu⁺(c-C₂H₄O) seems likely to contribute one band observed.

Spectroscopic Identification of the Carbyne Hydride Structure of the Dehydrogenation Product of Methane Activation by Osmium Cations

The present work explores the structures of species formed by dehydrogenation of methane (CH₄) and perdeuterated methane (CD₄) by the 5d transition metal cation osmium (Os⁺). Using infrared multiple photon dissociation (IRMPD) action spectroscopy and density functional theory (DFT), the structures of the [Os,C,2H]⁺ and [Os,C,2D]⁺ products are explored. This study complements previous work on the related species formed by dehydrogenation of methane by four other 5d transition metal cations (M⁺ = Ta⁺, W⁺, Ir⁺, and Pt⁺). Osmium cations are formed in a laser ablation source, react with methane pulsed into a reaction channel downstream, and the resulting products spectroscopically characterized through photofragmentation using the Free-Electron Laser for IntraCavity Experiments (FELICE) in the 300-1800 cm^-1 range. Photofragmentation was monitored by the loss of H₂/D₂. Comparison of the experimental spectra and DFT calculated spectra leads to identification of the ground state carbyne hydride, HOsCH⁺ (²A') as the species formed, as previously postulated theoretically. Further, a full description of the systematic spectroscopic shifts observed for deuterium labeling of these complexes, some of the smallest systems to be studied using IRMPD action spectroscopy, is achieved. A full rotational contour analysis explains the observed linewidths as well as the observation of doublet structures in several bands, consistent with previous observations for HIrCH⁺ (²A'). Graphical Abstract ᅟ.

Infrared Spectroscopy of Au+(CH4) n Complexes and Vibrationally-Enhanced C-H Activation Reactions

A combined spectroscopic and computational study of gas-phase Au⁺(CH₄)_n (n = 3–8) complexes reveals a strongly-bound linear Au⁺(CH₄)₂ core structure to which up to four additional ligands bind in a secondary coordination shell. Infrared resonance-enhanced photodissociation spectroscopy in the region of the CH₄a₁ and t₂ fundamental transitions reveals essentially free internal rotation of the core ligands about the H₄C–Au⁺–CH₄ axis, with sharp spectral features assigned by comparison with spectral simulations based on density functional theory. In separate experiments, vibrationally-enhanced dehydrogenation is observed when the t₂ vibrational normal mode in methane is excited prior to complexation. Clear infrared-induced enhancement is observed in the mass spectrum for peaks corresponding 4u below the mass of the Au⁺(CH₄)_n=2,3 complexes corresponding, presumably, to the loss of two H₂ molecules.

Submerged Barriers in the Ni(+) Assisted Decomposition of Propionaldehyde

The reaction dynamics of the Ni(+) mediated decarbonylation of propionaldehyde was assessed using the single photon initiated decomposition rearrangement reaction (SPIDRR) technique. The exothermic production of Ni(+)CO was temporally monitored and the associated rate constants, k(E), were extracted as a function of activating photon energy. In addition, the reaction potential energy surface was calculated at the UCCSD(T)/def2-TZVP//PBEPBE/cc-pVDZ level of theory to provide an atomistic description of the reaction profile. The decarbonylation of propionaldehyde can be understood as proceeding through parallel competitive reaction pathways that are initiated by Ni(+) insertion into either the C-C or C-H bond of the propionaldehyde carbonyl carbon. Both paths lead to the elimination of neutral ethane and are governed by submerged barriers. The lower energy sequence is a consecutive C-C/C-H addition process with a submerged barrier of 14 350 ± 600 cm(-1). The higher energy sequence is a consecutive C-H/C-C addition process with a submerged barrier of 15 400 ± 600 cm(-1). Both barriers were determined using RRKM calculations fit to the experimentally determined k(E) values. The measured energy difference between the two barriers agrees with the DFT computed difference in rate limiting transition-state energies, 18 413 and 19 495 cm(-1).

Gas-phase perspective on the thermodynamics and kinetics of heterogeneous catalysis

Gas-phase studies of small transition metal cluster cations provide thermochemistry of utility to surface science and heterogeneous catalysis.

Co(+)-assisted decomposition of h6-acetone and d6-acetone: acquisition of reaction rate constants and dynamics of the dissociative mechanism

Reaction rate constants have been acquired for the gaseous unimolecular decomposition reaction of the Co(+)(OC(CH(3))(2)) cluster ion and its deuterium labeled analog. Each rate constant is measured at a well resolved cluster internal energy within the range 12,300-16,100 cm(-1). The weighted, averaged kinetic isotope effect (KIE), k(H)/k(D) = 1.54 ± 0.05, is about three times smaller than the KIE measured for the rate-determining rate constants in the similar Ni(+)(OC(CH(3))(2)) decomposition reaction. These reactions likely follow the same oxidative addition-reductive elimination mechanism. Thus, this unexpected change in the KIE magnitudes is not due to differences in the dissociative reaction coordinates. Rather, we propose that the unique dissociation dynamics of these two similar systems is due to differences in the low-lying electronic structure of each transition metal ion.

Reaction rate constants and mechanistic detail of the Ni+ + butanone reaction

The unimolecular decomposition kinetics of the jet-cooled Ni(+)-butanone cluster ion has been monitored over a range of internal energies (16000-18800 cm⁻¹). First-order rate constants are acquired for the precursor ion dissociation into three product channels. The temporal growth of each fragment ion is selectively monitored in a custom instrument and yields similar valued rate constants at a common ion internal energy. The decomposition reaction is proposed to proceed along two parallel reaction coordinates. Each dissociative pathway is rate-limited by the initial Ni(+) oxidative addition into either the C-CH₃ or C-C₂H₅ σ-bond in the butanone molecule. Ratios of integrated product ion intensities as well as the measured rate constants are used to determine values for each σ-bond activation rate constant. The lowest energy measurement presented in this study occurs when the binary complex ion possesses an internal energy of 16000 cm⁻¹. Under this condition, the Ni(+) assisted decomposition of the butanone molecule is rate limited by k(act)(C-C₂H₅) = (0.92 ± 0.08) × 10⁵ s⁻¹ and k(act)(C-CH₃) = (0.37 ± 0.03) × 10⁵ s⁻¹. The relative magnitudes of the two rate constants reflect the greater probability for reaction to occur along the C-C₂H₅ σ-bond insertion pathway, consistent with thermodynamic arguments. DFT calculations at the B3LYP/6-311++G(d,p) level of theory suggest the most likely geometries and relative energies of the reactants, intermediates, and products.

Low-energy reaction rate constants for the Ni+-assisted decomposition of acetaldehyde: observation of C-H and C-C activation

Rate constants for the low-energy Ni(+)-assisted dissociative reaction of acetaldehyde have been measured under jet-cooled conditions in the gas phase. The rate constants are acquired through monitoring the time dependence of fragment Ni(+)CO formation. The decomposition of the precursor Ni(+)-acetaldehyde cluster ion proceeds via consecutive, parallel reaction coordinates that originate with the Ni(+)-assisted cleavage of either a C-C or an aldehyde C-H bond. The energies used to initiate these reactions are well below that required to cleave sigma-bonds in the isolated acetaldehyde molecule. Direct measurement of the reaction kinetics over a range of energies indicates that the rate-limiting step in the dissociative mechanism changes at cluster ion internal energies = 17,200 +/- 400 cm(-1). Arguments are presented that this energy marks the closure of the dissociative coordinate that initiates with C-H sigma-bond activation and thus provides a measure of the activation energy of this dissociative pathway.

Rate-limiting step in the low-energy unimolecular decomposition reaction of Ni+* acetone into Ni+CO + ethane

Rate constants for the low-energy Ni(+)-assisted C-C bond cleavage reaction of deuterium-labeled acetone have been acquired under jet-cooled conditions in the gas phase. The energies used to initiate the dissociative reactions of the precursor complex ion Ni(+)(d(6)-Ac) are well below that required to cleave C-C sigma-bonds in isolated organic molecules. The rate constants are compared to those acquired previously for the lighter Ni(+)(h(6)-Ac) isotope and result in a substantial kinetic isotope effect (k(H)/k(D) approximately 5.5). Arguments are made that implicate isomerization leading to C-C bond coupling as the rate-limiting step (not C-C sigma-bond activation) in the dissociative reaction.

The low-energy unimolecular reaction rate constants for the gas phase, Ni+-mediated dissociation of the C-C sigma bond in acetone

The time dependence of the gaseous unimolecular decomposition of the jet-cooled adduct ion, Ni+-OC(CH3)2, was monitored through selective detection of the Ni+CO fragment ion. Various resolved amounts of energy in the range 15600-18800 cm(-1) were supplied to initiate the dissociation reaction through absorption of laser photons by the title molecular complex. First-order rate constants, k(E), ranged from 113000 to 55000 s(-1) and decreased with decreasing amounts of internal excitation. The energy used to initiate the reaction is well below that required to fragment C-C sigma bonds and indicates the necessity of the Ni+ cation to induce bond activation and fragmentation. These measurements are carried out in a unique apparatus and represent the first direct kinetic study of such catalytic type reactions.

Reactions of atomic cations with methane: gas phase room-temperature kinetics and periodicities in reactivity

Reactions of methane have been measured with 59 atomic metal cations at room temperature in helium bath gas at 0.35 Torr using an inductively-coupled plasma/selected-ion flow tube (ICP/SIFT) tandem mass spectrometer. The atomic cations were produced at approximately 5500 K in an ICP source and allowed to decay radiatively and to thermalize by collisions with argon and helium atoms prior to reaction. Rate coefficients and product distributions are reported for the reactions of fourth-row atomic cations from K(+) to Se(+), of fifth-row atomic cations from Rb(+) to Te(+) (excluding Tc(+)), of sixth-row atomic cations from Cs(+) to Bi(+), and of the lanthanide cations from La(+) to Lu(+) (excluding Pm(+)). Two primary reaction channels were observed: C-H bond insertion with elimination of H(2), and CH(4) addition. The bimolecular H(2) elimination was observed in the reactions of CH(4) with As(+), Nb(+), and some sixth-row metal cations, i.e., Ta(+), W(+), Os(+), Ir(+), Pt(+); secondary and higher-order H(2) elimination was observed exclusively for Ta(+), W(+), and Ir(+) ions. All other transition-metal cations except Mn(+) and Re(+) were observed to react with CH(4) exclusively by addition, and up to two methane molecules were observed to add sequentially to most transition-metal ions. CH(4) addition was also observed for Ge(+), Se(+), La(+), Ce(+), and Gd(+) ions, while the other main-group and lanthanide cations did not react measurably with methane.

Gas-phase activation of methane by ligated transition-metal cations

Motivated by the search for ways of a more efficient usage of the large, unexploited resources of methane, recent progress in the gas-phase activation of methane by ligated transition-metal ions is discussed. Mass spectrometric experiments demonstrate that the ligands can crucially influence both reactivity and selectivity of transition-metal cations in bond-activation processes, and the most reactive species derive from combinations of transition metals with the electronegative elements fluorine, oxygen, and chlorine. Furthermore, the collected knowledge about intramolecular kinetic isotope effects associated with the activation of C-H(D) bonds of methane can be used to distinguish the nature of the bond activation as a mere hydrogen-abstraction, a metal-assisted mechanism or more complex reactions such as formation of insertion intermediates or sigma-bond metathesis.

Ligand and Substrate Effects in Gas-Phase Reactions of NiX+/RH Couples (X=F, Cl, Br, I; R=CH3, C2H5,nC3H7,nC4H9)

The reactions of small saturated hydrocarbons by gaseous nickel cations NiX+ (X=F, Cl, Br, I) are investigated by means of electrospray ionization mass spectrometry. The halide cations are obtained from solutions of the corresponding Ni(II) salts in water or methanol as solvents. NiF+ is the only Ni(II) halide complex that brings about thermal activation of methane. The branching ratios of the observed reactions with C2H6, C3H8, and nC4H10 are shifted systematically by changing the nature of both the ligand X and the substrate RH. In the elimination of HX (X=F, Cl, Br, I), the formal oxidation state of the metal ion appears to be conserved, and the importance of this reaction channel decreases in going from NiF+ to NiI+. A reversed trend is observed in the losses of small closed-shell neutral molecules, that is, H2, CH4 and C2H6, which dominate the gas-phase ion chemistry of NiI+/RH couples. Additionally, inner-sphere electron-transfer reactions take place for a few systems, that is, the delivery of hydride or methanide ions from the hydrocarbon to NiX+ in the course of which the hydrocarbon is converted to a carbenium ion and the cationic metal complex gives rise to a neutral RNiX molecule (R=H, CH3). This process gains importance with decreasing atomic number of the halides and with increasing the size of the alkane. Thus, it constitutes the major pathway in the reactions of NiF+ with propane and n-butane, whereas it is not observed for any of the NiI+/RH couples investigated. Concerning the regioselectivity of the reactions with propane and n-butane, heterolytic cleavage of secondary carbon--hydrogen bonds is clearly preferred compared to that of primary ones, as revealed by deuterium labeling studies. For the NiF+/C3H8 couple, the selectivity of the hydride transfer is as large as 360 in favor of the secondary positions. Though smaller, large preferences for the activation of secondary C--H bonds are also operative in homolytic bond activation of RH (R=nC3H7, nC4H9).