[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