vasp跑NVT thermostat统计原子displacement分布问题求助

Oliviaz · 发表于 Post on 2020-3-1 00:16:06

本帖最后由 Oliviaz 于 2020-3-1 18:08 编辑

用vasp跑Cu(111) slab的NVT系综，想看看在给定温度下表面原子的位移分布，但得到的分布却不是正态分布。想问问大家知不知道有哪些可能的出错原因，INCAR和相关图片见下。

备注：我想跑的是Andersen thermostat, 但是Nose-Hoover也试过了，都得不到正态分布。恳请大家不吝赐教，多谢大家！

get-it · 发表于 Post on 2020-3-1 07:59:44

只有当体系可以用简谐声子完全描述时，位移分布才是正态分布。而实际的相互作用都是非谐的。如果你一定要跑出正态分布，那就设一个很低的温度。

Oliviaz · 发表于 Post on 2020-3-1 15:51:59

get-it 发表于 2020-3-1 07:59
只有当体系可以用简谐声子完全描述时，位移分布才是正态分布。而实际的相互作用都是非谐的。如果你一定要跑 ...

你好，我这里设置的是300K，也不能出现正态分布，而且在九百甚至一千K下，都应该是正态分布的呢。

get-it · 发表于 Post on 2020-3-1 16:29:53

Oliviaz 发表于 2020-3-1 15:51
你好，我这里设置的是300K，也不能出现正态分布，而且在九百甚至一千K下，都应该是正态分布的呢。

请问，你这个“应该是正态分布的结论”是怎么得到的？

Oliviaz · 发表于 Post on 2020-3-1 18:00:36

本帖最后由 Oliviaz 于 2020-3-1 18:03 编辑

get-it 发表于 2020-3-1 16:29
请问，你这个“应该是正态分布的结论”是怎么得到的？

因为我用已有的势能面对不同温度的位移进行统计，产生的都是正态分布，包括我们组内其他的体系跑thermostat也得到了正态分布，所以我更怀疑是我跑的AIMD的问题~
图一就是我势能面的位移图了~

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] vasp跑NVT thermostat统计原子displacement分布问题求助