Stability#
We’re going to start simple here - let’s run a local relaxation (optimize the unit cell and positions) using a pre-trained EquiformerV2-31M-OMAT24-MP-sAlex checkpoint. This checkpoint has a few fun properties
It’s a relatively small (31M) parameter model
It was pre-trained on the OMat24 dataset, and then fine-tuned on the MPtrj and Alexandria datasets, so it should emit energies and forces that are consistent with the MP GGA (PBE/PBE+U) level of theory
Need to install fairchem-core or get UMA access or getting permissions/401 errors?
Install the necessary packages using pip, uv etc
Get access to any necessary huggingface gated models
Get and login to your Huggingface account
Request access to https://huggingface.co/facebook/UMA
Create a Huggingface token at https://huggingface.co/settings/tokens/ with the permission “Permissions: Read access to contents of all public gated repos you can access”
Add the token as an environment variable using
huggingface-cli login
or by setting the HF_TOKEN environment variable.
from __future__ import annotations
import pprint
from ase.build import bulk
from ase.optimize import LBFGS
from quacc.recipes.mlp.core import relax_job
# Make an Atoms object of a bulk Cu structure
atoms = bulk("Cu")
# Run a structure relaxation
result = relax_job(
atoms,
method="fairchem",
name_or_path="uma-s-1",
task_name="omat",
opt_params={"fmax": 1e-3, "optimizer": LBFGS},
)
/home/runner/work/_tool/Python/3.12.11/x64/lib/python3.12/site-packages/torchtnt/utils/version.py:12: UserWarning: pkg_resources is deprecated as an API. See https://setuptools.pypa.io/en/latest/pkg_resources.html. The pkg_resources package is slated for removal as early as 2025-11-30. Refrain from using this package or pin to Setuptools<81.
import pkg_resources
WARNING:root:device was not explicitly set, using device='cuda'.
INFO:root:Setting random seed to 41
INFO:root:Setting random seed to 41
INFO:quacc.runners.prep:Calculation will run at /home/runner/work/fairchem/fairchem/docs/inorganic_materials/examples_tutorials/tmp-quacc-2025-07-01-20-21-39-844066-59941
INFO:quacc.runners.prep:Calculation results stored at /home/runner/work/fairchem/fairchem/docs/inorganic_materials/examples_tutorials/quacc-2025-07-01-20-21-39-844066-59941
pprint.pprint(result)
{'atoms': Atoms(symbols='Cu', pbc=True, cell=[[0.0, 1.805, 1.805], [1.805, 0.0, 1.805], [1.805, 1.805, 0.0]]),
'converged': np.True_,
'dir_name': '/home/runner/work/fairchem/fairchem/docs/inorganic_materials/examples_tutorials/quacc-2025-07-01-20-21-39-844066-59941',
'input_atoms': {'atoms': Atoms(symbols='Cu', pbc=True, cell=[[0.0, 1.805, 1.805], [1.805, 0.0, 1.805], [1.805, 1.805, 0.0]]),
'structure_metadata': {'builder_meta': {'build_date': datetime.datetime(2025, 7, 1, 20, 21, 40, 570271, tzinfo=datetime.timezone.utc),
'emmet_version': '0.84.8',
'pymatgen_version': '2025.6.14'},
'chemsys': 'Cu',
'composition': Composition('Cu1'),
'composition_reduced': Composition('Cu1'),
'density': 8.971719813140734,
'density_atomic': 11.761470249999999,
'elements': [Element Cu],
'formula_anonymous': 'A',
'formula_pretty': 'Cu',
'nelements': 1,
'nsites': 1,
'symmetry': {'angle_tolerance': 5.0,
'crystal_system': <CrystalSystem.cubic: 'Cubic'>,
'number': 225,
'point_group': 'm-3m',
'symbol': 'Fm-3m',
'symprec': 0.1,
'version': '2.6.0'},
'volume': 11.761470249999999}},
'name': 'fairchem Relax',
'nid': 'runnervm1w6u6.fmwxo1vhichurb1wybuswyps5c.dx.internal.cloudapp.net',
'parameters': {'version': '2.3.0'},
'parameters_opt': {'fmax': 0.001,
'max_steps': 1000,
'maxstep': 0.2,
'optimizer': 'LBFGS',
'type': 'optimization'},
'quacc_version': '0.16.0',
'results': {'energy': -3.7579006783217954,
'forces': array([[0., 0., 0.]], dtype=float32),
'free_energy': -3.7579006783217954,
'stress': array([-6.3947155e-03, -6.3947132e-03, -6.3947104e-03, 9.1077640e-10,
-2.4789384e-11, -2.3730695e-09], dtype=float32)},
'structure_metadata': {'builder_meta': {'build_date': datetime.datetime(2025, 7, 1, 20, 21, 40, 580795, tzinfo=datetime.timezone.utc),
'emmet_version': '0.84.8',
'pymatgen_version': '2025.6.14'},
'chemsys': 'Cu',
'composition': Composition('Cu1'),
'composition_reduced': Composition('Cu1'),
'density': 8.971719813140734,
'density_atomic': 11.761470249999999,
'elements': [Element Cu],
'formula_anonymous': 'A',
'formula_pretty': 'Cu',
'nelements': 1,
'nsites': 1,
'symmetry': {'angle_tolerance': 5.0,
'crystal_system': <CrystalSystem.cubic: 'Cubic'>,
'number': 225,
'point_group': 'm-3m',
'symbol': 'Fm-3m',
'symprec': 0.1,
'version': '2.6.0'},
'volume': 11.761470249999999},
'trajectory': [Atoms(symbols='Cu', pbc=True, cell=[[0.0, 1.805, 1.805], [1.805, 0.0, 1.805], [1.805, 1.805, 0.0]], calculator=SinglePointCalculator(...))],
'trajectory_results': [{'energy': -3.7579006783217954,
'forces': array([[0., 0., 0.]]),
'free_energy': -3.7579006783217954,
'stress': array([-6.39471551e-03, -6.39471319e-03, -6.39471039e-03, 9.10776399e-10,
-2.47893841e-11, -2.37306952e-09])}]}
Congratulations; you ran your first relaxation using an OMat24-trained checkpoint and quacc
!