It is relatively easy but there are many little details and you have to do it yourself;
the file case.int has the atom or orbital/atom that you want to plot, the DOS of these are put in case.dos1ev, case.dos2ev (for non sp-cases) the 'case.dos1ev' has the first 7 cases of the 'case.int' 'case.dos2ev' has from 8 to 14, etc. it also creates the :dos2 file, this last file you can rename to Graf-DOS and you now can add the cases in 'case.dos*ev', for example; plot "BiRuO-U-afm.dos1evup" using 1:2 title "Total" w l lt 1 lw 2, \ "BiRuO-U-afm.dos1evup" using 1:5 title "Ru" w l lt 3 lw 2, \ "BiRuO-U-afm.dos1evup" using 1:6 title "Ru" w l lt 4 lw 2, \ "BiRuO-U-afm.dos2evup" using 1:($2+$3+$4) title "O1" w l lt 8 lw 2, \ "BiRuO-U-afm.dos1evup" using 1:($3+$4) title "Bi" w l lt 2 lw 2 pause -1 Now case.int is: Title -0.50 0.002 1.500 0.003 # EMIN, DE, EMAX, Gauss-broadening(>de) 11 N 0.000 # NUMBER OF DOS-CASES specified below, G/L/B broadeni ng (Ry) 0 1 total # atom, case=column in qtl-header, label ***dos1*** 1 1 Bi 2 1 Bi 3 1 Ru 8 1 Ru 3 1 Ru 8 1 Ru ***dos1*** 4 1 O1 ***dos2*** 5 1 O1 6 1 O1 7 1 O2 1 1 Atom1 p 2 1 Atom2 tot ***dos2*** Here O1 represents atoms 4, 5 and 6 and ends up in "BiRuO-U-afm.dos2evup" using 1:($2+$3+$4) means add columns 2, 3 and 4, the '$' is put to distinguish from a number. You can put ($2+3), which means add '3' to column '2' By the way, you plot with gnuplot Graf-DOS Hope to be clear enough!!! (but I doubt it) ________________________________ De: wien-boun...@zeus.theochem.tuwien.ac.at <wien-boun...@zeus.theochem.tuwien.ac.at> en nombre de Farshad Nejadsattari <fneja...@uottawa.ca> Enviado: miƩrcoles, 16 de julio de 2014 08:11 p.m. Para: wien@zeus.theochem.tuwien.ac.at Asunto: [Wien] specific DOS plots Dear Wien2k community, I am working on a compound with 20 independent atoms, 5 of which are iron and 15 Aluminium, my question is how to plot the density of states resultant from the contribution of all 5 iron atoms and not each of them individually. As it appears the options available will only allow me to calculate the total DOS(contributions from all the atoms) or the total DOS of each individual atom and not for instance the collective contribution of 5 of them. below is my input file configuration from w2web: please specify which PDOS you want to calculate You can specify: total (for plotting 'Total DOS') N (to select atom N) tot,s,p,d,... (to select a set of PDOS for previously selected atom N) The following PDOS is possible: (Header from Al13Fe4.qtl ): ATOM Fe1: 1 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Fe2: 2 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Fe3: 3 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Fe4: 4 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Fe5: 5 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al1: 6 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al2: 7 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al3: 8 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al4: 9 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al5: 10 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al6: 11 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al7: 12 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al8: 13 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al9: 14 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al10: 15 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al11: 16 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al12: 17 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al13: 18 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al14: 19 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f ATOM Al15: 20 tot,s,p,PX,PY,PZ,d,DZ2,DX2Y2,DXY,DXZ,DYZ,f can you please inform me on how to plot the density of states arising only from all 5 iron atoms and excluding the contributions from the aluminium atoms such that I can produce a plot for density of states of "Iron". your help and assistance is truly appreciated. with regards, Farshad Nejadsattari
_______________________________________________ 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