<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Thank you for the update</p>
<p><br>
</p>
<p><br>
</p>
<div id="x_Signature">
<div id="x_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>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Wien <wien-bounces@zeus.theochem.tuwien.ac.at> on behalf of Peter Blaha <pblaha@theochem.tuwien.ac.at><br>
<b>Sent:</b> Friday, May 5, 2017 3:28:28 PM<br>
<b>To:</b> A Mailing list for WIEN2k users<br>
<b>Subject:</b> Re: [Wien] Non-spin ploarized and spin polarized SOC for Gold</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">I have fixed the normalization problems for a non-spinpolarized SO
<br>
calculation of the optical properties.<br>
(It was just a factor of 2 in case.joint and a factor of sqrt(2) in the <br>
plasma frequency).<br>
<br>
The modified files<br>
<br>
opmain.f (SRC_optic) and<br>
<br>
readop.f, read_diag.f, joint.f (SRC_joint) are attached.<br>
<br>
Regards<br>
<br>
On 05/04/2017 08:14 AM, Osama Yassin wrote:<br>
> Dear Prof Blaha<br>
><br>
> Guten Morgen<br>
><br>
><br>
> Thank you for your reply. Before I answer your questions I have just<br>
> went back to Wien2k 14.2 and run non-spin polarized OPTIC+SOC. It gave<br>
> me what I expect. I followed the same procedure as I did with Wien2k 16.1.<br>
><br>
><br>
> Now I have doubt about how Wien2k 16.1 deals with non-spin polarized<br>
> "optic"+soc calculation. It may be related to how joint behaves or<br>
> something else!!.<br>
><br>
><br>
> May you please test what I said.<br>
><br>
><br>
> Below I answer your questions*:*<br>
><br>
> ------------------------------------------------------------------------<br>
><br>
><br>
> I'd expect that a spin-polarized Au calculation converges to a zero<br>
> moment solution, which should be identical to the non-spinpolarized one.<br>
><br>
> Do the spin-polarized+so and non-spinpolarized+so calculations give the<br>
> same DOS ??<br>
><br>
> My answer:<br>
><br>
> I did not do DOS calculations because my intention is to reproduce the<br>
> dielectric function of gold. I will check it.<br>
><br>
> Did you allow for symmetry break in the spin-polarized SO calculation ?.<br>
><br>
> Yes I did.<br>
><br>
> What are the results for spin-polarized and non-spinpolarized<br>
> calculations without SO ??.<br>
><br>
> Both work well and gave me the value of the plasma frequency within the<br>
> expected range for the case without soc.<br>
><br>
> Best regards,,<br>
><br>
> Am 03.05.2017 um 19:57 schrieb Osama Yassin:<br>
>> Dear colleagues<br>
>><br>
>><br>
>> With reference to my earlier post, which was concerned with the optical<br>
>> properties of gold with spin-orbit coupling, I have a query that I hope<br>
>> to find an answer for it.<br>
>><br>
>><br>
>> To calculate the dielectric properties of Gold (Au) I used spin<br>
>> polarized calculation and everything went fine. However, the non-spin<br>
>> polarized calculation did not give the same results as that of the spin<br>
>> polarized one. The later gives the expected values for the plasma<br>
>> frequency ~ 8.7 eV.<br>
>><br>
>><br>
>> My query what is wrong with the non-spin polarized Spin-orbit<br>
>> calculations when following the steps:<br>
>><br>
>><br>
>> 1. Generate the struct file with lattice constant 4.1011A (optimized)<br>
>> and space group No. 225 .<br>
>><br>
>><br>
>> 2- Set the Kmax to 9 and Gmax to 16.<br>
>><br>
>><br>
>> 3- Make dense k-mesh (31x31x31).<br>
>><br>
>> 4- Initialize the spin orbit calculation and set Emax to 4.5 ry.<br>
>><br>
>><br>
>> 5- The inso is :<br>
>><br>
>> WFFIL<br>
>> 4 0 0 llmax,ipr,kpot<br>
>> -10 1.5 Emin, Emax<br>
>> 0 0 1 h,k,l (direction of magnetization)<br>
>> 0 number of atoms with RLO<br>
>> 0 0 number of atoms without SO, atomnumbers<br>
>><br>
>><br>
>> Then the scf converged successfully.<br>
>><br>
>><br>
>> 6- When executing the OPTIC the following steps are followed<br>
>><br>
>><br>
>> 99999 1 number of k-points, first k-point<br>
>> -5.0 3.0 9999 Emin, Emax for matrix elements, NBvalMAX<br>
>> 1 number of choices (columns in *outmat): 2: hex or<br>
> tetrag. case<br>
>> 3 Re zz<br>
>> OFF ON/OFF writes MME to unit 4<br>
>><br>
>><br>
>> and<br>
>><br>
>><br>
>> 1 9999 9999 : LOWER,UPPER and (optional) UPPER-VAL<br>
>> BANDINDEX<br>
>> 0.0000 0.00100 4.0000 : EMIN DE EMAX FOR ENERGYGRID IN ryd<br>
>> eV : output units eV / ryd / cm-1<br>
>> 6 : SWITCH<br>
>> 1 : NUMBER OF COLUMNS<br>
>> 0.02 0.1 0.3 : BROADENING (FOR DRUDE MODEL - switch<br>
> 6,7 -<br>
>> ONLY)<br>
>><br>
>> then the switch set to 4.<br>
>><br>
>><br>
>> The results of the plasma frequency is 12.25 ev.<br>
>><br>
>><br>
>> Look forward to hearing from you...<br>
>><br>
>><br>
>> Yassin<br>
>><br>
>><br>
>><br>
>> Sent from Outlook <<a href="http://aka.ms/weboutlook">http://aka.ms/weboutlook</a>><br>
> Check out Outlook.com – free, personal email from Microsoft.<br>
> <<a href="http://aka.ms/weboutlook">http://aka.ms/weboutlook</a>><br>
> aka.ms<br>
> Take your email anywhere you go when you add your free, personal,<br>
> Outlook.com webmail to your Android, iPhone, or Windows mobile devices.<br>
> Send and receive messages with mobile mail from Outlook.com<br>
><br>
><br>
>><br>
>><br>
>> _______________________________________________<br>
>> Wien mailing list<br>
>> Wien@zeus.theochem.tuwien.ac.at<br>
>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> Wien -- A Mailing list for WIEN2k users<br>
> <<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>><br>
> zeus.theochem.tuwien.ac.at<br>
> A Mailing list for WIEN2k users. Please post questions, suggestions or<br>
> comments about WIEN2k ONLY in this list. Please follow the following<br>
> "Nettiquette" (depending ...<br>
><br>
><br>
>> SEARCH the MAILING-LIST at:<br>
> <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> wien - The Mail Archive<br>
> <<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a>><br>
> <a href="http://www.mail-archive.com">www.mail-archive.com</a><br>
> Messages by Thread [Wien] How to know the user id & password for<br>
> graphical interface of wien2k shamik chakrabarti. Re: [Wien] How to know<br>
> the user id & password for ...<br>
><br>
><br>
>><br>
><br>
> --<br>
> --------------------------------------------------------------------------<br>
> Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
> Phone: +43-1-58801-165300 FAX: +43-1-58801-165982<br>
> Email: blaha@theochem.tuwien.ac.at WIEN2k: <a href="http://www.wien2k.at">http://www.wien2k.at</a><br>
> WWW:<br>
> <a href="http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------">
http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------</a><br>
><br>
> IMC : Prof. Dr. P. Blaha: Computational Materials Science - Home of<br>
> WIEN2k<br>
> <<a href="http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------">http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------</a>><br>
> <a href="http://www.imc.tuwien.ac.at">www.imc.tuwien.ac.at</a><br>
> Homepage of Institute of Materials Chemistry<br>
><br>
><br>
><br>
> _______________________________________________<br>
> Wien mailing list<br>
> Wien@zeus.theochem.tuwien.ac.at<br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> Wien -- A Mailing list for WIEN2k users<br>
> <<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>><br>
> zeus.theochem.tuwien.ac.at<br>
> A Mailing list for WIEN2k users. Please post questions, suggestions or<br>
> comments about WIEN2k ONLY in this list. Please follow the following<br>
> "Nettiquette" (depending ...<br>
><br>
><br>
> SEARCH the MAILING-LIST at:<br>
> <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> wien - The Mail Archive<br>
> <<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a>><br>
> <a href="http://www.mail-archive.com">www.mail-archive.com</a><br>
> Messages by Thread [Wien] How to know the user id & password for<br>
> graphical interface of wien2k shamik chakrabarti. Re: [Wien] How to know<br>
> the user id & password for ...<br>
><br>
><br>
><br>
><br>
> _______________________________________________<br>
> Wien mailing list<br>
> Wien@zeus.theochem.tuwien.ac.at<br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.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">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
><br>
<br>
-- <br>
<br>
P.Blaha<br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982<br>
Email: blaha@theochem.tuwien.ac.at WIEN2k: <a href="http://www.wien2k.at">http://www.wien2k.at</a><br>
WWW: <a href="http://www.imc.tuwien.ac.at/TC_Blaha">http://www.imc.tuwien.ac.at/TC_Blaha</a><br>
--------------------------------------------------------------------------<br>
</div>
</span></font>
</body>
</html>