Discovery of Dramatically Improved Ammonia Synthesis Catalysts through Hierarchical High-Throughput Catalyst Screening of the Fe(211) Surface

Jon Fuller,Alessandro Fortunelli,William A. Goddard,Qi An
DOI: https://doi.org/10.1021/acs.chemmater.0c02701
IF: 10.508
2020-11-17
Chemistry of Materials
Abstract:In order to improve efficiency of ammonia synthesis using the Haber–Bosch (HB) process with Fe-based catalysts, we employed quantum mechanics (QM)-based hierarchical high-throughput catalyst screening (HHTCS) of 49 possible metal dopants. Here, we consider the Fe(211) surface (one of the two most active iron catalyst facets) to identify dopants that dramatically increase the turnover frequency (TOF) for HB synthesis. We found that under HB conditions, this surface reconstructs to form the Fe(211)R missing-row surface. Focusing on dopants with a strong preference for the subsurface site, we found that Co is the most promising candidate among the 49. We then examined the full reaction pathway on this Co-doped Fe(211)R surface, considering all 19 important 2 × 2 configurations and calculated the free-energy barriers (Δ<i>G</i><sup>∫</sup>) for all 12 important reaction steps. At 673 K and 20 atm, we find a decrease, δ(Δ<i>G</i><sup>∫</sup>) = −0.19 eV, in the overall reaction free-energy barrier for the Co-doped case. We then carried out kinetic Monte Carlo simulations for 60–120 min using 100 replicas with the full reaction path using rates from QM free-energy reaction barriers to predict that the TOF for the Co-doped surface increases by a factor of 2.8 with respect to the undoped Fe(211)R surface. Thus, the Co-doped Fe(211)R system could lower the extreme HB pressure of 200 atm to ∼40 atm at 773 K while maintaining the same TOF as that of undoped Fe(211)R. We conclude that Co dopants in the Fe catalyst could significantly improve the catalytic efficiency of ammonia synthesis under industrial conditions. This excellent performance of the Co-doped system is explained in terms of a surface spin analysis on the N<sub>2</sub>-bonded configurations that show how Co dopants shift the N<sub>2</sub> surface-binding mode. This demonstrates that metal surface spins can be used as quantitative descriptors to understand reaction energetics. This study demonstrates that the HHTCS kinetic analysis of the free-energy reaction path in terms of essential configurations can enable discovery of the salient barriers to overcome and best dopant candidates for further improvements.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.chemmater.0c02701?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.chemmater.0c02701</a>.Preferred positions of elements in the first and second layer for Fe-bcc(211)R; lowest-energy state for doping elements in Fe-bcc(211)R; convergence with respect to kMC steps for Co-doped and pure Fe-bcc(211)R surfaces; preferred positions of HHTCS elements; surface behavior of large elements; potential energy curves for hydrogenation transition states; potential energy curves for nitrogen dissociation transition states; and full reaction pathway for the undoped Fe-bcc(211)R surface (<a class="ext-link" href="/doi/suppl/10.1021/acs.chemmater.0c02701/suppl_file/cm0c02701_si_001.pdf">PDF</a>)Raw DFT data for HHTCS (<a class="ext-link" href="/doi/suppl/10.1021/acs.chemmater.0c02701/suppl_file/cm0c02701_si_002.xlsx">XLSX</a>)Raw data for the reaction energy path of the HB process on the Co-doped Fe-bcc(211)R surface (<a class="ext-link" href="/doi/suppl/10.1021/acs.chemmater.0c02701/suppl_file/cm0c02701_si_003.xlsx">XLSX</a>)This article has not yet been cited by other publications.
materials science, multidisciplinary,chemistry, physical
What problem does this paper attempt to address?