First experiment using lammps with a pretrained potential (mace mp0) about mace HOT 6 OPEN

Can you try 'atom_style atomic'?

from mace.

Can you try 'atom_style atomic'?

Thanks, I tried and It kinda works. I forgot to update the question.

If I understand correctly, atomic is the style to be used because it is useless to explicitly define bonds and angle between atoms in molecules using this kind of potentials. Maybe I am wrong but if MACE mp0 is trained on potential coming DFT calculations, so an ab initio method, after a number of iterations, water molecules should "form" only because of the potential, without the need of explicitly fixing distances, bonds, etc.

I tried the following two scripts and it works, the problem is that it is apparently too slow and it segfaults for lack of memory I think, but I'm just trying on a small laptop with 16 GB of RAM...

Lammps input

Initial configuration

from mace.

Does it work on your laptop for a smaller system? Try the smallest that is reasonable.

from mace.

Does it work on your laptop for a smaller system? Try the smallest that is reasonable.

I tried with a 64 molecules systems, so small enough I think but I have still some segfaults problems.

This is the initial configuration, it is a water system that should be equilibrated

https://dpaste.org/Jr81n

This is the error:

https://dpaste.org/JYVtB

This it the lammps input script:

https://dpaste.org/nEPBU

I tried with both domain decomposition and with no domain decomposition.

Can you try and see if I'm doing something wrong? Thank you.

from mace.

I've tried your example and reproduced similar errors. As best I can tell, the main issue is that you are using a gigantic timestep (1 ps). Try something closer to 1 fs.

from mace.

I've tried your example and reproduced similar errors. As best I can tell, the main issue is that you are using a gigantic timestep (1 ps). Try something closer to 1 fs.

Thanks, that was the problem. With a timestep of 0.001 (1 fs) it seems to work.

Another question:
I'm running it on a cluster with four Intel(R) Xeon(R) Gold 6252N CPU @ 2.30GHz cpus. I found that the simulation still crashes if I run it with more than one processor, for example with mpirun -np 2 lmp but I have to experiment more.

The best performances are obtained with running a job with just one mpi process but using a large number of threads (96 for example setting the env variable with export OMP_NUM_THREADS=96) without domain decomposition, just following the lammps mace documentation. With this setting it takes 25 seconds to compute 1 timestep for a system of 125 water molecules starting from an equilibrated configuration using the mace mp0 L1 model. Is it normal to get this results?

By the way the next task is trying to compile the GPU version to run on my modest Nvidia RTX 3060Ti 8 GB on my Arch Linux system :-)

from mace.