How to Gain Atomistic Insights on Reactions at the Water/Solid Interface?

Modeling Single-Atom Catalysis

Electronic structure calculations represent an essential complement of experiments to characterize single-atom catalysts (SACs), consisting of isolated metal atoms stabilized on a support, but also to predict new catalysts. However, simulating SACs with quantum chemistry approaches is not as simple as often assumed. In this work, the essential factors that characterize a reliable simulation of SACs activity are examined. The Perspective focuses on the importance of precise atomistic characterization of the active site, since even small changes in the metal atom's surroundings can result in large changes in reactivity. The dynamical behavior and stability of SACs under working conditions, as well as the importance of adopting appropriate methods to solve the Schrödinger equation for a quantitative evaluation of reaction energies are addressed. The Perspective also focuses on the relevance of the model adopted. For electrocatalysis this must include the effects of the solvent, the presence of electrolytes, the pH, and the external potential. Finally, it is discussed how the similarities between SACs and coordination compounds may result in reaction intermediates that usually are not observed on metal electrodes. When these aspects are not adequately considered, the predictive power of electronic structure calculations is quite limited.

Characterizing Interfaces by Voronoi Tessellation

The chemistry of interfaces differs markedly from that of the bulk. Calculation of interfacial properties depends strongly on the definition of the interface, which can lead to ambiguous results that vary between studies. There is a need for a method that can explicitly define the interfaces and boundaries in molecular systems. Voronoi tessellation offers an attractive solution to this problem through its ability to determine neighbors among specified groups of atoms. Here we discuss three cases where Voronoi tessellation combined with modeling of vibrational sum frequency generation (SFG) spectroscopy yields relevant insights: the breakdown of the air-water interface into clear and intuitive molecular layers, the study of the hydration shell in biological systems, and the acceleration of difficult spectral calculations where intermolecular vibrational couplings dominate. The utility of Voronoi tessellation has broad applications that extend beyond any single type of spectroscopy or system.

CO Electroreduction Mechanism on Single-Atom Zn (101) Surfaces: Pathway to C2 Products

Electrocatalytic reduction of carbon dioxide (CO₂RR) employs electricity to store renewable energy in the form of reduction products. The activity and selectivity of the reaction depend on the inherent properties of electrode materials. Single-atom alloys (SAAs) exhibit high atomic utilization efficiency and unique catalytic activity, making them promising alternatives to precious metal catalysts. In this study, density functional theory (DFT) was employed to predict stability and high catalytic activity of Cu/Zn (101) and Pd/Zn (101) catalysts in the electrochemical environment at the single-atom reaction site. The mechanism of C2 products (glyoxal, acetaldehyde, ethylene, and ethane) produced by electrochemical reduction on the surface was elucidated. The C-C coupling process occurs through the CO dimerization mechanism, and the formation of the *CHOCO intermediate proves beneficial, as it inhibits both HER and CO protonation. Furthermore, the synergistic effect between single atoms and Zn results in a distinct adsorption behavior of intermediates compared to traditional metals, giving SAAs unique selectivity towards the C2 mechanism. At lower voltages, the Zn (101) single-atom alloy demonstrates the most advantageous performance in generating ethane on the surface, while acetaldehyde and ethylene exhibit significant certain potential. These findings establish a theoretical foundation for the design of more efficient and selective carbon dioxide catalysts.

Evaluating Adsorbate-Solvent Interactions: Are Dispersion Corrections Necessary?

Incorporating solvent-adsorbate interactions is paramount in models of aqueous (electro)catalytic reactions. Although a number of techniques exist, they are either highly demanding in computational terms or inaccurate. Microsolvation offers a trade-off between accuracy and computational expenses. Here, we dissect a method to swiftly outline the first solvation shell of species adsorbed on transition-metal surfaces and assess their corresponding solvation energy. Interestingly, dispersion corrections are generally not needed in the model, but caution is to be exercised when water-water and water-adsorbate interactions are of similar magnitude.

An Electrostatically Embedded QM/MM Scheme for Electrified Interfaces

Atomistic modeling of electrified interfaces remains a major issue for detailed insights in electrocatalysis, corrosion, electrodeposition, batteries, and related devices such as pseudocapacitors. In these domains, the use of grand-canonical density functional theory (GC-DFT) in combination with implicit solvation models has become popular. GC-DFT can be conveniently applied not only to metallic surfaces but also to semiconducting oxides and sulfides and is, furthermore, sufficiently robust to achieve a consistent description of reaction pathways. However, the accuracy of implicit solvation models for solvation effects at interfaces is in general unknown. One promising way to overcome the limitations of implicit solvents is going toward hybrid quantum mechanical (QM)/molecular mechanics (MM) models. For capturing the electrochemical potential dependence, the key quantity is the capacitance, i.e., the relation between the surface charge and the electrochemical potential. In order to retrieve the electrochemical potential from a QM/MM hybrid scheme, an electrostatic embedding is required. Furthermore, the charge of the surface and of the solvent regions has to be strictly opposite in order to consistently simulate charge-neutral unit cells in MM and in QM. To achieve such a QM/MM scheme, we present the implementation of electrostatic embedding in the VASP code. This scheme is broadly applicable to any neutral or charged solid/liquid interface. Here, we demonstrate its use in the context of GC-DFT for the hydrogen evolution reaction (HER) over a noble-metal-free electrocatalyst, MoS₂. We investigate the effect of electrostatic embedding compared to the implicit solvent model for three contrasting active sites on MoS₂: (i) the sulfur vacancy defect, which is rather apolar; (ii) a Mo antisite defect, where the active site is a surface bound highly polar OH group; and (iii) a reconstructed edge site, which is generally believed to be responsible for most of the catalytic activity. According to our results, the electrostatic embedding leads to almost indistinguishable results compared to the implicit solvent for the apolar system but has a significant effect on polar sites. This demonstrates the reliability of the hybrid QM/MM, electrostatically embedded solvation model for electrified interfaces.

Thickness-dependent catalytic activity of hydrogen evolution based on single atomic catalyst of Pt above MXene

Hydrogen as the cleanest energy carrier is a promising alternative renewable resource to fossil fuels. There is an ever-increasing interest in exploring efficient and cost-effective approaches of hydrogen production. Recent experiments have shown that single platinum atom immobilized on the metal vacancies of MXenes allows a high-efficient hydrogen evolution reaction (HER). Here usingab initiocalculations, we design a series of substitutional Pt-doped Ti_{n+ 1}C_nT_x(Ti_{n+ 1}C_nT_x-Pt_SA) with different thicknesses and terminations (n= 1, 2 and 3, T_x= O, F and OH), and investigate the quantum-confinement effect on the HER catalytic performance. Surprisingly, we reveal a strong thickness effect of the MXene layer on the HER performance. Among the various surface-terminated derivatives, Ti₂CF₂-Pt_SAand Ti₂CH₂O₂-Pt_SAare found to be the best HER catalysts with the change of Gibbs free energy ΔG_H*∼ 0 eV, complying with the thermoneutral condition. Theab initiomolecular dynamics simulations reveal that Ti₂CF₂-Pt_SAand Ti₂CH₂O₂-Pt_SApossess a good thermodynamic stability. The present work shows that the HER catalytic activity of the MXene is not solely governed by the local environment of the surface such as Pt single atom. We point out the critical role of thickness control and surface decoration of substrate in achieving a high-performance HER catalytical activity.

Differences in solvation thermodynamics of oxygenates at Pt/Al2O3 perimeter versus Pt(111) terrace sites

A prominent role of water in aqueous-phase heterogeneous catalysis is to modify free energies; however, intuition about how is based largely on pure metal surfaces or even homogeneous solutions. Using multiscale modeling with explicit liquid water molecules, we show that the influence of water on the free energies of adsorbates at metal/support interfaces is different than that on pure metal surfaces. We specifically compute free energies of solvation for methanol and its constituents on a Pt/Al₂O₃ catalyst and compare the results to analogous values calculated on a pure Pt catalyst. We find that the more hydrophilic Pt/Al₂O₃ interface leads to smaller (more positive) free energies of solvation due to an increased entropy penalty resulting from the additional work necessary to disrupt the interfacial water structure and accommodate the interfacial species. The results will be of interest in other fields, including adsorption and proteins.

Materials Screening by the Descriptor G max(η): The Free-Energy Span Model in Electrocatalysis

To move from fossil-based energy resources to a society based on renewables, electrode materials free of precious noble metals are required to efficiently catalyze electrochemical processes in fuel cells, batteries, or electrolyzers. Materials screening operating at minimal computational cost is a powerful method to assess the performance of potential electrode compositions based on heuristic concepts. While the thermodynamic overpotential in combination with the volcano concept refers to the most popular descriptor-based analysis in the literature, this notion cannot reproduce experimental trends reasonably well. About two years ago, the concept of G _max(η), based on the idea of the free-energy span model, has been proposed as a universal approach for the screening of electrocatalysts. In contrast to other available descriptor-based methods, G _max(η) factors overpotential and kinetic effects by a dedicated evacuation scheme of adsorption free energies into an analysis of trends. In the present perspective, we discuss the application of G _max(η) to different electrocatalytic processes, including the oxygen evolution and reduction reactions, the nitrogen reduction reaction, and the selectivity problem of the competing oxygen evolution and peroxide formation reactions, and we outline the advantages of this screening approach over previous investigations.

Electrocatalytic Reduction of CO2 to C1 Compounds by Zn-Based Monatomic Alloys: A DFT Calculation

Electrocatalytic reduction of carbon dioxide to produce usable products and fuels such as alkanes, alkenes, and alcohols, is a very promising strategy. Recent experiments have witnessed great advances in precisely controlling the synthesis of single atom alloys (SAAs), which exhibit unique catalytic properties different from alloys and nanoparticles. However, only certain precious metals, such as Pd or Au, can achieve this transformation. Here, the density functional theory (DFT) calculations were performed to show that Zn-based SAAs are promising electrocatalysts for the reduction of CO2 to C1 hydrocarbons. We assume that CO2 reduction in Zn-based SAAs follows a two-step continuous reaction: first Zn reduces CO2 to CO, and then newly generated CO is captured by M and further reduced to C1 products such as methane or methanol. This work screens seven stable alloys from 16 SAAs (M = Fe, Ru, Os, Co, Rh, Ir, Ni, Pd, Pt, Cu, Ag, Au, V, Mo, Ti, Cr). Among them, Pd@Zn (101) and Cu@Zn (101) are promising catalysts for CO2 reduction. The reaction mechanisms of these two SAAs are discussed in detail. Both of them convert CO2 into methane via the same pathway. They are reduced by the pathway: *CO2 → *COOH → *CO + H2O; *CO → *CHO → *CH2O → *CH3O → *O + CH4 → *OH + CH4 → H2O + CH4. However, their potential determination steps are different, i.e., *CO2 → *COOH (ΔG = 0.70 eV) for Cu@Zn (101) and *CO → *CHO (ΔG = 0.72 eV) for Pd@Zn, respectively. This suggests that Zn-based SAAs can reduce CO2 to methane with a small overpotential. The solvation effect is simulated by the implicit solvation model, and it is found that H2O is beneficial to CO2 reduction. These computational results show an effective monatomic material to form hydrocarbons, which can stimulate experimental efforts to explore the use of SAAs to catalyze CO2 electrochemical reduction to hydrocarbons.

Beyond single-crystal surfaces: The GAL21 water/metal force field

Solvent effects are notoriously difficult to describe for metallic nanoparticles (NPs). Here, we introduce GAL21 which is the first pairwise additive force field that is specifically designed to modulate the near chemisorption energy of water as a function of the coordination numbers of the metallic atoms. We find a quadratic dependence to be most suitable for capturing the dependence of the adsorption energy of water on the generalized coordination number (GCN) of the metal atoms. GAL21 has been fitted against DFT adsorption energies for Cu, Ag, Au, Ni, Pd, Pt, and Co on 500 configurations and validated on about 3000 configurations for each metal, constructed on five surfaces with GCNs varying from 2.5 to 11.25. Depending on the metals, the root mean square deviation is found between 0.7 kcal mol−1 (Au) to 1.6 kcal mol−1 (Ni). Using GAL21, as implemented in the open-source code CP2K, we then evaluate the solvation energy of Au55 and Pt55 NPs in water using thermodynamic integration. The solvation free energy is found to be larger for Pt than for Au and systematically larger than 200 kcal mol−1, demonstrating the large impact of solvent on the surface energetics of NPs. Still, given that the amorphous NPs are both, the most stable and the most solvated ones, we do not predict a change in the preferred morphology between the gas-phase and in water. Finally, based on a linear regression on three sizes of NPs (from 38 to 147), the solvation energy for Au and Pt surface atoms is found to be −5.2 and −9.9 kcal mol−1, respectively.