| 查看: 3614 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有15人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 金狐多功能系统维护盘 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
30