| 查看: 1103 | 回复: 2 | |||
smartcheng铁虫 (小有名气)
|
[求助]
castep中限制原子只沿某个方向运动
|
|
我看前人的帖子说，只有修改.cell 文件，才能做到按方向性限制。如下所示，如果要限制第二个Al原子不能沿着X方向移动 4 Al 2 1.0000000000 0.0000000000 0.0000000000 5 Al 2 0.0000000000 0.0000000000 1.0000000000 这样怎么会达到限制x方向的目的，这里面的1和0都代表什么啊。 原帖： Q3：如何在CASTEP计算中限制某个原子的移动方向？ A3：CASTEP计算中只支持Fraction坐标，不支持Coordinate坐标。采用常规方式无法限制原子的移动方向，需要手动修改CASTEP的Cell文件，具体方式如下，在 %BLOCK POSITIONS_FRAC Al 0.0000000000 0.0000000000 0.0000000000 Al 0.2270171283 0.2094837802 0.5009850469 Al 0.5000000000 0.0000000000 0.5000000000 Al 0.5000000000 0.5000000000 0.0000000000 %ENDBLOCK POSITIONS_FRAC 中包括每个原子的分数坐标，请留意每个原子的编号，在下面的计算中会有所帮助。 在Cell文件中，关于原子坐标限制的部分如下： FIX_COM : false %BLOCK IONIC_CONSTRAINTS 1 Al 1 1.0000000000 0.0000000000 0.0000000000 2 Al 1 0.0000000000 1.0000000000 0.0000000000 3 Al 1 0.0000000000 0.0000000000 1.0000000000 4 Al 3 1.0000000000 0.0000000000 0.0000000000 5 Al 3 0.0000000000 1.0000000000 0.0000000000 6 Al 3 0.0000000000 0.0000000000 1.0000000000 7 Al 4 1.0000000000 0.0000000000 0.0000000000 8 Al 4 0.0000000000 1.0000000000 0.0000000000 9 Al 4 0.0000000000 0.0000000000 1.0000000000 %ENDBLOCK IONIC_CONSTRAINTS 其中包括每行编号（紫色），原子类型（红色），原子编号（绿色）和需要限制的分数坐标（黑色）。如果要限制第二个Al原子不能沿着X方向移动，那么在该部分的设置如下： 1 Al 1 1.0000000000 0.0000000000 0.0000000000 2 Al 1 0.0000000000 1.0000000000 0.0000000000 3 Al 1 0.0000000000 0.0000000000 1.0000000000 4 Al 2 1.0000000000 0.0000000000 0.0000000000 5 Al 2 0.0000000000 0.0000000000 1.0000000000 6 Al 3 1.0000000000 0.0000000000 0.0000000000 7 Al 3 0.0000000000 1.0000000000 0.0000000000 8 Al 3 0.0000000000 0.0000000000 1.0000000000 9 Al 4 1.0000000000 0.0000000000 0.0000000000 10 Al 4 0.0000000000 1.0000000000 0.0000000000 11 Al 4 0.0000000000 0.0000000000 1.0000000000 注意，编号要顺次后移，否则会出错，黑色部分分别对应于x、y、z三个方向，完成后，使用Run Files方式即可进行计算。 此外，在 %BLOCK CELL_CONSTRAINTS 1 2 3 4 5 6 %ENDBLOCK CELL_CONSTRAINTS 中可以限制晶胞向量的变化，第一行分别表示三个晶胞参数a、b、c， 第二行表示晶胞角度α、β、γ。上述结构表明参数可以被优化，如果需要加以限制，则需要改为0。 |
回复此楼» 收录本帖的淘帖专辑推荐
 计算MS等 |
» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 大家自由Ni原子是怎么计算，能使它收敛的CASTEP
已经有12人回复
- CASTEP原子布居分析的问题
已经有4人回复
- CASTEP超胞边界上的原子为什么不能单个被替换？
已经有8人回复
- 急！求助MoSe2的原子位置，要用castep做计算，但是找不到原子位置。
已经有4人回复
- 原子的分数坐标 用CASTEP手动建模用的
已经有4人回复
- 单原子能量的castep计算
已经有8人回复
- Castep 优化时，原子坐标位置的优化
已经有10人回复
- castep 计算界面能，几何优化后原子完全混乱排列什么原因？
已经有16人回复
- 求助：用castep 做掺杂计算时，掺杂原子怎么添加到原有晶格中？
已经有14人回复
- 请问96G的内存, 12核, 2.66GHz 的电脑配置, CASTEP能计算最多的原子数 ？
已经有23人回复
- CASTEP计算间隙式原子掺杂，怎么加掺杂的原子？
已经有5人回复
- 请问castep计算的时候如何来固定部分原子
已经有3人回复
- 大家用castep计算单原子O或N有没有出错？
已经有11人回复
- 【求助】在用Castep算出能带，怎么观察一种物质中各个原子的分态（S,P等）的态密度图
已经有9人回复
- 【求助】Castep中如何计算自由原子的能量？
已经有10人回复
- 【求助】castep原子电荷计算问题
已经有3人回复
- 【求助】castep建模 为什么原子之间化学键没有自动生成 要怎么添加呢
已经有8人回复
- 【求助】表面吸附原子的位置问题！
已经有9人回复
泽润东方
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1702.8
- 散金: 406
- 帖子: 470
- 在线: 222.7小时
- 虫号: 684187
- 注册: 2008-12-29
- 性别: GG
- 专业: 基础物理学
2楼2015-04-01 14:52:32
rollingveau
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 35.7
- 散金: 72
- 帖子: 128
- 在线: 72.8小时
- 虫号: 5051340
- 注册: 2016-09-23
- 性别: GG
- 专业: 理论和计算化学
3楼2019-12-16 10:13:29