PyIron+Sphinx中的自旋约束

Question

我想限制块体原子的自旋，同时让我的超级细胞的自由表面原子放松它们的磁矩。在PyIron+SPhinx中，是否有可能约束超级细胞中原子子集(不是所有原子)的自旋？

Answer 1

原则上是可以的：

from pyiron import Project
import numpy as np

spx = pr.create.job.Sphinx('spx')
spx.structure = pr.create.structure.bulk('Fe', a=2.83, cubic=True)
spx.structure.set_initial_magnetic_moments([2, 2])
spx.fix_spin_constraint = True
spx.structure.spin_constraint = np.array([True, False])
spx.calc_static()
spx.run()

简短说明：spx.fix_spin_constraint = True初始化属性spx.structure.spin_constraint，该属性在开头只包含所有原子的True。对于不受约束的原子，可以设置False。