Gold-doped graphene: A highly stable and active electrocatalysts for the oxygen reduction reaction

In addressing the growing need of renewable and sustainable energy resources, hydrogen-fuel-cells stand as one of the most promising routes to transform the current energy paradigm into one that integrally fulfills environmental sustainability. Nevertheless, accomplishing this technology at a large scale demands to surpass the efficiency and enhance the cost-effectiveness of platinum-based cathodes, which catalyze the oxygen reduction reaction (ORR). In this work, our first-principles calculations show that Au atoms incorporated into graphene di-vacancies form a highly stable and cost-effective electrocatalyst that is, at the same time, as or more (dependently of the dopant concentration) active toward ORR than the best-known Pt-based electrocatalysts. We reveal that partial passivation of defected-graphene by gold atoms reduces the reactivity of C dangling bonds and increases that of Au, thus optimizing them for catalyzing the ORR and yielding a system of high thermodynamic and electrochemical stabilities. We also demonstrate that the linear relation among the binding energies of the reaction intermediates assumed in computational high-throughput material screening does not hold, at least for this non-purely transition-metal material. We expect Au-doped graphene to finally overcome the cathode-related challenge hindering the realization of hydrogen-fuel cells as the leading means of powering transportation and portable devices

Friedel oscillations responsible for stacking fault of adatoms: The case of Mg(0001) and Be(0001)

We perform a first-principles study of Mg adatom and adislands on the Mg(0001) surface, and Be adatom on Be(0001), to obtain further insights into the previously reported energetic preference of the fcc faulty stacking of Mg monomers on Mg(0001). We first provide a viewpoint on how Friedel oscillations influence ionic relaxation on these surfaces. Our three-dimensional charge-density analysis demonstrates that Friedel oscillations have maxima which are more spatially localized than what one-dimensional average density or two-dimensional cross sectional plots could possibly inform: The well-known charge-density enhancement around the topmost surface layer of Mg(0001) is strongly localized at its fcc hollow sites. The charge accumulation at this site explains the energetically preferred stacking fault of the Mg monomer, dimer and trimer. Yet, larger islands prefer the normal hcp stacking. Surprisingly, the mechanism by which the fcc site becomes energetically more favorable is not that of enhancing the surface-adatom bonds but rather those between surface atoms. To confirm our conclusions, we analyze the stacking of Be adatom on Be(0001) - a surface also largely influenced by Friedel oscillations. We find, in fact, a much stronger effect: The charge enhancement at the fcc site is even larger and, consequently, the stacking-fault energy favoring the fcc site is quite large, 44 meV.Comment: Submitted to Physical Review

The perturbation energy: A missing key to understand the nobleness of bulk gold

The nobleness of gold surfaces has been appreciated since long before the beginning of recorded history. Yet, the origin of this phenomenon remains open because the so far existing explanations either incorrectly imply that silver should be the noblest metal or would fail to predict the dissolution of Au in aqua regia. Here, based on our analyses of oxygen adsorption, we advance that bulk gold\u27s unique resistance to oxidation is traced to the large energy cost associated with the perturbation its surfaces undergo upon adsorption of highly electronegative species. This fact is related to the almost totally filled d-band of Au and relativistic effects, but does not imply that the strength of the adsorbate-Au bond is weak. The magnitude of the structural and charge-density perturbation energy upon adsorption of atomic oxygen-which is largest for Au-is assessed from first-principles calculations and confirmed via a multiple regression analysis of the binding energy of oxygen on metal surfaces

Diffusion of the Cu monomer and dimer on Ag(111): Molecular dynamics simulations and density functional theory calculations

We present results of molecular dynamics (MD) simulations and density functional theory (DFT) calculations of the diffusion of Cu adatom and dimer on Ag(111). We have used potentials generated by the embedded-atom method for the MD simulations and pseudopotentials derived from the projected-augmented-wave method for the DFT calculations. The MD simulations (at three different temperatures: 300, 500, and 700 K) show that the diffusivity has an Arrhenius behavior. The effective energy barriers obtained from the Arrhenius plots are in excellent agreement with those extracted from scanning tunneling microscopy experiments. While the diffusion barrier for Cu monomers on Ag(111) is higher than that reported (both in experiment and theory) for Cu(111), the reverse holds for dimers [which, for Cu(111), has so far only been theoretically assessed]. In comparing our MD result with those for Cu islets on Cu(111), we conclude that the higher barriers for Cu monomers on Ag(111) results from the comparatively large Ag-Ag bond length, whereas for Cu dimers on Ag(111) the diffusivity is taken over and boosted by the competition in optimization of the Cu-Cu dimer bond and the five nearest-neighbor Cu-Ag bonds. Our DFT calculations confirm the relatively large barriers for the Cu monomer on Ag(111)-69 and 75 meV-compared to those on Cu(111) and hint a rationale for them. In the case of the Cu dimer, the relatively long Ag-Ag bond length makes available a diffusion route whose highest relevant energy barrier is only 72 meV and which is not favorable on Cu(111). This process, together with another involving an energy barrier of 83 meV, establishes the possibility of low-barrier intercell diffusion by purely zigzag mechanisms

High CO tolerance of Pt/Ru nano-catalyst: insight from first principles calculation

Density functional theory based calculations of the energetics of adsorption and diffusion of CO on Pt islets and on the Ru(0001) substrate show that CO has the lowest adsorption energy at the center of the islet, and its bonding increases as it moves to the edge of the island and further onto the substrate. Activation energy barriers for CO diffusion from the islet to the Ru surface are found to be lower than 0.3 eV making the process feasible and leading to the conclusion that this hydrogen oxidation catalyst is CO tolerant because of the spillover of CO from active Pt sites to the Ru substrate. We present the rationale for this effect using insights from detailed electronic structure calculations.Comment: 6 pages, 5 figure

First principles calculations of the electronic and geometric structure of Ag27Cu7Ag_{27}Cu_{7} nanoalloy

\emph{Ab initio} calculations of the structure and electronic density of states (DOS) of the perfect core-shell $Ag_{27}Cu_{7}$ nanoalloy attest to its $D_{5h}$ symmetry and confirm that it has only 6 non-equivalent (2 $Cu$ and 4 $Ag$) atoms. Analysis of bond-length, average formation energy, heat of formation of $Ag_{27}Cu_{7}$ and $L1_2$ $Ag-Cu$ alloys provide an explanation for the relative stability of the former with respect to the other nanoalloys in the same family. The HOMO-LUMO gap is found to be 0.77 eV, in agreement with previous results. Analysis of the DOS of $Ag_{27}Cu_{7}$, $L1_2$ $Ag-Cu$ alloys and related systems provides insight into the effects of low coordination, contraction/expansion and the presence of foreign atoms on the DOS of $Cu$ and $Ag$. While some characteristics of the DOS are reminiscent of those of the phonon-stable $L1_2$ $Ag-Cu$ alloys, the $Cu$ and $Ag$ states hybridize significantly in $Ag_{27}Cu_{7}$, compensating the $d$-band narrowing that each atom undergoes and hindering the dip in the DOS found in the bulk alloys. Charge density plots of $Ag_{27}Cu_{7}$ provide further insights into the relative strengths of the various interatomic bonds. Our results for the electronic and geometric structure of this nanoalloy can be explained in terms of length and strength hierarchies of the bonds, which may have implications also for the stability of alloy in any phase or size.Comment: 16 figure

Effect of dipolar interactions on the magnetization of a cubic array of nanomagnets

We investigated the effect of intermolecular dipolar interactions on a cubic 3D ensemble of 5X5X4=100 nanomagnets, each with spin $S = 5$. We employed the Landau-Lifshitz-Gilbert equation to solve for the magnetization $M(B)$ curves for several values of the damping constant $\alpha$, the induction sweep rate, the lattice constant $a$, the temperature $T$, and the magnetic anisotropy field $H_A$. We find that the smaller the $\alpha$, the stronger the maximum induction required to produce hysteresis. The shape of the hysteresis loops also depends on the damping constant. We find further that the system magnetizes and demagnetizes at decreasing magnetic field strengths with decreasing sweep rates, resulting in smaller hysteresis loops. Variations of $a$ within realistic values (1.5 nm - 2.5 nm) show that the dipolar interaction plays an important role in the magnetic hysteresis by controlling the relaxation process. The $T$ dependencies of $\alpha$ and of $M$ are presented and discussed with regard to recent experimental data on nanomagnets. $H_A$ enhances the size of the hysteresis loops for external fields parallel to the anisotropy axis, but decreases it for perpendicular external fields. Finally, we reproduce and test an $M(B)$ curve for a 2D-system [M. Kayali and W. Saslow, Phys. Rev. B {\bf 70}, 174404 (2004)]. We show that its hysteretic behavior is only weakly dependent on the shape anisotropy field and the sweep rate, but depends sensitively upon the dipolar interactions. Although in 3D systems, dipole-dipole interactions generally diminish the hysteresis, in 2D systems, they strongly enhance it. For both square 2D and rectangular 3D lattices with ${\bm B}||(\hat{\bm x}+\hat{\bm y})$, dipole-dipole interactions can cause large jumps in the magnetization.Comment: 15 pages 14 figures, submitted to Phys. Rev.

Effect of c(2x2)-CO overlayer on the phonons of Cu(001): a first principles study

We have examined the effect of a c(2x2) overlayer of CO on the surface phonons of the substrate, Cu(001), by applying the density functional perturbation theory with both the local (LDA) and the generalized-gradient (GGA) density approximations, through the Hedin-Lundqvist and the Perdew-Burke-Ernzerhof functionals, respectively. Our results (GGA) trace the Rayleigh wave softening detected by helium atom scattering (HAS) experiments to changes in the force constants between the substrate surface atoms brought about by CO chemisorption, resolving an ongoing debate on the subject. The calculated surface phonon dispersion curves document the changes in the polarization of some modes and show those of the modes originally along the $\bar{Y}$ direction of the clean surface Brillouin zone (SBZ) which are back-folded along the $\bar{\Delta}$ direction of the chemisorbed SBZ, to be particularly consequential. The vertical and shear horizontal section of $S_1$ in the SBZ of the clean surface, for example, is back-folded as a longitudinal-vertical mode, indicating thereby that $S_1$ $-$ predicted a long time back along $\bar{\Delta}$ for the clean surface $-$ may be indirectly assessed at $\bar{X}$ upon CO adsorption by standard planar scattering techniques. These findings further suggest that some of the energy losses detected by HAS along $\bar{\Delta}$, which were associated to multiphonon excitations of the adlayer frustrated translation mode, may actually correspond to the back-folded substrate surface modes

Vibrational dynamics of a c(2x2) phase induced by nitrogen adsorption on Cu(001)

Helium-atom scattering and density-functional perturbation theory (DFPT) calculations invoking the linear-response approximation and the pseudopotential approach have been used to study the vibrational dynamics of c(2x2)-like phases produced by nitrogen-ion implantation and subsequent annealing of Cu(001) surfaces. We find that, while the c(2x2) phonon dispersion relations are different from those of clean Cu(001), neither the acoustic nor the optical surface phonon mode energies measured along both [100] and [110] directions are dependent on N coverage once the c(2x2) pattern is formed. We show that the dispersion of the surface phonon modes is well reproduced with an analysis of the DFPT calculations of a nonstress relieved c(2x2) structure. A marked softening of a zone center optical mode is very apparent both experimentally and in the calculations. We show this softening arises largely because of interplanar Cu surface relaxations induced by N adsorption

Rational Catalyst Design Methodologies: Principles And Factors Affecting The Catalyst Design

The main focus of this chapter is the rational design of efficient and cost-effective electrocatalysts for the oxygen reduction reaction (ORR) on hydrogen fuel cell cathodes. The prohibitively high cost of the so-far best platinum-based ORR catalysts is one of the major obstacles for the large-scale application of hydrogen fuel cells. In addition, such catalysts have relatively low activity toward ORR and low electrochemical stability. Therefore, designing cost-effective, highly active, and stable ORR catalysts is of great importance to address the pressing need of renewable and sustainable energy resources. At present, the main research direction in such regard focuses on core-shell structures in which the shell is a monolayer of a catalytically active element (AE) and the core is an inexpensive metal substrate (MS). Therefore, it is critical to pin down the key factors determining the thermodynamic/electrochemical stability and catalytic activity of tailored catalysts such as structural match or mismatch between AE and MS and the hybridization between the AE and MS electronic states. As such, the rational selection of the AE and MS elements requires understanding of the tight relation between four variables: composition, electronic structure, binding energies of the ORR intermediates, and activity toward ORR . From our point of view, the least known link in the above chained dependence is the relation between electronic structure and binding energy. A number of models rationalizing that link (d-band center model, linear scaling between binding energies of different ORR intermediates, and relation between surface strain and its reactivity) have been shown to be insufficiently accurate to predict surface reactivity. Moreover, the electrochemical and thermodynamic stability of the AE/MS structures calls for a strong hybridization between the AE and MS electronic states, for which designing catalysts for practical applications requires fulfilling such demand. In this chapter, based on the wealth of information available from two decades of research, we will discuss the paths to efficiently design better ORR catalysts. Regarding the compromise of electrochemical stability (namely, most elements dissolve in the reaction environment) and ORR activity, we will discuss ways (a) to increase the dissolution potential of AE by depositing it on the proper MS; (b) to activate toward ORR some electrochemically stable elements (e.g., Au) through the AEMS interaction. The thermodynamic stability will also be rationalized in terms of the ratio between the strength of the AEAE, AEMS, and MSMS bonds, cohesive energy of AE and MS elements, and surface energies. Our approach to rational design of the ORR catalysts consists of (a) using the existing knowledge on the above-mentioned relations to preselect promising candidates for the ORR catalysts (a dozen, not hundreds or thousands); (b) performing first principles calculations to confirm and quantify the properties in question (in particular, thermodynamic and electrochemical stability and activity toward ORR), and narrow down the selection; (c) testing experimentally the systems found to be most promising in steps (a) and (b). We will use our recent results to illustrate the efficiency of this approach