Dear Lammps,
I am trying to minimize a configuration of ZrCu using the pair_style eam/alloy and suitable parameters for Zr and Cu (EAM potentials generated by Howard Sheng at GMU)
I would like that all forces of particles are very small, specifically
I want that the maxforce is < 1e-10.
However there is no way to minimize such configuration, I have tried with different minimization protocols, and procedures: cg, fire, damped dynamics.
Could you please give me an advice on what i am doing wrong?
I paste below my script.
Thank you very much for your kind attention,
kind regards.
###############################################
units metal
atom_style atomic
read_restart glass_50.restart.25502036
Interatomic potential information.
neigh_modify every 1000 delay 0 check no
pair_style eam/alloy
pair_coeff * * …/ZrCuAl.lammps.eam Zr Cu Al
timestep ${dt}
thermo 1000
thermo_style custom step temp pe ke fmax fnorm ndanger press pxx pyy pzz
thermo_modify flush yes format float %10.10g
variable maxforce equal abs(fmax)
fix HALT all halt 1 v_maxforce < 1e-10 error continue
fix ti2 all nve
fix ti1 all langevin 0.0 0.0 1.0 7774
run 100000
unfix ti1
unfix ti2
min_style fire
minimize 0.0 1e-10 1000000 1000000
min_style cg
minimize 0.0 1e-10 1000000 1000000
min_style fire
minimize 0.0 1e-10 1000000 1000000
min_style cg
minimize 0.0 1e-10 1000000 1000000
write_dump all custom final_20000.dump id type x y z fx fy fz modify append no sort id format line ‘%d %d %20.20g %20.20g %20.20g %20.20g %20.20g %20
.20g’
write_data final_20000.dat