| 查看: 4424 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
WJ-MS新虫 (小有名气)
|
[求助]
LAMMPS温度－步长曲线震荡
|
|
|
有几个问题想跟大家请教一下。 1 从网上得到的一个程序（见附件），跑了一下，发现体系的温度一直上升，但是在设定初始温度的时候不是设置的10K么，体系的温度是不是通过 “velocity all create 1000 87287 dist gaussian“ 这个命令设定阿？随后的温度变化又由什么决定的呢？ 2 想计算1000K下的性质，采用“velocity all create 1000 87287 dist gaussia“， 计算时，先fix npt运算50000步，再fix nvt运算500000步，体系的温度随时间的变化不是应该趋于1000K么，为什么我得到的温度曲线是震荡的呢？ 3 此外，RDF曲线，我看别人都是一条曲线，我得到的怎么是很多条曲线的样子？这个该怎么处理呢？具体看附件。 菜鸟级，希望大家多多指教！  温度时间步长.png  rdf.png |
回复此楼» 猜你喜欢
277求调剂
已经有20人回复
-
一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣
已经有10人回复
-
279学硕食品专业求调剂院校
已经有11人回复
-
调剂
已经有26人回复
-
生物学308求调剂
已经有4人回复
-
食品与营养（0955）271求调剂
已经有3人回复
-
290求调剂
已经有14人回复
-
335求调剂
已经有13人回复
-
279求调剂
已经有6人回复
-
322求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps的时间步长问题
已经有16人回复
- 在LAMMPS中怎么控制温度变化速率
已经有4人回复
- nvt系综，系统总能量震荡厉害
已经有9人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
- 【讨论】高温高压的体系需要注意的环节
已经有18人回复
【答案】应助回帖
★ ★ ★ ★ ★
WJ-MS: 金币+5, ★★★很有帮助 2013-05-28 16:17:52
WJ-MS: 金币+5, ★★★很有帮助 2013-05-28 16:17:52
| 据我了解，velocity all create 1000 87287 dist gaussian是一个系统初始化的命令，是根据你所提供的系统温度，提供系统动能，再根据动能为各原子随机分配速度，并不适用于恒温操作。fix temp/rescale可以实现温度在一定误差范围内进行波动，是最常用的控温方法。而我们所定义的npt或nvt是系综，主要针对系统的物理类型进行计算，调温时还是应该同时定义调温方法的，不能定义npt后将fix temp/rescale命令delete掉，我也问一个问题，就是你所说的这个软件是什么软件，可以发给我一下么，我之前一直是自己编程画图，很是麻烦，谢谢，邮箱393831415@qq.com |
mileyfan9楼2013-05-28 12:50:13
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
2楼2013-05-03 13:06:12
WJ-MS
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1559.2
- 散金: 100
- 帖子: 194
- 在线: 92.5小时
- 虫号: 1697144
- 注册: 2012-03-17
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
3楼2013-05-04 09:11:07
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
4楼2013-05-05 08:15:57
