| 10 | 1/1 | 返回列表 |
| 查看: 3305 | 回復: 9 | ||
ydxia鐵蟲 (小有名氣)
|
[求助]
nvt系綜,系統(tǒng)總能量震蕩厲害
|
|
本人剛開始接觸lammps,作為練手,我計算了一條聚合度為10的聚乙烯分子鏈在NVT系綜下的動力學過程。在materials studio里面建好模型后得到data文件,之后開始用lammps計算。下面是我的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 本來應該把我得到的溫度和系統(tǒng)能量隨時間的變化圖片加上來,但是總是添加失敗,只好文字描述一下了: 溫度在300上下震蕩,振幅很大;按理說系統(tǒng)的能量應該是有一個降低的過程,但是我的結算結果卻是能量自始至終都在劇烈震蕩,根本沒有降低的趨度。 我把run的時間由10000增加至1000000,還是同樣地現(xiàn)象,應該不是run的時間不夠。到底是什么原因呢?百思不得其解,大家?guī)臀页龀鲋饕獍,謝謝啦! [ Last edited by ydxia on 2011-10-18 at 21:22 ] |
鐵蟲 (小有名氣)
鐵桿木蟲 (正式寫手)

鐵蟲 (小有名氣)
鐵桿木蟲 (正式寫手)

鐵桿木蟲 (正式寫手)

鐵蟲 (小有名氣)
鐵蟲 (小有名氣)
鐵蟲 (小有名氣)
鐵桿木蟲 (正式寫手)

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305求調劑 +5 | RuiFairyrui 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 本科新能源科學與工程,一志愿華理能動285求調劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 一志愿華北電力大學能動專碩,293,求調劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
[考研] 085602 化工專碩 338分 求調劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 一志愿中南大學化學0703總分337求調劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 346求調劑 一志愿070303有機化學 +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-24 | 4/200 |
|
|
[考研] 299求調劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 266求調劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 314求調劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 0703化學一志愿南京師范大學303求調劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 324求調劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 336材料求調劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 349求調劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 環(huán)境專碩324分求調劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 293求調劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 材料調劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 086003食品工程求調劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研] 300求調劑,材料科學英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
[考研] 求調劑一志愿武漢理工大學材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|