计算化学公社

标题: 求助贴：Dmol3做汞吸附，吸附后汞原子的电荷计算结果异常偏高，参数设置是否存在问题 [打印本页]

作者
Author: FreesiaCC 时间: 2019-4-19 15:51
标题: 求助贴：Dmol3做汞吸附，吸附后汞原子的电荷计算结果异常偏高，参数设置是否存在问题
尝试计算Hg原子在V2O5/TiO2表面的吸附情况，得到吸附构型如下图，其中最上方的原子为Hg，所示键长为两个Hg-O键，算得Hg的吸附能为-32.79kJ/mol，population分析结果显示吸附后Hg原子的Mulliken atomic charges为0.237，Hirshfeld charge为0.2290。然而查阅文献中类似的Hg原子在V2O5/TiO2表面的吸附结构中Hg原子的电荷仅为0.07和0.019，虽然所用软件和参数有一定差别，但数值相差实在巨大，而O原子V原子的电荷与文献结果相差较小，是否我设置的某些参数并不适用于Hg原子的计算呢？以下为本次计算过程中设置的主要参数：MS软件Dmol3模块，GGA-PBE方法，使用TS方法进行DFT-D校正，选用spin unrestricted，Multiplicity为triplet，DNP基组，ECP core treatment，k点2*4*1，orbital cutoff取5.2，使用偶极校正，mulliken analysis选择atomic charge，Hirshfeld analysis选择charge。相关结果文件已上传附件，请各位老师同学批评指正。
另外，请问mulliken charge和Hirshfeld charge哪一个更适用于描述吸附前后电荷变化情况呢？谢谢大家
(, 下载次数 Times of downloads: 46) (, 下载次数 Times of downloads: 3)

(, 下载次数 Times of downloads: 4)

作者
Author: FreesiaCC 时间: 2019-4-22 12:53
希望各位老师同学有空的时候帮忙解答一下，万分感谢！！

作者
Author: sobereva 时间: 2019-4-23 01:50
Hirshfeld电荷物理意义比Mulliken明显更强，虽说有数量级偏小等的问题，但对于当前情况还是更值得使用

作者
Author: FreesiaCC 时间: 2019-4-23 09:51

sobereva 发表于 2019-4-23 01:50
Hirshfeld电荷物理意义比Mulliken明显更强，虽说有数量级偏小等的问题，但对于当前情况还是更值得使用

好的谢谢您，请问我在计算过程中设置的参数是否存在不妥呢，结果中V,O和文献中电荷仅相差在小数点后两位，而我细算得出的Hg原子电荷明显高于文献，由于吸附并不是很强，我认为计算得出的Hg电荷偏高了，是否我采用的某些参数不适用于描述Hg原子呢

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