| 5 | 1/1 | 返回列表 |
| 查看: 2705 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
xiaowu787木蟲 (正式寫手)
|
[交流]
【討論】lammps-怪異的group溫度 已有4人參與
|
||
|
說明一下,雙金屬體系,時(shí)間設(shè)為1fs,體系整體從0.05到300K升溫,發(fā)現(xiàn)group的溫度與體系的溫度對(duì)不上,并且處于模型最外端的group的溫度更是異常,比其他group的溫度高幾十度。 原來(lái)不小心把時(shí)間步長(zhǎng)取為10fs,進(jìn)行模擬同樣的模擬也沒有出現(xiàn)類似的問題。并且可以認(rèn)定的是勢(shì)函數(shù)是可靠的。不知是哪里的問題? 下面是數(shù)據(jù): 時(shí)間步長(zhǎng)為0.001 單位為metal fix 2 all nvt temp 0.05 300.0 0.05 第一個(gè)數(shù)是整個(gè)體系的溫度,后面五個(gè)數(shù)據(jù)分別是五個(gè)group的溫度(把長(zhǎng)方體模型均分為9個(gè)部分,取的五個(gè)group分別是從模型中心到模型一邊的邊界---相當(dāng)于中心對(duì)稱) 體系溫度 group1 group2 group3 group4 group5 0.32223265 0.35672397 0.35130176 0.43155934 0.49054604 0.59470303 14.681225 17.707719 18.038659 19.401222 18.578224 23.696852 17.908917 22.595134 22.381331 22.859198 22.735383 29.888281 21.236465 26.594786 27.373716 27.549981 25.78319 33.537153 24.658653 32.005349 29.70771 29.110636 29.385247 40.579523 28.866911 35.896002 36.528466 36.668734 37.459447 45.291189 33.150806 40.469357 40.669588 40.651187 42.624405 55.096068 38.69363 47.161347 47.656391 48.248048 49.56008 65.557493 43.565413 57.637558 52.837027 55.77367 52.481708 71.654626 295.43779 314.76351 305.85202 324.41969 308.7731 386.26485 299.80219 313.80556 316.32717 320.05518 311.6842 389.04083 [ Last edited by xiaowu787 on 2010-10-31 at 09:16 ] |
木蟲 (正式寫手)
木蟲 (著名寫手)
木蟲 (著名寫手)
|
你的九個(gè)部分原子數(shù)相等么? A 如果原子數(shù)相等,且溫度是對(duì)稱的,那就意味著九個(gè)部分的平均溫度都高于總平均溫度。這個(gè)情況應(yīng)是不可能出現(xiàn)的。所以溫度可能不是對(duì)稱的。 B 假如你九個(gè)部分原子數(shù)不相等,則情況比較復(fù)雜。 實(shí)際上,你在MD當(dāng)中,分成九個(gè)部分,即使原子數(shù)相等,且初始結(jié)構(gòu)對(duì)稱,但你的初速度的方向和大小都是隨機(jī)的,尤其是元胞中的中間部分經(jīng)過運(yùn)動(dòng)變化后,也不會(huì)是完全不偏不倚的,所以溫度一般不會(huì)完全對(duì)稱。一邊高一邊低的可能性不能排除,這與你的邊界條件也有關(guān)系。 [ Last edited by 老虎大王 on 2010-10-31 at 12:16 ] |
木蟲 (正式寫手)
|
group temp0 region temp0 47528 atoms in group temp0 group temp1 region temp1 47534 atoms in group temp1 group temp2 region temp2 47525 atoms in group temp2 group temp3 region temp3 47518 atoms in group temp3 group temp4 region temp4 48582 atoms in group temp4 謝謝您的答復(fù),現(xiàn)在還有一點(diǎn)不明白,溫度不是統(tǒng)計(jì)平均值嗎?為什么還與group原子數(shù)目有關(guān)?希望您能再指點(diǎn)一下? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工328分調(diào)劑 +8 | 。,。,。,。i 2026-03-23 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +4 | 兔兔兔111223314 2026-03-29 | 4/200 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過六級(jí)。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|