24小时热门版块排行榜    

CyRhmU.jpeg
查看: 15218  |  回复: 88
【奖励】 本帖被评价49次,作者灯塔守望者增加金币 40.9

[资源] 【原创】vaspkit正在更新中

相信很多人都为vasp计算结果的处理感到很头痛,网上虽然有很多处理vasp计算结果的软件或者脚本,不过要么只针对于某一个问题,要么使用不太友好。为此,我弄了一个软件包 (vaspkit),把结构可视化,提取态密度,能带,电荷,自旋密度等数据等功能整合到一个软件包里,方便大家的使用。由于时间有限,我只花了两天时间做得,所以vaspkit工具包比较粗糙,而且有很多bug, 目前只实现最基本的功能。由于vaspkit 目前处于测试阶段,要不断的更新,所以不上传到论坛了。感兴趣的朋友可以到http://abinitio.workstag.com/topics/185下载。希望大家能提出宝贵的意见以及需要增加的功能。同时,我们希望对vaspkit感兴趣的朋友一起加入并不断地增加功能,让vaspkit能成为一个很强大的数据处理工具。在此,特别感谢论坛版主csfn, foxhunter, 语过添情以及蓝等朋友的支持。

当前版本:0.14 (2009-12-3).
1,修正若干bug;
2,增加输出xsf格式结构文件;
3,增加输出DOS积分;
4,增加提取指定原子指定轨道(s,p,d)的能带数据.


主要功能有:
1,把POSCAR 或者CONTCAR转化为cif或xsf格式。
2,提取电荷,自旋密度,使用vesta软件可视化。
3,提取总,投影态密度数据;
4,提取总,投影能带数据。


期待增加功能:
1,增加计算差分电荷数据;
2,xsf格式电荷密度。

文件输出格式:
1,结构可视化文件输出文件为poscar.cif 或contcar.cif, 用MS, VESTA或其它可视化软件软件打开;
2,总态密度输出文件为tdos.dat, 偏态密度文件为pdos.dat, 使用origin或gnuplot绘制。
3,能带数据输出文件名为band-x.dat (x=up, dn). 同样使用origin或gnuplot绘制。
4,  电荷和自旋密度输出文件名分别为chgden.dat 和spnden.dat, 用VESTA打开。注意同时把OUTCAR放进本地文件夹。


注意:最好在完整的vasp计算作业目录下运行,因为要用到OUTCAR读取元素符号,CHG读电荷密度,DOSCAR读取态密度,PROCAR读取能带数据,还有CONTCAR和POSCAR文件。请确保以上文件完整。

vasptip软件运行主页面为:
Choose the problem to solve:                        
=============== Structure options ===================
1: POSCAR                                            
2: CONTCAR                                          
                                                      
=============== DOS options ======================
11: Total DOS                                       
12: Projected DOS                                      
13: T2g-Eg DOS                           

=============== Energy Band options =================
21: Band                                             
22: Projected Band                     
                                                      
=============== Density options ====================
31: Charge Density                                   
32: Spin Density                                    
                                                      
=============== MISC options ======================
                                                      
0: Quit                                             
------------>>

[ Last edited by 灯塔守望者 on 2009-12-4 at 18:02 ]

记得09年我花了一周时间写了vaspkit, 那时候还对fortran不熟悉,之后由于忙碌的缘故,vaspkit再没有更新。一些同仁曾写email问我vaspkit是否有新版本,鉴于一些朋友的鼓励和我的需要,我打算近两个月更新一下vaspkit, 增加一些新的功能。借此抛砖引玉,希望更多对此有兴趣的同仁能够参与其中,vaspkit主页为http://vaspkit.sourceforge.net

vaspkit-0.20 预览版功能如下:
+---------------------------------------------------+
|>>>   Postprocesing tool for VASP                  |
|      Current version: 0.20 (18.Nov.2013)          |
|...................................................|
|>>>   Written by: Vei WANG (wangvei@live.com)      |
|>>>   Code URL: http://vaspkit.sourceforge.net     |
+---------------------------------------------------+

Choose the problem to solve:
=============== Structural options ==================
1: pos2cif or pos2xsf (unfinished)
2: Building supercell (unfinished)
3: EOS fitting

=============== Electronic options ==================
11: Total DOS
12: Projected DOS
21: Band Structure
22: Projected Band Structure
31: Charge Density
32: Charge Density Difference (unfinished)
33: Spin density
34: Spin Up/Down Density (unfinished)
41: Work Function

=============== Optical options =====================
51: Linear Optics

=============== MISC options ========================

0: Quit
------------>>

[ Last edited by 灯塔守望者 on 2013-11-18 at 21:39 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)


2楼2009-11-25 14:26:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★ ★
yjmaxpayne(金币+5,VIP+0):谢谢 11-25 21:22


[ Last edited by wuli8 on 2009-11-25 at 14:43 ]
3楼2009-11-25 14:41:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freer007

木虫 (正式写手)


强! 支持楼主。
4楼2009-11-25 14:59:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)



小木虫(金币+0.5):给个红包,谢谢回帖交流
At line 86 of file structinfo.f90
Fortran runtime error: Bad real number in item 4 of list input
5楼2009-11-25 15:36:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjmaxpayne(金币+1,VIP+0):谢谢 11-25 21:23
为什么不用脚本写呢?我觉得vts的那些脚本其实挺好的,修改野比较方便,而且不用编译。

其实我觉得自己写脚本比较可靠,fortran的文本处理能力有限,不如perl来得方便。等我脚本写的比较稳定了,再放上来

我现在在看楼主的代码
6楼2009-11-25 15:41:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

benhuzhou

铁杆木虫 (小有名气)


★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjmaxpayne(金币+1,VIP+0):谢谢 11-25 21:23
非常感谢!大力支持!
我用的是vasp5.2,我分析态密度的时候出现这样的错误,同样分析能带等其它的性质时,好像读的时候老是出错,可能是格式的问题。想问楼主写的程序主要是针对vasp哪个版本的。
At line 32 of file writetdos.f90
Fortran runtime error: Bad integer for item 5 in list input
7楼2009-11-25 16:00:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice_rain

铁杆木虫 (正式写手)



小木虫(金币+0.5):给个红包,谢谢回帖交流
支持,如果能做出Medea 的功能和效果就更好了
8楼2009-11-25 16:30:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★★★ 五星级,优秀推荐

good
9楼2009-11-25 16:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjmaxpayne(金币+1,VIP+0): 11-25 21:23
我使用的是vasp 4.6,为什么没有出现大家遇到的问题啊?你们再核查一下,谢谢。还有,如果你想以后使用fortran编计算程序的话,不推荐使用脚本。我最主要的目的是如何用fortran开发一个完整的计算软件包。支持各位热心的朋友,把你们的好东西也拿出来分享。

[ Last edited by 灯塔守望者 on 2009-11-25 at 20:21 ]
10楼2009-11-25 18:16:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)


★★★★★ 五星级,优秀推荐

鼓励楼主,并高度统一5楼的说法!
我也建议大家使用shell script 和awk进行数据处理!
大家可以一起开发!
11楼2009-11-25 18:43:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangdelin0000

木虫 (著名写手)


强烈支持师兄,
12楼2009-11-25 19:41:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★★★ 五星级,优秀推荐

恭喜 兄弟 终于小试牛刀啦 希望能继续奉献大作
13楼2009-11-25 22:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
我用ASE处理VASP,比较方便,而且ASE处理SIESTA,ABINIT就和VASP是一样的.
不过需要python语言基础
14楼2009-11-25 23:59:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangle_xq

银虫 (小有名气)


支持,回头发送测试结果。
不知道有没有考虑组织一个开发团队。我有兴趣参与。打造中国人自己的VASP处理软件。
引用回帖:
Originally posted by 灯塔守望者 at 2009-11-25 18:16:
我使用的是vasp 4.6,为什么没有出现大家遇到的问题啊?你们再核查一下,谢谢。还有,如果你想以后使用fortran编计算程序的话,不推荐使用脚本。我最主要的目的是如何用fortran开发一个完整的计算软件包。支持各 ...

15楼2009-11-26 00:43:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

这挺好,感谢灯塔兄辛苦工作并与大家分享成果!!也期望更多的朋友贡献更多精华资源!
16楼2009-11-26 01:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★★★ 五星级,优秀推荐

感谢楼主。非常好!!
17楼2009-11-26 02:19:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)


楼主辛苦了,感谢
18楼2009-11-26 12:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

benhuzhou

铁杆木虫 (小有名气)


我又试了很久,没有成功。和4楼出现的错误时类似的,不知道什么原因。
用的是vasp5.2,我分析态密度的时候出现这样的错误,同样分析能带等其它的性质时,好像读的时候老是出错,可能是格式的问题。想问楼主写的程序主要是针对vasp哪个版本的。
At line 32 of file writetdos.f90
Fortran runtime error: Bad integer for item 5 in list input,麻烦大家多多帮助,不知道是哪里格式有问题。
下面是DOSCAR这个文件的具体内容:

    1   1   1   0
  0.1697933E+02  0.2884996E-09  0.2884996E-09  0.2884996E-09  0.5000000E-15
  1.000000000000000E-004
  CAR
fcc Au                                 
     27.81778758     -6.32977860  301      6.82147560      1.00000000
     -6.330  0.0000E+00  0.0000E+00
     -6.216  0.0000E+00  0.0000E+00
     -6.102  0.0000E+00  0.0000E+00
     -5.988  0.0000E+00  0.0000E+00
     -5.874  0.0000E+00  0.0000E+00
     -5.761  0.0000E+00  0.0000E+00
     -5.647  0.0000E+00  0.0000E+00
     -5.533  0.0000E+00  0.0000E+00
     -5.419  0.2819E-14  0.3209E-15
     -5.305  0.1149E-11  0.1311E-12
     -5.192  0.2492E-09  0.2849E-10
     -5.078  0.2881E-07  0.3307E-08
     -4.964  0.1782E-05  0.2062E-06
     -4.850  0.5930E-04  0.6956E-05
     -4.736  0.1068E-02  0.1286E-03
     -4.622  0.1051E-01  0.1325E-02
     -4.509  0.5724E-01  0.7841E-02
     -4.395  0.1761E+00  0.2789E-01
     -4.281  0.3181E+00  0.6409E-01
     -4.167  0.3645E+00  0.1056E+00
     -4.053  0.3080E+00  0.1406E+00
     -3.939  0.2382E+00  0.1678E+00
     -3.826  0.1965E+00  0.1901E+00
     -3.712  0.1724E+00  0.2097E+00
     -3.598  0.1581E+00  0.2277E+00
     -3.484  0.1500E+00  0.2448E+00
     -3.370  0.1382E+00  0.2605E+00
     -3.256  0.1357E+00  0.2760E+00
     -3.143  0.1352E+00  0.2914E+00
     -3.029  0.1249E+00  0.3056E+00
     -2.915  0.1272E+00  0.3201E+00
     -2.801  0.1337E+00  0.3353E+00
     -2.687  0.1289E+00  0.3500E+00
     -2.574  0.1310E+00  0.3649E+00
     -2.460  0.1495E+00  0.3819E+00
     -2.346  0.1731E+00  0.4016E+00
     -2.232  0.2469E+00  0.4297E+00
     -2.118  0.4524E+00  0.4812E+00
     -2.004  0.7318E+00  0.5645E+00
     -1.891  0.9280E+00  0.6701E+00
     -1.777  0.1058E+01  0.7906E+00
     -1.663  0.1220E+01  0.9294E+00
     -1.549  0.1388E+01  0.1087E+01
     -1.435  0.1444E+01  0.1252E+01
     -1.321  0.1364E+01  0.1407E+01
     -1.208  0.1311E+01  0.1556E+01
     -1.094  0.1315E+01  0.1706E+01
     -0.980  0.1310E+01  0.1855E+01
     -0.866  0.1391E+01  0.2013E+01
     -0.752  0.1619E+01  0.2197E+01
     -0.639  0.1820E+01  0.2405E+01
     -0.525  0.1789E+01  0.2608E+01
     -0.411  0.1525E+01  0.2782E+01
     -0.297  0.1215E+01  0.2920E+01
     -0.183  0.1051E+01  0.3040E+01
     -0.069  0.1091E+01  0.3164E+01
      0.044  0.1275E+01  0.3309E+01
      0.158  0.1503E+01  0.3480E+01
      0.272  0.1886E+01  0.3695E+01
      0.386  0.2648E+01  0.3996E+01
      0.500  0.3439E+01  0.4388E+01
      0.614  0.3553E+01  0.4792E+01
      0.727  0.2971E+01  0.5130E+01
      0.841  0.2153E+01  0.5375E+01
      0.955  0.1450E+01  0.5540E+01
      1.069  0.1104E+01  0.5666E+01
      1.183  0.1089E+01  0.5790E+01
      1.297  0.1170E+01  0.5923E+01
      1.410  0.1151E+01  0.6054E+01
      1.524  0.1043E+01  0.6173E+01
      1.638  0.1059E+01  0.6293E+01
      1.752  0.1500E+01  0.6464E+01
      1.866  0.2424E+01  0.6740E+01
      1.979  0.3154E+01  0.7099E+01
      2.093  0.3033E+01  0.7445E+01
      2.207  0.2390E+01  0.7717E+01
      2.321  0.1787E+01  0.7920E+01
      2.435  0.1371E+01  0.8076E+01
      2.549  0.1178E+01  0.8210E+01
      2.662  0.1175E+01  0.8344E+01
      2.776  0.1293E+01  0.8491E+01
      2.890  0.1645E+01  0.8678E+01
      3.004  0.2333E+01  0.8944E+01
      3.118  0.2918E+01  0.9276E+01
      3.232  0.2902E+01  0.9606E+01
      3.345  0.2725E+01  0.9916E+01
      3.459  0.2707E+01  0.1022E+02
      3.573  0.2175E+01  0.1047E+02
      3.687  0.1106E+01  0.1060E+02
      3.801  0.3414E+00  0.1064E+02
      3.914  0.1083E+00  0.1065E+02
      4.028  0.9483E-01  0.1066E+02
      4.142  0.9400E-01  0.1067E+02
      4.256  0.7483E-01  0.1068E+02
      4.370  0.7383E-01  0.1069E+02
      4.484  0.1096E+00  0.1070E+02
      4.597  0.1611E+00  0.1072E+02
      4.711  0.1883E+00  0.1074E+02
      4.825  0.1758E+00  0.1076E+02
      4.939  0.1529E+00  0.1078E+02
      5.053  0.1224E+00  0.1079E+02
      5.167  0.8375E-01  0.1080E+02
      5.280  0.6655E-01  0.1081E+02
      5.394  0.6229E-01  0.1082E+02
      5.508  0.6033E-01  0.1082E+02
      5.622  0.8435E-01  0.1083E+02
      5.736  0.1117E+00  0.1084E+02
      5.850  0.9502E-01  0.1086E+02
      5.963  0.6426E-01  0.1086E+02
      6.077  0.6331E-01  0.1087E+02
      6.191  0.7586E-01  0.1088E+02
      6.305  0.9236E-01  0.1089E+02
      6.419  0.1343E+00  0.1090E+02
      6.532  0.2069E+00  0.1093E+02
      6.646  0.2638E+00  0.1096E+02
      6.760  0.2529E+00  0.1099E+02
      6.874  0.2088E+00  0.1101E+02
      6.988  0.1777E+00  0.1103E+02
      7.102  0.1476E+00  0.1105E+02
      7.215  0.1225E+00  0.1106E+02
      7.329  0.1247E+00  0.1108E+02
      7.443  0.1341E+00  0.1109E+02
      7.557  0.1223E+00  0.1110E+02
      7.671  0.1050E+00  0.1112E+02
      7.785  0.1095E+00  0.1113E+02
      7.898  0.1253E+00  0.1114E+02
      8.012  0.1182E+00  0.1116E+02
      8.126  0.1272E+00  0.1117E+02
      8.240  0.2564E+00  0.1120E+02
      8.354  0.5483E+00  0.1126E+02
      8.468  0.8008E+00  0.1135E+02
      8.581  0.7502E+00  0.1144E+02
      8.695  0.5236E+00  0.1150E+02
      8.809  0.3628E+00  0.1154E+02
      8.923  0.2695E+00  0.1157E+02
      9.037  0.2149E+00  0.1160E+02
      9.150  0.2384E+00  0.1162E+02
      9.264  0.3569E+00  0.1166E+02
      9.378  0.5275E+00  0.1172E+02
      9.492  0.6388E+00  0.1180E+02
      9.606  0.6216E+00  0.1187E+02
      9.720  0.5593E+00  0.1193E+02
      9.833  0.5335E+00  0.1199E+02
      9.947  0.5626E+00  0.1206E+02
     10.061  0.7455E+00  0.1214E+02
     10.175  0.1070E+01  0.1226E+02
     10.289  0.1241E+01  0.1240E+02
     10.403  0.1091E+01  0.1253E+02
     10.516  0.8271E+00  0.1262E+02
     10.630  0.6591E+00  0.1270E+02
     10.744  0.5964E+00  0.1276E+02
     10.858  0.5847E+00  0.1283E+02
     10.972  0.6183E+00  0.1290E+02
     11.085  0.6933E+00  0.1298E+02
     11.199  0.7032E+00  0.1306E+02
     11.313  0.5720E+00  0.1313E+02
     11.427  0.4254E+00  0.1317E+02
     11.541  0.3726E+00  0.1322E+02
     11.655  0.3885E+00  0.1326E+02
     11.768  0.4582E+00  0.1331E+02
     11.882  0.5053E+00  0.1337E+02
     11.996  0.4514E+00  0.1342E+02
     12.110  0.3577E+00  0.1346E+02
     12.224  0.2981E+00  0.1350E+02
     12.338  0.2908E+00  0.1353E+02
     12.451  0.3153E+00  0.1357E+02
     12.565  0.3431E+00  0.1360E+02
     12.679  0.4287E+00  0.1365E+02
     12.793  0.6421E+00  0.1373E+02
     12.907  0.8238E+00  0.1382E+02
     13.021  0.7304E+00  0.1390E+02
     13.134  0.4457E+00  0.1395E+02
     13.248  0.2247E+00  0.1398E+02
     13.362  0.1478E+00  0.1400E+02
     13.476  0.1350E+00  0.1401E+02
     13.590  0.1280E+00  0.1403E+02
     13.703  0.1321E+00  0.1404E+02
     13.817  0.1266E+00  0.1406E+02
     13.931  0.1138E+00  0.1407E+02
     14.045  0.1222E+00  0.1408E+02
     14.159  0.1390E+00  0.1410E+02
     14.273  0.1480E+00  0.1412E+02
     14.386  0.1589E+00  0.1413E+02
     14.500  0.1743E+00  0.1415E+02
     14.614  0.1844E+00  0.1417E+02
     14.728  0.2339E+00  0.1420E+02
     14.842  0.4313E+00  0.1425E+02
     14.956  0.7605E+00  0.1434E+02
     15.069  0.9764E+00  0.1445E+02
     15.183  0.9474E+00  0.1456E+02
     15.297  0.8048E+00  0.1465E+02
     15.411  0.6862E+00  0.1472E+02
     15.525  0.6656E+00  0.1480E+02
     15.638  0.7821E+00  0.1489E+02
     15.752  0.9358E+00  0.1500E+02
     15.866  0.9057E+00  0.1510E+02
     15.980  0.6724E+00  0.1518E+02
     16.094  0.4464E+00  0.1523E+02
     16.208  0.3358E+00  0.1526E+02
     16.321  0.3005E+00  0.1530E+02
     16.435  0.2716E+00  0.1533E+02
     16.549  0.2476E+00  0.1536E+02
     16.663  0.2590E+00  0.1539E+02
     16.777  0.2694E+00  0.1542E+02
     16.891  0.2319E+00  0.1544E+02
     17.004  0.1814E+00  0.1547E+02
     17.118  0.1967E+00  0.1549E+02
     17.232  0.2800E+00  0.1552E+02
     17.346  0.3541E+00  0.1556E+02
     17.460  0.3359E+00  0.1560E+02
     17.574  0.2522E+00  0.1563E+02
     17.687  0.1938E+00  0.1565E+02
     17.801  0.1866E+00  0.1567E+02
     17.915  0.2217E+00  0.1570E+02
     18.029  0.2669E+00  0.1573E+02
     18.143  0.3020E+00  0.1576E+02
     18.256  0.3211E+00  0.1580E+02
     18.370  0.2975E+00  0.1583E+02
     18.484  0.2543E+00  0.1586E+02
     18.598  0.2489E+00  0.1589E+02
     18.712  0.2407E+00  0.1592E+02
     18.826  0.1749E+00  0.1594E+02
     18.939  0.1181E+00  0.1595E+02
     19.053  0.1132E+00  0.1596E+02
     19.167  0.1256E+00  0.1598E+02
     19.281  0.1539E+00  0.1599E+02
     19.395  0.1970E+00  0.1602E+02
     19.509  0.2004E+00  0.1604E+02
     19.622  0.1576E+00  0.1606E+02
     19.736  0.1235E+00  0.1607E+02
     19.850  0.1206E+00  0.1608E+02
     19.964  0.1527E+00  0.1610E+02
     20.078  0.2510E+00  0.1613E+02
     20.191  0.4276E+00  0.1618E+02
     20.305  0.6204E+00  0.1625E+02
     20.419  0.6812E+00  0.1633E+02
     20.533  0.5239E+00  0.1639E+02
     20.647  0.3054E+00  0.1642E+02
     20.761  0.2215E+00  0.1645E+02
     20.874  0.3155E+00  0.1648E+02
     20.988  0.5010E+00  0.1654E+02
     21.102  0.6177E+00  0.1661E+02
     21.216  0.5759E+00  0.1668E+02
     21.330  0.4409E+00  0.1673E+02
     21.444  0.3298E+00  0.1676E+02
     21.557  0.2996E+00  0.1680E+02
     21.671  0.3279E+00  0.1683E+02
     21.785  0.3637E+00  0.1688E+02
     21.899  0.3927E+00  0.1692E+02
     22.013  0.4378E+00  0.1697E+02
     22.127  0.4503E+00  0.1702E+02
     22.240  0.3772E+00  0.1706E+02
     22.354  0.3065E+00  0.1710E+02
     22.468  0.2889E+00  0.1713E+02
     22.582  0.2657E+00  0.1716E+02
     22.696  0.2208E+00  0.1719E+02
     22.809  0.1764E+00  0.1721E+02
     22.923  0.1642E+00  0.1723E+02
     23.037  0.2244E+00  0.1725E+02
     23.151  0.3184E+00  0.1729E+02
     23.265  0.3531E+00  0.1733E+02
     23.379  0.3150E+00  0.1736E+02
     23.492  0.2483E+00  0.1739E+02
     23.606  0.1891E+00  0.1741E+02
     23.720  0.1680E+00  0.1743E+02
     23.834  0.1750E+00  0.1745E+02
     23.948  0.1739E+00  0.1747E+02
     24.062  0.1752E+00  0.1749E+02
     24.175  0.1990E+00  0.1752E+02
     24.289  0.2053E+00  0.1754E+02
     24.403  0.1887E+00  0.1756E+02
     24.517  0.2143E+00  0.1758E+02
     24.631  0.3115E+00  0.1762E+02
     24.745  0.4416E+00  0.1767E+02
     24.858  0.5263E+00  0.1773E+02
     24.972  0.4817E+00  0.1779E+02
     25.086  0.3380E+00  0.1782E+02
     25.200  0.2201E+00  0.1785E+02
     25.314  0.1814E+00  0.1787E+02
     25.427  0.1823E+00  0.1789E+02
     25.541  0.1901E+00  0.1791E+02
     25.655  0.2135E+00  0.1794E+02
     25.769  0.2300E+00  0.1796E+02
     25.883  0.1888E+00  0.1798E+02
     25.997  0.1015E+00  0.1800E+02
     26.110  0.3287E-01  0.1800E+02
     26.224  0.6094E-02  0.1800E+02
     26.338  0.6288E-03  0.1800E+02
     26.452  0.3548E-04  0.1800E+02
     26.566  0.1082E-05  0.1800E+02
     26.680  0.1772E-07  0.1800E+02
     26.793  0.1548E-09  0.1800E+02
     26.907  0.7195E-12  0.1800E+02
     27.021  0.1775E-14  0.1800E+02
     27.135  0.0000E+00  0.1800E+02
     27.249  0.0000E+00  0.1800E+02
     27.362  0.0000E+00  0.1800E+02
     27.476  0.0000E+00  0.1800E+02
     27.590  0.0000E+00  0.1800E+02
     27.704  0.0000E+00  0.1800E+02
     27.818  0.0000E+00  0.1800E+02
     27.81778758     -6.32977860  301      6.82147560      1.00000000
     -6.330  0.0000E+00  0.0000E+00  0.0000E+00
     -6.216  0.0000E+00  0.0000E+00  0.0000E+00
     -6.102  0.0000E+00  0.0000E+00  0.0000E+00
     -5.988  0.0000E+00  0.0000E+00  0.0000E+00
     -5.874  0.0000E+00  0.0000E+00  0.0000E+00
     -5.761  0.0000E+00  0.0000E+00  0.0000E+00
     -5.647  0.0000E+00  0.0000E+00  0.0000E+00
     -5.533  0.0000E+00  0.0000E+00  0.0000E+00
     -5.419  0.1929E-14  0.5964E-18  0.4001E-20
     -5.305  0.7863E-12  0.2989E-15  0.2351E-17
     -5.192  0.1705E-09  0.7583E-13  0.6357E-15
     -5.078  0.1971E-07  0.1036E-10  0.9460E-13
     -4.964  0.1219E-05  0.7693E-09  0.7840E-11
     -4.850  0.4056E-04  0.3142E-07  0.3664E-09
     -4.736  0.7305E-03  0.7167E-06  0.9796E-08
     -4.622  0.7186E-02  0.9302E-05  0.1527E-06
     -4.509  0.3911E-01  0.7038E-04  0.1430E-05
     -4.395  0.1202E+00  0.3214E-03  0.8395E-05
     -4.281  0.2168E+00  0.9343E-03  0.3273E-04
     -4.167  0.2476E+00  0.1869E-02  0.9102E-04
     -4.053  0.2079E+00  0.2841E-02  0.1956E-03
     -3.939  0.1593E+00  0.3660E-02  0.3530E-03
     -3.826  0.1299E+00  0.4338E-02  0.5659E-03
     -3.712  0.1127E+00  0.4864E-02  0.8352E-03
     -3.598  0.1021E+00  0.5428E-02  0.1199E-02
     -3.484  0.9566E-01  0.5952E-02  0.1629E-02
     -3.370  0.8695E-01  0.6189E-02  0.2141E-02
     -3.256  0.8399E-01  0.6874E-02  0.2925E-02
     -3.143  0.8251E-01  0.7395E-02  0.3736E-02
     -3.029  0.7479E-01  0.7324E-02  0.4740E-02
     -2.915  0.7435E-01  0.8085E-02  0.6383E-02
     -2.801  0.7655E-01  0.8906E-02  0.8291E-02
     -2.687  0.7172E-01  0.8752E-02  0.1065E-01
     -2.574  0.6988E-01  0.9117E-02  0.1448E-01
     -2.460  0.7611E-01  0.1047E-01  0.2125E-01
     -2.346  0.8055E-01  0.1090E-01  0.3642E-01
     -2.232  0.9310E-01  0.1115E-01  0.8714E-01
     -2.118  0.1342E+00  0.1268E-01  0.2184E+00
     -2.004  0.1790E+00  0.1311E-01  0.4118E+00
     -1.891  0.1829E+00  0.1202E-01  0.5792E+00
     -1.777  0.1555E+00  0.1265E-01  0.7192E+00
     -1.663  0.1319E+00  0.1579E-01  0.8826E+00
     -1.549  0.1141E+00  0.1904E-01  0.1055E+01
     -1.435  0.8905E-01  0.1899E-01  0.1145E+01
     -1.321  0.6156E-01  0.1625E-01  0.1119E+01
     -1.208  0.4623E-01  0.1530E-01  0.1094E+01
19楼2009-11-26 16:13:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chimegreen

木虫 (正式写手)


强人啊
此贴要顶
谢谢lz
20楼2009-11-26 16:20:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaiwddf

金虫 (小有名气)


★★★★★ 五星级,优秀推荐

如此伟大的人真不多见,真诚的赞一个
21楼2009-11-26 18:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqpie

木虫 (职业作家)


★★★★★ 五星级,优秀推荐

楼主强人!@@
22楼2009-11-26 22:27:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fox-fox

木虫 (知名作家)


★★★★★ 五星级,优秀推荐

学习了。
24楼2009-11-27 12:47:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddang100

银虫 (小有名气)


强烈支持
26楼2009-11-28 17:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
修正了很多bug,请及时更新.谢谢!
27楼2009-12-04 17:59:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yunml

铜虫 (小有名气)


引用回帖:
Originally posted by wangle_xq at 2009-11-26 00:43:
支持,回头发送测试结果。
不知道有没有考虑组织一个开发团队。我有兴趣参与。打造中国人自己的VASP处理软件。


希望哪天有中国人自己写的VASP,当然还要超过它。处理软件,可以交给三弟做。
28楼2009-12-05 09:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
下一步目标是写简单的第一性原理程序!希望有兴趣的朋友多交流,再次谢谢大家的热心讨论
29楼2009-12-05 10:22:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天路堂主

银虫 (初入文坛)


引用回帖:
Originally posted by 灯塔守望者 at 2009-11-25 13:32:29:
相信很多人都为vasp计算结果的处理感到很头痛,网上虽然有很多处理vasp计算结果的软件或者脚本,不过要么只针对于某一个问题,要么使用不太友好。为此,我弄了一个软件包 (vaspkit),把结构可视化,提取态密度,能 ...

为什么我得出的pdos.dat文件只有一个?后面的把前面的数据都覆盖掉了,怎么样得到每个原子的pdos.dat文件,比如说像1pdos.dat,2pdos.dat,3pdos.dat等等?
30楼2010-12-03 15:30:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awmc2008

至尊木虫 (文坛精英)


★★★★★ 五星级,优秀推荐

学习学习,谢谢楼主。
31楼2010-12-03 21:44:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowenzhi16

铁杆木虫 (正式写手)


★★★★★ 五星级,优秀推荐

感谢灯塔给我们的光明。这东西是好东西。不知道后期的功能实现了没有。
32楼2010-12-04 08:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)


★★★★★ 五星级,优秀推荐

不错,下来玩玩
33楼2010-12-04 10:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kgdu

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

谢谢了,支持
34楼2010-12-04 13:22:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangtr

金虫 (小有名气)


★★★★★ 五星级,优秀推荐

zan ~~
35楼2010-12-05 17:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddang100

银虫 (小有名气)


★★★★★ 五星级,优秀推荐

不错,大家努力
36楼2010-12-05 19:00:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youzhizhe(金币+1):期待你的更新哦,抓紧哈…… 2010-12-05 22:50:38
谢谢各位同行的支持,过段时间,我打算更新一下这个程序。以后多交流
37楼2010-12-05 19:04:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkleof

银虫 (正式写手)


引用回帖:
Originally posted by 灯塔守望者 at 2009-11-25 13:32:29:
相信很多人都为vasp计算结果的处理感到很头痛,网上虽然有很多处理vasp计算结果的软件或者脚本,不过要么只针对于某一个问题,要么使用不太友好。为此,我弄了一个软件包 (vaspkit),把结构可视化,提取态密度,能 ...

你是vaspkit的原作者?我用过vaspkit-0.14,发现一些问题,修改过其中一些模块的代码。其中能带处理的部分在体系比较大,比如nbd=600,nkd=300的情况下会出错,这个问题的根源好像是数组定义超出了fortran里面最大容许值,但是一些超胞的计算里nbd确实很大。还有从PROCAR里读取数据的部分好像原来的代码不对,或者可能和vasp5的格式不匹配,结果读取出来的轨道占据因子occ不正确,还有不知道为什么设定读取的列数为4,是设计时仅考虑包括p轨道的情况?在有d轨道的情况下读取的数据列应该为9啊。另外在处理CHGCAR文件时会由于文件中的元素符号一行出现错误,这个应该是由于vasp5的格式改变引起的了,不知道在新版里面有没有更新。
vaspkit是一个不错的程序包,lz能够开发并公布出来给大家用是非常热心的给大家做了贡献,希望能够不断的开发和完善。
38楼2010-12-06 16:56:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by nkleof at 2010-12-06 17:56:09:

你是vaspkit的原作者?我用过vaspkit-0.14,发现一些问题,修改过其中一些模块的代码。其中能带处理的部分在体系比较大,比如nbd=600,nkd=300的情况下会出错,这个问题的根源好像是数组定义超出了fortran里面最 ...

谢谢和欢迎你的提议和修改。原来编写这个程序是基于vasp 4.6, vasp 5.2的格式与4.6在CHGCAR文件格式等方面略有不同。当时花了一周时间在边学fortran边写vaspkit的情况下弄得,程序的确有一些粗糙。以至于升级变得很麻烦,我想把一部分重写一下,采用module传递变量,以后增加电荷密度的xrysden格式,以及差分电荷密度,自旋向上向下等电荷密度等等(先简单后复杂)。最近实在有些忙,欢迎其他朋友提出宝贵意见,谢谢!

[ Last edited by 灯塔守望者 on 2010-12-6 at 20:09 ]
39楼2010-12-06 19:07:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

za5010

铁虫 (初入文坛)


★★★ 三星级,支持鼓励

支持楼主,鼓励原创!
拿来借鉴,有需要还是自己写来得顺手...
40楼2010-12-06 21:30:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowenzhi16

铁杆木虫 (正式写手)


运用VASPKIT确实不错,但是在能带输出时候似乎总等到不能输出上下自旋能带的,上自旋没有数据输出,文件为空,下自旋有数据但是似乎不对。提示为:Floating point exception.不知道什么原因。多谢灯塔守望的不懈努力。
41楼2010-12-29 15:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddang100

银虫 (小有名气)


用的时候每次都要吧所有的文件拷到vaspkit0.14文件夹中吗?能不能直接在工作目录下运行呢,否则挺麻烦的

[ Last edited by ddang100 on 2011-3-4 at 17:31 ]
42楼2011-03-04 17:05:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by ddang100 at 2011-03-04 18:05:47:
用的时候每次都要吧所有的文件拷到vaspkit0.14文件夹中吗?能不能直接在工作目录下运行呢,否则挺麻烦的

[ Last edited by ddang100 on 2011-3-4 at 17:31 ]

没有必要啊,你把vaspkit放到~/bin/
每次直接在你的作业下运行~/bin/vaspkit,不用拷贝你的数据。
43楼2011-03-05 13:28:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by xiaowenzhi16 at 2010-12-29 16:48:45:
运用VASPKIT确实不错,但是在能带输出时候似乎总等到不能输出上下自旋能带的,上自旋没有数据输出,文件为空,下自旋有数据但是似乎不对。提示为:Floating point exception.不知道什么原因。多谢灯塔守望的不懈努 ...

这个程序的确有一些bug,现在事情太多了,实在想重新写,但总是力不从心。希望有朋友参与就好了。
44楼2011-03-05 13:29:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天路堂主

银虫 (初入文坛)


★★★ 三星级,支持鼓励

引用回帖:
Originally posted by 灯塔守望者 at 2009-11-25 13:32:29:
相信很多人都为vasp计算结果的处理感到很头痛,网上虽然有很多处理vasp计算结果的软件或者脚本,不过要么只针对于某一个问题,要么使用不太友好。为此,我弄了一个软件包 (vaspkit),把结构可视化,提取态密度,能 ...

为什么我用vaspkit得到的各原子的pdos加起来不等于总的DOS呢?求指教!
46楼2011-04-25 15:15:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nfsex

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

你好,楼主
     我首先输入vaspkit
             然后我输入ls
                    出现:At line 51 of file mainmenu.f90
Fortran runtime error: Bad integer for item 1 in list input
  这是为什么呢,谢谢楼主,我用的是vasp5.2
47楼2011-04-26 11:45:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

llg1987

金虫 (知名作家)


★★★★★ 五星级,优秀推荐

感谢楼主啊,好资源啊
48楼2011-04-27 09:00:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

llg1987

金虫 (知名作家)


0.8

谢谢楼主提供的好资源
49楼2011-04-27 09:01:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
wuli823楼
2009-11-26 22:51   回复  
 
zdhlover25楼
2009-11-27 22:27   回复  
 
sg1840892645楼
2011-03-05 18:30   回复  
五星好评  
hnuzhoulin50楼
2011-05-04 09:59   回复  
五星好评  顶一下,感谢分享!
相关版块跳转 我要订阅楼主 灯塔守望者 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见