 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1758 | 回復: 11 | ||||
| 本帖產生 2 個 1ST強帖 ，點擊這里進行查看 | ||||
輕風龍騎士新蟲 (初入文壇)
|
[交流]
【求助】siesta 做NOSE動力學出錯，請高手指教 已有4人參與
|
|||
| 在做nose的溫控動力學的時候，我設置的初始和目標溫度分別是2700K和2600K，做降溫的動力學模擬，但是我的systemlable.MDE文件中怎么顯示的每一步的溫度有問題�。涸诘谌降臅r候就已經到了2620K左右，之后溫度更低最低只有1400K左右，然后又開始升高，最高高達3000K以上，這是為什么？是不是我那個參數錯了，請高手指教 |
回復此樓» 收錄本帖的淘帖專輯推薦
 siesta問題 |
» 猜你喜歡
考研求調劑
已經有4人回復
-
085601專碩，總分342求調劑，地區(qū)不限
已經有4人回復
-
085601材料工程專碩求調劑
已經有5人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調劑
已經有6人回復
-
能源材料化學課題組招收碩士研究生8-10名
已經有6人回復
-
296求調劑
已經有5人回復
-
0703化學調劑 ，六級已過，有科研經歷
已經有8人回復
-
268求調劑
已經有8人回復
-
材料，紡織，生物（0856、0710），化學招生啦
已經有7人回復
-
301求調劑
已經有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求助萊克多巴胺的動力學特征
已經有5人回復
- 麻煩哪位高手幫忙翻譯一下，后半段實在是搞不懂，我急著做PPT的，謝謝了。
已經有1人回復
- 請問考博調劑的事，過來的前輩請指教�。�1!
已經有26人回復
- 出錯 warning#2 #180
已經有7人回復
- 請教高手指點迷津－固體磷核磁表征固體酸酸強度
已經有5人回復
- 關于Discover動力學計算的一個問題
已經有26人回復
- 請教：orca 中用CCSD(T)計算單點能出錯--提示內存不足 該怎么辦？
已經有10人回復
- 請教：管式積分反應器做動力學時，反應速率如何計算？
已經有10人回復
- 原來編譯的siesta突然不能用，重新編譯沒錯，運行出錯
已經有13人回復
- 【求助】siesta進行ordern計算出錯求助
已經有4人回復
- 【求助】誰做ZSM-5分子篩，請高手指教
已經有24人回復
- 【求助】安裝siesta3.0出錯：segmentation violation signal raised。請高手指教
已經有8人回復
- 【求助】用charmm軟件做分子動力學模擬是出現問題請指教
已經有4人回復
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強帖: 3
- 應助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
2樓2011-04-12 23:19:52
輕風龍騎士
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 393.8
- 帖子: 46
- 在線: 48小時
- 蟲號: 947650
- 注冊: 2010-01-24
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2011-04-13 12:52:16
wangsong1016
銅蟲 (初入文壇)
- 1ST強帖: 2
- 應助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時
- 蟲號: 408613
- 注冊: 2007-06-21
- 性別: GG
- 專業(yè): 物理
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
youzhizhe(金幣+2, 1ST強帖+1): 謝謝交流。 2011-05-02 19:56:53
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
youzhizhe(金幣+2, 1ST強帖+1): 謝謝交流。 2011-05-02 19:56:53
|
我沒有用過nose模塊，不過我在用verletMD計算時出現過類似的溫度起伏，甚至是劇烈起伏的問題。 我想可能有如下的原因給你參考。 首先宏觀體系的溫度就如二樓所說是個統計平均值，實際上是對分子熱運動的一個平均描述。 到了微觀層次，按照能量均分定理，每個原子，每個自由度能量是1/2kBT,簡化考慮就可以采用1/2kBT=1/2mv**2。意思是在一維無轉動的情況下，這一個自由度的能量對應于粒子平動的動能。在程序中就是這么計算得到的溫度。輸出文件中有一個systemlable.XV里面有各個原子不同方向的速度。通過公式就可以把各個速度對應的能量轉化為溫度。 你在輸入文件中給定了初始的溫度，程序就給各個原子不同的速度讓其震動起來。但是原子一旦發(fā)生震動，勢能必然變化，那么各個原子的動能（速度）也隨之變化，因為總能量要保持守恒。這樣就導致了體系溫度的變化。這種漲落是不可避免的，不算是錯誤，但是可能與我們想要的所謂恒定溫度情況有較大的偏差。 |
4樓2011-04-16 22:25:16
輕風龍騎士
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 393.8
- 帖子: 46
- 在線: 48小時
- 蟲號: 947650
- 注冊: 2010-01-24
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
5樓2011-05-02 14:59:18
wangsong1016
銅蟲 (初入文壇)
- 1ST強帖: 2
- 應助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時
- 蟲號: 408613
- 注冊: 2007-06-21
- 性別: GG
- 專業(yè): 物理
★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
youzhizhe(金幣+2, 1ST強帖+1): 謝謝交流。 2011-05-03 10:48:12
小木蟲(金幣+0.5):給個紅包，謝謝回帖
youzhizhe(金幣+2, 1ST強帖+1): 謝謝交流。 2011-05-03 10:48:12
|
恒溫的動力學過程我沒有計算過，所以對過程不太了解。至于10%是否合適我覺得也很難界定，這根據你的體系的大小不同應該是不一樣的。體系越小，漲落就越大。計算50個原子時漲落應該大于500個原子，兩個體系都用10%應該是不合適的。我這里也不敢給你定量的結論，很不好意思。 還有一個是我認為你所說的相變應該是宏觀的一個過程，某個溫度下的穩(wěn)定相應該是在固體的宏觀尺度范圍內進行考慮，但是看固體的微觀結構，還是應該有很劇烈的漲落的。計算到百十個原子的時候，無論漲落是否很大應該都很難反映出宏觀固體狀態(tài)。 我有個想法就是如果你想做一個恒溫的計算，可否考慮計算時間長一點。比如你要一個100K的溫度，你讓體系自己運行，跑個2000或者3000個MD過程，然后再看一下他的溫度變化范圍，有可能就是在你100K的溫度上下有較大的波動。但是你把3000個溫度取個平均值試試，也許就是100K左右的溫度。溫度本來就是一個較長時間對能量統計平均值，如果做了這么長那個時間的模擬，溫度平均值能夠是你所希望得到的那個值，我覺得應該是比較有說服力了。 |
6樓2011-05-03 10:34:49
輕風龍騎士
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 393.8
- 帖子: 46
- 在線: 48小時
- 蟲號: 947650
- 注冊: 2010-01-24
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
7樓2011-05-04 11:17:45
wangsong1016
銅蟲 (初入文壇)
- 1ST強帖: 2
- 應助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時
- 蟲號: 408613
- 注冊: 2007-06-21
- 性別: GG
- 專業(yè): 物理
8樓2011-05-05 09:29:59
輕風龍騎士
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 393.8
- 帖子: 46
- 在線: 48小時
- 蟲號: 947650
- 注冊: 2010-01-24
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
9樓2011-05-05 23:41:13
wangsong1016
銅蟲 (初入文壇)
- 1ST強帖: 2
- 應助: 0 (幼兒園)
- 金幣: 169
- 散金: 4
- 帖子: 45
- 在線: 26.3小時
- 蟲號: 408613
- 注冊: 2007-06-21
- 性別: GG
- 專業(yè): 物理
10樓2011-05-06 13:25:50
相關版塊跳轉
我要訂閱樓主 輕風龍騎士 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學調劑 ，六級已過，有科研經歷 +8 | 曦熙兮 2026-03-15 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 301求調劑 +5 | yy要上岸呀 2026-03-17 | 5/250 |
|
|
|
[考研] 293求調劑 +7 | zjl的號 2026-03-16 | 12/600 |
|
|
|
[考研] 304求調劑 +8 | 小熊joy 2026-03-14 | 8/400 |
|
|
|
[考研] 考研化學學碩調劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 化學工程321分求調劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
|
[考研] 290求調劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 材料與化工專碩調劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 302求調劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +7 | 邱gl 2026-03-11 | 8/400 |
|
|
|
[考研] 278求調劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 化學調劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 277材料科學與工程080500求調劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 294求調劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 【0703化學調劑】-一志愿華中師范大學-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[考研] 0856材料與化工301求調劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 333求調劑 +3 | 152697 2026-03-12 | 4/200 |
|