abnormal SHG coefficent

Questions and doubts about features of non linear optic in Yambo (yamb_nl)

Moderators: Davide Sangalli, claudio, myrta gruning

Post Reply
Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

abnormal SHG coefficent

Post by Dean » Wed Aug 30, 2023 11:24 am

Dear all,
I wan to calculate SHG coefficent by RT-BSE method. But, I get some abnormal SHG coefficent datas, such as 10^5 pm/v.
I find that the polarization is closely determained by the two parameters, NLtime and NLstep. Please see the figure.
My question is how to set these two parameters, NLtime and NLstep and then get the correct SHG coefficents?
Thanks in advance!
You do not have the required permissions to view the files attached to this post.
Dr. Yimin Ding
Soochow University, China.

User avatar
claudio
Posts: 438
Joined: Tue Mar 31, 2009 11:33 pm
Contact:

Re: abnormal SHG coefficent

Post by claudio » Fri Sep 01, 2023 5:45 pm

Dear

I see from the plot of your polarization, that there is a problem. The polarization should oscillate around zero and not go down.
Please make some test:
1) do you get the correct response in independent particle approximation?
2) how much smearing (dephasing) did you put in your calculation?

for the paramter NLTime and NLtstep,
try to reduce NLstep = 0.01 fs and to use a better integrator NLintegrator= "CRANKNIC"

let me know

best
Claudio
Claudio Attaccalite
[CNRS/ Aix-Marseille Université/ CINaM laborarory / TSN department
Campus de Luminy – Case 913
13288 MARSEILLE Cedex 09
web site: http://www.attaccalite.com

Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

Re: abnormal SHG coefficent

Post by Dean » Mon Sep 04, 2023 7:01 am

Dear Claudio,
Thanks for your reply.
This issue was solved by changing the version of Yambo from "Version 5.1.0 Revision 21422 Hash (prev commit) fde6e2a07" to "Version 5.1.0 Revision 22561 Hash (prev commit) 785b4fd6f".
Is ther any reference or test to set these two parameters, NLtime and NLstep and get the convergent SHG coefficents?
Best,
Dr. Yimin Ding
Soochow University, China.

User avatar
claudio
Posts: 438
Joined: Tue Mar 31, 2009 11:33 pm
Contact:

Re: abnormal SHG coefficent

Post by claudio » Tue Sep 05, 2023 9:10 am

Dear Yimin

usually the default parameters are fine.
My advice is to use the CRANCKNIC integrator.
NLtime is fixed automatically from the smearing parameter, it is inversely proportional to the smearing,
Small smearing equal long runs.
you can make some test to decrease it and speed up calculation if you see that the result does not change

best
Claudio
Claudio Attaccalite
[CNRS/ Aix-Marseille Université/ CINaM laborarory / TSN department
Campus de Luminy – Case 913
13288 MARSEILLE Cedex 09
web site: http://www.attaccalite.com

Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

Re: abnormal SHG coefficent

Post by Dean » Wed Sep 06, 2023 1:29 am

Dear Claudio,
Thanks for your reply.
I use the CRANCKNIC integrator.
I will make some test.
Best,
Dr. Yimin Ding
Soochow University, China.

Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

Re: abnormal SHG coefficent

Post by Dean » Fri Jan 26, 2024 3:57 am

Dear Claudio,
I got again abnormal SHG coefficent for another 2D materials, see below. This material show a small bandgap of about 0.5 eV.
I have used the latest version of Yambo(5.2.1). The input and output files are attahched. Please refer it.
Please give me some suggestions.
Thanks in advance!
You do not have the required permissions to view the files attached to this post.
Dr. Yimin Ding
Soochow University, China.

User avatar
claudio
Posts: 438
Joined: Tue Mar 31, 2009 11:33 pm
Contact:

Re: abnormal SHG coefficent

Post by claudio » Mon Jan 29, 2024 1:54 pm

Dear Yimin

I tried to run your calculation with slightly less converged parameters in such a way to have a fast result.

I found a correct X2, in the figure you can find the results IPA vs TD-SES for X2_yyy and X2_xyy

in attachment I put all input files, DFT and Yambo. (I used different pseudo and scissor operator, I did not have the GW results.)

please have a look and try to reproduce it, then you can increase cutoff, kpoints etc..

best
Claudio

ps use the last Yambo version on github, it is alwayer better
You do not have the required permissions to view the files attached to this post.
Claudio Attaccalite
[CNRS/ Aix-Marseille Université/ CINaM laborarory / TSN department
Campus de Luminy – Case 913
13288 MARSEILLE Cedex 09
web site: http://www.attaccalite.com

Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

Re: abnormal SHG coefficent

Post by Dean » Thu Feb 01, 2024 11:00 am

Dear Claudio,
Thanks for your help. I get resonable results by using your input files.
By comparing the input parameters, I think the abnormal SHG coefficent are caused by these parmeters, i.e., "NoDiagSC, ElecTemp= 0.0, BoseTemp=0 " in collisions and SHG calculations. I set these parmeters, because Yambo reports "metallic system", while it is a semiconductor for my system. So, I guess these parmeters, i.e., "NoDiagSC, ElecTemp= 0.0, BoseTemp=0 " will result in abnormal NL results?
Moreover, in QE runs, if I delete these parameters for a semicouductor, i.e.,"occupations='smearing', smearing='mv', degauss=0.0001," the following Yambo calculations shows correct outputs.
I wish these tips may be helpful for other users.
Best,
Dr. Yimin Ding
Soochow University, China.

Dean
Posts: 98
Joined: Thu Oct 10, 2019 7:03 am

Re: abnormal SHG coefficent

Post by Dean » Sat Feb 03, 2024 2:30 am

Dear Claudio,
After some tests, I find some unaccountable issues.
I can get resonable SHG coefficents, while the k-point of QE is set as 18*18*1, 15*15*1 or 12*12*1.....
But, abnormal SHG coefficents are produced with the same input files for Yambo, while the k-point of QE is set as 19*19*1, 20*20*1 or 21*21*1 .....
Please see the attached files.
Any suggestion is helpful. Thanks in advvance.
Best,
You do not have the required permissions to view the files attached to this post.
Dr. Yimin Ding
Soochow University, China.

Post Reply