<div dir="ltr">Dear all,<div><br></div><div>I am doing some calculations on the HSE06 (25% EXX) band gap of wurtzite ZnO with experimental structure, and find that I cannot get the usual value of ~2.4 eV. My value is just ~2.1 eV.</div><div>Maybe my norm-conserving pseudopotentials (NC-PPs) of Zn and O are not very good. So I also tried the NC-PPs for Zn and O from the pslibrary.0.3.1. However, the gap is even smaller (~1.8 eV). (All PPs are based on PBE xc.)</div><div>I have also done different tests with different parameters, containing nbnd, nqx1, nqx2, nqx3, k-mesh. But I still cannot get 2.4 eV for the HSE gap.</div><div><br></div><div>So could anyone provide a guidance on how to obtain the usual HSE band gap of 2.4 eV for wurtzite ZnO ? The input parameters are as follows, and the output file calculated with the NC-PPs from the pslibrary.0.3.1 is also attached. Thank you in advance!</div><div><br></div><div>My input parameters for the scf run:</div><div><div>&CONTROL</div><div> calculation = 'scf',</div><div> prefix = 'ZnO',</div><div> pseudo_dir = './',</div><div>! nstep = 200,</div><div> verbosity = 'high',</div><div> wf_collect = .true.</div><div> etot_conv_thr = 1.0D-5,</div><div> forc_conv_thr = 1.0D-4,</div><div> restart_mode = 'from_scratch',</div><div> outdir = './temp_out',</div><div>/</div><div>&SYSTEM</div><div> ibrav = 0,</div><div>! celldm(1) = 1.88972688,</div><div> nat = 4,</div><div> ntyp = 2,</div><div> ecutwfc = 60,</div><div>! ecutrho = 240,</div><div>! nosym = .true.</div><div> nbnd = 36,</div><div>! occupations = 'smearing',</div><div>! degauss = 0.000367,</div><div>! smearing = 'gaussian',</div><div> input_dft = 'hse',</div><div> exx_fraction = 0.25,</div><div>! screening_parameter = 0.106</div><div>! exxdiv_treatment = gygi-baldereschi</div><div> nqx1 = 3, nqx2 = 3, nqx3 = 2,</div><div>! nspin = 2</div><div>! starting_magnetization(1) = 1.0,</div><div>! starting_magnetization(2) = -1.0,</div><div>! lda_plus_u = .true.</div><div>! hubbard_U(5) = 2.59,</div><div>! hubbard_alpha(2) = 0.d0,</div><div>! hubbard_U(2) = 3,</div><div>/</div><div>&ELECTRONS</div><div>! electron_maxstep = 100,</div><div> mixing_mode = 'plain',</div><div> mixing_beta = 0.7,</div><div> conv_thr = 1.D-8,</div><div>! startingwfc ='atomic+random',</div><div>! startingpot ='atomic',</div><div>/</div><div>&IONS</div><div> ion_dynamics = 'bfgs',</div><div> pot_extrapolation = 'atomic',</div><div> upscale = 100,</div><div> trust_radius_ini = 0.5,</div><div>/</div><div>&CELL</div><div> cell_dynamics = 'bfgs',</div><div> press_conv_thr = 0.5,</div><div> cell_dofree = 'all',</div><div>/</div><div><br></div><div>ATOMIC_SPECIES</div><div> Zn 65.38 Zn.pbe-nc.UPF</div><div> O 15.999 O.pbe-nc.UPF</div><div><br></div><div>CELL_PARAMETERS (angstrom)</div><div> 3.249 0.000 0.000</div><div> -1.625 2.814 0.000</div><div> 0.000 0.000 5.205</div><div><br></div><div>ATOMIC_POSITIONS (crystal)</div><div>Zn 0.333333343 0.666666687 0.000000000</div><div>Zn 0.666666627 0.333333313 0.500000000</div><div>O 0.333333343 0.666666687 0.382600009</div><div>O 0.666666627 0.333333313 0.882600009</div><div><br></div><div>K_POINTS {automatic}</div><div> 9 9 6 0 0 0</div></div><div><div><br></div><div>Regards,</div><div>Hengxin Tan, PhD</div><div>Department of Physics</div><div>Tsinghua University.</div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><pre style="color:rgb(0,0,0);line-height:16.8px;word-wrap:break-word">Tan Hengxin
Department of physics, THU.
Beijing 100084, China
Office: B403,New Science Building
E-mail:<a href="mailto:E-mail%3Athx13@mails.tsinghua.edu.cn" target="_blank">t</a><a href="mailto:anhx90@gmail.com" target="_blank">anhx90@gmail.com</a></pre></div></div></div></div></div></div></div></div>
</div></div>