<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=utf-8">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">

<p>Dear Prof. Blaha and colleagues,,,</p>

<p><br>

</p>

<p>I 'm trying to reproducing the dielectric function of Gold as reported by <span>

<span itemtype="http://schema.org/Person" itemprop="author" class="nowrap"><span itemprop="name">Kathrin Glantschnig</span><sup>1,2,3</sup></span> and

<span itemtype="http://schema.org/Person" itemprop="author" class="nowrap"><span itemprop="name">Claudia Ambrosch-Drax in:</span></span></span></p>

<p><span><span itemtype="http://schema.org/Person" itemprop="author" class="nowrap"><span itemprop="name"><br>

</span></span></span></p>

<p><span><span itemtype="http://schema.org/Person" itemprop="author" class="nowrap"><span itemprop="name"><a href="http://iopscience.iop.org/article/10.1088/1367-2630/12/10/103048/meta" class="OWAAutoLink" id="LPlnk560647" previewremoved="true">http://iopscience.iop.org/article/10.1088/1367-2630/12/10/103048/meta</a><br>

</span></span></span></p>

<p><span><span itemtype="http://schema.org/Person" itemprop="author" class="nowrap"><span itemprop="name">and from which I quote:</span></span></span></p>

<p><span><span itemtype="http://schema.org/Person" itemprop="author" class="nowrap"><span itemprop="name"></span></span></span></p>

<div>"Relativistic effects also influence the plasma frequency changing from 9.44 to 8.96 eV when taking scalar-relativistic terms into account. The inclusion of SO interaction leads to an additional decrease to 8.81 eV. It now agrees well with a value of 8.83 eV

 extracted from experimental data by Cooper <em>et al</em> [<a class="cite hastip tooltipstered" href="http://iopscience.iop.org/article/10.1088/1367-2630/12/10/103048/meta#nj349811bib18" id="fnref-nj349811bib18">18</a>]."<br>

<br>

1- I use Wien2k 16.1 (after fixing the bug related to lapw2 in the DFT+U calculations) to calculate the optical constants of bulk Gold.<br>

2- I followed the user guide (i.e. increasing the maximum energy and setting Rkmax to a higher value ~11). and got the expected value of Wp when using using GGA alone.<br>

2- Upon inclusion of the spin-orbit coupling I got a value ~ 12.6 eV!!!. The inclusion of SOC is expected to reduce the value of Wp but here the opposite happened.

<br>

3- By switching off the SOC and I again got the expected value of Wp when GGA alone is used.<br>

<br>

My question:<br>

<br>

1- Why Wp increases when SOC is included in the optical calculations?. What is wrong<br>

<br>

- the lattice parameter I used is 4.1011 Angstrom.<br>

- a dense mesh of 41x41x41 i used for the optical calculations..<br>

<br>

<br>

O A Yassin<br>

<br>

<br>

</div>

<br>

<p></p>

<p><br>

</p>

<p><br>

</p>

<p><br>

</p>

<p><br>

</p>

<div id="Signature">

<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">

Sent from <a href="http://aka.ms/weboutlook" id="LPNoLP">Outlook</a></div>

</div>

</div>

</body>

</html>