计算化学公社

标题: 关于SOS计算超极化率的问题 [打印本页]

作者
Author: cornerj 时间: 2018-6-26 17:42
标题: 关于SOS计算超极化率的问题
最近做第一超极化率的计算，碰到几个小问题，有点不明白，特来请教。
对于六十个原子左右（含CNHF）的共轭配合物体系，如果利用SOS方法计算其第一超极化率，选取CAM-B3LYP 泛函，按照培训讲义推荐的基组选用aug-cc-pVDZ。但是计算的体系中如若包含一个金属原子 In，似乎aug-cc-pVDZ并没有定义它，该换成什么基组？另外，对于不同的分子，结合Multiwfn计算出的β在偶极矩方向上的分量β||，有的是负值有的是正值比如-1000a.u.和300a.u.，这个正负号表示的是什么含义？如果是要比较其大小，是比较数值还是绝对值大小？

作者
Author: sobereva 时间: 2018-6-26 22:55
用aug-cc-pVDZ-PP
谈谈赝势基组的选用
http://sobereva.com/373

根据超极化率导数形式的定义理解。比绝对值

PS：Multiwfn 3.5的SOS功能有bug，如果用的3.5，改用最新的3.6(dev)或3.5之前的版本

作者
Author: cornerj 时间: 2018-6-27 21:21

sobereva 发表于 2018-6-26 22:55
用aug-cc-pVDZ-PP
谈谈赝势基组的选用
http://sobereva.com/373

哦哦，好的谢谢老师！
我正是用的Multiwfn3.5版本计算的β值，那我重新计算一遍。
另，发现计算过程中生成S0S.txt文件这一步骤极度耗费时间，对于正在计算的体系，笔记本要跑三四十小时甚至更久才能完成，除了增大settings.ini中的nthreads，不知有没有其他加快生成文件速度的方法？

作者
Author: sobereva 时间: 2018-6-28 02:07

cornerj 发表于 2018-6-27 21:21
哦哦，好的谢谢老师！
我正是用的Multiwfn3.5版本计算的β值，那我重新计算一遍。
另，发现计算过程中 ...

我刚刚更新了主页上的3.6(dev)版，这个版本在生成SOS.txt的过程中把非常耗时的一步（产生分子轨道间的跃迁偶极矩）的耗时降低了甚巨，可以试试。如果还慢，那就只能用更小点基组，或者把IOp(9/40=x)的x设更小点（但不能小于3），代码已经没法进一步优化了。

作者
Author: cornerj 时间: 2018-6-28 23:30

sobereva 发表于 2018-6-28 02:07
我刚刚更新了主页上的3.6(dev)版，这个版本在生成SOS.txt的过程中把非常耗时的一步（产生分子轨道间的 ...

十分感谢老师的指点！今天尝试了下用3.6（dev）版重新生成S0S.txt文件，速度快的飞起！太爽了！
大致耗时仅仅半小时至一小时而已，极度省时了！
麻烦再问点问题，对比了一下3.5版和3.6（dev）版得到的数据，对于同样的.fch和.log文件，体系的极化率α数值是可重现的，一模一样，
但β|| 的数值却有较大差别，如有略微增大几百a.u.的,也有暴减的，从之前+300a.u.变为-20a.u,对这个变化很是疑惑，不知是何原因？

作者
Author: sobereva 时间: 2018-6-29 00:12

cornerj 发表于 2018-6-28 23:30
十分感谢老师的指点！今天尝试了下用3.6（dev）版重新生成S0S.txt文件，速度快的飞起！太爽了！
大致耗 ...

3.5版的alpha没问题，beta和delta是错的

作者
Author: cornerj 时间: 2018-6-29 19:40

sobereva 发表于 2018-6-29 00:12
3.5版的alpha没问题，beta和delta是错的

明白了！谢谢老师！

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