First-principles-based multiscale modelling of heterogeneous catalysis

Kinetic Monte Carlo Simulations of Low-Temperature NH3-SCR over Cu-Exchanged Chabazite

Cu-exchanged chabazite (Cu-CHA) is widely used for ammonia assisted selective catalytic reduction of nitrogen oxides (NH3-SCR). The Cu+ ions are at low temperatures solvated by NH3 forming mobile [Cu(NH3)2]+ complexes. The dynamic behaviour of the complexes is critical as O2 adsorption requires a pair of complexes to form a [Cu2(NH3)4O2]2+ peroxo-species over which NO couples with NH3. Here we introduce a first principles-based kinetic Monte Carlo approach to explore the effect of the Al-distribution on the reaction kinetics of NH3-SCR over Cu-CHA. The method allows us to scrutinize the interplay between the pairing of [Cu(NH3)2]+ complexes and the reaction landscape for the NH3-SCR reaction over the peroxo-complex. The Al-distribution affects the stability of the [Cu(NH3)2]+ pairs as well as the kinetic parameters of the SCR-reaction. The turn-over frequency is determined by the stability of the [Cu(NH3)2]+ pairs and the relative strength of NO and NH3 adsorption once a pair is present. The results establish the hierarchy of effects that influences the performance of Cu-CHA over NH3-SCR and provide a computational basis for further development of the Cu-CHA material.

Analytical approach to structural chemistry origins of mechanical, acoustical and thermal properties

Crystalline matters with periodically arranged atoms found wide applications in modern science and technology. To facilitate the design of new materials and the advancement of existing ones, accurate and efficient models without relying too much on known inputs for predicting the functionalities are essential. Here, we propose an analytical approach for such a purpose, with only the knowledge of the structural chemistry of crystals. Based on the electrostatic interaction between periodically arranged atoms, the 1st, 2nd and 3rd derivatives of interatomic potential, respectively, enable a prediction of ten kinds in total of mechanical, acoustical and thermal properties. Over a thousand measurements are collected from ∼500 literatures, this results in the symmetric mean percentage error (SMPE) within ±25% and the symmetric mean absolute percentage error (SMAPE) ranging from 22%∼74% across all properties predicted, which further enables a revelation of bond characteristics as the most important but implicit origin for functionalities.

Utilizing a synergistic strategy that combines electromagnetic and chemical enhancement to analyze the SERS effect of the Fe3O4@GO@Ag on PAHs detection

A comprehensive understanding of the enhancement mechanism of the substrate material is crucial to ensure the repeatability and functionality of SERS detection technology. Therefore, this study introduces a theoretical analysis method that integrates electromagnetic and chemical enhancement to achieve a comprehensive understanding of the SERS effect on the magnetic composite substrate. The visual model is employed in this study to comprehensively analyze and illustrate the electric field enhancement and optical effects of composite substrate materials. The study also elucidated the adsorption and charge transfer between the substrate material and target molecules. Based on this theory, Fe3O4@GO@Ag material was prepared and used to detect hydrophobic organic molecules such as polycyclic aromatic hydrocarbons (PAHs), with a concentration as low as 0.5 nM. This study comprehensively analyzed the SERS enhancement effect of the composite substrate for the first time, and prepared a magnetic composite substrate material for the detection of hydrophobic organic molecules, opening up a new avenue for theoretical guidance and experimental exploration in SERS detection and analysis.

Small basis set density functional theory method for cost-efficient, large-scale condensed matter simulations

We present an efficient first-principles based method geared toward reliably predicting the structures of solid materials across the Periodic Table. To this end, we use a density functional theory baseline with a compact, near-minimal min+s basis set, yielding low computational costs and memory demands. Since the use of such a small basis set leads to systematic errors in chemical bond lengths, we develop a linear pairwise correction, available for elements Z = 1-86 (excluding the lanthanide series), parameterized for use with the Perdew-Burke-Ernzerhof exchange-correlation functional. We demonstrate the reliability of this corrected approach for equilibrium volumes across the Periodic Table and the transferability to differently coordinated environments and multi-elemental crystals. We examine relative energies, forces, and stresses in geometry optimizations and molecular dynamics simulations.

Simulating Structural Dynamics of Metal Catalysts under Operative Conditions

Structural reconstructions of metal catalysts have been recognized as common phenomena during catalytic reactions, which play a key role in their activities in heterogeneous catalysis. Precisely identifying the structures under the operative conditions becomes a prerequisite to establish a reliable structure-activity relationship and further rationalize the design of metal catalysts. However, real-time capture of the structural variations of catalysts at the atomic level with high-temporal resolution is a grand challenge for present in situ characterizations. During the past decade, significant progress has been made in theory to couple the structures with the reaction conditions to reproduce the experimental observations and predict the adsorbate-induced changes of catalysts in composition, morphology, size, etc. Modeling the dynamic correlation between the structure and activity of the metal catalysts brings us advanced knowledge of heterogeneous catalysis and becomes indispensable for accurate evaluation of the performance of metal catalysts.

i-PI 3.0: A flexible and efficient framework for advanced atomistic simulations

Atomic-scale simulations have progressed tremendously over the past decade, largely thanks to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD, and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows for deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities.

Curing and Cross-Linking Processes in the Poly(3,3-bis-azidomethyl oxetane)-tetrahydrofuran/Toluene Diisocyanate/Trimethylolpropane System: A Density Functional Theory and Accelerated ReaxFF Molecular Dynamics Investigation

The physical and chemical properties of solid propellant are influenced by the composition and structure of the binder, with its network structure being formed through curing and cross-linking reactions. Therefore, understanding the mechanisms of these reactions is crucial. In this study, we investigated the curing and cross-linking mechanisms of poly(3,3-bis-azidomethyl oxetane)-tetrahydrofuran (PBT), toluene diisocyanate (TDI), and trimethylolpropane (TMP) using a combination of density functional theory (DFT) calculations and accelerated ReaxFF molecular dynamics (MD) simulations. DFT calculations revealed that the steric effect of the -CH3 group in TDI exerts a significant influence on the curing reaction between TDI and PBT. Additionally, in the cross-linking process, the energy barrier for TDI reacting with TMP was found to be much lower than that for TDI reacting with the PBT-TDI intermediate. Subsequently, we conducted competing reaction processes of TMP/TDI-PBT-TDI cross-linking and TDI-PBT-TDI self-cross-linking using accelerated MD simulations within the fitted ReaxFF framework. The results showed that the successful frequency of TMP/TDI-PBT-TDI cross-linking was substantially higher than that of TDI-PBT-TDI self-cross-linking, consistent with the energy barrier results from DFT calculations. These findings deepen our understanding of the curing and cross-linking mechanisms of the PBT system, providing valuable insights for the optimization and design of solid propellants.

Advances in the Dynamics of Adsorbate Diffusion on Metal Surfaces: Focus on Hydrogen and Oxygen

Adsorbates on metal surfaces are typically formed from the dissociative chemisorption of molecules occurring at gas-solid interfaces. These adsorbed species exhibit unique diffusion behaviors on metal surfaces, which are influenced by their translational energy. They play crucial roles in various fields, including heterogeneous catalysis and corrosion. This review examines recent theoretical advancements in understanding the diffusion dynamics of adsorbates on metal surfaces, with a specific emphasis on hydrogen and oxygen atoms. The diffusion processes of adsorbates on metal surfaces involve two energy transfer mechanisms: surface phonons and electron-hole pair excitations. This review also surveys new theoretical methods, including the characterization of the electron-hole pair excitation within electronic friction models, the acceleration of quantum chemistry calculations through machine learning, and the treatment of atomic nuclear motion from both quantum mechanical and classical perspectives. Furthermore, this review offers valuable insights into how energy transfer, nuclear quantum effects, supercell sizes, and the topography of potential energy surfaces impact the diffusion behavior of hydrogen and oxygen species on metal surfaces. Lastly, some preliminary research proposals are presented.

On-the-fly kinetic Monte Carlo simulations with neural network potentials for surface diffusion and reaction

We develop an adaptive scheme in the kinetic Monte Carlo simulations, where the adsorption and activation energies of all elementary steps, including the effects of other adsorbates, are evaluated "on-the-fly" by employing the neural network potentials. The configurations and energies evaluated during the simulations are stored for reuse when the same configurations are sampled in a later step. The present scheme is applied to hydrogen adsorption and diffusion on the Pd(111) and Pt(111) surfaces and the CO oxidation reaction on the Pt(111) surface. The effects of interactions between adsorbates, i.e., adsorbate-adsorbate lateral interactions, are examined in detail by comparing the simulations without considering lateral interactions. This study demonstrates the importance of lateral interactions in surface diffusion and reactions and the potential of our scheme for applications in a wide variety of heterogeneous catalytic reactions.

Design Principle of Molybdenum-Based Metal Nitrides for Lattice Nitrogen-Mediated Ammonia Production

Chemical looping ammonia synthesis (CLAS) is a promising technology for reducing the high energy consumption of the conventional ammonia synthesis process. However, the comprehensive understanding of reaction mechanisms and rational design of novel nitrogen carriers has not been achieved due to the high complexity of catalyst structures and the unrevealed relationship between electronic structure and intrinsic activity. Herein, we propose a multistage strategy to establish the connection between catalyst intrinsic activity and microscopic electronic structure fingerprints using density functional theory computational energetics as bridges and apply it to the rational design of metal nitride catalysts for lattice nitrogen-mediated ammonia production. Molybdenum-based nitride catalysts with well-defined structures are employed as prototypes to elucidate the decoupled effects of electronic and geometrical features. The electron-transfer and spin polarization characteristics of the magnetic metals are constructed as descriptors to disclose the atomic-scale causes of intrinsic activity. Based on this design strategy, it is demonstrated that Ni3Mo3N catalysts possess the highest lattice nitrogen-mediated ammonia synthesis activity. This work reveals the structure-activity relationship of metal nitrides for CLAS and provides a multistage perspective on catalyst rational design.

Spontaneous Charge Separation at the Metal-Water Interface

Reactions at the metal-water interface are essential in a range of fundamental and technological processes. Using Density Functional Theory calculations, we demonstrate that water substantially affects the adsorption of H and O2 on Cu(111), Ag(111), Au(111), Pd(111) and Pt(111). In water, H is found to undergo a spontaneous charge separation, where a proton desorbs to the water solution while an electron is donated to the surface. The reaction is exothermic over Au and Pt and associated with low barriers. The process is facile also over Pd, albeit slightly endothermic. For O2, water is found to increase the metal-to-adsorbate charge transfer, enhancing the adsorption energy and O-O bond length as compared to the adsorption in the absence of water. The magnitudes of the effects are system dependent, which implies that calculations should treat water explicitly. The results elucidate previous experimental results and highlights the importance of charge-transfer effects at the metal-water interface; both to describe the potential energy landscape, and to account for alternative reaction routes in the presence of water.

Three-Center Tight-Binding Together with Multipolar Auxiliary Functions

We present an ab initio tight-binding method that allows to improve on the effective potential and minimal basis approximations employed in semiempirical calculations. Three-center expansions are used to evaluate the zeroth-order Hamiltonian matrix elements and repulsive energy terms in the spirit of the Horsfield method. Self-consistency is handled by expanding atomic orbital products in an auxiliary basis following the work of Giese and York, combined with a two-center expansion of the exchange-correlation kernels. Together with nonminimal main basis sets (double-ζ plus polarization), we show that the resulting method trades a modest amount of accuracy for a significant gain in speed, compared to that of numerical atomic orbital density functional theory, in calculations on small molecules, bulk compounds, and metal nanoclusters.

Discovering High Entropy Alloy Electrocatalysts in Vast Composition Spaces with Multiobjective Optimization

High entropy alloys (HEAs) are a highly promising class of materials for electrocatalysis as their unique active site distributions break the scaling relations that limit the activity of conventional transition metal catalysts. Existing Bayesian optimization (BO)-based virtual screening approaches focus on catalytic activity as the sole objective and correspondingly tend to identify promising materials that are unlikely to be entropically stabilized. Here, we overcome this limitation with a multiobjective BO framework for HEAs that simultaneously targets activity, cost-effectiveness, and entropic stabilization. With diversity-guided batch selection further boosting its data efficiency, the framework readily identifies numerous promising candidates for the oxygen reduction reaction that strike the balance between all three objectives in hitherto unchartered HEA design spaces comprising up to 10 elements.

Computational chemistry for water-splitting electrocatalysis

Electrocatalytic water splitting driven by renewable electricity has attracted great interest in recent years for producing hydrogen with high-purity. However, the practical applications of this technology are limited by the development of electrocatalysts with high activity, low cost, and long durability. In the search for new electrocatalysts, computational chemistry has made outstanding contributions by providing fundamental laws that govern the electron behavior and enabling predictions of electrocatalyst performance. This review delves into theoretical studies on electrochemical water-splitting processes. Firstly, we introduce the fundamentals of electrochemical water electrolysis and subsequently discuss the current advancements in computational methods and models for electrocatalytic water splitting. Additionally, a comprehensive overview of benchmark descriptors is provided to aid in understanding intrinsic catalytic performance for water-splitting electrocatalysts. Finally, we critically evaluate the remaining challenges within this field.

What Makes On-Chip Microdevices Stand Out in Electrocatalysis?

Clean and sustainable energy conversion and storage through electrochemistry shows great promise as an alternative to traditional fuel or fossil-consumption energy systems. With regards to practical and high-efficient electrochemistry application, the rational design of active sites and the accurate description of mechanism remain a challenge. Toward this end, in this Perspective, a unique on-chip micro/nano device coupling nanofabrication and low-dimensional electrochemical materials is presented, in which material structure analysis, field-effect regulation, in situ monitoring, and simulation modeling are highlighted. The critical mechanisms that influence electrochemical response are discussed, and how on-chip micro/nano device distinguishes itself is emphasized. The key challenges and opportunities of on-chip electrochemical platforms are also provided through the Perspective.

A Data and DFT-Driven Framework for Predicting the Microstructure of Submonolayer Inverse Metal Oxide on Metal Catalysts

Metal oxides on metal (inverse) catalysts can selectively drive many important reactions. However, understanding the active site under experimentally relevant conditions is lacking. Herein, we introduce a computational framework for predicting atomic models of stable inverse catalysts and demonstrate it for WOx on Pt(553) and a Pt79 nanoparticle at variable WOx coverages. An evolutionary algorithm identifies a small (5%) subset of promising atomic configurations on which DFT simulations are performed. We predict a maximum coverage of ∼50% WOx on Pt(553), consisting of small clusters (tetramers and pentamers), which preferentially reside on the terrace, with their oxygen atoms interacting with the Pt step sites. Consistently, WOx does not lie on curved and undercoordinated metal sites of Pt nanoparticles. The oxide clusters prefer a partially reduced oxidation state. Theoretical EXAFS spectra for select configurations provide insights into interpreting experimental spectra of inverse catalysts. The framework applies to other catalysts.

Reference-Quality Free Energy Barriers in Catalysis from Machine Learning Thermodynamic Perturbation Theory

For the first time, we report calculations of the free energies of activation of cracking and isomerization reactions of alkenes that combine several different electronic structure methods with molecular dynamics simulations. We demonstrate that the use of a high level of theory (here Random Phase Approximation-RPA) is necessary to bridge the gap between experimental and computed values. These transformations, catalyzed by zeolites and proceeding via cationic intermediates and transition states, are building blocks of many chemical transformations for valorization of long chain paraffins originating, e.g., from plastic waste, vegetable oils, Fischer-Tropsch waxes or crude oils. Compared with the free energy barriers computed at the PBE+D2 production level of theory via constrained ab initio molecular dynamics, the barriers computed at the RPA level by the application of Machine Learning thermodynamic Perturbation Theory (MLPT) show a significant decrease for isomerization reaction and an increase of a similar magnitude for cracking, yielding an unprecedented agreement with the results obtained by experiments and kinetic modeling.

Automated MUltiscale simulation environment

Multiscale techniques integrating detailed atomistic information on materials and reactions to predict the performance of heterogeneous catalytic full-scale reactors have been suggested but lack seamless implementation. The largest challenges in the multiscale modeling of reactors can be grouped into two main categories: catalytic complexity and the difference between time and length scales of chemical and transport phenomena. Here we introduce the Automated MUltiscale Simulation Environment AMUSE, a workflow that starts from Density Functional Theory (DFT) data, automates the analysis of the reaction networks through graph theory, prepares it for microkinetic modeling, and subsequently integrates the results into a standard open-source Computational Fluid Dynamics (CFD) code. We demonstrate the capabilities of AMUSE by applying it to the unimolecular iso-propanol dehydrogenation reaction and then, increasing the complexity, to the pre-commercial Pd/In2O3 catalyst employed for the CO2 hydrogenation to methanol. The results show that AMUSE allows the computational investigation of heterogeneous catalytic reactions in a comprehensive way, providing essential information for catalyst design from the atomistic to the reactor scale level.

Visualization of Self-Assembly and Hydration of a β-Hairpin through Integrated Small and Wide-Angle Neutron Scattering

Fundamental understanding of the structure and assembly of nanoscale building blocks is crucial for the development of novel biomaterials with defined architectures and function. However, accessing self-consistent structural information across multiple length scales is challenging. This limits opportunities to exploit atomic scale interactions to achieve emergent macroscale properties. In this work we present an integrative small- and wide-angle neutron scattering approach coupled with computational modeling to reveal the multiscale structure of hierarchically self-assembled β hairpins in aqueous solution across 4 orders of magnitude in length scale from 0.1 Å to 300 nm. Our results demonstrate the power of this self-consistent cross-length scale approach and allows us to model both the large-scale self-assembly and small-scale hairpin hydration of the model β hairpin CLN025. Using this combination of techniques, we map the hydrophobic/hydrophilic character of this model self-assembled biomolecular surface with atomic resolution. These results have important implications for the multiscale investigation of aqueous peptides and proteins, for the prediction of ligand binding and molecular associations for drug design, and for understanding the self-assembly of peptides and proteins for functional biomaterials.

Estimating Free Energy Barriers for Heterogeneous Catalytic Reactions with Machine Learning Potentials and Umbrella Integration

Predicting the rate constants of elementary reaction steps is key for the computational modeling of catalytic processes. Within transition state theory (TST), this requires an accurate estimation of the corresponding free energy barriers. While sophisticated methods for estimating free energy differences exist, these typically require extensive (biased) molecular dynamics simulations that are computationally prohibitive with the first-principles electronic structure methods that are typically used in catalysis research. In this contribution, we show that machine-learning (ML) interatomic potentials can be trained in an automated iterative workflow to perform such free energy calculations at a much reduced computational cost as compared to a direct density functional theory (DFT) based evaluation. For the decomposition of CHO on Rh(111), we find that thermal effects are substantial and lead to a decrease in the free energy barrier, which can be vanishingly small, depending on the DFT functional used. This is in stark contrast to previously reported estimates based on a harmonic TST approximation, which predicted an increase in the barrier at elevated temperatures. Since CHO is the reactant of the putative rate limiting reaction step in syngas conversion on Rh(111) and essential for the selectivity toward oxygenates containing multiple carbon atoms (C2+ oxygenates), our results call into question the reported mechanism established by microkinetic models.

Automated, Consistent, and Even-Handed Selection of Active Orbital Spaces for Quantum Embedding

A widely used strategy to reduce the computational cost of quantum-chemical calculations is to partition the system into an active subsystem, which is the focus of the computational efforts, and an environment that is treated at a lower computational level. The system partitioning is mostly based on localized molecular orbitals. When reaction paths or energy differences are to be calculated, it is crucial to keep the orbital space consistent for all structures. Inconsistencies in orbital space can lead to unpredictable errors on the potential energy surface. While successful strategies to ensure this consistency have been established for organic and even metal-organic systems, these methods often fail for metal clusters or nanoparticles with a high density of near-degenerate and delocalized molecular orbitals. However, such systems are highly relevant for catalysis. Accurate yet feasible quantum-mechanical ab initio calculations are therefore highly desired. In this work, we present an approach based on the subsystem projected atomic orbital decomposition algorithm that allows us to ensure automated and consistent partitioning even for systems with delocalized and near-degenerate molecular orbitals and demonstrate the validity of this method for the binding energies of small molecules on transition-metal clusters.

Advanced electrocatalytic redox processes for environmental remediation of halogenated organic water pollutants

Halogenated organic compounds are widespread, and decades of heavy use have resulted in global bioaccumulation and contamination of the environment, including water sources. Here, we introduce the most common halogenated organic water pollutants, their classification by type of halogen (fluorine, chlorine, or bromine), important policies and regulations, main applications, and environmental and human health risks. Remediation techniques are outlined with particular emphasis on carbon-halogen bond strengths. Aqueous advanced redox processes are discussed, highlighting mechanistic details, including electrochemical oxidations and reductions of the water-oxygen system, and thermodynamic potentials, protonation states, and lifetimes of radicals and reactive oxygen species in aqueous electrolytes at different pH conditions. The state of the art of aqueous advanced redox processes for brominated, chlorinated, and fluorinated organic compounds is presented, along with reported mechanisms for aqueous destruction of select PFAS (per- and polyfluoroalkyl substances). Future research directions for aqueous electrocatalytic destruction of organohalogens are identified, emphasizing the crucial need for developing a quantitative mechanistic understanding of degradation pathways, the improvement of analytical detection methods for organohalogens and transient species during advanced redox processes, and the development of new catalysts and processes that are globally scalable.

Automated Generation of Microkinetics for Heterogeneously Catalyzed Reactions Considering Correlated Uncertainties

The study presents an ab-initio based framework for the automated construction of microkinetic mechanisms considering correlated uncertainties in all energetic parameters and estimation routines. 2000 unique microkinetic models were generated within the uncertainty space of the BEEF-vdW functional for the oxidation reactions of representative exhaust gas emissions from stoichiometric combustion engines over Pt(111) and compared to experiments through multiscale modeling. The ensemble of simulations stresses the importance of considering uncertainties. Within this set of first-principles-based models, it is possible to identify a microkinetic mechanism that agrees with experimental data. This mechanism can be traced back to a single exchange-correlation functional, and it suggests that Pt(111) could be the active site for the oxidation of light hydrocarbons. The study provides a universal framework for the automated construction of reaction mechanisms with correlated uncertainty quantification, enabling a DFT-constrained microkinetic model optimization for other heterogeneously catalyzed systems.

CO2 Hydrogenation to Gasoline and Aromatics: Mechanistic and Predictive Insights from DFT, DRIFTS and Machine Learning

The emission of CO2 from fossil fuels is the largest driver of global climate change. To realize the target of a carbon-neutrality by 2050, CO2 capture and utilization is crucial. The efficient conversion of CO2 to C5+ gasoline and aromatics remains elusive mainly due to CO2 thermodynamic stability and the high energy barrier of the C-C coupling step. Herein, advances in mechanistic understanding via Diffuse Reflectance Infrared Fourier Transform Spectroscopy (DRIFTS), density functional theory (DFT), and microkinetic modeling are discussed. It further emphasizes the power of machine learning (ML) to accelerate the search for optimal catalysts. A significant effort has been invested into this field of research with volumes of experimental and characterization data, this study discusses how they can be used as input features for machine learning prediction in a bid to better understand catalytic properties capable of accelerating breakthroughs in the process.

Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations

Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.

A robust sparse identification method for nonlinear dynamic systems affected by non-stationary noise

In the field of science and engineering, identifying the nonlinear dynamics of systems from data is a significant yet challenging task. In practice, the collected data are often contaminated by noise, which often severely reduce the accuracy of the identification results. To address the issue of inaccurate identification induced by non-stationary noise in data, this paper proposes a method called weighted ℓ1-regularized and insensitive loss function-based sparse identification of dynamics. Specifically, the robust identification problem is formulated using a sparse identification mathematical model that takes into account the presence of non-stationary noise in a quantitative manner. Then, a novel weighted ℓ1-regularized and insensitive loss function is proposed to account for the nature of non-stationary noise. Compared to traditional loss functions like least squares and least absolute deviation, the proposed method can mitigate the adverse effects of non-stationary noise and better promote the sparsity of results, thereby enhancing the accuracy of identification. Third, to overcome the non-smooth nature of the objective function induced by the inclusion of loss and regularization terms, a smooth approximation of the non-smooth objective function is presented, and the alternating direction multiplier method is utilized to develop an efficient optimization algorithm. Finally, the robustness of the proposed method is verified by extensive experiments under different types of nonlinear dynamical systems. Compared to some state-of-the-art methods, the proposed method achieves better identification accuracy.

Iterative multiscale and multi-physics computations for operando catalyst nanostructure elucidation and kinetic modeling

Modern heterogeneous catalysis has benefitted immensely from computational predictions of catalyst structure and its evolution under reaction conditions, first-principles mechanistic investigations, and detailed kinetic modeling, which are rungs on a multiscale workflow. Establishing connections across these rungs and integration with experiments have been challenging. Here, operando catalyst structure prediction techniques using density functional theory simulations and ab initio thermodynamics calculations, molecular dynamics, and machine learning techniques are presented. Surface structure characterization by computational spectroscopic and machine learning techniques is then discussed. Hierarchical approaches in kinetic parameter estimation involving semi-empirical, data-driven, and first-principles calculations and detailed kinetic modeling via mean-field microkinetic modeling and kinetic Monte Carlo simulations are discussed along with methods and the need for uncertainty quantification. With these as the background, this article proposes a bottom-up hierarchical and closed loop modeling framework incorporating consistency checks and iterative refinements at each level and across levels.

Accurate quantification of the stability of the perylene-tetracarboxylic dianhydride on Au(111) molecule-surface interface

Studying inorganic/organic hybrid systems is a stepping stone towards the design of increasingly complex interfaces. A predictive understanding requires robust experimental and theoretical tools to foster trust in the obtained results. The adsorption energy is particularly challenging in this respect, since experimental methods are scarce and the results have large uncertainties even for the most widely studied systems. Here we combine temperature-programmed desorption (TPD), single-molecule atomic force microscopy (AFM), and nonlocal density-functional theory (DFT) calculations, to accurately characterize the stability of a widely studied interface consisting of perylene-tetracarboxylic dianhydride (PTCDA) molecules on Au(111). This network of methods lets us firmly establish the adsorption energy of PTCDA/Au(111) via TPD (1.74 ± 0.10 eV) and single-molecule AFM (2.00 ± 0.25 eV) experiments which agree within error bars, exemplifying how implicit replicability in a research design can benefit the investigation of complex materials properties.

Designing main-group catalysts for low-temperature methane combustion by ozone

The catalytic combustion of methane at a low temperature is becoming increasingly key to controlling unburned CH₄ emissions from natural gas vehicles and power plants, although the low activity of benchmark platinum-group-metal catalysts hinders its broad application. Based on automated reaction route mapping, we explore main-group elements catalysts containing Si and Al for low-temperature CH₄ combustion with ozone. Computational screening of the active site predicts that strong Brønsted acid sites are promising for methane combustion. We experimentally demonstrate that catalysts containing strong Bronsted acid sites exhibit improved CH₄ conversion at 250 °C, correlating with the theoretical predictions. The main-group catalyst (proton-type beta zeolite) delivered a reaction rate that is 442 times higher than that of a benchmark catalyst (5 wt% Pd-loaded Al₂O₃) at 190 °C and exhibits higher tolerance to steam and SO₂. Our strategy demonstrates the rational design of earth-abundant catalysts based on automated reaction route mapping.

Site Communication in Direct Formation of H2O2 over Single-Atom Pd@Au Nanoparticles

Single atom alloy catalysts offer possibilities to obtain turnover frequencies and selectivities unattainable by their monometallic counterparts. One example is direct formation of H₂O₂ from O₂ and H₂ over Pd embedded in Au hosts. Here, a first-principles-based kinetic Monte Carlo approach is developed to investigate the catalytic performance of Pd embedded in Au nanoparticles in an aqueous solution. The simulations reveal an efficient site separation where Pd monomers act as active centers for H₂ dissociation, whereas H₂O₂ is formed over undercoordinated Au sites. After dissociation, atomic H may undergo an exothermic redox reaction, forming a hydronium ion in the solution and a negative charge on the surface. H₂O₂ is preferably formed from reactions between dissolved H⁺ and oxygen species on the Au surface. The simulations show that tuning the nanoparticle composition and reaction conditions can enhance the selectivity toward H₂O₂. The outlined approach is general and applicable for a range of different hydrogenation reactions over single atom alloy nanoparticles.

An automated reaction route mapping for the reaction of NO and active species on Ag4 clusters in zeolites

A computational investigation of the catalytic reaction on multinuclear sites is very challenging. Here, using an automated reaction route mapping method, the single-component artificial force induced reaction (SC-AFIR) algorithm, the catalytic reaction of NO and OH/OOH species over the Ag₄²⁺ cluster in a zeolite is investigated. The results of the reaction route mapping for H₂ + O₂ reveal that OH and OOH species are formed over the Ag₄²⁺ cluster via an activation barrier lower than that of OH formation from H₂O dissociation. Then, reaction route mapping is performed to examine the reactivity of the OH and OOH species with NO molecules over the Ag₄²⁺ cluster, resulting in the facile reaction path of HONO formation. With the aid of the automated reaction route mapping, the promotion effect of H₂ addition on the SCR reaction was computationally proposed (boosting the formation of OH and OOH species). In addition, the present study emphasizes that automated reaction route mapping is a powerful tool to elucidate the complicated reaction pathway on multi-nuclear clusters.

Data-Centric Heterogeneous Catalysis: Identifying Rules and Materials Genes of Alkane Selective Oxidation

Artificial intelligence (AI) can accelerate catalyst design by identifying key physicochemical descriptive parameters correlated with the underlying processes triggering, favoring, or hindering the performance. In analogy to genes in biology, these parameters might be called "materials genes" of heterogeneous catalysis. However, widely used AI methods require big data, and only the smallest part of the available data meets the quality requirement for data-efficient AI. Here, we use rigorous experimental procedures, designed to consistently take into account the kinetics of the catalyst active states formation, to measure 55 physicochemical parameters as well as the reactivity of 12 catalysts toward ethane, propane, and n-butane oxidation reactions. These materials are based on vanadium or manganese redox-active elements and present diverse phase compositions, crystallinities, and catalytic behaviors. By applying the sure-independence-screening-and-sparsifying-operator symbolic-regression approach to the consistent data set, we identify nonlinear property-function relationships depending on several key parameters and reflecting the intricate interplay of processes that govern the formation of olefins and oxygenates: local transport, site isolation, surface redox activity, adsorption, and the material dynamical restructuring under reaction conditions. These processes are captured by parameters derived from N₂ adsorption, X-ray photoelectron spectroscopy (XPS), and near-ambient-pressure in situ XPS. The data-centric approach indicates the most relevant characterization techniques to be used for catalyst design and provides "rules" on how the catalyst properties may be tuned in order to achieve the desired performance.

Graphene at Liquid Copper Catalysts: Atomic-Scale Agreement of Experimental and First-Principles Adsorption Height

Liquid metal catalysts have recently attracted attention for synthesizing high-quality 2D materials facilitated via the catalysts' perfectly smooth surface. However, the microscopic catalytic processes occurring at the surface are still largely unclear because liquid metals escape the accessibility of traditional experimental and computational surface science approaches. Hence, numerous controversies are found regarding different applications, with graphene (Gr) growth on liquid copper (Cu) as a prominent prototype. In this work, novel in situ and in silico techniques are employed to achieve an atomic-level characterization of the graphene adsorption height above liquid Cu, reaching quantitative agreement within 0.1 Å between experiment and theory. The results are obtained via in situ synchrotron X-ray reflectivity (XRR) measurements over wide-range q-vectors and large-scale molecular dynamics simulations based on efficient machine-learning (ML) potentials trained to first-principles density functional theory (DFT) data. The computational insight is demonstrated to be robust against inherent DFT errors and reveals the nature of graphene binding to be highly comparable at liquid Cu and solid Cu(111). Transporting the predictive first-principles quality via ML potentials to the scales required for liquid metal catalysis thus provides a powerful approach to reach microscopic understanding, analogous to the established computational approaches for catalysis at solid surfaces.

Mesokinetics as a Tool Bridging the Microscopic-to-Macroscopic Transition to Rationalize Catalyst Design

Heterogeneous catalysis is the workhorse of the chemical industry, and a heterogeneous catalyst possesses numerous active sites working together to drive the conversion of reactants to desirable products. Over the decades, much focus has been placed on identifying the factors affecting the active sites to gain deep insights into the structure-performance relationship, which in turn guides the design and preparation of more active, selective, and stable catalysts. However, the molecular-level interplay between active sites and catalytic function still remains qualitative or semiquantitative, ascribed to the difficulty and uncertainty in elucidating the nature of active sites for its controllable manipulation. Hence, bridging the microscopic properties of active sites and the macroscopic catalytic performance, that is, microscopic-to-macroscopic transition, to afford a quantitative description is intriguing yet challenging, and progress toward this promises to revolutionize catalyst design and preparation.In this Account, we propose mesokinetics modeling, for the first time enabling a quantitative description of active site characteristics and the related mechanistic information, as a versatile tool to guide rational catalyst design. Exemplified by a pseudo-zero-order reaction, the kinetics derivation from the Pt particle size-sensitive catalytic activity and size-insensitive activation energy suggests only one type of surface site as the dominant active site, in which the Pt(111) with almost unchanged turnover frequency (TOF₁₁₁) is further identified as the dominating active site. Such a method has been extended to identify and quantify the number (N_i) of active sites for various thermo-, electro-, and photocatalysts in chemical synthesis, hydrogen generation, environment application, etc. Then, the kinetics derivation from the kinetic compensation effects suggests a thermodynamic balance between the activation entropy and enthalpy, which exhibit linear dependences on Pt charge. Accordingly, the Pt charge can serve as a catalytic descriptor for its quantitative determination of TOF_i. This strategy has been further applied to Pt-catalyzed CO oxidation with nonzero-order reaction characteristic by taking the site coverages of surface species into consideration.Hence, substituting the above statistical correlations of N_i and TOF_i into the rate equation R = ∑N_i × TOF_i offers the mesokinetics model, which can precisely predict catalytic function and screen catalysts. Finally, based on the disentanglement of the factors underlying Pt electronic structures, a de novo strategy, from the interfacial charge distribution to reaction mechanism, kinetics, and thermodynamics parameters of the rate-determining step, and ultimately catalytic performance, is developed to map the unified mechanistic and kinetics picture of reaction. Overall, the mesokinetics not only demonstrates much potential to elucidate the quantitative interplay between active sites and catalytic activity but also provides a new research direction in kinetics analysis to rationalize catalyst design.

Graph-Driven Reaction Discovery: Progress, Challenges, and Future Opportunities

Graph-based descriptors, such as bond-order matrices and adjacency matrices, offer a simple and compact way of categorizing molecular structures; furthermore, such descriptors can be readily used to catalog chemical reactions (i.e., bond-making and -breaking). As such, a number of graph-based methodologies have been developed with the goal of automating the process of generating chemical reaction network models describing the possible mechanistic chemistry in a given set of reactant species. Here, we outline the evolution of these graph-based reaction discovery schemes, with particular emphasis on more recent methods incorporating graph-based methods with semiempirical and ab initio electronic structure calculations, minimum-energy path refinements, and transition state searches. Using representative examples from homogeneous catalysis and interstellar chemistry, we highlight how these schemes increasingly act as "virtual reaction vessels" for interrogating mechanistic questions. Finally, we highlight where challenges remain, including issues of chemical accuracy and calculation speeds, as well as the inherent challenge of dealing with the vast size of accessible chemical reaction space.

High-Throughput Screening of Bicationic Redox Materials for Chemical Looping Ammonia Synthesis

Ammonia recently has gained increasing attention as a carrier for the efficient and safe usage of hydrogen to further advance the hydrogen economy. However, there is a pressing need to develop new ammonia synthesis techniques to overcome the problem of intense energy consumption associated with the widely used Haber-Bosch process. Chemical looping ammonia synthesis (CLAS) is a promising approach to tackle this problem, but the ideal redox materials to drive these chemical looping processes are yet to be discovered. Here, by mining the well-established MP database, the reaction free energies for CLAS involving 1699 bicationic inorganic redox pairs are screened to comprehensively investigate their potentials as efficient redox materials in four different CLAS schemes. A state-of-the-art machine learning strategy is further deployed to significantly widen the chemical space for discovering the promising redox materials from more than half a million candidates. Most importantly, using the three-step H2 O-CL as an example, a new metric is introduced to determine bicationic redox pairs that are "cooperatively enhanced" compared to their corresponding monocationic counterparts. It is found that bicationic compounds containing a combination of alkali/alkaline-earth metals and transition metal (TM)/post-TM/metalloid elements are compounds that are particularly promising in this respect.

Toward Computing Accurate Free Energies in Heterogeneous Catalysis: a Case Study for Adsorbed Isobutene in H-ZSM-5

Herein, we propose a novel computational protocol that enables calculating free energies with improved accuracy by combining the best available techniques for enthalpy and entropy calculation. While the entropy is described by enhanced sampling molecular dynamics techniques, the energy is calculated using ab initio methods. We apply the method to assess the stability of isobutene adsorption intermediates in the zeolite H-SSZ-13, a prototypical problem that is computationally extremely challenging in terms of calculating enthalpy and entropy. We find that at typical operating conditions for zeolite catalysis (400 °C), the physisorbed π-complex, and not the tertiary carbenium ion as often reported, is the most stable intermediate. This method paves the way for sampling-based techniques to calculate the accurate free energies in a broad range of chemistry-related disciplines, thus presenting a big step forward toward predictive modeling.

Stability of oxidized states of free-standing and ceria-supported PtOx particles

Nanostructured materials based on CeO2 and Pt play a fundamental role in catalyst design. However, their characterization is often challenging due to their structural complexity and the tendency of the materials to change under reaction conditions. In this work, we combine calculations based on the density functional theory, a machine-learning assisted global optimization method (GOFEE), and ab initio thermodynamics to characterize stable oxidation states of ceria-supported PtyOx particles in different environments. The collection of global minima structures for different stoichiometries resulting from the global optimisation effort is used to assess the effect of temperature, oxygen pressure, and support interactions on the phase diagrams, oxidation states, and geometries of the PtyOx particles. We thus identify favoured structural motifs and O:Pt ratios, revealing that oxidized states of free-standing and ceria-supported platinum particles are more stable than reduced ones under a wide range of conditions. These results indicate that studies rationalizing activity of ceria-supported Pt clusters must consider oxidized states, and that previous understanding of such materials obtained only with fully reduced Pt clusters may be incomplete.