计算化学公社

标题: vasp跑NVT thermostat统计原子displacement分布问题求助 [打印本页]

作者
Author: Oliviaz 时间: 2020-3-1 00:16
标题: vasp跑NVT thermostat统计原子displacement分布问题求助
本帖最后由 Oliviaz 于 2020-3-1 18:08 编辑

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

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

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

作者
Author: Oliviaz 时间: 2020-3-1 15:51

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

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

作者
Author: get-it 时间: 2020-3-1 16:29

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

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

作者
Author: Oliviaz 时间: 2020-3-1 18:00
本帖最后由 Oliviaz 于 2020-3-1 18:03 编辑

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

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

欢迎光临计算化学公社 (http://ccc.keinsci.com/)