justinzhouke

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線: 6.5小時
蟲號: 3513863
注冊: 2014-11-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[交流] PWscf(QE)的結(jié)構(gòu)優(yōu)化問題已有3人參與

我在用vc-relax做結(jié)構(gòu)優(yōu)化的時候，還沒有達到應(yīng)力收斂條件的時候就退出了

   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

回復(fù)此樓

» 猜你喜歡

一志愿華東師范大學(xué)有機化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! 已經(jīng)有4人回復(fù)
298求調(diào)劑已經(jīng)有4人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有11人回復(fù)
0856求調(diào)劑已經(jīng)有7人回復(fù)
356求調(diào)劑已經(jīng)有3人回復(fù)
332求調(diào)劑已經(jīng)有7人回復(fù)
348求調(diào)劑已經(jīng)有4人回復(fù)
305求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有13人回復(fù)
調(diào)劑310 已經(jīng)有9人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

pwscf結(jié)構(gòu)優(yōu)化報錯問題求助已經(jīng)有0人回復(fù)
第一次運用QE計算貌似出現(xiàn)問題，求解答已經(jīng)有12人回復(fù)
QE（pwscf）計算計算聲子譜問題已經(jīng)有1人回復(fù)
Error in routine read_namelists (1): reading namelist ions 求助已經(jīng)有3人回復(fù)
quantum-espresso（原先叫pwscf）安裝問題已經(jīng)有8人回復(fù)
QE_GPU加速的使用范圍已經(jīng)有1人回復(fù)
QE PWscf pw-gpu 計算速度終止條件收斂速度已經(jīng)有0人回復(fù)
Quantum-ESPRESSO 5.0.3，Windows 32-bit，串行版已經(jīng)有55人回復(fù)
PWscf終端運行報錯，菜鳥求助已經(jīng)有4人回復(fù)
有關(guān)QE(Pwscf） PHonon編譯計算的一些基本問題已經(jīng)有10人回復(fù)
QE(Pwscf) 非共線磁性計算已經(jīng)有0人回復(fù)
pwscf-5.0得到的聲子譜很亂是怎么回事？已經(jīng)有8人回復(fù)
Problem - nscf - c_bands: eigenvalues not converged 已經(jīng)有20人回復(fù)
pwscf中average.x的使用已經(jīng)有0人回復(fù)
緊急求助：運行pwscf出現(xiàn)對稱性操作警告已經(jīng)有43人回復(fù)
求專家解惑：pwscf并行計算cpu分配問題已經(jīng)有3人回復(fù)
pwscf 安裝測試出問題。。求救已經(jīng)有2人回復(fù)
pwscf畫能帶圖求助已經(jīng)有6人回復(fù)
pwscf計算integrated charge出錯及不收斂已經(jīng)有13人回復(fù)

whoisyourdady

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

justinzhouke

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線: 6.5小時
蟲號: 3513863
注冊: 2014-11-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

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后，還是出現(xiàn)

   The maximum number of steps has been reached.

其中1.0D-5的部分輸出結(jié)果為：

   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.

贊一下

回復(fù)此樓

whoisyourdady

3樓2015-08-28 10:47:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個回答

liqizuiyang

木蟲 (著名寫手)

應(yīng)助: 59 (初中生)
金幣: 6859.5
散金: 261
紅花: 94
帖子: 2093
在線: 708.9小時
蟲號: 2366686
注冊: 2013-03-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

press_conv_thr = 5.0D-4
===============

0.5巴的壓強……這個值低了，難怪不收斂。
默認的是0.5，個人感覺0.1不能再低了。

贊一下

回復(fù)此樓

2樓2015-08-28 09:10:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liqizuiyang

木蟲 (著名寫手)

應(yīng)助: 59 (初中生)
金幣: 6859.5
散金: 261
紅花: 94
帖子: 2093
在線: 708.9小時
蟲號: 2366686
注冊: 2013-03-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

3樓: Originally posted by justinzhouke at 2015-08-28 10:47:04
你好，
我把press_conv_thr 改為1.0D-4和1.0D-5后，還是出現(xiàn)

   The maximum number of steps has been reached.

其中1.0D-5的部分輸出結(jié)果為：

   Forces acting on atoms (Ry/au):

   atom ...

這里的單位是kBar，不是Bar。直接寫0.1，不要換算。

贊一下

回復(fù)此樓

4樓2015-08-28 11:18:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

justinzhouke

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線: 6.5小時
蟲號: 3513863
注冊: 2014-11-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

4樓: Originally posted by liqizuiyang at 2015-08-28 11:18:43
這里的單位是kBar，不是Bar。直接寫0.1，不要換算。...

我改為0.1后，優(yōu)化后的結(jié)果為：

   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

沒有出現(xiàn) The maximum number of steps has been reached.可以沒有達到應(yīng)力為0.
當(dāng)我改為0.01后還是出現(xiàn) 沒有出現(xiàn) The maximum number of steps has been reached.

贊一下

回復(fù)此樓

whoisyourdady

5樓2015-08-28 11:51:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 309求調(diào)劑 +6	誰不是少年 2026-03-29	6/300	2026-03-29 08:11 by fmesaito
[考研] 調(diào)劑考研 +3	王杰一 2026-03-29	3/150	2026-03-29 08:09 by fmesaito
[考研] 291求調(diào)劑 +3	Y-cap 2026-03-29	3/150	2026-03-29 07:34 by 無際的草原
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9	天天好運來上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 295材料工程專碩求調(diào)劑 +7	1428151015 2026-03-27	7/350	2026-03-28 19:58 by S240
[考研] 275求調(diào)劑 +10	Micky11223 2026-03-25	14/700	2026-03-28 15:48 by Micky11223
[考研] 0703化學(xué)求調(diào)劑 +9	奶油草莓. 2026-03-22	10/500	2026-03-28 13:30 by 唐沐兒
[考研] 調(diào)劑 +3	好好讀書。 2026-03-28	3/150	2026-03-28 12:04 by 王保杰33
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3	熠然好運氣 2026-03-27	3/150	2026-03-28 04:11 by fmesaito
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 086000調(diào)劑 +3	7901117076 2026-03-26	3/150	2026-03-27 21:34 by Jianing_Mi
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 求調(diào)劑一志愿本科北科大化學(xué) 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17	RichLi_ 2026-03-25	17/850	2026-03-26 19:44 by plmuchong
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 機械學(xué)碩總分317求調(diào)劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[有機交流] 有機合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 材料專碩331求調(diào)劑 +4	鮮當(dāng)牛 2026-03-24	4/200	2026-03-24 15:58 by JourneyLucky