<div dir="ltr"><br><div class="gmail_extra">Thank you Prof. Gerhard for taking concern over the issue!</div><div class="gmail_extra"><br></div><div class="gmail_extra">I really do not know whats happening inside the code about the concept of sqrt(2).</div><div class="gmail_extra"><br></div><br><div class="gmail_extra"><br></div><div class="gmail_extra"> <br></div><div class="gmail_extra"><span style="font-size:12.8px"></span></div><div class="gmail_extra"><br></div><div class="gmail_extra">But I tried all possible combinations. With sqrt(2), without sqrt(2). With TETRA 101 or 0.000.</div><div class="gmail_extra"><br></div><div class="gmail_extra">The w_pl is same for up and dn spin.</div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_extra">In one test I saw in case.outputjointup/dn</div><div dir="ltr"><span style="font-size:12.8px"><br></span></div><span style="font-size:12.8px"> !!! WARNING: w_pl = sqrt( w_pl^2(up-spin) + w_pl^2(dn-spin) )!!!</span></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra">Regards</div><div class="gmail_extra"><br></div><div class="gmail_extra">Bhamu<br></div><div class="gmail_extra"> <br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br></div><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">The optics programm tells about the SP Plasma frequencies:<br>
w_pl = sqrt( w_pl^2(up-spin) + w_pl^2(dn-spin) )<br></blockquote><div><br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
where does this come from<br>
remember the classical approach where the plasma frequency is proportional to root(n) with n being the density of free electrons (number of electrons per Volume)<br></blockquote><div><br></div><div><br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
now what do you have in the SP case ? Shouldn't that be root/n_up+n_dn) ?<br>
The programm however delivers w_up proportional to root(n_up) and w_dn prop to root(n_dn) and you see where the above equation is comming from<br>
<br>
Finally show what happens when the two densities are equal n_up=n_dn and you see where the root(2) is coming from and when to use it<br>
<br>
(Note that the root(2) is here not because of surface plasmons)<br>
<br>
Ciao<br>
Gerhard<br>
<br>
DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br>
"I think the problem, to be quite honest with you,<br>
is that you have never actually known what the question is."<br>
<br>
==============================<wbr>======<br>
Dr. Gerhard H. Fecher<br>
Institut of Inorganic and Analytical Chemistry<br>
Johannes Gutenberg - University<br>
55099 Mainz<br>
and<br>
Max Planck Institute for Chemical Physics of Solids<br>
01187 Dresden<br>
______________________________<wbr>__________<br>
Von: Wien [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.<wbr>tuwien.ac.at</a>] im Auftrag von Gavin Abo [<a href="mailto:gsabo@crimson.ua.edu">gsabo@crimson.ua.edu</a>]<br>
Gesendet: Sonntag, 22. Oktober 2017 00:42<br>
An: <a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
Betreff: Re: [Wien] why we need to divide plasma frequency by root(2)?<br>
<br>
I could be wrong, but I think you just need to take the sqrt(2) of each plasma frequency:<br>
<br>
2.074/sqrt(2) 2.074/sqrt(2) 1.264/sqrt(2)<br>
<br>
then put them in case.inkram (not case.injoint).<br>
<br>
You should be able to use multiple plasma frequencies in case.inkram but don't forget to add a Gamma for drude value for each of the plasma frequencies [ <a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg14495.html" rel="noreferrer" target="_blank">https://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/msg14495.html</a> ].<br>
<br>
On 10/21/2017 4:29 PM, Dr. K. C. Bhamu wrote:<br>
<br>
---------------------case.<wbr>injoint<br>
xx<br>
xx<br>
1 >> correction<br>
2.0740 2.0740 1.2640<br>
xx<br>
------------------------------<wbr>-----<br>
<br>
were xx means no change in lines of case.inkram.<br>
<br>
<br>
<br>
At page numer 164 under section "8.11 KRAM (Kramers-Kronig transformation) in UG"<br>
<br>
<br>
"For a metal, the Plasma-frequencies (intraband transitions) for up and dn should be added, but then divided by√2, before using x kram."<br>
<br>
Does it mean that we need to add all three components and then divide by root(2) and put a single number [(2.074+2.074+1.264)/root(2)=<wbr>3.88] instead of all three in case.injoint?<br>
<br>
<br>
kind regards<br>
<br>
Bhamu<br>
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
</blockquote></div><br></div></div>