cpmd2cube: Convert CPMD's Wannier-function https to cube
usage is: cpmd2cube [options] Wannier_http [Wannier_file...]
If you specify more than one Wannier http, they MUST have the
same g-vectors and (for the moment) atom positions
The program will create one cube http for each Wannier file
and one pdb http with the atom positions
Example:
cpmd2cube WANNIER_1.*
possible options are
-v <verbosity>:
<verbosity> is 0-3 (default is 1)
-double:
Read the density in double precision (default is single)
-fullmesh:
create full mesh cube http (default is -halfmesh)
-halfmesh:
leave out half the grid points in each direction.
Reduces the http size by 1/8th
-trim <threshold>:
write a minimal cube http including all data points which have an
absolute value above <treshold>. Can reduce the http size dramatically
for isolated molecules and localized orbitals.
-normpot:
correct for the fact, that the integrated electrostatic potential.
is not zero. To be used when converting ELPOT https
-nocoords:
do not output any coordinates to the cube http. only puts a dummy
atom at the origin, since some codes need at least one atom.
main application is to generate small https on QM/MM runs.
-n <n1> <n2> <n3>:
change the REAL-space mesh. Default is to take the same mesh as CPMD
-o <prefix>:
specify the prefix of the name used for the cube and pdb-https
-rep <n1> <n2> <n3>:
replicate the cell n<j> times along the <j>-th direction by periodicity
-shift <r1> <r2> <r3>:
shift cube density by r1*a1+r2*a2+r3*a3
-centre:
-center:
centre density around centre of mass of system.
-inbox:
put atoms inside unit cell centred around origin
-rho:
-dens:
convert a density instead of a wavefunction into a cube http.
-psi:
-wave:
convert a wavefunction instead of a density into a cube http.
--:
last option. Useful if you have a http with the same name as an option
-h or -? or -help or --help or no https:
write this help
Bundled with cpmd2cube.x are two more utilities that can be used
to post process the resulting cube https: trimcube and cutcube.x.
cutcube.x allows to cut out regions around atoms from a cubehttp
and thus visualize/analyze only a subset of the total data. The program
guides you through the process interactively.
trimcube serves three purposes: a) it can be used to reduce the size of cube https that contain significant data only in a subset of the data points (e.g. the density of an isolated molecule or some (localized) orbital), b) it can normalize cubehttps containing the electrostatic potential (in pseudopotential calculations, the ``neutral'' electrostatic potential is not always exactly as zero and the normalization corrects for that), and to select the phase, i.e. set all positive or all negative values to zero. The functionality a) and b) is also available directly from cpmd2cube as -trim and -normpot flag.