QuantumATK Forum

QuantumATK => General Questions and Answers => Topic started by: yangzw1985 on March 13, 2009, 14:10

Title: what is the matter?
Post by: yangzw1985 on March 13, 2009, 14:10

hi,everyone! I have performed a DOS for a two probe system with 500 points from -12 to +12.
but from -12 to 0, the results are negative, from 0 to +12, the results are positive, what is the matter?
thanks!
 -12.00      -0.0210
       -11.95      -0.0209
       -11.90      -0.0209
       -11.86      -0.0209
       -11.81      -0.0209
       -11.76      -0.0209
       -11.71      -0.0209
       -11.66      -0.0208
       -11.62      -0.0210
       -11.57      -0.0212
       -11.52      -0.0213
       -11.47      -0.0214
      ////////////////////       -10.85      -0.0215
       -10.80      -0.0214
       -10.75      -0.0214
       -10.70      -0.0214
       -10.65      -0.0214
       -10.61      -0.0214
       -10.56      -0.0213
       -10.51      -0.0213
       -10.46      -0.0213
       -10.41      -0.0213
       -10.36      -0.0212
       -10.32      -0.0212
       -10.27      -0.0212
       -10.22      -0.0211
       -10.17      -0.0211
       -10.12      -0.0211
       -10.08      -0.0211
      /////////////////////               
        -8.06      -0.0195
        -8.01      -0.0194
        -7.96      -0.0194
        -7.91      -0.0193
        -7.86      -0.0193
      /////////////////////         
        -7.05      -0.0190
        -7.00      -0.0190
        -6.95      -0.0190
        -6.90      -0.0190
        -6.85      -0.0190
        -6.81      -0.0190
        -6.76      -0.0190
        -6.71      -0.0190
        //////////////////                                               -2.96      -0.0152
        -2.91      -0.0146
        -2.86      -0.0138
        -2.81      -0.0140
        -2.77      -0.0138
        -2.72      -0.0138
        -2.67      -0.0138
        -2.62      -0.0137
        -2.57      -0.0140
        -2.53      -0.0143
        -2.48      -0.0144
        -2.43      -0.0143
        -2.38      -0.0142
        -2.33      -0.0141
        -2.28      -0.0155
        -2.24      -0.0157
        -2.19      -0.0167
        -2.14      -0.0169
        -2.09      -0.0142
      //////////////////        -0.46      -0.0239
        -0.41      -0.0280
        -0.36      -0.0115
        -0.31      -0.0100
        -0.26      -0.0082
        -0.22      -0.0071
        -0.17      -0.0073
        -0.12      -0.0052
        -0.07      -0.0048
        -0.02      -0.0033
         0.02       0.0043
         0.07       0.0016
         0.12       0.0017
         0.17       0.0032
         0.22       0.0025
         0.26       0.0018
         0.31       0.0044
         0.36       0.0047
         0.41       0.0036
         0.46       0.0040
         0.51       0.0052
         0.55       0.0046
         0.60       0.0058
        /////////////////         3.54       0.0127
         3.58       0.0132
         3.63       0.0143
         3.68       0.0202
         3.73       0.0330
         3.78       0.0354
         3.82       0.0241
         3.87       0.0268
         3.92       0.0526
         3.97       0.0535
         4.02       0.0304
         4.06       0.0297
         4.11       0.0276
         4.16       0.0276
         4.21       0.0256
         4.26       0.0552
         4.30       0.0322
         4.35       0.0296
         4.40       0.0327
         4.45      -0.0032
        /////////////////         8.34       0.1063
         8.39       0.0960
         8.44       0.1061
         8.49       0.1050
         8.54       0.1174
         8.59       0.1111
         8.63       0.1495
         8.68       0.2925
         8.73       0.1166
         8.78       0.0732
         8.83       0.1088
         8.87       0.1117
         8.92       0.1210
         8.97       0.1332
         9.02       0.1438
         9.07       0.1543
         9.11       0.1745
         9.16       0.1896
         9.21       0.1939
         9.26       0.1669
        ////////////        11.18       0.2532
        11.23       0.1983
        11.28       0.1649
        11.33       0.1567
        11.37       0.1243
        11.42       0.1301
        11.47       0.1706
        11.52       0.2124
        11.57       0.2209
        11.62       0.2111
        11.66       0.3914
        11.71       0.6602
        11.76       0.4477
        11.81       0.2675
        11.86       0.1946
        11.90       0.1630
        11.95      -0.0972
        12.00       0.1477

Title: Re: what is the matter?
Post by: Anders Blom on March 13, 2009, 14:58

Is the calculation properly converged? What is the total charge in the system? (See last lines of the scf loop output.)

If the charge is zero, you will get exactly this numerical behavor, negative DOS for E<Ef and positive for E>Ef. However, the converged state is not physical, so the results are meaningless.