UT theoretical chemistry code forum

New scripts under vtstscripts are:

1) nebef.pl : Gets the energy and force for each images in a neb chain
2) nebsetup.pl : Sets up a linear neb chain from two POSCARs
3) make_mode.pl : (Under construction). Makes a MODECAR for the lanczos method used by VASP.

Scripts made obsolete by this:

1) intPOS
2) nebforce.pl
3) nebcontrun.pl
4) pclean (alias)

As with all new scripts, they should be used with extreme caution to begin with and ALL problems should be reported back asap.

Andri

The following scripts were modified:

neb_setup.pl -> nebsetup.pl
make_mode.pl -> modemake.pl
chenergy -> /dev/null

nebresults.pl: changed the output format, and called the internal scripts via $Bin/intscript.pl.

I added three scripts: avgCHGCAR.pl, avgPARCHG.pl and diffCHGCAR.pl, which either average or subtract CHGCAR/PARCHG files.

The two types of files are very similar (CHGCAR vs PARCHG) except in how many entries they have in each line for the charge density on each FFT grid point.

For them to work, the ENAUG of each file must be the same. This is typically not a problem, unless, for example, in one system you have Si and P but in the second only Si. Since P has a large ENAUG, the one with only Si will have a different one and this affects the FFT grid, if I recall correctly.

These are crude, but may serve as a good starting point for further work.

Blas

I've change these script names, using the prefix chg:

avgCHGCAR.pl -> chgavg.pl
difffCHGCAR.pl -> chgdiff.pl
avgPARCHG.pl -> chgparavg.pl

I've also set the perl binary to be the general call which works on all platforms.

Thanks again for these scripts, I think they will be very useful for visualizing charge density differences. The output from these scripts should work with vaspview.