NaN in GW numerical calculation

You can find here problems arising when using old releases of Yambo (< 5.0). Issues as parallelization strategy, performance issues and other technical aspects.

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

Locked
zhouym
Posts: 1
Joined: Thu Nov 10, 2022 3:32 am
Location: Kaifeng, Henan 475001, China

NaN in GW numerical calculation

Post by zhouym » Sat Nov 19, 2022 6:33 am

Dear developers,

I'm using Yambo 5.1.1 to perform the example of hBN. When I use the GW of PPA approximation : the values of (E-Eo) and Sc(Eo) are NaN for each K-point and each band.

o-gw_PPA.qp:
---------------------------------------------------------------------------------------------------------------------------------------------------------------
# K-point Band Eo [eV] E-Eo [eV] Sc|Eo [eV]
#
7 8 -0.411870 NaN NaN
7 9 3.877799 NaN NaN
---------------------------------------------------------------------------------------------------------------------------------------------------------------
r-gw_PPA_HF_and_locXC_gw0_dyson_em1d_ppa:
---------------------------------------------------------------------------------------------------------------------------------------------------------------
QP [eV] @ K [7] (iku): 0.000000 -0.500000 0.000000
B=8 Eo= -0.41 E= NaN E-Eo= NaN Re(Z)= NaN Im(Z)= NaN nlXC=-19.12909 lXC=-16.10656 So= NaN
B=9 Eo= 3.88 E= NaN E-Eo= NaN Re(Z)= NaN Im(Z)= NaN nlXC= -5.53657 lXC=-10.66983 So= NaN
---------------------------------------------------------------------------------------------------------------------------------------------------------------

There was no error message in the process of software compilation and calculation, and there was no problem when running on the virtual machine given by the official website. The problem of input file can be excluded. Can you give me some advice?

Best,

Zhou
You do not have the required permissions to view the files attached to this post.

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

Re: NaN in GW numerical calculation

Post by Daniele Varsano » Wed Nov 23, 2022 8:12 am

Dear Zhou,

it seems a compilation problem, you can check if the NaN is already present in the ndb.dipoles and/or ndb.pp filesby using the grep command.
You can try to recompile your code using the flag --enable-int-linalg in the configure process, even if it seems that the linear algebra is not the problem as the exchnage part of the self-energy is correctly computed.

If this does not solve the problem, please report here.

Please add your full name and affiliation in your signature, this is a rule of the forum, you can do it once for all by filling the signature in the user profile.

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/

Locked