Lost atoms in lammps

Dear Rubén,

It’s because you don’t have any bonds defined between O and H atoms. You need to use the “create_bonds” command. The “set” command doesn’t create bonds.

Michal