Joshua J. Melko

Activation of Methane by FeO+: Determining Reaction Pathways through Temperature-Dependent Kinetics and Statistical Modeling

The temperature dependences of the rate constants and product branching ratios for the reactions of FeO(+) with CH4 and CD4 have been measured from 123 to 700 K. The 300 K rate constants are 9.5 × 10(-11) and 5.1 × 10(-11) cm(3) s(-1) for the CH4 and CD4 reactions, respectively. At low temperatures, the Fe(+) + CH3OH/CD3OD product channel dominates, while at higher temperatures, FeOH(+)/FeOD(+) + CH3/CD3 becomes the majority channel. The data were found to connect well with previous experiments at higher translational energies. The kinetics were simulated using a statistical adiabatic channel model (vibrations are adiabatic during approach of the reactants), which reproduced the experimental data of both reactions well over the extended temperature and energy ranges. Stationary point energies along the reaction pathway determined by ab initio calculations seemed to be only approximate and were allowed to vary in the statistical model. The model shows a crossing from the ground-state sextet surface to the excited quartet surface with large efficiency, indicating that both states are involved. The reaction bottleneck for the reaction is found to be the quartet barrier, for CH4 modeled as -22 kJ mol(-1) relative to the sextet reactants. Contrary to previous rationalizations, neither less favorable spin-crossing at increased energies nor the opening of additional reaction channels is needed to explain the temperature dependence of the product branching fractions. It is found that a proper treatment of state-specific rotations is crucial. The modeled energy for the FeOH(+) + CH3 channel (-1 kJ mol(-1)) agrees with the experimental thermochemical value, while the modeled energy of the Fe(+) + CH3OH channel (-10 kJ mol(-1)) corresponds to the quartet iron product, provided that spin-switching near the products is inefficient. Alternative possibilities for spin switching during the reaction are considered. The modeling provides unique insight into the reaction mechanisms as well as energetic benchmarks for the reaction surface.

Effect of Charge and Composition on the Structural Fluxionality and Stability of Nine Atom Tin−Bismuth Zintl Analogues

Synergistic studies of bismuth doped tin clusters combining photoelectron spectra with first principles theoretical investigations establish that highly charged Zintl ions, observed in the condensed phase, can be stabilized as isolated gas phase clusters through atomic substitution that preserves the overall electron count but reduces the net charge and thereby avoids instability because of coulomb repulsion. Mass spectrometry studies reveal that Sn(8)Bi(-), Sn(7)Bi(2)(-), and Sn(6)Bi(3)(-) exhibit higher abundances than neighboring species, and photoelectron spectroscopy show that all of these heteroatomic gas phase Zintl analogues (GPZAs) have high adiabatic electron detachment energies. Sn(6)Bi(3)(-) is found to be a particularly stable cluster, having a large highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap. Theoretical calculations demonstrate that the Sn(6)Bi(3)(-) cluster is isoelectronic with the well know Sn(9)(-4) Zintl ion; however, the fluxionality reported for Sn(9)(-4) is suppressed by substituting Sn atoms with Bi atoms. Thus, while the electronic stability of the clusters is dominated by electron count, the size and position of the atoms affects the dynamics of the cluster as well. Substitution with Bi enlarges the cage compared with Sn(9)(-4) making it favorable for endohedral doping, findings which suggest that these cages may find use for building blocks of cluster assembled materials.

Temperature dependence of the OH- + CH3I reaction kinetics. experimental and simulation studies and atomic-level dynamics

Direct dynamics simulations and selected ion flow tube (SIFT) experiments were performed to study the kinetics and dynamics of the OH(-) + CH3I reaction versus temperature. This work complements previous direct dynamics simulation and molecular beam ion imaging experiments of this reaction versus reaction collision energy (Xie et al. J. Phys. Chem. A 2013, 117, 7162). The simulations and experiments are in quite good agreement. Both identify the SN2, OH(-) + CH3I → CH3OH + I(-), and proton transfer, OH(-) + CH3I → CH2I(-) + H2O, reactions as having nearly equal importance. In the experiments, the SN2 pathway constitutes 0.64 ± 0.05, 0.56 ± 0.05, 0.51 ± 0.05, and 0.46 ± 0.05 of the total reaction at 210, 300, 400, and 500 K, respectively. For the simulations this fraction is 0.56 ± 0.06, 0.55 ± 0.04, and 0.50 ± 0.05 at 300, 400, and 500 K, respectively. The experimental total reaction rate constant is (2.3 ± 0.6) × 10(-9), (1.7 ± 0.4) × 10(-9), (1.9 ± 0.5) × 10(-9), and (1.8 ± 0.5) × 10(-9) cm(3) s(-1) at 210, 300, 400, and 500 K, respectively, which is approximately 25% smaller than the collision capture value. The simulation values for this rate constant are (1.7 ± 0.2) × 10(-9), (1.8 ± 0.1) × 10(-9), and (1.6 ± 0.1) × 10(-9) cm(3)s(-1) at 300, 400, and 500 K. From the simulations, direct rebound and stripping mechanisms as well as multiple indirect mechanisms are identified as the atomic-level reaction mechanisms for both the SN2 and proton-transfer pathways. For the SN2 reaction the direct and indirect mechanisms have nearly equal probabilities; the direct mechanisms are slightly more probable, and direct rebound is more important than direct stripping. For the proton-transfer pathway the indirect mechanisms are more important than the direct mechanisms, and stripping is significantly more important than rebound for the latter. Calculations were performed with the OH(-) quantum number J equal to 0, 3, and 6 to investigate the effect of OH(-) rotational excitation on the OH(-) + CH3I reaction dynamics. The overall reaction probability and the probabilities for the SN2 and proton-transfer pathways have little dependence on J. Possible effects on the atomistic mechanisms were investigated for the SN2 pathway and the probability of the direct rebound mechanism increased with J. However, the other atomistic mechanisms were not appreciably affected by J.

Further Insight into the Reaction FeO+ + H2 → Fe+ + H2O: Temperature Dependent Kinetics, Isotope Effects, and Statistical Modeling

The reactions of FeO(+) with H2, D2, and HD were studied in detail from 170 to 670 K by employing a variable temperature selected ion flow tube apparatus. High level electronic structure calculations were performed and compared to previous theoretical treatments. Statistical modeling of the temperature and isotope dependent rate constants was found to reproduce all data, suggesting the reaction could be well explained by efficient crossing from the sextet to quartet surface, with a rigid near thermoneutral barrier accounting for both the inefficiency and strong negative temperature dependence of the reactions over the measured range of thermal energies. The modeling equally well reproduced earlier guided ion beam results up to translational temperatures of about 4000 K.

Exploring the reactions of Fe+ and FeO+ with NO and NO2.

We report for the first time temperature dependences (from 300 to 600 K) of the reactions of Fe(+) and FeO(+) with NO and NO(2). Both ions react quickly with NO(2), and their rate constants have weak negative temperature dependences. The former is consistent with the calculated energy profile along the Fe(+) + NO(2) reaction coordinate. Ground state Fe(+) reacts with NO(2) to produce only FeO(+), while FeO(+) reacts with NO(2) to produce NO(+) exclusively. Certain source conditions produce excited Fe(+), as evidenced by production of primary NO(+), which is endothermic with the ground state by 0.35 eV. The room temperature rate constants are in agreement with previous values. For the reactions of Fe(+) and FeO(+) with NO, we find an upper limit of <1.0 × 10(-12) cm(3) s(-1) for both rate constants, in contrast to a previous report of a rate constant of ∼1.7 × 10(-11) cm(3) s(-1) for Fe(+) + NO. Because this is an endothermic process, the prior report interpreted the reaction as a termolecular process involving two NO molecules; instead, we show that the previous results were likely due to an NO(2) impurity. Implications for other metal cation reactions which have been speculated to occur by the termolecular mechanism are discussed.

Iron cation catalyzed reduction of N2O by CO: Gas-phase temperature dependent kinetics.

The ion-molecule reactions Fe(+) + N2O → FeO(+) + N2 and FeO(+) + CO → Fe(+) + CO2, which catalyze the reaction CO + N2O → CO2 + N2, have been studied over the temperature range 120-700 K using a variable temperature selected ion flow tube apparatus. Values of the rate constants for the former two reactions were experimentally derived as k2 (10(-11) cm(3) s(-1)) = 2.0(±0.3) (T/300)(-1.5(±0.2)) + 6.3(±0.9) exp(-515(±77)/T) and k3 (10(-10) cm(3) s(-1)) = 3.1(±0.1) (T/300)(-0.9(±0.1)). Characterizing the energy parameters of the reactions by density functional theory at the B3LYP/TZVP level, the rate constants are modeled, accounting for the intermediate formation of complexes. The reactions are characterized by nonstatistical intrinsic dynamics and rotation-dependent competition between forward and backward fluxes. For Fe(+) + N2O, sextet-quartet switching of the potential energy surfaces is quantified. The rate constant for the clustering reaction FeO(+) + N2O + He → FeO(N2O)(+) + He was also measured, being k4 (10(-27) cm(6) s(-1)) = 1.1(±0.1) (T/300)(-2.5(±0.1)) in the low pressure limit, and analyzed in terms of unimolecular rate theory.

Combined Experimental and Theoretical Study of AlnX (n = 1−6; X = As, Sb) Clusters: Evidence of Aromaticity and the Jellium Model

The electronic structure of Al(n)X (n = 1-6; X = As, Sb) clusters has been investigated using a synergistic approach combining negative ion photoelectron spectroscopy and first principles electronic structure calculations. It is shown that Al(3)X and Al(5)X exhibit enhanced energetic stability, as evidenced from calculated removal and embedding energies as well as chemical stability manifested through a large gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). However, the stabilities of these species are derived from different mechanisms. Al(3)As and Al(3)Sb, with HOMO-LUMO gaps of 1.86 and 1.73 eV, respectively, are shown to have planar geometries where the p orbitals combine to form one pi and two sigma aromatic orbitals reminiscent of conventional all-metal aromatic species. Al(5)As and Al(5)Sb, with 20 valence electrons, possess a closed electronic shell (1s(2), 1p(6), 1d(10), 2s(2)) within a jellium framework and have HOMO-LUMO gaps of 1.12 and 1.17 eV, respectively.

A novel technique for measurement of thermal rate constants and temperature dependences of dissociative recombination: CO2+, CF3+, N2O+, C7H8+, C7H7+, C6H6+, C6H5+, C5H6+, C4H4+, and C3H3+

A novel technique using a flowing afterglow-Langmuir probe apparatus for measurement of temperature dependences of rate constants for dissociative recombination (DR) is presented. Low (~10(11) cm(-3)) concentrations of a neutral precursor are added to a noble gas∕electron afterglow plasma thermalized at 300-500 K. Charge exchange yields one or many cation species, each of which may undergo DR. Relative ion concentrations are monitored at a fixed reaction time while the initial plasma density is varied between 10(9) and 10(10) cm(-3). Modeling of the decrease in concentration of each cation relative to the non-recombining noble gas cation yields the rate constant for DR. The technique is applied to several species (O2(+), CO2(+), CF3(+), N2O(+)) with previously determined 300 K values, showing excellent agreement. The measurements of those species are extended to 500 K, with good agreement to literature values where they exist. Measurements are also made for a range of CnHm(+) (C7H7(+), C7H8(+), C5H6(+), C4H4(+), C6H5(+), C3H3(+), and C6H6(+)) derived from benzene and toluene neutral precursors. CnHm(+) DR rate constants vary from 8-12 × 10(-7) cm(3) s(-1) at 300 K with temperature dependences of approximately T(-0.7). Where prior measurements exist these results are in agreement, with the exception of C3H3(+) where the present results disagree with a previously reported flat temperature dependence.

The applicability of three-dimensional aromaticity in BiSnn- Zintl analogues

Three-dimensional aromaticity is shown to play a role in the stability of deltahedral Zintl clusters and here we examine the connection between aromaticity and stability. In order to gain further insight, we have studied Zintl analogs comprised of bismuth doped tin clusters with photoelectron spectroscopy and theoretical methods. To assign aromaticity, we examine the ring currents induced around the cage by using the nucleus independent chemical shift. In the current study, BiSn(4)(-) is a stable cluster and fits aromatic criteria, while BiSn(5)(-) is found to fit antiaromatic criteria and has reduced stability. The more stable clusters exhibit an aromatic character which originates from weakly interacting s-states and bonding orbitals parallel to the surface of the cluster, while nonbonding lone pairs perpendicular to the surface of the cluster account for antiaromaticity and reduced stability. The effect of three-dimensional aromaticity on the electronic structure does not result in degeneracies, so the resulting variations in stability are smaller than those seen in conventional aromaticity.

Temperature Dependences for the Reactions of O2- and O- with N and O Atoms in a Selected-Ion Flow Tube Instrument

Rate constants for the reactions of O2(-) and O(-) with N and O atoms have been measured for the first time as a function of temperature from 173 to 500 K for O(-) reactions and 173 to 400 K for O2(-) reactions. Room temperature rate constants for O2(-) reacting with N and O are 3.1 × 10(-10) and 1.7 × 10(-10) cm(3) s(-1), respectively, and the corresponding O(-) rate constants are 1.7 × 10(-10) and 1.5 × 10(-10) cm(3) s(-1), in good agreement with previous values. Temperature dependences are about T(-1.7) for both O2(-) reactions and T(-0.6) and T(-1.3) for the reactions of O(-) with N and O, respectively. Branching for the O2(-) reaction with N is found to predominantly form O(-) (>85%) in contrast to previous measurements, which reported NO2 + e(-) as the main channel. Calculations point to the present results being correct. The potential energy surface for this reaction was calculated using density functional theory, coupled cluster with singles, doubles (triples), complete active space self-consistent field, and complete active space second-order perturbation methods and is found to be quite complex, with agreement between the calculated surface and the observed kinetic data only possible through the inclusion of dynamical correlation.