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

FreesiaCC · 发表于 Post on 2019-4-19 15:51:26

尝试计算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哪一个更适用于描述吸附前后电荷变化情况呢？谢谢大家

opt.txt (236.18 KB, 下载次数 Times of downloads: 3)

population.txt (51.62 KB, 下载次数 Times of downloads: 4)

FreesiaCC · 发表于 Post on 2019-4-22 12:53:37

希望各位老师同学有空的时候帮忙解答一下，万分感谢！！

sobereva · 发表于 Post on 2019-4-23 01:50:31

Hirshfeld电荷物理意义比Mulliken明显更强，虽说有数量级偏小等的问题，但对于当前情况还是更值得使用

FreesiaCC · 发表于 Post on 2019-4-23 09:51:53

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

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

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

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

浏览过的版块