Reaction mechanism and kinetics for ammonia synthesis on the Fe(111) Surface
Abstract
The Haber-Bosch industrial process for synthesis of ammonia (NH_3) from hydrogen and nitrogen produces the millions of tons of ammonia gas annually needed to produce nitrates for fertilizers required to feed the earth's growing populations. This process has been optimized extensively, but it still uses enormous amounts of energy (2% of the world's supply), making it essential to dramatically improve its efficiency. To provide guidelines to accelerate this improvement, we used quantum mechanics to predict reaction mechanisms and kinetics for NH_3 synthesis on Fe(111)—the best Fe single crystal surface for NH_3 synthesis. We predicted the free energies of all reaction barriers for all steps in the mechanism and built these results into a kinetic Monte Carlo model for predicting steady state catalytic rates to compare with single-crystal experiments at 673 K and 20 atm. We find excellent agreement with a predicted turnover frequency (TOF) of 17.7 s^(-1) per 2 × 2 site (5.3 × 10^(-9) mol/cm^2/sec) compared to TOF = 10 s^(-1) per site from experiment.
Additional Information
© 2018 American Chemical Society. Received: December 25, 2017; Published: April 27, 2018. This work was supported by the U.S. Department of Energy (USDOE), Office of Energy Efficiency and Renewable Energy (EERE), Advanced Manufacturing Office Next Generation R&D Projects under contract no. DE-AC07-05ID14517 (program manager Dickson Ozokwelu, in collaboration with Idaho National Laboratories, Rebecca Fushimi). This project was initiated with modest support from the U.S. Department of Energy's Advanced Research Projects Agency-Energy (ARPA-E) under contract number DE-AR0000552 (Patrick McGrath) aimed at examining novel methods for accelerating catalytic reactions. A.F. gratefully acknowledges financial support from a Short-Term Mission (STM) funded by Italian Consiglio Nazionale delle Ricerche (CNR). We thank Dr. Tao Cheng for help with the ER reaction for N2 dissociation. Many of the calculations were carried out on a GPU-cluster provided by DURIP (Cliff Bedford, program manager). This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562. Author Contributions: J.Q. and Q.A. contributed equally to this work. The authors declare no competing financial interest.Attached Files
Supplemental Material - ja7b13409_si_001.xlsx
Supplemental Material - ja7b13409_si_002.pdf
Files
Name | Size | Download all |
---|---|---|
md5:b1ee0b049b3096ba6707d43a54d93d61
|
1.3 MB | Preview Download |
md5:10e0e8ab76e367fd71c13ffff57b44e9
|
56.4 kB | Download |
Additional details
- Eprint ID
- 86078
- Resolver ID
- CaltechAUTHORS:20180427-105456274
- Department of Energy (DOE)
- DE-AC07-05ID14517
- Advanced Research Projects Agency-Energy (ARPA-E)
- DE-AR0000552
- Consiglio Nazionale delle Ricerche (CNR)
- NSF
- ACI-1548562
- Created
-
2018-04-27Created from EPrint's datestamp field
- Updated
-
2021-11-15Created from EPrint's last_modified field
- Other Numbering System Name
- WAG
- Other Numbering System Identifier
- 1282