[gpaw-users] Intraband contributions in dielectric function calculations

[Sharmila Shirodkar]

Dear All

I have looked through the pair.py and chi0.py files in gpaw/response
folders.
The condition for adding intraband contributions to the dielectric or
response matrix
is given by

optical_limit = not self.no_optical_limit and np.allclose(q_c, 0.0) ------
(1)

Implying that the intraband contributions are added only when q_c = (0,0,0).
Hence the dielectric matrix only contains the interband transitions.
However, the paper "Plasmons in metallic monolayer and bilayer transition
metal dichalcogenides",
Andersen et al, PHYSICAL REVIEW B 88, 155128 (2013) does report plasmon
dispersion with intraband contributions at non-zero q_c.

I also changed the condition (1) to

optical_limit = not self.no_optical_limit and self.intraband ----- (2)

And made appropriate changes to chi0.py and pair.py files to accommodate
condition (2), as given below.

A) In chi0.py
    Line 139
    replace  if np.allclose(q_c, 0.0):
    with   if self.intraband:

    Line 167
    replace optical_limit = not self.no_optical_limit and np.allclose(q_c,
0.0)
    with   optical_limit = not self.no_optical_limit and self.intraband

B) In pair.py
    Line 757
    replace optical_limit = not disable_optical_limit and np.allclose(q_c,
0.0)
    with      optical_limit = not disable_optical_limit and intraband
#SHARMILA

    Line 916
    replace assert np.allclose(pd.kd.bzk_kc[0], 0.0)
    with assert np.allclose(0.0, 0.0)

However, I was unable to reproduce the results of the above mentioned paper.
I observe that potential of the eigenmodes have a curvature (implying that
charge is not conserved, rho is non-zero (0.005) in vacuum).
Please find the plots of potential and charge density (rho) of the plasmon
eigenmodes attached to this email.

Any advice will be appreciated.

Best Regards
Sharmila Shirodkar

Postdoctoral Research Fellow
SEAS, Harvard University
Cambridge, MA 02138
USA
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listserv.fysik.dtu.dk/pipermail/gpaw-users/attachments/20160323/c5b553cd/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: potential.png
Type: image/png
Size: 206175 bytes
Desc: not available
URL: <http://listserv.fysik.dtu.dk/pipermail/gpaw-users/attachments/20160323/c5b553cd/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rho.png
Type: image/png
Size: 215841 bytes
Desc: not available
URL: <http://listserv.fysik.dtu.dk/pipermail/gpaw-users/attachments/20160323/c5b553cd/attachment-0003.png>