I think in case.injoint number of columns should be 7, base on your input in case.inop. <br><br><div id="ymail_android_signature"><a id="ymail_android_signature_link" href="https://go.onelink.me/107872968?pid=InProduct&c=Global_Internal_YGrowth_AndroidEmailSig__AndroidUsers&af_wl=ym&af_sub1=Internal&af_sub2=Global_YGrowth&af_sub3=EmailSignature">Sent from Yahoo Mail on Android</a></div> <br> <blockquote style="margin: 0 0 20px 0;"> <div style="font-family:Roboto, sans-serif; color:#6D00F6;"> <div>On Sat, May 4, 2019 at 12:15 PM, tran@theochem.tuwien.ac.at</div><div><tran@theochem.tuwien.ac.at> wrote:</div> </div> <div style="padding: 10px 0 0 20px; margin: 10px 0 0 0; border-left: 1px solid #6D00F6;"> Hi,<br clear="none"><br clear="none">To which extent do the DFT and experimental spectra differ?<br clear="none">Is there still some very vague qualitative resemblance?<br clear="none"><br clear="none">In practice DFT is an approximate method.<br clear="none">That means that disagreement with experiment is<br clear="none">always expected.<br clear="none"><br clear="none">F. Tran<br clear="none"><br clear="none">On Saturday 2019-05-04 12:30, KRISH wrote:<br clear="none"><br clear="none">>Date: Sat, 4 May 2019 12:30:12<br clear="none">>From: KRISH <<a shape="rect" ymailto="mailto:krish.odc@gmail.com" href="mailto:krish.odc@gmail.com">krish.odc@gmail.com</a>><br clear="none">>Reply-To: A Mailing list for WIEN2k users <<a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br clear="none">>To: <a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a><br clear="none">>Subject: [Wien] XMCD spectra<div class="yqt1355163137 yQTDBase" id="yqtfd58381"><br clear="none">><br clear="none">>Dear Wien2k Users,<br clear="none">>I am new to Wien2k. I am using WIEN2k 18.2 in Linux. I tried to calculated XMCD spectra for US compound. <br clear="none">>I have used the following steps to calculate the XMCD spectra for US (NaCl structure)<br clear="none">><br clear="none">>1. cif2cell us.cif<br clear="none">>2. init_lapw<br clear="none">>3. runsp_lapw<br clear="none">>4. delete brodyn files<br clear="none">>5. initso_lapw<br clear="none">>       Magnetic axis 1 1 1<br clear="none">>       No RLOs were used (Since optics does not work well with RLOs)<br clear="none">>6. us.in2 -TOT was changed to FERMI<br clear="none">>7. The system is metal -TETRA 101 inserted to us.in2 file<br clear="none">>8. IPRINT option was changed from 0 to 1 in 'us.inc' and 'us.inc_so' files<br clear="none">>9. Created us.indm, us.indmc, and us.inorb files (LSDA+U)<br clear="none">>10. runsp_lapw -so -orb -dm<br clear="none">>11. Using w2web -selected <x lapw2 -fermi –up><br clear="none">>12. Edited the us.inop file<br clear="none">><br clear="none">>us.inop<br clear="none">><br clear="none">>9999 1       number of k-points, first k-point<br clear="none">>-10.0 10.0 9999 Emin, Emax for matrix elements, NBvalMAX<br clear="none">>XMCD 2 M45      optional line: for XMCD of 2nd atom and M45 spectrum<br clear="none">>7             number of choices (columns in *outmat): 2: hex or tetrag. case<br clear="none">>1             Re xx<br clear="none">>2             Re yy<br clear="none">>3             Re zz<br clear="none">>4             Re xy<br clear="none">>7             Im xy<br clear="none">>8             Im xz<br clear="none">>9             Im yz<br clear="none">>OFF           ON/OFF   writes MME to unit 4 <br clear="none">><br clear="none">>13.  Using w2web -selected <x optic -up><br clear="none">>14. Edited us.injoint<br clear="none">><br clear="none">>us.injoint<br clear="none">>  <br clear="none">>1  9999 9999              : LOWER,UPPER and (optional) UPPER-VAL BANDINDEX<br clear="none">>-0.0000    0.00100   1.0000 : EMIN DE EMAX FOR ENERGYGRID IN ryd<br clear="none">>eV                            : output units  eV / ryd  / cm-1<br clear="none">>XMCD                          : omitt these 4 lines for non-XMCD<br clear="none">>-255.34 -268.49               : core energies in Ry (grep :3D case.scfcup) 3d and 3d*(scfcup)<br clear="none">>1.6  0.6                      : core-hole broadening (eV) for both core states<br clear="none">>0.1                           : spectrometer broadening (eV)<br clear="none">>     6                        : SWITCH<br clear="none">>     4                        : SWITCH<br clear="none">>     2                        : NUMBER OF COLUMNS<br clear="none">>   0.1  0.1  0.3              : BROADENING (FOR DRUDE MODEL - switch 6,7 -<br clear="none">>ONLY)<br clear="none">><br clear="none">>SWITCH:<br clear="none">><br clear="none">>   0...JOINTDOS FOR EACH BAND COMBINATION<br clear="none">>   1...JOINTDOS AS SUM OVER ALL BAND COMBINATIONS<br clear="none">>   2...DOS FOR EACH BAND<br clear="none">>   3...DOS AS SUM OVER ALL BANDS<br clear="none">>   4...Im(EPSILON)<br clear="none">>   5...Im(EPSILON) for each band combination<br clear="none">>   6...INTRABAND contributions<br clear="none">>   7...INTRABAND contributions including band analysis<br clear="none">><br clear="none">>15. Using w2web -selected <x joint -up><br clear="none">><br clear="none">> I obtained an XMCD (M4,5) spectra which does not match with the experiment. <br clear="none">>Am I missing anything here? I do not understand why the calculated XMCD spectra does not match with experiment. <br clear="none">>Any help will be highly appreciated. <br clear="none">><br clear="none">><br clear="none">>Thank you.<br clear="none">>KRISH<br clear="none">><br clear="none">><br clear="none">><br clear="none">><br clear="none">></div> </div> </blockquote>