on RPA for metallic systems
Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan
-
- Posts: 4
- Joined: Wed Jun 23, 2010 5:34 pm
on RPA for metallic systems
Dear Yambo developers
I’m trying to reproduce results on Cu bulk at RPA level, published in PRB 64 195125 (A.Marini, G.Onida, R.Del Sole). I’m using a PP with 17 electrons without NLCC and my starting tests seem quite satisfactory. However I have not understood well how Yambo calculates the different contributions to the dielectric function yet.
I would like to know what is the meaning of the two entries in the variable DrudeWXd that Yambo should use to calculate the Drude term of the dielectric function.
According to the short description of Yambo Documentation, it is the Drude Plasmon.
1) Why are there two entries?
2) Why does Yambo need of the Drude plasma frequency? Doesn’t it calculate this frequency solving the equation 17 of the PRB 64 195125?
4) Does Yambo calculate the intraband contribution to the dielectric function by the equation 23 of the PRB 64 195125? Is this equation suitable for metallic surfaces also ?
Thank you in advance for your answers.
Best Regards.
Luca Sementa.
IPCF - CNR
Via G.Moruzzi 1, 56124 Pisa
I’m trying to reproduce results on Cu bulk at RPA level, published in PRB 64 195125 (A.Marini, G.Onida, R.Del Sole). I’m using a PP with 17 electrons without NLCC and my starting tests seem quite satisfactory. However I have not understood well how Yambo calculates the different contributions to the dielectric function yet.
I would like to know what is the meaning of the two entries in the variable DrudeWXd that Yambo should use to calculate the Drude term of the dielectric function.
According to the short description of Yambo Documentation, it is the Drude Plasmon.
1) Why are there two entries?
2) Why does Yambo need of the Drude plasma frequency? Doesn’t it calculate this frequency solving the equation 17 of the PRB 64 195125?
4) Does Yambo calculate the intraband contribution to the dielectric function by the equation 23 of the PRB 64 195125? Is this equation suitable for metallic surfaces also ?
Thank you in advance for your answers.
Best Regards.
Luca Sementa.
IPCF - CNR
Via G.Moruzzi 1, 56124 Pisa
- andrea marini
- Posts: 325
- Joined: Mon Mar 16, 2009 4:27 pm
- Contact:
Re: on RPA for metallic systems
The two fields represent the real and imaginary part of the Drude frequency. The imaginary part corresponds to the width of the peak appearing in the EELS.Luca Sementa wrote:However I have not understood well how Yambo calculates the different contributions to the dielectric function yet.
I would like to know what is the meaning of the two entries in the variable DrudeWXd that Yambo should use to calculate the Drude term of the dielectric function.
According to the short description of Yambo Documentation, it is the Drude Plasmon.
1) Why are there two entries?
Beacause, as written in the paragraph following Eq.(17) we had to use 16000 random k-points distributed very near the Fermi surface to converge Eq.(17). For this reason Yambo uses an user-defined Drude frequency plugged in a jellium-like analytic expression.Luca Sementa wrote: 2) Why does Yambo need of the Drude plasma frequency? Doesn’t it calculate this frequency solving the equation 17 of the PRB 64 195125?
Yes. You can find more information in my PhD thesis, pages 50-55 and pages 137-138.Luca Sementa wrote: 4) Does Yambo calculate the intraband contribution to the dielectric function by the equation 23 of the PRB 64 195125? Is this equation suitable for metallic surfaces also ?
In case you need more help do not hesitate in posting again here in the forum and I will try top help you.
Andrea
Andrea MARINI
Istituto di Struttura della Materia, CNR, (Italy)
Istituto di Struttura della Materia, CNR, (Italy)
-
- Posts: 4
- Joined: Wed Jun 23, 2010 5:34 pm
Re: on RPA for metallic systems
Dear Andrea
Thank you very much for your reply.
Sincerly.
Luca Sementa
IPCF-CNR
Via G.Moruzzi 56124 Pisa
Thank you very much for your reply.
Sincerly.
Luca Sementa
IPCF-CNR
Via G.Moruzzi 56124 Pisa
-
- Posts: 6
- Joined: Mon Feb 11, 2013 3:36 pm
Intraband transition
Dear All,
Hoe can i add intraband transition to optical properties?
thanks in advance
Hoe can i add intraband transition to optical properties?
thanks in advance
Masoud
Physics Department, Faculty of Science
Razi University
Physics Department, Faculty of Science
Razi University
- Daniele Varsano
- Posts: 3955
- Joined: Tue Mar 17, 2009 2:23 pm
- Contact:
Re:Intraband transitions
Dear m.shahrokhi,
fist of all let me ask you to fill your signature with your affiliation. This is a rule of the forum.
Next, in the following please post you question in the right place of the forum as it can be useful to other users and it's easier to find it. The forum it is divided in subforums divided in topics. I will move your question in the linear response subforum.
About your question, this has been answered here:
viewtopic.php?f=13&t=200&p=1146&hilit=intraband#p1146
Cheers,
Daniele
fist of all let me ask you to fill your signature with your affiliation. This is a rule of the forum.
Next, in the following please post you question in the right place of the forum as it can be useful to other users and it's easier to find it. The forum it is divided in subforums divided in topics. I will move your question in the linear response subforum.
About your question, this has been answered here:
viewtopic.php?f=13&t=200&p=1146&hilit=intraband#p1146
Cheers,
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/
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/
-
- Posts: 6
- Joined: Mon Feb 11, 2013 3:36 pm
Re: on RPA for metallic systems
Dear Daniele
Thank you very much for your quick reply.
Sincerely
Thank you very much for your quick reply.
Sincerely
Masoud
Physics Department, Faculty of Science
Razi University
Physics Department, Faculty of Science
Razi University
-
- Posts: 6
- Joined: Mon Feb 11, 2013 3:36 pm
Re: on RPA for metallic systems
Dear Daniele
How do i find the plasma frequency for a metal with yambo?
Should I use the experimental plasma frequency for metals?
How do i find the plasma frequency for a metal with yambo?
Should I use the experimental plasma frequency for metals?
Masoud
Physics Department, Faculty of Science
Razi University
Physics Department, Faculty of Science
Razi University
- Daniele Varsano
- Posts: 3955
- Joined: Tue Mar 17, 2009 2:23 pm
- Contact:
Re: on RPA for metallic systems
Dear Masoud,
yes this is the easiest way.
Fore more details please look at the Andrea Marini answer:
viewtopic.php?f=13&t=200#p1148
to Luca Sementa question.
Where to find the Drude frequency, I do not have an answer but they should be in textbook ad literature.
Cheers,
Daniele
yes this is the easiest way.
Fore more details please look at the Andrea Marini answer:
viewtopic.php?f=13&t=200#p1148
to Luca Sementa question.
Where to find the Drude frequency, I do not have an answer but they should be in textbook ad literature.
Cheers,
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/
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/
-
- Posts: 91
- Joined: Tue Sep 09, 2014 7:57 pm
Re: on RPA for metallic systems
Does the intraband contributions being included in the calculation for a metallic system even if the Drude term in the input file NOT intuduced?
Are the Drude energy and the damping values entered merely a parameter? If not, does the convergence guaranteed?
Are the Drude energy and the damping values entered merely a parameter? If not, does the convergence guaranteed?
Burak Ozdemir
Post-doc,
University of Nantes, France
Post-doc,
University of Nantes, France
- Daniele Varsano
- Posts: 3955
- Joined: Tue Mar 17, 2009 2:23 pm
- Contact:
Re: on RPA for metallic systems
Dear Burak,
the intraband contribution is *not* taken into account if not specified in the input.
When specified, in enters via the Drude model, so using empirical parameters. Ab-initio inclusion, as you argued, is very hard to converge, you can see an example in the Andrea Marini PhD thesis.
Best,
Daniele
the intraband contribution is *not* taken into account if not specified in the input.
When specified, in enters via the Drude model, so using empirical parameters. Ab-initio inclusion, as you argued, is very hard to converge, you can see an example in the Andrea Marini PhD thesis.
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/
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/