24小时热门版块排行榜    

查看: 529  |  回复: 6
当前主题已经存档。

liuhuan2017

新虫 (小有名气)

[交流] 【求助】关于晶格优化

我对GaAs进行晶格优化,in文件如下,GaAs的实验值是10.68,这个优化出来的结果是7.28.为什么结果相差这么多呢?
ionmov  3
ntime  15
dilatmx  1.3
ecutsm  0.8
prtgeo 1
kptopt 1         
nshiftk 4
shiftk  0.5 0.5 0.5  
        0.5 0.0 0.0
        0.0 0.5 0.0
        0.0 0.0 0.5

ngkpt  6 6 6        
acell 3*10.68         
rprim  0.0  0.5  0.5  
       0.5  0.0  0.5   
       0.5  0.5  0.0
ntypat 2         
znucl 31 33        
natom 2         
typat 1 2        
xred                           
   0.0  0.0  0.0  
   1/4  1/4  1/4  
ecut 8.0         
nstep 10         
toldfe 1.0d-6
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjmaxpayne(金币+1,VIP+0):谢谢 11-18 20:38
ecut 8.0  ?
有点小吧:)

注:我没用过abinit,如果说错了别见怪;)
2楼2009-11-18 20:33:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuchenwf(金币+1,VIP+0):谢谢回帖交流 11-19 20:52
ecut和k-point没收敛。。
3楼2009-11-18 20:41:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guifan

木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuchenwf(金币+2,VIP+0):谢谢回帖交流 11-19 20:52
建议先对主要参数进行测试
参数选好了,才能保证结果的准确性
楼上说的两个参数要测试收敛才行
4楼2009-11-18 22:18:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuchenwf(金币+2,VIP+0):谢谢回帖交流 11-19 20:52
ecut 太小,  k点适当再多一些。
我觉得是不是你用的赝势有问题啊?
再有,你需要给个压力,如果是常压,就不用给了。
真理是一点点接近的!
5楼2009-11-19 17:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy1123

木虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zdhlover(金币+2,VIP+0):欢迎常来交流讨论 11-19 23:56
ecut太小,ngkpt较低,另外一种可能是优化结果给出的是原胞的acell
我在计算MgS-225时出现过这个情况,输出的是原胞的acell
你可以在MS上建晶胞后导出体系的原胞,以原胞的参数做输入,计算结果是一样的
6楼2009-11-20 08:50:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuan2017

新虫 (小有名气)

果然算出来的是原胞的acell啊,谢谢虫友们的关注
7楼2009-11-21 10:12:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuhuan2017 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见