| 查看: 4374 | 回复: 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 |
回复此楼» 猜你喜欢
环境调剂
已经有4人回复
-
272求调剂
已经有9人回复
-
290分材料工程085601求调剂 数二英一
已经有5人回复
-
0856材料求调剂
已经有13人回复
-
材料学硕318求调剂
已经有10人回复
-
材料270求调剂
已经有7人回复
-
0856求调剂285
已经有11人回复
-
085600材料工程一志愿中科大总分312求调剂
已经有11人回复
-
材料284求调剂，一志愿郑州大学英一数二专硕
已经有15人回复
-
理学，工学，农学调剂，少走弯路，这里欢迎您！
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps的时间步长问题
已经有16人回复
- 在LAMMPS中怎么控制温度变化速率
已经有4人回复
- nvt系综，系统总能量震荡厉害
已经有9人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
- 【讨论】高温高压的体系需要注意的环节
已经有18人回复
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
WJ-MS2楼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
WJ-MS
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1559.2
- 散金: 100
- 帖子: 194
- 在线: 92.5小时
- 虫号: 1697144
- 注册: 2012-03-17
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2013-05-06 20:39:09
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
6楼2013-05-07 08:44:12
WJ-MS
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1559.2
- 散金: 100
- 帖子: 194
- 在线: 92.5小时
- 虫号: 1697144
- 注册: 2012-03-17
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
7楼2013-05-07 09:20:03
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
8楼2013-05-07 11:00:33
【答案】应助回帖
★ ★ ★ ★ ★
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 |
» 本帖已获得的红花（最新10朵）
9楼2013-05-28 12:50:13
WJ-MS
新虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1559.2
- 散金: 100
- 帖子: 194
- 在线: 92.5小时
- 虫号: 1697144
- 注册: 2012-03-17
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
10楼2013-05-28 16:21:33
5