[DFTB-Plus-User] work function calculation for a metallic slab
Ravi Kumar
ravikumar381709 at gmail.com
Tue Jun 17 11:46:02 CEST 2025
Dear all,
I am trying to calculate the work function of a tungsten slab using an
*ElectrostaticPotential* calculation (an example input file is given
below). However, the XY-planer averaged potential does not appear to be
correct when compared with VASP results.
I am not sure whether the issue lies in the calculation setup or in the
postprocessing script. I would greatly appreciate any guidance to help
identify and resolve the problem.
Any help or suggestions would be highly appreciated.
Geometry = VASPFormat {
W
1.00000000000000
9.4950000000000010 0.0000000000000000 0.0000000000000000
0.0000000000000000 9.4950000000000010 0.0000000000000000
0.0000000000000000 0.0000000000000000 23.1649999999999991
W
27
Direct
0.0000000000000000 0.0000000000000000 0.4316857327865335
0.3333333333333357 0.0000000000000000 0.4316857327865335
0.6666666666666643 0.0000000000000000 0.4316857327865335
0.0000000000000000 0.3333333333333357 0.4316857327865335
0.3333333333333357 0.3333333333333357 0.4316857327865335
0.6666666666666643 0.3333333333333357 0.4316857327865335
0.0000000000000000 0.6666666666666643 0.4316857327865335
0.3333333333333357 0.6666666666666643 0.4316857327865335
0.6666666666666643 0.6666666666666643 0.4316857327865335
0.1666666666666643 0.1666666666666643 0.5000000000000000
0.5000000000000000 0.1666666666666643 0.5000000000000000
0.8333333333333357 0.1666666666666643 0.5000000000000000
0.1666666666666643 0.5000000000000000 0.5000000000000000
0.5000000000000000 0.5000000000000000 0.5000000000000000
0.8333333333333357 0.5000000000000000 0.5000000000000000
0.1666666666666643 0.8333333333333357 0.5000000000000000
0.5000000000000000 0.8333333333333357 0.5000000000000000
0.8333333333333357 0.8333333333333357 0.5000000000000000
0.0000000000000000 0.0000000000000000 0.5683142672134665
0.3333333333333357 0.0000000000000000 0.5683142672134665
0.6666666666666643 0.0000000000000000 0.5683142672134665
0.0000000000000000 0.3333333333333357 0.5683142672134665
0.3333333333333357 0.3333333333333357 0.5683142672134665
0.6666666666666643 0.3333333333333357 0.5683142672134665
0.0000000000000000 0.6666666666666643 0.5683142672134665
0.3333333333333357 0.6666666666666643 0.5683142672134665
0.6666666666666643 0.6666666666666643 0.5683142672134665
}
Hamiltonian = DFTB {
SCC = Yes
SCCTolerance = 1e-6
MaxSCCIterations = 1000
SlaterKosterFiles = Type2FileNames {
Prefix = "/usr/local/slakos/Voss/rPBE/"
Separator = "-"
Suffix = ".skf"
LowerCaseTypeName = No
}
MaxAngularMomentum = {
W = "d"
}
KPointsAndWeights = SupercellFolding {
12 0 0
0 12 0
0 0 1
0.0 0.0 0.0
}
}
Analysis = {
AtomResolvedEnergies = Yes
ElectrostaticPotential = {
Grid [F] = {
Spacing = {1.0 1.0 1.0}
Origin = {0 0 0}
GridPoints = {100 100 100}
}
}
}
Options = {
WriteAutotestTag = Yes
}
ParserOptions = {
ParserVersion = 5
}
Thanks and best regards,
Ravi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.zfn.uni-bremen.de/pipermail/dftb-plus-user/attachments/20250617/3f3c8d2b/attachment.htm>
More information about the DFTB-Plus-User
mailing list