 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 2078 | 回復(fù): 7 | ||
xwu10銅蟲 (小有名氣)
|
[求助]
SIESTA 計算石墨烯，模型增大體系不收斂 已有1人參與
|
|
各位好 我在用SIESTA進行石墨烯的簡單計算時，遇到了一個問題。我用MS進行石墨烯體系的建模，當(dāng)石墨烯體系較小時（小于30個原子），SIESTA計算可以很好地收斂；然而當(dāng)我增大石墨烯體系時（其他參數(shù)不變），去發(fā)現(xiàn)體系始終無法收斂，這個時候不管我修改何種參數(shù)，體系都無法收斂。我想問的是，SIESTA針對大體系計算難道就會出現(xiàn)收斂性問題嗎？這是軟件本身的bug還是什么其他原因？各位是否有可靠建議？（我試過修改MixingWeight和NumberPulay參數(shù)，以及Kgrid，結(jié)果依舊不能收斂） 兩種體系的fdf文件如下，也供初學(xué)者學(xué)習(xí)參考 #graphene-small supercell SystemName graphene SystemLabel graphene NumberOfAtoms 30 NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 6 C %endblock ChemicalSpeciesLabel %block PAO.Basis # Define Basis set C 2 # Species label, number of l-shells n=2 0 1 # n, l, Nzeta 4.088 1.000 n=2 1 1 # n, l, Nzeta, Polarization, NzetaPol 4.870 1.000 %endblock PAO.Basis XC.functional GGA LatticeConstant 1.0 Ang %block LatticeVectors 20.00000000 0.0000000000 0.0000000000 0.000000000 15.0000000000 0.0000000000 0.000000000 0.0000000000 20.0000000000 %endblock LatticeVectors AtomicCoordinatesFormat Ang %block AtomicCoordinatesAndAtomicSpecies 1.7 2.13 1.23 1 1.7 2.841 0 1 1.7 4.261 0 1 1.7 4.971 1.23 1 1.7 6.391 1.23 1 1.7 7.101 0 1 1.7 8.522 0 1 1.7 9.232 1.23 1 1.7 10.652 1.23 1 1.7 11.362 0 1 1.7 2.13 3.69 1 1.7 2.841 2.46 1 1.7 4.261 2.46 1 1.7 4.971 3.69 1 1.7 6.391 3.69 1 1.7 7.101 2.46 1 1.7 8.522 2.46 1 1.7 9.232 3.69 1 1.7 10.652 3.69 1 1.7 11.362 2.46 1 1.7 2.13 6.15 1 1.7 2.841 4.92 1 1.7 4.261 4.92 1 1.7 4.971 6.15 1 1.7 6.391 6.15 1 1.7 7.101 4.92 1 1.7 8.522 4.92 1 1.7 9.232 6.15 1 1.7 10.652 6.15 1 1.7 11.362 4.92 1 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 10 0 0.0 0 0 10 0.5 %endblock Kgrid_Monkhorst_Pack MeshCutoff 500.0 Ry MaxSCFIterations 1000 DM.MixingWeight 0.05 DM.NumberPulay 10 DM.Tolerance 0.0001 WriteCoorXmol .true. BandLinesScale ReciprocalLatticeVectors %block BandLines 1 0.0 0.0 0.0 \Gamma 100 0.0 0.0 0.5 X %endblock BandLines SolutionMethod diagon TS.WriteHS .true. #graphene-large supercell SystemName graphene SystemLabel graphene NumberOfAtoms 50 NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 6 C %endblock ChemicalSpeciesLabel %block PAO.Basis # Define Basis set C 2 # Species label, number of l-shells n=2 0 1 # n, l, Nzeta 4.088 1.000 n=2 1 1 # n, l, Nzeta, Polarization, NzetaPol 4.870 1.000 %endblock PAO.Basis XC.functional GGA XC.authors PBE LatticeConstant 1.0 Ang %block LatticeVectors 20.00000000 0.0000000000 0.0000000000 0.000000000 15.0000000000 0.0000000000 0.000000000 0.0000000000 20 %endblock LatticeVectors AtomicCoordinatesFormat Ang %block AtomicCoordinatesAndAtomicSpecies 1.7 2.13 1.23 1 1.7 2.841 0 1 1.7 4.261 0 1 1.7 4.971 1.23 1 1.7 6.391 1.23 1 1.7 7.101 0 1 1.7 8.522 0 1 1.7 9.232 1.23 1 1.7 10.652 1.23 1 1.7 11.362 0 1 1.7 2.13 3.69 1 1.7 2.841 2.46 1 1.7 4.261 2.46 1 1.7 4.971 3.69 1 1.7 6.391 3.69 1 1.7 7.101 2.46 1 1.7 8.522 2.46 1 1.7 9.232 3.69 1 1.7 10.652 3.69 1 1.7 11.362 2.46 1 1.7 2.13 6.15 1 1.7 2.841 4.92 1 1.7 4.261 4.92 1 1.7 4.971 6.15 1 1.7 6.391 6.15 1 1.7 7.101 4.92 1 1.7 8.522 4.92 1 1.7 9.232 6.15 1 1.7 10.652 6.15 1 1.7 11.362 4.92 1 1.7 2.13 8.61 1 1.7 2.841 7.38 1 1.7 4.261 7.38 1 1.7 4.971 8.61 1 1.7 6.391 8.61 1 1.7 7.101 7.38 1 1.7 8.522 7.38 1 1.7 9.232 8.61 1 1.7 10.652 8.61 1 1.7 11.362 7.38 1 1.7 2.13 11.07 1 1.7 2.841 9.84 1 1.7 4.261 9.84 1 1.7 4.971 11.07 1 1.7 6.391 11.07 1 1.7 7.101 9.84 1 1.7 8.522 9.84 1 1.7 9.232 11.07 1 1.7 10.652 11.07 1 1.7 11.362 9.84 1 %endblock AtomicCoordinatesAndAtomicSpecies %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 10 0 0.0 0 0 10 0.5 %endblock Kgrid_Monkhorst_Pack MeshCutoff 500.0 Ry MaxSCFIterations 1000 DM.MixingWeight 0.1 DM.NumberPulay 10 DM.Tolerance 0.0001 WriteCoorXmol .true. BandLinesScale ReciprocalLatticeVectors %block BandLines 1 0.0 0.0 0.0 \Gamma 100 0.0 0.0 0.5 X %endblock BandLines SolutionMethod diagon TS.WriteHS .true. |
回復(fù)此樓» 猜你喜歡
一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑
已經(jīng)有3人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料學(xué)碩297已過四六級求調(diào)劑推薦
已經(jīng)有4人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有14人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有8人回復(fù)
-
材料080500調(diào)劑求收留
已經(jīng)有5人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有7人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
085410人工智能專碩317求調(diào)劑（0854都可以）
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 靜態(tài)計算，加U計算不收斂！當(dāng)U>=7時，就不收斂了，當(dāng)U小于7時很快收斂！
已經(jīng)有8人回復(fù)
- 計算機程序設(shè)計藝術(shù)（第三版）1~3卷
已經(jīng)有46人回復(fù)
- 石墨烯轉(zhuǎn)移膜層不連續(xù)
已經(jīng)有6人回復(fù)
- 請教如何加快siesta的收斂速度
已經(jīng)有9人回復(fù)
- 雙層石墨烯的間距多大啊 優(yōu)化后 謝謝
已經(jīng)有7人回復(fù)
- X-crysden安裝各種問題
已經(jīng)有6人回復(fù)
- VASP NEB 不收斂問題
已經(jīng)有25人回復(fù)
- VASP 計算石墨烯與石墨烯帶的區(qū)別
已經(jīng)有11人回復(fù)
- Siesta從官網(wǎng)下載贗勢計算出錯
已經(jīng)有7人回復(fù)
- VASP計算表面吸附 但是能量總是不收斂/。。。
已經(jīng)有9人回復(fù)
- 求助 高斯計算過程中一直不收斂 scf(maxcyc=500,xqc)
已經(jīng)有9人回復(fù)
- 算例第一個時間步就不收斂？
已經(jīng)有14人回復(fù)
- 用Dmol3結(jié)構(gòu)優(yōu)化時不收斂
已經(jīng)有5人回復(fù)
- 被積函數(shù)很快收斂到0，為什么積分不收斂。
已經(jīng)有3人回復(fù)
- siesta結(jié)構(gòu)優(yōu)化收斂標準
已經(jīng)有5人回復(fù)
- 結(jié)構(gòu)優(yōu)化 siesta 和 vasp
已經(jīng)有7人回復(fù)
- siesta-3.1 計算石墨烯條帶不收斂，不知是不是我的輸入有問題。
已經(jīng)有24人回復(fù)
- SIESTA計算能級不準怎么解決
已經(jīng)有10人回復(fù)
- 求助關(guān)于用vasp計算二維石墨烯納米帶所得到的能帶結(jié)果與報道對不上
已經(jīng)有7人回復(fù)
- vasp計算石墨烯能帶時，dirac點是打開的，帶隙不為零，與靜態(tài)的DOSCAR也對不上
已經(jīng)有10人回復(fù)
- 【求助】請問用siesta算輸運怎么建模啊 我是搞石墨烯
已經(jīng)有7人回復(fù)
- 【求助】siesta中transiesta模塊的收斂標準是啥~~
已經(jīng)有20人回復(fù)
- 【求助】siesta中transiesta模塊的收斂問題
已經(jīng)有13人回復(fù)
- 【求助】siesta中關(guān)于幾何結(jié)構(gòu)優(yōu)化的問題【完結(jié)】
已經(jīng)有17人回復(fù)
xwu10
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 599.6
- 散金: 228
- 帖子: 239
- 在線: 273小時
- 蟲號: 1366976
- 注冊: 2011-08-12
- 專業(yè): 金屬材料表面科學(xué)與工程
2樓2015-01-26 14:55:53
銳利的碎片
木蟲 (正式寫手)
star watcher
- 1ST強帖: 1
- 應(yīng)助: 136 (高中生)
- 金幣: 3637.1
- 散金: 252
- 紅花: 22
- 帖子: 988
- 在線: 1224.9小時
- 蟲號: 961933
- 注冊: 2010-03-05
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2015-01-26 21:09:01
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
4樓2015-01-26 21:42:12
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2015-01-26 21:45:10
guohuazhong
至尊木蟲 (職業(yè)作家)
- 1ST強帖: 1
- 應(yīng)助: 52 (初中生)
- 貴賓: 0.032
- 金幣: 35131.5
- 散金: 5
- 紅花: 19
- 沙發(fā): 1
- 帖子: 3924
- 在線: 885.6小時
- 蟲號: 550098
- 注冊: 2008-04-22
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
6樓2015-01-27 20:44:34
xwu10
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 599.6
- 散金: 228
- 帖子: 239
- 在線: 273小時
- 蟲號: 1366976
- 注冊: 2011-08-12
- 專業(yè): 金屬材料表面科學(xué)與工程
7樓2015-01-28 09:44:28
8樓2016-05-12 22:58:07
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xwu10 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-19 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 329求調(diào)劑 +6 | 想上學(xué)吖吖 2026-03-19 | 6/300 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +8 | 墨墨漠 2026-03-18 | 8/400 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +8 | 步川酷紫123 2026-03-13 | 8/400 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +20 | rare12345 2026-03-18 | 20/1000 |
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3 | Ncdx123456 2026-03-19 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|