Materials Project Documentation
Return to materialsproject.org
  • Introduction
  • Frequently Asked Questions (FAQ)
    • Glossary of Terms
  • Changes and Updates
    • Database Versions
    • Website Changelog
  • Documentation Credit
  • Community
    • Getting Help
    • Getting Involved
      • Contributor Guide
      • Potential Collaborators
      • MP Community Software Ecosystem
    • Community Resources
    • Code of Conduct
  • Services
    • MPContribs
  • Methodology
    • Materials Methodology
      • Overview
      • Calculation Details
        • GGA/GGA+U Calculations
          • Parameters and Convergence
          • Hubbard U Values
          • Pseudo-potentials
        • r2SCAN Calculations
          • Parameters and Convergence
          • Pseudopotentials
      • Thermodynamic Stability
        • Energy Corrections
          • Anion and GGA/GGA+U Mixing
          • GGA/GGA+U/r2SCAN Mixing
        • Phase Diagrams (PDs)
        • Chemical Potential Diagrams (CPDs)
        • Finite Temperature Estimation
      • Electronic Structure
      • Phonon Dispersion
      • Diffraction Patterns
      • Aqueous Stability (Pourbaix)
      • Magnetic Properties
      • Elastic Constants
      • Piezoelectric Constants
      • Dielectric Constants
      • Equations of State (EOS)
      • X-ray Absorption Spectra (XAS)
      • Surface Energies
      • Grain Boundaries
      • Charge Density
      • Suggested Substrates
      • Related Materials
      • Optical absorption spectra
      • Alloys
    • Molecules Methodology
      • Overview
      • Calculation Details
      • Atomic Partial Charges
      • Atomic Partial Spins
      • Bonding
      • Metal Coordination and Binding
      • Natural Atomic and Molecular Orbitals
      • Redox and Electrochemical Properties
      • Molecular Thermodynamics
      • Vibrational Properties
      • Legacy Data
    • MOF Methodology
      • Calculation Parameters
        • DFT Parameters
        • Density Functionals
        • Pseudopotentials
        • DFT Workflow
  • Apps
    • Explore and Search Apps
      • Materials Explorer
        • Tutorial
      • Molecules Explorer
        • Tutorial
        • Legacy Data
      • Battery Explorer
        • Background
        • Tutorial
      • Synthesis Explorer
        • Background
        • Tutorial
      • Catalysis Explorer
        • Tutorial
      • MOF Explorer
        • Downloading the Data
        • Structure Details
          • QMOF IDs
          • Structure Sources
          • Finding MOFs by Common Name
          • Structural Fidelity
        • Property Definitions
          • SMILES, MOFid, and MOFkey
          • Pore Geometry
          • Topology
          • Electronic Structure
          • Population Analyses and Bond Orders
          • Symmetry
        • Version History
        • How to Cite
    • Analysis Apps
      • Phase Diagram
        • Background
        • Tutorials
        • FAQ
      • Pourbaix Diagram
        • Background
        • Tutorial
        • FAQ
      • Crystal Toolkit
        • Background
        • Tutorial
        • FAQ
      • Reaction Calculator
      • Interface Reactions
    • Characterization Apps
      • X-ray Absorption Spectra (XAS)
    • Explore Contributed Data
  • Downloading Data
    • How do I download the Materials Project database?
    • Using the API
      • Getting Started
      • Querying Data
      • Tips for Large Downloads
      • Examples
      • Advanced Usage
    • Differences between new and legacy API
    • Query and Download Contributed Data
    • AWS OpenData
  • Uploading Data
    • Contribute Data
  • Data Production
    • Data Workflows
    • Data Builders
Powered by GitBook
On this page
  • Introduction
  • Formalism
  • DFT Parameters
  • References

Was this helpful?

Edit on GitHub
Export as PDF
  1. Methodology
  2. Materials Methodology

Surface Energies

PreviousX-ray Absorption Spectra (XAS)NextGrain Boundaries

Last updated 2 years ago

Was this helpful?

Introduction

Surface energy is a measure of the energy change associated with the breaking of intermolecular bonds in a bulk material to create a surface. In thermodynamically stable materials, the creation of a surface will always increase energy, otherwise there would be a thermodynamic driving force to create surfaces and the material would sublimate. In theory, surface energy is equal to half of the energy of cohesion (the energy needed to break all of the bonds required to form two new surfaces). However, this perfect cleaving of surfaces is rarely achieved. In reality, surfaces often rearrange and/or react with their surroundings to passivate or adsorb molecules or atoms to lower their surface energy from the theoretical cohesive energy value.

Formalism

Surface energy is calculated using a slab model where a supercell of a crystal is oriented such that a given facet of interest is created and then exposed to vacuum by removing atoms from the supercell. If we are interested in creating a surface with the plane (hkl) exposed, lattice vector transformations are performed on the supercell with lattice vectors a and b parallel to the exposed plane (hkl) and lattice vector c as close to perpendicular to the exposed plane as is feasbile. This new unit cell is referred to as the oriented unit cell. The atoms in the oriented unit cell must also be shifted in the c direction in order to expose all possible symmetrically distinct atomic terminations. This algorithm for generating slabs is implemented in pymatgen .

The surface energy γhklσ\gamma^{\sigma}_{hkl}γhklσ​ of facet (hkl) of a slab model is calculated as:

γhklσ=Eslabhkl,σ−Ebulkhkl⋅nslab2⋅Aslab\gamma^{\sigma}_{hkl} = \frac{E^{hkl,\sigma}_{slab} - E^{hkl}_{bulk} \cdot n_{slab}}{2 \cdot A_{slab}}γhklσ​=2⋅Aslab​Eslabhkl,σ​−Ebulkhkl​⋅nslab​​

where Eslabhkl,σE^{hkl,\sigma}_{slab}Eslabhkl,σ​​is the total energy of the slab with termination σ\sigmaσ​, EbulkhklE^{hkl}_{bulk}Ebulkhkl​​is the per atom total energy of the bulk oriented unit cell, nslabn_{slab}nslab​ is the total number of atoms in the slab and AAA is the surface area of the slab. The bulk oriented unit cell's atomic positions as well as its volume are relaxed, whereas in the slab model, only the atomic positions are relaxed.

DFT Parameters

All DFT calculations are performed in the Vienna Ab-initio Simulation Package (VASP) with the projector augmented wave (PAW) method. Exchange correlation effects are modeled using the Perdew-Berke-Ernzerhof (PBE) generalized gradient approximation (GGA) funcitonal. All calculations are spin polarized using a plane wave cutoff energy of 400eV. Full details can be found in .

References

[1]: Ong, S. P. et al. Python Materials Genomics (pymatgen): A robust, open-source python library for materials analysis. Computational Materials Science 68, 314–319 (2013)

[2]: Tran, R., Xu, Z., Radhakrishnan, B. et al. Surface energies of elemental crystals. Sci Data 3, 160080 (2016)

[1]
[2]