Wa Ode Nur Fitriah Rajaelo, Harry H. Halim, M. Fadhlan Anshor, Yuki Yamada, Syahrul Ramadhani, Fikri A. Mubarok, Septia E. M. Putra, Atthar L. Ivansyah, Triati D. K. Wungu, Suprijadi, Yoshitada Morikawa, Fahdzi Muttaqien
The discrepancy in the catalytic behaviors in ultrahigh vacuum vs realistic gas environments, known as “pressure gap”, has been a long-standing challenge in the field of surface science. The challenge lies in the limited knowledge of the surface morphology that might undergo reconstruction as a result of exposure to gas and adsorbates. In this study, we address the challenge by applying a set of multiscale simulations that included density functional theory (DFT), machine learning-based molecular dynamics (MLMD), and microkinetic modeling (MKM) simulations, taking the case of hydrogenation of CO2 into formate on the Cu(100) surface. MLMD simulations allow direct observation of hydrogen-induced reconstruction of the Cu(100) surface in the form of a rotated and line-shifted structure even at room temperature. This insight into the surface dynamics is used to provide more realistic structures of the active sites under experimental conditions. Our calculations reveal that this reconstruction significantly influences the CO2 hydrogenation rate on the Cu(100) surface. In particular, the line-shifted reconstruction on the Cu(100) surface brings the energy barrier and the hydrogenation rate much closer to those on the Cu(111) surface, matching the experimentally observed insensitivity of the Cu surface facets to the CO2 hydrogenation rate. On the other hand, the hydrogenation on the pristine Cu(100) surface exhibits a lower barrier compared to the Cu(111) surface, resulting in significant overestimation of the reaction rate compared to the experiment. This study demonstrates that adopting more realistic active sites from in situ observation of simulated catalytic reactions paves the way for a better connection between theory and experiment in the catalysis field. © 2026 The Authors. Published by American Chemical Society
Master Program in Computational Science, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Jl. Ganesa 10, Bandung, 40132, Indonesia; Department of Precision Engineering, Graduate School of Engineering, The University of Osaka, 2-1 Yamada-Oka, Osaka, Suita, 565-0871, Japan; Innovative Catalysis Science Division, Institute for Open and Transdisciplinary Research Initiatives (ICS-OTRI), The University of Osaka, Osaka, Suita, 565-0871, Japan; Engineering Physics, Institut Teknologi Sumatera, Jl. Terusan Ryacudu, Way Hui, Jati Agung, South Lampung, 35365, Indonesia; Faculty of Mathematics and Natural Sciences, Analytic Chemistry Research Group, Institut Teknologi Bandung, Jl. Ganesa 10, Bandung, 40132, Indonesia; Nuclear and Biophysics Research Group, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Jl. Ganesa 10, Bandung, 40132, Indonesia; Research Center for Nanoscience and Nanotechnology, Institut Teknologi Bandung, Jl. Ganesha 10, Bandung, 40132, Indonesia; Instrumentation and Computational Physics Research Group, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Jl. Ganesa 10, Bandung, 40132, Indonesia; Research Center for Precision Engineering, Graduate School of Engineering, The University of Osaka, 2-1 Yamadaoka, Osaka, Suita, 565-0871, Japan