No Arabic abstract
Density Functional Theory (DFT) calculations have been widely used to predict the activity of catalysts based on the free energies of reaction intermediates. The incorporation of the state of the catalyst surface under the electrochemical operating conditions while constructing the free energy diagram is crucial, without which even trends in activity predictions could be imprecisely captured. Surface Pourbaix diagrams indicate the surface state as a function of the pH and the potential. In this work, we utilize error-estimation capabilities within the BEEF-vdW exchange correlation functional as an ensemble approach to propagate the uncertainty associated with the adsorption energetics in the construction of Pourbaix diagrams. Within this approach, surface-transition phase boundaries are no longer sharp and are therefore associated with a finite width. We determine the surface phase diagram for several transition metals under reaction conditions and electrode potentials relevant for the Oxygen Reduction Reaction (ORR). We observe that our surface phase predictions for most predominant species are in good agreement with cyclic voltammetry experiments and prior DFT studies. We use the OH$^*$ intermediate for comparing adsorption characteristics on Pt(111), Pt(100), Pd(111), Ir(111), Rh(111), and Ru(0001) since it has been shown to have a higher prediction efficiency relative to O$^*$, and find the trend Ru>Rh>Ir>Pt>Pd for (111) metal facets, where Ru binds OH$^*$ the strongest. We robustly predict the likely surface phase as a function of reaction conditions by associating c-values to quantifying the confidence in predictions within the Pourbaix diagram. We define a confidence quantifying metric using which certain experimentally observed surface phases and peak assignments can be better rationalized.
The success of descriptor-based material design relies on eliminating bad candidates and keeping good candidates for further investigation. While DFT has been widely successfully for the former, often times good candidates are lost due to the uncertainty associated with the DFT-predicted material properties. Uncertainty associated with DFT predictions has gained prominence and has led to the development of exchange correlation functionals that have built-in error estimation capability. In this work, we demonstrate the use of built-in error estimation capabilities within the BEEF-vdW exchange correlation functional for quantifying the uncertainty associated with the magnetic ground state of solids. We demonstrate this approach by calculating the uncertainty estimate for the energy difference between the different magnetic states of solids and compare them against a range of GGA exchange correlation functionals as is done in many first principles calculations of materials. We show that this estimate reasonably bounds the range of values obtained with the different GGA functionals. The estimate is determined as a post-processing step and thus provides a computationally robust and systematic approach to estimating uncertainty associated with predictions of magnetic ground states. We define a confidence value (c-value) that incorporates all calculated magnetic states in order to quantify the concurrence of the prediction at the GGA level and argue that predictions of magnetic ground states from GGA level DFT is incomplete without an accompanying c-value. We demonstrate the utility of this method using a case study of Li and Na-ion cathode materials and the c-value metric correctly identifies that GGA level DFT will have low predictability for NaFePO$_4$F.
Doped transition-metal dichalcogenides monolayers exhibit exciting magnetic properties for the benefit of two-dimensional spintronic devices. Using density functional theory (DFT) incorporating Hubbard-type of correction (DFT$+U$) to account for the electronic correlation, we study the magnetocrystalline anisotropy energy (MAE) characterizing Mn-doped MS$_2$ (M=Mo, W) monolayers. A single isolated Mn dopant exhibits a large perpendicular magnetic anisotropy of 35 meV (8 meV) in the case of Mn-doped WS$_2$ (MoS$_2$) monolayer. This value originates from the Mn in-plane orbitals degeneracy lifting due to the spin-orbit coupling. In pairwise doping, the magnetization easy axis changes to the in-plane direction with a weak MAE compared to single Mn doping. Our results suggest that diluted Mn-doped MS$_2$ monolayers, where the Mn dopants are well separated, could potentially be a candidate for the realization of ultimate nanomagnet units.
This work explores the use of joint density-functional theory, a new form of density-functional theory for the ab initio description of electronic systems in thermodynamic equilibrium with a liquid environment, to describe electrochemical systems. After reviewing the physics of the underlying fundamental electrochemical concepts, we identify the mapping between commonly measured electrochemical observables and microscopically computable quantities within an, in principle, exact theoretical framework. We then introduce a simple, computationally efficient approximate functional which we find to be quite successful in capturing a priori basic electrochemical phenomena, including the capacitive Stern and diffusive Gouy-Chapman regions in the electrochemical double layer, quantitative values for interfacial capacitance, and electrochemical potentials of zero charge for a series of metals. We explore surface charging with applied potential and are able to place our ab initio results directly on the scale associated with the Standard Hydrogen Electrode (SHE). Finally, we provide explicit details for implementation within standard density-functional theory software packages at negligible computational cost over standard calculations carried out within vacuum environments.
The metal-insulator transition observed in the In/Si(111)-4x1 reconstruction is studied by means of ab initio calculations of a simplified model of the surface. Different surface bands are identified and classified according to their origin and their response to several structural distortions. We support the, recently proposed [New J. of Phys. 7 (2005) 100], combination of a shear and a Peierls distortions as the origin of the metal-insulator transition. Our results also seem to favor an electronic driving force for the transition.
Transition-metal dichalcogenide IrTe2 has attracted attention because of striped lattice, charge ordering and superconductivity. We have investigated the surface structure of IrTe2, using low energy electron diffraction (LEED) and scanning tunneling microscopy (STM). A complex striped lattice modulations as a function of temperature is observed, which shows hysteresis between cooling and warming. While the bulk 5x1 and 8x1 phases appear at high temperatures, the surface ground state has the 6x1 phase, not seen in the bulk, and the surface transition temperatures are distinct from the bulk. The broken symmetry at the surface creates a quite different phase diagram, with the coexistence of several periodicities resembling devils staircase behavior.