计算化学公社

标题: 使用Multiwfn计算超极化率密度的分量问题 [打印本页]

作者
Author: 瑶瑶 时间: 2020-6-6 13:43
标题: 使用Multiwfn计算超极化率密度的分量问题
各位老师好。在用Multiwfn计算第一超极化率密度时，只关注最重要的一个分量，比如βzzz。但是，各个分量中如果xyy是最大的，该怎么办呢。

作者
Author: sobereva 时间: 2020-6-6 13:52
旋转体系，让X轴冲着偶极矩方向，或者说让体系的偶极矩正好平行于X方向，重算一遍（记得加上nosymm），看看届时XXX分量是不是最大的，如果是的话就研究XXX

作者
Author: 瑶瑶 时间: 2020-6-6 13:54
谢谢老师的回复

作者
Author: 瑶瑶 时间: 2020-6-6 18:58

sobereva 发表于 2020-6-6 13:52
旋转体系，让X轴冲着偶极矩方向，或者说让体系的偶极矩正好平行于X方向，重算一遍（记得加上nosymm），看看 ...

上图就是X轴冲着偶极矩方向的，计算时加了nosymm，xyy分量最大。刚看了一个文献，算的是βyyz分量的-zρyy。这样算是不是在y轴和z都加外电场呢。

作者
Author: doublezhang 时间: 2020-6-7 20:33

瑶瑶发表于 2020-6-6 18:58
上图就是X轴冲着偶极矩方向的，计算时加了nosymm，xyy分量最大。刚看了一个文献，算的是βyyz分量的-zρy ...

有很多几十年前年前的文献讲怎么用有限场算超极化率的分量的，感觉很复杂，不仅仅是几个外电场的事

作者
Author: 瑶瑶 时间: 2020-6-7 21:44
算的第一超极化率密度。调整好坐标之后，就可得到想要的结果了。

作者
Author: sobereva 时间: 2020-6-8 18:57

doublezhang 发表于 2020-6-7 20:33
有很多几十年前年前的文献讲怎么用有限场算超极化率的分量的，感觉很复杂，不仅仅是几个外电场的事

如今Gaussian都有现成的polar关键词。完全不需要自己用有限场方法计算
看
使用Multiwfn分析Gaussian的极化率、超极化率的输出
http://sobereva.com/231

另外还可以用SOS来算
使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率
http://sobereva.com/232

作者
Author: doublezhang 时间: 2020-6-9 09:33

sobereva 发表于 2020-6-8 18:57
如今Gaussian都有现成的polar关键词。完全不需要自己用有限场方法计算
看
使用Multiwfn分析Gaussian的 ...

如果是算静态超极化率的话，有限场方法会快一点，就是体系大的时候加了电场scf有可能不收敛。
谢谢社长，社长这两篇文章我早就日日拜读，参悟出来的大道都够我们全组同学达到毕业要求啦

作者
Author: sobereva 时间: 2020-6-9 15:28

doublezhang 发表于 2020-6-9 09:33
如果是算静态超极化率的话，有限场方法会快一点，就是体系大的时候加了电场scf有可能不收敛。
谢谢社长 ...

有解析导数就完全没有必要自己用有限场，后者只会更慢、更麻烦，而且精度还低（尤其是纯数值导数，到了高阶时累计误差容易很大），除非仅仅是需要(超)极化率某个分量。
有解析导数却还用有限场的基本都是没搞清楚什么叫解析导数和数值导数。

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