| 查看: 1826 | 回复: 7 | |||
| 【悬赏金币】回答本帖问题，作者yszheng将赠送您 22 个金币 | |||
[求助]
Siesta 分子动力学模拟，10x10x10的盒子里将氢分子旋转90度，控温效果为什么会不同？
|
|||
在10x10x10的盒子里，做单个H2分子动力学模拟，目标温度为30k。在input.fdf文件里只改变了H2的坐标（键长不变，仅仅旋转一定角度），控温效果却完全不同: 两个氢原子沿着y方向，温度在8到31 k范围内波动；旋转为沿着z方向，温度在28到31 k范围内波动。我用了不同的热浴方法，结果差不多。单个水分子也存在这种情况。有哪位大佬知道是什么原因？             发自小木虫Android客户端 |
回复此楼» 猜你喜欢
281求调剂（0805）
已经有8人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有3人回复
-
268求调剂
已经有3人回复
-
311求调剂
已经有10人回复
-
被我言中：新模板不强调格式了，假专家开始管格式了
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3508.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
2楼2021-05-12 00:42:15
3楼2021-05-14 09:09:54
4楼2021-05-14 09:14:39
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3508.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
5楼2021-05-14 09:53:30
|
k点密度不影响能量和力。但旋转前后力是不一样的。第一张是diffy（两个氢原子仅y不同），第二张是diff-xyz（两个氢原子xyz都不同）（键长一致）。diff-y的两个氢原子y方向受力要比diff-xyz的大许多。此外我还用vasp测试，氢分子旋转前后温度波动也是有很大差异的！    发自小木虫Android客户端 |
6楼2021-05-14 15:45:44
7楼2021-05-14 15:46:41
8楼2021-05-14 15:54:18
20