版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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)化的時候，還沒有達(dá)到應(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ù)此樓

» 猜你喜歡

0856求調(diào)劑已經(jīng)有13人回復(fù)
0703化學(xué) 已經(jīng)有11人回復(fù)
329求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有8人回復(fù)
0856材料化工調(diào)劑總分330 已經(jīng)有12人回復(fù)
309求調(diào)劑已經(jīng)有3人回復(fù)
070300化學(xué)354求調(diào)劑已經(jīng)有4人回復(fù)
330分求調(diào)劑已經(jīng)有3人回復(fù)
299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)

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

pwscf結(jié)構(gòu)優(yōu)化報錯問題求助已經(jīng)有0人回復(fù)
第一次運(yùn)用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終端運(yùn)行報錯，菜鳥求助已經(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ù)
緊急求助：運(yùn)行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的回帖

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巴的壓強(qiáng)……這個值低了，難怪不收斂。
默認(rèn)的是0.5，個人感覺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小時
蟲號: 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巴的壓強(qiáng)……這個值低了，難怪不收斂。
默認(rèn)的是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的回帖

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.可以沒有達(dá)到應(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的回帖

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, 給個紅包，謝謝回帖

引用回帖:

5樓: Originally posted by justinzhouke at 2015-08-28 11:51:53
我改為0.1后，優(yōu)化后的結(jié)果為：

   Total force =    0.000000    Total SCF correction =    0.000000

   entering subroutine stress ...

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

為什么要用這么低的壓強(qiáng)呢？我見過的文獻(xiàn)中有用0.5kBar的。

贊一下

回復(fù)此樓

6樓2015-08-28 13:44:54

已閱回復(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é)和

引用回帖:

6樓: Originally posted by liqizuiyang at 2015-08-28 13:44:54
為什么要用這么低的壓強(qiáng)呢？我見過的文獻(xiàn)中有用0.5kBar的。...

壓強(qiáng)用的高，收斂時候應(yīng)力還是很大

贊一下

回復(fù)此樓

whoisyourdady

7樓2015-08-28 15:03:46

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

茶味果凍

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 265.8
散金: 5
帖子: 78
在線: 110小時
蟲號: 2793749
注冊: 2013-11-11
專業(yè): 電化學(xué)分析

★
小木蟲: 金幣+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結(jié)果的力這個具體是什么意思啊 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

贊一下

回復(fù)此樓

8樓2015-08-29 15:00:35

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

yyyu200

金蟲 (小有名氣)

應(yīng)助: 27 (小學(xué)生)
金幣: 713
散金: 790
紅花: 4
帖子: 225
在線: 129.5小時
蟲號: 1760414
注冊: 2012-04-17
專業(yè): 物理學(xué)I

多大應(yīng)力算大

發(fā)自小木蟲Android客戶端

回復(fù)此樓

9樓2015-09-08 07:15:39

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 justinzhouke 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

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

[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +4	來財aa 2026-03-25	4/200	2026-03-28 07:41 by 棒棒球手
[考研] 085701環(huán)境工程求調(diào)劑 +9	多久上課 2026-03-27	9/450	2026-03-28 03:58 by fmesaito
[有機(jī)交流] 高溫高壓反應(yīng)求助 10+4	chibby 2026-03-25	4/200	2026-03-27 21:08 by BT20230424
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 考研調(diào)劑 +10	呼呼？~+123456 2026-03-24	10/500	2026-03-27 11:46 by wangjy2002
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 351求調(diào)劑 +4	麥克阿磊 2026-03-24	4/200	2026-03-27 00:32 by wxiongid
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 281求調(diào)劑 +3	亞克西good 2026-03-26	5/250	2026-03-26 19:48 by 不吃魚的貓
[考研] 340求調(diào)劑 +3	Amber00 2026-03-26	3/150	2026-03-26 18:57 by 不吃魚的貓
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 0854電子信息求調(diào)劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 070300，一志愿北航320求調(diào)劑 +3	Jerry0216 2026-03-22	5/250	2026-03-23 09:16 by 。。堂堂

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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