Reactive Force Field Development for Propane Dehydrogenation on Platinum Surfaces.

Propane dehydrogenation (PDH) is an on-purpose catalytic technology to produce propylene from propane that operates at high temperatures, 773-973 K. Several key industry players have been active in developing new catalysts and processes with improved carbon footprint and economics, where Pt-based catalysts have played a central role. The optimization of these catalytic systems through computational and atomistic simulations requires large-scale models that account for their reactivity and dynamic properties. To address this challenge, we developed a new reactive ReaxFF force field (2023-Pt/C/H) that enables large-scale simulations of PDH reactions catalyzed on Pt surfaces. The optimization of force-field parameters relies on a large training set of density functional theory (DFT) calculations of Pt-catalyzed PDH mechanism, including geometries, adsorption and relative energies of reaction intermediates, and key C-H and C-C bond-breaking/forming reaction steps on the Pt(111) surface. The internal validation supports the accuracy of the developed 2023-Pt/C/H force-field parameters, resulting in mean absolute errors (MAE) against DFT data of 14 and 12 kJ mol-1 for relative energies of intermediates and energy barriers, respectively. We demonstrated the applicability of the 2023-Pt/C/H force field with reactive molecular dynamics simulations of propane on different Pt surface topologies and temperatures. The simulations successfully model the formation of propene in the gas phase as well as competitive, unproductive reactions such as deep dehydrogenation and C-C bond cleavage that produce H, C1 and C2 adsorbed species responsible of catalytic deactivation of Pt surface. Results show the following reactivity order: Pt(111) < Pt(100) < Pt(211), and that for the stepped Pt(211) surface, propane activation occurs on low-coordinated Pt atoms at the steps. The measured selectivity as a function of surface topology follows the same trend as activity, the Pt(211) facet being the most selective. The 2023-Pt/C/H reactive force field can also describe the increase of reactivity with the temperature. From these simulations, we were able to estimate the Arrhenius activation energy, 73 kJ mol-1, whose value is close to those reported experimentally for PDH catalyzed by large, supported Pt nanoparticles . The newly developed 2023-Pt/C/H reactive force field can be used in subsequent investigations of different Pt topologies and of collective effects such as temperature, propane pressure, or H surface coverage.

Bringing Selectivity in H/D Exchange Reactions Catalyzed by Metal Nanoparticles through Modulation of the Metal and the Ligand Shell.

Ru and Rh nanoparticles catalyze the selective H/D exchange in phosphines using D2 as the deuterium source. The position of the deuterium incorporation is determined by the structure of the P-based substrates, while activity depends on the nature of the metal, the properties of the stabilizing agents, and the type of the substituent on phosphorus. The appropriate catalyst can thus be selected either for the exclusive H/D exchange in aromatic rings or also for alkyl substituents. The selectivity observed in each case provides relevant information on the coordination mode of the ligand. Density functional theory calculations provide insights into the H/D exchange mechanism and reveal a strong influence of the phosphine structure on the selectivity. The isotope exchange proceeds via C-H bond activation at nanoparticle edges. Phosphines with strong coordination through the phosphorus atom such as PPh3 or PPh2Me show preferred deuteration at ortho positions of aromatic rings and at the methyl substituents. This selectivity is observed because the corresponding C-H moieties can interact with the nanoparticle surface while the phosphine is P-coordinated, and the C-H activation results in stable metallacyclic intermediates. For weakly coordinating phosphines such as P(o-tolyl)3, the interaction with the nanoparticle can occur directly through phosphine substituents, and then, other deuteration patterns are observed.

N═N Bond Cleavage by Tantalum Hydride Complexes: Mechanistic Insights and Reactivity.

The reaction of [TaCpRX4] (CpR = η5-C5Me5, η5-C5H4SiMe3, η5-C5HMe4; X = Cl, Br) with SiH3Ph resulted in the formation of the dinuclear hydride tantalum(IV) compounds [(TaCpRX2)2(µ-H)2], structurally identified by single-crystal X-ray analyses. These species react with azobenzene to give the mononuclear imide complex [TaCpRX2(NPh)] along with the release of molecular hydrogen. Analogous reactions between the [{Ta(η5-C5Me5)X2}2(µ-H)2] derivatives and the cyclic diazo reagent benzo[c]cinnoline afford the biphenyl-bridged (phenylimido)tantalum complexes [{Ta(η5-C5Me5)X2}2(µ-NC6H4C6H4N)] along with the release of molecular hydrogen. When the compounds [(TaCpRX2)2(µ-H)2] (CpR = η5-C5H4SiMe3, η5-C5HMe4; X = Cl, Br) were employed, we were able to trap the side-on-bound diazo derivatives [(TaCpRX)2{µ-(η2,η2-NC6H4C6H4N)}] (CpR = η5-C5H4SiMe3, η5-C5HMe4; X = Cl, Br) as intermediates in the N═N bond cleavage process. DFT calculations provide insights into the N═N cleavage mechanism, in which the ditantalum(IV) fragment can promote two-electron reductions of the N═N bond at two different metal-metal bond splitting stages.

A Bridging bis-Allyl Titanium Complex: Mechanistic Insights into the Electronic Structure and Reactivity.

Treatment of the dinuclear compound [{Ti(η5-C5Me5)Cl2}2(µ-O)] with allylmagnesium chloride provides the formation of the allyltitanium(III) derivative [{Ti(η5-C5Me5)(µ-C3H5)}2(µ-O)] (1), structurally identified by single-crystal X-ray analysis. Density functional theory (DFT) calculations confirm that the electronic structure of 1 is a singlet state, and the molecular orbital analysis, along with the short Ti-Ti distance, reveal the presence of a metal-metal single bond between the two Ti(III) centers. Complex 1 reacts rapidly with organic azides, RN3 (R = Ph, SiMe3), to yield the allyl µ-imido derivatives [{Ti(η5-C5Me5)(CH2CH═CH2)2}2(µ-NR)(µ-O)] [R = Ph(2), SiMe3(3)] along with molecular nitrogen release. Reaction of 2 and 3 with H2 leads to the µ-imido propyl species [{Ti(η5-C5Me5)(CH2CH2CH3)2}2(µ-NR)(µ-O)] [R = Ph(4), SiMe3(5)]. Theoretical calculations were used to gain insight into the hydrogenation mechanism of complex 3 and rationalize the lower reactivity of 2. Initially, the µ-imido bridging group in these complexes activates the H2 molecule via addition to the Ti-N bonds. Subsequently, the titanium hydride intermediates induce a change in hapticity of the allyl ligands, and the nucleophilic attack of the hydride to the allyl groups leads to metallacyclopropane intermediates. Finally, the proton transfer from the amido group to the metallacyclopropane moieties affords the propyl complexes 4 and 5.