| 查看: 3291 | 回复: 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 ] |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 金狐多功能系统维护盘 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
- 专业: 高分子物理与高分子物理化
2楼2011-10-18 21:22:35
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
dinghuachen
铁杆木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 7181.6
- 散金: 199
- 红花: 3
- 帖子: 342
- 在线: 516.3小时
- 虫号: 1079849
- 注册: 2010-08-22
- 性别: GG
- 专业: 高分子物理与高分子物理化
6楼2011-10-20 10:44:45
ydxia
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1448.9
- 散金: 60
- 红花: 1
- 帖子: 145
- 在线: 62.3小时
- 虫号: 574179
- 注册: 2008-06-15
- 专业: 高分子物理与高分子物理化
7楼2011-10-21 04:51:25
ydxia
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1448.9
- 散金: 60
- 红花: 1
- 帖子: 145
- 在线: 62.3小时
- 虫号: 574179
- 注册: 2008-06-15
- 专业: 高分子物理与高分子物理化
8楼2011-10-21 04:52:29
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
- 专业: 高分子物理与高分子物理化
10楼2011-10-21 11:13:14
