<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Cameron,<div class=""><br class=""></div><div class="">I don’t think that this strange behavior is due to mixing_beta or conv_thr, so you can remove those flags from the input. Moreover, the vacuum space you use (~27 A) is a bit large, you can slightly reduce it to speedup the calculation.</div><div class=""><br class=""></div><div class="">As far as I can see you should: 1) check if the lattice parameter is correct, for some materials it can happen that if you use a lattice parameter quite different from the equilibrium one the band structure undergoes more or less significant changes; 2) check your k-point path, are you sure you gave the correct k-points in the correct units? Maybe seeing the input for the nscf calculation would be useful as well.</div><div class=""><br class=""></div><div class="">Finally, even though I do not think that this might affect your band structure calculation, if you include one or more ultrasoft pseudo potentials in you calculations, as you do, typical values of ecutrho, depending on the atomic element and on the pseudo potential, can range from ~ 6*ecutwfc to even 12*ecutwfc, so your choice 300 Ry ~ 4.3 * 70 Ry might produce not fully converged results. </div><div class=""><br class=""></div><div class="">Last, if I were you, I would replace the norm conserving pseudo potential for S with an ultrasoft pseudo potential, which allows you to use much lower values of ecutwfc.</div><div class=""><br class=""></div><div class="">I attach a band structure of monolayer MoS2 that has been done using QE 5.2.0, where the direct gap at K (~1.74 eV) is evident.</div><div class=""><br class=""></div><div class="">Giovanni</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><img apple-inline="yes" id="05D51B95-C7DE-4894-875B-B546075AC775" height="214" width="320" apple-width="yes" apple-height="yes" src="cid:EE569A64-97F2-4918-BC90-3E05B0A5AE5B" class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 24 Sep 2015, at 08:20, Cameron Foss <<a href="mailto:cjfoss@umass.edu" class="">cjfoss@umass.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Hello,<div class=""><br class=""></div><div class="">I am attempting to obtain the band structure of monolayer MoS2, however my results so far do not agree well with what has been documented for monolayer (ML) MoS2. I consistently get an indirect band gap for ML MoS2 caused by the 'Q-valley' minimum seen in bulk MoS2.</div><div class=""><br class=""></div><div class=""><div class="">I have tested the code with interplanar distances of 15 and 30 angstroms to eliminate any interplanar interactions. There have been several papers that report the band structure for ML MoS2 from first principles, they have used functionals from LDA, GGA, and an HSE hybrid functional. My best success has been with GGA-PW91(see input file below). Note in all my simulations I have measured a band gap of ~1.9eV which is in agreement with reported values.</div><div class=""><br class=""></div><div class="">I am uncertain as to what the issue could be? I have tested 3 different PPs. </div></div><div class=""><br class=""></div><div class="">%%%%%%%%%%%%%%%%%%%%</div><div class="">Details on simulations:</div><div class=""><br class=""></div><div class=""> Specifically, my simulations continue to show that the 'Q-valley' is slightly lower (16meV~30meV) than the minimum at the K point, suggesting an indirect band gap. However ML mos2 has been widely shown to have a direct band gap at the K point where this 'Q-valley' along the conduction band is suppressed above the K-point minimum.</div><div class=""><br class=""></div><div class="">I run a scf (pw.x) calculation first, followed by an nscf (pw.x) run, and a bands calculation with pw.x (and of course bands.x to retrieve the energy values). I specify 16 bands for nbnd, kpoints along symmetry paths, and a MP grid of 8x8x1 with a 1 1 1 offset. </div><div class=""><br class=""></div><div class="">Input file: </div><div class=""><br class=""></div><div class=""><div class="">&control</div><div class=""> calculation='scf' !nscf, then bands</div><div class=""> restart_mode='from_scratch',</div><div class=""> !pseudo_dir='directory where pseudopotentials are stored/',</div><div class=""> !outdir='directory where large files are written/'</div><div class=""> pseudo_dir='/home/cameron/QE/espresso-5.1/pseudo/',</div><div class=""> outdir='/home/cameron/QE/espresso-5.1/2dout'</div><div class=""> prefix='mos2-gga',</div><div class=""> /</div><div class=""> &system</div><div class=""> ibrav=4, celldm(1)=5.8364, celldm(3)=10,</div><div class=""> nat=3, ntyp=2, ecutwfc =70, ecutrho=300 ! specified nbnd=16</div><div class=""> /</div><div class=""> &electrons</div><div class=""> conv_thr = 1.0d-15 ! reduced to 1.0d--12 for nscf, bands</div><div class=""> mixing_beta = 0.5 ! increased to 0.7 for nscf, bands</div><div class=""> /</div><div class="">ATOMIC_SPECIES</div><div class=""> Mo 95.94 Mo.pw91-n-van.UPF</div><div class=""> S 32.065 S.pw91-n-mt.UPF</div><div class="">ATOMIC_POSITIONS bohr</div><div class="">S 0.00002 -0.00001 0.000000000</div><div class="">Mo 2.91822 1.68481 2.929886569</div><div class="">S 0.00002 -0.00001 5.859870374</div><div class="">K_POINTS automatic</div><div class=""> 8 8 1 1 1 1</div></div><div class=""><br class=""></div><div class="">I ran into some convergence issues with all the bands with this conv_thr in the nscf and bands calculations, so i had to reduce it to 1.0d-12 for those calculations. I also used a mixing beta of 0.7 for both the nscf and bands calculations. I am uncertain as to whether changing mixing_beta's between calculations is of importance or not as well, but I do not think this is the case.</div><div class=""><br class=""></div><div class="">Best regards,</div><div class="">Cameron</div></div>
_______________________________________________<br class="">Pw_forum mailing list<br class=""><a href="mailto:Pw_forum@pwscf.org" class="">Pw_forum@pwscf.org</a><br class="">http://pwscf.org/mailman/listinfo/pw_forum</div></blockquote></div><br class=""><div apple-content-edited="true" class="">
-- <br class=""><br class="">Giovanni Cantele, PhD<br class="">CNR-SPIN<br class="">c/o Dipartimento di Fisica<br class="">Universita' di Napoli "Federico II"<br class="">Complesso Universitario M. S. Angelo - Ed. 6<br class="">Via Cintia, I-80126, Napoli, Italy<br class="">e-mail: <a href="mailto:giovanni.cantele@spin.cnr.it" class="">giovanni.cantele@spin.cnr.it</a><br class="">Phone: +39 081 676910<br class="">Skype contact: giocan74<br class=""><br class="">ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" class="">http://www.researcherid.com/rid/A-1951-2009</a><br class="">Web page: <a href="http://people.na.infn.it/~cantele" class="">http://people.na.infn.it/~cantele</a><br class="">
</div>
<br class=""></div></body></html>