mat-surface-energy

Calculate surface energies for (hkl) planes and construct equilibrium Wulff shapes.

144|21|Updated Jan 8, 2026
One-click install
npx skills add https://github.com/learningmatter-mit/AtomisticSkills --skill mat-surface-energy
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: mat-surface-energy
Source: https://github.com/learningmatter-mit/AtomisticSkills/tree/main/.agents/skills/mat-surface-energy
Command: npx skills add https://github.com/learningmatter-mit/AtomisticSkills --skill mat-surface-energy

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires pymatgen, ase, matplotlib, numpy, json, pathlib, and includes scripts (resource) components.

What problem does it solve?

It helps you quantify the surface energy of specific crystallographic planes (hkl) and use those values to predict the equilibrium crystal morphology (Wulff shape) for a material.

Core Features & Use Cases

  • Plane-resolved surface energy: Relaxes a bulk reference structure and multiple oriented slabs for different (hkl) terminations to compute γ for each plane.
  • MLIP/Model Context Protocol integration: Uses MLIP relaxations via MCP tools (e.g., MatGL/CHGNet) with consistent settings between bulk and slabs.
  • Automatic Wulff construction: Converts computed plane energies into a Wulff shape plot and associated geometric properties for comparison across materials or potential models.
  • Typical use case: Compare how different MLIPs (or different levels of theory) change predicted equilibrium nanoparticle shapes by recalculating γ and regenerating the Wulff diagram.

Quick Start

Run the surface workflow for your bulk CIF by creating slabs up to a chosen Miller index, relaxing bulk and slabs with the same MLIP settings, computing surface energies, and generating the Wulff shape image.

Frequently Asked Questions about mat-surface-energy

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I calculate surface energy for different crystallographic planes?

Calculating surface energy involves relaxing a bulk reference structure and multiple oriented slabs for different (hkl) terminations to compute the energy per area. This workflow requires consistent relaxation settings between bulk and slab models.

What is a Wulff shape and when do I need it for nanoparticle morphology studies?

A Wulff shape is an equilibrium crystal polyhedron constructed from plane-resolved surface energies. You need Wulff construction when predicting equilibrium nanoparticle morphology and catalyst shapes from specific crystallographic plane energies.

Can I use pymatgen to generate slabs and compute surface energies with MLIP relaxations?

Yes, pymatgen surface utilities generate oriented slabs and MLIP-based structural relaxations run via MCP tools. The workflow computes surface energies by keeping the slab cell area fixed during structural relaxation.

How do I construct a Wulff shape from relaxed bulk and slab energies?

You construct a Wulff shape by post-processing computed surface energies for each (hkl) plane. The workflow converts these plane-energy pairs into a Wulff polyhedron plot and calculates associated geometric properties.

Does this workflow support comparing different MLIPs for equilibrium crystal morphology?

Yes, the workflow supports comparing how different MLIPs change predicted equilibrium nanoparticle shapes. By recalculating surface energies and regenerating the Wulff diagram with consistent settings, you can compare results across potential models.

What are the limitations of using fixed slab cell area for structure relaxation?

Using a fixed slab cell area during structure relaxation ensures that the computed surface energy remains normalized to a constant geometric area. This constraint prevents the slab from undergoing lateral relaxation, which may limit accuracy for surfaces with significant in-plane stress.