卡开发发 发表于 2024-11-29 17:27 谢谢卡卡老师解答，脚本没问题，我坐标轴的问题，按照您的3的操作做了，现在已经成功导出了！ |
Flase 发表于 2024-11-29 17:22 好的我试试，谢谢！ |
|
1、脚本没啥问题啊。只是因为你切表面的时候c轴放到x，所以(26.xxx 0 0)就是c轴这是完全合理的，这应该不会影响到结果，只是计算后的朝向看上去比较怪而已。 2、之所以导出到vesta没问题的原因是因为导出后会不一样是因为存到cif的时候是按照abc长度和夹角存的，然后导入到vesta的时候vesta自身读取规则是按照a沿着x罢了。 3、要解决这个问题或许应该在切面的时候让a沿着x，在build->cleave surface的对话框有个options选项卡，然后让u along x, v in xy plane，你可以这样试试看再用脚本是否正常。或者已经切好面你可以在build->lattice parameters对话框的advance选项卡进行类似设置。 4、除了上面的方案，你也可以导出cif用ase的gui进行原子固定，安装好ase之后，你可以用ase gui Ag001.cif启动ase的图形界面并加载这个cif结构，可以通过框选的方式选中多个原子然后在菜单栏tools->constraints对话框点击fix按钮，最后存成Ag.vasp或者POSCAR_Ag即可。 |
|
vaspkit 402 1-选择POSCAR 3-按高度选择 0空格x：指x高度下所有原子固定，高度可以从vesta看 1或者2-选择坐标系 all或者1-3：固定所有坐标系，也可以自己选择 最后会输出一个POSCAR-FIX文件 |
snljty2 发表于 2024-11-29 16:19 谢谢大佬！ |
|
这种事用五分钟写个小脚本就行了，下面链接是我随手写的，啥都没多想，变量名和实现方式都有点儿丑。 把原来的POSCAR重命名为POSCAR.old，用Multiwfn载入POSCAR.old，输入gi然后按两下回车。GaussView打开生成的POSCAR.gjf，用鼠标各种方式选取好要固定的原子，点Tools -> Atom Selection，把里面的原子编号复制出来，粘贴到这个脚本开头的fixed_indices后面的双引号里面，保存，然后确定脚本目录下有刚才的POSCAR.old，直接运行这个脚本就能得到设置好了固定原子标志的POSCAR。 
set_POSCAR_optflag.py
(1.03 KB, 下载次数 Times of downloads: 7)
|
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 |
chever 发表于 2024-11-29 15:46 嗯嗯，都在推荐vaspkit，那我得下载了！ |
| 试试vaspkit |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-18 07:40 , Processed in 0.488674 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites