mirror of https://github.com/abinit/abinit.git
81 lines
1.9 KiB
Plaintext
81 lines
1.9 KiB
Plaintext
# Crystalline silicon
|
|
# Same physical test as in tutorial t31.in, but here
|
|
# test the DOS using the tetrahedron method
|
|
#
|
|
|
|
ndtset 2
|
|
|
|
#Dataset 1 : SCF calculation
|
|
ngkpt1 2 2 2 # Use a rather loose grid for SCF calculation
|
|
toldfe1 1.0d-10
|
|
prtden1 1
|
|
|
|
#Dataset 2 : DOS calculation
|
|
iscf2 -3
|
|
getden2 1
|
|
getwfk2 1
|
|
prtdos2 2
|
|
dosdeltae 0.00005 # This is also the default value for prtdos2=2,
|
|
# but the presence of dosdeltae is there for checking purposes
|
|
ngkpt2 4 4 4 # This is not sufficient for a nice DOS.
|
|
# One should use ngkpt2 8 8 8
|
|
tolwfr2 1.0d-16
|
|
shiftk2 0.0 0.0 0.0 # Use a non-shifted k point grid, giving a much
|
|
0.0 0.5 0.5 # better-looking DOS graph than the shifted one,
|
|
0.5 0.0 0.5 # because tetrahedra with all four same energies
|
|
0.5 0.5 0.0 # are avoided
|
|
|
|
|
|
#Definition of the unit cell
|
|
acell 3*10.18
|
|
rprim 0.0 0.5 0.5
|
|
0.5 0.0 0.5
|
|
0.5 0.5 0.0
|
|
|
|
#Definition of the atom types
|
|
ntypat 1
|
|
znucl 14
|
|
|
|
#Definition of the atoms
|
|
natom 2
|
|
typat 1 1
|
|
xred 0.0 0.0 0.0
|
|
1/4 1/4 1/4
|
|
|
|
#Definition of the planewave basis set
|
|
ecut 8.0
|
|
|
|
#Definition of the k-point grid
|
|
kptopt 1
|
|
nshiftk 4
|
|
shiftk 0.5 0.5 0.5
|
|
0.5 0.0 0.0
|
|
0.0 0.5 0.0
|
|
0.0 0.0 0.5
|
|
|
|
#Definition of the SCF procedure
|
|
nstep 10
|
|
diemac 12.0
|
|
|
|
|
|
pp_dirpath "$ABI_PSPDIR"
|
|
pseudos "PseudosTM_pwteter/14si.pspnc"
|
|
|
|
#%%<BEGIN TEST_INFO>
|
|
#%% [setup]
|
|
#%% executable = abinit
|
|
#%% [files]
|
|
#%% files_to_test =
|
|
#%% t46.abo, tolnlines = 0, tolabs = 0.000e+00, tolrel = 0.000e+00;
|
|
#%% t46o_DS2_DOS, tolnlines = 0, tolabs = 0.000e+00, tolrel = 0.000e+00
|
|
#%% [paral_info]
|
|
#%% max_nprocs = 2
|
|
#%% [extra_info]
|
|
#%% authors = Unknown
|
|
#%% keywords =
|
|
#%% description =
|
|
#%% 2 Silicon atoms, in a diamond structure.
|
|
#%% Test the computation of the DOS, using the linear tetrahedron method.
|
|
#%% topics = ElecDOS
|
|
#%%<END TEST_INFO>
|