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

登錄注冊

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

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

返回列表

didi5158

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 482.4
帖子: 250
在線: 89.9小時(shí)
蟲號: 273449
注冊: 2006-08-20
性別: GG
專業(yè): 無機(jī)材料化學(xué)

[交流] 【求助】優(yōu)化的結(jié)果，scf計(jì)算該結(jié)構(gòu)，force為什么仍然很大？已有7人參與

vc-relax計(jì)算晶體的晶格常數(shù)，用計(jì)算所得到的結(jié)構(gòu)再進(jìn)行自洽計(jì)算，而所得結(jié)果中force，stress，total energy都差別較大，這個(gè)是為什么？請大家指教。

vc-relax計(jì)算結(jié)果：
---------------------------------------------------------------------------------------------------------
total energy             = -588.70418379 Ry
Harris-Foulkes estimate = -588.70418379 Ry
estimated scf accuracy <       3.9E-10 Ry

Total force =    0.000201    Total SCF correction =    0.000059
SCF correction compared to forces is too large, reduce conv_thr

entering subroutine stress ...
      total stress  (Ry/bohr**3)                (kbar)    P= 0.06
0.00000076 0.00000000 0.00000000       0.11    0.00    0.00
0.00000000 0.00000076 0.00000000       0.00    0.11    0.00
0.00000000 0.00000000  -0.00000024       0.00    0.00    -0.03

bfgs converged in  14 scf cycles and  12 bfgs steps
End of BFGS Geometry Optimization

Final enthalpy = -588.7041837859 Ry
new unit-cell volume = 657.99457 a.u.^3 ( 97.50473 Ang^3 )

CELL_PARAMETERS (alat)
0.781116722 0.000000000 0.000000000
  -0.390558361 0.676466924 0.000000000
0.000000000 0.000000000 1.486696131

ATOMIC_POSITIONS (crystal)
Li    0.333333333 0.666666667 0.250000000
Li    -0.333333333  -0.666666667  -0.250000000
F       0.420499391  -0.010919477 0.122948838
F       0.010919477 0.431418866 0.122948838
F    -0.431418866  -0.420499390 0.122948838
F       0.010919477  -0.420499390 0.377051161
F    -0.431418866  -0.010919477 0.377051161
F       0.420499390 0.431418866 0.377051161
F    -0.420499390 0.010919478  -0.122948838
F    -0.010919477  -0.431418866  -0.122948838
F       0.431418866 0.420499390  -0.122948838
F    -0.010919477 0.420499390 0.622948838
F       0.431418866 0.010919477 0.622948838
F    -0.420499390  -0.431418866 0.622948838
Al    0.666666667 0.333333333 0.250000000
Al    -0.666666667  -0.333333333  -0.250000000
Ca    0.000000000 0.000000000 0.000000000
Ca    0.000000000 0.000000000 0.500000000

----------------------------------------------------------------------------------------------------------

用上面所得到結(jié)構(gòu)，scf計(jì)算結(jié)果：

----------------------------------------------------------------------------------------------------------
total energy             = -571.49242692 Ry
Harris-Foulkes estimate = -571.49242692 Ry
estimated scf accuracy <       2.4E-11 Ry

Total force =    0.262955    Total SCF correction =    0.000025

entering subroutine stress ...

      total stress  (Ry/bohr**3)                (kbar)    P=-4722.80
  -0.03210576 0.00000000 0.00000000    -4722.92    0.00    0.00
0.00000000  -0.03210576 0.00000000       0.00  -4722.92    0.00
0.00000000 0.00000000  -0.03210340       0.00    0.00  -4722.57
----------------------------------------------------------------------------------------------------------

回復(fù)此樓

» 猜你喜歡

0856調(diào)劑，是學(xué)校就去已經(jīng)有6人回復(fù)
288求調(diào)劑已經(jīng)有13人回復(fù)
288求調(diào)劑，一志愿華南理工大學(xué)071005 已經(jīng)有5人回復(fù)
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有5人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有12人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有13人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有11人回復(fù)
復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有11人回復(fù)
321求調(diào)劑已經(jīng)有10人回復(fù)

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

為什么采用CASTEP計(jì)算CaTiO3的光學(xué)性質(zhì)其中的靜態(tài)介電常數(shù)與實(shí)驗(yàn)值相差很大已經(jīng)有14人回復(fù)
求助結(jié)構(gòu)解析。。。。已經(jīng)有16人回復(fù)
關(guān)于碳材料的內(nèi)部結(jié)構(gòu) 已經(jīng)有5人回復(fù)
vasp計(jì)算過渡態(tài)問題已經(jīng)有13人回復(fù)
優(yōu)化和性質(zhì)計(jì)算成功，但是分析結(jié)果時(shí)鍵結(jié)構(gòu)圖看不了，怎么辦？已經(jīng)有3人回復(fù)
幾何優(yōu)化已經(jīng)有7人回復(fù)
結(jié)構(gòu)收斂與scf收斂的關(guān)系已經(jīng)有7人回復(fù)
【求助】VASP中，“結(jié)構(gòu)優(yōu)化”和“電子自洽計(jì)算”的關(guān)系？已經(jīng)有8人回復(fù)
【求助】用gaussian view打開優(yōu)化好的結(jié)構(gòu)后,將斷裂鍵連上,是否影響后面計(jì)算已經(jīng)有9人回復(fù)
【求助】castep結(jié)構(gòu)優(yōu)化與性質(zhì)計(jì)算已經(jīng)有5人回復(fù)
【求助】為什么優(yōu)化結(jié)果總是有虛頻已經(jīng)有5人回復(fù)

1樓 2010-09-21 10:50:11

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

goldenfisher

金蟲 (著名寫手)

1ST強(qiáng)帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時(shí)
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zxzj05(金幣+1):謝謝回帖交流 2010-09-21 16:52:39

是否scf和relax用的某些參數(shù)不一樣，比如k點(diǎn)，比如截?cái)嗟取?br /> 這就證明你的relax并沒有收斂，還應(yīng)該按照更精確的參數(shù)算一算

贊一下(3人)

回復(fù)此樓

2樓2010-09-21 16:35:58

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

dxcharlary

專家顧問 (職業(yè)作家)

專家經(jīng)驗(yàn): +83
1ST強(qiáng)帖: 2
應(yīng)助: 251 (大學(xué)生)
金幣: 10418.6
散金: 155
紅花: 65
帖子: 3137
在線: 757.4小時(shí)
蟲號: 936718
注冊: 2010-01-04
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 第一性原理

★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zxzj05(金幣+1):謝謝回帖交流 2010-09-22 10:47:42

一般只有在K點(diǎn)和，ecut設(shè)置不一樣的時(shí)候·會出現(xiàn)這樣的情況。仔細(xì)的對比一下，兩個(gè)輸入文件。

贊一下(2人)

回復(fù)此樓

3樓2010-09-22 10:03:40

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

qphll

金蟲 (正式寫手)

應(yīng)助: 18 (小學(xué)生)
金幣: 2708.7
散金: 2294
紅花: 8
帖子: 554
在線: 227.1小時(shí)
蟲號: 89654
注冊: 2005-08-29
性別: GG
專業(yè): 碳素材料與超硬材料

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

學(xué)習(xí)了.

回復(fù)此樓

Life, Love, Laugh.

4樓2010-09-24 09:09:20

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

gemucai

木蟲 (正式寫手)

應(yīng)助: 16 (小學(xué)生)
金幣: 1774.6
散金: 157
紅花: 4
帖子: 693
在線: 201.7小時(shí)
蟲號: 816478
注冊: 2009-07-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖
uuv2010(金幣+1): 多謝提示 2011-11-17 11:39:38

一般就算K點(diǎn)和截?cái)嗄芡耆恢乱矔霈F(xiàn)很大差別。所以每次我做完vc-relax都要重新算一遍自洽計(jì)算。

贊一下(2人)

回復(fù)此樓

5樓2011-11-17 10:53:48

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

hanyanli0475

銅蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線: 130.5小時(shí)
蟲號: 1160645
注冊: 2010-12-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

引用回帖:

1樓: Originally posted by didi5158 at 2010-09-21 10:50:11:
vc-relax計(jì)算晶體的晶格常數(shù)，用計(jì)算所得到的結(jié)構(gòu)再進(jìn)行自洽計(jì)算，而所得結(jié)果中force，stress，total energy都差別較大，這個(gè)是為什么？請大家指教。

vc-relax計(jì)算結(jié)果：
---------------------------------- ...

個(gè)人建議應(yīng)該進(jìn)行收斂測試后（ecut k都要測試），在一定收斂精度范圍內(nèi)，再去做relax和scf，并且保證兩個(gè)輸入文件是一樣的，這種情況應(yīng)該就不會發(fā)生。

贊一下(1人)

回復(fù)此樓

什么時(shí)候能算完�。。。�

6樓2011-11-17 17:03:55

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

zqclyyq

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 296.7
帖子: 40
在線: 36.3小時(shí)
蟲號: 3105748
注冊: 2014-04-01
專業(yè): 無機(jī)非金屬材料

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

引用回帖:

5樓: Originally posted by gemucai at 2011-11-17 10:53:48
一般就算K點(diǎn)和截?cái)嗄芡耆恢乱矔霈F(xiàn)很大差別。所以每次我做完vc-relax都要重新算一遍自洽計(jì)算。

vc-relax最后它自己不是會進(jìn)行一次scf嗎，如果我把vc最后得到的原子位置拿出來進(jìn)行一次scf，這兩次scf的區(qū)別在哪里，老師你知道嗎?

贊一下

回復(fù)此樓

7樓2015-12-22 11:10:15

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

siestafonda

至尊木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 20250.2
散金: 5
紅花: 6
帖子: 1486
在線: 131.5小時(shí)
蟲號: 958984
注冊: 2010-03-01
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

One possibility is that a fixed set of G-vector is selected using the initial cell based on the specified `ecutwfc`. This G-vector set is kept fixed during a vc-relax run (which can modify effectively the stored G-vectors through changes in the unit cell). As such, at the end of vc-relax, the planewave cutoff can differ from direct scf calculation if the unit cell undergoes sizable strain. One solution to this problem is to use the approach by M. Bernasconi et al, [J. Phys. Chem. Solids 56, 501 (1995), doi:10.1016/0022-3697(94)00228-2]. Relevant QE input parameters include: ecfixed, qcutz, and q2sigma see: https://www.quantum-espresso.org/Doc/INPUT_PW.html#idm406

贊一下(1人)

回復(fù)此樓

8樓2022-06-19 23:18:42

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

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

返回列表

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

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

[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +11	曦熙兮 2026-03-15	11/550	2026-03-19 17:13 by yywzz
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +8	zlingli 2026-03-13	8/400	2026-03-19 16:32 by 輕松不少隨
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 085700資源與環(huán)境308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086

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

[交流] 【求助】優(yōu)化的結(jié)果，scf計(jì)算該結(jié)構(gòu)，force為什么仍然很大？ 已有7人參與

» 猜你喜歡

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

[交流] 【求助】優(yōu)化的結(jié)果，scf計(jì)算該結(jié)構(gòu)，force為什么仍然很大？已有7人參與

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