 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3309 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
ydxia鐵蟲 (小有名氣)
|
[求助]
nvt系綜，系統(tǒng)總能量震蕩厲害
|
|
|
本人剛開始接觸lammps，作為練手，我計(jì)算了一條聚合度為10的聚乙烯分子鏈在NVT系綜下的動(dòng)力學(xué)過程。在materials studio里面建好模型后得到data文件，之后開始用lammps計(jì)算。下面是我的in文件： dimension 3 boundary p p p units real atom_style full #-------------------------------------- pair_style lj/charmm/coul/long 8.0 10.0 10.0 bond_style harmonic angle_style harmonic dihedral_style harmonic #improper_style cvff kspace_style pppm 0.0001 #-------------------------------------- read_data data.PE #-------------------------------------- neighbor 2.0 bin neigh_modify delay 0 every 1 check yes #-------------------------------------- min_style cg minimize 1.0e-4 1.0e-6 100 1000 # Perform an energy minimization of the system, by iteratively adjusting atom coordinates. #--------------------------------------- velocity all create 300.0 4928459 rot yes dist gaussian fix 1 all nvt temp 300.0 300.0 100.0 #-------------------------------------- timestep 2.0 # Set the timestep size for subsequent molecular dynamics simulations. #-------------------------------------- #compute 1 all gyration #-------------------------------------- dump 1 all custom 100 dump.PE id type x y z vx vy vz #-------------------------------------- thermo_style custom step temp press etotal pe ke emol ebond eangle edihed eimp epair evdwl ecoul elong etail thermo 100 #--------------------------------------- run 10000 本來應(yīng)該把我得到的溫度和系統(tǒng)能量隨時(shí)間的變化圖片加上來，但是總是添加失敗，只好文字描述一下了： 溫度在300上下震蕩，振幅很大；按理說系統(tǒng)的能量應(yīng)該是有一個(gè)降低的過程，但是我的結(jié)算結(jié)果卻是能量自始至終都在劇烈震蕩，根本沒有降低的趨度。 我把run的時(shí)間由10000增加至1000000，還是同樣地現(xiàn)象，應(yīng)該不是run的時(shí)間不夠。到底是什么原因呢？百思不得其解，大家?guī)臀页龀鲋饕獍桑x謝啦！ [ Last edited by ydxia on 2011-10-18 at 21:22 ] |
回復(fù)此樓» 猜你喜歡
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑
已經(jīng)有8人回復(fù)
-
349求調(diào)劑
已經(jīng)有6人回復(fù)
-
329求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有9人回復(fù)
-
327求調(diào)劑
已經(jīng)有5人回復(fù)
-
求收留
已經(jīng)有5人回復(fù)
-
375求調(diào)劑
已經(jīng)有4人回復(fù)
-
340求調(diào)劑
已經(jīng)有6人回復(fù)
-
298求調(diào)劑
已經(jīng)有3人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 金狐多功能系統(tǒng)維護(hù)盤 2012 簡(jiǎn)體中文版 （一盤在手裝機(jī)無憂）【轉(zhuǎn)載】
已經(jīng)有4人回復(fù)
- 美國Envia Systems實(shí)現(xiàn)400Wh/kg能量密度的鋰離子充電電池,Li2MnO3-LiMO2體系
已經(jīng)有15人回復(fù)
- 請(qǐng)問稻米食品方向的研究比較厲害的國家有？
已經(jīng)有22人回復(fù)
- 【求助】VASP做分子動(dòng)力學(xué)模擬NVT系綜下為啥內(nèi)能竟然是升高的？
已經(jīng)有8人回復(fù)
- 為什么要進(jìn)行能量最小化？
已經(jīng)有7人回復(fù)
- NVT系綜的英語全稱怎么寫？
已經(jīng)有3人回復(fù)
- 【求助】怎樣分析分子間相互作用體系的能量強(qiáng)弱
已經(jīng)有4人回復(fù)
- 【討論】lammps中怎樣檢驗(yàn)?zāi)M系統(tǒng)在NVT系綜下弛豫是否達(dá)到平衡？
已經(jīng)有6人回復(fù)
ydxia
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1448.9
- 散金: 60
- 紅花: 1
- 帖子: 145
- 在線: 62.3小時(shí)
- 蟲號(hào): 574179
- 注冊(cè): 2008-06-15
- 專業(yè): 高分子物理與高分子物理化
7樓2011-10-21 04:51:25
dinghuachen
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 7181.6
- 散金: 199
- 紅花: 3
- 帖子: 342
- 在線: 516.3小時(shí)
- 蟲號(hào): 1079849
- 注冊(cè): 2010-08-22
- 性別: GG
- 專業(yè): 高分子物理與高分子物理化
3樓2011-10-19 16:19:31
ydxia
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1448.9
- 散金: 60
- 紅花: 1
- 帖子: 145
- 在線: 62.3小時(shí)
- 蟲號(hào): 574179
- 注冊(cè): 2008-06-15
- 專業(yè): 高分子物理與高分子物理化
4樓2011-10-19 22:05:50
dinghuachen
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 7181.6
- 散金: 199
- 紅花: 3
- 帖子: 342
- 在線: 516.3小時(shí)
- 蟲號(hào): 1079849
- 注冊(cè): 2010-08-22
- 性別: GG
- 專業(yè): 高分子物理與高分子物理化
5樓2011-10-20 09:39:56
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +5 | 小懶蟲不懶了 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 309求調(diào)劑 +6 | 誰不是少年 2026-03-29 | 6/300 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +9 | 奶油草莓. 2026-03-22 | 10/500 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
|
[考研]
材料學(xué)碩，求調(diào)劑
 6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|