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

xiongxiong5712

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 543
散金: 12
紅花: 2
帖子: 257
在線: 267.7小時(shí)
蟲號(hào): 706467
注冊(cè): 2009-02-23
專業(yè): 半導(dǎo)體材料

[求助] 原子受力沒有達(dá)到收斂標(biāo)準(zhǔn)計(jì)算就結(jié)束了已有1人參與

用siesta進(jìn)行結(jié)構(gòu)優(yōu)化，由于是六方晶系，首先通過(guò)擬合得到了能量最低對(duì)應(yīng)的體積，MD.TypeOfRun設(shè)置為CG，然后再保持這個(gè)體積不變，優(yōu)化晶格參數(shù)，MD.TypeOfRun 設(shè)置為Broyden，具體的fdf文件如下：
SystemName       bulk
SystemLabel       bulk

NumberOfAtoms    14
NumberOfSpecies    2

%block ChemicalSpeciesLabel
1 14  Si
2 7 N
%endblock ChemicalSpeciesLabel

LatticeConstant    1 Ang

%block LatticeVectors
   6.571864087507342    -3.794267500000001    0.000000000000000
   0.000000000000000    7.588534999999999    0.000000000000000
   0.000000000000000    0.000000000000000    2.906249000000000
%endblock LatticeVectors

%block kgrid_Monkhorst_Pack
5 0 0 0
0 5 0 0
0 0 12 0
%endblock kgrid_Monkhorst_Pack

AtomicCoordinatesFormat    Fractional
%block AtomicCoordinatesAndAtomicSpecies
0.3323000000000000 0.0314000000000000 0.2500000000000000 2
-0.0314000000000000 0.3009000000000000 0.2500000000000000 2
-0.3009000000000000 -0.3323000000000000 0.2500000000000000 2
-0.3323000000000000 -0.0314000000000000 0.7500000000000000 2
0.0314000000000000 -0.3009000000000000 0.7500000000000000 2
0.3009000000000000 0.3323000000000000 0.7500000000000000 2
0.3333333333333330 0.6666666666666670 0.2500000000000000 2
-0.3333333333333330 -0.6666666666666670 0.7500000000000000 2
0.1733000000000000 -0.2306000000000000 0.2500000000000000 1
0.2306000000000000 0.4039000000000000 0.2500000000000000 1
-0.4039000000000000 -0.1733000000000000 0.2500000000000000 1
-0.1733000000000000 0.2306000000000000 0.7500000000000000 1
-0.2306000000000000 -0.4039000000000000 0.7500000000000000 1
0.4039000000000000 0.1733000000000000 0.7500000000000000 1
%endblock AtomicCoordinatesAndAtomicSpecies

xc.functional       LDA
xc.authors          CA
MeshCutoff          100 Ry

PAO.BasisSize       DZP
PAO.EnergyShift    75 mev

MD.MaxForceTol    0.05 eV/Ang

DM.MixingWeight    0.01
DM.NumberPulay    10
MaxSCFIterations    200

SolutionMethod    diagon

MD.TypeOfRun       Broyden
MD.NumCGsteps       500

MD.VariableCell       .true.
MD.ConstantVolume    .true.
MD.RelaxCellOnly    .true.

優(yōu)化很快， Begin CG move =  1的時(shí)候就結(jié)束了，dMax達(dá)到了收斂標(biāo)準(zhǔn)，但是原子受力沒有收斂，結(jié)果是：
siesta: Atomic forces (eV/Ang):
----------------------------------------
Tot -0.000014 -0.000041 -0.000144
----------------------------------------
Max 0.722244
Res 0.314921 sqrt( Sum f_i^2 / 3N )
----------------------------------------
Max 0.722244 constrained

可見原子受力沒有達(dá)到0.05以下，而計(jì)算就正常結(jié)束了。請(qǐng)教大家，這是什么原因引起的啊？是不是我MD.TypeOfRun 設(shè)置為Broyden的原因呢？用這種方法優(yōu)化的時(shí)候應(yīng)該注意什么？有勞各位幫幫忙，多謝了。

[ Last edited by xiongxiong5712 on 2013-5-31 at 21:46 ]

回復(fù)此樓

» 猜你喜歡

323求調(diào)劑已經(jīng)有6人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有4人回復(fù)
352求調(diào)劑已經(jīng)有3人回復(fù)
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)

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

燃燒計(jì)算，先收斂后發(fā)散已經(jīng)有12人回復(fù)
關(guān)于回流與計(jì)算結(jié)果收斂的討論已經(jīng)有11人回復(fù)
ansys接觸計(jì)算不收斂問(wèn)題已經(jīng)有16人回復(fù)
pwscf計(jì)算integrated charge出錯(cuò)及不收斂已經(jīng)有13人回復(fù)
怎樣判定fluent計(jì)算是否收斂已經(jīng)有4人回復(fù)
NEB計(jì)算擴(kuò)散路徑的力收斂精度設(shè)置問(wèn)題已經(jīng)有16人回復(fù)
Castep計(jì)算無(wú)法收斂，請(qǐng)高人指教已經(jīng)有17人回復(fù)
加溶劑后計(jì)算，不收斂，請(qǐng)高手指點(diǎn) 已經(jīng)有5人回復(fù)
CASTEP計(jì)算摻雜不收斂已經(jīng)有5人回復(fù)
計(jì)算時(shí)候總是不能收斂如何解決？已經(jīng)有5人回復(fù)
靜態(tài)計(jì)算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？已經(jīng)有4人回復(fù)
ELK計(jì)算嚴(yán)重不收斂～已經(jīng)有10人回復(fù)
vasp計(jì)算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計(jì)算不收斂已經(jīng)有3人回復(fù)
IRC計(jì)算沒收斂已經(jīng)有9人回復(fù)
DMol3計(jì)算頻率不收斂已經(jīng)有4人回復(fù)
【求助】用Orca計(jì)算一個(gè)片簇的單點(diǎn)能，不收斂。已經(jīng)有7人回復(fù)
【求助】計(jì)算不收斂怎么辦？已經(jīng)有4人回復(fù)
【求助】請(qǐng)教PWSCF晶體建模和自洽計(jì)算，不收斂已經(jīng)有21人回復(fù)

1樓 2013-05-31 16:23:13

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

newton3915

鐵桿木蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 4742.3
散金: 40
紅花: 4
帖子: 418
在線: 391.3小時(shí)
蟲號(hào): 609189
注冊(cè): 2008-09-22
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by xiongxiong5712 at 2013-06-01 15:10:48
自己頂一下。麻煩各位大俠百忙中幫我看看，或者提示一下也行。我先說(shuō)說(shuō)我的看法：是不是因?yàn)镸D.TypeOfRun 不是CG，計(jì)算不會(huì)進(jìn)行CGmove的循環(huán)，所以在Begin CG move = 1的時(shí)候就結(jié)束了，是否收斂也不是看force里的M ...

LZ問(wèn)題解決了嗎，我也是遇到同樣的問(wèn)題？

贊一下

回復(fù)此樓

3樓2014-12-15 04:22:57

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

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

xiongxiong5712

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 543
散金: 12
紅花: 2
帖子: 257
在線: 267.7小時(shí)
蟲號(hào): 706467
注冊(cè): 2009-02-23
專業(yè): 半導(dǎo)體材料

自己頂一下。麻煩各位大俠百忙中幫我看看，或者提示一下也行。我先說(shuō)說(shuō)我的看法：是不是因?yàn)镸D.TypeOfRun 不是CG，計(jì)算不會(huì)進(jìn)行CGmove的循環(huán)，所以在Begin CG move =  1的時(shí)候就結(jié)束了，是否收斂也不是看force里的Max，而是有一個(gè)地方給出了每個(gè)原子的受力：
siesta: Atomic forces (eV/Ang):
siesta:    1 -0.368915 -0.160954 -0.000062
siesta:    2 0.323778 -0.238986 -0.000043
siesta:    3 0.045044 0.399972 0.000016
siesta:    4 0.368938 0.160918 0.000009
siesta:    5 -0.323802 0.238974 -0.000020
siesta:    6 -0.045044 -0.399962 -0.000040
siesta:    7 0.000000 -0.000020 -0.000041
siesta:    8 -0.000020 0.000009 0.000023
siesta:    9 0.450554 -0.573922 -0.000001
siesta:    10 0.271759 0.677123 0.000003
siesta:    11 -0.722307 -0.103215 -0.000009
siesta:    12 -0.450542 0.573917 0.000002
siesta:    13 -0.271775 -0.677130 0.000005
siesta:    14 0.722304 0.103217 0.000011
siesta: ----------------------------------------
siesta: Tot -0.000027 -0.000058 -0.000148
最后一行是總的受力，是不是看這個(gè)地方收斂即可？（CG方法沒有給出這一結(jié)果）。不知道我猜測(cè)的對(duì)不對(duì)，請(qǐng)大家指正。

贊一下

回復(fù)此樓

2樓2013-06-01 15:10:48

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

chengxiajie

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 43.5
散金: 10
帖子: 7
在線: 7.3小時(shí)
蟲號(hào): 4127351
注冊(cè): 2015-10-08

【答案】應(yīng)助回帖

你的問(wèn)題解決了嘛？你看一下是不是你的MD.NumCGsteps設(shè)置的太小了

贊一下

回復(fù)此樓

4樓2016-05-10 22:10:45

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

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

[考研] 277材料科學(xué)與工程080500求調(diào)劑 +7	自由煎餅果子 2026-03-16	7/350	2026-03-22 22:40 by ACS Nano——
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	4/200	2026-03-22 15:12 by 朱芷琳
[考研] 生物學(xué)調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +3	酥酥魚.. 2026-03-21	3/150	2026-03-22 11:31 by 楊楊楊紫
[考研] 354求調(diào)劑 +7	Tyoumou 2026-03-18	10/500	2026-03-22 11:11 by 人來(lái)盛
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 求助 +5	夢(mèng)里的無(wú)言 2026-03-21	6/300	2026-03-21 17:51 by 學(xué)員8dgXkO
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 22408 344分求調(diào)劑一志愿華電計(jì)算機(jī)技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶 2026-03-17	9/450	2026-03-20 14:24 by 無(wú)懈可擊111
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao