计算化学公社

标题: 求教ELF绘图问题 [打印本页]

作者
Author: emma_9821 时间: 2017-9-30 09:14
标题: 求教ELF绘图问题
sob老师，您好：
两个Eu原子用SDD基组，默认SDD（6d，10f）正常结束后，用multiwfn在XZ平面内绘制ELF图，此时需输入Y值，当Y=0.1和Y=3是得到不同的ELF， Y=0.1时得到的ELF更细致，能看到离核近处电子定域情况（右图）。而Y=3时（左图）离核近处则是无明显的定域离域情况，请问Y值选取的物理意义是什么呢？如果是3个原子，在用multiwfn时，选择某个分子平面，则无法调节垂直于该平面方向的值（无该选项），若想调节，该怎么办呢？

(, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 27)

作者
Author: sobereva 时间: 2017-9-30 14:59
绘制的是一个截面，显然得让这个截面取在最能反映体系电子结构特征的那个面上。比如Eu2这个体系，如果只绘制一个面的话，那么最能展现体系特征的那个面显然就是沿着键轴方向并且切开分子轴的面，即Y=0的XY面。而Y=3则偏离分子轴挺远了。

三个原子的时候通常就在Multiwfn里通过三个原子来定义作图平面。

作者
Author: emma_9821 时间: 2017-10-9 12:41
(, 下载次数 Times of downloads: 22) 谢谢老师的解答！还遇到一个问题：
三个Eu原子，同样的优化好的结构，用Gaussian自带的SDD关键字，默认Standard basis: SDD (6D, 10F)，再用stuttgart的网站http://www.tc.uni-koeln.de/PP/clickpse.en.html的小核赝势ECP28MWB和赝势基组，看了一下输出文件，读入的是（5d, 7f），得到的能量，前者是-2130.6990546，后者是-2131.018341，我该用哪个基组呢？
还有，我用两种基组得到的wfn文件选用3个原子所在平面绘制ELF图，直接用SDD关键字和用 stuttgart这两种情况得到的图也不一样，上图是直接用SDD关键字得到的ELF图（更细致），下图这个是用stuttgart的网站小核赝势ECP28MWB和赝势基组得到的ELF，如何解释这两者不同呢？或者是我在做图时，可还要调整什么参数吗？望老师能指点迷津，万分感谢！
(, 下载次数 Times of downloads: 25)

作者
Author: sobereva 时间: 2017-10-9 16:22
你在高斯里用gfinput打印一下赝势基组参数和stuttgart网站上进行对比来了解差异的原因，如果基组定义都相同，那就是笛卡尔和球谐函数的差异所致，用哪个能量都可以，所有计算保持一致即可。

光从图上不好说什么差异导致的，感觉基组定义差了很多。反正你用stuttgart网站上的是绝对没问题的

作者
Author: emma_9821 时间: 2017-10-9 17:04

emma_9821 发表于 2017-10-9 12:41
谢谢老师的解答！还遇到一个问题：
三个Eu原子，同样的优化好的结构，用Gaussian自带的SDD关键字，默认Sta ...

嗯，非常感谢！

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