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

登录注册

应《网络安全法》要求，自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

回复此楼

» 猜你喜欢

求调剂推荐已经有7人回复
289 分105500药学专硕求调剂(找B区学校) 已经有4人回复
0854求调剂已经有13人回复
初试324 中药学一志愿天中医求调剂已经有3人回复
药学求调剂已经有14人回复
327求调剂已经有27人回复
急需调剂已经有5人回复
300求调剂已经有7人回复
271求调剂已经有33人回复
273求调剂已经有8人回复

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

pwscf结构优化报错问题求助已经有0人回复
第一次运用QE计算貌似出现问题，求解答已经有12人回复
QE（pwscf）计算计算声子谱问题已经有1人回复
Error in routine read_namelists (1): reading namelist ions 求助已经有3人回复
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 个回答

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

[考研] 291 求调剂 +40	化工2026届毕业� 2026-04-09	41/2050	2026-04-15 21:54 by noqvsozv
[考研] 一志愿A区211，22408 321求调剂 +6	随心所欲☆ 2026-04-15	7/350	2026-04-15 21:45 by lbsjt
[考研] 085801电气专硕272求调剂 +19	电气李 2026-04-13	21/1050	2026-04-15 13:37 by 黑科技矿业
[考研] 245求调剂 +6	冰糖橘?汽水 2026-04-13	10/500	2026-04-14 10:49 by jyl0317
[考研] 085404 298分求调剂 +11	呼啦呼啦呼呼呼 2026-04-10	12/600	2026-04-14 08:38 by wfj257
[考研] 302求调剂 +10	易！? 2026-04-13	10/500	2026-04-13 19:04 by lbsjt
[考研] 一志愿211 0703化学 346分求调剂 +26	土豆er? 2026-04-09	29/1450	2026-04-13 15:15 by 独醉梦孤城
[考研] 346分，工科0854求调剂，专硕 +6	moser233 2026-04-12	7/350	2026-04-12 22:11 by fqwang
[考研] 296求调剂 +14	汪！？！ 2026-04-10	16/800	2026-04-12 10:48 by zhouyuwinner
[考研] 求调剂，一志愿大连理工大学354分 +5	雨声余生 2026-04-11	6/300	2026-04-11 16:12 by 雨声余生
[考研] 求调剂 +6	archer.. 2026-04-09	8/400	2026-04-11 10:55 by zhq0425
[考研] 288求调剂 +15	代fish 2026-04-09	16/800	2026-04-11 10:26 by wwj2530616
[考研] 281求调剂 +11	觉得好的吧 2026-04-10	11/550	2026-04-11 09:35 by 逆水乘风
[考研] 本科211 工科085400 280分求调剂可跨专业 +11	LZH（等待调剂中 2026-04-10	11/550	2026-04-11 08:39 by zhq0425
[考研] 考研调剂 +26	硕星赴 2026-04-09	27/1350	2026-04-10 22:24 by 猪会飞
[考研] 中科院总分315求调剂 +8	lallalh 2026-04-09	8/400	2026-04-10 19:30 by dick_runner
[考研] 初试261 +3	Asht少 2026-04-10	6/300	2026-04-10 16:38 by Asht少
[考研] 085601初试330分找调剂 +10	流心奶黄包l 2026-04-09	10/500	2026-04-10 08:14 by Sammy2
[考研] 332，085601求调剂 +12	ydfyh 2026-04-09	14/700	2026-04-09 17:28 by wp06
[考研] 化学工程与技术专业一志愿哈工程 291分B区国家级大创负责人有一作论文 +13	Emmy~ 2026-04-09	13/650	2026-04-09 14:47 by only周

24小时热门版块排行榜

justinzhouke

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

» 猜你喜欢

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

justinzhouke

liqizuiyang

justinzhouke

liqizuiyang

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