BSE calculate

Concerns any physical issues arising during the setup step (-i option). This includes problems with symmetries, k/q-point sets, and so on. For technical problems (running in parallel, etc), refer to the Technical forum.

Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, Daniele Varsano

Post Reply
Abeshipour
Posts: 11
Joined: Mon Mar 22, 2021 10:45 am

BSE calculate

Post by Abeshipour » Mon May 17, 2021 6:10 pm

Dear Yambo developers,

I'm trying to calculate the dielectric function of NaBO2 bulk with BSE according to the "Calculating optical spectra including excitonic effects: a step-by-step guide' tutorial on the website, however...

I first calculated optical spectra at the independent particle level for bulk NaBO2 according to "First steps... step4" and find that The highest peak of the imaginary part of the dielectric function of this material is about 6 eV.
But when I do calculations BSE for this system and plot the imaginary part of the dielectric function according to the output file, o-3D_BSE.eps_q1_diago_bse in the same direction x, the highest peak of this part is about
1 eV. That is, it gets smaller.
It is also non-convergent, meaning that for every value I give to the input data, it is drawn the same and does not change...

I tried very hard to understand my mistake but I could not. can you help me?

I put the input files of BSE in this section.
You do not have the required permissions to view the files attached to this post.
Aref. Abeshipour
MSc student at chamran university of Ahvaz, Iran
https://scu.ac.ir/

User avatar
Daniele Varsano
Posts: 3773
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculate

Post by Daniele Varsano » Wed May 19, 2021 7:49 am

Dear Abeshipur,
please sign your posts with your full name and affiliation. This is a rule of the forum, you can do once for all by filling your signature in you forum profile.
Next please, post your question in the appropriate subform, your question is not related with initialisation.

Coming to your question, can you explain better what is the problem?
What I can see is that you applied a scissor-shift of this kind:

Code: Select all

% KfnQP_E 
 2.50000 | 1.40000 | 0.950000 |        # [EXTQP BSK BSS] E parameters  (c/v) eV|adim|adim
% 
How did you choose these parameters?
The peak position of the spectrum depends on the scissor and stretching you choose as you are applying a rigid shit to the bands plus a stretching which is rather large.

Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Abeshipour
Posts: 11
Joined: Mon Mar 22, 2021 10:45 am

Re: BSE calculate

Post by Abeshipour » Wed May 19, 2021 10:53 am

Dear Daniele

Thanks for your reply, I Will correct my profile info ASAP.



I've asked that, in my BSE calculation, the absorption peak of dielectric function is smaller than the peak which calculated by RPA.

Is this correct or not?
If it's not correct, how can I fix it?


About your second question, I have to say:


I got the way for choosing KfnQP_E from this link

"viewtopic.php?t=1378"

according to that, I found out that to select this parameter, I should calculate the slope of the linearly QP energies with respect to KS energies.
Which they are column 3+ column 4, and column 3 in the o.qp output file.


I put Figure 1 again for two approximations. for better comparison

I also put the QP energies with respect to KS energies. To express your opinion about parameter KfnQP_E
You do not have the required permissions to view the files attached to this post.
Aref. Abeshipour
MSc student at chamran university of Ahvaz, Iran
https://scu.ac.ir/

User avatar
Daniele Varsano
Posts: 3773
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculate

Post by Daniele Varsano » Wed May 19, 2021 11:14 am

Dear Abeshipour,
your results sound strange, are the RPA and BSE calculated with the same number of bands?
Can you post your BSE output and report as well as the output/report of your GW calculations?

Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Abeshipour
Posts: 11
Joined: Mon Mar 22, 2021 10:45 am

Re: BSE calculate

Post by Abeshipour » Wed May 19, 2021 12:47 pm

Dear Daniele

Thanks for quick reply,

Yes, I do not know where the problem is
Please help me...

I put the outputs.

Thank you
You do not have the required permissions to view the files attached to this post.
Aref. Abeshipour
MSc student at chamran university of Ahvaz, Iran
https://scu.ac.ir/

User avatar
Daniele Varsano
Posts: 3773
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculate

Post by Daniele Varsano » Wed May 19, 2021 1:00 pm

Dear Href,
the spectra when compared with the Independent particle absorption looks ok.
The intensity comparison is most probably due to the different number of bands included in the two calculations.

In the GW calculation I would check the convergence of bands and blocks in the screening as a total of 70 bands with 48 occupied states seems to me it is very low.

Best,
Daniele
You do not have the required permissions to view the files attached to this post.
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Abeshipour
Posts: 11
Joined: Mon Mar 22, 2021 10:45 am

Re: BSE calculate

Post by Abeshipour » Wed May 19, 2021 1:30 pm

Dear Daniele

Thanks for your helpful comments.

I'll check the GW calculations again. It would be great if let me know your insights about the GW calculations and it's problem.

Is there any necessity to increase the number of bands or blocks in order to achieve a better convergence?


Best regards
Aref
Aref. Abeshipour
MSc student at chamran university of Ahvaz, Iran
https://scu.ac.ir/

User avatar
Daniele Varsano
Posts: 3773
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculate

Post by Daniele Varsano » Wed May 19, 2021 1:52 pm

Dear Aref,

convergence parameters are system dependent and should carefully checked.
There is not a precise rule to know in advance how many empty band are needed, but at a first glance the proportion of the bands you used with respect the occupied bands seems to me very low.

Best,
Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Abeshipour
Posts: 11
Joined: Mon Mar 22, 2021 10:45 am

Re: BSE calculate

Post by Abeshipour » Fri Sep 17, 2021 11:38 am

Dear Daniele

Thanks to you and your colleagues for all of these insightful comments/answers/recommendations, which help me to learn Yambo very well.

However, I am often doubtful about the convergence of the "BSEBands" in the BSE part.
I have calculated the "BSEBands" up to several intervals. Based on the attached convergence plot it would be great if you answer the question below:

Is the selection of "BSEBands" in the interval of (34-64) is an appropriate selection or should be charged to a larger one?

Best regards
Aref
You do not have the required permissions to view the files attached to this post.
Aref. Abeshipour
MSc student at chamran university of Ahvaz, Iran
https://scu.ac.ir/

User avatar
Daniele Varsano
Posts: 3773
Joined: Tue Mar 17, 2009 2:23 pm
Contact:

Re: BSE calculate

Post by Daniele Varsano » Fri Sep 17, 2021 12:31 pm

Dear Aref,
this depends on which part of the spectrum you are interested in and to what level of accuracy you need.
The first peak is perfectly converged, the second intense peak is is probably also converged, I suggest you to use more frequencies (BEnSteps) in order to better judge the quality of the convergence as it does not add computational effort.

Best,

Daniele
Dr. Daniele Varsano
S3-CNR Institute of Nanoscience and MaX Center, Italy
MaX - Materials design at the Exascale
http://www.nano.cnr.it
http://www.max-centre.eu/

Post Reply