应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【讨论】lammps-怪异的group温度
51/1返回列表
查看: 2446  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-10-30 17:11:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

453695608

银虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主你好，你这体系温度是怎么弄出来的啊，还有你这是多少步升多少温度啊
赞一下 回复此楼
学习更多的知识，为了未来
10楼2015-10-11 11:13:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

老虎大王

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你把11个group温度都搞出来看看再说。
赞一下(1人) 回复此楼
2楼2010-10-30 19:31:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)
★
zh1987hs(金币+1): 鼓励交流~ 2010-11-02 20:24:36
引用回帖:
Originally posted by 老虎大王 at 2010-10-30 19:31:54:
你把11个group温度都搞出来看看再说。

整个模型是均分的9个部分，9个部分结构完全一样相当于是个超晶胞，取了包含中心部分的group和一边的四个group,没有输出另一边的。不过就这感觉已经不正常了，两边的温度应该是一样吧，不会出现一边高一边低吧。后面的模拟就是认为两边的温度是一样的，这有问题吗？
赞一下(1人) 回复此楼
3楼2010-10-31 09:28:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+6):谢谢 2010-11-01 09:32:08
你的九个部分原子数相等么？

A 如果原子数相等，且温度是对称的，那就意味着九个部分的平均温度都高于总平均温度。这个情况应是不可能出现的。所以温度可能不是对称的。
B 假如你九个部分原子数不相等，则情况比较复杂。


实际上，你在MD当中，分成九个部分，即使原子数相等，且初始结构对称，但你的初速度的方向和大小都是随机的，尤其是元胞中的中间部分经过运动变化后，也不会是完全不偏不倚的，所以温度一般不会完全对称。一边高一边低的可能性不能排除，这与你的边界条件也有关系。

[ Last edited by 老虎大王 on 2010-10-31 at 12:16 ]
赞一下(2人) 回复此楼
4楼2010-10-31 12:12:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭