24小时热门版块排行榜

返回列表

【奖励】本帖被评价12次，作者xujc1983增加金币 9.5 个

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xujc1983

木虫 (著名写手)

1ST强帖: 1
应助: 71 (初中生)
金币: 2065.3
帖子: 1096
在线: 575小时
虫号: 285496

[资源] 【原创】通过CASTEP产生POSCAR文件【已搜索“CASTEP”+“POSCAR”，无重复】

前几天回帖说过自己产生POSCAR文件是通过CASTEP完成的，具体过程如下：
1.  在CASTEP中把模型建好。
2.  模块中的CASTEP，计算，单击“Run”旁边的“Files...”按钮，然后点击“Save Files”按钮，弹出对话框，选“No”，将相关文件保存到project中，会出现类似计算出现的下一级文件夹。
3.  从电脑中找到此文件夹，找到***.cell文件。此文件为隐藏文件，***为MS中模型的名字。有用的为前几行，以BCC-Fe单胞为例，产生的Fe.cell文件开头内容如下：
%BLOCK LATTICE_CART
   2.866400000000000    0.000000000000000    0.000000000000000
   0.000000000000000    2.866400000000000    0.000000000000000
   0.000000000000000    0.000000000000000    2.866400000000000
%ENDBLOCK LATTICE_CART

%BLOCK POSITIONS_FRAC
Fe 0.0000000000000000 0.0000000000000000 0.0000000000000000
Fe 0.5000000000000000 0.5000000000000001 0.5000000000000001
%ENDBLOCK POSITIONS_FRAC
......
4.  书写POSCAR文件。
第一行，自定义体系名，随便自己写了，以自己方便为好
第2-5行为晶格缩放因子（晶格常数）和晶格矢量，可以有两种写法
第一种：
1.0000
   2.866400000000000    0.000000000000000    0.000000000000000
   0.000000000000000    2.866400000000000    0.000000000000000
   0.000000000000000    0.000000000000000    2.866400000000000
第二种：
2.866400001.0000
   1.000000000000000    0.000000000000000    0.000000000000000
   0.000000000000000    1.000000000000000    0.000000000000000
   0.000000000000000    0.000000000000000    1.000000000000000
推荐第二种写法
第6行，原子数目。应该按照种类分别写，后面坐标以及POTCAR中的顺序应与此处相同。
第7行，给出的坐标类型。“cartesian”和“direct”之分，分别是笛卡尔坐标和分数坐标，通常用第2中，分数坐标。
第8行及其以后，按照第6行原子种类排序的各个原子的分数坐标。
注：如果要考虑，部分坐标驰豫的话，需要在第7行前加入一行，“Selective dynamics”，并在各个原子坐标后，以“T”，“F”标记那个坐标需要驰豫
BCC－Fe产生POSCAR如下，具体书写格式可参考VASP Guide
BCC-Fe
2.866400000000000
   1.0000000000000000 0.0000000000000000 0.0000000000000000
   0.0000000000000000 1.0000000000000000 0.0000000000000000
   0.0000000000000000 0.0000000000000000 1.0000000000000000
2
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.5000000000000000  0.5000000000000000  0.5000000000000000

[ Last edited by xujc1983 on 2009-6-9 at 19:40 ]

回复此楼