| 10 | 1/1 | 返回列表 |
| 查看: 1529 | 回復: 9 | |||
[交流]
【求助】SIESTA計算總能等自洽計算時力準不準阿 已有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)化后的受力,已經達收斂標準了 我用這個結構,再做scf計算 發(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 可以看出,自洽計算時的輸出的力和優(yōu)化后輸出的力不一樣阿 現(xiàn)在的問題:是不是 SIESTA計算總能等自洽計算時力算不準呢 |
木蟲 (著名寫手)

|
這就是自洽最后的結果: 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 ] |
木蟲 (著名寫手)

木蟲 (著名寫手)

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

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專碩調劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 求調劑,一志愿:南京航空航天大學大學 ,080500材料科學與工程學碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 求調劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] A區(qū)線材料學調劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 一志愿北京化工大學0703化學318分,有科研經歷,求調劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
[考研] 081700化工學碩調劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|
[考研] 0703化學調劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 286求調劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 材料考研調劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 328求調劑,英語六級551,有科研經歷 +4 | 生物工程調劑 2026-03-16 | 12/600 |
|
|
[考研] 化學工程321分求調劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 278求調劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 326求調劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 機械專碩325,尋找調劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 318求調劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 288求調劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
[考研] 復試調劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|