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

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

24小時(shí)熱門(mén)版塊排行榜

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

返回列表

justinzhouke

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線(xiàn): 6.5小時(shí)
蟲(chóng)號(hào): 3513863
注冊(cè): 2014-11-02
性別: GG
專(zhuān)業(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ù)此樓

» 猜你喜歡

數(shù)一英一271專(zhuān)碩（085401）求調(diào)劑，可跨已經(jīng)有8人回復(fù)
295求調(diào)劑已經(jīng)有4人回復(fù)
0856求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有10人回復(fù)
299求調(diào)劑已經(jīng)有6人回復(fù)
一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑已經(jīng)有6人回復(fù)
289求調(diào)劑已經(jīng)有12人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑，求導(dǎo)師收已經(jīng)有10人回復(fù)
各位老師好，我的一志愿為北京科技大學(xué)085601材料專(zhuān)碩已經(jīng)有5人回復(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ò)，菜鳥(niǎo)求助已經(jīng)有4人回復(fù)
有關(guān)QE(Pwscf） PHonon編譯計(jì)算的一些基本問(wèn)題已經(jīng)有10人回復(fù)
QE(Pwscf) 非共線(xiàn)磁性計(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ì)稱(chēng)性操作警告已經(jīng)有43人回復(fù)
求專(zhuān)家解惑：pwscf并行計(jì)算cpu分配問(wèn)題已經(jīng)有3人回復(fù)
pwscf 安裝測(cè)試出問(wèn)題。。求救已經(jīng)有2人回復(fù)
pwscf畫(huà)能帶圖求助已經(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

木蟲(chóng) (著名寫(xiě)手)

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

★
小木蟲(chóng): 金幣+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

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線(xiàn): 6.5小時(shí)
蟲(chóng)號(hào): 3513863
注冊(cè): 2014-11-02
性別: GG
專(zhuān)業(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的回帖

liqizuiyang

木蟲(chóng) (著名寫(xiě)手)

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

★
小木蟲(chóng): 金幣+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。直接寫(xiě)0.1，不要換算。

贊一下

回復(fù)此樓

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

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

justinzhouke

新蟲(chóng) (初入文壇)

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

引用回帖:

4樓: Originally posted by liqizuiyang at 2015-08-28 11:18:43
這里的單位是kBar，不是Bar。直接寫(xiě)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的回帖

liqizuiyang

木蟲(chóng) (著名寫(xiě)手)

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

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

引用回帖:

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)呢？我見(jiàn)過(guò)的文獻(xiàn)中有用0.5kBar的。

贊一下

回復(fù)此樓

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

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

justinzhouke

新蟲(chóng) (初入文壇)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

whoisyourdady

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

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

茶味果凍

新蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 265.8
散金: 5
帖子: 78
在線(xiàn): 110小時(shí)
蟲(chóng)號(hào): 2793749
注冊(cè): 2013-11-11
專(zhuān)業(yè): 電化學(xué)分析

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

樓主我問(wèn)個(gè)問(wèn)題，像力
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
這怎么解釋?zhuān)瑸樯妒侨齻€(gè)值，還有怎么檢查各個(gè)原子上的力，下面是我計(jì)算relax結(jié)果的力這個(gè)具體是什么意思啊 3個(gè)值
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

金蟲(chóng) (小有名氣)

應(yīng)助: 27 (小學(xué)生)
金幣: 713
散金: 790
紅花: 4
帖子: 225
在線(xiàn): 129.5小時(shí)
蟲(chóng)號(hào): 1760414
注冊(cè): 2012-04-17
專(zhuān)業(yè): 物理學(xué)I

多大應(yīng)力算大

發(fā)自小木蟲(chóng)Android客戶(hù)端

回復(fù)此樓

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

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

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

返回列表

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

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

[考研] 083000學(xué)碩274求調(diào)劑 +8	Li李魚(yú) 2026-03-26	8/400	2026-03-28 20:33 by 加油向未來(lái)啊
[考研] 295材料工程專(zhuān)碩求調(diào)劑 +7	1428151015 2026-03-27	7/350	2026-03-28 19:58 by S240
[考研] 調(diào)劑 +3	好好讀書(shū)。 2026-03-28	3/150	2026-03-28 12:04 by 王保杰33
[考研] 材料與化工考研調(diào)劑 +17	孅華 2026-03-22	17/850	2026-03-28 08:35 by WYUMater
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 266求調(diào)劑 +11	陽(yáng)陽(yáng)哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 求調(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
[考研] 生物學(xué) 296 求調(diào)劑 +4	朵朵- 2026-03-26	6/300	2026-03-26 19:01 by 不吃魚(yú)的貓
[考研] 340求調(diào)劑 +3	Amber00 2026-03-26	3/150	2026-03-26 18:57 by 不吃魚(yú)的貓
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 0854電子信息求調(diào)劑 +7	α____ 2026-03-22	9/450	2026-03-25 13:37 by α____
[考研] 086003食品工程求調(diào)劑 +6	淼淼111 2026-03-24	6/300	2026-03-25 10:29 by 3Strings
[考研] B區(qū)考研調(diào)劑 +4	yqdszhdap－ 2026-03-22	5/250	2026-03-25 08:51 by baoball
[考研] 求調(diào)劑 +6	研研，接電話(huà) 2026-03-24	7/350	2026-03-24 17:01 by barlinike
[考研] 333求調(diào)劑 +3	ALULU4408 2026-03-23	3/150	2026-03-23 19:04 by macy2011
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋?xiě)?/a>
[考研] 306求調(diào)劑 +5	來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320

24小時(shí)熱門(mén)版塊排行榜

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

» 猜你喜歡

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

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

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