24小时热门版块排行榜    

Znn3bq.jpeg
查看: 779  |  回复: 0

caqhero

金虫 (正式写手)

[求助] abinit 几何结构优化无法收敛

如题,最近在用abinit对单层石墨烯进行计算,几何结构优化无法收敛!想在此讨教解决方案。请各位大侠不吝赐教。谢谢!现贴上我的输入文件。
# AGNR
#ndtset 10
#ngkpt 1 1 4
#ngkpt+ 0 0 2
#Definition of the unit cell
acell    19.29 10.0 4.26258  Angstrom #设置晶格常数a=b=c为9.65285 a.u.
angdeg    90 90 90

#Definition of the atom types
ntypat 2          # There is only one type of atom
znucl 6 1         # The keyword "znucl" refers to the atomic number of the
                  # possible type(s) of atom. The pseudopotential(s)
                  # mentioned in the "files" file must correspond
                  # to the type(s) of atom. Here, the only type is Hydrogen.
                        

#Definition of the atoms
natom 18           # There are two atoms
typat 1 1 1 1 1 1 1
      1 1 1 1 1 1 1
      2 2 2 2
                  # They both are of type 1, that is, Hydrogen

xred             # This keyword indicates that the location of the atoms
                  # will follow, one triplet of number for each atom
                  # Triplet giving the cartesian coordinates of atom 1, in Bohr
                  # Triplet giving the cartesian coordinates of atom 2, in Bohr
0.30863122  0.5         0.11111112
0.30863122  0.5         0.22222223
0.37242081  0.5         0.27777778
0.37242081  0.5         0.05555556
0.43621039  0.5         0.11111112
0.43621039  0.5         0.22222223
0.50000001  0.5         0.27777778
0.50000001  0.5         0.05555556
0.56378962  0.5         0.11111112
0.56378962  0.5         0.22222223
0.62757918  0.5         0.27777778
0.62757918  0.5         0.05555556
0.69136879  0.5         0.11111112
0.69136879  0.5         0.22222223
0.25920175  0.5         0.06806207  
0.25920175  0.5         0.26527127  
0.74079826  0.5         0.06806207  
0.74079826  0.5         0.26527127  

kptopt 1 #设置生成k点的方法,这里表明有ngkpt和nshiftk来确定k点网格的大小

ngkpt 1 1 4 #设置对布里渊区进行8x8x8网格的划分

nshiftk 1 #对划分得到的k点按下面的偏移量进行平移

shiftk 0.0 0.0 0.5

#Definition of the SCF procedure
ionmov  3             # Use the modified Broyden algorithm
ntime  500            # Maximum number of Broyden "timesteps"
tolmxf  5.0d-4        # Stopping criterion for the geometry optimization :
ecut 40.0
iscf 7 #自洽迭代中的算法,这里采用cg方法优化。
toldff 5.0d-5 #自洽迭代迭过程总能收敛的标准
diemac 3.0 #体系宏观的介电常数
nstep 500 #最大自洽迭代步数
#getwfk -1
getden -1
#Definition of the outfile information
prtwf 0 #不输出波函数文件。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 caqhero 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290求调剂 +13 luoziheng 2026-04-10 13/650 2026-04-11 23:15 by labixiaoqiao
[考研] 291分调剂 +5 上岸小莹加油 2026-04-09 6/300 2026-04-11 21:06 by 逆水乘风
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 300分求调剂 (085501机械专硕,本科扬大) +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 吉大计算机技术331分,英语六级,求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 278求调剂 +27 范婷娜 2026-04-07 31/1550 2026-04-09 20:49 by zhouxiaoyu
[考研] 一志愿华工085600 331分 +6 天下ww 2026-04-09 6/300 2026-04-09 18:59 by l_paradox
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 288求调剂 +15 没有答案_ 2026-04-05 15/750 2026-04-09 10:22 by 5268321
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
信息提示
请填处理意见