I happened to stumble upon something today that I thought was unusual.  The problem is illustrated in a simple example: if you take the following basic cubic Si input file and, say, the Si.pz-vbc.UPF pseudopotential from the QE website and do a test for convergence with respect to ectutwfc.<br>

<br>If we look at the 4 bands that are calculated we see the following results:<br><br>At ecutwfc=40 Ry:   -5.1740   7.0370   7.0370   7.0370<br>    ecutwfc=100Ry:   -5.1746   7.0369   7.0369   7.0369<br>    ecutwfc=200Ry:    -5.1746   7.0369   7.0369   7.0369<br>

Now try <br>   ecutwfc=2000Ry:       -5.2862   6.9066   6.9066  10.2399<br><br>Now why do the values change? If I look at the output file for the 2000 Ry case I see that there is a negative starting charge:<br>
<br>Initial potential from superposition of free atoms<br>     Check: negative starting charge=   -0.057614<br><br>But still, why is this? The usual answer for a negative starting charge is to increase the wavefunction cutoff, although I suspect that's not the problem in this case..... <br>
<br>So besides the using of a 2000 Ry cutoff for silicon, what else is wrong here?<br><br>Best,<br>Brad<br>UC Berkeley<br><br><br>&control<br>   prefix = 'si'<br>   calculation = 'scf'<br>   restart_mode = 'from_scratch'<br>

   wf_collect = .false.<br>   tstress = .true.<br>   tprnfor = .true.<br>   outdir = './'<br>   wfcdir = './'<br>   pseudo_dir = './'<br>/<br>&system<br>   ibrav = 0<br>   celldm(1) = 10.2612<br>

   nat = 2<br>   ntyp = 1<br>   nbnd = 4<br>   ecutwfc = 2000.0<br>/<br>&electrons<br>   electron_maxstep = 100<br>   conv_thr = 1.0d-10<br>   mixing_beta = 0.7<br>   diago_full_acc = .true.<br>/<br>CELL_PARAMETERS cubic<br>

      0.000000000  0.500000000  0.500000000<br>      0.500000000  0.000000000  0.500000000<br>      0.500000000  0.500000000  0.000000000<br>ATOMIC_SPECIES<br>  Si  28.086  Si.pz-vbc.UPF<br>ATOMIC_POSITIONS crystal<br>  Si  -0.125000000  -0.125000000  -0.125000000<br>

  Si  0.125000000  0.125000000  0.125000000<br>K_POINTS automatic<br>1 1 1 0 0 0 <br><br>