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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » PWscf(QE)的结构优化问题
91/1返回列表
查看: 4319  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

justinzhouke

新虫 (初入文坛)

[交流] PWscf(QE)的结构优化问题 已有3人参与

我在用vc-relax做结构优化的时候，还没有达到应力收敛条件的时候就退出了

     The maximum number of steps has been reached.

这个怎么处理啊。
求大神帮忙。

............
    Forces acting on atoms (Ry/au):

     atom    1 type  1   force =     0.00000000    0.00000000    0.00000000
     atom    2 type  2   force =     0.00000000    0.00000000    0.00000000

     Total force =     0.000000     Total SCF correction =     0.000000


     entering subroutine stress ...

          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.02
  -0.00000013  -0.00000000   0.00000000         -0.02     -0.00      0.00
   0.00000000  -0.00000013  -0.00000000          0.00     -0.02     -0.00
   0.00000000  -0.00000000  -0.00000013          0.00     -0.00     -0.02


     The maximum number of steps has been reached.

     End of BFGS Geometry Optimization
Begin final coordinates
     new unit-cell volume =    337.13304 a.u.^3 (    49.95796 Ang^3 )

CELL_PARAMETERS (alat= 11.69699086)
   0.000000000   0.472261950   0.472261950
   0.472261950  -0.000000000   0.472261950
   0.472261950   0.472261950  -0.000000000

ATOMIC_POSITIONS (crystal)
In      -0.000000000  -0.000000000   0.000000000
As       0.250000000   0.250000000   0.250000000
End final coordinates
............

输入文件如下：
&control
        calculation = 'vc-relax',
          restart_mode= 'from_scratch',
          pseudo_dir  = './',
          outdir      = './tmp',
        prefix='inas',
        wf_collect = .true.
        tstress = .true.
        tprnfor = .true.
        etot_conv_thr = 1.0D-4
        forc_conv_thr = 1.0D-3
        nstep = 20000
/
&system
    ibrav=  0,
    celldm(1)=11.6969908603;
    nat=  2,
    ntyp= 2,
    ecutwfc =48,
/
&electrons
    electron_maxstep=2500
    mixing_beta = 0.7
    conv_thr =  5.0d-12
/
&IONS
ion_dynamics = 'bfgs'
/
&CELL
cell_dynamics = 'bfgs'
press=0.0
press_conv_thr = 5.0D-4
/
ATOMIC_SPECIES
In  111.818  In.pz-bhs.UPF
As   74.920  As.pz-bhs.UPF

ATOMIC_POSITIONS crystal
In 0.00 0.00 0.00
As 0.25 0.25 0.25
K_POINTS automatic
    9 9 9 0 0 0
CELL_PARAMETERS alat
   0.000000000   0.5   0.5
   0.5  -0.000000000   0.5
   0.5   0.5   0.000000000
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
whoisyourdady
1楼 2015-08-27 22:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
press_conv_thr = 5.0D-4
===============
0.5巴的压强……这个值低了，难怪不收敛。
默认的是0.5，个人感觉0.1不能再低了。
赞一下 回复此楼
2楼2015-08-28 09:10:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

justinzhouke

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by liqizuiyang at 2015-08-28 09:10:56
press_conv_thr = 5.0D-4
===============
0.5巴的压强……这个值低了，难怪不收敛。
默认的是0.5，个人感觉0.1不能再低了。

你好，
我把press_conv_thr 改为1.0D-4和1.0D-5后，还是出现

     The maximum number of steps has been reached.

其中1.0D-5的部分输出结果为：

     Forces acting on atoms (Ry/au):

     atom    1 type  1   force =     0.00000000    0.00000000    0.00000000
     atom    2 type  2   force =     0.00000000    0.00000000    0.00000000

     Total force =     0.000000     Total SCF correction =     0.000000


     entering subroutine stress ...

          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.02
  -0.00000013  -0.00000000   0.00000000         -0.02     -0.00      0.00
   0.00000000  -0.00000013  -0.00000000          0.00     -0.02     -0.00
   0.00000000  -0.00000000  -0.00000013          0.00     -0.00     -0.02


     The maximum number of steps has been reached.
赞一下 回复此楼
whoisyourdady
3楼2015-08-28 10:47:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by justinzhouke at 2015-08-28 10:47:04
你好，
我把press_conv_thr 改为1.0D-4和1.0D-5后，还是出现

     The maximum number of steps has been reached.

其中1.0D-5的部分输出结果为：

     Forces acting on atoms (Ry/au):

     atom   ...

这里的单位是kBar，不是Bar。直接写0.1，不要换算。
赞一下 回复此楼
4楼2015-08-28 11:18:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

justinzhouke

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by liqizuiyang at 2015-08-28 11:18:43
这里的单位是kBar，不是Bar。直接写0.1，不要换算。...

我改为0.1后，优化后的结果为：

     Total force =     0.000000     Total SCF correction =     0.000000


     entering subroutine stress ...

          total   stress  (Ry/bohr**3)                   (kbar)     P=   -0.05
  -0.00000031  -0.00000000   0.00000000         -0.05     -0.00      0.00
  -0.00000000  -0.00000031   0.00000000         -0.00     -0.05      0.00
   0.00000000  -0.00000000  -0.00000031          0.00     -0.00     -0.05


     bfgs converged in   6 scf cycles and   5 bfgs steps
     (criteria: energy < 0.10E-03, force < 0.10E-02, cell < 0.10E+00)

     End of BFGS Geometry Optimization

     Final enthalpy =     -16.9554247603 Ry
Begin final coordinates
     new unit-cell volume =    337.14584 a.u.^3 (    49.95986 Ang^3 )

CELL_PARAMETERS (alat= 11.69699086)
   0.000000000   0.472267929   0.472267929
   0.472267929  -0.000000000   0.472267929
   0.472267929   0.472267929  -0.000000000

ATOMIC_POSITIONS (crystal)
In       0.000000000   0.000000000  -0.000000000
As       0.250000000   0.250000000   0.250000000
End final coordinates

没有出现 The maximum number of steps has been reached.可以没有达到应力为0.
当我改为0.01后还是出现 没有出现 The maximum number of steps has been reached.
赞一下 回复此楼
whoisyourdady
5楼2015-08-28 11:51:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
5楼: Originally posted by justinzhouke at 2015-08-28 11:51:53
我改为0.1后，优化后的结果为：

     Total force =     0.000000     Total SCF correction =     0.000000


     entering subroutine stress ...

          total   stress  (Ry/bohr**3)            ...

为什么要用这么低的压强呢？我见过的文献中有用0.5kBar的。
赞一下 回复此楼
6楼2015-08-28 13:44:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

justinzhouke

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by liqizuiyang at 2015-08-28 13:44:54
为什么要用这么低的压强呢？我见过的文献中有用0.5kBar的。...

压强用的高，收敛时候应力还是很大
赞一下 回复此楼
whoisyourdady
7楼2015-08-28 15:03:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

茶味果冻

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主我问个问题，像力
Forces acting on atoms (Ry/au):

      atom    1 type  1   force =     0.00000000    0.00000000    0.00000000
      atom    2 type  2   force =     0.00000000    0.00000000    0.00000000
这怎么解释，为啥是三个值，还有怎么检查各个原子上的力，下面是我计算relax结果的力这个具体是什么意思啊 3个值
Forces acting on atoms (au):
   Ti  -0.602690E-03  0.299794E-02 -0.103070E+00
   Ti   0.114534E-02 -0.280239E-02 -0.235120E-01
   Ti   0.154898E-02  0.768104E-02  0.230297E-01
   Ti  -0.305820E-02 -0.297500E-02  0.888280E-01
赞一下 回复此楼
8楼2015-08-29 15:00:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyyu200

金虫 (小有名气)
多大应力算大

发自小木虫Android客户端
回复此楼
9楼2015-09-08 07:15:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 justinzhouke 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
关闭
请填处理意见
关闭