luminescence calculation

Run-time issues concerning Yambo that are not covered in the above forums.

Moderators: myrta gruning, andrea marini, Daniele Varsano, Conor Hogan

SAM
Posts: 19
Joined: Tue Jan 30, 2024 3:36 pm

luminescence calculation

Post by SAM » Wed Apr 17, 2024 6:29 pm

dear user,
I have a problem with the luminescence calculation. By following the tutorials the last step causes me a problem the step where you have to obtain the luminescence from the yambo_ph -excph o command. when I compile it I get the error saying that the el-ph data is not found. While I executed the other previous steps successfully.

my second concern about the q2r.x calculation of quantum espresso. In the tutorials on the calculation of luminescence assisted by phonon, after the calculation of phonon it is necessary to calculate q2r.x, matdyn. however, before executing q2r you must have generated the dyn0 as well as the dyn* obtained from the phonon calculation. my question is how to generate the dyn0 in order to do the q2r.x calculation.

SAM
PhD Student at university of Lome

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: luminescence calculation

Post by claudio » Fri Apr 19, 2024 4:20 pm

Dear SAM

I repeated all the steps of the tutorial and it works perfectly.
I used QE 7.2 and the Yambo version of the tutorial, and compiled them with gfortran.
Some additional comments:

1) in the original tutorial q2r.x is never used, please check the script run_dvscf.sh.
q2r.x is used only if you want to use the double-grid. I advice you to run the tutorial without double-grid
and then if you want you can add it.

2) I added the YamboDatabases including electron-phonon in the tutorial web page, you can download them if you want.
https://www.yambo-code.eu/wiki/index.ph ... minescence

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

SAM
Posts: 19
Joined: Tue Jan 30, 2024 3:36 pm

Re: luminescence calculation

Post by SAM » Mon Apr 22, 2024 4:31 pm

Dear claudio,
Thank you very much for your answer. The difficulty I encountered was that when I compiled the bash ./rundvscf: I encountered the error: ypp_ph: unrecognized option '-noopenmp'.
the difficulty compared to the q2r.x calculation was based on the phonon-assisted photoluminescence tutorials. It was necessary to calculate matdyn.x. to generate the Q.modes.
my second difficulty lies in the calculation of finite difference after obtaining the dipoles. is this done externally?

thanks
SAM
PhD student at university of Lome/Togo

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: luminescence calculation

Post by claudio » Wed Apr 24, 2024 1:03 pm

Dear SAM

1) remove the option "-noopenmp" form the script run_dvscf.sh, but this is strage, if you do ypp_ph --help
do you see this option?

2) for the second tutorial
https://www.yambo-code.eu/wiki/index.ph ... placements
I will check it and let you 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

SAM
Posts: 19
Joined: Tue Jan 30, 2024 3:36 pm

Re: luminescence calculation

Post by SAM » Wed Apr 24, 2024 1:14 pm

Dear Claudio,

Thank you for answering me .
when I compile ypp_ph I get *nompi*

I have a complaint about the tutorials on phonon-assisted luminescence. Is it possible to leave us a bash script which automates all the calculation up to the calculation of finite difference, as is done in the tutorials on luminescence. Because it seems a bit complicated to me. Maybe looking at the automated bash will help me.
THANKS
SAM
PhD student at University of Lome/Togo

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: luminescence calculation

Post by claudio » Thu Apr 25, 2024 2:12 pm

Dear SAM

for the tutorial at finite differences, I fixed the wiki webpage, now you can download all inputs and run calculations
for phonons as explained in the tutorial, then run q2r.x

https://www.yambo-code.eu/wiki/index.ph ... transition

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

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: luminescence calculation

Post by claudio » Thu Apr 25, 2024 2:13 pm

You can remove also -nompi from the input,

regarding the automatic script I will have a look and let you 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

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: luminescence calculation

Post by claudio » Thu Apr 25, 2024 2:15 pm

Anyway at the end of the tutorial all inputs are provided


https://www.yambo-code.eu/wiki/index.ph ... s_for_h-BN
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

SAM
Posts: 19
Joined: Tue Jan 30, 2024 3:36 pm

Re: luminescence calculation

Post by SAM » Thu Apr 25, 2024 4:14 pm

Dear Claudio,

thank you very much for adding more detailed information to the tutorials, I think it's a bit clear for me.
My second question was: how is the difference calculated? Is this done externally?
is it possible to use yambopy to do this?
We will have a total of 12 dipoles with which we must calculate the finished difference.


SAM
PhD student at university of Lome/Togo

User avatar
claudio
Posts: 458
Joined: Tue Mar 31, 2009 11:33 pm
Location: Marseille
Contact:

Re: luminescence calculation

Post by claudio » Sat Apr 27, 2024 5:13 pm

Dear SAM

first of all I added some additional script in the wiki

https://www.yambo-code.eu/wiki/index.ph ... placements

the derivative of dipoles are calculated in this way:
1) the lowest dipole are zero because the transition are indirect (in the unit cell)
2) when you displace atoms along a phonon mode you activate a dipole, this mean that the dipole that before was zero, not it is different from zero.
Therefore the derivative of the diple id d Dip/ d u = [ Dip(diplaced) - 0.0 ] / du
at the end you can take as derivative the value you obtain when you displace atoms. You should divide for "du" but it is just a multiplicative factor.

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

Post Reply