版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

justinzhouke

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 3.5
红花: 1
帖子: 24
在线: 6.5小时
虫号: 3513863
注册: 2014-11-02
性别: GG
专业: 凝聚态物性I:结构、力学和

[交流] 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

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

pwscf结构优化报错问题求助已经有0人回复
quantum-espresso（原先叫pwscf）安装问题已经有8人回复
QE_GPU加速的使用范围已经有1人回复
QE PWscf pw-gpu 计算速度终止条件收敛速度已经有0人回复
Quantum-ESPRESSO 5.0.3，Windows 32-bit，串行版已经有55人回复
PWscf终端运行报错，菜鸟求助已经有4人回复
有关QE(Pwscf） PHonon编译计算的一些基本问题已经有10人回复
QE(Pwscf) 非共线磁性计算已经有0人回复
pwscf-5.0得到的声子谱很乱是怎么回事？已经有8人回复
Problem - nscf - c_bands: eigenvalues not converged 已经有20人回复
pwscf中average.x的使用已经有0人回复
紧急求助：运行pwscf出现对称性操作警告已经有43人回复
求专家解惑：pwscf并行计算cpu分配问题已经有3人回复
pwscf 安装测试出问题。。求救已经有2人回复
pwscf画能带图求助已经有6人回复
pwscf计算integrated charge出错及不收敛已经有13人回复

whoisyourdady

1楼 2015-08-27 22:18:15

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

justinzhouke

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 3.5
红花: 1
帖子: 24
在线: 6.5小时
虫号: 3513863
注册: 2014-11-02
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

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的回帖

查看全部 9 个回答

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+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

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 3.5
红花: 1
帖子: 24
在线: 6.5小时
虫号: 3513863
注册: 2014-11-02
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

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

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+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的回帖

查看全部 9 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 基金系统什么内容也没有 30+3	winsaint 2026-08-27	6/300	2026-08-28 01:46 by alongwaytogo
[基金申请] 投票: 有多少人是今天查系统知道结果的？ +15	爱看书的可乐 2026-08-26	17/850	2026-08-28 01:23 by xxniao123
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	XLGfUIdOmEAO 2026-08-27	3/150	2026-08-28 00:49 by 3RAFyYXYSwGS
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7	苏东坡二世 2026-08-23	7/350	2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6	lyz123lyz 2026-08-27	7/350	2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中，共勉 +10	eulota 2026-08-26	10/500	2026-08-27 17:23 by lqllinqiaoli
[基金申请] 怎么看青基中了没有啊 +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 北京莱茵编辑
[基金申请] 系统查不到 +10	董八千 2026-08-26	10/500	2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5	huang1991js 2026-08-26	5/250	2026-08-26 16:09 by Equinoxhua
[基金申请] 为什么国自然不能直接公布 +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
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6	wittyboy 2026-08-26	14/700	2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3	paperzjh 2026-08-26	3/150	2026-08-26 10:41 by LemmonTr
[基金申请] 牛来！米来！面来！ +8	beefly 2026-08-26	8/400	2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了！？ +6	chengyan1220 2026-08-23	6/300	2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9	图啦图啦 2026-08-24	10/500	2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11	fansofjerry 2026-08-21	11/550	2026-08-25 11:04 by 孤独的英雄6
[基金申请] 应该是下周三26日公布了吧？ +4	哈哈蛤？ 2026-08-21	4/200	2026-08-21 10:58 by Vivilian

24小时热门版块排行榜

justinzhouke

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

justinzhouke

liqizuiyang

justinzhouke

liqizuiyang

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