Toward Fast and Reliable Potential Energy Surfaces for Metallic Pt Clusters by Hierarchical Delta Neural Networks.

Data-driven machine learning force fields (MLFs) are more and more popular in atomistic simulations and exploit machine learning methods to predict energies and forces for unknown structures based on the knowledge learned from an existing reference database. The latter usually comes from density functional theory calculations. One main drawback of MLFs is that physical laws are not incorporated in the machine learning models, and instead, MLFs are designed to be very flexible to simulate complex quantum chemistry potential energy surface (PES). In general, MLFs have poor transferability, and hence, a very large trainset is required to span all the target feature space to get a reliable MLF. This procedure becomes more troublesome when the PES is complicated, with a large number of degrees of freedom, in which building a large database is inevitable and very expensive, especially when accurate but costly exchange-correlation functionals have to be used. In this manuscript, we exploit a high-dimensional neural network potential (HDNNP) on Pt clusters of sizes from 6 to 20 as one example. Our standard level of energy calculation is DFT GGA (PBE) using a plane wave basis set. We introduce an approximate but fast level with the PBE functional and a minimal atomic orbital basis set, and then, a more accurate but expensive level, using a hybrid functional or nonlocal vdW functional and a plane wave basis set, is reliably predicted by learning the difference with HDNNP. The results show that such a differential approach (named ΔHDNNP) can deliver very accurate predictions (error <10 meV/atom) in reference to converged basis set energies as well as more accurate but expensive xc functionals. The overall speedup can be as large as 900 for a 20 atom Pt cluster. More importantly, ΔHDNNP shows much better transferability due to the intrinsic smoothness of the delta potential energy surface, and accordingly, one can use much smaller trainset data to obtain better accuracy than the conventional HDNNP. A multilayer ΔHDNNP is thus proposed to obtain very accurate predictions versus expensive nonlocal vdW functional calculations in which the required trainset is further reduced. The approach can be easily generalized to any other machine learning methods and opens a path to study the structure and dynamics of Pt clusters and nanoparticles

Tuning the Hydrogenation Selectivity of an Unsaturated Aldehyde via Single-Atom Alloy Catalysts

Selective hydrogenation of α,β-unsaturated aldehydes to produce unsaturated alcohols remains a challenge in catalysis. Here, we explore, on the basis of first-principles simulations, single-atom alloy (SAA) catalysts on copper as a class of catalytic materials to enhance the selectivity for C═O bond hydrogenation in unsaturated aldehydes by controlling the binding strength of the C═C and C═O bonds. We show that on SAA of early transition metals such as Ti, Zr, and Hf, the C═O binding mode of acrolein is favored but the strong binding renders subsequent hydrogenation and desorption impossible. On SAA of late-transition metals, on the other hand, the C═C binding mode is favored and C═C bond hydrogenation follows, resulting in the production of undesired saturated aldehydes. Mid-transition metals (Cr and Mn) in Cu(111) appear as the optimal systems, since they favor acrolein adsorption via the C═O bond but with a moderate binding strength, compatible with catalysis. Additionally, acrolein migration from the C═O to the C═C binding mode, which would open the low energy path for C═C bond hydrogenation, is prevented by a large barrier for this process. SAA of Cr in Cu appears as an optimal candidate, and kinetic simulations show that the selectivity for propenol formation is controlled by preventing the acrolein migration from the more stable C═O to the less stable C═C binding mode and subsequent H-migration and by the formation of the O-H bond from the monohydrogenated intermediate. Dilute alloy catalysts therefore enable tuning the binding strength of intermediates and transition states, opening control of catalytic activity and selectivity

Rh single atoms on TiO2 dynamically respond to reaction conditions by adapting their site.

Single-atom catalysts are widely investigated heterogeneous catalysts; however, the identification of the local environment of single atoms under experimental conditions, as well as operando characterization of their structural changes during catalytic reactions are still challenging. Here, the preferred local coordination of Rh single atoms is investigated on TiO2 during calcination in O2, reduction in H2, CO adsorption, and reverse water gas shift (RWGS) reaction conditions. Theoretical and experimental studies clearly demonstrate that Rh single atoms adapt their local coordination and reactivity in response to various redox conditions. Single-atom catalysts hence do not have static local coordinations, but can switch from inactive to active structure under reaction conditions, hence explaining some conflicting literature accounts. The combination of approaches also elucidates the structure of the catalytic active site during reverse water gas shift. This insight on the real nature of the active site is key for the design of high-performance catalysts

H and CO Co-Induced Roughening of Cu Surface in CO2 Electroreduction Conditions

The dynamic restructuring of Cu has been observed under electrochemical conditions, and it has been hypothesized to underlie the unique reactivity of Cu toward CO2 electroreduction. Roughening is one of the key surface phenomena for Cu activation, whereby numerous atomic vacancies and adatoms form. However, the atomic structure of such surface motifs in the presence of relevant adsorbates has remained elusive. Here, we explore the chemical space of Cu surface restructuring under coverage of CO and H in realistic electroreduction conditions, by combining grand canonical DFT and global optimization techniques, from which we construct a potential-dependent grand canonical ensemble representation. The regime of intermediate and mixed CO and H coverage─where structures exhibit some elevated surface Cu─is thermodynamically unfavorable yet kinetically inevitable. Therefore, we develop a quasi-kinetic Monte Carlo simulation to track the system's evolution during a simulated cathodic scan. We reveal the evolution path of the system across coverage space and identify the accessible metastable structures formed along the way. Chemical bonding analysis is performed on the metastable structures with elevated Cu*CO species to understand their formation mechanism. By molecular dynamics simulations and free energy calculations, the surface chemistry of the Cu*CO species is explored, and we identify plausible mechanisms via which the Cu*CO species may diffuse or dimerize. This work provides rich atomistic insights into the phenomenon of surface roughening and the structure of involved species. It also features generalizable methods to explore the chemical space of restructuring surfaces with mixed adsorbates and their nonequilibrium evolution