Author Topic: Negative numbers in density of states  (Read 6419 times)

0 Members and 1 Guest are viewing this topic.

Offline Roc

  • QuantumATK Guru
  • ****
  • Posts: 104
  • Reputation: 0
    • View Profile
Negative numbers in density of states
« on: January 7, 2009, 02:41 »
Dear everyone,

Error happens in ATK caculation, as follows,


Traceback (most recent call last):
  File "AumolecularAuscript.py", line 291, in ?
    mulliken_population = calculateMullikenPopulation(self_consistent_calculation = scf)
NLTypeError: The population must contain non-negative numbers only.[/b]

that‘s to say, some negetive numbers happen in DOS, the .dat document is that,

# Density Of States
# -----------------------------------------------------------------------------
# Energy (eV)   DOS (1/eV)
        -4.00       0.0001
        -3.90       0.0002
        -3.80       0.0003
          ................
         3.60      -0.0012
         3.70      -0.0012
         3.80      -0.0016
         3.90      -0.0017
         4.00      -0.0016

how to deal with it?
Thanks!

Moderator edit: Changed subject to improve searching
« Last Edit: January 8, 2009, 10:29 by Anders Blom »

Offline Nordland

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 812
  • Reputation: 18
    • View Profile
Re: Another error
« Reply #1 on: January 7, 2009, 14:41 »
This one is a bit nasty - I guess it should not be able to happen, but in some way the convergence has gone haywire, but still managed to find a fix point.
You need to improve the quality of your twoprobe calculation by getting a more physical correct convergence.