<html><body>
<div style="font-size: 16px; font-family: Tahoma; color: rgb(0, 0, 0); font-weight: 400; font-style: normal;">Dear Jose,</div>
<div style="font-size: 16px; font-family: Tahoma; color: rgb(0, 0, 0); font-weight: 400; font-style: normal;">Usually a "single-valued" question generates a single-valued answer.</div>
<div style="font-size: 16px; font-family: Tahoma; color: rgb(0, 0, 0); font-weight: 400; font-style: normal;">But as you see, yours was not so.</div>
<div style="font-size: 16px; font-family: Tahoma; color: rgb(0, 0, 0); font-weight: 400; font-style: normal;"> </div>
<div style="font-size: 16px; font-family: Tahoma; color: rgb(0, 0, 0); font-weight: 400; font-style: normal;">Bests,</div>
<div style="font-size: 16px; font-family: Tahoma; color: rgb(0, 0, 0); font-weight: 400; font-style: normal;">Mahmoud</div>
<div style="font-size: 16px; font-family: Tahoma; color: rgb(0, 0, 0); font-weight: 400; font-style: normal;"> </div>
<div style="font-family: Tahoma; color: rgb(0, 0, 0); font-size: 16px; font-weight: 400; font-style: normal;">
<blockquote class="as91pFk8AX8sDPOUp66RUy8IoCRjefxNG8oLk8fg" style="PADDING-RIGHT: 0px; PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #000000 2px solid; MARGIN-RIGHT: 0px">From: Kazume NISHIDATE <nisidate@iwate-u.ac.jp><br />
To: José Xavier <jxln_karate@yahoo.com.br><br />
Cc: Quantum ESPRESSO users Forum <users@lists.quantum-espresso.org><br />
Date: Thu, 4 Nov 2021 15:06:17 +0900<br />
Subject: Re: [QE-users] occupation in biological molecules and band gap<br />
<div dir="ltr">
<div dir="ltr">
<div>> <span class="gmail-Apple-converted-space"> </span>would like to ask if you could tell me why this question generates so different answers.<span class="gmail-Apple-converted-space"> </span></div>
<div>It is partly because you did not provide the detailed information on your calculation.</div>
<br />
I supposed it is a crystal composed of molecules and you applied some<br />
k-point mesh to incorporate its periodicity. I thought that the system<br />
might be in a semiconducting state since you were interested in the<br />
gap value. In that case, you should apply the 'smearing' method.<br />
<br />
If you are calculating an isolated molecule in vacuume with only the<br />
gamma for the k-point, then the 'fixed' method is a good choice since<br />
the occupancy can be defined explicitly at the gamma.<br />
<br />
- note<br />
The 'nbnd' is the number of band. You should include sufficient empty<br />
bands to draw a nice band/dos structure.
<div> </div>
<div> </div>
<div>
<div>
<div class="gmail_signature" data-smartmail="gmail_signature" dir="ltr">
<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></div>
</div>
</div>
</div>
</div>
<div class="gmail_quote">
<div class="gmail_attr" dir="ltr">2021年11月4日(木) 14:09 José Xavier <<a href="mailto:jxln_karate@yahoo.com.br">jxln_karate@yahoo.com.br</a>>:</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 />
<br />
Thank you for your explanation.<br />
I'll try both occupancies to see if the results of bandgap and DOS are very different. When I got the values, I can post the differences here.<br />
<br />
I would like to ask if you could tell me why this question generates so different answers. As just like in the tutorials I saw, each one gave me a different answer.<br />
<br />
Best wishes,<br />
José Xavier<br />
<br />
<br />
Em terça-feira, 2 de novembro de 2021 00:41:19 BRT, Mahmoud Payami Shabestari <<a href="mailto:mpayami@aeoi.org.ir" target="_blank">mpayami@aeoi.org.ir</a>> escreveu:<br />
<br />
<br />
<br />
<br />
<br />
<br />
Dear Xavier,<br />
<br />
<br />
<br />
> So, I can use the Fixed occupation and, if the energy levels show a small bandgap, I introduce the >Smearing. Is it right?<br />
<br />
<br />
<br />
Yes, exactly.<br />
<br />
<br />
<br />
<br />
>If yes, I just have to add the nbnd function in the input file?<br />
<br />
<br />
<br />
Yes. It suffices to put nbnd equal to N/2 + 1 (N is number of electrons).<br />
<br />
<br />
<br />
Mahmoud<br />
<br />
> From: xavier neto via users <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br />
><br />
> To: Mahmoud Payami Shabestari <<a href="mailto:mpayami@aeoi.org.ir" target="_blank">mpayami@aeoi.org.ir</a>>, Quantum ESPRESSO users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br />
> Date: Mon, 1 Nov 2021 21:04:31 -0300<br />
> Subject: [QE-users] occupation in biological molecules and band gap<br />
> <br />
> Dear Mahmoud Payami,<br />
><br />
> Thank you for your answer. It helps me a lot.<br />
><br />
> So, I can use the Fixed occupation and, if the energy levels show a small bandgap, I introduce the Smearing. Is it right?<br />
><br />
> <br />
><br />
> About the second question. I've read that the result of HOMO - LUMO, obtained after the SCF calculation, can give the bandgap. Can I obtain the bandgap in that way? If yes, I just have to add the nbnd function in the input file?<br />
><br />
> On 11/1/21 02:04, Mahmoud Payami Shabestari wrote:<br />
><br />
>> <br />
>> Dear José Xavier,<br />
>><br />
>> Hi.<br />
>><br />
>> Smearing for the occupations is used when the levels near HOMO are so close (commonly in metals near Fermi level) that in consecutive iterations they may exchange positions and lead to nonconvergence. For semiconductors with 3-5 eV gap I think it won't help. <br />
>><br />
>> Bests<br />
>><br />
>> <br />
>> Mahmoud Payami<br />
>><br />
>> NSTRI, AEOI, Tehran, Iran<br />
>><br />
>> <br />
>><br />
>><br />
>> <br />
>>> From: José Xavier via users <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br />
>>> To: "<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>" <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br />
>>> Date: Mon, 1 Nov 2021 00:18:55 +0000 (UTC)<br />
>>> Subject: [QE-users] occupation in biological molecules and band gap<br />
>>> <br />
>>> Dear all,<br />
>>><br />
>>> I'm new in QE. The group that I'm working with usually calculates the properties of biologically relevant molecules, like amino acids, hormones, neurotransmitters, and drugs, and I would like to introduce the QE code to the group. Their previous papers have shown that the bandgap of these crystals is something between 3-5 eV, which I understand to be in the semiconductor range.<br />
>>><br />
>>> I've watched/read some tutorials about QE, and some of these showed that the occupation of semiconductors should be "smearing", other ones that it should be "fixed", but no one showed calculations with this kind of molecules. What value should I use for the SCF and electronic calculations of the crystals that I'm going to work on (amino acids, hormones, neurotransmitters, and drugs)?<br />
>>><br />
>>> *Are there any differences when the molecule has a metal, like the Heme group?<br />
>>><br />
>>> Besides, I would like to ask if there is a calculation that I could perform to obtain the bandgap of these molecules. Is it only possible if I create the band structure figure and calculate the difference between the peaks in VB and CB?<br />
>>><br />
>>> Thank you for your help,<br />
>>><br />
>>> Sincerely,<br />
>>> José Xavier<br />
>>><br />
>>> _______________________________________________<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><br />
>>><br />
>><br />
>><br />
><br />
<br />
</blockquote>
</div>
</div>
</blockquote>
</div>
</body></html>