In silico characterization of nanoparticles.

Nanoparticles (NPs) make for intriguing heterogeneous catalysts due to their large active surface area and excellent and often size-dependent catalytic properties that emerge from a multitude of chemically different surface reaction sites. NP catalysts are, in principle, also highly tunable: even small changes to the NP size or surface facet composition, doping with heteroatoms, or changes of the supporting material can significantly alter their physicochemical properties. Because synthesis of size- and shape-controlled NP catalysts is challenging, the ability to computationally predict the most favorable NP structures for a catalytic reaction of interest is an in-demand skill that can help accelerate and streamline the material optimization process. Fundamentally, simulations of NP model systems present unique challenges to computational scientists. Not only must considerable methodological hurdles be overcome in performing calculations with hundreds to thousands of atoms while retaining appropriate accuracy to be able to probe the desired properties. Also, the data generated by simulations of NPs are typically more complex than data from simulations of, for example, single crystal surface models, and therefore often require different data analysis strategies. To this end, the present work aims to review analytical methods and data analysis strategies that have proven useful in extracting thermodynamic trends from NP simulations.

Atomistic Studies on Water-Induced Lithium Corrosion.

It is well known that lithium reacts violently with water under the release of molecular hydrogen and the formation of lithium hydroxide. In this work, the initial mechanisms for the surface reactions of metallic lithium with water from the gas phase were investigated by means of periodic density functional theory calculations. For this purpose, adsorption/absorption structures and diffusion and dissociation processes of hydrogen, OH, and H2 O on low-index metallic lithium surfaces were investigated. Through thermodynamic and kinetic considerations, negatively charged centers on the surface were identified as the origin of hydrogen formation. The strikingly low reaction barriers for the reaction at these centers implied a self-supporting effect of hydrogen evolution and the associated lithium degradation.

Assessment of the Accuracy of Density Functionals for Calculating Oxygen Reduction Reaction on Nitrogen-Doped Graphene.

Experimental studies of the oxygen reduction reaction (ORR) at nitrogen-doped graphene electrodes have reported a remarkably low overpotential, on the order of 0.5 V, similar to Pt-based electrodes. Theoretical calculations using density functional theory have lent support to this claim. However, other measurements have indicated that transition metal impurities are actually responsible for the ORR activity, thereby raising questions about the reliability of both the experiments and the calculations. To assess the accuracy of the theoretical calculations, various generalized gradient approximation (GGA), meta-GGA, and hybrid functionals are employed here and calibrated against high-level wave-function-based coupled-cluster calculations (CCSD(T)) of the overpotential as well as self-interaction corrected density functional calculations and published quantum Monte Carlo calculations of O adatom binding to graphene. The PBE0 and HSE06 hybrid functionals are found to give more accurate results than the GGA and meta-GGA functionals, as would be expected, and for a low dopant concentration, 3.1%, the overpotential is calculated to be 1.0 V. The GGA and meta-GGA functionals give a lower estimate by as much as 0.4 V. When the dopant concentration is doubled, the overpotential calculated with hybrid functionals decreases, while it increases in GGA functional calculations. The opposite trends result from different potential-determining steps, the *OOH species being of central importance in the hybrid functional calculations, while the reduction of *O determines the overpotential obtained in GGA and meta-GGA calculations. The results presented here are mainly based on calculations of periodic representations of the system, but a comparison is also made with molecular flake models that are found to give erratic results due to finite size effects and geometric distortions during energy minimization. The presence of the electrolyte has not been taken into account explicitly in the calculations presented here but is estimated to be important for definitive calculations of the overpotential.

Simulations of the Oxidation and Degradation of Platinum Electrocatalysts.

Improved understanding of the fundamental processes leading to degradation of platinum nanoparticle electrocatalysts is essential to the continued advancement of their catalytic activity and stability. To this end, the oxidation of platinum nanoparticles is simulated using a ReaxFF reactive force field within a grand-canonical Monte Carlo scheme. 2-4 nm cuboctahedral particles serve as model systems, for which electrochemical potential-dependent phase diagrams are constructed from the thermodynamically most stable oxide structures, including solvation and thermochemical contributions. Calculations in this study suggest that surface oxide structures should become thermodynamically stable at voltages around 0.80-0.85 V versus standard hydrogen electrode, which corresponds to typical fuel cell operating conditions. The potential presence of a surface oxide during catalysis is usually not accounted for in theoretical studies of Pt electrocatalysts. Beyond 1.1 V, fragmentation of the catalyst particles into [Pt6 O8 ]4- clusters is observed. Density functional theory calculations confirm that [Pt6 O8 ]4- is indeed stable and hydrophilic. These results suggest that the formation of [Pt6 O8 ]4- may play an important role in platinum catalyst degradation as well as the electromotoric transport of Pt2+/4+ ions in fuel cells.

Microscopic Properties of Na and Li-A First Principle Study of Metal Battery Anode Materials.

Using density functional theory, we studied the bulk and surface properties of Li and Na electrodes on an atomistic level. To get a better understanding of the initial stages of surface growth phenomena (and thus dendrite formation), various self-diffusion mechanisms were studied. For this purpose, dedicated diffusion pathways on the surfaces of Na and Li were investigated within the terrace-step-kink (TSK) model utilizing nudged elastic band calculations. We were able to prove that the mere investigation of terrace self-diffusion on the respective low-index surfaces does not provide a possible descriptor for dendritic growth. Finally, we provide an initial view of the surface growth behavior of both alkali metals as well as provide a basis for experimental investigations and theoretical long-scale kinetic Monte Carlo simulations.

First principles studies of self-diffusion processes on metallic lithium surfaces.

Due to the theoretical high specific capacity (3860 mAh/g) and the low standard electrode potential (-3.040 V vs. standard hydrogen electrode), rechargeable lithium metal batteries are considered as excellent energy storage systems. Unfortunately, security concerns related to dendrite formation during charge/discharge cycles still hinder the commercial use of Li metal-based batteries. Using density functional theory, we have studied the bulk and surface properties of metallic lithium at an atomistic level. In this process, bcc Li(100) proved to be the most stable metallic lithium surface. Subsequently, possible self-diffusion mechanisms on perfect and imperfect Li(100) surfaces were examined. For this purpose, nudged elastic band calculations were performed to characterize the respective diffusion processes and to determine the relevant pre-exponential factors and activation barriers. On the basis of the acquired data, it became possible to derive activation temperatures and reaction rates for the respective processes, which are useful for experimental verification as well as for the implementation in long-scale kinetic Monte Carlo simulations.

Molybdenum Doping Augments Platinum-Copper Oxygen Reduction Electrocatalyst.

Improving the efficiency of Pt-based oxygen reduction reaction (ORR) catalysts while also reducing costs remains an important challenge in energy research. To this end, we synthesized highly stable and active carbon-supported Mo-doped PtCu (Mo-PtCu/C) nanoparticles (NPs) from readily available precursors in a facile one-pot reaction. Mo-PtCu/C displays two-to-fourfold-higher ORR half-cell kinetics than reference PtCu/C and Pt/C materials, a trend that was confirmed in proof-of-concept experiments by using a H2 /O2 microlaminar fuel cell. This Mo-induced activity increase mirrors observations for Mo-PtNi/C NPs and possibly suggests an emerging trend. Electrochemical-accelerated stability tests revealed that dealloying was greatly reduced in Mo-PtCu/C in contrast to the binary alloys PtCu/C and PtMo/C. Supporting DFT studies suggested that the exceptional stability of Mo-PtCu could be attributed to oxidative resistance of the Mo-doped atoms. Furthermore, our calculations revealed that oxygen could induce segregation of Mo to the catalytic surface, at which it effected beneficial changes to the surface oxygen adsorption energetics in the context of the Sabatier principle.

Growth of Stable Surface Oxides on Pt(111) at Near-Ambient Pressures.

Detailed knowledge of the structure and degree of oxidation of platinum surfaces under operando conditions is essential for understanding catalytic performance. However, experimental investigations of platinum surface oxides have been hampered by technical limitations, preventing in situ investigations at relevant pressures. As a result, the time-dependent evolution of oxide formation has only received superficial treatment. In addition, the amorphous structures of many surface oxides have hindered realistic theoretical studies. Using near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) we show that a time scale of hours (t≥4 h) is required for the formation of platinum surface oxides. These experimental observations are consistent with ReaxFF grand canonical Monte Carlo (ReaxFF-GCMC) calculations, predicting the structures and coverages of stable, amorphous surface oxides at temperatures between 430-680 K and an O2 partial pressure of 1 mbar.

Surface Buckling and Subsurface Oxygen: Atomistic Insights into the Surface Oxidation of Pt(111).

Platinum is a catalyst of choice in scientific investigations and technological applications, which are both often carried out in the presence of oxygen. Thus, a fundamental understanding of platinum's (electro)catalytic behavior requires a detailed knowledge of the structure and degree of oxidation of platinum surfaces in operando. ReaxFF reactive force field calculations of the surface energies for structures with up to one monolayer of oxygen on Pt(111) reveal four stable surface phases characterized by pure adsorbate, high- and low-coverage buckled, and subsurface-oxygen structures, respectively. These structures and temperature programmed desorption (TPD) spectra simulated from them compare favorably with and complement published scanning tunneling microscopy (STM) and TPD experiments. The surface buckling and subsurface oxygen observed here influence the surface oxidation process, and are expected to impact the (electro)catalytic properties of partially oxidized Pt(111) surfaces.

Development of a ReaxFF potential for Pt-O systems describing the energetics and dynamics of Pt-oxide formation.

ReaxFF force field parameters describing Pt-Pt and Pt-O interactions have been developed and tested. The Pt-Pt parameters are shown to accurately account for the chemical nature, atomic structures and other materials properties of bulk platinum phases, low and high-index platinum surfaces and nanoclusters. The Pt-O parameters reliably describe bulk platinum oxides, as well as oxygen adsorption and oxide formation on Pt(111) terraces and the {111} and {100} steps connecting them. Good agreement between the force field and both density functional theory (DFT) calculations and experimental observations is demonstrated in the relative surface free energies of high symmetry Pt-O surface phases as a function of the oxygen chemical potential, making ReaxFF an ideal tool for more detailed investigations of more complex Pt-O surface structures. Validation for its application to studies of the kinetics and dynamics of surface oxide formation in the context of either molecular dynamics (MD) or Monte Carlo simulations are provided in part by a two-part investigation of oxygen diffusion on Pt(111), in which nudged elastic band (NEB) calculations and MD simulations are used to characterize diffusion processes and to determine the relevant diffusion coefficients and barriers. Finally, the power of the ReaxFF reactive force field approach in addressing surface structures well beyond the reach of routine DFT calculations is exhibited in a brief proof-of-concept study of oxygen adsorbate displacement within ordered overlayers.

Nanoscale-faceting of metal surfaces induced by adsorbates.

Using density functional theory and thermodynamic considerations, adsorbate-induced faceting of high-index metal surfaces such as Ir(210) and Re(112 1) has been studied. Focusing on these two systems we first discuss the adsorption behaviour of oxygen and nitrogen on the various surfaces relevant for the faceting, and afterwards use these energies to evaluate the stability of substrates and facets in the presence of oxygen and nitrogen. The faceting phase diagrams of Ir(210) and Re(112 1) show that both adsorbates enhance the anisotropy in surface free energy, finally causing nanofacets to become the thermodynamically favourable surface structure. We also generated analogous electrochemical phase diagrams for both surfaces in contact with an oxygen- or nitrogen-containing electrolyte and found that the same nanofacets should also become stable at positive electrode potentials. Thus, our calculations not only reproduce the experimentally observed surface faceting under UHV conditions, but also predict facet formation under electrochemical conditions.

Two-dimensional assembly of magnetic binuclear complexes: a scanning tunneling microscopy study.

Mono- and binuclear metal-organic compounds bearing long alkyl chains were synthesized and studied at the liquid/graphite interface using scanning tunneling microscopy. Two different lamellar surface patterns as well as a star like structure were obtained driven by van der Waals interactions of the alkyl chains and weak hydrogen bonds of the phenoxy moieties. In the case of the star like structure solvent molecules (1,2,4-trichlorobenzene) are supposed to play an important role for the stabilization of the created pattern. Magnetic investigation of the bulk material by a superconducting quantum interference device magnetometer revealed magnetic moments up to 1.7 mu(B) (NiCo) and most likely antiferromagnetic coupling between the two metals within a single complex. The presented two-dimensional crystallization of the binuclear complexes may provide an easy access to new designable materials in molecular electronics.