Skip to content

Commit 9906d67

Browse files
ElliottKasoarDomantas Kuryla
authored andcommitted
Set missing spin and charge for Orb (ddmms#384)
1 parent 78d2f6f commit 9906d67

2 files changed

Lines changed: 8 additions & 0 deletions

File tree

ml_peg/calcs/molecular_crystal/CPOSS209/calc_CPOSS209.py

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -70,6 +70,9 @@ def test_lattice_energy(mlip: tuple[str, Any]) -> None:
7070

7171
# Read crystal structure
7272
solid = read(crystal_path, index=0)
73+
# Set default charge and spin
74+
solid.info.setdefault("charge", 0)
75+
solid.info.setdefault("spin", 1)
7376
solid.calc = copy(calc)
7477
solid.get_potential_energy()
7578

@@ -126,6 +129,8 @@ def test_lattice_energy(mlip: tuple[str, Any]) -> None:
126129

127130
# Read gas phases
128131
molecule = read(molecule_path, index=0)
132+
molecule.info.setdefault("charge", 0)
133+
molecule.info.setdefault("spin", 1)
129134
molecule.calc = copy(calc)
130135
molecule.get_potential_energy()
131136

ml_peg/calcs/physicality/diatomics/calc_diatomics.py

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -171,6 +171,9 @@ def run_diatomics(model_name: str, model) -> None:
171171
[element1, element2],
172172
positions=[(0.0, 0.0, 0.0), (0.0, 0.0, float(distance))],
173173
)
174+
# Set default charge and spin
175+
atoms.info.setdefault("charge", 0)
176+
atoms.info.setdefault("spin", 1)
174177

175178
atoms.calc = calc
176179
energy = float(atoms.get_potential_energy())

0 commit comments

Comments
 (0)