| 5 | 1/1 | 返回列表 |
| 查看: 2704 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
xiaowu787木蟲 (正式寫手)
|
[交流]
【討論】lammps-怪異的group溫度 已有4人參與
|
||
|
說明一下,雙金屬體系,時(shí)間設(shè)為1fs,體系整體從0.05到300K升溫,發(fā)現(xiàn)group的溫度與體系的溫度對不上,并且處于模型最外端的group的溫度更是異常,比其他group的溫度高幾十度。 原來不小心把時(shí)間步長取為10fs,進(jìn)行模擬同樣的模擬也沒有出現(xiàn)類似的問題。并且可以認(rèn)定的是勢函數(shù)是可靠的。不知是哪里的問題? 下面是數(shù)據(jù): 時(shí)間步長為0.001 單位為metal fix 2 all nvt temp 0.05 300.0 0.05 第一個數(shù)是整個體系的溫度,后面五個數(shù)據(jù)分別是五個group的溫度(把長方體模型均分為9個部分,取的五個group分別是從模型中心到模型一邊的邊界---相當(dāng)于中心對稱) 體系溫度 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 ] |
木蟲 (著名寫手)
|
你的九個部分原子數(shù)相等么? A 如果原子數(shù)相等,且溫度是對稱的,那就意味著九個部分的平均溫度都高于總平均溫度。這個情況應(yīng)是不可能出現(xiàn)的。所以溫度可能不是對稱的。 B 假如你九個部分原子數(shù)不相等,則情況比較復(fù)雜。 實(shí)際上,你在MD當(dāng)中,分成九個部分,即使原子數(shù)相等,且初始結(jié)構(gòu)對稱,但你的初速度的方向和大小都是隨機(jī)的,尤其是元胞中的中間部分經(jīng)過運(yùn)動變化后,也不會是完全不偏不倚的,所以溫度一般不會完全對稱。一邊高一邊低的可能性不能排除,這與你的邊界條件也有關(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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +5 | 零八# 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 08工科求調(diào)劑286 +3 | tgs_001 2026-03-28 | 3/150 |
|
|
[基金申請] 面上5B能上會嗎? +4 | redcom 2026-03-29 | 4/200 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過六級。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +6 | minht 2026-03-29 | 6/300 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +8 | 麻辣魷魚 2026-03-27 | 9/450 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。。 +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|