Hallo to everybody.
I'm facing a strange phenomenon. I was trying to compute RPA absorption spectrum using GW energies, so I put OMEGAMAX=20 and NEDOS=2001 in my GW input file (see related topic: 'Physics questions / RPA calculation for GW ' a help also there will be appreciated).
I did not obtain the RPA spectrum but not only that! I did not even obtained the same GW corrections I computed for an exactly equivalent calculation without NEDOS nor OMEGAMAX, although these parameters should not enter in the calculation of sigma.
Here you can find INCAR and part of the OUTCAR from the two calculations:
================================
== calculation with NEDOS and OMEGAMAX ==
================================
INCAR:
NPAR = 16
ALGO = GW0
LSPECTRAL = True
LOPTICS = True
LRPA = True
PRECFOCK = FAST
PREC = High
ENCUT = 704.00 eV
NOMEGA = 100
NBANDS = 560
ENCUTGW = 704.00 eV
CSHIFT = 0.1
NBANDSGW = 48
LMAXFOCKAE = 4
NKRED = 2
ODDONLYGW = True
MAXMEM = 4000
NEDOS = 2001
OMEGAMAX = 20
CSHIFT = 0.1
ISMEAR = 0 ; SIGMA = 0.05
OUTCAR:
k-point 1 : 0.0000 0.0000 0.0000
band No. KS-energies QP-energies sigma(KS) V_xc(KS) V^pw_x(r,r') Z occupation
1 -16.8994 -17.9330 -22.6790 -21.2405 -30.2115 0.7185 2.0000
2 -4.2271 -5.6216 -24.1207 -22.2912 -30.7202 0.7622 2.0000
3 -4.2271 -5.6216 -24.1207 -22.2912 -30.7202 0.7622 2.0000
4 -4.2271 -5.6216 -24.1207 -22.2912 -30.7202 0.7622 2.0000
5 -3.0956 -4.5970 -25.3654 -23.4590 -32.4407 0.7875 2.0000
6 -3.0956 -4.5970 -25.3654 -23.4590 -32.4407 0.7875 2.0000
7 2.4616 1.0541 -23.9387 -22.1911 -27.0544 0.8054 2.0000
8 2.4616 1.0541 -23.9387 -22.1911 -27.0544 0.8054 2.0000
9 2.4616 1.0541 -23.9387 -22.1911 -27.0544 0.8054 2.0000
================================
== calculation without NEDOS and OMEGAMAX ==
================================
INCAR:
NPAR = 16
ALGO = GW0
LSPECTRAL = True
LOPTICS = True
LRPA = True
PRECFOCK = FAST
PREC = High
ENCUT = 704.00 eV
NOMEGA = 100
NBANDS = 560
ENCUTGW = 704.00 eV
CSHIFT = 0.1
NBANDSGW = 48
LMAXFOCKAE = 4
NKRED = 2
ODDONLYGW = True
MAXMEM = 4000
ISMEAR = 0 ; SIGMA = 0.05
OUTCAR:
k-point 1 : 0.0000 0.0000 0.0000
band No. KS-energies QP-energies sigma(KS) V_xc(KS) V^pw_x(r,r') Z occupation
1 -16.8994 -17.3244 -21.8446 -21.2405 -30.2115 0.7036 2.0000
2 -4.2271 -4.9498 -23.2762 -22.2912 -30.7202 0.7337 2.0000
3 -4.2271 -4.9498 -23.2762 -22.2912 -30.7202 0.7337 2.0000
4 -4.2271 -4.9498 -23.2762 -22.2912 -30.7202 0.7337 2.0000
5 -3.0956 -3.8983 -24.5191 -23.4590 -32.4407 0.7571 2.0000
6 -3.0956 -3.8983 -24.5191 -23.4590 -32.4407 0.7571 2.0000
7 2.4616 1.7613 -23.0983 -22.1911 -27.0544 0.7719 2.0000
8 2.4616 1.7613 -23.0983 -22.1911 -27.0544 0.7719 2.0000
9 2.4616 1.7613 -23.0983 -22.1911 -27.0544 0.7719 2.0000
Is this a bug or I am messing up things?
Thanks for any help!
<span class='smallblacktext'>[ Edited ]</span>
NEODS changes GW results
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 13
- Joined: Wed Nov 20, 2013 5:01 pm
- License Nr.: 5-311
NEODS changes GW results
Last edited by sponza on Wed May 21, 2014 1:19 pm, edited 1 time in total.
-- -- --
LORENZO SPONZA
LORENZO SPONZA
-
- Global Moderator
- Posts: 1817
- Joined: Mon Nov 18, 2019 11:00 am
Re: NEODS changes GW results
Hi,
We're sorry that we didn’t answer your question. This does not live up to the quality of support that we aim to provide. The team has since expanded. If we can still help with your problem, please ask again in a new post, linking to this one, and we will answer as quickly as possible.
Best wishes,
VASP