<div dir="ltr"><div>Dear Pascal Boulet:</div><div><br></div><div>Thank you for the information.😃 <br>It's all about the energy gap.<br>And still it is a problem we are tackling. <br></div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><br></span></div><div dir="ltr"><span style="color:rgb(0,0,0);font-family:monospace;font-size:18px">best regards</span><br style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><span style="color:rgb(0,0,0);font-family:monospace;font-size:18px">kazume NISHIDATE</span><br style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><span style="color:rgb(0,0,0);font-family:monospace;font-size:18px">敬具　西館数芽</span><br style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><br style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><span style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><a href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a></span><br style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><span style="color:rgb(0,0,0);font-family:monospace;font-size:18px"><a href="mailto:kazume.nishidate@gmail.com" target="_blank">kazume.nishidate@gmail.com</a></span><br></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">2021年11月5日(金) 18:57 pboulet <<a href="mailto:pascal.boulet@univ-amu.fr">pascal.boulet@univ-amu.fr</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div style="word-wrap:break-word;line-break:after-white-space">Dear all,<div><br></div><div>Regarding the ‘scissor’ operator, when I was a PhD student working on molecular systems with programs using Slater types of orbitals (ADF), I sometimes had to use the so-called ‘level-shift’ trick to facilitate the SCF convergence (it avoids the swapping of occupied/unoccupied orbitals from one step to the next when the gap is small). If I am right, correct me if I am wrong, this level-shift trick ressembles a lot the ‘scissor’ operator. </div><div><br></div><div>Best,</div><div>Pascal</div><div><br></div><div><br><div>
<div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;word-wrap:break-word;line-break:after-white-space"><div dir="auto" style="text-align:start;text-indent:0px;word-wrap:break-word;line-break:after-white-space"><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><span style="font-size:16pt;font-family:Mistral">Pascal Boulet</span></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font face="Lucida Handwriting"><span style="font-size:21px">—</span></font></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><em style="color:rgb(34,187,234);font-family:verdana;font-size:11px;line-height:22px">Professor in computational materials chemistry - DEPARTMENT OF CHEMISTRY</em></div><div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Tahoma,Verdana,Arial,Helvetica,sans-serif;font-size:15px;line-height:22px"><span style="font-size:11px;font-family:verdana">University of Aix-Marseille </span><span style="font-size:11px;font-family:verdana">- Avenue Escadrille Normandie Niemen - F-13013 Marseille - FRANCE</span></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none;font-family:Tahoma,Verdana,Arial,Helvetica,sans-serif;font-size:15px;line-height:22px"><span style="font-size:11px;font-family:verdana">Tél: +33(0)4 13 55 18 10 - Fax : +33(0)4 13 55 18 50</span></div><div style="color:rgb(0,0,0);letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><font face="verdana"><span style="font-size:12px;line-height:22px">Email : </span></font><font color="#22bbea" face="verdana"><span style="font-size:12px;line-height:20px"><a href="mailto:pascal.boulet@univ-amu.fr" target="_blank">pascal.boulet@univ-amu.fr</a></span></font></div><div style="font-size:14px"><br></div></div></div></div></div><br></div><br></div><br></div><br></div><br></div><br></div><br><br>
</div>
<div><br><blockquote type="cite"><div>Le 5 nov. 2021 à 02:57, Kazume NISHIDATE <<a href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a>> a écrit :</div><br><div><div dir="ltr"><div><br></div><div><br>> So, a new calculation of the bandgap was performed with HSE to obtain<br>> the difference between Egap HSE and Egap PBE.<br><br>Ok, they used the 'scissor' operator to correct the PBE dielectric constant.<br><br>> If I remember the explanation, this value is added in the optical<br>> results to push all the results obtained with PBE toward the value<br>> that could be found using the HSE functional.<br><br>And they concluded that the correction was small despite the<br>underestimation tendency of the E_g_PBE (left column, p4).<br><br>> What do you suggest to me? What value should I use instead of the<br>> 2x2x2?<br><br>The 2x2x2 k-mesh must be used for the 'scf' calculation.  And then you </div><div>should try to evaluate the band structure along the specific k-point<br>path with the 'bands' (non-scf) calculation. This is the procedure they</div><div>used to evaluate the bands.<br><br>Have you ever done a band structure calculation of simple<br>semiconductor crystal? You should once try to draw a band structure </div><div>of Si (diamond structure) so that you can understand the procedure.<br><br></div><div><br></div><div><div dir="ltr"><div dir="ltr"><span style="font-family:monospace;font-size:18px"><br></span></div><div dir="ltr"><span style="font-family:monospace;font-size:18px">best regards</span><br style="font-family:monospace;font-size:18px"><span style="font-family:monospace;font-size:18px">kazume NISHIDATE</span><br style="font-family:monospace;font-size:18px"><span style="font-family:monospace;font-size:18px">敬具　西館数芽</span><br style="font-family:monospace;font-size:18px"><br style="font-family:monospace;font-size:18px"><span style="font-family:monospace;font-size:18px"><a href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a></span><br style="font-family:monospace;font-size:18px"><span style="font-family:monospace;font-size:18px"><a href="mailto:kazume.nishidate@gmail.com" target="_blank">kazume.nishidate@gmail.com</a></span><br></div></div></div><br><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">2021年11月5日(金) 9:23 José Xavier <<a href="mailto:jxln_karate@yahoo.com.br" target="_blank">jxln_karate@yahoo.com.br</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Dear, <br>
About the 'scissor', it was used because a reviewer asked to also perform the calculations with a hybrid functional. So, a new calculation of the bandgap was performed with HSE to obtain the difference between Egap HSE and Egap PBE. If I remember the explanation, this value is added in the optical results to push all the results obtained with PBE toward the value that could be found using the HSE functional. So, 3.8 eV is the bandgap with PBE, while the HSE result is 4.7 eV.<br>
<br>
>In other wards, you should correctly plot the band structure along typical k-points to<br>
find the 'reasonable' band gaps.<br>
<br>
What do you suggest to me? What value should I use instead of the 2x2x2?<br>
<br>
<br>
Best Wishes<br>
José Xavier<br>
<br>
Em quinta-feira, 4 de novembro de 2021 07:27:48 BRT, Kazume NISHIDATE <<a href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a>> escreveu: <br>
<br>
<br>
<br>
<br>
<br>
Dear Stefano:<br>
> For isolated molecules, smearing helps dumping the convergence<br>
> oscillations arising from HOMO-LUMO level crossings.<br>
<br>
I did not know it.<br>
Thank you for the information.<br>
<br>
<br>
Dear José Xavier,<br>
> The band gap for Dopamine was 3.8 eV, and the other papers published<br>
>  by them showed values between 3-5 eV. That is why I have asked<br>
<br>
According to the paper you mentioned, they utlized the 'scissor'<br>
correction on the band gap,<br>
<br>
 E_g_scissor = E_g_HSE - E_g_PBE<br>
<br>
This should be the possible reason of the discrepancies.<br>
<br>
It is somewhat empirical way, but surprisingly, it looks me work for some<br>
transition oxide materials, according to their paper. As far as I know,<br>
the 'scissor' method is not well recognized in the DFT community.<br>
#This is the first time I have ever heard of.<br>
<br>
And there is one other point. The valence band maximum (VBM), which<br>
corresponds to the HOMO, is located at the Y point in the BZ (see<br>
Fig.2 of the paper).  To estimate correctly the gap from the DOS<br>
analysis, it must include the contribution from the Y point.  The<br>
2x2x2 mesh you mentioned is apparently insufficient.  In other wards,<br>
you should correctly plot the band structure along typical k-points to<br>
find the 'reasonable' band gaps.<br>
<br>
<br>
<br>
best regards<br>
kazume NISHIDATE<br>
敬具　西館数芽<br>
<br>
<a href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a><br>
<a href="mailto:kazume.nishidate@gmail.com" target="_blank">kazume.nishidate@gmail.com</a><br>
<br>
<br>
<br>
</blockquote></div></div>
_______________________________________________<br>Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</a>)<br>users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></blockquote></div><br></div></div>_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>