Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory (DFT) calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Evaluation of Sb-Nd and Te-Nd phases within the U-Zr fuel matrix and their interactions with HT9 alloy

Antimony (Sb) and tellurium (Te) were investigated as potential additives for U-10Zr (wt.%) metallic fuel to limit the fuel-cladding chemical interaction (FCCI) with HT-9 alloy. Neodymium (Nd) was utilized to simulate the formation of lanthanide-based solid fission products which are known to play a detrimental role in FCCI. Fuel alloys of U-Zr-Sb-Nd and U-Zr-Te-Nd were evaluated in their annealed condition and compared against their as-cast conditions. Isothermal diffusion couple experiments were performed between U-Zr-Nd, U-Zr-Sb-Nd, and U-Zr-Te-Nd against the cladding alloy HT9 to evaluate the effectiveness of the additives to stabilize Nd within the fuel alloys, as well as investigate the interaction regions that form between the different fuel alloys and HT9. Further, SbNd and Sb 3 Nd 4 , and TeNd are found to be the primary neodymium-based phases formed in the U-Zr-Sb-Nd and U-Zr-Te-Nd alloys, respectively. The zirconium-based phase, Zr 2 Sb, is also found to form within the former alloy. All phases were found to remain stable through the diffusion experiments and exhibited no interaction with HT9 constituent elements. Preferential interaction between Nd with additivities Te and Sb compared to constituting elements in HT9 was further verified based on density functional theory (DFT) calculated enthalpy of mixing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhancing perovskite electrocatalysis through synergistic functionalization of B-site cation for efficient water splitting

The family of perovskite oxides is a promising class of catalysts for diverse energy conversion processes including water splitting. In this work, a facile two-step manipulation (in-situ exsolution and post-sulfurization) strategy was proposed and applied to LaCo 0.2 Fe 0.8 O 3 (LCF) perovskite parent, through which, the electronic state, spatial immersion and intrinsic activity of B-site cobalt (Co) were stepwise tuned at nanoscale proximity accordingly (i.e., lattice Co ions segregated Co 0 → embedded CoS 2 ). Impressively, the as-prepared catalyst (S-LCF) obtains an emergent oxygen deficient microstructure seamlessly pinned with uniformly distributed CoS 2 nanoparticles (NPs), which demonstrates enhanced performance toward both oxygen evolution reaction (OER) and hydrogen evolution reaction (HER), and shows good stability in overall water splitting. The density functional theory (DFT) calculations illustrate the optimized metal-oxygen covalency and hydrogen adsorption Gibbs free energy (Δ G H *) on S-LCF, which further buttresses the prominence of our B-site cation engineering tactics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phonon Spectroscopy in Antimony and Tellurium Oxides

α-Sb 2 O 3 (senarmontite), β-Sb 2 O 3 (valentinite), and α-TeO 2 (paratellurite) are compounds with pronounced stereochemically active Sb and Te lone pairs. The vibrational and lattice properties of each have been previously studied but often lead to incomplete or unreliable results due to modes being inactive in infrared or Raman spectroscopy. Here in this paper, we present a study of the relationship between bonding and lattice dynamics of these compounds. Mössbauer spectroscopy is used to study the structure of Sb in α-Sb 2 O 3 and β-Sb 2 O 3 , whereas the vibrational modes of Sb and Te for each oxide are investigated using nuclear inelastic scattering, and further information on O vibrational modes is obtained using inelastic neutron scattering. Additionally, vibrational frequencies obtained by density functional theory (DFT) calculations are compared with experimental results in order to assess the validity of the utilized functional. Good agreement was found between DFT-calculated and experimental density of phonon states with a 7% scaling factor. The Sb–O–Sb wagging mode of α-Sb 2 O 3 whose frequency was not clear in most previous studies is experimentally observed for the first time at ~340 cm -1 . Softer lattice vibrational modes occur in orthorhombic β-Sb 2 O 3 compared to cubic α-Sb 2 O 3 , indicating that the antimony bonds are weakened upon transforming from the molecular α phase to the layer-chained β structure. The resulting vibrational entropy increase of 0.45 ± 0.1 kB/Sb 2 O 3 at 880 K accounts for about half of the α–β transition entropy. The comparison of experimental and theoretical approaches presented here provides a detailed picture of the lattice dynamics in these oxides beyond the zone center and shows that the accuracy of DFT is sufficient for future calculations of similar material structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phonon spectrum in the spin-Peierls phase of CuGeO 3

CuGeO 3 has long been studied as a prototypical example of the spin-Peierls transition in a 𝑆 = 1/2 Heisenberg chain. Despite intensive investigation of this quasi-one-dimensional material, systematic measurements and calculations of the phonon excitations in the dimerized phase have not to date been possible, leaving certain aspects of the spin-Peierls phenomenon unresolved. We perform state-of-the-art density functional theory (DFT) calculations to compute the electronic structure and phonon dynamics in the low-temperature dimerized phase. We also perform high-resolution neutron spectroscopy to measure the full phonon spectrum over multiple Brillouin zones. We find excellent agreement between our numerical and experimental results that extend to all measurement temperatures. Notable features of our phonon spectra include a number of steeply dispersive modes, nonmonotonic dispersion features, and specific phonon anticrossings, which we relate to the mode eigenvectors. By calculating the magnetic interactions within DFT and studying the effects of different phonon modes on the superexchange paths, we discuss the possibility of observing spin-phonon hybridization effects in experiments performed both in and out of equilibrium.

density functional theory↗

Identifying the Molecular Edge Termination of Exfoliated Hexagonal Boron Nitride Nanosheets with Solid-State NMR Spectroscopy and Plane-Wave DFT Calculations

Hexagonal boron nitride nanosheets (h-BNNS), the isoelectronic analog to graphene, have received interest over the past decade due to their high thermal oxidative resistance, high bandgap, catalytic activity, and low cost. The functional groups that terminate boron and nitrogen zigzag and/or armchair edges directly affect their chemical, physical, and electronic properties. However, an understanding of the molecular edge termination present in h-BNNS is lacking. Here, high-resolution magic-angle spinning (MAS) solid-state NMR (SSNMR) spectroscopy, and plane-wave density-functional theory (DFT) calculations are used to determine the molecular edge termination in exfoliated h-BNNS. 1 H → 11 B cross-polarization MAS (CPMAS) SSNMR spectra of h-BNNS revealed multiple hydroxyl/oxygen coordinated boron edge sites that were not detectable in direct excitation experiments. A dynamic nuclear polarization (DNP)-enhanced 1 H → 15 N CPMAS spectrum of h-BNNS displayed four distinct 15 N resonances while a 2D 1 H{ 14 N} dipolar-HMQC spectrum acquired with fast MAS revealed three distinct 14 N environments. Plane-wave DFT calculations were used to construct model edge structures and predict the corresponding 11 B, 14 N and 15 N SSNMR spectra. Comparison of the experimental and predicted SSNMR spectra confirms that zigzag and armchair edges with both amine and boron hydroxide/oxide termination are present. The detailed characterization of h-BNNS molecular edge termination will prove useful for many material science applications. The techniques outlined here should also be applicable to understand the molecular edge terminations in other 2D materials

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DFT calculations and theory do not support enantiospecificity in NMR J-coupling constants

Distinguishing the enantiomers of small organic molecules is an industrially relevant problem with important implications for the health of the population. In a recent publication, Bouchard and co-workers have suggested that large differences in indirect spin-spin (J) coupling constants between enantiomers are possible. A close inspection of their work revealed significant flaws in their density functional theory (DFT) calculations and that the reported effects disappear with appropriate care. We thus conclude that enantiospecificity in spin-spin coupling constants has not be demonstrated either experimentally or theoretically.

Perras, Frédéric A. [Ames Laboratory (AMES), Ames,↗

CH 4 steam reforming on Pt + Pd/Al 2 O 3 monolith: impact of Mn 0.5 Fe 2.5 O 4 spinel addition

The performance of Pt + Pd/Al 2 O 3 monolith catalysts was analyzed under a feed (H 2 + CO + CH 4 + NO + O 2 + H 2 O + CO 2 ) simulating the exhaust gas of stoichiometric natural gas vehicles (NGVs). The flow reactor results show enhanced low temperature (<400 °C) CH 4 oxidation activity under cyclic lean/rich feed modulation and with Mn 0.5 Fe 2.5 O 4 (MFO) spinel addition to the platinum group metal (PGM) catalyst. In the absence of the MFO spinel, the PGM catalyst can achieve full CH 4 conversion, but the presence of spinel in the catalyst limits the high temperature (>500 °C) CH 4 conversion. The depletion of gas phase O 2 , as indicated by the detection of products CO and H 2 in the reactor effluent at higher temperature and CH 4 conversion, portends the emergence of steam reforming of CH 4 (SRM). To investigate the cause of the detrimental impact of spinel under the high temperature, net rich regime, we show through a combination of SRM experiments and post-reaction catalyst characterization that SRM is negatively impacted by the accumulation of Mn and Fe spinel species on the PGM catalyst. This is due to the migration of the base metal species and potential encapsulation and blockage of PGM active sites. The latter is supported by density functional theory (DFT) calculations that indicate favorable metal oxide decoration of active Pt step sites. Lastly, a zoned catalyst design is developed that mitigates the detrimental metals migration while exploiting the methane conversion enhancement afforded by spinel oxides under oxygen rich conditions.

03 NATURAL GAS↗

Hydrogen trapping energetics at BCC iron-helium interfaces

Density functional theory (DFT) calculations have been performed to assess the trapping and segregation strength of hydrogen (H) to helium (He) interfaces in BCC (Fe) as a function of He density and surface orientation. The He density ranges from 0 to 2 He/V which corresponds to equilibrium bubbles that are expected to form in the structural components in a fusion reactor environment. The DFT calculations consist of a slab of BCC metal and an initial lattice of FCC He. Here, the report on the binding energy of H to these interfaces as well as the migration barriers into and away from the surface, which together provide information on the de-trapping energy for H at He bubbles. It was found that the binding energy of H to He-Fe interfaces decreases with increasing He density. Migration barriers are also modified due to the He density, which are sensitive to the surface orientation as well as the diffusion pathway. These results provide valuable first-principles energetics that are necessary for higher-order, mesoscale models, and provide insight into the extent to which He bubbles may trap tritium in fusion structural materials.

42 ENGINEERING↗

Nd, SbNd and Sb 3 Nd 4 and their interactions with the cladding alloy HT9

Lanthanide fission products, such as neodymium, formed during the irradiation of metallic fuels are known to cause deleterious effects from chemical interactions occurring at the fuel-cladding interface; a phenomenon known as fuel-cladding chemical interaction (FCCI). The use of fuel-based additives that bind with the lanthanide elements within the fuel meat, alleviating their interactions at the fuel-cladding interface, is one potential method proposed to mitigate the FCCI phenomenon and extend the burnup potential of such metallic fuel systems. In this study, antimony (Sb) is evaluated as one such additive, and neodymium (Nd) is used to represent the lanthanides. A Sb-Nd alloy is fabricated, which consists of two intermetallic phases, SbNd and Sb 3 Nd 4 . Isothermal diffusion couple experiments are carried out at 675 °C for 24 h between Nd and the Fe-12Cr based HT9 cladding alloy. The results are compared against similar diffusion couple experiments carried out between the Sb-Nd alloy and also HT9. Inter-diffusion between Nd/HT9 diffusion couple is characterized and the phases Fe 17 Nd 5 and Fe 17 Nd 2 were found to form whereas no interactions were observed in the Sb-Nd/HT9 diffusion couple. The lack of compound forming tendencies between Fe, the primary alloying constituent of HT9, and SbNd is elucidated through density functional theory (DFT) calculations on enthalpy of mixing, aligning well with the experimental observations. Lastly, the strong binding of Nd with Sb appears to be favorable, alleviating Nd interactions with HT9 constituent elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Composition-Dependent Reconstructions of Titanium Oxide Clusters and Cu(111) Support via Cluster-Adatom Interactions

Mass-selected cluster deposition and scanning tunneling microscopy (STM) were used to investigate the surface structure and thermal stability of Ti 3 O y (y = 5, 6) clusters on a Cu(111) surface. STM imaging shows that the substoichiometric Ti 3 O 5 clusters aggregate into large dendritic assemblies at step edges, while the stoichiometric Ti 3 O 6 clusters are highly dispersed as single or small groupings of clusters. Cluster deposition also results in extensive reconstruction of the Cu(111) surface. The Ti 3 O 5 /Cu(111) surface exhibits one Cu-layer deep pits that are ringed by cluster assemblies that nucleate and grow at the perimeter step edges. By comparison, the Ti 3 O 6 /Cu(111) surface consists of irregular step edges and small one Cu-layer islands that are decorated with clusters. The restructuring of the Cu(111) surface is attributed to strong Cu-cluster interactions that release Cu adatoms from the step edges. This mechanism is supported by density functional theory (DFT) calculations, which predict that the Cu adatom attachment to the clusters is energetically favorable. Annealing the surfaces at elevated temperatures (400-700 K) results in higher densities of irregular steps, pits and isolated Cu islands as well as increased cluster aggregation. In conclusion, this work highlights adatom-cluster interactions as another manifestation of strong cluster-support interactions that can significantly alter the nature of the oxide-metal interface.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

S K-edge XAS of Cu Ⅱ , Cu Ⅰ , and Zn Ⅱ oxidized Dithiolene complexes: Covalent contributions to structure and the Jahn-Teller effect

We report reduced dithiolene ligands are bound to high valent Mo centers in the active site of the oxotransferase family of enzymes. Related model complexes have been studied with great insight by Prof. Holm and his colleagues. This study focuses on the other limit of dithiolene chemistry: an investigation of the 2-electron oxidized dithiolene bound to low-valent late transition metal (TM) ions (Zn II , Cu I , and Cu II ). The bonding descriptions of the oxidized dithiolene [N,N-dimethyl piperazine 2,3-dithione (Me 2 Dt 0 )] complexes are probed using S K-edge X-ray absorption spectroscopy (XAS) and the results are correlated to density functional theory (DFT) calculations. These experimentally supported calculations are then extended to explain the different geometric structures of the three complexes. The Zn II (Me 2 Dt 0 ) 2 complex has only ligand-ligand repulsion so it is stabilized at the D 2d symmetry limit. The Cu I (Me 2 Dt 0 ) 2 complex has additional weak backbonding thus distorts somewhat from D 2d toward D 2h symmetry. The Cu II (Me 2 Dt 0 ) 2 complex has a strong σ donor bond that leads to both a large Jahn-Teller stabilization to D 2h and an additional covalent contribution to the geometry. The combined strong stabilization results in the square planar, D 2h structure. This study quantifies the competition between the ligand-ligand repulsion and the change in electronic structures in determining the final geometric structures of the oxidized dithiolene complexes, and provides quantitative insights into the Jahn-Teller stabilization energy and its origin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational design of passivants for CdTe grain boundaries

CdTe is the second most-widely deployed photovoltaic (PV) material due to its high efficiency and low manufacturing costs. Currently, polycrystalline CdTe (poly-CdTe) has a record efficiency of similar to 22%, which is still well below the theoretical limit (similar to 30%). Polycrystalline CdTe films that have incorporated Cl or Se show higher efficiency, possibly due to the segregation of these ions to the grain boundaries (GBs), where they may passivate the dangling bonds. However, the efficiency enhancement mechanisms of passivants in CdTe GBs, and the feasibility of employing alternative passivants, have not been well explored. Here, we present a systematic computational study of CdTe GBs with potential passivants, namely S, P, As, Se, and Sb on Te sites, and Na, Mg, Al, Sc, Cu, and Zn on Cd sites. Density functional theory (DFT) calculations were performed on GB dislocation core structures derived from scanning transmission electron microscopy (STEM) images of a model GB (bicrystal). We computed the segregation thermodynamics, electronic density of states, and charge variations near doped CdTe GBs. We find that segregation of impurities to GBs is thermodynamically favorable. For a Te-terminated core, Se, S, and P on Te sites effectively reduce midgap states. For both Cd-and Te-terminated dislocation cores, Sc and Al reduce midgap states when substituted for Cd atoms. Overall, the greatest improvement was achieved with co-doping, i.e. simultaneously substituting Te with Se and substituting Cd with Cu or Al. The elimination of midgap states is predicted to increase the photovoltaic efficiency of CdTe by reducing the recombination at grain boundaries.

14 SOLAR ENERGY↗

Non-reactive facet specific adsorption as a route to remediation of chlorinated organic contaminants

The present work quantifies metal-contaminant interactions between palladium substrates and three salient chlorinated organic contaminants, namely trichloroethylene 1,3,5-trichlorobenzene (TCB), and 3,3′,4,4′-tetrachlorobiphenyl (PCB77). Given that Pd is one of the conventional catalytically active materials known for contaminant removal, maximizing catalytic efficiency through optimal adsorption dynamics reduces the cost of remediation of contaminants that are persistent water pollutants chronically affecting public health. Adsorption efficiency analyses from all-atom molecular dynamics (MD) simulations advance the understanding of reaction mechanisms available from density functional theory (DFT) calculations to an extractable feature scale that can fit the parametric design of supported metal catalytic systems and feed into high throughput catalyst selection. Data on residence time, site-specific adsorption, binding energies, packing geometries, orientation profiles, and the effect of adsorbate size show the anomalous behaviour of organic contaminant adsorption on the undercoordinated {110} surface as compared to the {111} and {100} surfaces. The intermolecular interaction within contaminants from molecular dynamics simulation exhibits refreshing results than ordinary single molecule density functional theory calculation. Since complete adsorption and dechlorination is an essential step for chlorinated organic contaminant remediation pathways, the presented profiles provide essential information for designing efficient remediation systems through facet-controlled palladium nanoparticles.

Guo, Hao↗

Doping Effects on the Ductility of a Lightweight Refractory High-Entropy Alloy: Grain Boundary and Bulk Lattice Aspects

Doping elements in small amounts often segregate to grain boundaries (GBs) in alloys and can significantly impact mechanical properties and performance. Refractory high-entropy alloys (RHEAs) are known for their poor ductility, especially at low temperatures. Promoting GB cohesion through segregation can be an effective approach to mitigate embrittlement. Here, in this study, first-principles density functional theory (DFT) calculations were performed to examine the effects of important interstitial dopants (O, B, C, and N) and substitutional dopants (Cr, Y, La, and Ce) on the Σ5(310) [001] tilt GB of a lightweight RHEA Nb 32.5 Ti 27.5 Mo 22.5 Ta 12.5 Hf 2.5 Zr 2.5 . The DFT calculations reveal that certain dopants, such as B, C, Cr, N, and Y, exhibit favorable GB strengthening effects by improving bonding interactions with the bulk alloy. The impact of doping on the ductility parameter of the bulk lattice, defined as the ratio of surface energy to unstable stacking fault energy for the {110} <111> slip system, was also studied; and the results show that doping reduces the intrinsic ductility of the alloy, decreasing the D-parameter from 2.90 to 2.55, depending on the specific dopant. The present findings provide a foundational understanding at atomic level of the effect of representative dopants on mechanical properties of RHEAs and can be used to guide future alloy design for improved mechanical properties.

deformation charge density↗

Ternary PtIrNi Catalysts for Efficient Electrochemical Ammonia Oxidation

Due to high energy density, carbon-free feature, and easiness to be liquefied, ammonia (NH 3 ) has proved to be an effective alternative to hydrogen in low temperature fuel cells via its direct the ammonia oxidation reaction (AOR) for electricity generation. However, the kinetically sluggish AOR has prohibitively hindered the attractive direct ammonia fuel cell (DAFC) applications. Here we report an efficient AOR catalyst, in which ternary PtIrNi alloy nanoparticles well dispersed on a binary composite support consisting of porous silicon dioxide (SiO 2 ) and carboxyl-functionalized carbon nanotube (PtIrNi/SiO 2 -CNT-COOH) through a sonochemical-assisted synthesis strategy. The PtIrNi alloy nanoparticles, with the aid of abundant OHad provided by porous SiO 2 , and the improved electrical conductivity by CNTs, exhibit remarkable catalytic activity for the AOR in alkaline media. It is evidenced by a lower onset potential (~0.40 V vs. RHE) at room temperature, than that of commercial PtIr/C (ca. 0.43 V vs. RHE). Increasing NH 3 concentrations and operation temperatures significantly improve AOR performance of this catalyst. Specifically, AOR activity of the optimal PtIrNi nanoparticle catalyst can be significantly enhanced by elevating the temperature to 80ºC, with a much lower onset potential (~0.32 V vs. RHE), indicating that DAFC can be operated at higher temperature for increased performance. Constant-potential density functional theory (DFT) calculations showed that the Pt-Ir ensembles on {100}-terminated surfaces serve as the active site. Importantly, the introduction of Ni raises the center energy of the density of states projected onto the group d-orbitals of surface sites and thus lowers the theoretical onset potential for *NH 2 dehydrogenation to *NH when compared to Pt and Pt 3 Ir alloy.

25 ENERGY STORAGE↗

Rb 4 CuSb 2 Cl 11 and Rb 2 In 0.91(0.2) Sb 0.09 Cl 5 ·H 2 O: Wide Band Gap 0D Metal Halide Semiconductors

Herein, we report the discovery, structural and photophysical characterization of a new zero-dimensional (0D) lead-free all-inorganic halide, Rb 4 CuSb 2 Cl 11 , which adopts a new structure type. Single-crystal X-ray diffraction (SCXRD) shows that the structure consists of isolated, distorted seesaw [SbCl 4 ] − and trigonal planar [CuCl 3 ] 2− units, separated by Rb + cations that provide charge balance. Optoelectronic measurements and density functional theory (DFT) calculations indicate an indirect band gap of 2.89 eV, making it a candidate for wide-bandgap optoelectronic applications. Electrical resistivity was measured at 1.29 × 10 10 Ω·cm, and the trap-state density (n trap ) was found to be 7.44 × 10 10 cm −3 . Attempts to synthesize substitution analogs of Rb 4 CuSb 2 Cl 11 led to the synthesis of Rb 2 In 0.91(0.2) Sb 0.09 Cl 5 ·H 2 O, which was erroneously reported as Rb 2 SbCl 5 O in a previous study. Rb 2 In 0.91(0.2) Sb 0.09 Cl 5 ·H 2 O adopts a vacancyordered perovskite structure and exhibits broad-band yellow emission under UV excitation. The measured photoluminescence quantum yield (PLQY) for Rb 2 In 0.91 (0.2)Sb 0.09 Cl 5 ·H 2 O is 18.2%. These findings add to the growing class of quaternary metal halides with multiple cation and anion compositions, expanding the chemical phase space for the discovery of new materials with functional properties.

Crystal structure↗

A First-Principles Approach to Modeling Surface Site Stabilities on Multimetallic Catalysts

The study of multimetallic alloys and the multitude of possible surface compositions have sparked a tremendous interest in engineering low-cost materials with high activity and selectivity in heterogeneous catalysis. Multimetallic systems provide complementary functionalities and an unprecedented tunability when designing catalyst formulations. However, due to their immense structural and compositional complexity, the investigation and identification of an optimal catalyst is a tedious and time-consuming process, both experimentally and theoretically. Therefore, theoretical design principles are highly desirable to accelerate the screening of catalyst structures across the vast compositional space. In this paper, we introduce a simple and general model for predicting the site stability of multimetallic surfaces and nanoparticles, which is based on physical principles. The model requires only a small set of density functional theory (DFT) calculations of metal atom binding energies on monometallic and dilute alloy surface slabs to optimize the parameters in the simple model. The resulting model allows for the quantification of the stability of any particular atom site in any conceivable chemical environment across a wide range of morphologies, sizes, and arrangements by interpolating the derived parameters from a monometallic system to a completely diluted alloyed system. Herein, we demonstrate the robustness of the model across an extensive data set of transition metal alloy surfaces and 147-atoms cuboctahedral nanoparticles (NPs) composed of IrRhRu and PtPdRu. In conclusion, our approach yields mean absolute errors of ≈0.15 (IrRhRu), 0.20 (PtPdRu), 0.19 (IrRhRu NP), and 0.26 (PtPdRu NP) eV relative to site binding energies calculated using DFT.

36 MATERIALS SCIENCE↗