Search NASA⌕ Search

SEARCH · Search NASA

Results for “GGA”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Tight-binding study of beryllium

In this work, we have applied the Naval Research Laboratory Tight-binding (NRL-TB) Method to the hcp alkaline-earth metal Beryllium. This approach consists of fitting to energy band and total energy data generated by the Linearized Augmented Plane Wave (LAPW) for the fcc, bcc, sc, and hcp structures as a function of volume. First, we found that using the Generalized Gradient Approximation (GGA) for the input data secures a better agreement to experimental volumes than the LDA. Second, we found that including small volumes to the fit is needed to get TB parameters that are transferable enough to successfully predict the energies of structures that were not fitted, and to perform molecular dynamics (MD) simulations. The non-orthogonal Hamiltonian obtained, is successful in correctly predicting the hcp lattice as the ground state of Be and finding eight other structures, that were not fitted to LAPW, positioned at higher energies. In addition, the NRL-TB produces accurate energy bands, densities of states, elastic constants, the Bain path and several quantities derived from MD simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theoretical predictive screening of noble-metal-containing M 3 AuC 2 (M = Ti, V, and Cr) MAX phases

Using first-principles density functional theory (DFT), the ground state physical properties of the newly synthesized noble-metal-containing Ti 3 AuC 2 MAX phase have been investigated. The effect of transition-element Ti replacement (V and Cr) on physical properties, including structural, elastic, electronic, thermal, and optical are presented. The optimized lattice parameters of Ti 3 AuC 2 MAX phase are in good agreement with the experimental values and decrease when we replace Ti with V and then Cr. The magnetic properties of M 3 AuC 2 are predicted by using generalized gradient approximation (GGA) incorporating onsite electron correlation parameter (Hubbard U parameter) within DFT. Out of three studied MAX phases, V 3 AuC 2 and Cr 3 AuC 2 phases exhibit magnetism. The thermodynamical stability of M 3 AuC 2 (M = Ti, V, and Cr) phases is discussed by the formation enthalpy with respect to their most competing phases. The electronic structure reveals that these phases have metallic nature and are electrically anisotropic. The M-element replacement has an effect on the bonding properties of M 3 AuC 2 . The bonding between M-C is in the order of Cr 3 AuC 2 > V 3 AuC 2 > Ti 3 AuC 2 according to their peak positions and heights of density of states in the occupied site, which is also confirmed from the charge density distribution. The elastic properties indicate that M 3 AuC 2 phases are ductile, machinable, less stiff, and better resistant to thermal shock and are elastically anisotropic. To understand the properties of M 3 AuC 2 for the extreme environment, we employed a quasi-harmonic Debye model at the pressure and temperature range of 0–50 GPa and 0–1600 K, respectively. Additionally, we evaluated thermal conductivity and melting temperature to further understand the potential of these MAX phases in coating applications at elevated temperature. In relation to optical properties, these MAX phases have a reasonable absorption coefficient in the visible as well as in the ultra-violet regions. The reflectivity of M 3 AuC 2 phases is polarization dependent and of 45 % against the visible region, which reveals its potential as coating material to minimize solar heating.

312 MAX phases↗
DOE OSTI2021

Structure sensitivity of ammonia electro-oxidation on transition metal surfaces: A first-principles study

Ammonia electro-oxidation is a promising catalytic reaction for application in alkaline fuel cells. In this paper we study the trends in this reaction on the (1 0 0) facet of eight fcc transition metals: Au, Ag, Cu, Pt, Pd, Ni, Ir, and Rh. We calculate from first-principles (DFT-GGA-PW91) the energetics for two mechanisms: (i) the N + N mechanism and (ii) the Gerischer-Mauerer mechanism. The onset potentials for both mechanisms are provided for both (1 0 0) and (1 1 1) facets. We also calculate the activation energies for the N-N bond-making non-Faradaic events on Cu, Pt, Pd, Ni, Ir, and Rh. We discover that N* – which is considered a poison for Pt(1 1 1) – is a reactive intermediate on Pt(1 0 0), allowing for a facile dimerization to N2 on Pt(1 0 0) at high potentials. We draw similar analyses for the structure sensitivity on the other metals, and we conclude with some general principles to guide catalyst design for this reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The microwave spectra of the conformers of $\mathcal{n}$-butyl nitrate

We report the microwave spectrum of n-butyl nitrate was recorded in the 5 to 20 GHz frequency range using broadband chirp and narrowband pulse excitation molecular jet Fourier transform microwave spectrometers. A quantum chemistry structural analysis yielded thirteen stable conformers. Among them, the five most energetically stable conformers were observed in the experimental spectra. The most stable conformer features a butyl chain with an anti-gauche-anti conformation (AGA) where the γ-carbon atom is about 64° out of the nitrate plane. For this conformer, spectra of all 13 C and 15 N minor isotopologues could be measured. The conformer with a straight butyl chain (AAA), and three other conformers (GAA, GGA, and AGG) were also observed. Accurate rotational constants, centrifugal distortion constants, and 14 N nuclear quadrupole coupling constants could be deduced and compared to the theoretical values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Small bimetallic clusters Agn-1M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15): Density functional theory and genetic algorithm

We investigated the effect of size and composition on the properties of bimetallic nanoclusters. The geometric structures, stabilities, and electronic properties of size-selected Ag n-1 M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15) bimetallic nanoclusters are systematically analyzed using spin-polarized density functional theory (DFT) within the generalized gradient approximation (GGA). We determine the most stable geometries for these clusters using a genetic algorithm (GA) in combination with DFT. Our results show that doping pure silver clusters with an M atom (transition metal), referred to as a “guest atom”, increases the stability as compared to pure Ag n (n = 3, 9, 15) clusters. The results for various properties including formation energy per atom, electronic structure, magnetic moments, and vibrational density of states (VDOS) are evaluated as a function of both size and composition of the system. The adsorption of selected bimetallic clusters on hydroxylated alumina substrate shows weak binding and minor changes in geometric properties except for Ag 8 Pt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Anionic Oxygen Redox in the High-Lithium Material Li 8 SnO 6

Lithiated transition metal oxide cathode materials that combine oxygen-anion redox [O 2 – /(O 2 ) n– ] with cation redox offer substantial capacities and higher voltage than cathodes without oxygen redox, but oxygen-anion redox is usually accompanied by the formation of peroxo or superoxo species that may cause oxygen release. To clarify the relationship between the oxygen release and peroxide and superoxide dimers without the complication of transition-metal redox, we performed density functional theory calculations to study the layered Li 8 SnO 6 cathode because it has only oxygen-anion redox during delithiation. Features with O–O distances of ~1.5 Å (corresponding to the formation of a peroxide) were observed at 75% lithium concentration, and further delithiation to 62.5% lithium concentration shortened the O–O bond distance to ~1.3 Å (corresponding to the formation of superoxide). The solely anionic oxygen redox stabilizes the delithiated structure of Li 8–x SnO 6 , and it does not cause structural disorder or release oxygen gas. Ab initio molecular dynamics calculations in the supercell at 300 K indicate that oxygen dimerization may occur at 75% lithium concentration. On the (0001) surface with a high lithium concentration (98%), oxygen dimerization is unfavorable. The calculations by HSE06 show that oxygen redox in Li 8 SnO 6 offers an initial voltage (vs Li/Li + ) greater than 4.2 V. The GGA+U calculations indicate that the Li + -polaron in Li 95 Sn 12 O 72 migrates with a migration barrier of only 0.43 eV. Here, this work also quantifies a solely anionic oxygen redox mechanism in Li 8 SnO 6 , and it provides a deeper understanding of anionic oxygen redox in Li-excess oxide cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Calculating and Characterizing the Charge Distributions in Solids

Accurate estimation of the partial atomic charges on metal centers is useful for understanding electronic and catalytic properties of materials. However, different methods of calculating these charges may give quite different results; this issue has been more widely studied for molecules than for solids. Here we study the charges on the metal centers of a test set of 18 solids containing transition metals by using density functional theory with several density functionals (PBE, PBE+U, TPSS, revTPSS, HLE17, revM06-L, B3LYP, B3LYP*, and other exchange-modified B3LYP functionals) and four charge models (Bader, Hirshfeld, CM5, and DDEC6). The test set contains 12 systems with nonmagnetic metal centers (eight metal oxides (MO 2 ), two metal sulfides (MS 2 ), and two metal selenides (MSe 2 )) and six ferromagnetic transition metal complexes. Our study shows that, among the four types of charges, Bader charges are the highest and Hirshfeld charges are the lowest for all the systems, regardless of the functional being used. The CM5 charges are bigger than DDEC6 charges for MX 2 with M = Ti or Mo and X = S or Se, but for the other 14 cases they are lower. We found that the most of the systems are sensitive to the Hubbard U parameters in PBE+U and to the percentage X of Hartree–Fock exchange in exchange-modified B3LYP; as we increase U or X, the charges on the metal atoms in MX 2 increase steadily. Testing different density functionals shows charges calculated with higher Hubbard U parameters in PBE+U are comparable to B3LYP (with 20% Hartree–Fock exchange). Among four meta-GGA functionals studied, the charges with HLE17 have the closest agreement with B3LYP. In conclusion, the variation of charges with choice of charge model is greater than the variation with choice of density functional.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Norm-Conserving Pseudopotentials and Basis Sets to Explore Actinide Chemistry in Complex Environments

We constructed and tested a new set of norm-conserving pseudopotentials, as well as companion Gaussian basis sets, for the actinide series (Ac–Lr) using the Goedecker, Teter and Hutter (GTH) formalism within the generalized gradient approximation (GGA) functional for exchange-correlation of Perdew, Burke and Ernzerhof (PBE). This fills a critical void in the literature for studying the chemistry of 5f–block elements in condensed phase. The accuracy and reliability of the newly parameterized An–GTH pseudopotentials and basis sets, a variety of tests on actinide-containing molecules were carried out and compared to all-electron and available experimental results. The new pseudopotentials include both medium ( [Xe]4f 14 ) and large core ( [Xe]4f 14 5d 10 ) options that successfully reproduced structure and energetics, particularly redox processes. The medium core size set, in particular, reproduced all electron calculations over multiple oxidation states from 0 to VII, whereas the large core set was suitable only for the early series elements and low oxidation states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: A Comparison Between Energy - and Variance-based Approaches

An efficient implementation of geometrical derivatives at the Hartree–Fock (HF) and current-density functional theory (CDFT) levels is presented for the study of molecular structure in strong magnetic fields. The required integral derivatives are constructed using a hybrid McMurchie–Davidson and Rys quadrature approach, which combines the amenability of the former to the evaluation of derivative integrals with the efficiency of the latter for basis sets with high angular momentum. In addition to its application to evaluating derivatives of four-center integrals, this approach is also applied to gradients using the resolution-of-the-identity approximation, enabling efficient optimization of molecular structure for many-electron systems under a strong magnetic field. The CDFT contributions have been implemented for a wide range of density functionals up to and including the meta-GGA level with current-density dependent contributions and (range-separated) hybrids for the first time. Illustrative applications are presented to the OH and benzene molecules, revealing the rich and complex chemistry induced by the presence of an external magnetic field. Challenges for geometry optimization in strong fields are highlighted, along with the requirement for careful analysis of the resulting electronic structure at each stationary point. The importance of correlation effects is examined by comparison of results at the HF and CDFT levels. The present implementation of molecular gradients at the CDFT level provides a cost-effective approach to the study of molecular structure under strong magnetic fields, opening up many new possibilities for the study of chemistry in this regime.

Comparison↗
DOE OSTI2022

Accurate Free Energies for Complex Condensed-Phase Reactions Using an Artificial Neural Network Corrected DFTB/MM Methodology

Semiempirical methods like density functional tight-binding (DFTB) allow extensive phase space sampling, making it possible to generate free energy surfaces of complex reactions in condensed-phase environments. Such a high efficiency often comes at the cost of reduced accuracy, which may be improved by developing a specific reaction parametrization (SRP) for the particular molecular system. Thiol–disulfide exchange is a nucleophilic substitution reaction that occurs in a large class of proteins. Its proper description requires a high-level ab initio method, while DFT-GAA and hybrid functionals were shown to be inadequate, and so is DFTB due to its DFT-GGA descent. We develop an SRP for thiol–disulfide exchange based on an artificial neural network (ANN) implementation in the DFTB+ software and compare its performance to that of a standard SRP approach applied to DFTB. Furthermore, as an application, we use both new DFTB-SRP as components of a QM/MM scheme to investigate thiol–disulfide exchange in two molecular complexes: a solvated model system and a blood protein. Demonstrating the strengths of the methodology, highly accurate free energy surfaces are generated at a low cost, as the augmentation of DFTB with an ANN only adds a small computational overhead.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

DeePKS + ABACUS as a Bridge between Expensive Quantum Mechanical Models and Machine Learning Potentials

Recently, the development of machine learning (ML) potentials has made it possible to perform large-scale and long-time molecular simulations with the accuracy of quantum mechanical (QM) models. However, for different levels of QM methods, such as density functional theory (DFT) at the meta-GGA level and/or with exact exchange, quantum Monte Carlo, etc., generating a sufficient amount of data for training an ML potential has remained computationally challenging due to their high cost. In this work, we demonstrate that this issue can be largely alleviated with Deep Kohn–Sham (DeePKS), an ML-based DFT model. DeePKS employs a computationally efficient neural network-based functional model to construct a correction term added upon a cheap DFT model. Upon training, DeePKS offers closely matched energies and forces compared with high-level QM method, but the number of training data required is orders of magnitude less than that required for training a reliable ML potential. As such, DeePKS can serve as a bridge between expensive QM models and ML potentials: one can generate a decent amount of high-accuracy QM data to train a DeePKS model and then use the DeePKS model to label a much larger amount of configurations to train an ML potential. Further, this scheme for periodic systems is implemented in a DFT package ABACUS, which is open source and ready for use in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Critical Assessment of the Thermodynamics of Vacancy Formation in Fe 2 O 3 Using Hybrid Density Functional Theory

Fe 2 O 3 hematite is a technologically important material with applications in energy storage as well as being the key phase formed in the rust of iron-based materials. Despite this central importance, there is much that is still unknown regarding the properties of defects in Fe 2 O 3 and consequently oxide growth. In this work, using screened hybrid density functional theory (HSE06), we consider the thermodynamics of vacancies in Fe 2 O 3 , considering the effects of ionic and electronic chemical potentials on both iron and oxygen vacancy formation. We find that, in the oxygen-rich limit, iron vacancies are easier to form, though the difference in formation energy between the two vacancies is only about 1 eV. In contrast, in the Fe-rich limit, oxygen vacancies have an extremely low formation energy, only 0.07 eV at mid-gap, and would spontaneously form as the Fermi level is reduced, while Fe vacancies require over 5 eV to form. Consistent with experiment, this indicates that Fe 2 O 3 is relatively easily reduced but not oxidized. However, the theoretical picture is very different when considering other exchange–correlation functionals (GGA + U or SCAN), emphasizing the critical role of the exchange–correlation functional in describing this system accurately.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Quantifying Temperature Dependence of Electronic Band Gaps and Optical Properties in SnO 2 and SnO via First-Principles Simulations

Tin metal oxides SnO x ( x = 1, 2) have gained interest as gas-sensing materials. For their applications as high-temperature sensors, however, a better understanding of their temperature dependence sensing responses is needed. Here, we comparatively quantify the temperature-dependent electronic band gaps and optical properties of SnO 2 and SnO using first-principles calculations. Without considering the temperature effect, SnO 2 and SnO are predicted to have direct and indirect band gaps of 2.18 and 1.66 eV, respectively, at the PBE + U -GGA level that we employed. The temperature effect on the electronic and optical properties is captured by taking account of the electron–phonon interaction. Band gap renormalization with temperature is calculated via the Allen–Heine–Cardona theory. For both oxides, we find a monotonic decrease in the electronic band gap such that renormalization at zero point (0 K) is ~–0.17 and ~–0.52 eV at 1000 K. These results are also analyzed by employing an analytical equation that helps characterize the band gap shift with temperature. In addition, the optical properties at finite temperatures are simulated using the frozen-phonon method that combines electron–phonon coupling with the momentum matrix. As temperature increases, the optical property spectra are smoothed because of the smearing effect, which diminishes optical constants at shorter wavelengths. Our results are of interest for high-temperature functional materials in applications of optical detection.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Trends in Formic Acid Electro-Oxidation on Transition Metals Alloyed with Platinum and Palladium

Direct formic acid fuel cells (DFAFCs) have emerged as a promising power source to meet increased demands for alternative energy sources in the transportation and portable energy storage sectors. Furthermore, these fuel cells utilize formic acid (FA), a nontoxic and carbon-neutral fuel when produced from biomass or via CO 2 reduction. Despite the promise of DFAFCs, the best monometallic catalysts, platinum and palladium, are poisoned by CO through the indirect oxidation pathway and require large overpotentials. By alloying Pt and Pd with other metals, we aim to improve both the activity and selectivity of these catalysts. Here, we present a systematic density functional theory (GGA-PW91) study of FA electro-oxidation (FAO) on the (111) facet of bimetallic Pt (Pt 3 M) or Pd (Pd 3 M) catalysts (M = Au, Ag, Cu, Pt, Pd, Ir, Rh, Ru, or Re) to evaluate the catalytic performance of these surfaces. For each surface, we calculate free energy diagrams and onset potentials of three key reaction mechanisms: direct oxidation of FA via carboxyl (COOH), direct oxidation of FA via formate (HCOO), and the indirect oxidation of FA that first forms CO en route to full oxidation to CO 2 . We then display the trends in the form of phase diagrams that compare the activity of the calculated surfaces against regions of high activity using the free energies of adsorbed CO and OH as descriptors, enabling high-throughput screening and design of improved catalysts, particularly those alloying Pt or Pd with Ir, Ru, or Re.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗