Nonzero temperature from QE database

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
Weiqing Zhou
Posts: 48
Joined: Thu Sep 06, 2018 7:57 am
Location: Wuhan University

Nonzero temperature from QE database

Post by Weiqing Zhou » Sun Dec 09, 2018 1:47 pm

Hi all,

After initialization (yambo), I monitor the report file r_setup and find that the temperature from QE output sometimes is nonzero, always 0.02585 eV (300K)

[RD./SAVE//ns.db1]------------------------------------------
Bands : 150
K-points : 1214
G-vectors [RL space]: 67815
Components [wavefunctions]: 8548
Symmetries [spatial]: 4
Spinor components : 1
Spin polarizations : 1
Temperature [ev]: 0.02585
Electrons : 28.00000
WF G-vectors : 10323
Max atoms/species : 2
No. of atom species : 3

-----------------------------------------------------------------------
And in this case, Default of ElecTemp in all following yambo input will be 0.025 eV, too. I think this change maybe make some difference in my calculation
compared to the 0K one.
I also find nonzero temperature only occur when I set ' smearing = 'mp' ' in PWscf.
Why the temperature of the DFT ground state is not zero ? Does it matter ?

Weiqing
Weiqing Zhou
Phd student
Wuhan University
Wuhan, Hubei, 430072

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

Re: Nonzero temperature from QE database

Post by Daniele Varsano » Mon Dec 10, 2018 9:26 am

Dear Weiqing,
if you use a smearing then you will have a finite electronic temperature (it affects the occupation numbers). In presence of a gap this should not affect your results. In any case, if your system is not metallic, you can set 0K by adding in your input files (also for the setup) the variable:

Code: Select all

ElecTemp= 0.0 eV
You can also set an arbitrary temperature in case of a metallic system but 0K it is not allowed.

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