<div dir="ltr"><div>Hi Coralie,<br></div><div>Could you attach the output file? To see the forces and the geometry after the first move. Also, the pseudopotential could be useful to someone can run your input.</div><div><br></div><div>Your parameters are good, for paw pseudopotentials 50 Ry is an standart. <br></div><div>One thing you could change is the beta factor and the mixing mode, sometimes those are the key, an afternoon playing with it could solve the problem.</div><div>Another thing is the initial geometry, the methane molecule is very close to the slab... 1.6A is like an hydrogen bond, for methane I expect longer distance for a physisorption. If this is the problem, your forces in the first SCF step should be very big for methane and near slab atoms...</div><div><br></div><div>I do not have so much experience with magnetic systems, but I see in your output, the magnetization changes... maybe the initial magnetization is bad behaved and converges to weird density.</div><div><br></div><div>Regards<br></div></div>