Dear Lammps Users,
I am trying to fix pressure and temperature at 300 K for an ionic crystal: Li10GeP2S12. However, the temperature rises to about 650 K and cannot be fixed at 300K. As a new user the system does not let me to upload my files.
Here is my input file:
--------------------------Initialization
units metal
dimension 3
boundary p p p
atom_style charge
timestep 0.0001
atom_modify map array sort 0 0.0
neighbor 3.0 bin
#————————————————————
# System definition
read_data Li.data
replicate 1 1 1
#————————————————————
#Simulation setting
pair_style hybrid morse 6 coul/long 8.9
pair_coeff 1 4 morse 1.066825 1.531394 2.349741
pair_coeff 2 4 morse 2.087340 1.945525 2.104578
pair_coeff 3 4 morse 3.939695 1.766784 1.951558
pair_coeff 1 1 coul/long
pair_coeff 1 2 coul/long
pair_coeff 1 3 coul/long
pair_coeff 2 2 coul/long
pair_coeff 2 3 coul/long
pair_coeff 3 3 coul/long
pair_coeff 4 4 coul/long
kspace_style ewald 1.0e-6
#____________________
# Definition Li10GeP2S12 group
group Li type 1
group Ge type 2
group P type 3
group S type 4
# ---------- Define Settings ---------------------
compute eng all pe/atom
compute eatoms all reduce sum c_eng
#————————————————————
#Minimization
fix 1 all box/relax iso 0.0 vmax 0.001
thermo 400
thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms vol temp
min_style cg
minimize 1.0e-25 1.0e-25 10000 10000
unfix 1
#-------------------
# Equilibrium Setup
velocity all create 300 1357 mom yes rot yes
fix 2 all npt temp 300.0 300.0 100 iso 0.0 0.0 1000
dump 1 all atom 100 Equilibrium.atom
run 20000
unfix 2
undump 1
reset_timestep 0
#-----------------------
Could you please let me know where the problem is?
Many thanks in advance,
Marsusi,
Amirkabir University