24小时热门版块排行榜    

查看: 6432  |  回复: 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 科研

» 猜你喜欢

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

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

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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等)。
至于收敛性测试,要注意收敛条件是
引用回帖:
0.1meV/原子

,而不是
引用回帖:
0.1meV

在优化几何结构的时候,建议不要使用ISIF=3选项!速度会非常慢。在同样的计算条件下,建议用loop.sh方案。
一般来说,先在1x1x1体系下在实验测得的晶格常数基础上用ISIF=7优化得到体系的晶格常数。
然后作等比例放大,得到nxnxn体系的POSCAR,然后用循环静态计算这个晶格常数附近几个值对应的体系的能量,用murn方程拟合得到体系的平衡态晶格常数。
9楼2012-05-17 00:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

xueht987

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
uuv2010: 金币+1, 多谢指导 2012-05-14 14:14:10
Illusionist: 金币+20, ★★★很有帮助 2012-05-14 14:16:09
首先打破对称性是可以的,只不过计算可能会多费点时间;
至于替位掺杂,要考虑一下掺杂原子与被掺杂原子是否是等价的,若果是非等价掺杂要考虑电荷补偿的问题,相关参数是NELECT,具体关于这方面之前本版有很多的相关帖子可以参考,自己搜一下。
2楼2012-05-14 14:13:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

引用回帖:
2楼: Originally posted by xueht987 at 2012-05-14 14:13:29:
首先打破对称性是可以的,只不过计算可能会多费点时间;
至于替位掺杂,要考虑一下掺杂原子与被掺杂原子是否是等价的,若果是非等价掺杂要考虑电荷补偿的问题,相关参数是NELECT,具体关于这方面之前本版有很多的 ...

打破对称性之后是否需要提高ENCUT?因为我打破之后,发现晶格常数变大了许多·······3.18变到了3.21左右
3楼2012-05-14 14:15:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueht987

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
贺仪: 金币+2, 谢谢指教! 2012-05-14 14:29:32
Illusionist: 金币+10, ★★★很有帮助, 谢谢,我再看看··有什么问题 2012-05-14 14:36:36
引用回帖:
3楼: Originally posted by Illusionist at 2012-05-14 14:15:57:
打破对称性之后是否需要提高ENCUT?因为我打破之后,发现晶格常数变大了许多·······3.18变到了3.21左右

只要是收敛的,我个人觉得不需要;至于膨胀可与实验值比较一下,一般在1.5%内都可以接受,甚至更大。
4楼2012-05-14 14:25:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +16 Lanmanbaby 2026-08-09 26/1300 2026-08-11 17:16 by 天神眷顾
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-11 3/150 2026-08-11 17:15 by CLKTSgjM0Uqj
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +6 布布和一二 2026-08-10 7/350 2026-08-11 16:57 by Ai19911123
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +7 Tide man 2026-08-10 8/400 2026-08-11 15:49 by 铁帽子农民
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 6/300 2026-08-11 13:19 by JACS/Angew
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
信息提示
请填处理意见