24小时热门版块排行榜    

查看: 10656  |  回复: 8
【悬赏金币】回答本帖问题,作者yszheng将赠送您 22 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yszheng

新虫 (初入文坛)

[求助] Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?

在10x10x10的盒子里,做单个H2分子动力学模拟,目标温度为30k。在input.fdf文件里只改变了H2的坐标(键长不变,仅仅旋转一定角度),控温效果却完全不同: 两个氢原子沿着y方向,温度在8到31 k范围内波动;旋转为沿着z方向,温度在28到31 k范围内波动。我用了不同的热浴方法,结果差不多。单个水分子也存在这种情况。有哪位大佬知道是什么原因?

Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-1


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-2


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-3


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-4


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-5


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-6


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-7


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-8


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-9


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-10


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-11


Siesta 分子动力学模拟,10x10x10的盒子里将氢分子旋转90度,控温效果为什么会不同?-12


发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yszheng

新虫 (初入文坛)

不太明白“基组试一下用平滑的限制势产生”

发自小木虫Android客户端
8楼2021-05-14 15:54:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

暂时似乎难判断。是否尝试过对分子平移很小的一个长度,看看是否有差异?如果仍然差异明显我觉得可能是密度网格或其他精度不足引起的。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2021-05-12 00:42:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yszheng

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 卡开发发 at 2021-05-12 00:42:15
暂时似乎难判断。是否尝试过对分子平移很小的一个长度,看看是否有差异?如果仍然差异明显我觉得可能是密度网格或其他精度不足引起的。

必须有15金币才能发帖。我测试过k点,cutoff等,对结果基本没有影响!

发自小木虫Android客户端
3楼2021-05-14 09:09:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yszheng

新虫 (初入文坛)

我还跑过2个H2和2个H2O分子的,也会出现上面的情况

发自小木虫Android客户端
4楼2021-05-14 09:14:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-10-01 21:07 by Equinoxhua
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 IDs3scOF0tjC 2026-09-28 8/400 2026-10-01 19:08 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 CfXuS1rDhLYN 2026-09-28 4/200 2026-10-01 09:23 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见