 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 2680 | 回復: 9 | |||
xiaowu787木蟲 (正式寫手)
|
[交流]
【討論】lammps-怪異的group溫度 已有4人參與
|
|
說明一下，雙金屬體系，時間設為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分別是從模型中心到模型一邊的邊界---相當于中心對稱） 體系溫度 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 ] |
回復此樓» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有5人回復
-
281求調(diào)劑（0805）
已經(jīng)有16人回復
-
304求調(diào)劑
已經(jīng)有6人回復
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復
-
求材料調(diào)劑
已經(jīng)有8人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復
-
一志愿華中科技大學，080502，354分求調(diào)劑
已經(jīng)有4人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有6人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- MS生成Lammps的data文件，文件中原子的類型和順序并不是自己想要的，請問如何修改？
已經(jīng)有8人回復
- 求助LAMMPS關于MSD計算的問題
已經(jīng)有17人回復
- LAMMPS溫度－步長曲線震蕩
已經(jīng)有11人回復
- lammps的in文件
已經(jīng)有4人回復
- lammps deposit 粒子初始速度為0？
已經(jīng)有12人回復
- LAMMPS中meam勢無法調(diào)用
已經(jīng)有10人回復
- 怎么樣開始LAMMPS的學習呢
已經(jīng)有18人回復
- lammps的時間步長問題
已經(jīng)有16人回復
- LAMMPS Error:lost atoms
已經(jīng)有7人回復
- 在LAMMPS中怎么建立圖示的模型
已經(jīng)有6人回復
- 【求助】LAMMPS運行錯誤：Out of range atoms - cannot compute PPPM
已經(jīng)有7人回復
- 【求助】fedora中LAMMPS編譯問題
已經(jīng)有4人回復
- 【分享】嘗試lammps中, 分享中...
已經(jīng)有15人回復
- 【討論】lammps中利用velocity命令無法控制住速度問題
已經(jīng)有11人回復
- 【討論】lammps中region和group的問題
已經(jīng)有4人回復
- 【討論】lammps中關于fix ave/spatial命令的輸出問題
已經(jīng)有4人回復
- 【討論】lammps弛豫過程離奇的快
已經(jīng)有4人回復
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應助: 26 (小學生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時
- 蟲號: 659094
- 注冊: 2008-11-21
- 專業(yè): 金屬結(jié)構材料
2樓2010-10-30 19:31:54
xiaowu787
木蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 1406.9
- 散金: 95
- 紅花: 1
- 帖子: 550
- 在線: 430.3小時
- 蟲號: 326816
- 注冊: 2007-03-18
- 性別: GG
3樓2010-10-31 09:28:07
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應助: 26 (小學生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時
- 蟲號: 659094
- 注冊: 2008-11-21
- 專業(yè): 金屬結(jié)構材料
★ ★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
ghcacj(金幣+6):謝謝 2010-11-01 09:32:08
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
ghcacj(金幣+6):謝謝 2010-11-01 09:32:08
|
你的九個部分原子數(shù)相等么？ A 如果原子數(shù)相等，且溫度是對稱的，那就意味著九個部分的平均溫度都高于總平均溫度。這個情況應是不可能出現(xiàn)的。所以溫度可能不是對稱的。 B 假如你九個部分原子數(shù)不相等，則情況比較復雜。 實際上，你在MD當中，分成九個部分，即使原子數(shù)相等，且初始結(jié)構對稱，但你的初速度的方向和大小都是隨機的，尤其是元胞中的中間部分經(jīng)過運動變化后，也不會是完全不偏不倚的，所以溫度一般不會完全對稱。一邊高一邊低的可能性不能排除，這與你的邊界條件也有關系。 [ Last edited by 老虎大王 on 2010-10-31 at 12:16 ] |
4樓2010-10-31 12:12:26
xiaowu787
木蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 1406.9
- 散金: 95
- 紅花: 1
- 帖子: 550
- 在線: 430.3小時
- 蟲號: 326816
- 注冊: 2007-03-18
- 性別: GG
★
zh1987hs(金幣+1):鼓勵交流 2010-11-02 20:24:48
zh1987hs(金幣+1):鼓勵交流 2010-11-02 20:24:48
|
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 謝謝您的答復，現(xiàn)在還有一點不明白，溫度不是統(tǒng)計平均值嗎？為什么還與group原子數(shù)目有關？希望您能再指點一下？ |
5樓2010-11-02 17:44:29
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應助: 26 (小學生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時
- 蟲號: 659094
- 注冊: 2008-11-21
- 專業(yè): 金屬結(jié)構材料
★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zh1987hs(金幣+5):謝謝老虎兄 2010-11-02 20:25:19
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
zh1987hs(金幣+5):謝謝老虎兄 2010-11-02 20:25:19
|
我的意思不是說溫度與原子數(shù)目有關，而是當你考慮總平均溫度（即所有原子的總動能對總原子數(shù)的商 Ek(tot)/N），與各部分平均溫度的平均（1/9*(Ek1/N1+Ek2/N2+...)），它們兩個之間的關系時，這個關系（是相等，還是大于，還是小于，這個應該是與各部分原子數(shù)是否相等有關的 （上述兩式是個大致的意思，不是溫度的確切表達，沒有寫上有關的系數(shù)） 總之我基本上斷定你兩邊的溫度不會是對稱的。 [ Last edited by 老虎大王 on 2010-11-2 at 19:32 ] |
6樓2010-11-02 19:23:22
xiaowu787
木蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 1406.9
- 散金: 95
- 紅花: 1
- 帖子: 550
- 在線: 430.3小時
- 蟲號: 326816
- 注冊: 2007-03-18
- 性別: GG
★
zh1987hs(金幣+1):鼓勵交流 2010-11-03 09:56:57
zh1987hs(金幣+1):鼓勵交流 2010-11-03 09:56:57
|
模型結(jié)構本來是均勻的，不過分成九個部分時人為造成group原子數(shù)目的不同，原來想當然的認為各部分動能平均是一樣，那么溫度也應該一樣。 下面是分為7個部分（模型小了一些）的結(jié)果（注：時間步長不小心設成了10fs--主要區(qū)別) group temp0 region temp0 55180 atoms in group temp group temp1 region temp1 55177 atoms in group temp1 group temp2 region temp2 55124 atoms in group temp2 group temp3 region temp3 56290 atoms in group temp3 可惜沒有記錄之前的，第一行數(shù)據(jù)是在300K下馳豫后的結(jié)果，溫度倒是沒有表現(xiàn)出很大的波動 體系溫度 group1 group2 group3 group4 300.69665 302.18029 301.84022 301.07516 298.5104 313.3621 393.59714 302.1085 300.1085 297.70397 334.20323 497.64952 306.07788 300.52586 298.067 359.62189 603.66492 316.23812 303.46189 298.99988 387.76138 707.79523 330.64041 304.5327 299.22969 417.99158 816.07707 349.605 311.21851 301.95479 450.30771 918.12848 372.49825 318.89365 306.08764 483.03501 1026.523 397.37457 328.27153 307.04527 517.06064 1128.7104 426.83542 338.56959 311.56416 553.46222 1235.7885 455.4556 352.70137 320.88119 590.08605 1335.7681 488.29725 366.44313 327.79355 627.75478 1448.4179 520.09948 386.86076 336.64048 667.25754 1549.6889 551.42325 402.66075 347.6617 707.2067 1653.0293 578.8763 420.80881 360.50564 748.4945 1769.3749 617.00291 442.83465 374.4264 788.56642 1877.2389 658.98013 461.85488 390.24975 [ Last edited by xiaowu787 on 2010-11-2 at 22:32 ] |
7樓2010-11-02 22:22:30
老虎大王
木蟲 (著名寫手)
- 模擬EPI: 17
- 應助: 26 (小學生)
- 貴賓: 0.17
- 金幣: 4774.1
- 散金: 8
- 紅花: 42
- 帖子: 1361
- 在線: 215.2小時
- 蟲號: 659094
- 注冊: 2008-11-21
- 專業(yè): 金屬結(jié)構材料
8樓2010-11-03 07:52:50
奈落之夜
木蟲 (正式寫手)
- 應助: 6 (幼兒園)
- 金幣: 1890.6
- 紅花: 1
- 帖子: 393
- 在線: 73.9小時
- 蟲號: 2550025
- 注冊: 2013-07-17
- 專業(yè): 凝聚態(tài)物性I:結(jié)構、力學和
9樓2014-02-27 10:35:29
10樓2015-10-11 11:13:38
相關版塊跳轉(zhuǎn)
我要訂閱樓主 xiaowu787 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 0703化學求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 0703 物理化學調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 學碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考博] 福州大學楊黃浩課題組招收2026年專業(yè)學位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|