<p>
Dear experts,
</p>
<p>
Focusing on spin-polarized system, I still have trouble on how to set variables for both scf and nscf calculations. And this trouble comes from BSE based calculation. Here my previous settings on my system:
</p>
<p>
&system<br>
degauss = 7.3498618000d-03<br>
ibrav = 0<br>
nat = 4<br>
nspin = 2<br>
ntyp = 2<br>
ecutwfc = 8.0000000000d+01<br>
occupations = 'smearing'<br>
smearing = 'cold'<br>
starting_magnetization(1) = 1.0000000000d-01<br>
starting_magnetization(2) = 1.0000000000d-01<br>
/<br>
&electrons<br>
conv_thr = 1.0d-12<br>
diagonalization = 'cg'<br>
diago_full_acc = .true.<br>
/<br>
ATOMIC_SPECIES<br>
C 12.01 C_ONCV_PBE_fr.upf<br>
N 14.01 N_ONCV_PBE_fr.upf<br>
ATOMIC_POSITIONS crystal<br>
N 0.0000000000 0.0000000000 0.5124828468 <br>
C 0.0000000000 0.0000000000 0.7579976664 <br>
C 0.0000000000 0.0000000000 0.0125047498 <br>
C 0.0000000000 0.0000000000 0.2670147420 <br>
CELL_PARAMETERS angstrom<br>
8.1142534611 0.0000000000 0.0000000000<br>
-4.0571267305 7.0271496300 0.0000000000<br>
0.0000000000 0.0000000000 5.0158372896<br>
K_POINTS {automatic}<br>
4 4 5 0 0 0
</p>
<p>
While I carried out BSE related calculation based on shifted k-grids (sigma,kernel and absorption), the warning occurs: WARNING: occupations outside of range [0,1] in wfn<br>
The expert of BSE response like follow:
</p>
<p>
<span class="ke-content-forecolor" style="color:#E53333;">Note that the input you take from elsewhere (e.g. your QE input) may not be suitable for GW calculations, plus it happens sometimes that those online massive-calculation project will use a uniform set of parameters for many materials, and to give roughly more or less correct results. But the accuracy in these calculations may not work for GW at all.</span>
</p>
<p>
<br>
</p>
<p>
<span class="ke-content-forecolor" style="color:#E53333;">For example, your case is a semiconductor, then you do NOT need to use smearing at all. You should use fixed occupation in this case. Also, you use the cold smearing, that's why you get occupations larger than 1.0 (you may consult the original literature). If you use Gaussian or Fermi-Dirac, this won't happen. But again, to run BerkeleyGW correctly, if you are sure your case is a semiconductor, you should use fixed occupation. (</span><span class="ke-content-forecolor" style="color:#E53333;font-family:Arial;font-size:14px;">qeinputgenerator</span><span class="ke-content-forecolor" style="color:#E53333;font-family:Arial;font-size:14px;"> for large-scale massive calculations use smearing and large smearing makes sense because the automatic code would not know if this material is metal or insulator or semiconductor beforehand. But if you already know from your test runs, you should use the best setup for your production run.</span><span class="ke-content-forecolor" style="color:#E53333;">)</span>
</p>
<p>
After this, I constrained the occupation like this:
</p>
<p>
&system<br>
ibrav = 0<br>
nat = 4<br>
nspin = 2<br>
ntyp = 2<br>
ecutwfc = 80<br>
occupations = 'fixed'<br>
tot_magnetization = 1<br>
/<br>
&electrons<br>
conv_thr = 1.0d-10<br>
mixing_mode = 'plain'<br>
mixing_beta = 0.7<br>
mixing_ndim = 8<br>
diagonalization = 'david'<br>
diago_david_ndim = 4<br>
diago_full_acc = .true.<br>
/<br>
ATOMIC_SPECIES<br>
C 12.01 C_ONCV_PBE_fr.upf<br>
N 14.01 N_ONCV_PBE_fr.upf<br>
ATOMIC_POSITIONS crystal<br>
N 0.0000000000 0.0000000000 0.5124828468 <br>
C 0.0000000000 0.0000000000 0.7579976664 <br>
C 0.0000000000 0.0000000000 0.0125047498 <br>
C 0.0000000000 0.0000000000 0.2670147420 <br>
CELL_PARAMETERS angstrom<br>
8.1142534611 0.0000000000 0.0000000000<br>
-4.0571267305 7.0271496300 0.0000000000<br>
0.0000000000 0.0000000000 5.0158372896<br>
K_POINTS {automatic}<br>
4 4 5 0 0 0
</p>
<p>
But this time when I ran the first epsilon step, the error occurs:
</p>
<p>
<span style="color:#9E0F30;">Reading header of WFN<br>
Highest occupied band (unshifted grid) = 9<br>
Valence max (unshifted grid) = .948300 eV<br>
Conduction min (unshifted grid) = -3.444951 eV<br>
Middle energy (unshifted grid) = -1.248325 eV<br>
<br>
WFN ifmin/ifmax fields are inconsistent:<br>
- there is a valence state above the middle energy<br>
- there is a conduction state below the middle energy<br>
Possible causes are:<br>
(1) Your k-point sampling is too coarse and cannot resolve the Fermi energy. Try to carefully inspect your mean-field energies, and consider using a finer k-grid.<br>
(2) You are using eqp.dat and the QP corrections change the character of some states from valence<->conduction. In this case, you should use another mean-field theory that gives the same ground state as your GW calculation.<br>
(3) You are running inteqp, but you are either shifting the Fermi energy or using restricted transformation.<br>
ERROR: WFN ifmin/ifmax fields are inconsistent</span>
</p>
<p>
It is very strange that v max is higher than middle energy and c min lower than it. The BSE expert suggested that I should check the qe outputs.
</p>
<p>
<span class="ke-content-forecolor" style="color:#E53333;">The error says your lowest conduction band has lower energy than your highest occupied band, which cannot be true. You should first check your DFT results to see whether it makes sense or not before proceeding with the GW calculations. Note that even if you have to use smearing for spin calculation, use nearly zero smearing in generating WFNs. This is required to get the correct Fermi level in BerkeleyGW.</span>
</p>
<p>
</p><p>
<span class="ke-content-forecolor" style="color:#000000;">All I could figure out is the setting is not corrected. For my system, What should I do? Thank you.</span>
</p>
<p>
<span class="ke-content-forecolor" style="color:#000000;">Sincerely,</span>
</p>
<p>
<span class="ke-content-forecolor" style="color:#000000;">Kan Zhang</span>
</p>
<br><br><hr>
<table border="0">
<tbody>
<tr style=";height:110px">
<td style="padding: 0px 3px" valign="center" width="100">
<img alt="" src="http://www.sysu.edu.cn/2012/images/sysu-mail-logo.png" style="width: 90px; height: 90px;">
</td>
<td style="padding: 0px " width="600">
<p>
<span style="font-size:10px;font-family:宋体;color:#7F7F7F">本邮件及其附件含有发送给特定个人和用于特定目的的信息。如果您不是预期的收件人，请立即删除本邮件并通知发件人。严禁任何非预期的收件人使用、传播、分发或复制本邮件或其附件。</span><br>
<span style="font-size:10px;font-family:Arial,sans-serif;color:#7F7F7F">This email and its attachments may contain confidential information intended for a specific individual and purpose. If you are not the intended recipient, you should delete this email and notify the sender immediately. Any use, dissemination, distribution, or copying of this email or its attachments by persons other than the intended recipient(s), is strictly prohibited.</span>
</p>
</td>
</tr>
</tbody>
</table>