Skip to content

Commit

Permalink
Update chapters in Basic usage
Browse files Browse the repository at this point in the history
Adapt to changed output, set current parser version in input, make
run.sh scripts complete.
  • Loading branch information
aradi committed Nov 30, 2018
1 parent 4faea9e commit ab8d67d
Show file tree
Hide file tree
Showing 9 changed files with 143 additions and 118 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -14,8 +14,8 @@ Geometry = GenFormat {
}

Hamiltonian = DFTB {
SCC = Yes
SCCTolerance = 1e-5
Scc = Yes
SccTolerance = 1e-5
SlaterKosterFiles = Type2FileNames {
Prefix = "../../../slakos/mio-ext/"
Separator = "-"
Expand All @@ -25,9 +25,6 @@ Hamiltonian = DFTB {
Ti = "d"
O = "p"
}
Filling = Fermi {
Temperature [Kelvin] = 0.0
}
KPointsAndWeights = SupercellFolding {
4 0 0
0 4 0
Expand All @@ -52,6 +49,6 @@ Analysis {
}

ParserOptions {
ParserVersion = 5
ParserVersion = 6
}

7 changes: 7 additions & 0 deletions docs/_archives/recipes/basics/bandstruct/1_density/run.sh
Original file line number Diff line number Diff line change
@@ -1,3 +1,10 @@
#!/bin/bash
set -e
dftb+ > output
dp_dos band.out dos_total.dat >> output
dp_dos -w dos_ti.1.out dos_ti.s.dat >> output
dp_dos -w dos_ti.2.out dos_ti.p.dat >> output
dp_dos -w dos_ti.3.out dos_ti.d.dat >> output
dp_dos -w dos_o.1.out dos_o.s.dat >> output
dp_dos -w dos_o.2.out dos_o.p.dat >> output

9 changes: 3 additions & 6 deletions docs/_archives/recipes/basics/bandstruct/2_bands/dftb_in.hsd
Original file line number Diff line number Diff line change
Expand Up @@ -14,9 +14,9 @@ Geometry = GenFormat {
}

Hamiltonian = DFTB {
SCC = Yes
Scc = Yes
ReadInitialCharges = Yes
MaxSCCIterations = 1
MaxSccIterations = 1
SlaterKosterFiles = Type2FileNames {
Prefix = "../../../slakos/mio-ext/"
Separator = "-"
Expand All @@ -26,9 +26,6 @@ Hamiltonian = DFTB {
Ti = "d"
O = "p"
}
Filling = Fermi {
Temperature [Kelvin] = 0.0
}
KPointsAndWeights = Klines {
1 0.5 0.5 -0.5 # Z
20 0.0 0.0 0.0 # G
Expand All @@ -38,5 +35,5 @@ Hamiltonian = DFTB {
}

ParserOptions {
ParserVersion = 5
ParserVersion = 6
}
2 changes: 2 additions & 0 deletions docs/_archives/recipes/basics/bandstruct/2_bands/run.sh
Original file line number Diff line number Diff line change
Expand Up @@ -2,3 +2,5 @@
set -e
cp ../1_density/charges.bin .
dftb+ > output
dp_bands band.out band >> output

3 changes: 0 additions & 3 deletions docs/_archives/recipes/basics/firstcalc/dftb_in.hsd
Original file line number Diff line number Diff line change
Expand Up @@ -26,9 +26,6 @@ Hamiltonian = DFTB {
O = "p"
H = "s"
}
Filling = Fermi {
Temperature [Kelvin] = 0.0
}
}

Options {}
Expand Down
38 changes: 18 additions & 20 deletions docs/basics/bandstruct.rst
Original file line number Diff line number Diff line change
Expand Up @@ -17,20 +17,20 @@ you will need the Slater-Koster sets `mio` and `tiorg`. The sample input files
assume that the necessary Slater-Koster files have been copied into a
subdirectory `mio-ext`.

The examples here are based on DFTB+ version 17.1, the input/output
files in other versions may slightly differ from those shown here.


Introduction
============

The calculation of the band structure for a periodic system consists of two
steps. First for self-consistent (SCC) calculations, the charges in the system
must be calculated using a converged k-point sampling. Then, keeping the
obtained self-consistent charges fixed, the one-electron levels must be
calculated for k-points chosen along the specific lines in k-space of the chosen
band structure. These are usually between high symmetry points in the Brillouin
zone of that unit cell.
steps.

* First for self-consistent (SCC) calculations, the charges in the system must
be calculated using a converged k-point sampling.

* Then, keeping the obtained self-consistent charges fixed, the one-electron
levels must be calculated for k-points chosen along the specific lines in
k-space of the chosen band structure. These are usually between high symmetry
points in the Brillouin zone of that unit cell.


Creating the proper input charges
Expand Down Expand Up @@ -84,9 +84,6 @@ good k-point sampling. A sample `dftb_in.hsd` input looks like::
Ti = "d"
O = "p"
}
Filling = Fermi {
Temperature [Kelvin] = 0.0
}
KPointsAndWeights = SupercellFolding {
4 0 0
0 4 0
Expand All @@ -111,7 +108,7 @@ good k-point sampling. A sample `dftb_in.hsd` input looks like::
}
ParserOptions {
ParserVersion = 5
ParserVersion = 6
}

In the input above, the coordinates have been specified in relative (fractional)
Expand All @@ -123,7 +120,7 @@ line of the geometry specification::
6 F
:

The k-points are been generated automatically using the ``SupercellFolding``
The k-points are generated automatically using the ``SupercellFolding``
method, which enables among others the generation of Monkhorst-Pack schemes. In
the current example, a k-point set equivalent to the Monkhorst-Pack scheme
:math:`4 \times 4 \times 4` has been chosen (For details how to specify the
Expand All @@ -141,7 +138,7 @@ an accuracy in the range of 1e-3 eV for the total energy. Also, by specifying a
smaller SCC tolerance than the chosen one (1e-5) you can check that converging
the charges more precisely does not significantly decrease the total energy. We
note in passing that these settings provide well converged results for the total
energy in the current example, by in principal may not provide converged values
energy in the current example, but in principal may not provide converged values
for other properties. One should, in principal, test the convergence of any
evaluated properties with respect to the calculation parameters.

Expand Down Expand Up @@ -259,8 +256,9 @@ to be changed slightly::

# ...

Note: only the relevant part of the input are shown, here. See the
`Introduction`_ how to obtain the archive with the full input.
Note: only the relevant parts of the input are shown, here. See the
:ref:`sec-introduction` section on how to obtain the archive with the full
input.

The input is (must be) almost the same as in the previous case, with only a few
adaptions:
Expand Down Expand Up @@ -303,9 +301,9 @@ adaptions:
point only, as demonstrated above for the Z-point.

Running DFTB+ with the input above, the eigenlevel spectrum is calculated at
the required k-points. The results are written to the file `detailed.out` and
in more readable format to `band.out`. You can use the script `dp_bands` from
the `dptools` package to convert this file into NXY format. By issuing::
the required k-points. The results are written to the file `band.out`. You can
use the script `dp_bands` from the `dptools` package to convert this file into
XNY format. By issuing::

dp_bands band.out band

Expand Down
Loading

0 comments on commit ab8d67d

Please sign in to comment.