reproducing data in restarting CPMD simulation

From NWChem

Viewed 1224 times, With a total of 0 Posts
Jump to: navigation, search

Just Got Here
Threads 1
Posts 1
Hi, all.
I have an issue with reproducing data (temperture) in restarting CPMD simulation. I ran the CPMD simulation using PSPW(psedopotential planewave). It is NVE simulation, and I didn’t use any thermostat method. Because I set the number of time steps too less, the simulation was completed before the temperature of the system was equilibrated. So, I restarted program to resume the simulation. To make sure the restart had been correctly done, I checked the .xyz and .energy files- the ones at the last time step in the original simulation and the others at the first time step of restarted simulation. xyz files are exactly same, and kinetic/ potential energies look pretty much same. But there is big discrepancy in temperature: 1230K vs 150K.

I have two questions regarding this restarted simulation. First, what should I do to reproduce the temperature correctly. Is there any command to be added in input file for that? Second, To my knowledge, the temperature is calculated based on the velocities of atoms. In my restarted-simulation, atom velocities were correctly read. So, I’m thinking the temperature also should be same. Would you briefly explain why I got the different temperature in a restarted simulation. As a student, any comments on that would be greatly helpful and appreciated.

I attach the input file for your reference.

geometry units an noautosym noautoz print
 # geometry info  

ngrid 32 32 32
cutoff 29
ewald_ncut 8
xc xperdew91

VECTORS input.movecs
DENSITY total edensity

fake_mass 400.0
time_step 5
xc xperdew91
loop 1 2000
ion_motion_filename ionmotion
emotion_filename energy

set pspw_dplot:iteration_list 1000
task pspw energy
task pspw Car-Parrinello

Forum >> NWChem's corner >> General Topics

Who's here now Members 0 Guests 1 Bots/Crawler 0

AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC