A general-purpose machine learning Pt interatomic potential for an accurate description of bulk, surfaces, and nanoparticles.

A Gaussian approximation machine learning interatomic potential for platinum is presented. It has been trained on density-functional theory (DFT) data computed for bulk, surfaces, and nanostructured platinum, in particular nanoparticles. Across the range of tested properties, which include bulk elasticity, surface energetics, and nanoparticle stability, this potential shows excellent transferability and agreement with DFT, providing state-of-the-art accuracy at a low computational cost. We showcase the possibilities for modeling of Pt systems enabled by this potential with two examples: the pressure-temperature phase diagram of Pt calculated using nested sampling and a study of the spontaneous crystallization of a large Pt nanoparticle based on classical dynamics simulations over several nanoseconds.

Reassignment of magic numbers for icosahedral Au clusters: 310, 564, 928 and 1426.

Icosahedral Au clusters with three and four shells of atoms are found to deviate significantly from the commonly assumed Mackay structures. By introducing additional atoms in the surface shell and creating a vacancy in the center of the cluster, the calculated energy per atom can be lowered significantly, according to several different descriptions of the interatomic interaction. Analogous icosahedral structures with five and six shells of atoms are generated using the same structural motifs and are similarly found to be more stable than Mackay icosahedra. The lowest energy per atom obtained here is for clusters containing 310, 564, 928 and 1426 atoms, as compared with the commonly assumed magic numbers of 309, 561, 923 and 1415. Some of the vertices in the optimized clusters have a hexagonal ring of atoms, rather than a pentagon, with the vertex atom missing. An inner shell atom in some cases moves outwards by more than an Ångström into the surface shell at such a vertex site. This feature, as well as the wide distribution of nearest-neighbor distances in the surface layer, can strongly influence the properties of icosahedral clusters, for example catalytic activity. The structural optimization is initially carried out using the GOUST method with atomic forces estimated with the EMT empirical potential function, but the atomic coordinates are then refined by minimization using electron density functional theory (DFT) or Gaussian approximation potential (GAP). A single energy barrier is found to separate the Mackay icosahedron from a lower energy structure where a string of atoms moves outwards in a concerted manner from the center so as to create a central vacancy while placing an additional atom in the surface shell.

Amorphization mechanism of SrIrO3 electrocatalyst: How oxygen redox initiates ionic diffusion and structural reorganization.

The use of renewable electricity to prepare materials and fuels from abundant molecules offers a tantalizing opportunity to address concerns over energy and materials sustainability. The oxygen evolution reaction (OER) is integral to nearly all material and fuel electrosyntheses. However, very little is known about the structural evolution of the OER electrocatalyst, especially the amorphous layer that forms from the crystalline structure. Here, we investigate the interfacial transformation of the SrIrO3 OER electrocatalyst. The SrIrO3 amorphization is initiated by the lattice oxygen redox, a step that allows Sr2+ to diffuse and O2- to reorganize the SrIrO3 structure. This activation turns SrIrO3 into a highly disordered Ir octahedral network with Ir square-planar motif. The final Sr y IrO x exhibits a greater degree of disorder than IrO x made from other processing methods. Our results demonstrate that the structural reorganization facilitated by coupled ionic diffusions is essential to the disordered structure of the SrIrO3 electrocatalyst.

All-electron ab initio Bethe-Salpeter equation approach to neutral excitations in molecules with numeric atom-centered orbitals.

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as "Thiel's set." Literature reference data based on Gaussian-type orbitals are reproduced to about one millielectron-volt precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted "tier2+aug2" is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

Sr3CrN3: A New Electride with a Partially Filled d-Shell Transition Metal.

Electrides are ionic crystals in which the electrons prefer to occupy free space, serving as anions. Because the electrons prefer to be in the pockets, channels, or layers to the atomic orbitals around the nuclei, it has been challenging to find electrides with partially filled d-shell transition metals, since an unoccupied d-shell provides an energetically favorable location for the electrons to occupy. We recently predicted the existence of electrides with partially filled d-shells using high-throughput computational screening. Here, we provide experimental support using X-ray absorption spectroscopy and X-ray and neutron diffraction to show that Sr3CrN3 is indeed an electride despite its partial d-shell configuration. Our findings indicate that Sr3CrN3 is the first known electride with a partially filled d-shell transition metal, in agreement with theory, which significantly broadens the criteria for the search for new electride materials.

Measurements of Oxygen Electroadsorption Energies and Oxygen Evolution Reaction on RuO2(110): A Discussion of the Sabatier Principle and Its Role in Electrocatalysis.

We report the hydroxide (OHad) and oxide (Oad) experimental electroadsorption free energies, their dependences on pH, and their correlations to the oxygen evolution reaction (OER) electrocatalysis on RuO2(110) surface. The Sabatier principle predicts that catalyst is most active when the intermediate stabilization is moderate, not too strong such that the bound intermediate disrupts the subsequent catalytic cycle, nor too weak such that the surface is ineffective. For decades, researchers have used this concept to rationalize the activity trend of many OER electrocatalysts including RuO2, which is among the state-of-the-art OER catalysts. In this article, we report an experimental assessment of the Sabatier principle by comparing the oxygen electroadsorption energy to the OER electrocatalysis for the first time on RuO2. We find that the OHad and Oad electroadsorption energies on RuO2(110) depend on pH and obey the scaling relation. However, we did not observe a direct correlation between the OHad and Oad electroadsorption energies and the OER activity in the comparative analysis that includes both RuO2(110) and IrO2(110). Our result raises a question of whether the Sabatier principle can describe highly active electrocatalysts, where the kinetic aspects may influence the electrocatalysis more strongly than the electroadsorption energy, which captures only the thermodynamics of the intermediates and not yet kinetics.

Influence of Surface Adsorption on the Oxygen Evolution Reaction on IrO2(110).

A catalyst functions by stabilizing reaction intermediates, usually through surface adsorption. In the oxygen evolution reaction (OER), surface oxygen adsorption plays an indispensable role in the electrocatalysis. The relationship between the adsorption energetics and OER kinetics, however, has not yet been experimentally measured. Herein we report an experimental relationship between the adsorption of surface oxygen and the kinetics of the OER on IrO2(110) epitaxially grown on a TiO2(110) single crystal. The high quality of the IrO2 film grown using molecular-beam epitaxy affords the ability to extract the surface oxygen adsorption and its impact on the OER. By examining a series of electrolytes, we find that the adsorption energy changes linearly with pH, which we attribute to the electrified interfacial water. We support this hypothesis by showing that an electrolyte salt modification can lead to an adsorption energy shift. The dependence of the adsorption energy on pH has implications for the OER kinetics, but it is not the only factor; the dependence of the OER electrocatalysis on pH stipulates two OER mechanisms, one operating in acidic solution and another operating in alkaline solution. Our work points to the subtle adsorption-kinetics relationship in the OER and highlights the importance of the interfacial electrified interaction in electrocatalyst design.