| 查看: 3596 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
ydxia铁虫 (小有名气)
|
[求助]
nvt系综，系统总能量震荡厉害
|
|
|
本人刚开始接触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 本来应该把我得到的温度和系统能量随时间的变化图片加上来，但是总是添加失败，只好文字描述一下了： 温度在300上下震荡，振幅很大；按理说系统的能量应该是有一个降低的过程，但是我的结算结果却是能量自始至终都在剧烈震荡，根本没有降低的趋度。 我把run的时间由10000增加至1000000，还是同样地现象，应该不是run的时间不够。到底是什么原因呢？百思不得其解，大家帮我出出主意吧，谢谢啦！ [ Last edited by ydxia on 2011-10-18 at 21:22 ] |
回复此楼» 猜你喜欢
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有11人回复
-
有没有大神帮我看看基金代码
已经有9人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
放榜前的不淡定
已经有18人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 金狐多功能系统维护盘 2012 简体中文版 （一盘在手装机无忧）【转载】
已经有4人回复
- 美国Envia Systems实现400Wh/kg能量密度的锂离子充电电池,Li2MnO3-LiMO2体系
已经有15人回复
- 请问稻米食品方向的研究比较厉害的国家有？
已经有22人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- 为什么要进行能量最小化？
已经有7人回复
- NVT系综的英语全称怎么写？
已经有3人回复
- 【求助】怎样分析分子间相互作用体系的能量强弱
已经有4人回复
- 【讨论】lammps中怎样检验模拟系统在NVT系综下弛豫是否达到平衡？
已经有6人回复
ydxia
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1448.9
- 散金: 60
- 红花: 1
- 帖子: 145
- 在线: 62.3小时
- 虫号: 574179
- 注册: 2008-06-15
- 专业: 高分子物理与高分子物理化
9楼2011-10-21 04:56:18
dinghuachen
铁杆木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 7181.6
- 散金: 199
- 红花: 3
- 帖子: 342
- 在线: 516.3小时
- 虫号: 1079849
- 注册: 2010-08-22
- 性别: GG
- 专业: 高分子物理与高分子物理化
3楼2011-10-19 16:19:31
ydxia
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1448.9
- 散金: 60
- 红花: 1
- 帖子: 145
- 在线: 62.3小时
- 虫号: 574179
- 注册: 2008-06-15
- 专业: 高分子物理与高分子物理化
4楼2011-10-19 22:05:50
dinghuachen
铁杆木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 7181.6
- 散金: 199
- 红花: 3
- 帖子: 342
- 在线: 516.3小时
- 虫号: 1079849
- 注册: 2010-08-22
- 性别: GG
- 专业: 高分子物理与高分子物理化
5楼2011-10-20 09:39:56
20