| 查看: 2734 | 回复: 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 ] |
回复此楼» 猜你喜欢
322求调剂
已经有12人回复
-
电子信息270求调剂
已经有16人回复
-
272分材料子求调剂
已经有50人回复
-
295分求调剂
已经有13人回复
-
339求调剂
已经有9人回复
-
一志愿哈工大 085600 277 12材科基求调剂
已经有28人回复
-
调剂求收留
已经有28人回复
-
290调剂生物0860
已经有26人回复
-
2本，初试303，0860求调剂
已经有7人回复
-
材料工程281还有调剂机会吗
已经有28人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
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
- 专业: 金属结构材料
2楼2010-10-30 19:31:54
老虎大王
木虫 (著名写手)
- 模拟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
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1406.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.3小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
★
zh1987hs(金币+1):鼓励交流 2010-11-02 20:24:48
zh1987hs(金币+1):鼓励交流 2010-11-02 20:24:48
|
group temp0 region temp0 47528 atoms in group temp0 group temp1 region temp1 47534 atoms in group temp1 group temp2 region temp2 47525 atoms in group temp2 group temp3 region temp3 47518 atoms in group temp3 group temp4 region temp4 48582 atoms in group temp4 谢谢您的答复，现在还有一点不明白，温度不是统计平均值吗？为什么还与group原子数目有关？希望您能再指点一下？ |
5楼2010-11-02 17:44:29
