LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] How to raise the temperature after box relaxation at 0K?
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] How to raise the temperature after box relaxation at 0K?


From: Meral Sharkas <eng_meral.sh10@...16...>
Date: Tue, 22 Aug 2017 20:59:20 +0000 (UTC)

Dear Team,

I wander if you can help me with my question.

I am trying to relax my box at 0K to decrease pressure as much as I can as a start then I want to run my equilibration at 300 K and constant NPT ensemble.

When I ran the "fix/Npt" , I got Nans.

I suppose that the problem was "raising the temperature", so if that is the case, how can I raise the temperature to 300 K for equilibration?


timestep 0.0001
thermo 100
velocity      all create 0 7928459 dist gaussian
fix 1 all box/relax iso 0.0

#minimization

min_style cg
minimize 1e-25 1e-25 5000 100000
run 1000

#Equilibration
fix                     1 all npt temp 0 300.0 100 iso 0.0 0.0 1000
run 100000

Thank you in advance

Regards,