应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 如何做vasp掺杂
162/2返回列表上一页12
查看: 5451  |  回复: 15
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

condensed

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Illusionist: 金币+5, ★有帮助, 108个原子，这是折算的晶格常数 2012-05-17 08:35:43
liliangfang: 金币+1, 谢谢交流 2012-05-17 08:54:31
引用回帖:
3楼: Originally posted by Illusionist at 2012-05-14 14:15:57:
打破对称性之后是否需要提高ENCUT？因为我打破之后，发现晶格常数变大了许多·······3.18变到了3.21左右

发现晶格常数变大了许多·······3.18变到了3.21左右
'''''''''''''''''''''''''''''''''''''''''''''''''''''

胞（原子数）这么小（少）啊，可靠么
赞一下 回复此楼
11楼2012-05-17 00:34:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
fzx2008: 金币+3, 谢谢指教 2012-05-17 14:57:51
引用回帖:
9楼: Originally posted by dxf_dxf_dxf at 2012-05-17 00:07:44:
掺杂的问题，建议把K格点设成只有Gamma点，再用更大的晶格(例如4x4x4,5x5x5,6x6x6等)。
至于收敛性测试，要注意收敛条件是，而不是
在优化几何结构的时候，建议不要使用ISIF=3选项！速度会非常慢。在同样的计算 ...

先在1x1x1体系下用ISIF=7优化得到粗略的晶格常数
然后做ENCUT和KPOINTS的收敛性测试
按满足收敛条件的ENCUT和KPOINTS值，计算不同晶格常数下体系的能量
用murn.f拟合得到平衡态晶格常数。
把1x1x1体系按比例放大，得到2x2x2/3x3x3/4x4x4等体系。
对3x3x3体系：
KPOINTS设成
CODE:
K-Points
0
Monkhorst Pack
1 1 1
0  0  0

然后进行ENCUT收敛性测试，记录不同ENCUT情况下能量和cpu time
然后在4x4x4x体系下做同样的事情，记录不同ENCUT情况下能量和cpu time
再重复至5x5x5,6x6x6,7x7x7
计算时间超过48小时的任务就可以忽略了。
综合考虑性能和精度的影响，确定最后选用的掺杂体系，例如我就选用4x4x4体系。
在此体系下，加入掺杂原子，按照之前收敛性测试的结果设置ENCUT
接下来就开始几何优化了。这里一般要求保持原胞的形状和体积，不过你也可以先优化原胞的体积，一般不要优化原胞的形状。不过这也说不定，例如极高压情况中体系的晶格会改变(例如从体心立方变成六角密堆，所以原胞形状是否改变跟你所做的具体工作有关。)
优化好原胞以后再做离子迟豫。你可以自己写脚本计算，也可以让vasp自己算。
赞一下(1人) 回复此楼
12楼2012-05-17 13:44:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖

★
fzx2008: 金币+1, 谢谢指教 2012-05-17 14:58:02
引用回帖:
11楼: Originally posted by condensed at 2012-05-17 00:34:06:
发现晶格常数变大了许多·······3.18变到了3.21左右
'''''''''''''''''''''''''''''''''''''''''''''''''''''

胞（原子数）这么小（少）啊，可靠么

3.18变到3.21，变化达到9%。一般来说，LDA近似和GGA近似计算得到的晶格常数和实验结果的差别在5%左右。
赞一下 回复此楼
13楼2012-05-17 13:49:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanjiayu0116

金虫 (初入文坛)
引用回帖:
5楼: Originally posted by 贺仪 at 2012-05-14 14:32:59
首先对你的数据：
2）ENCUT
280   4093.144030
300   4093.144030
320   4093.144030
340   4093.144030
表示怀疑

其次对称性没必要打破，设置isif=3就可以了

楼主可不可以把你测试ENCUT的脚本给我一个
赞一下 回复此楼
14楼2013-03-26 16:11:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

随便注册

银虫 (初入文坛)
LZ，我也想要你的ENCUT测试脚本。当伸手党了，不好意思，非常感谢您！
tom4587@hotmail.com
赞一下 回复此楼
15楼2013-03-29 08:16:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

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

这个电荷补偿一定要设置吗？本身不就是电中性的吗
赞一下 回复此楼
希望有一天成为大神
16楼2013-10-21 15:58:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Illusionist 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭