| 5 | 1/1 | 返回列表 |
| 查看: 2703 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
xiaowu787木蟲 (正式寫手)
|
[交流]
【討論】lammps-怪異的group溫度 已有4人參與
|
||
|
說明一下,雙金屬體系,時間設(shè)為1fs,體系整體從0.05到300K升溫,發(fā)現(xiàn)group的溫度與體系的溫度對不上,并且處于模型最外端的group的溫度更是異常,比其他group的溫度高幾十度。 原來不小心把時間步長取為10fs,進行模擬同樣的模擬也沒有出現(xiàn)類似的問題。并且可以認定的是勢函數(shù)是可靠的。不知是哪里的問題? 下面是數(shù)據(jù): 時間步長為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ù)雜。 實際上,你在MD當(dāng)中,分成九個部分,即使原子數(shù)相等,且初始結(jié)構(gòu)對稱,但你的初速度的方向和大小都是隨機的,尤其是元胞中的中間部分經(jīng)過運動變化后,也不會是完全不偏不倚的,所以溫度一般不會完全對稱。一邊高一邊低的可能性不能排除,這與你的邊界條件也有關(guān)系。 [ Last edited by 老虎大王 on 2010-10-31 at 12:16 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 各位老師您好:本人初試372分 +4 | jj涌77 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +6 | 顧九笙要謙虛 2026-03-24 | 6/300 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +17 | 陌の森林 2026-03-18 | 17/850 |
|
|
[考研] 總分293求調(diào)劑 +3 | 加一一九 2026-03-25 | 3/150 |
|
|
[考研] 298求調(diào)劑 +3 | Rong$$$ 2026-03-18 | 3/150 |
|
|
[考研] 食品專碩 一志愿雙一流 328 +3 | xiaom99 2026-03-21 | 4/200 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +6 | NIFFFfff 2026-03-20 | 6/300 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 求材料,環(huán)境專業(yè)調(diào)劑 +3 | 18567500178 2026-03-18 | 3/150 |
|
|
[考研] 335分 | 材料與化工專碩 | GPA 4.07 | 有科研經(jīng)歷 +4 | cccchenso 2026-03-23 | 4/200 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|