| 查看: 6514 | 回复: 15 | ||||
Illusionist银虫 (正式写手)
|
[求助]
如何做vasp掺杂
|
|||
|
本人想用vasp做掺杂计算，刚开始学 1.直接建立3x3x3的超晶格模型 1）关于K的选择 3x3x3 -1019.397709 4x4x4 -1019.397873 所以选择了3x3x3 2）ENCUT 280 4093.144030 300 4093.144030 320 4093.144030 340 4093.144030 这里我最后选择了500ENCUT去进行结构优化。 2.结构优化INCAR： SYSTEM=SnO2 ENCUT=500 ISMEAR=-5 GGA=91 VOSKOWN=1 ISYM=0 ISIF=3 NSW=100 PTIOM=0.1 IBRION=2 这里我选择打破对称性不知道是否合理 ISYM=0设置前后影响挺大，起码晶格常数会变化。LREAL=auto设置后，晶格常数也会有明显变化，请问一下一般如何处理这样的问题 另外我想在这一步结束后掺入替位的Sb，但是Sb有+3和+5价之分，该如何设定价态呢 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP | 科研 |
» 猜你喜欢
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于TiO2掺杂DOS图的问题
已经有12人回复
- vasp掺杂的问题讨论
已经有9人回复
- 看看别人仅用VASP做DOS计算发好文章（两篇PRL）
已经有30人回复
- 掺杂求高手指点：vasp 不同掺杂浓度模型设置。
已经有8人回复
- vasp 结构优化问题
已经有14人回复
- 如何处理在体结构中掺杂小原子（如B，N，C等）并做结构弛豫？
已经有39人回复
- vasp计算中如何比较掺杂前后结合能
已经有13人回复
- vasp计算掺杂体系的形成能
已经有7人回复
- 【求助】设置铁磁和反铁磁结果一样???
已经有5人回复
- 【求助】如何设定带分数电荷的H，用于钝化表面的悬挂键
已经有10人回复
7楼2012-05-14 20:57:00
5楼2012-05-14 14:32:59
dxf_dxf_dxf
银虫 (小有名气)
- 1ST强帖: 4
- 应助: 15 (小学生)
- 金币: 269.1
- 红花: 3
- 帖子: 104
- 在线: 29.1小时
- 虫号: 1449737
- 注册: 2011-10-19
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★
fzx2008: 金币+3, 谢谢指教 2012-05-17 14:57:51
fzx2008: 金币+3, 谢谢指教 2012-05-17 14:57:51
|
先在1x1x1体系下用ISIF=7优化得到粗略的晶格常数 然后做ENCUT和KPOINTS的收敛性测试 按满足收敛条件的ENCUT和KPOINTS值，计算不同晶格常数下体系的能量 用murn.f拟合得到平衡态晶格常数。 把1x1x1体系按比例放大，得到2x2x2/3x3x3/4x4x4等体系。 对3x3x3体系： KPOINTS设成 然后进行ENCUT收敛性测试，记录不同ENCUT情况下能量和cpu time 然后在4x4x4x体系下做同样的事情，记录不同ENCUT情况下能量和cpu time 再重复至5x5x5,6x6x6,7x7x7 计算时间超过48小时的任务就可以忽略了。 综合考虑性能和精度的影响，确定最后选用的掺杂体系，例如我就选用4x4x4体系。 在此体系下，加入掺杂原子，按照之前收敛性测试的结果设置ENCUT 接下来就开始几何优化了。这里一般要求保持原胞的形状和体积，不过你也可以先优化原胞的体积，一般不要优化原胞的形状。不过这也说不定，例如极高压情况中体系的晶格会改变(例如从体心立方变成六角密堆，所以原胞形状是否改变跟你所做的具体工作有关。) 优化好原胞以后再做离子迟豫。你可以自己写脚本计算，也可以让vasp自己算。 |
12楼2012-05-17 13:44:14
xueht987
木虫 (正式写手)
- 应助: 112 (高中生)
- 贵宾: 0.081
- 金币: 3514.2
- 散金: 106
- 红花: 6
- 帖子: 593
- 在线: 273.4小时
- 虫号: 1059317
- 注册: 2010-07-17
- 性别: GG
- 专业: 金属材料表面科学与工程
2楼2012-05-14 14:13:29
Illusionist
银虫 (正式写手)
- 1ST强帖: 1
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
3楼2012-05-14 14:15:57
xueht987
木虫 (正式写手)
- 应助: 112 (高中生)
- 贵宾: 0.081
- 金币: 3514.2
- 散金: 106
- 红花: 6
- 帖子: 593
- 在线: 273.4小时
- 虫号: 1059317
- 注册: 2010-07-17
- 性别: GG
- 专业: 金属材料表面科学与工程
4楼2012-05-14 14:25:10
Illusionist
银虫 (正式写手)
- 1ST强帖: 1
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
6楼2012-05-14 14:40:13
liu108
木虫 (正式写手)
- 应助: 125 (高中生)
- 金币: 1753.7
- 红花: 6
- 帖子: 362
- 在线: 122.9小时
- 虫号: 1210273
- 注册: 2011-02-23
- 性别: GG
- 专业: 理论和计算化学
8楼2012-05-14 21:03:38
dxf_dxf_dxf
银虫 (小有名气)
- 1ST强帖: 4
- 应助: 15 (小学生)
- 金币: 269.1
- 红花: 3
- 帖子: 104
- 在线: 29.1小时
- 虫号: 1449737
- 注册: 2011-10-19
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教 2012-05-17 08:38:22
Illusionist: 金币+20, ★★★很有帮助, Gamma计算什么，结构优化还是收敛测试？对于K空间的gamma点对应的应该是是空间的晶格边界还是中心，掺杂的时候最好不要放在边间还是中心？ 2012-05-17 08:48:09
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指教 2012-05-17 08:38:22
Illusionist: 金币+20, ★★★很有帮助, Gamma计算什么，结构优化还是收敛测试？对于K空间的gamma点对应的应该是是空间的晶格边界还是中心，掺杂的时候最好不要放在边间还是中心？ 2012-05-17 08:48:09
|
掺杂的问题，建议把K格点设成只有Gamma点，再用更大的晶格(例如4x4x4,5x5x5,6x6x6等)。 至于收敛性测试，要注意收敛条件是 ，而不是 在优化几何结构的时候，建议不要使用ISIF=3选项！速度会非常慢。在同样的计算条件下，建议用loop.sh方案。 一般来说，先在1x1x1体系下在实验测得的晶格常数基础上用ISIF=7优化得到体系的晶格常数。 然后作等比例放大，得到nxnxn体系的POSCAR，然后用循环静态计算这个晶格常数附近几个值对应的体系的能量，用murn方程拟合得到体系的平衡态晶格常数。 |
9楼2012-05-17 00:07:44
dxf_dxf_dxf
银虫 (小有名气)
- 1ST强帖: 4
- 应助: 15 (小学生)
- 金币: 269.1
- 红花: 3
- 帖子: 104
- 在线: 29.1小时
- 虫号: 1449737
- 注册: 2011-10-19
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2012-05-17 00:09:49
