计算化学公社

标题: forcite跑动力学后黄铁矿吸附的水集中在一部分晶面是怎么回事？求助分子动力... [打印本页]

作者
Author: hgkjqx 时间: 2019-6-21 09:01
标题: forcite跑动力学后黄铁矿吸附的水集中在一部分晶面是怎么回事？求助分子动力...
如图，在几何优化后跑动力学之前水层均匀分布在黄铁矿表面，跑完动力学后有一部分水在一小部分晶面上向真空层方向跑，形成水柱，晶面是40*40埃，真空层100埃，UFF力场，黄铁矿和水是在castep中优化的，黄铁矿的电荷用的是dmol3中计算的mulliken电荷，晶面应该每个原子是一样的，为什么出现了作用不一致的情况呢？

作者
Author: sobereva 时间: 2019-6-21 10:49
Mulliken电荷严重不适合分子动力学目的
仔细看此文相关讨论
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

作者
Author: hgkjqx 时间: 2019-6-21 14:30

sobereva 发表于 2019-6-21 10:49
Mulliken电荷严重不适合分子动力学目的
仔细看此文相关讨论
RESP拟合静电势电荷的原理以及在Multiwfn中的 ...

请问bader电荷呢？

作者
Author: hgkjqx 时间: 2019-6-21 14:34
是看美国一个教授的论文这样操作的，因为这样做出来得到的接触角与实际差别较小，所以用的这个电荷

作者
Author: sobereva 时间: 2019-6-21 22:47

hgkjqx 发表于 2019-6-21 14:30
请问bader电荷呢？

Bader电荷更是糟糕透顶，结果毫无意义
仔细看了《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）里面静电势重现性的对比自然就明白了

作者
Author: sobereva 时间: 2019-6-21 22:47

hgkjqx 发表于 2019-6-21 14:34
是看美国一个教授的论文这样操作的，因为这样做出来得到的接触角与实际差别较小，所以用的这个电荷

别盲目效仿文献
就连JACS上的某些文章都不乏低级荒唐的理论计算角度的错误

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