24小时热门版块排行榜    

查看: 769  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
信息提示
请填处理意见