24小时热门版块排行榜    

北京石油化工学院2025年硕士研究生招生接受调剂公告
查看: 17954  |  回复: 165
【奖励】 本帖被评价127次,作者souledge增加金币 100
本帖产生 1 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 强势增强CASTEP输出的工具:check2xsf (Windows版)

相信用过Materials Studio下CASTEP的人都知道,CASTEP的图形输出,也就是所谓的post-processing是基于MS提供的功能的,除了例如能带之类的点-线图可以轻易的导出到Origin或者MatLAB之类做出比较专业的图象,其他的例如电荷密度之类的,导出重画一是麻烦,二是不能得到别人论文中那样漂亮的图片。几个免费做结构图的软件,比如Molekel,XCrySDen等,做图的质量基本是秒杀MS的模型图的。
更更郁闷的,是如果在Linux下直接用RunCASTEP.sh运行计算,计算后的结果很难再导回到Materials Studio的图形界面,给结果分析造成了困难。
跟别的基于DFT的平面波赝势软件相比,CASTEP的输出种类也的确偏少,总是缺乏别人那么花样繁多的结果。比如该帖中的问题:http://muchong.com/bbs/viewthread.php?tid=3701889&fpage=3,取出HOMO或者LUMO的电荷密度,就无法直接使用MS来解决。
这也怪不得CASTEP,毕竟CASTEP是商业软件,其核心结果文档(check文件)的格式并没有公开,所以很难存在什么软件轻易读取它的结果。
事实上,这些问题早就得到了很大程度上的解决了。因为,CASTEP至少在英国学术界内,是开源的。我们这里要感谢英国剑桥大学卡文迪实实验室凝聚态物理组的Michael Rutter,他为我们提供了一个工具check2xsf:http://www.tcm.phy.cam.ac.uk/~mjr/vis/。
该工具是开源的,在Unix/Linux下非常容易编译,只需要修改Makefile中一两个参数就可以直接make了。但是,我们用MS更多是在Windows下。
这两天使用Visual Studio,总算成功在Windows平台下成功编译了check2xsf,最大的困难已经不是修改编译时的参数,而是要修改几处C代码,因为原始的代码使用的是POSIX UNIX的标准C库函数,放在Windows下显然是不能使用的~编译好的程序在附件中。
下面来做个例子,说明怎么使用这个软件。
使用ZnO作为例子,ZnO作为直接带隙半导体材料,有着明确的HOMO和LUMO,并且位于Gamma点上。我们本次任务就是提取HOMO和LUMO在Gamma点上的电荷密度。首先,对ZnO进行计算,计算时一定要选择回传完整check文件:

然后,计算完成,到存储计算文件的文件夹中,显示所有隐藏文件,把扩展名为check的文件取消隐藏属性。
这时,需要再做两个检查:1,HOMO和LUMO分别是第几条能带;2,第几个k点是Gamma点。
对于问题1,我们可以看能带图。

对于问题2,我们既可以看扩展名为castep的文件,又可以使用命令
CODE:
check2xsf -l --cell ZnO.check ZnO_k.cell

,然后查看ZnO_k.cell。不管哪个方法,都可以看到如下的结果:

这说明了,Gamma点是第48号k点。
然后,就可以提取数据了。数据保存为cube格式,因为这个格式几乎可以用任何晶体作图软件打开,不仅仅限于check2xsf的作者推荐的XCrySDen和Jmol。我在这里使用的是V_Sim,这个软件有Windows版,可以到官方网站随便下载。不会使用V_Sim又想学习的,可以到该帖看我以前制作的教程:http://muchong.com/bbs/viewthread.php?tid=4923069&fpage=1
首先,是HOMO的结果。显示3D的等值面图,可以看到,如同Abinit区的那个帖子,我们不但得到了HOMO的等值面,还可以对不同位相使用不同颜色。所以,特此更正一下,这里输出的其实不是Gamma点处HOMO能级的电荷密度,而是该点的波函。我们知道,波函是有位相的,但是电荷密度是,是不可能有位相并且只有非负值的。更有特色的一点在于,输出波函,在Abinit和PWscf中都不支持,或者只能用wannier方程计算才能得到有位相区分的图(VASP不清楚,暂时保留意见)。

然后,是对某个截面的波函分布的平面图:

再再然后,是Gamma点处LUMO的波函:

详细的check2xsf的使用命令如下:
引用回帖:
Usage: check2xsf [-aAbBckmRsSvx] [--FORMAT] infile [outfile]

-a           rotate as though outputing in abc format
-A           accumulate (sum) requested bands
-b[=range]   include bands (as psi)
-B[=range]   include bands (as densities)
-c           include charge density
-d           read also a .chdiff file, constructing its name from
               the .cell or .check file given
-H           shift atoms by half a grid cell
-k=range     include given k-points (default 1) for bands
-l           list k-points in .cell output, not MP parameters
-L           produce (incorrect) left-handed abc output
-m[=a,b,c]   assume input is molecule, not crystal, and move by
               given nos of grid cells, or move automatically
-R           don't rescale densities
-s           include spin densities
-S=range     include given spins (0 or 1) for bands
-t=(x1,y1,z1)(x2,y2,z2)[(x3,y3,z3)]
             rotate coords so 1st vector becomes 2nd, using third
             vector as axis, else axis is perpendicular to others
-T=(x1,y1,z1)(x2,y2,z2)[(x3,y3,z3)]
             ditto, but first vec expressed in absolute coords
-v           be verbose (may be repeated)
-x           expand rhombohedral cell to cubic/hexagonal automatically
-x=(x1,y1,z1)(x2,y2,z2)(x3,y3,z3)
             re-express in new basis given in terms of old
-X=(x1,y1,z1)(x2,y2,z2)(x3,y3,z3)
             re-express in new basis given in absolute terms
-3           if swapping axes to convert lhs to rhs, third axis
               is special & first two swapped. Else first is special.

FORMAT is one of: xsf       XCrySDen (default)
                  cell      CASTEP .cell, cartesian and fractional
                  cell_abc                abc and fractional
                  cell_abs                cartesian and absolute
                  cell_abc_abs            abc and absolute
                  cml       Chemical Markup Language
                  cube      Gaussian cube
                  dx        OpenDX
                  fdf       Flexible Data Format (Siesta) (beta)
                  null      Discard output
                  pdb       PDB
                  pdbn      PDB with atoms numbered
                  vasp      VASP output
                  xplor     Xplor
                  xyz       XYZ

range specifies band numbers as "a,b-c,d"
-b and -B are mutually exclusive. Only one of the x and t options may be given

Input files ending .pdb are assumed to be in pdb format, otherwise
automatic detection of .cell or .check input. Compatible with
.check files from CASTEP 3.0 to 5.5 (and perhaps beyond).

Version 1.04, (c) MJ Rutter 2007, 2008, 2011 licenced under the GPL v2

回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : check2xsf.rar
  • 2013-01-10 14:41:48, 115.45 K

» 收录本帖的淘帖专辑推荐

纳米技术与能源及模拟 精华网帖收集 万卷阁之搜神记 VASP
第一原理资料汇编 第一性原理计算-tg Material Studio 科研技巧
MS第一性原理 资源 计算模拟-MS-CASTEP 第一性原理计算辅助工具
物理前沿 MS资料 生活休闲 软件及使用
VASP 计算资料 @个人收集@ 科研必备知识
量化软件 差分电荷密度 Photochemistry 杂七杂八
MS学习笔记 good resources 计算化学 第一性原理
MS学习 我的分享 田凯 藏宝阁
晶体软件 计算化学相关 电子书 模拟计算
第一原理

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paopaotu326

金虫 (正式写手)


送红花一朵
楼主好,我的系统也是win7 64位的,打开cmd后,进入命令文件的文件夹,输入命令,结果提示: no Input file specified.
如图片所示:

https://pan.baidu.com/s/1_IA6Yqs-YG3X3-3imwxAwg

咋回事?
152楼2018-04-10 11:40:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paopaotu326

金虫 (正式写手)


还有,v-sim好像只能在win32下运行,有没有其他推荐呢?可以在win64下运行的
153楼2018-04-10 11:41:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
paopaotu326151楼
2018-04-10 11:19   回复  
五星好评  顶一下,感谢分享!
相关版块跳转 我要订阅楼主 souledge 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 公办高校,不坐班,工作日没课时去外地也要请假吗? +4 河西夜郎 2025-04-12 4/200 2025-04-15 22:59 by 豫椒
[考研] 工科274求调剂 +4 madankuaipao 2025-04-15 5/250 2025-04-15 22:46 by kaoyan的家
[考研] 青岛科技大学环境学院环境科学专业二次调剂名额 +5 lixin1994 2025-04-09 5/250 2025-04-15 10:21 by 是西瓜熟了
[考研] 0856,总分250,英一数二,一志愿哈工大,四级已过,校二等奖学金,校优秀团员 +3 仲s2 2025-04-13 3/150 2025-04-15 09:17 by taosheng130
[考研] 复试求调剂 +5 yjlお 2025-04-11 6/300 2025-04-14 23:11 by 白白12138
[考研] 工学271求老师收留 +4 qjj0914 2025-04-12 8/400 2025-04-14 10:39 by edgexjy
[论文投稿] 论文二审这种审稿意见返修后还需要再发审稿人审一次吗? 5+4 qweasd12345 2025-04-12 5/250 2025-04-14 09:33 by 北京莱茵润色
[考研] 求调剂 +5 benben0216 2025-04-13 16/800 2025-04-13 22:03 by zhouyuwinner
[考研] 专硕数二英二材料专业085601分数308本科西南交通大学 +4 上岸成功丁 2025-04-13 4/200 2025-04-13 21:36 by zglnasyy
[考研] 已录取,此帖作废,感谢各位关注,祝好! +5 carpr 2025-04-13 6/300 2025-04-13 19:50 by carpr
[考研] 241求调剂 +3 叶叶287 2025-04-11 5/250 2025-04-13 10:20 by 小小毅prince
[考研] 250求调剂 +6 仲s2 2025-04-12 10/500 2025-04-13 09:08 by 仲s2
[考研] 调剂材料306 +28 小曌 2025-04-10 31/1550 2025-04-12 05:28 by 星海S~
[考研] 复试调剂 +24 邹gv 2025-04-10 31/1550 2025-04-11 21:32 by zyh001021
[考研] 08工学求调剂 +4 范俊鹏 2025-04-11 4/200 2025-04-11 12:05 by 褚迎松小小
[考研] 260求调剂 +3 cccmmmy 2025-04-11 8/400 2025-04-11 10:00 by wjmyutian
[考研] 331求调剂 +5 灰咸鱼 2025-04-10 8/400 2025-04-11 09:24 by 徐小七七
[考研] 工科085500,总分275 +5 jsha01 2025-04-10 5/250 2025-04-11 09:17 by 云边小卖部。
[考研] 求085400的调剂 +4 enter!!! 2025-04-09 7/350 2025-04-10 23:56 by enter!!!
[考研] 282求调剂 +3 谦谦念念Dpp 2025-04-09 3/150 2025-04-09 17:27 by ??⒍
信息提示
请填处理意见
婢跺秷鐦拫鍐ㄥ婢跺秷鐦拫鍐ㄥ
婢跺秷鐦拫鍐ㄥ娴滃瞼娣惍锟�閼惧嘲褰囨径宥堢槸 & 鐠嬪啫澧忛幐鍥у础閿涳拷瀵邦喕淇婇幍顐$閹碉拷
閸熷棗濮熼崥鍫滅稊閸熷棗濮熼崥鍫滅稊
閸熷棗濮熼崥鍫滅稊娴滃瞼娣惍锟�閸熷棗濮熼崥鍫滅稊閸溿劏顕�瀵邦喕淇婇幍顐$閹碉拷(婢跺洦鏁為敍姘櫌閸斺€虫値娴o拷)
鐎涳附婀抽獮鑼舵彛鐎涳附婀抽獮鑼舵彛
鐎涳附婀抽獮鑼舵彛娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷娴滃棜袙閺堚偓閺傛澘顒熼張顖氬叡鐠愶拷
鐠佺儤鏋冩潏鍛嚤鐠佺儤鏋冩潏鍛嚤
鐠佺儤鏋冩潏鍛嚤娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷濞h濮為懓浣哥瑎閿涘苯鍘ょ拹鐟版尒鐠囷拷
鐢牐绁幏娑樺珯鐢牐绁幏娑樺珯