Vasp transition state theory tools

Moderator: moderators

flyhigh
Posts: 3
Joined: Wed May 25, 2016 5:43 am

Hello,

I calculated d-band center of Pt surface with this, and I got the results below.
Only Emax is different in the case (1) and the case (2).
I do not understand the meaning of the last message,"Can't take sqrt of -415.985 at /bin/vtstscripts/dosanalyze.pl line 292."
Why did it make problems to increase Emax?
Logically, it does not seem to have problems.
In other words, how can I decide the range of Emin and Emax?
Average Energy (band center) is affected by the range of Emin and Emax too much.

Thanks.

case (1)
\$ dosanalyze.pl e=-7,0 d

Integrate from -7eV to 0eV.
WARNING: For DOS0, ALL bands are analyzed
WARNING: band center has shifted by more than 10% compared to full integrated range

Total States: 69.43606769
Average Energy (band center): -2.25169355987311
Standard Deviation: 2.48325612734491

case (2)
\$ dosanalyze.pl e=-7,1 d

Integrate from -7eV to 1eV.
WARNING: For DOS0, ALL bands are analyzed
Can't take sqrt of -415.985 at /bin/vtstscripts/dosanalyze.pl line 292.
graeme
Posts: 1999
Joined: Tue Apr 26, 2005 4:25 am
Contact:

If you can post the directory with the data, I'll try to reproduce the problem and find a solution.
flyhigh
Posts: 3
Joined: Wed May 25, 2016 5:43 am

[quote="graeme"]If you can post the directory with the data, I'll try to reproduce the problem and find a solution.[/quote]

I checked the script "dosanalyze.pl ", and I found the error message is related to the variance.

290 \$eneavg = \$dstsum/\$dossum;
291 \$enevar = (\$dssqsum*\$dossum - \$dstsum*\$dstsum)/(\$dossum*\$dossum);
292 \$enestd = sqrt(\$enevar);

In case I used ISPIN=2 in INCAR file, can this problem be caused by this?
graeme