# Example: NaCl: Bulk Crystal¶

Download NaCl.run

A bulk crystal computation for Sodium Chloride (common salt), with a subsequent DOS analysis, using a Restart facility to use the results from a preceding calculation.

The BAND input follows slightly different conventions from the ADF input, for historical reasons.

The COMMENT keyword allows users to provide some information about the run which may be of use later. Usually a brief summary of the run is given here.

Numerical integration precision is controlled with the key Accuracy. The accuracy for integrals over the Brillouin Zone is set by the KSpace key. The latter should, generally, take as value an odd number (3, 5...) to invoke the more accurate quadratic tetrahedron integration procedure. For even values it will revert to the linear tetrahedron method, which is almost always inferior in accuracy.

The key Angstroms specifies that geometric data, such as lattice constants are in Angstrom units.

Since there are 3 data records in the Lattice block, the calculation will assume 3-dimensional periodicity, with lattice vectors as indicated. Note that lattice vectors are undefined up to linear combinations among themselves. Internally, the program will recombine the input vectors so as to minimize the size of the actually used vectors.

The input line Coordinates Fractional means that atomic positions are input as coefficients in terms of the lattice vectors, rather than as absolute (Cartesian) coordinate values.

For each of the atoms in the calculations, Na and Cl here, there must be data blocks to specify various items. First, their positions in the crystal unit cell (key Atoms). Second, the single isolated atom computation that will serve as start-up (Dirac). Third, any Slater-type orbital basis functions (BasisFunctions) for that atom. Fourth, the fit functions (FitFunctions) for the calculation of the Coulomb potential and. The third item (BasisFunctions) is optional and not present in this example.

The basis set, numerical atomic orbitals that are computed by the Herman-Skillman type subprogram DIRAC, and the fit sets that one has to insert into the BAND input files can be taken from the database you find in the directory $ADFHOME/atomicdata/Band. In contrast to the shown example, a more comfortable and automatized way of specifying the basis and fit sets is the keyblock BasisDefaults. After the first calculation, a single point calculation of bulk NaCl, is finished the standard result file from BAND, the RUNKF file, is saved under the name ‘NaCl.runkf’. This file will be used in the follow-up calculation of Density-of-States properties. $ADFBIN/band << eor
Title Title NaCl (from neutral atoms)

Comment
Technical
Hybrid K space integration (3D)
Low real space integration accuracy
Fractional coordinates
Lengths in Angstrom
Parameters Dirac procedure
Features
Lattice   : 3D
Unit cell : 2 atoms
Basis     : NO w/ core
Options   : Save restart file
End

Accuracy 3.5
Kspace 3

Units
Length Angstrom
End

Lattice
0    2.75   2.75
2.75  0     2.75
2.75  2.75   0
End

Coordinates Fractional
Atoms
Na 0.0 0.0 0.0
Cl 0.5 0.5 0.5
End

AtomType Na
Dirac   Na
4 1
RMin    1E-4
RMax    60
VALENCE
1 0
2 0
2 1
3 0 1.0
SubEnd

FitFunctions
1 0 18.9
2 0 30.3
2 0 15.5
3 0 14.9
3 0  8.9
4 0  7.8
4 0  5.1
4 0  3.3
5 0  2.8
5 0  1.9
5 0  1.3
2 1 14.3
3 1  9.9
4 1  6.7
4 1  3.4
5 1  2.4
5 1  1.3
3 2 10.5
4 2  5.4
5 2  3.0
5 2  1.3
4 3  5.8
5 3  1.7
5 4  2.0
SubEnd
End

AtomType Cl
Dirac  Cl
5 3
VALENCE
1 0
2 0
2 1
3 0
3 1  5.0
SubEnd

FitFunctions
1 0 29.1
2 0 49.5
2 0 26.1
3 0 25.8
3 0 15.8
4 0 14.2
4 0  9.4
4 0  6.2
5 0  5.4
5 0  3.8
5 0  2.6

2 1 21.2
3 1 16.5
4 1 12.4
4 1  6.8
5 1  5.1
5 1  3.1

3 2 16.6
4 2  9.4
5 2  5.5
5 2  2.6

4 3  8.7
5 3  3.3

5 4  4.0
SubEnd
End

End Input
eor

mv RUNKF NaCl.runkf

rm Points


The next run has largely the same input and provides a restart of the previous run.

The key DOS in the block Restart tells the program to pick up the indicated file as restart file and to use it for DOS analysis purposes.

The DOS key block details the energy grid (and range) and the file to write the data to. The optional keys GROSSPOPULATIONS and OverlapPopulations invoke the computation of, respectively, gross populations and overlap populations (i.e. for each of these the density-of-states values in the user-defined energy grid).

\$ADFBIN/band << eor
Title Title NaCl (from neutral atoms)   DOS analysis (restart)

Comment
Technical
Hybrid K space integration (3D)
Low real space integration accuracy
Fractional coordinates
Lengths in Angstrom
Parameters Dirac procedure
Features
Lattice   : 3D
Unit cell : 2 atoms
Basis     : NO w/ core
Options   : Use restart file for DOS
Analysis: DOS, PDOS, COOP
End

Restart
File NaCl.runkf
DOS
End

Accuracy 3.5
Kspace 3

SCF
Iterations 15
End

Units
Length Angstrom
End

Lattice
0    2.75   2.75
2.75  0     2.75
2.75  2.75   0
End

Coordinates Fractional
Atoms
Na 0.0 0.0 0.0
Cl 0.5 0.5 0.5
End

DOS
File NaCl.dos
Energies 1000
Min -0.5
Max  0.5
End

GROSSPOPULATIONS
FRAG 1
FRAG 2
SUM
1 0
2 0
ENDSUM
End

OVERLAPPOPULATIONS
LEFT
FRAG 1
RIGHT
FRAG 2
LEFT
1 0
1 1
RIGHT
2 0
2 1
End

AtomType Na
Dirac   Na
4 1
RMin    1E-4
RMax    60
VALENCE
1 0
2 0
2 1
3 0 1.0
SubEnd

FitFunctions
1 0 18.9
2 0 30.3
2 0 15.5
3 0 14.9
3 0  8.9
4 0  7.8
4 0  5.1
4 0  3.3
5 0  2.8
5 0  1.9
5 0  1.3
2 1 14.3
3 1  9.9
4 1  6.7
4 1  3.4
5 1  2.4
5 1  1.3
3 2 10.5
4 2  5.4
5 2  3.0
5 2  1.3
4 3  5.8
5 3  1.7
5 4  2.0
SubEnd
End

AtomType Cl
Dirac  Cl
5 3
VALENCE
1 0
2 0
2 1
3 0
3 1  5.0
SubEnd

FitFunctions
1 0 29.1
2 0 49.5
2 0 26.1
3 0 25.8
3 0 15.8
4 0 14.2
4 0  9.4
4 0  6.2
5 0  5.4
5 0  3.8
5 0  2.6

2 1 21.2
3 1 16.5
4 1 12.4
4 1  6.8
5 1  5.1
5 1  3.1

3 2 16.6
4 2  9.4
5 2  5.5
5 2  2.6

4 3  8.7
5 3  3.3

5 4  4.0
SubEnd
End

End Input
eor


Finally, we copy the contents of the DOS result file to standard output

echo Contents of DOS file
cat NaCl.dos