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

c_chi

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 17 (小學(xué)生)
金幣: 2449.7
散金: 228
紅花: 2
帖子: 551
在線: 493.6小時(shí)
蟲(chóng)號(hào): 1253512
注冊(cè): 2011-04-02
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] siesta基組優(yōu)化問(wèn)題已有1人參與

和PW方法相比，LCAO方法有時(shí)候的計(jì)算準(zhǔn)確度確實(shí)有所不夠，特別是對(duì)于過(guò)渡金屬和貴金屬感覺(jué)不太容易計(jì)算，之前用siesta優(yōu)化一個(gè)Beta相的MnO2，發(fā)現(xiàn)結(jié)果基本是錯(cuò)的。問(wèn)了些人，也在論壇里翻了其他人的帖子，說(shuō)是可以通過(guò)優(yōu)化基組來(lái)達(dá)到PW的準(zhǔn)確度（工具：/Util/Optimizer目錄下的simplex），但是還是不怎么明白。
http://www.gaoyang168.com/bbs/viewthread.php?tid=2486484&fpage=1&target=blank
http://www.gaoyang168.com/bbs/viewthread.php?tid=2485045&fpage=1&target=self&page=1
http://www.gaoyang168.com/bbs/viewthread.php?tid=3146982
也看了高手提供的一篇這方面的文章，看了下，沒(méi)能全看懂。
Optimal strictly localized basis sets for noble metal surfaces.（PHYSICAL REVIEW B79, 075441 2009）

這里希望有人能夠解一下疑惑：
1詳細(xì).解釋下基組，大家經(jīng)驗(yàn)來(lái)看什么情況下使用什么基組可以達(dá)到可行的精度？
2.解釋下siesta基組優(yōu)化的具體方法和技巧，或者提供這方面的資料？或者可不可以參考高斯的基組來(lái)設(shè)置siesta的基組？
3.我在siesta計(jì)算過(guò)渡金屬氧化物時(shí)，感覺(jué)scf比較難以收斂，如果優(yōu)化基組后可能更難，計(jì)算量更大，我一般的做法是加大DM.NumberPulay并減小DM.MixingWeight，如下。請(qǐng)問(wèn)我的做法合理嗎？還有什么其他方法？
DM.NumberPulay 8
DM.MixingWeight 0.1
4. 最后對(duì)于贗勢(shì)一般是用軟件提供的贗勢(shì)，這里也能夠優(yōu)化來(lái)提高準(zhǔn)確度？需要怎么做？
先謝謝了！��！

回復(fù)此樓

» 猜你喜歡

一志愿大連理工大學(xué)材料求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工調(diào)劑一志愿大連海事085600，349 已經(jīng)有6人回復(fù)
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程調(diào)劑已經(jīng)有11人回復(fù)
311求調(diào)劑一志愿合肥工業(yè)大學(xué) 已經(jīng)有4人回復(fù)
一志愿食品科學(xué)與工程083200求調(diào)劑已經(jīng)有4人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有10人回復(fù)
一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利｜求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩已經(jīng)有4人回復(fù)
085600材料與化工調(diào)劑已經(jīng)有15人回復(fù)

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

基組or結(jié)構(gòu)問(wèn)題，B3lyp優(yōu)化成功，B3PW91/WB97XD優(yōu)化全出錯(cuò) 已經(jīng)有16人回復(fù)
請(qǐng)教如何加快siesta的收斂速度已經(jīng)有9人回復(fù)
Siesta從官網(wǎng)下載贗勢(shì)計(jì)算出錯(cuò) 已經(jīng)有7人回復(fù)
求助siesta收斂問(wèn)題已經(jīng)有3人回復(fù)
贗試基組輸入出錯(cuò)，求教已經(jīng)有5人回復(fù)
siesta結(jié)構(gòu)優(yōu)化收斂標(biāo)準(zhǔn) 已經(jīng)有5人回復(fù)
SIESTA計(jì)算能級(jí)不準(zhǔn)怎么解決已經(jīng)有10人回復(fù)
求助siesta優(yōu)化問(wèn)題已經(jīng)有3人回復(fù)
siesta 計(jì)算電子占據(jù)數(shù)的問(wèn)題已經(jīng)有3人回復(fù)
求助贗勢(shì)基組的問(wèn)題已經(jīng)有9人回復(fù)
關(guān)于SIESTA中如何計(jì)算功函數(shù)的一種做法探討已經(jīng)有8人回復(fù)
手動(dòng)寫(xiě)入原子軌道基組的原則和技巧已經(jīng)有5人回復(fù)
SIESTA幾何優(yōu)化不收斂是不是就說(shuō)明了這個(gè)結(jié)構(gòu)是不穩(wěn)定的? 已經(jīng)有4人回復(fù)
siesta計(jì)算弱相互作用VDW 已經(jīng)有17人回復(fù)
【求助】siesta 優(yōu)化晶胞已經(jīng)有3人回復(fù)
【分享】SIESTA中Bader analysis 已經(jīng)有43人回復(fù)
【求助】siesta 計(jì)算DOS應(yīng)該注意什么問(wèn)題已經(jīng)有3人回復(fù)
【討論】關(guān)于如何優(yōu)化大體系的問(wèn)題已經(jīng)有9人回復(fù)
【求助】siesta多個(gè)/種分子優(yōu)化的問(wèn)題已經(jīng)有7人回復(fù)
【求助】siesta輸入文件中贗勢(shì)文件的作用？已經(jīng)有8人回復(fù)
【求助】誰(shuí)能具體說(shuō)一下各個(gè)基組的優(yōu)缺點(diǎn) 已經(jīng)有5人回復(fù)
【求助】siesta中關(guān)于幾何結(jié)構(gòu)優(yōu)化的問(wèn)題【完結(jié)】已經(jīng)有17人回復(fù)
【求助】求gamess 輸入使用贗勢(shì)基組與分裂基組的例子已經(jīng)有6人回復(fù)

1樓 2014-09-09 14:35:20

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

c_chi

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 17 (小學(xué)生)
金幣: 2449.7
散金: 228
紅花: 2
帖子: 551
在線: 493.6小時(shí)
蟲(chóng)號(hào): 1253512
注冊(cè): 2011-04-02
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

2樓: Originally posted by liqizuiyang at 2014-09-09 16:01:22
推薦樓主去看下openmx的幫助文檔和項(xiàng)目組發(fā)的文章，我記得有篇PRB文章專(zhuān)門(mén)講基組優(yōu)化的。雖然方法不同，但都是LCAO，應(yīng)該有相通的部分。

LCAO使用原子軌道作為基組，原子軌道可分為徑向函數(shù)和球諧函數(shù)兩部分。球 ...

謝謝，我去研究下opemx的文檔和文章。

贊一下

回復(fù)此樓

3樓2014-09-10 08:30:24

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

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

liqizuiyang

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

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

★ ★ ★
franch: 金幣+3, 謝謝回帖交流。。。 2014-09-22 22:29:11

推薦樓主去看下openmx的幫助文檔和項(xiàng)目組發(fā)的文章，我記得有篇PRB文章專(zhuān)門(mén)講基組優(yōu)化的。雖然方法不同，但都是LCAO，應(yīng)該有相通的部分。

LCAO使用原子軌道作為基組，原子軌道可分為徑向函數(shù)和球諧函數(shù)兩部分。球諧函數(shù)有確定的表達(dá)式，徑向函數(shù)有slater型，gauss型和數(shù)值型。數(shù)值型就是把徑向函數(shù)分成一系列格點(diǎn)，存儲(chǔ)格點(diǎn)上的函數(shù)值，計(jì)算時(shí)再用插值法還原。基組優(yōu)化實(shí)際上就是調(diào)整格點(diǎn)上的值。

贊一下

回復(fù)此樓

2樓2014-09-09 16:01:22

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

c_chi

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 17 (小學(xué)生)
金幣: 2449.7
散金: 228
紅花: 2
帖子: 551
在線: 493.6小時(shí)
蟲(chóng)號(hào): 1253512
注冊(cè): 2011-04-02
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

問(wèn)題還是沒(méi)解決，我主要是想優(yōu)化MnO2的基組，我看了下例子不是很明白，優(yōu)化參數(shù)和原子軌道有什么關(guān)系嗎？能解釋下下面的內(nèi)容：
1.有個(gè)如下的VASRS文件，應(yīng)該是優(yōu)化的閾值，但是怎么�。坑惺裁匆�(guī)則還不是很明白？為什么要這樣取這些軌道？
s_z1_rc_O 2.5  8.0
spl_s_O 0.05  0.5  0.15
p_z1_rc_O 2.5  8.0
spl_p_O 0.05  0.5
d_z1_rc_O 2.5  8.0
s_z1_rc_H 2.0  8.0
spl_s_H 0.05  0.8
p_z1_rc_H 2.0  8.0
2.同時(shí)運(yùn)行的輸入文件有個(gè)基組的block，這個(gè)和上面的優(yōu)化閾值有什么關(guān)系嗎？謝謝了�。�！
%block PAO.Basis                                  # Define Basis set
O    3
n=2 0 2  S $spl_s_O
$s_z1_rc_O 0.0
1.000 1.000
n=2 1 2  S $spl_p_O
$p_z1_rc_O 0.0
1.000 1.000
n=3 2 1
$d_z1_rc_O
1.000
H    2
n=1 0 2  S $spl_s_H
$s_z1_rc_H  0.0
1.000 1.000
n=2 1 1
$p_z1_rc_H
1.000
%endblock PAO.Basis

贊一下

回復(fù)此樓

4樓2014-09-11 14:45:11

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

c_chi

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 17 (小學(xué)生)
金幣: 2449.7
散金: 228
紅花: 2
帖子: 551
在線: 493.6小時(shí)
蟲(chóng)號(hào): 1253512
注冊(cè): 2011-04-02
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

沒(méi)有蟲(chóng)子對(duì)這個(gè)有研究�。�

贊一下

回復(fù)此樓

5樓2014-09-17 14:04:48

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

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

[考研] 288資源與環(huán)境專(zhuān)碩求調(diào)劑，不限專(zhuān)業(yè)，有學(xué)上就行 +11	lllllos 2026-03-30	11/550	2026-03-30 21:19 by 研究僧導(dǎo)導(dǎo)
[考研] 359求調(diào)劑 +5	王了個(gè)楠 2026-03-25	5/250	2026-03-30 19:36 by 源_2020
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +5	12.27 2026-03-27	5/250	2026-03-30 19:31 by 源_2020
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[考研] 322求調(diào)劑 +10	宋明欣 2026-03-27	10/500	2026-03-30 18:47 by 544594351
[考研] 298求調(diào)劑 +3	什么是胖頭魚(yú) 2026-03-30	5/250	2026-03-30 14:41 by 青海小西牛
[碩博家園] 求調(diào)劑有機(jī)化學(xué)考研356分 +10	Nadiums 2026-03-25	11/550	2026-03-30 11:56 by yjolah
[考研] 070300一志愿211，312分求調(diào)劑院校 +4	小黃鴨寶 2026-03-30	4/200	2026-03-30 11:49 by 棒棒球手
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8	柳溪邊 2026-03-26	8/400	2026-03-29 20:42 by 唐沐兒
[考研] 298求調(diào)劑 +3	種圣賜 2026-03-29	3/150	2026-03-29 12:06 by longlotian
[考研] 085600，專(zhuān)業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 320分，材料與化工專(zhuān)業(yè)，求調(diào)劑 +9	一定上岸aaa 2026-03-27	13/650	2026-03-28 15:00 by 神馬都不懂
[考研] 一志愿廈門(mén)大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +10	y7czhao 2026-03-26	10/500	2026-03-28 14:23 by 唐沐兒
[考研] 286求調(diào)劑 +12	PolarBear11 2026-03-26	12/600	2026-03-28 12:14 by zllcz
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 一志愿211院校 344分東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5	丶風(fēng)雪夜歸人丶 2026-03-26	8/400	2026-03-27 19:22 by 丶風(fēng)雪夜歸人丶
[考研] 考研調(diào)劑 +10	呼呼？~+123456 2026-03-24	10/500	2026-03-27 11:46 by wangjy2002
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu