NaN in quasiparticle energies calculation with yambo-5.1.2

Concerns issues with computing quasiparticle corrections to the DFT eigenvalues - i.e., the self-energy within the GW approximation (-g n), or considering the Hartree-Fock exchange only (-x)

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

Post Reply
Gao Siyuan
Posts: 4
Joined: Tue Aug 29, 2023 3:32 pm

NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Gao Siyuan » Fri Sep 01, 2023 4:45 am

Dear developers,
I'm following the tutorial to obtain the quasi-particle band structure with yambo-5.1.2. There are on errors presented in my log files, but the values of (E-Eo) and Sc(Eo) are NaN when I set "DysSolver" to "n" or "s".
The HF quasi-particle correction calculations are normal. But X_Q_1 in ndb.pp_fragment_1 file are all "NaN".
I attached my input and report files. How can I fix it?
Thanks for your generous help!

Best,
Siyuan
You do not have the required permissions to view the files attached to this post.
Gao Siyuan
XMU, Xiamen, China

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

Re: NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Daniele Varsano » Thu Sep 07, 2023 5:11 pm

Dear Gao,
we will try to reproduce your problem, I suspect that there could be something related to the compilation, can you post also your config.log?

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/

User avatar
Davide Sangalli
Posts: 614
Joined: Tue May 29, 2012 4:49 pm
Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy
Contact:

Re: NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Davide Sangalli » Fri Sep 08, 2023 9:39 am

Dear Siyuan,
what about the ndb.dipoles. Does it contain NaNs as well?

Best,
D.
Davide Sangalli, PhD
CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
https://sites.google.com/view/davidesangalli
http://www.max-centre.eu/

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

Re: NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Daniele Varsano » Fri Sep 08, 2023 11:43 am

Dear Gao,
Can you download the new release (5.2) and see if the problem persists.

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/

Gao Siyuan
Posts: 4
Joined: Tue Aug 29, 2023 3:32 pm

Re: NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Gao Siyuan » Sat Sep 09, 2023 7:07 am

Daniele Varsano wrote: Thu Sep 07, 2023 5:11 pm Dear Gao,
we will try to reproduce your problem, I suspect that there could be something related to the compilation, can you post also your config.log?

Best,
Daniele
Dear Daniele,
Thanks for reply, here's config.log. And I'll try v5.2.
You do not have the required permissions to view the files attached to this post.
Gao Siyuan
XMU, Xiamen, China

Gao Siyuan
Posts: 4
Joined: Tue Aug 29, 2023 3:32 pm

Re: NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Gao Siyuan » Sat Sep 09, 2023 7:17 am

Davide Sangalli wrote: Fri Sep 08, 2023 9:39 am Dear Siyuan,
what about the ndb.dipoles. Does it contain NaNs as well?

Best,
D.
Dear Davide
Thanks for your reply. The ndb.dipoles does contained several NaNs.
Gao Siyuan
XMU, Xiamen, China

Gao Siyuan
Posts: 4
Joined: Tue Aug 29, 2023 3:32 pm

Re: NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Gao Siyuan » Thu Oct 05, 2023 9:10 am

I‘ve run a tutorial hBN again using Yambo5.2.0, but it still cause NaN.
Gao Siyuan
XMU, Xiamen, China

User avatar
Davide Sangalli
Posts: 614
Joined: Tue May 29, 2012 4:49 pm
Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy
Contact:

Re: NaN in quasiparticle energies calculation with yambo-5.1.2

Post by Davide Sangalli » Thu Oct 05, 2023 9:57 am

The ndb.dipoles does contained several NaNs.
So the issue is surely there.
The operation done to generate the dipoles are simple linear algebra operations.
When NaNs are there, usually it is due to a problem with the linking of blas/lapack libraries.

Can you try to recompile yambo using the flag

Code: Select all

--enable-int-linalg
when running the configure?

It will download and compile the blas and lapack library together with the code.

Best,
D.
Davide Sangalli, PhD
CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
https://sites.google.com/view/davidesangalli
http://www.max-centre.eu/

Post Reply