| 5 | 1/1 | 返回列表 |
| 查看: 3306 | 回復(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ī)臀页龀鲋饕獍,謝謝啦! [ Last edited by ydxia on 2011-10-18 at 21:22 ] |
鐵蟲 (小有名氣)
鐵桿木蟲 (正式寫手)

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。。 +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|