next up previous contents
Next: 10. Examples Up: 2. Detailed description of Previous: 8. Analysis, Properties and   Contents

Subsections


9. Utility Programs


1. afminput

This program helps to determine an input for clmcopy, which is a program which copies spin-up densities of atom 1 to spin-down densities of atom 2 and vice versa in an anti-ferromagnetic arrangement. It uses case.clmup and case.clmdn to find out how and which atoms (and Fourier coefficients) must be interchanged. Therefore "exact" antiferromagnetic densities are necessary, generated either by using the appropriate starting densities from dstart (after you put AFM configuration in case.inst) or after a very well converged SCF calculation. Please see the comments in sect. 4.5.4 on how to proceed in detail for AFM calculations.

There is no other input file, but you may be ask to decide, whether two densities are "equal" or not.


1. Execution

The program afminput is executed by invoking the command:
afminput afminput.def or x afminput

2. Dimensioning parameters

The following parameters are used:

NCOM number of LM components in the density (in param.inc)
NDIG number of significant digits for "equivalency" of density components
NDIG1 as above for Fouriercoefficients
EMIN exponent of a density component, which will be considered as "zero" in the program


2. clmcopy

This program generates the spin-dn density (case.clmdn) from a given spin-up density (case.clmup) according to rules in case.inclmcopy (generated earlier by afminput) for an AFM calculation. Please see the comments in sect. 4.5.4 on how to proceed in detail for AFM calculations.


1. Execution

The program clmcopy is executed by invoking the command:
clmcopy clmcopy.def or x clmcopy

2. Dimensioning parameters

The following parameters are used in param.inc:

NCOM number of LM components in the density
NRAD number of radial mesh points
NSYM number of symmetryoperations


3. Input

An example is given below:


---------------- top of file: case.inclmcopy -----------------------
   3                   NUMBER of ATOMS to CHANGE
   1   2               INTERCHANGE these ATOMS
   0                   NUMBER of LM to CHANGE SIGN
                       LM list
   3   3               INTERCHANGE these ATOMS
   4                   NUMBER of LM to CHANGE SIGN
 2 2 4 2 6 2 6 6       LM list
   4   4               INTERCHANGE these ATOMS
   6                   NUMBER of LM to CHANGE SIGN
 1 0 3 0 3 2 5 0 5 2 5 4    LM list
   3                   NUMBER of K-VECTOR-types to change sign
 x   y   z
odd 0   0   -1.0
odd 0   odd -1.0
evenodd 0   -1.0

Interpretive comments on this file are as follows:

line 1:
free format
NATOM   Number of atoms for which rules for copying the density will be defined
line 2:
free format
N1, N2   Interchange spin-up and dn densities of atoms N1 and N2
line 3:
free format
NLM   Number of LM values, for which you have to change the sign when swapping up and dn-densities
line 4:
free format
L,M   NLM pairs of L and M, for which you have to change the sign when swapping up and dn-densities

Lines 2-4 have to be repeated NATOM times.

line 5:
free format
NKK   number of K-vector types, for which you have to change the sign when swapping up and dn-densities
line 6:
free format
    empty line (just comments)
line 7:
free format
x,y,z,factor    
  x,y,z character of the K-vector component. Can be even, odd or 0.
  factor factor with which one multiplies the swaped up and dn Fouriercomponents (usually -1.)

Line 7 has to be repeated NKK times.

3. reformat

To produce a surface plot of the electron density using rhoplot_lapw (which is an interface to gnuplot), data from the file case.rho created by lapw5 must be converted using reformat

The sources of the program reformat.c are supplied in SRC_reformat.

4. hex2rhomb and rhomb_in5

hex2rhomb interactively converts the positions of an atom from hexagonal to rhombohedral coordinates (needed in case.struct).

rhomb_in5 interactively helps to generate input case.in5 for rhombohedral systems. It defines a plane as needed in the input file when you specify 3 atoms of that plane.

The sources of these programs are supplied in SRC_trig.


5. eosfit

Small program to calculate the Equation of States (EOS; Equilibrium volume $V_0$, Bulk modulus $B_0$ and it's derivative $B_0'$, see Murnaghan 1944). It relies on the file case.analysis created from w2web using Total Energy and Volume.

The sources are supplied in SRC_eosfit.


6. spacegroup

This program was contributed by:

\framebox {
\parbox[c]{12cm}{
Vaclav Petricek\\
Institute of Physics\\
A...
...ilinglist. If necessary, we will communicate
the problem to the authors.}
}
}


Interactive program to generate equivalent positions for a given spacegroup and lattice. The program is also used internally from w2web to generate positions when selecting spacegroups in the StructGen $^{\mbox{\textsc{TM}}}$.


7. analyse

This program was contributed by:

\framebox {
\parbox[c]{12cm}{
Hartmut Enkisch \\
Institute of Physics E1b\...
...ilinglist. If necessary, we will communicate
the problem to the authors.}
}
}


analyse reads data from case.scf and writes them into analyse.dat. For each iteration it writes one line containing the parameters specified in the command line. analyse is meant to work in conjunction with the scfmonitor script (see 5.2.5).


8. StructGen $^{\mbox{\textsc{TM}}}$ of w2web

The new StructGen $^{\mbox{\textsc{TM}}}$ helps to generate the master input file case.struct. It has the following additional features:


9. Visualization


1. BALSAC

balsac (Build and Analyze Lattices, Surfaces and Clusters) was written by Klaus Hermann (Fritz-Haber Institut, Berlin). It provides high quality postscript files. In SRC_balsac-utils we provide the following interface programs to convert from WIEN2k to balsac:

For an example see figure 3.1 For scientific questions concerning BALSAC please contact Klaus Hermann at hermann@FHI-Berlin.MPG.DE

Balsac is available from:

Garching Innovation GmbH, Mrs. M. Pasecky Koeniginstr. 19, D-80539 Munich, Germany
Tel.: +49 89 288279, Fax.: +49 89 21081593
e-mail: 100627.1572@compuserve.com


2. XCrysDen

XCrysDen (Kokalj 1999) is a render and analysis package. It has the following features:

Figure 9.1: 3D electron density in TiC generated with XCrysDen
\begin{figure}
\begin{center}
\leavevmode
\epsfig {figure=figs/tic/xcrysden, height=50mm} \end{center}\end{figure}

XCrysDen is available from:

Tone Kokalj
Jozef Stefan Institute, Dept. of Physical and Organic Chemistry
Jamova 39, SI-1000 Ljubljana, Slovenia
Tel.: +386 61 177 3520, Fax: +386 61 177 3811
Tone.Kokalj@ijs.si


next up previous contents
Next: 10. Examples Up: 2. Detailed description of Previous: 8. Analysis, Properties and   Contents
Dieter Kvasnicka
2001-12-05