Example: NaCl: Bulk Crystal¶
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
Radial 2000
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
Radial 2000
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