[Wien] Wien2k OPTIC crashing for XMCD calculation
Santu Baidya
santubaidya2009 at gmail.com
Fri May 15 09:01:11 CEST 2015
Dear Prof. Blaha,Gavin Abo andXavier Rocquefelte*,*
Thank you for your suggestions. As you see from my previous
mail I first did GGA+U+SO calculation for XMCD of Co L23 edge.
It did not work and I moved to GGA+U cal.
I first did scf calculation under GGA+U+SO. then ran these commands as i
see from forum.....
a) runsp -so -dm -orb -c -s lapw1 -e lcore
b) x kgen -so (for dense grid)
c) x lapw1 -up/dn -orb -p(for eigenvalues)
d) Edit case.inop with xtra line "XMCD 1 L23" as it is first atom
e) x optic -so -up
Then the problem comes with error.......
running OPTIC in parallel mode
[1] 24537
forrtl: severe (64): input conversion error, unit 35, file
/users/santu/ProjectCoxOy/AFM-A-CoO/GGA+U+SO/xmcd/try/CoO/CoO.corewfup
Image PC Routine Line Source
opticc 00000000004E9ABD Unknown Unknown Unknown
opticc 00000000004E85C5 Unknown Unknown Unknown
opticc 0000000000492599 Unknown Unknown Unknown
opticc 0000000000451DCA Unknown Unknown Unknown
opticc 00000000004515C0 Unknown Unknown Unknown
opticc 000000000047126B Unknown Unknown Unknown
opticc 000000000046E60E Unknown Unknown Unknown
opticc 000000000043BEFA cor_mat_ 220
sph-UPcor_tmp.f
opticc 000000000041D3A5 MAIN__ 460 opmain.f
opticc 00000000004036AC Unknown Unknown Unknown
libc.so.6 00002AF40916AC36 Unknown Unknown Unknown
opticc 0000000000403589 Unknown Unknown Unknown
[1] + Done ( cd $PWD; $t $exe ${def}_${loop}.def;
rm -f .lock_$lockfile[$p] ) >> ...
[1] 24543
OPTIC - ERROR
[1] + Done ( cd $PWD; $t $exe ${def}_${loop}.def;
rm -f .lock_$lockfile[$p] ) >> ...
** OPTIC crashed!
0.036u 0.024s 0:05.26 0.9% 0+0k 0+1016io 0pf+0w
error: command /home/santu/wien2k/14/opticcpara -up -c -so upoptic.def
failed
That is why I moved to normal GGA+U calculation and got again problem in
optic and posted in forum.
I would like to mention here when I did not add line "XMCD 1 L23" in
case.inop optic program ran without any problem.
So please suggest me to know the problem here and solve it.
Thankin you,
Santu Baidya
On 15 May 2015 at 11:59, Peter Blaha <pblaha at theochem.tuwien.ac.at> wrote:
> You need spin-orbit coupling for xmcd
>
> Am 14.05.2015 um 20:58 schrieb Santu Baidya:
>
>> Dear Prof. Blaha and Wien2k users,
>>
>> I am using wien2k version WIEN2k_14.1 (Release 8/9/2014) for
>> calculating xmcd spectra of Co L23 edge for 8 atoms system. So I generated
>> case.struct file. The did normal
>> GGA+U spin polarized calculation using "runsp_lapw -orb -p".
>>
>> Then I ran few commands as mentioned in userguide and wien2k forum:
>> 1) x kgen -p (for denser grid)
>> 2) x lapw1 -orb -up/dn -p (for eigenvalue calculations)
>> 3) x lapw2 -fermi -up/dn -p (for fermi energy calculation)
>> 4) editing case.inop file with additional line "XMCD 1 L23"
>> 5) x optic -orb -up -p
>>
>> Before optic it ran very well. But after running optic it crashes with
>> error:
>>
>> running OPTIC in parallel mode
>> [1] 13265
>> [1] + Done ( cd $PWD; $t $exe
>> ${def}_${loop}.def; rm -f .lock_$lockfile[$p] ) >> ...
>> [1] 13273
>> [1] + Done ( cd $PWD; $t $exe
>> ${def}_${loop}.def; rm -f .lock_$lockfile[$p] ) >> ...
>> [1] 13281
>> [1] + Done ( cd $PWD; $t $exe
>> ${def}_${loop}.def; rm -f .lock_$lockfile[$p] ) >> ...
>> [1] 13289
>> [1] + Done ( cd $PWD; $t $exe
>> ${def}_${loop}.def; rm -f .lock_$lockfile[$p] ) >> ...
>> [1] 13295
>> [1] + Done ( cd $PWD; $t $exe
>> ${def}_${loop}.def; rm -f .lock_$lockfile[$p] ) >> ...
>> ** OPTIC crashed!
>> 0.016u 0.032s 0:05.26 0.7%0+0k 0+840io 0pf+0w
>> error: command /home/santu/wien2k/14/opticcpara -up -c upoptic.def
>> failed
>>
>> After looking into other such error in wen2k forum I do not find its
>> solution which led me to report it in the forum and asking for solution if
>> known. My system has only 8
>> atoms which is not big. So why should it crash !!!!
>> Even I tried in GGA+U+SO calculation I got same problem which led me to
>> do normal GGA+U calculation for trial. But it gives me same "optic crash".
>>
>>
>> I would like to request you to help me to solve this problem in running
>> optic program. Is wien2k optic program is unstable for XMCD calculation ?!
>>
>>
>> Thanking you in advance,
>>
>> Santu Baidya
>> University of Duisburg
>> Germany
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> Wien mailing list
>> Wien at 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
>>
>>
> --
> -----------------------------------------
> Peter Blaha
> Inst. Materials Chemistry, TU Vienna
> Getreidemarkt 9, A-1060 Vienna, Austria
> Tel: +43-1-5880115671
> Fax: +43-1-5880115698
> email: pblaha at theochem.tuwien.ac.at
> -----------------------------------------
>
> _______________________________________________
> Wien mailing list
> Wien at 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20150515/d30cbf90/attachment.html>
More information about the Wien
mailing list