當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

justinzhouke

新蟲 (初入文壇)

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

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

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

   The maximum number of steps has been reached.

這個(gè)怎么處理啊。
求大神幫忙。

............
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ù)此樓

» 猜你喜歡

309求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)

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

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

whoisyourdady

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

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

liqizuiyang

木蟲 (著名寫手)

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

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

引用回帖:

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

查看全部 9 個(gè)回答

liqizuiyang

木蟲 (著名寫手)

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

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

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

0.5巴的壓強(qiáng)……這個(gè)值低了，難怪不收斂。
默認(rèn)的是0.5，個(gè)人感覺(jué)0.1不能再低了。

贊一下

回復(fù)此樓

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

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

justinzhouke

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線: 6.5小時(shí)
蟲號(hào): 3513863
注冊(cè): 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巴的壓強(qiáng)……這個(gè)值低了，難怪不收斂。
默認(rèn)的是0.5，個(gè)人感覺(jué)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的回帖

justinzhouke

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線: 6.5小時(shí)
蟲號(hào): 3513863
注冊(cè): 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

沒(méi)有出現(xiàn) The maximum number of steps has been reached.可以沒(méi)有達(dá)到應(yīng)力為0.
當(dāng)我改為0.01后還是出現(xiàn) 沒(méi)有出現(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ā)表

[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +10	jhybd 2026-03-23	15/750	2026-03-29 21:49 by nanaliuyun
[考研] 337求調(diào)劑 +4	《樹》 2026-03-29	4/200	2026-03-29 20:52 by 唐沐兒
[考研] 297求調(diào)劑 +11	田洪有 2026-03-26	11/550	2026-03-29 13:14 by peike
[考研] 348求調(diào)劑 +5	小懶蟲不懶了 2026-03-28	5/250	2026-03-29 10:34 by 唐沐兒
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 11408軟件工程求調(diào)劑 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 21:50 by zhq0425
[考研] 317求調(diào)劑 +6	十閑wx 2026-03-24	6/300	2026-03-28 13:27 by Iveryant
[考研] 292求調(diào)劑 +14	鵝鵝鵝額額額額?/a> 2026-03-25	15/750	2026-03-28 08:45 by WYUMater
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 化學(xué)調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 272求調(diào)劑 +7	腳滑的守法公民 2026-03-27	7/350	2026-03-27 17:23 by laoshidan
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 調(diào)劑求收留 +7	果然有我 2026-03-26	7/350	2026-03-27 00:26 by wxiongid
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt