Evolution of Ferromagnetic and Antiferromagnetic States in Iron Nitride Clusters FenN and FenN2 (n = 1-10)

Quantum hardware calculations of the activation and dissociation of nitrogen on iron clusters and surfaces

Catalytic processes are the cornerstone of chemical industry, and catalytic conversion of nitrogen to ammonia remains one of the largest industrial processes implemented. Rational design of catalysts and catalytic reactions largely depends on approximate computational chemistry methods, such as density functional theory, which, however, suffer from limited accuracy, especially for strongly-correlated materials. Rigorous ab initio methods which account for static and dynamic electron correlation, while arbitrarily accurate for small systems, are generally too expensive to be applied to modelling of catalytic cycles, due to prohibitive time and space computational complexity with respect to the size of the active space. Recent advances in quantum computing give hope for enabling access to accurate ab initio methods at scale. Herein, we present a prototype hybrid quantum-classical workflow for modeling chemical reactions on surfaces, applied to proof-of-concept models of activation and dissociation of nitrogen on small Fe clusters and a single-layer (221) iron surface. First, we determined the structures of species present in the catalytic cycle at DFT level and studied their electronic structure using CASSCF. We show that it is possible to decouple the half-filled Fe-3d band from the Fe-N and N-N bond orbitals, thereby reducing the active space significantly. Subsequently, we translated the CASSCF wavefunctions into corresponding qubit quantum states, using the Adaptive Variational Quantum Eigensolver, and estimated their energies using a state vector simulator, H1-1E quantum emulator and (for selected systems) H1-1 quantum computer. We demonstrated that if a sufficiently small active orbital space is chosen, ground state energies obtained with classical methods and with the quantum computer are in reasonable agreement. We argue that once quantum computing methods are scaled up so that larger active spaces are accessible, they can offer a tremendous practical advantage to the computational catalysis community.

Nitrogen Reduction to Ammonia on a Fe16 Nanocluster: A Computational Study of Catalysis

The sequence of elementary steps leading to reductive ammonia formation from N2 and H2 catalyzed by a Fe16 cluster is studied using generalized gradient approximation density functional theory and an all-electron basis set of triple-ζ quality. The computational methods are validated by comparison to experimental data such as binding energies where possible. First, the associative and dissociative attachment of N2 to Fe16 is considered, followed by exploration of the pathways leading to distal (Fe16-N-NH2) and enzymatic (NFe16-NH2) formation of an amino group. Next, the pathways leading to NH3 formation in both distal and enzymatic cases are examined. Two mechanisms for NH3 detachment have been discovered. An interesting peculiarity of the pathways is that they often proceed with total spin fluctuations, which are related to the rupture and formation of bonds on the surface of the catalyst over the course of the reactions. The reaction Fe16 + N2 + 2H2 → Fe16NH + NH3 is found to be exothermic by 1.02 eV (93.8 kJ/mol).

Plasma-Assisted Dinitrogen Activation on Small Cobalt Clusters: Co4 N9 + with Enhanced Stability

We have performed a study on the accommodation of nitrogen doping toward superatomic states of transition metal clusters. By reacting cobalt clusters with N₂ in the presence of plasma radiation, a large number of odd-nitrogen clusters were observed, typically Co₃ N_2m-1 ⁺ (m=1-5) and Co₄ N_2m-1 ⁺ (m=1-6) series, showing N≡N bond cleavage in the mild plasma atmosphere. Interestingly, the Co₃ N₇ ⁺ , Co₄ N₉ ⁺ , and Co₅ N₉ ⁺ clusters exhibit prominent mass abundances. First-principles calculation results elucidate the stability of the diverse cobalt nitride clusters and find unique stability of Co₄ N₉ ⁺ with a swallow-kite structure of which four coordinated N₂ molecules causes a significantly enlarged HOMO-LUMO gap, while the single N atom doping gives rise to superatomic states of 1S² 1P³ ||1D⁰ . We reveal an efficient dinitrogen activation strategy by reacting multiple N₂ molecules with cobalt clusters under a plasma atmosphere.