<div>Hi Gabriele ,<br></div><div><br></div><div>as a quick workaround make the atoms you want to fix rigid ions.</div><div><br></div><div>Seems to me that the fixed ion routines (fix_atom_control) only zero the forces on particles, but the forces on cores and shells are separate thingies, so an extra bit of logic needs to go in. This only really matters if you wanted to have the shells relax, but the cores remain fixed - otherwise changing to rigid ions causes no problems.</div><div><br></div><div>Cheers,</div><div><br></div><div>Matt</div>