24小时热门版块排行榜    

查看: 768  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
信息提示
请填处理意见