版塊導(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í)熱門版塊排行榜

返回列表

njuxyh

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1802
散金: 12
紅花: 3
帖子: 1320
在線: 771.5小時(shí)
蟲號(hào): 455673
注冊: 2007-11-10
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】SIESTA計(jì)算總能等自洽計(jì)算時(shí)力準(zhǔn)不準(zhǔn)阿已有3人參與

siesta: Atomic forces (eV/Ang):
   1 0.002280 0.001885 0.000011
   2 -0.000823 -0.002001 0.000006
   3 -0.002785 0.000413 -0.000004
   4 -0.003841 -0.003846 0.000010
   5 -0.003414 0.002833 -0.000021
   6 0.002270 0.000616 0.000007
   7 -0.000137 0.002821 0.000011
   8 0.000639 0.000494 -0.000006
   9 -0.004456 0.004070 0.000019
上面是優(yōu)化后的受力，已經(jīng)達(dá)收斂標(biāo)準(zhǔn)了
我用這個(gè)結(jié)構(gòu)，再做scf計(jì)算
發(fā)現(xiàn)siesta: Atomic forces (eV/Ang):
1 0.002875 -0.084867 0.035258
2 -0.000382 0.039908 0.042227
3 -0.003277 0.002771 0.036465
4 -0.001384 -0.004595 0.035799
5 -0.001858 0.005840 0.035947
6 0.006556 -0.000592 0.036087
7 0.001250 -0.038206 0.042518
8 0.000911 0.087377 0.035318
9 -0.003271 -0.001110 0.034860
  10 0.002253 -0.013611 0.034873
  11 -0.001831 0.004953 0.035918
  12 0.001754 -0.005937 0.035983
可以看出，自洽計(jì)算時(shí)的輸出的力和優(yōu)化后輸出的力不一樣阿
現(xiàn)在的問題：是不是 SIESTA計(jì)算總能等自洽計(jì)算時(shí)力算不準(zhǔn)呢

回復(fù)此樓

» 猜你喜歡

材料學(xué)碩297已過四六級(jí)求調(diào)劑推薦已經(jīng)有9人回復(fù)
一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有5人回復(fù)
316求調(diào)劑已經(jīng)有4人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有8人回復(fù)
329求調(diào)劑已經(jīng)有8人回復(fù)
350求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工 322求調(diào)劑已經(jīng)有3人回復(fù)
一志愿武理材料工程348求調(diào)劑已經(jīng)有4人回復(fù)
265求調(diào)劑已經(jīng)有13人回復(fù)
290求調(diào)劑已經(jīng)有5人回復(fù)

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

求助Siesta計(jì)算熟手，想學(xué)習(xí)此軟件已經(jīng)有8人回復(fù)
siesta 計(jì)算電子占據(jù)數(shù)的問題已經(jīng)有3人回復(fù)
vasp自洽計(jì)算K點(diǎn)不匹配問題，急��！已經(jīng)有15人回復(fù)
靜態(tài)自洽計(jì)算遇到的問題已經(jīng)有5人回復(fù)
關(guān)于SIESTA中如何計(jì)算功函數(shù)的一種做法探討已經(jīng)有8人回復(fù)
【求助】siesta進(jìn)行ordern計(jì)算出錯(cuò)求助已經(jīng)有4人回復(fù)
【求助】我在用ATK做自洽計(jì)算的時(shí)候，出現(xiàn)如下的錯(cuò)誤：已經(jīng)有3人回復(fù)
【求助】VASP中，“結(jié)構(gòu)優(yōu)化”和“電子自洽計(jì)算”的關(guān)系？已經(jīng)有8人回復(fù)
【求助】結(jié)構(gòu)優(yōu)化之后的靜態(tài)自洽計(jì)算已經(jīng)有12人回復(fù)
【討論】自洽計(jì)算與非自洽計(jì)算的含義、所做的工作已經(jīng)有10人回復(fù)
【原創(chuàng)】并行自洽計(jì)算后得到和節(jié)點(diǎn)數(shù)一樣多的波函數(shù)文件wfc 已經(jīng)有4人回復(fù)
【求助】靜態(tài)自洽計(jì)算讀WAVECAR有用嗎？已經(jīng)有9人回復(fù)
【求助】請(qǐng)教PWSCF晶體建模和自洽計(jì)算，不收斂已經(jīng)有21人回復(fù)
【求助】態(tài)密度計(jì)算究竟是自洽計(jì)算，還是非自洽計(jì)算？已經(jīng)有24人回復(fù)

1樓 2010-11-05 15:13:37

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

zhangguangping

木蟲 (著名寫手)

1ST強(qiáng)帖: 23
應(yīng)助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時(shí)
蟲號(hào): 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

★
zxzj05(金幣+1):感謝交流！ 2010-11-05 21:04:27

引用回帖:

Originally posted by njuxyh at 2010-11-05 08:13:37:
siesta: Atomic forces (eV/Ang):
   1 0.002280 0.001885 0.000011
   2 -0.000823 -0.002001 0.000006
   3 -0.002785 0.000413 -0.000004
   4 -0.003841 -0.003846 ...

如果你新的計(jì)算SCF自洽結(jié)束的話，應(yīng)該差別不是很大。

贊一下(1人)

回復(fù)此樓

弘德明志博學(xué)篤行

2樓2010-11-05 17:18:21

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

njuxyh

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1802
散金: 12
紅花: 3
帖子: 1320
在線: 771.5小時(shí)
蟲號(hào): 455673
注冊: 2007-11-10
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

這就是自洽最后的結(jié)果：
1 0.002875 -0.084867 0.035258
2 -0.000382 0.039908 0.042227
3 -0.003277 0.002771 0.036465
4 -0.001384 -0.004595 0.035799
5 -0.001858 0.005840 0.035947
6 0.006556 -0.000592 0.036087
7 0.001250 -0.038206 0.042518
8 0.000911 0.087377 0.035318
9 -0.003271 -0.001110 0.034860
  10 0.002253 -0.013611 0.034873
  11 -0.001831 0.004953 0.035918
  12 0.001754 -0.005937 0.035983

本文來自: 小木蟲論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=2565072&page=1#pid1801284

[ Last edited by njuxyh on 2010-11-5 at 17:56 ]

贊一下

回復(fù)此樓

3樓2010-11-05 17:55:28

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

zhangguangping

木蟲 (著名寫手)

1ST強(qiáng)帖: 23
應(yīng)助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時(shí)
蟲號(hào): 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

★
zxzj05(金幣+1):感謝交流！ 2010-11-05 21:04:44
njuxyh(金幣+5): 2010-11-08 10:12:43

引用回帖:

Originally posted by njuxyh at 2010-11-05 10:55:28:
這就是自洽最后的結(jié)果：
1 0.002875 -0.084867 0.035258
2 -0.000382 0.039908 0.042227
3 -0.003277 0.002771 0.036465
4 -0.001384 -0.004595 0.035799
5 - ...

好像這個(gè)就是傳說中的不同的初始DM猜測得到不同的結(jié)果吧。我在siesta-3.0-b和siesta-3.0-rc2之間遇到過這個(gè)問題。

贊一下(1人)

回復(fù)此樓

弘德明志博學(xué)篤行

4樓2010-11-05 19:46:34

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

njuxyh

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1802
散金: 12
紅花: 3
帖子: 1320
在線: 771.5小時(shí)
蟲號(hào): 455673
注冊: 2007-11-10
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

你是說版本不同造成的?

贊一下

回復(fù)此樓

5樓2010-11-05 20:19:37

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

zhangguangping

木蟲 (著名寫手)

1ST強(qiáng)帖: 23
應(yīng)助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時(shí)
蟲號(hào): 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

★
zxzj05(金幣+1):感謝交流！ 2010-11-05 21:05:04

引用回帖:

Originally posted by njuxyh at 2010-11-05 13:19:37:
你是說版本不同造成的?

我還沒遇到同一個(gè)版本中你的情況，但是我遇到兩個(gè)版本中是有這個(gè)問題的。但是他們同一個(gè)初始條件能得到一致的結(jié)果，但是讓他們自己猜測，結(jié)果就不同了。

贊一下(1人)

回復(fù)此樓

弘德明志博學(xué)篤行

6樓2010-11-05 20:31:52

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

njuxyh

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1802
散金: 12
紅花: 3
帖子: 1320
在線: 771.5小時(shí)
蟲號(hào): 455673
注冊: 2007-11-10
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

謝謝，不知道這些差別會(huì)不會(huì)影響我得總能計(jì)算？因?yàn)楦杏X后者的原子受力比較大，好想不是最穩(wěn)定的結(jié)構(gòu)？

贊一下

回復(fù)此樓

7樓2010-11-06 10:06:43

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

zhangguangping

木蟲 (著名寫手)

1ST強(qiáng)帖: 23
應(yīng)助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時(shí)
蟲號(hào): 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
cenwanglai(金幣+2):謝謝回復(fù)！ 2010-11-08 12:48:01

引用回帖:

Originally posted by njuxyh at 2010-11-06 03:06:43:
謝謝，不知道這些差別會(huì)不會(huì)影響我得總能計(jì)算？因?yàn)楦杏X后者的原子受力比較大，好想不是最穩(wěn)定的結(jié)構(gòu)？

即使是同一個(gè)版本，在用siesta的時(shí)候如果收斂了，再次運(yùn)算的話，如果用上一個(gè)DM，會(huì)在1到2個(gè)SCF之內(nèi)自洽。但是如果是做輸運(yùn)用transiesta的話，即使是收斂過的DM，如果再次用，還會(huì)發(fā)現(xiàn)不收斂，并且經(jīng)過的SCF比較多。這個(gè)確實(shí)是一個(gè)問題。但是你的問題我在使用siesta的時(shí)候沒有遇到過。尤其是你的幾何結(jié)構(gòu)是按照絕對(duì)的收斂條件來的。按照這個(gè)結(jié)構(gòu)再去算一次，即使是從頭猜測DM，只要收斂了，當(dāng)然DM的收斂標(biāo)準(zhǔn)一定要和優(yōu)化的時(shí)候用的一樣的收斂標(biāo)準(zhǔn)。那么應(yīng)該收斂后的受力應(yīng)該不會(huì)差別很大，最起碼還是滿足幾何結(jié)構(gòu)優(yōu)化的時(shí)候的那個(gè)條件的。你的這個(gè)問題可以寫信給siesta-l@uam.es，這個(gè)是siesta的一個(gè)討論區(qū)，但是需要注冊才能給這個(gè)地址發(fā)郵件。有可能你的問題，在這個(gè)能找到答案。

贊一下(3人)

回復(fù)此樓

弘德明志博學(xué)篤行

8樓2010-11-08 10:23:43

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

njuxyh

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 1802
散金: 12
紅花: 3
帖子: 1320
在線: 771.5小時(shí)
蟲號(hào): 455673
注冊: 2007-11-10
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

謝謝，學(xué)習(xí)了

回復(fù)此樓

9樓2010-11-09 10:20:16

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

295399157

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 86
帖子: 3
在線: 7.2小時(shí)
蟲號(hào): 1138517
注冊: 2010-11-04
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

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

我用atk來算一個(gè)碳管的透射譜，算了幾天后出現(xiàn)一個(gè)這樣的錯(cuò)誤命令提示“Exceed maximum self-consistent ”。請(qǐng)問哪位高人知道這該怎么辦？

贊一下(1人)

回復(fù)此樓

10樓2010-11-27 22:38:25

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

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

返回列表

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

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

[考研] 一志愿武理材料工程348求調(diào)劑 +3	￣^￣゜汗 2026-03-19	4/200	2026-03-20 21:01 by zhukairuo
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	13/650	2026-03-20 21:01 by 無際的草原
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4	一瓶苯甲酸 2026-03-14	4/200	2026-03-20 20:36 by fen_rao
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3	123阿標(biāo) 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 304求調(diào)劑 +5	曼殊2266 2026-03-18	5/250	2026-03-20 09:00 by ZHANG0tao
[考研] 307求調(diào)劑 +9	冷笙123 2026-03-17	9/450	2026-03-19 22:44 by 學(xué)員8dgXkO
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	6/300	2026-03-19 17:18 by fei626-918
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3	@taotao 2026-03-19	3/150	2026-03-19 14:07 by peike
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi

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

[交流] 【求助】SIESTA計(jì)算總能等自洽計(jì)算時(shí)力準(zhǔn)不準(zhǔn)阿 已有3人參與

» 猜你喜歡

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

[交流] 【求助】SIESTA計(jì)算總能等自洽計(jì)算時(shí)力準(zhǔn)不準(zhǔn)阿已有3人參與