Re: [Wien] XMCD spectra
I think in case.injoint number of columns should be 7, base on your input in case.inop. Sent from Yahoo Mail on Android On Sat, May 4, 2019 at 12:15 PM, t...@theochem.tuwien.ac.at wrote: Hi, To which extent do the DFT and experimental spectra differ? Is there still some very vague qualitative resemblance? In practice DFT is an approximate method. That means that disagreement with experiment is always expected. F. Tran On Saturday 2019-05-04 12:30, KRISH wrote: >Date: Sat, 4 May 2019 12:30:12 >From: KRISH >Reply-To: A Mailing list for WIEN2k users >To: wien@zeus.theochem.tuwien.ac.at >Subject: [Wien] XMCD spectra > >Dear Wien2k Users, >I am new to Wien2k. I am using WIEN2k 18.2 in Linux. I tried to calculated >XMCD spectra for US compound. >I have used the following steps to calculate the XMCD spectra for US (NaCl >structure) > >1. cif2cell us.cif >2. init_lapw >3. runsp_lapw >4. delete brodyn files >5. initso_lapw > Magnetic axis 1 1 1 > No RLOs were used (Since optics does not work well with RLOs) >6. us.in2 -TOT was changed to FERMI >7. The system is metal -TETRA 101 inserted to us.in2 file >8. IPRINT option was changed from 0 to 1 in 'us.inc' and 'us.inc_so' files >9. Created us.indm, us.indmc, and us.inorb files (LSDA+U) >10. runsp_lapw -so -orb -dm >11. Using w2web -selected >12. Edited the us.inop file > >us.inop > > 1 number of k-points, first k-point >-10.0 10.0 Emin, Emax for matrix elements, NBvalMAX >XMCD 2 M45 optional line: for XMCD of 2nd atom and M45 spectrum >7 number of choices (columns in *outmat): 2: hex or tetrag. case >1 Re xx >2 Re yy >3 Re zz >4 Re xy >7 Im xy >8 Im xz >9 Im yz >OFF ON/OFF writes MME to unit 4 > >13. Using w2web -selected >14. Edited us.injoint > >us.injoint > >1 : LOWER,UPPER and (optional) UPPER-VAL BANDINDEX >-0. 0.00100 1. : EMIN DE EMAX FOR ENERGYGRID IN ryd >eV : output units eV / ryd / cm-1 >XMCD : omitt these 4 lines for non-XMCD >-255.34 -268.49 : core energies in Ry (grep :3D case.scfcup) 3d >and 3d*(scfcup) >1.6 0.6 : core-hole broadening (eV) for both core states >0.1 : spectrometer broadening (eV) > 6 : SWITCH > 4 : SWITCH > 2 : NUMBER OF COLUMNS > 0.1 0.1 0.3 : BROADENING (FOR DRUDE MODEL - switch 6,7 - >ONLY) > >SWITCH: > > 0...JOINTDOS FOR EACH BAND COMBINATION > 1...JOINTDOS AS SUM OVER ALL BAND COMBINATIONS > 2...DOS FOR EACH BAND > 3...DOS AS SUM OVER ALL BANDS > 4...Im(EPSILON) > 5...Im(EPSILON) for each band combination > 6...INTRABAND contributions > 7...INTRABAND contributions including band analysis > >15. Using w2web -selected > > I obtained an XMCD (M4,5) spectra which does not match with the experiment. >Am I missing anything here? I do not understand why the calculated XMCD >spectra does not match with experiment. >Any help will be highly appreciated. > > >Thank you. >KRISH > > > > > ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
Re: [Wien] XMCD spectra
Hi, To which extent do the DFT and experimental spectra differ? Is there still some very vague qualitative resemblance? In practice DFT is an approximate method. That means that disagreement with experiment is always expected. F. Tran On Saturday 2019-05-04 12:30, KRISH wrote: Date: Sat, 4 May 2019 12:30:12 From: KRISH Reply-To: A Mailing list for WIEN2k users To: wien@zeus.theochem.tuwien.ac.at Subject: [Wien] XMCD spectra Dear Wien2k Users, I am new to Wien2k. I am using WIEN2k 18.2 in Linux. I tried to calculated XMCD spectra for US compound. I have used the following steps to calculate the XMCD spectra for US (NaCl structure) 1. cif2cell us.cif 2. init_lapw 3. runsp_lapw 4. delete brodyn files 5. initso_lapw Magnetic axis 1 1 1 No RLOs were used (Since optics does not work well with RLOs) 6. us.in2 -TOT was changed to FERMI 7. The system is metal -TETRA 101 inserted to us.in2 file 8. IPRINT option was changed from 0 to 1 in 'us.inc' and 'us.inc_so' files 9. Created us.indm, us.indmc, and us.inorb files (LSDA+U) 10. runsp_lapw -so -orb -dm 11. Using w2web -selected 12. Edited the us.inop file us.inop 1 number of k-points, first k-point -10.0 10.0 Emin, Emax for matrix elements, NBvalMAX XMCD 2 M45 optional line: for XMCD of 2nd atom and M45 spectrum 7 number of choices (columns in *outmat): 2: hex or tetrag. case 1 Re xx 2 Re yy 3 Re zz 4 Re xy 7 Im xy 8 Im xz 9 Im yz OFF ON/OFF writes MME to unit 4 13. Using w2web -selected 14. Edited us.injoint us.injoint 1 : LOWER,UPPER and (optional) UPPER-VAL BANDINDEX -0. 0.00100 1. : EMIN DE EMAX FOR ENERGYGRID IN ryd eV : output units eV / ryd / cm-1 XMCD : omitt these 4 lines for non-XMCD -255.34 -268.49 : core energies in Ry (grep :3D case.scfcup) 3d and 3d*(scfcup) 1.6 0.6 : core-hole broadening (eV) for both core states 0.1 : spectrometer broadening (eV) 6 : SWITCH 4 : SWITCH 2 : NUMBER OF COLUMNS 0.1 0.1 0.3 : BROADENING (FOR DRUDE MODEL - switch 6,7 - ONLY) SWITCH: 0...JOINTDOS FOR EACH BAND COMBINATION 1...JOINTDOS AS SUM OVER ALL BAND COMBINATIONS 2...DOS FOR EACH BAND 3...DOS AS SUM OVER ALL BANDS 4...Im(EPSILON) 5...Im(EPSILON) for each band combination 6...INTRABAND contributions 7...INTRABAND contributions including band analysis 15. Using w2web -selected I obtained an XMCD (M4,5) spectra which does not match with the experiment. Am I missing anything here? I do not understand why the calculated XMCD spectra does not match with experiment. Any help will be highly appreciated. Thank you. KRISH ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
[Wien] XMCD spectra
Dear Wien2k Users, I am new to Wien2k. I am using WIEN2k 18.2 in Linux. I tried to calculated XMCD spectra for US compound. I have used the following steps to calculate the XMCD spectra for US (NaCl structure) 1. cif2cell us.cif 2. init_lapw 3. runsp_lapw 4. delete brodyn files 5. initso_lapw Magnetic axis 1 1 1 No RLOs were used (Since optics does not work well with RLOs) 6. us.in2 -TOT was changed to FERMI 7. The system is metal -TETRA 101 inserted to us.in2 file 8. IPRINT option was changed from 0 to 1 in 'us.inc' and 'us.inc_so' files 9. Created us.indm, us.indmc, and us.inorb files (LSDA+U) 10. runsp_lapw -so -orb -dm 11. Using w2web -selected 12. Edited the us.inop file *us.inop* * 1 number of k-points, first k-point-10.0 10.0 Emin, Emax for matrix elements, NBvalMAXXMCD 2 M45 optional line: for XMCD of 2nd atom and M45 spectrum7 number of choices (columns in *outmat): 2: hex or tetrag. case1 Re xx2 Re yy3 Re zz4 Re xy7 Im xy8 Im xz9 Im yzOFF ON/OFF writes MME to unit 4 * 13. Using w2web -selected 14. Edited us.injoint *us.injoint* *1 : LOWER,UPPER and (optional) UPPER-VAL BANDINDEX-0.0.00100 1. : EMIN DE EMAX FOR ENERGYGRID IN rydeV: output units eV / ryd / cm-1XMCD : omitt these 4 lines for non-XMCD-255.34 -268.49 : core energies in Ry (grep :3D case.scfcup) 3d and 3d*(scfcup)1.6 0.6 : core-hole broadening (eV) for both core states0.1 : spectrometer broadening (eV) 6: SWITCH 4: SWITCH 2: NUMBER OF COLUMNS 0.1 0.1 0.3 : BROADENING (FOR DRUDE MODEL - switch 6,7 -ONLY)SWITCH: 0...JOINTDOS FOR EACH BAND COMBINATION 1...JOINTDOS AS SUM OVER ALL BAND COMBINATIONS 2...DOS FOR EACH BAND 3...DOS AS SUM OVER ALL BANDS 4...Im(EPSILON) 5...Im(EPSILON) for each band combination 6...INTRABAND contributions 7...INTRABAND contributions including band analysis* 15. Using w2web -selected I obtained an XMCD (M4,5) spectra which does not match with the experiment. Am I missing anything here? I do not understand why the calculated XMCD spectra does not match with experiment. Any help will be highly appreciated. Thank you. KRISH ___ Wien mailing list Wien@zeus.theochem.tuwien.ac.at http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
