<div dir="ltr"><div><div><div>Dear Prof. Peter, <span class="gmail-sender gmail-italic">Gerhard and Gavin</span></div><div><span class="gmail-sender gmail-italic"><br></span></div><div><span class="gmail-sender gmail-italic">The calculations I am doing is well optimised for a binary compound and bulk modules and all relevant output parameters are in the good acceptable range of the already reported values. So I think the initialized and scf parameters are well defined in my calculations. Now onward I doped the case with transition metal and optimised the doped structure to get band, doss and optical properties. <br></span></div><div><span class="gmail-sender gmail-italic">   <br></span></div><div><span class="gmail-sender gmail-italic">After doping I am getting metal-like behavior (for both spins) so there should be the role of w_pl in optical properties in intraband transitions. <br></span></div><div><span class="gmail-sender gmail-italic"><br></span></div><div><span class="gmail-sender gmail-italic">The present case is regarding the optical properties of metals.</span></div><div><span class="gmail-sender gmail-italic"><br></span></div><div><span class="gmail-sender gmail-italic"><br></span></div><div><br></div><div>The formula we see in the case.outputjount(up/dn) is:<br><br>   !!! WARNING: w_pl = sqrt( w_pl^2(up-spin) + w_pl^2(dn-spin) )!!!<br><br><br></div>Just below this the  program prints        w_p_<span style="color:rgb(0,0,255)">xx</span>    in [eV]<div><br></div><div>in one of my case, I am getting this value. <br></div><div><br></div><div>       1.7843<br></div><div><br></div><div>So from my previous queries and reply from Prof. <span class="gmail-sender gmail-italic">Gerhard</span> and Gavin, what I understand is;</div><div><br></div><div><br></div><div>If the case is -sp(so) [preset case] then, w_pl  should be sqrt( w_pl_<span style="color:rgb(0,0,255)">xx</span>^2(up-spin) + w_pl_<span style="color:rgb(0,0,255)">xx</span>^2(dn-spin) )</div><div><br></div><div><br></div><div>So the value printed above [1.7843] is outputted as w_pl_<span style="color:rgb(0,0,255)">xx</span>^2(up-spin) by the program and same for -dn spin where the outputted w_pl_xx 2.0058 [eV] in case.outputjoinudn is printed as w_pl_<span style="color:rgb(0,0,255)">xx</span>^2(dn-spin).</div><div><br></div><div><br></div><div>Now to get the final w_pl we need to do a small mathematical operation</div><div><br></div><div><br></div><div>w_pl =sqrt( 1.7843^2 + 2.0058^2 )= 2.6845 eV.</div><div><br></div><div>Now we set this value for kram and proceed for the calculation with the interband transition (option 4).</div><div><br></div><div><br></div><div>Please make any comment if I am doing any mistake in the procedure or it is the correct one.  <br></div><div><br></div><div>Prof. Peter also suggested putting <a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg14472.html">TETRA 101 for metal</a> cases with the additional statement that "It will not matter too much". So I did all calculations with TETRA 0.000 only.</div><div><br></div><br></div>regards<br><br></div>Bhamu<br></div>