[gpaw-users] Båndstrukturer med MGGA i GPAW

[Jens Jørgen Mortensen]

On 21-02-2012 17:42, Mie Andersen wrote:
>
> Hej Jens Jørgen
>
> Der er et problem med at udregne båndstrukturer med MGGA. Jeg har 
> vedhæftet et par slides, der illustrerer de resultater man får ved en 
> normal båndstruktursberegning og en ”workaround”, hvor jeg bruger et 
> tæt k punkt grid i den ene retning uden at fastholde densiteten.
>
> Problemet er formentlig at den kinetiske energi densitet ikke 
> fastholdes, den vil derfor blive dannet på ny ud fra en ikke-uniform k 
> punkt sampling, når man genstarter beregningen med fastholdt densitet 
> og nye k punkter.
>
> Som en midlertidig ”løsning” foreslår Lara at tilføje følgende i 
> scf.py under def run(self, wfs, hamiltonian, density, occupations): 
> assert not (self.fixdensity and hamiltonian.xc.type == 'MGGA'),"Can't 
> use fixdensity with mgga".
>
> Jeg vil gerne forsøge at finde en egentlig løsning på problemet. Hvis 
> du har nogle forslag til hvordan det kunne gøres, må du endelig sige 
> til, da jeg ikke har den store erfaring i det. Lara stopper desværre 
> her i Århus i næste uge.
>

Thank you for your explanation of the problem - I think I understand now 
what needs to be done to fix it.  I have opened a ticket:

     https://trac.fysik.dtu.dk/projects/gpaw/ticket/228

Jens Jørgen

> Mvh. Mie Andersen
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://listserv.fysik.dtu.dk/pipermail/gpaw-users/attachments/20120222/ecaf24a6/attachment.html