24小时热门版块排行榜    

Znn3bq.jpeg
查看: 5861  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +5 江淮北月 2026-04-05 5/250 2026-04-10 16:26 by 高维春
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 生物与医药273求调剂 +18 荔题南墙 2026-04-05 19/950 2026-04-10 08:14 by kangsm
[考研] 087100初试311求调剂 +3 任雅琴 2026-04-09 3/150 2026-04-09 22:42 by lbsjt
[考研] 085400 328分 求调剂 +8 喂你一个大橙子 2026-04-09 12/600 2026-04-09 19:20 by zl8213662
[考研] 材料299专硕求调剂 +10 +21 2026-04-09 10/500 2026-04-09 17:34 by 1753564080
[硕博家园] 新一代电子信息294求调剂 不挑学校 +5 Ytyt11 2026-04-09 6/300 2026-04-09 14:40 by Ytyt11
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 求助071001调剂!!! +7 黄守松 2026-04-05 8/400 2026-04-09 09:07 by 徐良白眉大侠
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 283分求调剂 +14 试试看呗 2026-04-04 14/700 2026-04-08 07:03 by lijunpoly
[考研] 304求调剂(085602,过四级,一志愿985) +25 化工人999 2026-04-04 26/1300 2026-04-07 22:06 by hemengdong
[考研] 318求调剂 +5 李青山山山 2026-04-07 5/250 2026-04-07 18:24 by 蓝云思雨
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干 +3 1021075758 2026-04-06 4/200 2026-04-07 09:23 by 蓝云思雨
[考研] 考研调剂 +7 15615482637 2026-04-04 7/350 2026-04-06 22:56 by chenzhimin
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
信息提示
请填处理意见