Dear Lammps helpers,
I have been facing a “killed” without any warnings while I was trying to find the Melting point of Al. I faced this problem at temp 900 and 1000 C. I would appreciate assistance with why I get this error. Please find below my input script.
# Calculating the melting point of Al by LAMMPS
# *********************************** Initialization Section ****************************************
units metal
dimension 3
boundary p p p
# *********************************** Creating Atoms ************************************************
atom_style atomic
lattice fcc 4
region box block 0 1 0 1 0 1 units lattice
create_box 1 box
create_atoms 1 box
replicate 6 6 6
# ************************************* Interatomic Potential Settings ******************************
pair_style eam/fs
pair_coeff * * Al1.eam.fs Al
# ************************************* Run Minimization *******************************************
timestep 0.001
thermo 100
velocity all create 900.0 4928459 rot yes dist gaussian
fix 1 all box/relax iso 0.0 vmax 0.001
thermo_style custom step pe lx ly lz press pxx pyy pzz
min_style cg
minimize 1e-25 1e-25 5000 10000
run 1000
unfix 1
write_restart restart.boxmin
# ************************************* Applying temperature ***************************************
fix 2 all npt temp 1000.0 1000.0 100.0 iso 0.0 1.0 100.0
thermo_style custom step atoms temp press vol lx ly lz density time elaplong etotal enthalpy elapsed
dump 1 all xyz 10000 file.*.xyz
run 1000000
write_restart restart.equil
# ************************************* Output Section *********************************************
variable length1 equal "lx/6"
variable length2 equal "ly/6"
variable length3 equal "lz/6"
print "Lattice constant a (Angstoms) = \{length1\};"
print "Lattice constant b \(Angstoms\) = {length2};"
print "Lattice constant c (Angstoms) = ${length3};"
Meral