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

Was this helpful?

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

Bonding

How chemical bonds are determined in MPcules

PreviousAtomic Partial SpinsNextMetal Coordination and Binding

Last updated 1 year ago

Was this helpful?

Especially when relying on the chemical definition of a molecule (see ), it is important to define the bonds in a molecule. Bonds can include covalent bonds, meaning that electrons are shared between multiple atoms, or other interactions between atoms like ionic bonds, hydrogen bonds, and coordinate bonds.

In MPcules, we currently determine molecular bonding in three ways. The simplest way relies on the cheminformatics toolkit [1] and the defined in pymatgen. This method relies purely on valence- and distance-based heuristics, meaning that it can be used on any molecular structure, without any specific electronic structure calculations. Because of this, we rely on this OpenBabel/pymatgen method when .

In addition to the purely heuristic OpenBabel/pymatgen method, we also use the method of Spotte-Smith, Blau, et al. [2] and natural bonding orbital (NBO) analysis [3, 4]. The Spotte-Smith-Blau method begins with the heuristic bonds defined by OpenBabel and pymatgen. Then, applying the critic2 tool [5] we identify additional bonds as the critical points of the optimized electron density from a DFT calculation. More specifically, if there are any critical points between atoms with a field strength greater than 0.02 (in atomic units) where the distance between atoms is < 2.5 Å, we say that those atoms are bonded.

NBO reports bonds based on electron sharing in hybrid orbitals between atoms (that is, covalent bonds). In addition to these bonds that are directly output by NBO, we can infer electrostatic bonds via orbital interactions. Specifically, to identify coordinate bonds between metals and nonmetals from NBO, we examine NBO's second-order perturbation theory analysis. If there are interactions between nonmetal lone pair orbitals and metal lone vacant or anti-Rydberg orbitals where the distance between the metal and the nonmetal is < 3.0 Å and the perturbation energy for the orbital interaction is ≥ 3.0 kcal/mol, then we say that there is a bond between the metal and the nonmetal.

In both the Spotte-Smith-Blau method based on critical point analysis and the NBO method based partially on orbital interactions, the cutoff values (in terms of interatomic distance, field strength, and perturbation energy) were determined heuristically by closely analyzing the NBO outputs for a modest, quasi-random set of molecules from MPcules.

References

  1. O'Boyle, N.M., Banck, M., James, C.A., Morley, C., Vandermeersch, T. and Hutchison, G.R., 2011. Open Babel: An open chemical toolbox. Journal of cheminformatics, 3(1), pp.1-14.

  2. Spotte-Smith, E.W.C., Blau, S.M., Xie, X., Patel, H.D., Wen, M., Wood, B., Dwaraknath, S. and Persson, K.A., 2021. Quantum chemical calculations of lithium-ion battery electrolyte and interphase species. Scientific data, 8(1), p.203.

  3. Glendening, E.D., Badenhoop, J.K., Reed, A.E., Carpenter, J.E., Bohmann, J.A., Morales, C.M., Karafiloglou, P., Landis, C.R., Weinhold, F., 2018. NBO 7.0. Theoretical Chemistry Institute, University of Wisconsin, Madison.

  4. Glendening, E.D., Landis, C.R. and Weinhold, F., 2012. Natural bond orbital methods. Wiley interdisciplinary reviews: computational molecular science, 2(1), pp.1-42.\

  5. Otero-de-la-Roza, A., Johnson, E.R. and Luaña, V., 2014. Critic2: A program for real-space analysis of quantum chemical interactions in solids. Computer Physics Communications, 185(3), pp.1007-1018.

Molecules Methodology - Overview
OpenBabel
metal_edge_extender utility
defining molecules