How do I confirm that the Molecular dynamics run right or wrong?

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
utpalsarkar
Newbie
Newbie
Posts: 2
Joined: Thu Oct 30, 2008 9:00 pm
License Nr.: 90, upgraded to paw
Location: France

How do I confirm that the Molecular dynamics run right or wrong?

#1 Post by utpalsarkar » Tue Apr 28, 2009 7:50 am

Dear VASP user,
I am doing a NVE molecular dynamics simulation of tungsten. At the end of the simulation when I am seeing the total free energy (E) in the OSZICAR file, it is fluctuating and decreasing contineously. Since the total energy is not constant, is my simulation correct or wrong?
My INCAR file lokks like:
---------------------
ALGO = Fast
ENMAX=478.93 eV
MAXMIX = -45
IBRION = 0
NSW = 200
NBLOCK = 1
POTIM = 3.0
TEBEG = 1200
ISYM = 0
SMASS = -3
LREAL =Auto
LCHARG =.TRUE.
NELM = 40
NELMIN = 4
NELMDL = 0
PREC = LOW
ISTART = 0
ICHARG = 0
ISMEAR = 1
SIGMA=0.1
--------------------
Any kind of help is highly appreciable.
regards
Utpal
Last edited by utpalsarkar on Tue Apr 28, 2009 7:50 am, edited 1 time in total.

panda

How do I confirm that the Molecular dynamics run right or wrong?

#2 Post by panda » Fri May 01, 2009 7:11 pm

Try PREC = MED or HIGH, also you could just set the energy cutoff at a flat 500, you could play around with IBRION = 2 or 3 because your geometry might not be relaxed. After obtaining the best structure you could try using a Nose-Hover Thermostat to see if this makes any difference. Good Luck!
Last edited by panda on Fri May 01, 2009 7:11 pm, edited 1 time in total.

panda

How do I confirm that the Molecular dynamics run right or wrong?

#3 Post by panda » Fri May 01, 2009 7:13 pm

Try PREC = MED or HIGH, also you could just set the energy cutoff at a flat 500, you could play around with IBRION = 2 or 3 because your geometry might not be relaxed. Also, LREAL = Auto might need to be set to LREAL = .FALSE. After obtaining the best structure you could try using a Nose-Hover Thermostat to see if this makes any difference. Reading through the user's manual will help address specific questions for your research. Good Luck!
Last edited by panda on Fri May 01, 2009 7:13 pm, edited 1 time in total.

Post Reply