| 查看: 2757 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xiaowu787木虫 (正式写手)
|
[交流]
【讨论】lammps-怪异的group温度 已有4人参与
|
||
|
说明一下，双金属体系，时间设为1fs，体系整体从0.05到300K升温，发现group的温度与体系的温度对不上，并且处于模型最外端的group的温度更是异常，比其他group的温度高几十度。 原来不小心把时间步长取为10fs，进行模拟同样的模拟也没有出现类似的问题。并且可以认定的是势函数是可靠的。不知是哪里的问题？ 下面是数据： 时间步长为0.001 单位为metal fix 2 all nvt temp 0.05 300.0 0.05 第一个数是整个体系的温度，后面五个数据分别是五个group的温度（把长方体模型均分为9个部分,取的五个group分别是从模型中心到模型一边的边界---相当于中心对称） 体系温度 group1 group2 group3 group4 group5 0.32223265 0.35672397 0.35130176 0.43155934 0.49054604 0.59470303 14.681225 17.707719 18.038659 19.401222 18.578224 23.696852 17.908917 22.595134 22.381331 22.859198 22.735383 29.888281 21.236465 26.594786 27.373716 27.549981 25.78319 33.537153 24.658653 32.005349 29.70771 29.110636 29.385247 40.579523 28.866911 35.896002 36.528466 36.668734 37.459447 45.291189 33.150806 40.469357 40.669588 40.651187 42.624405 55.096068 38.69363 47.161347 47.656391 48.248048 49.56008 65.557493 43.565413 57.637558 52.837027 55.77367 52.481708 71.654626 295.43779 314.76351 305.85202 324.41969 308.7731 386.26485 299.80219 313.80556 316.32717 320.05518 311.6842 389.04083 [ Last edited by xiaowu787 on 2010-10-31 at 09:16 ] |
回复此楼» 猜你喜欢
收到复试调剂但是去不了
已经有8人回复
-
260求调剂
已经有5人回复
-
急需调剂
已经有10人回复
-
申博/考博
已经有4人回复
-
化工学硕294分，求导师收留
已经有37人回复
-
求调剂
已经有11人回复
-
一志愿华中农业071010，320求调剂
已经有19人回复
-
304求调剂
已经有7人回复
-
求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师
已经有3人回复
-
二苯甲酮酸类衍生物
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MS生成Lammps的data文件，文件中原子的类型和顺序并不是自己想要的，请问如何修改？
已经有8人回复
- 求助LAMMPS关于MSD计算的问题
已经有17人回复
- LAMMPS温度－步长曲线震荡
已经有11人回复
- lammps的in文件
已经有4人回复
- lammps deposit 粒子初始速度为0？
已经有12人回复
- LAMMPS中meam势无法调用
已经有10人回复
- 怎么样开始LAMMPS的学习呢
已经有18人回复
- lammps的时间步长问题
已经有16人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 在LAMMPS中怎么建立图示的模型
已经有6人回复
- 【求助】LAMMPS运行错误：Out of range atoms - cannot compute PPPM
已经有7人回复
- 【求助】fedora中LAMMPS编译问题
已经有4人回复
- 【分享】尝试lammps中, 分享中...
已经有15人回复
- 【讨论】lammps中利用velocity命令无法控制住速度问题
已经有11人回复
- 【讨论】lammps中region和group的问题
已经有4人回复
- 【讨论】lammps中关于fix ave/spatial命令的输出问题
已经有4人回复
- 【讨论】lammps弛豫过程离奇的快
已经有4人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+5):谢谢老虎兄 2010-11-02 20:25:19
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+5):谢谢老虎兄 2010-11-02 20:25:19
|
我的意思不是说温度与原子数目有关，而是当你考虑总平均温度（即所有原子的总动能对总原子数的商 Ek(tot)/N），与各部分平均温度的平均（1/9*(Ek1/N1+Ek2/N2+...)），它们两个之间的关系时，这个关系（是相等，还是大于，还是小于，这个应该是与各部分原子数是否相等有关的 （上述两式是个大致的意思，不是温度的确切表达，没有写上有关的系数） 总之我基本上断定你两边的温度不会是对称的。 [ Last edited by 老虎大王 on 2010-11-2 at 19:32 ] |
6楼2010-11-02 19:23:22
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2010-10-30 19:31:54
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.3小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
3楼2010-10-31 09:28:07
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+6):谢谢 2010-11-01 09:32:08
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+6):谢谢 2010-11-01 09:32:08
|
你的九个部分原子数相等么？ A 如果原子数相等，且温度是对称的，那就意味着九个部分的平均温度都高于总平均温度。这个情况应是不可能出现的。所以温度可能不是对称的。 B 假如你九个部分原子数不相等，则情况比较复杂。 实际上，你在MD当中，分成九个部分，即使原子数相等，且初始结构对称，但你的初速度的方向和大小都是随机的，尤其是元胞中的中间部分经过运动变化后，也不会是完全不偏不倚的，所以温度一般不会完全对称。一边高一边低的可能性不能排除，这与你的边界条件也有关系。 [ Last edited by 老虎大王 on 2010-10-31 at 12:16 ] |
4楼2010-10-31 12:12:26
20