doublezhang 发表于 2020-6-9 09:33 有解析导数就完全没有必要自己用有限场,后者只会更慢、更麻烦,而且精度还低(尤其是纯数值导数,到了高阶时累计误差容易很大),除非仅仅是需要(超)极化率某个分量。 有解析导数却还用有限场的基本都是没搞清楚什么叫解析导数和数值导数。 |
sobereva 发表于 2020-6-8 18:57 如果是算静态超极化率的话,有限场方法会快一点,就是体系大的时候加了电场scf有可能不收敛。 谢谢社长,社长这两篇文章我早就日日拜读,参悟出来的大道都够我们全组同学达到毕业要求啦 ![]() |
doublezhang 发表于 2020-6-7 20:33 如今Gaussian都有现成的polar关键词。完全不需要自己用有限场方法计算 看 使用Multiwfn分析Gaussian的极化率、超极化率的输出 http://sobereva.com/231 另外还可以用SOS来算 使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率 http://sobereva.com/232 |
| 算的第一超极化率密度。调整好坐标之后,就可得到想要的结果了。 |
瑶瑶 发表于 2020-6-6 18:58 有很多几十年前年前的文献讲怎么用有限场算超极化率的分量的,感觉很复杂,不仅仅是几个外电场的事 |
sobereva 发表于 2020-6-6 13:52 上图就是X轴冲着偶极矩方向的,计算时加了nosymm,xyy分量最大。刚看了一个文献,算的是βyyz分量的-zρyy。这样算是不是在y轴和z都加外电场呢。 |
| 谢谢老师的回复 |
| 旋转体系,让X轴冲着偶极矩方向,或者说让体系的偶极矩正好平行于X方向,重算一遍(记得加上nosymm),看看届时XXX分量是不是最大的,如果是的话就研究XXX |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-19 18:04 , Processed in 0.174910 second(s), 25 queries , Gzip On.